YAP regulates cell mechanics by controlling focal adhesion assembly

Abstract

Hippo effectors YAP/TAZ act as on-off mechanosensing switches by sensing modifications in extracellular matrix (ECM) composition and mechanics. The regulation of their activity has been described by a hierarchical model in which elements of Hippo pathway are under the control of focal adhesions (FAs). Here we unveil the molecular mechanism by which cell spreading and RhoA GTPase activity control FA formation through YAP to stabilize the anchorage of the actin cytoskeleton to the cell membrane. This mechanism requires YAP co-transcriptional function and involves the activation of genes encoding for integrins and FA docking proteins. Tuning YAP transcriptional activity leads to the modification of cell mechanics, force development and adhesion strength, and determines cell shape, migration and differentiation. These results provide new insights into the mechanism of YAP mechanosensing activity and qualify this Hippo effector as the key determinant of cell mechanics in response to ECM cues.

ARTICLE

Received 8 Nov 2016 | Accepted 10 Mar 2017 | Published 15 May 2017

DOI: 10.1038/ncomms15321 OPEN

YAP regulates cell mechanics by controlling focal adhesion assembly

Giorgia Nardone1,*, Jorge Oliver-De La Cruz1,*, Jan Vrbsky1, Cecilia Martini1, Jan Pribyl2, Petr Skldal2, Martin Pel1,3, Guido Caluori1, Stefania Pagliari1, Fabiana Martino1,3, Zuzana Maceckova4, Marian Hajduch4, Andres Sanz-Garcia5,6, Nicola Maria Pugno7,8,9, Gorazd Bernard Stokin1 & Giancarlo Forte1,10

Hippo effectors YAP/TAZ act as onoff mechanosensing switches by sensing modications in extracellular matrix (ECM) composition and mechanics. The regulation of their activity has been described by a hierarchical model in which elements of Hippo pathway are under the control of focal adhesions (FAs). Here we unveil the molecular mechanism by which cell spreading and RhoA GTPase activity control FA formation through YAP to stabilize the anchorage of the actin cytoskeleton to the cell membrane. This mechanism requires YAP co-transcriptional function and involves the activation of genes encoding for integrins and FA docking proteins. Tuning YAP transcriptional activity leads to the modication of cell mechanics, force development and adhesion strength, and determines cell shape, migration and differentiation. These results provide new insights into the mechanism of YAP mechanosensing activity and qualify this Hippo effector as the key determinant of cell mechanics in response to ECM cues.

1 International Clinical Research Center (ICRC), St Annes University Hospital, CZ-65691 Brno, Czech Republic. 2 CEITEC MU, Masaryk University, CZ-65691 Brno, Czech Republic. 3 Faculty of Medicine, Department of Biology, Masaryk University, CZ-62500 Brno, Czech Republic. 4 Faculty of Medicine and Dentistry, Institute of Molecular and Translational Medicine, Palacky University Olomouc, Hnvotnsk 1333/5, CZ-77515 Olomouc, Czech Republic. 5 Faculty of Pharmacy, Division of Pharmaceutical Biosciences, University of Helsinki, Viikinkaari 5 E, FI-00014 Helsinki, Finland. 6 Institute of Advanced Biomedical Engineering and Science, Tokyo Womens Medical University, 8-1 Kawada-Cho, Shinjuku-Ku, Tokyo JP-162-8666, Japan. 7 Laboratory of Bio-inspired and Graphene Nanomechanics, Department of Civil, Environmental and Mechanical Engineering, University of Trento, I-38123 Trento, Italy. 8 Ket-Lab, Italian Space Agency, Via del Politecnico snc, 00133 Rome, Italy. 9 School of Engineering and Materials Science, Queen Mary University of London, Mile End Road, UK-E1 4NS London, UK. 10 Department of Biomaterials Science, Institute of Dentistry, University of Turku, FI-20014 Turku, Finland. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to G.F. (e-mail: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 8:15321 | DOI: 10.1038/ncomms15321 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15321

Cells are in constant isometric tension with the extracellular matrix (ECM), an equilibrium of forces needed to ensure to adopt the shape and volume suited to exert

their function1,2. On a larger scale, this condition keeps organ functionality, while changes in the mechanical balance between the cells and the surrounding milieu result in tissue malfunctioning or malignant transformation3,4.

The ability of cells to perceive ECM mechanics and spread is associated to Hippo pathway effectors Yes-associated protein (YAP) and WW domain-containing transcription regulator protein 1 (WWTR1 or TAZ) shuttling to the nucleus to exert their co-transcriptional activity5,6. By binding to cell- and context-specic transcription factors, YAP/TAZ contribute to ECM remodelling79.

Focal adhesions (FAs), the main hub for cell mechanosensing, act as a bridge between integrin-ECM connection and the cytoskeleton10. Changes in the signals propagated through FAs have been reported in malignant cells and are essential for tumour cell spreading11.

YAP/TAZ nuclear activity is correlated to the stability of actin cytoskeleton and cell tension, as controlled by myosin light chain II and Rho GTPase pathways1214. Integrin-FA signalling has been recently suggested to control Hippo pathway by phosphorylating large tumour suppressor (LATS) kinases through Src15. These results predicted a hierarchical mechanism by which Hippo effectors behave as downstream sensors of ECM mechanics through integrin-FA signalling and by perceiving cytoskeleton stability. Here we describe the molecular basis of the crosstalk among the different cell mechanosensing systems and propose a model by which YAP directly regulates FA assembly and cell mechanics.

ResultsCell area controls YAP shuttling regardless of FA assembly. Considering recent evidence suggests possible interplay between Hippo pathway and FAs1517, we investigated the correlation between YAP nuclear localization and the presence of FAs. To this end, we cultured adipose tissue-derived mesenchymal stem cells (AD-MSCs) onto bronectin (FN)-coated elastically supported surfaces of different stiffness (28 and 1.5 kPa) or onto glass surfaces coated either with FN or poly-L-lysine (PLL). FN coating onto the stiff surface (28 kPa) promotes FA assembly, whereas the exposure of cells to PLL abrogates FA formation regardless of substrate stiffness. In addition, FN is not able to foster FA arrangement on soft (1.5 kPa) surfaces. Interestingly, in all the conditions in which FA formation was prevented (PLL or soft surfaces), YAP was excluded from the nucleus and cell area was signicantly reduced as compared with the controls, consistent with the absence of FAs in the cell18 (n 10,

3 technical replicates, MannWhitney test, ****Po0.0001 and **Po0,001; Fig. 1a). Remarkably, when single cell area was controlled by growing AD-MSCs onto FN-coated micropatterns forcing the cell to acquire given sizes (300, 1,024, 2,025 and 10,000 mm2), YAP nuclear expression correlated with the presence of vinculin spikes. In fact, both YAP nucleus/ cytoplasm ratio and the number of FAs increased steadily with the size of the cell. Moreover, the expression of FA genes vinculin (VCL), zyxin (ZYX), talin 1 (TLN1) and talin 2 (TLN2) was found increased in cells having larger (10,000 mm2) as compared with smaller surface areas (1,024 mm2). Spread cells also displayed enhanced RHOA expression as compared with the conned controls. These results suggested that the expression of FA-related genes and the formation of the multiprotein complexes were correlated to YAP nuclear localization and controlled by cell spreading (Fig. 1b).

Given the correlation between FA presence and nuclear YAP, we asked whether YAP nuclear accumulation was determined by cell size or by the formation of FAs. Therefore, we overexpressed two main FA componentsvinculin and zyxinin AD-MSCs and detected no signicant change in nuclear YAP and cell surface area (Fig. 1c). Moreover, we designed single cell tools based on surface micropatterning by which cell area could be tailored independently of the availability of FN spots. By this means, cells were induced to maintain the same surface area while contacting a different number of adhesion sites. This experimental setting allowed us to nely tune FA number regardless of cell size. Within the same experiments, we were also able to compare cells having different cell area but constant adhesion area, as mediated by FN distribution (Fig. 1d). As expected, when the adhesion area was reduced (4,900 versus 1,024 mm2), the number of FAs dramatically dropped (28438 versus 8821, means.d.) and the ratio between nuclear and cytoplasmic YAP decreased accordingly (5.190.88 versus 2.991.44). On the contrary, when cell area was kept constant (4,900 mm2) but the adhesion area decreased (1,000 and 450 mm2), no signicant effect on YAP nuclear localization was detected (5.200.88 versus4.660.94 or 4.321.08, respectively), independently of the reduction in FA number (from 283.938.14 to 207.020.30 and 210.736.80, respectively; Fig. 1e). When the adhesion area was reduced to 254 mm2 in cells preserving constant area (4,900 mm2), no FA spike organization was detected, while YAP nuclear localization was not signicantly affected (n 10, 3 technical

replicates, Dunns test, ****Po0.0001, ***Po0.001 and **Po0.01). Lower adhesion areas (176 mm2) did not allow cell adhesion. These results pointed at YAP nuclear accumulation as a function of cell area rather than FA formation, as controlled by ECM adhesion. Moreover, they raised the possibility that the acquisition of cell area could directly affect YAP localization and function, independently of FA assembly.

YAP tunes FAs and cell stiffness in response to cell area. To explore the connection among cell area, YAP and FAs, we knocked down YAP in AD-MSCs by short hairpin RNA (shYAP) and controlled cell area through FN-coated micropatterns. The same experiment was performed with cells in which YAP paralogue protein TAZ (shTAZ) was silenced, as cells grown onto PLL (lacking FAs) displayed a signicant reduction in TAZ nuclear localization as compared with those cultured on FN (Supplementary Fig. 1). The selective silencing of YAP or TAZ was conrmed in quantitative PCR analysis (Supplementary Fig. 2a).

Although control cells (CTR FN) were able to acquire the given size and shape over FN micropatterns, shYAP cells displayed impaired FA maturation and acquired a bulging phenotype with deranged F-actin cytoskeleton protruding in spikes. On the contrary, TAZ-deprived cells preserved a correct cytoskeleton arrangement and developed FAs, while losing their ability to spread on the given surface: (Fig. 2a, Supplementary Fig. 2b and Supplementary Movies 13). Given these results, we decided to impair both YAP and TAZ nuclear localization in the cells, by culturing them onto N-cadherin micropatterns. The use of this experimental condition was suggested by the evidence that the concomitant silencing of both YAP and TAZ did not allow the cells to adhere.

Consistent with the model of cadherincatenin system sequestering both YAP and TAZ at the cellcell interaction site19,20, AD-MSCs grown onto N-cadherin-coated surface and stained with an antibody recognizing both YAP and TAZ displayed no nuclear localization of the proteins. In such conditions, FA formation was completely abrogated, cytoskeleton

2 NATURE COMMUNICATIONS | 8:15321 | DOI: 10.1038/ncomms15321 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15321 ARTICLE

a

b

28 kPa 1.5 kPa FN PLL 300

1,024

Total FA area (m)

Fold regulation

(10,000 mversus 1,024 m)

FA number

YAPF-actinvinculinDAPI

YAP

vinculin

****

**

Cell area (m ) Cell area (m )

2,025

10,000

YAP nuclear/cyto

ratio

****

**

10 m

20 m

Cell area (m )

c

DAPI YAP F-actin

YAP nuclear/cyto

ratio

YAP nuclear/cyto

ratio

CTR Zyxin CTR Vinculin

1

2

vinculin

zyxin

1

1 2

2

Cell area (m)

Cell area (m)

1

2

d e

DAPI vinculin F-actin FN DAPI YAP F-actin FN

***

***

FA number

Total FA area (m)

YAP nuclear/cyto

ratio

****

****

**

****

Figure 1 | YAP nuclear shuttling is determined by cell area. (a) Confocal images of AD-MSCs cultured onto elastically supported surfaces (ESS) with 28,1.5 kPa or bronectin- and poly-L-lysine-coated glass slides. Cells were stained with the indicated antibodies: anti-vinculin (green), anti-YAP (red). F-actin was decorated with Alexa Fluor 647 Phalloidin (white) and nuclei were counterstained with 40,6-diamidino-2-phenylindole (DAPI). Graphs: quantication of the number of FAs per cell (FA number), cell area and YAP nuclear/cytoplasmic ratio in AD-MSCs cultured on ESS, bronectin or poly-L-lysine(n 10, 3 technical replicates, MannWhitney test, ****Po0.0001 and **Po0,001). (b) Confocal analysis of single AD-MSCs cultured onto

micropatterned glass slides with 300, 1,024, 2,025 and 10,000 mm2 areas and coated with bronectin. Cells were stained with anti-YAP (red), and anti-vinculin (green). Graphs: quantication of FA number, total FA area and YAP nuclear/cytoplasmic ratio in single AD-MSCs grown onto bronectin-coated surfaces with increasing adhesion areas. Bottom right: Quantitative RTPCR analysis of vinculin (VCL), zyxin (ZYX), talin1 (TLN1), talin2 (TLN2), RHOA genes in AD-MSCs cultured onto 10,000 versus 1,024 mm2. The results are expressed as mean fold regulation obtained in 2 independent experiments and the bar indicates the s.d. (c) Confocal images and analysis of AD-MSCs transfected with vinculin-venus or RFP-zyxin and stained for YAP. F-actin was decorated with Alexa Fluor 647 Phalloidin (white). Graphs: quantication of YAP nucleus/cytoplasm ratio and cell area in transfected cells. (d) Left: Annotation of micropattern properties (cell area, adhesion area), schematic side view of single cell onto bronectin-coated micropatterned glass slides, schematic top view of bronectin distribution in the micropatterns. Fibronectin-covered area is indicated in blue. Center: confocal images of single AD-MSCs grown onto micropatterns stained with anti-vinculin (green). Right: confocal images of single AD-MSCs grown onto micropatterns stained with anti-YAP (green). In both panels cells are stained with anti-bronectin (white), Alexa Fluor 546 Phalloidin and DAPI. (e) Graphs: quantication of FA number, total FA area and YAP nuclear/cytoplasmic ratio in single cells cultured onto bronectin-coated surfaces. All error bars are s.d. (n 10, 3 technical

replicates, ****Po0.0001, *** Po0.001 and **Po0.01, as calculated by KruskalWallis test followed by post hoc Dunns test for multiple comparison).

was disrupted and cells failed to acquire the given shape and size (Fig. 2a and Supplementary Fig. 2c). FAs and cytoskeleton are required to control cell mechanics in response to ECM cues10,21. Therefore, we decided to measure the stiffness of shYAP and

shTAZ AD-MSCs by atomic force microscopy (AFM). The analysis of single cell Youngs Modulus claried that both shYAP and shTAZ cells were signicantly less stiff than the control and comparable to cells lacking FAs (Fig. 2b). Elasticity value

NATURE COMMUNICATIONS | 8:15321 | DOI: 10.1038/ncomms15321 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15321

a FNCTR

F-actin vinculin DAPI

shYAP shTAZ N-cadherin

CTR FN

d

CTR-N-cadherin

300 1,024 2,025 10,000

YAP

YAP

YAP

TAZ

TAZ

TAZ

shYAP FN shTAZ FN

b c

Youngs modulus (Pa)

15,000

10,000

5,000

0

shYAP

shTAZ

****

e AD-MSCs CAL51

DAPI YAP mutant

WT WT

PDZ

304

vinculin

CTR

S127A

DAPI vinculin

AD-MSCs

PDZ

304

S127A

WT

S89A

S127A

WT

S89A

DAPI TAZ mutant vinculin

DAPI YAP mutant vinculin F-actin

DAPI YAP mutant vinculin F-actin

CAL51 CTR

S127A

DAPI vinculin

Figure 2 | YAP regulates FA formation and cell stiffness in response to cell area in AD-MSCs. (a) Left: confocal analysis of individual cells grown onto 300, 1,024, 2,025, 10,000 mm2 micropatterns coated with either bronectin or N-cadherin. Right: three-dimensional (3D) reconstruction of individual cells grown onto 10,000 mm2. Cells were stained with: anti-vinculin (green) Alexa Fluor 546 Phalloidin and 40,6-diamidino-2-phenylindole (DAPI). Cells were infected with viral particles encoding for either control short hairpin RNA (CTR), short-hairpin RNA targeting YAP or TAZ (shYAP, shTAZ, respectively).

(b) Graph: Youngs modulus analysis of single AD-MSCs, shYAP and shTAZ onto micropatterns coated with FN or PLL. Values are shown as medianmin/max (n 24, 3 technical replicates, KruskalWallis test followed by Dunns test, ****Po0.0001). (c) Venn diagram showing the overlap

between genes signicantly regulated in shYAP and shTAZ cells as obtained by quantitative PCR arrays (see Supplementary Fig. 3a). (d) Schematic structure of YAP (S127A, DPDZ) and TAZ (S89A, D304) mutants used for the transfections. (e) Confocal images of focal adhesion distribution, as stained by anti-vinculin (red) and DAPI in AD-MSCs (left) and CAL51 (centre) transfected with the indicated YAP and TAZ mutants. Right: Representative image analysis of vinculin expression in AD-MSCs and CAL51 transfected or not with S127A YAP mutant. Transfected AD-MSCs were stained with anti-vinculin antibody (red). Transfected CAL51 cells were stained with anti-vinculin (red), Alexa Fluor 647 Phalloidin (white) and DAPI. Transfected cells are identied by GFP expression and highlighted by dashed lines.

for control cells grown onto FN-coated glass surfaces was 11,5181,786.5 Pa, whereas values calculated for shYAP, shTAZ on FN, and control cells cultured onto PLL-coated glass were 371200.35, 1,20570.5 and 755.7358.45 Pa, respectively (n 24, 3 technical replicates, Dunns test,

****Po0.0001). Owing to cell detaching from the surface,

it was not possible to measure neither cells in which both YAP and TAZ were depleted nor cells grown onto N-cadherin micropatterns.

As YAP- and TAZ-silenced cells displayed different phenotypes, we analysed the transcriptional activity of the two proteins as well as their impact on cell differentiation potential.

4 NATURE COMMUNICATIONS | 8:15321 | DOI: 10.1038/ncomms15321 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15321 ARTICLE

The paralogue proteins exhibited different transcriptional proles, like shown by PCR array analysis (Fig. 2c and Supplementary Fig. 3a,b, the full gene list is provided in Supplementary Data 1) and distinct adipogenic differentiation potential as compared to the control cells and among themselves (Supplementary Fig. 3b). Altogether, the results indicated thatbesides acting as mechanosensors by perceiving ECM mechanics, similar to that previously shown7YAP and TAZ are responsible of regulating cell stiffness in AD-MSCs. Moreover, the data suggested that YAP activity was exerted by controlling the assembly of FAs.

YAP/TAZ are known to exert their function as transcriptional co-activators mainly by binding TEAD transcription factor22. Therefore, we transiently overexpressed in two different cell types (AD-MSCs and Cal51) constructs coding for mutants of YAP and TAZ having a nuclear (YAPS127A, TAZS89A) or cytoplasmic (YAPDPDZ, TAZD304) localization2325. YAPS127A and TAZS89A mutants cannot be phosphorylated by LATS1/2 kinases and thus are predominantly located in cell nucleus and transcriptionally active26. On the other hand, YAPDPDZ and TAZD304 mutants lack the PDZ domain, which is required for nuclear shuttling22 and thus are sequestered to the cytoplasm (Fig. 2d). Staining transfected cells for FA protein vinculin highlighted an increase in the signal and an altered distribution of FA spikes induced by the transcriptionally active YAP mutant (YAPS127A). On the contrary, no change in FA expression was triggered by TAZ constitutively active mutant (TAZS89A) and the cytoplasmic mutants of both proteins (YAPDPDZ and TAZD304) (Fig. 2e).

Overall, these results suggested that the nuclear presence and the co-transcriptional activity of YAP were requested to promote FA formation.

YAP directly regulates the expression of FA-related genes. To understand whether YAP nuclear activity determines a direct activation of FA genes, we performed chromatin immunoprecipitation followed by deep sequencing (ChIP-Seq) analysis of YAP-specic DNA binding sites in CAL51 cells. A total of 7,278 peaks accounting for 4,321 unique hits were identied, this number being consistent with previous analysis performed on different cell types9,27. A very limited amount of YAP-binding sites were promoters, most of them being intergenic or intronic (Fig. 3a,b). The list of all targets identied is reported in Supplementary Data 2, whereas the methodology adopted and quality metrics parameters are reported in Supplementary Fig. 4. As YAP does not possess any ability to bind DNA per se, we analysed the transcription factor binding motifs known to target the sites identied by ChIP-Seq. The analysis yielded a number of potential partners involved in YAP co-transcriptional activity (Supplementary Data 3), with some of them being well-known to interact with YAP (TEAD1, p73, RUNX1, RUNX2 and PPARG), while others being newly related to YAP1, such as REST, RREB1 and ESR1 (Fig. 3c). Molecular function annotation of the identied hits claried most of the genes targeted by YAP are involved in cell development, cell growth and proliferation and in the acquisition of cell morphology, roles which have been previously associated to YAP activity (Fig. 3d)27. More interestingly, among the functions revealed by the analysis, those related to the formation and number of FAs displayed high scores (P 1.85 10 4 and P 1.88 10 3, respectively, n 3,

Fishers exact test), while also categories correlated to cytoskeleton organization and formation were represented (P 6.15 10 4) (Supplementary Data 4). The ability of YAP

to target FA genes was corroborated by comparing our original data with the results published by an independent research group in MDA-MB-231 breast cancer line27 (P 1.3 10 4; Easy

score; Supplementary Fig. 5a).

YAP controls FA integrity downstream of RhoA GTPase pathway. Given the evidence that YAP can directly target FA-related genes, we generated CAL51 stable mutant cell lines lacking YAP by CRISPR/Cas9 technology. The characterization of the cell lines is shown in Supplementary Fig. 5be. As predicted by ChIP-Seq analysis, YAP mutant clones underwent a switch in the expression of genes involved in FA formation and cellECM interaction (Fig. 4a). Remarkably, the expression of known targets of YAP like CTGF was consistently reduced in mutant cells. More importantly, the expression prole of integrin subunits was found altered in YAP-depleted cells. ITGA1, ITGA4 and ITGAV, which were also identied as direct targets of YAP in ChIP-Seq analysis, were consistently altered in YAP mutant cells (Fig. 4b). In addition, changes in ITGA5, ITGB3 and ITGB1 were conrmed at the protein level (Fig. 4c) and suggested a higher degree of complexity other than the direct gene targeting by YAP. These changes were accompanied by a clear reduction in key FA structural proteins vinculin and zyxin, and the displacement of talin (Fig. 4d,e and Supplementary Fig. 6a). Similar results were obtained in HEK293 YAP-depleted cells (Supplementary Fig. 6b).

As a result of the inability of developing FA spikes, cells lacking YAP encountered a dramatic shift in cell morphology with colonies failing to spread and acquire their regular shape (Supplementary Fig. 6c). Cell shape is determined by the tight interplay between FAs and cytoskeleton, with docking proteins Ezrin/Radixin/Moesin (ERM) and Enabled/vasodilator-stimulated phosphoprotein (VASP) being responsible for the stabilization of the interphase21,2831. Mutant cells displayed a signicant reduction in ERM expression and its phosphorylated form together with decreased phosphorylation on S157 site of VASP, responsible for its localization to the FAs. No change in VASP phosphorylation at S239 site was detected, a modication needed for cytoskeleton stabilization29. As expected for cells lacking FA spikes, a marked reduction in FAK autophosphorylation on Tyr397 site was found in YAP mutant cells, indicating a lack of integrin heterodimerization, a step needed for FA assembly32 (Fig. 4e).

The integrity of the interphase between FAs and cytoskeleton is required for the cell to establish a functional connection with the ECM during morphogenesis and migration33. As expected, mutant cells lost the ability to spread over the ECM, as measured by cell contact area (AWT 873.1331.0 mm2

versus ACRIPR/CAS9:YAP1: 254.4153.8 mm2), became softer

(EWT 1320.62409.78 Pa versus A

CRIPR/CAS9:YAP1:

405.28143.29 Pa) (n 24, 3 technical replicates, *Po0.001

and **Po0.0001, Welchs t-test) and failed to migrate in two-dimensional and three-dimensional conditions (two technical replicates, Dunns test, ****Po0.001 and **Po0.01)

(Fig. 4f,g). This effect was accompanied in YAP-decient cells by the modulation of genes encoding for proteins involved in the degradation of ECM (Fig. 4g and Supplementary Movie 4).

Rho/ROCK axis is considered to be the main mechanosensing pathway perceiving ECM mechanics34,35, being able to control not only FA formation36 but also actin polymerization dynamics by keeping colin severing protein in the phosphorylated, inactive form30. As Rho/ROCK have been found to control YAP shuttling in different reports37 and no impairment in RhoA, ROCK1 and ROCK2 gene expression was induced by YAP-depleted cells, we next asked whether these effects are connected.

YAP-decient cells displayed a rearrangement in F-actin cytoskeleton in a similar manner to what we described in shYAP AD-MSCs (Fig. 4h), but no evident changes in the tubulin network (Supplementary Fig. 6d). Western blot analysis of the ratio between F- and G-actin claried no impairment in actin bundle polymerization occurred in cells lacking YAP (Fig. 4i). The cells had unaltered ROCK2 levels and accumulated the

NATURE COMMUNICATIONS | 8:15321 | DOI: 10.1038/ncomms15321 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15321

a b

Sample1 Sample 2 Sample3 Sample input

Non-coding

0%

Promoter

1%

TTS

1% 3UTR

1%

5UTR

0% EXON

1%

3 UTR

83

3

Intergenic

49%

Intron

47%

5 UTR Exon

Intron Intergenic

Non-coding

Promoter TTS

66 3,443 3,564

24

28

67

5 kb

5 kb

d Cellular development Cellular growth and proliferation

Cell morphology Cellular assembly and organization

Cellular function and maintenance

Cellular movement Cell-to-cell signalling and interaction

Cell death and survival

Cell cycle Small molecule biochemistry

Cell signalling Post-translational modication

Free radical scavenging Molecular transport

Protein synthesis DNA replication, recombination and repair

Lipid metabiolism Carbohydrate metabolism

Cellular compromise Drug metabolism

0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5

c

log10 (P value)

Figure 3 | YAP directly targets FA genes in CAL51 cells. (a) Heatmap representing the consistent enrichment of YAP binding sites in three independent preparations (A1, A2, A3) of CAL51 cells. (b) Graph: location annotation of YAP DNA binding sites as obtained by ChIP-Seek web tool. Table: location annotation of YAP DNA binding sites reported as absolute values. (c) word cloud representation of the most signicantly represented transcription factors known to bind the sequences identied as YAP targets. Font size is inversely correlated to log10(P value). (d) Table: Gene functional annotation of

the 7,278 hits identied as YAP targets on immunoprecipitated chromatin with their respective P-value.

phosphorylated inactive form of actin severing protein colin (Fig. 4j). These results claried that the lack of FA structural and docking proteins prevented the anchorage of cell cytoskeleton to the membrane, and that Rho/ROCK role in actin dynamics was independent of YAP activity.

On the other hand, we asked whether YAP nuclear shuttling was required for Rho/ROCK acknowledged control over FA formation. The treatment with pharmacological inhibitors interfering with Rho/ROCK pathway at different levels, (Y27362, LIMKi3 and cytochalasin D) caused YAP shuttling to the cytoplasm, triggered cytoskeleton disassembly and the disappearance of FAs (Supplementary Fig. 7a,b). On the contrary, the overexpression of RhoA-Q63L constitutively active mutant increased YAP nuclear localization (Fig. 4k) and resulted in a marked increase in FA signal together with the formation of stress bres. Rho hyperactive mutant failed to induce the formation of FAs in YAP-deprived lines (Fig. 4l). These data claried that YAP was necessary downstream of Rho/ROCK to govern FA formation and cellmatrix interaction.

Cell mechanics is determined by YAP-induced FA assembly. Given the plethora of DNA targets hit by YAP and connected to FA assembly, we asked what are the key genes targeted by YAP that trigger the formation of FAs, thus determining cell spreading and mechanics. Therefore, we complemented YAP mutant cells with single components of the FA complex with a structural function that were depleted or mismatched in mutant cells

(vinculin, zyxin, talin, aV-, b1- and b3-integrins). As a control for proteins that were not affected, we also included a9 integrin. As expected, although YAP transduction in mutant cells restored some degree of vinculin and zyxin expression, single FA proteins such as vinculin, zyxin, talin, a9-, aV-, b1- and b3-integrin failed to induce any change in cell ability to interact with the ECM, as measured by cell surface area and cell morphology (n 26,

3 technical replicates, Dunns test, ****Po0.0001). (Fig. 5a and Supplementary Fig. 8). In previous reports, aVb1 and aVb3 integrinstwo major subsets of integrins displaced in YAP-depleted cellswere associated to cell mechanics, adhesion and migration38,39. Therefore, we co-transfected mutant cells with aVb1 and aVb3-integrins. No signicant increase in mutant cell surface area was induced by aVb1, whereas a signicant and consistent restoration was obtained with aVb3 transfection (AaVb3 751.5.1316.8 mm2 versus A

CRIPR/CAS9:YAP1:

346.6161.1 mm2), similar to the values obtained by re-expressing S127A YAP (AS127A 903.8404.8 mm2) constitutively active

mutant and wild type cells (AWT 823.1355.2 mm2; Fig. 5b,c)

(n 3 independent experiments, Dunns test, ****Po0.0001).

As a result, YAP-defective cells transfected with aVb3 integrins developed a higher surface energy (gaVb3 : (3.71.3) 10 3 mm 2

versus gCRIPR/CAS9:YAP1: (4.21.6) 10 4 J mm 2, n 3 inde

pendent experiments, GamesHowell test *Po0.05 and **Po0.01) and became signicantly stiffer than the mutant cells (gaVb3: 952.8248.5 Pa versus gCRIPR/CAS9:YAP1: 533.980.9 Pa, n 3

independent experiments, GamesHowell test *Po0.05 and ***Po0.001) when measured in AFM.

6 NATURE COMMUNICATIONS | 8:15321 | DOI: 10.1038/ncomms15321 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15321 ARTICLE

DiscussionThe search for the processes by which YAP/TAZ exert their acknowledged mechanosensing activity recently led to a model in which two different upstream pathways partially overlap with Hippo pathway to control their nuclear shuttling and thus regulate their role as transcriptional co-activators. Small GTPase Rho indirectly controls YAP/TAZ nuclear localization by promoting the formation of actin bundles and stress bres in

response to the spreading of the cell over ECM12,40. Moreover, YAP/TAZ upstream kinase LATS is targeted by integrin-FAK signalling through Src15. These results depicted a role for Hippo effectors as a hub for different signalling pathways sensing ECM composition and mechanics through cytoskeleton stability. Here we demonstrate that: (1) YAP nuclear localization is controlled by cell area through Rho/ROCK activation and independent of FA formation; and (2) YAP transcriptionally controls FA formation

a b c

Integrin 4

Integrin 5 Integrin V

Integrin 1 Integrin 3

Integrin 4 Integrin 5

Vinculin Zyxin

Talin GAPDH

100 kDa 150 kDa

100 kDa 100 kDa 100 kDa 250 kDa

75 kDa

75 kDa

250 kDa

37 kDa

Fold regulation

150 kDa

d

DAPI vinculin Detail

e f

75 kDa 75 kDa

50 kDa

p-ERM

ERM VASP (157)

p-VASP (239) VASP

FAK p-FAK (397)

GAPDH

Contact area (m)

2,500 2,000 1,500 1,000

500

0

Youngs modulus (Pa)

2,500 2,000 1,500 1,000

500

0

*** ****

CRISPR/Cas9:Yap1

50 kDa

50 kDa

100 kDa

100 kDa

37 kDa

g h i

DAPI F-actin Detail

DAPI EGFP-RhoA-Q63L vinculin F-actin

G-actin

F-actin

G-actin

F-actin

G-actin

F-actin

CAL51 WT

c3c2a

2D migration 3D migration

CRISPR/Cas9:Yap1 CAL51 WT

** ****

G-actin/F-actin

ratio

50 m

CRISPR/Cas9:Yap1

CTR c2a c3

j

ROCK2 p-cofilin cofilin GAPDH

150 kDa

20 kDa

20 kDa

37 kDa

k l

DAPI YAP

CRISPR/Cas9:Yap1CAL51 WT

EGFP-RHoA-Q63L

F-actin

NATURE COMMUNICATIONS | 8:15321 | DOI: 10.1038/ncomms15321 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15321

YAP displayed a specic DNA-binding signature for genes encoding for proteins involved in cytoskeleton stabilization (CAPZA1, CAPZA2, ABRA, ACTRT1, CKAP4, CKAP4 and TLN2) and connection to the membrane (PCDH15, CTNNA2, CTNNA3, CTNNB1, CTNND2, FAT3, PCDH7, PCDH17, PCDH10, PCDH20, CDH8, CDH5, CDH2, CDH20 and CDH26).Among YAP-targeted genes, those related to dystrophin/ sarcoglycan complex (SGCD, SGCG, ZSC1, SNTB1 and IMPG2) have also been identied in the fetal heart13 as YAP targets and directly related to the stabilization of the cytoskeleton.

Interestingly, YAP depletion resulted in the rearrangement of integrin expression pattern. Integrin expression prole modulates cell migration in response to ECM cues, with aVb1 being responsible for random cell migration, whereas aVb3 induce persistent directional cell migration43. The latter has also been shown to increase cell stiffness and is currently being targeted by etaracizumab monoclonal antibody in clinical studies aimed at hindering cancer propagation44. Although not being detected by our analysis as a direct target of YAP, b3 integrin is among the genes displaced by YAP knockout and, when re-expressed in combination with its partner aV, restored cell mechanical properties and the ability to spread. This evidence reinforces the indication of YAP as a target for cancer molecular therapy45.

The evidence that YAP-depleted cells encountered a signicant reduction in their stiffness qualied YAP as a pivotal component of the mechanosensing apparatus and depicted a model in which YAP co-transcriptional activity is required for the synthesis of docking and stabilizing proteins allowing cell tension generation.Not surprisingly, the absence of adhesion complexes in YAP mutant cells resulted in their impaired ability to get a grip on the surrounding ECM. Indeed, cells lacking YAP displayed lower adhesion energy and reduced contact area when interacting with elements of the ECM. As a proof-of-concept, mutant cells failed to migrate and invade the surrounding matrix, when challenged in two-dimensional and three-dimensional assays.

YAP has been historically described to exert overlapping/ redundant functions with its paralog protein TAZ37. Although TAZ depletion performed similarly in terms of cell mechanics in our cell models, it did not affect FA dynamics, this being an exclusive feature of YAP. Moreover, the mild observed upregulation of TAZ in YAP-silenced cells was incapable to counteract FA loss46. The apparent discrepancy between YAP and TAZ activities could be explained by the slight but consistent differences in the transcriptional and functional signatures of YAP and TAZ reported here and by others46,27. Solving this puzzle deserves further investigations.

Figure 4 | YAP regulates focal adhesion/cytoskeleton integrity downstream of RhoA. (a) Heatmap representing the RNA expression levels by normalized DCt values obtained in two qRTPCR independent replicates of YAP mutant clone C3 as compared with CAL51 WT. (b) Graph: barplot representing mean values.d. of up- and downregulated genes with a fold change higher than 2.0 in YAP mutant clone C3 as compared to CAL51 WT.

(c) Western blot analysis of the indicated focal adhesion (FA) proteins; n 3. (d) Representative confocal images of WT CAL51 and YAP mutant clones

C2a and C3, showing vinculin staining (green). (e) Protein analysis of the indicated cytoskeleton- and FA-associated proteins in WT CAL51 and YAP mutant clones; n 3. (f) Graphs: quantication of cell contact area and Youngs Modulus in WT CAL51 and YAP-defective clone C3. Values are shown as

medianmin/max (n 24, 3 technical replicates, Welchs t-test, ***Po0.001 and **** Po0.0001,). (g) Box: Representative brighteld image of WT

CAL51 and C3 clone migrated through ECM-coated transwell membrane (8 mm) and stained with crystal violet after 24 h, quantication is shown in the upper graph. The data represent the mean values.d. (12 random elds, 2 technical replicates, **Po0.01 and ***Po0.001, KruskalWallis test followed by post hoc Dunns test). Bottom: qRTPCR array analysis of genes involved in ECM remodelling in YAP mutant cells (C3) as compared with the control.

The threshold was set at 2.0. Right: Representative brighteld image of CAL51 WT and YAP mutant cells C3 grown in three-dimensional (3D) Matrigel for 120 h. (h) Representative confocal images showing the arrangement of F-actin in WT Cal51 and mutant clones C2a and C3 stained with Alexa Fluor 546 Phalloidin. (i) Representative western blot quantication of G-actin/F-actin ratio in WT CAL51 and YAP mutant clones. (j) Representative western blot analysis of the indicated proteins; n 3. (k) Representative confocal image of WT CAL51 cell transfected with EGFP-RhoA-Q63L and stained with anti-YAP

antibody (red), Alexa Fluor 647 Phalloidin (white). Image analysis shows the nuclear/cytoplasmic distribution of YAP protein in non-transfected and transfected cells. (l) Confocal images of WT CAL51 and mutant clone C3 transfected with EGFP-RhoA-Q63L and stained with anti-vinculin antibody (red), Alexa Fluor 647 Phalloidin (white).

and cytoskeleton stability, in turn determining cell mechanics and the degree of cell adhesion to the ECM.

Through a single-cell tool based on micropatterns able to control FA formation independently of cell size41, we claried that YAP nuclear presence in AD-MSCs does not depend on the formation of FAs, but it is rather controlled by the area of the cell itself. Cell spreading over ECM directly controls RHOA gene expression and activation and this event is crucial to stabilize cell cytoskeleton by inhibiting the F-actin severing protein colin through phosphorylation by LIMK35. As conrmed in previous studies12, YAP acts downstream of Rho/ROCK and is sensitive to F-actin bundle polymerization. Conversely, the integrity of cytoskeleton actin bundles was not affected per se by YAP depletion, as shown by the unaltered ratio between F- and G-actin, and the sustained phosphorylation of colin severing protein, thus indicating that Rho/ROCK role in actin dynamics was independent of YAP activity.

The activation of Rho GTPase pathway is also required to foster the maturation of FAs in cells in contact with the ECM36. The anchorage of cell cytoskeleton to the FAs is of paramount importance to control cell shape and mechanics, with docking proteins such as vinculin, talin, zyxin, ERM and VASP being primarily involved in coupling cytoskeleton to integrins21,2831. The lack of YAP in the cell caused the complete absence of the typical FA spikes highlighting the inability of YAP mutant cells to interact with the surrounding ECM. As a control, cells over-expressing the constitutively active and nuclear form of YAP (S127A YAP) showed an increased FA formation. YAP-depleted cells showed a reduced expression of FA stabilizing proteins ERM, vinculin and zyxin, and showed altered localization of actin docking protein VASP. The decreased expression and phosphorylation of ERM, a set of proteins providing structural linkage between transmembrane components and actin laments42 together with the decreased phosphorylation on S157 site of VASPa protein entitled to stabilize the interaction among FA docking proteins and actin cytoskeleton29accounted for the failure in promoting the anchorage between the cytoskeleton and FAs.

Our ChIP-Seq analysis claries that the stability of FA-cytoskeleton interphase is directly under the control of YAP: YAP co-transcriptional function, as determined by cell area through Rho/ROCK-mediated cytoskeleton assembly, directly targets DNA binding sites responsible for the production of a number of proteins involved in FA and cytoskeleton assembly and crucial to cell-ECM interaction. These results, obtained in breast cancer CAL51 cells, are conrmed in other cell types by an unbiased meta-analysis performed on ChIP-Seq data previously obtained by other research groups9,27.

8 NATURE COMMUNICATIONS | 8:15321 | DOI: 10.1038/ncomms15321 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15321 ARTICLE

a

b

****

1,400

1,200

1,000

800

600

400

200

DAPI

F-actin

CRISPR/CAS9:YAP1

cell area (m)

V integrin

GFP

0 YAP vinculin

zyxin

V integrin

9 integrin

1 integrin

3 integrin

V integrin

3 integrin

1 integrin

c

d

DAPI vinculin

** *

** *

*

**** ****

***

*

Cell area (m)

Surface energy (Jm)

Young's modulus (Pa)

V integrin 3 integrin

e

FA

genes

TEAD

YAP

TEAD

CFL1 CFL1

Cell spreading

CFL1 CFL1

Cell spreading

Actin stabilization

LIMK

LIMK

ZYX

ZYX

IT

ROCK

ROCK

YAP

VASP

VASP

V

V

V

V

V

RHO

RHO

TLN VCL

VCL

V

V

TLN

VCL

E

ERM

E

ERM

IT IT

IT

Figure 5 | YAP controls cell biophysical properties and the interaction with ECM through aVb3 integrin. (a) Graph: YAP mutant cell contact area. YAP-depleted cells were transfected as to express the indicated proteins and the contact area of transfected cells was evaluated and compared to the GFP-transfected cells. (n 26, 3 technical replicates, ****Po0.0001, as calculated by KruskalWallis test followed by post hoc Dunns test). Barplot

represents mean valuesd. (b) Confocal representative images of YAP-deprived cells transfected with either aVb1 or aVb3 integrins. Cell F-actin was stained with Alexa Fluor 647 Phalloidin (white) and nuclei counterstained with 40,6-diamidino-2-phenylindole (DAPI) (blue). (c) Confocal images of CAL51 mutant cells transfected with either aV (green), b3 (red) or co-transfected with aVb3 integrins (co-stained by red and green). Double transfected cell is highlighted by a dashed line. Cells were decorated with vinculin (white) and nuclei counterstained with DAPI (blue). (d) Graphs: comparison of the indicated biophysical properties of WT CAL51 (grey) and YAP mutant clone C3 (pale grey), transfected or not to overexpress the reported proteins. Cell area: n 10, 3 technical replicates, ****Po0.0001, as calculated by Kruskal-Wallis test followed by post hoc Dunns test. The barplot represents mean

valuesd. The surface energy is given in logarithmic scale. Youngs Modulus, surface energy: n 10, 3 technical replicates, *Po0.05, **Po0.01 and

***Po0.001 as calculated by Welchs ANOVA test followed by post hoc Games-Howell test. Values are shown as medianmin/max. (e) Model proposing YAP control over focal adhesion (FA) assembly. Left: Cell area controls RhoA activity to phosphorylate colin (CFL1) and maintain the integrity of actin ber cytoskeleton, which triggers YAP nuclear shuttling. In the nucleus, YAP promotes the transcription of genes encoding for proteins participating in FA assembly. Right: when YAP is lost, Integrin subunit remodelling occurs, docking proteins like vinculin, zyxin, ERM and VASP are removed from the interphase between FAs and F-actin, and the anchorage of cytoskeleton to the membrane is disrupted. All error bars are s.d.

Altogether, the data provided emphasize the role of YAP in controlling cell mechanics and describe a model in which cell spreading over the surrounding ECM modulates YAP transcriptional activity through Rho/ROCK axis to determine FA-cytoskeleton remodelling (Fig. 5e). When taken together with recent reports describing YAP as a regulator of ECM mechanics8, our ndings point at YAP as the master regulator of cellECM interaction.

Methods

Cell culture and differentiation. ASC52telo, hTERT immortalized AD-MSC cells (ATCC SCRC-4000) were purchased from American Type Culture Collection

(ATCC, Manassas, USA). HEK293 cells were kindly provided by Dr V. Pekarik (Department of Physiology, Masaryk University, Brno, Czech Republic). CAL51 cells were a gift of Dr L. Krej (Department of Biology, Masaryk University). Cells were cultured in DMEM medium 4.5 g l 1 Glucose (DMEM high Glucose, Lonza,

Basel, Switzerland) supplemented with 10% fetal bovine serum, 2 mM L-glutamine and 100 U ml 1 penicillin/streptomycin. YAP and TAZ knockdown was performed in AD-MSCs by using a lentiviral vector short hairpin RNA, (shYAP, shTAZ; Santa Cruz Biotechnology, Dallas, USA), according to the manufacturers instructions. Control cells were transfected with Control short hairpin RNA Lentiviral Particles (Santa Cruz Biotechnology). Briey, 50% conuent AD-MSCs were cultured on 24-well plates and incubated with retroviral particles (multiplicity of infection: 1) in the presence of 5 mg ml 1 polybrene (Santa Cruz Biotechnology).

One day after infection, the medium was replaced with fresh medium and cells were allowed to recover for 48 h. Infected cells were selected by adding 4 mg ml 1

NATURE COMMUNICATIONS | 8:15321 | DOI: 10.1038/ncomms15321 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15321

puromycin (Santa Cruz Biotechnology) to the culture medium until the experiments were performed. AD-MSC adipogenic differentiation was assessed by culturing cells in complete differentiation medium (Thermo Fisher Scientic, Waltham, USA), according to the manufacturers instructions. Three days after induction, lipid droplets were stained using AdipoRed, counterstained with 40,6-diamidino-2-phenylindole (Thermo Fisher Scientic) and visualized by confocal microscopy. Nine independent elds were acquired at low magnication and the ratio of cells containing lipid droplets was calculated.

For the inhibition of Rho/ROCK pathway, cell were treated with Y27632(10 mM, Sigma-Aldrich, St. Louis, USA), LIMKi3 (20 mM, Tocris Bioscience, Bristol,

UK), cytochalasin D (1 mM, Sigma-Aldrich), or the corresponding amount of DMSO diluted in complete media for 24 h and then xed and stained as described afterwards.

Cell transfection. Cell transfection and co-transfection were performed by using Lipofectamine 3000 (Thermo Fisher Scientic). The plasmids pEGFP-C3-hYAP1 (17843, gift from Marius Sudol), 2 FLAGhYAP1 (17791, gift from Marius Sudol),

2 FLAGhYAP1-S127A (17790, gift from Marius Sudol), pLX304-YAP1_PDZ

(59147, gift from William Hahn), 3 Flag pCMV5-TOPO TAZ wild type (WT)

(24809, gift from Jeff Wrana), 3 Flag pCMV5-TOPO TAZ (S89A) (24815,

gift from Jeff Wrana), 3 Flag pCMV5-TOPO TAZ (q304) (24810, gift from Jeff

Wrana), pcDNA3-EGFP (13031, gift from Doug Golenbock), pcDNA3-EGFPRhoA-Q63L (12968, gift from Gary Bokoch), EGFP-N3 Integrin alpha9 (Plasmid 13600, gift from Dean Sheppard), mEos2-Alpha-V-Integrin-25 (Plasmid 57345, gift from Michael Davidson and Catherine Galbraith), mCherry-Integrin-Beta1-N-18 (Plasmid 55064, gift from Michael Davidson), b3-integrin-YFP (Plasmid 26653, gift from Jonathan Jones) and 8 GTIIC-lux (Plasmid 34615, gift from Stefano

Piccolo) were obtained from Addgene. AD-MSCs and CAL51 were seeded into eight-well chamber slides. After 24 h, a preincubated mixture containing 20 ml of

Opti-MEM, 0.2 ng of DNA, 0.6 ml of Lipofectamine 3000 and 0.4 ml P3000 reagent was added in each well. After 3 h, the medium was replaced and cells allowed to recover for 48 h before analyses.

For the AFM measurements, cells were transfected, sorted for the corresponding uorescent channel and positive cells were seeded on FN-coated dishes and measured after 48 h similar to that described.

Luciferase assays analysis was quantied with Berthold CENTRO LB 960

Microplate Luminometer (Berthold Technologies GmbH, Austria).

Generation of YAP mutant CAL51 lines. YAP-decient Cal51 lines were produced by CRISPR/Cas9 technology. Guiding RNA was designed to hit exon 1 of YAP1 gene, which is common to all nine YAP1 splicing variants (Supplementary Fig. 5b). Two sets of complementary single-stranded DNA oligonucleotides (YAP1_R1: 50-CACCgtgcacgatctgatgcc-30, YAP1_R2: 50-AAACccgggcatcagatcgtg cac-30, YAP1_F1: 50-CACCGcatcagatcgtgcacgt-30, YAP1_F2 50-AAACcggacgtgcacg atctgatgC-30) were cloned into pSpCas9(BB)-2A-GFP (PX458) (Addhene 48138, gift from Feng Zhang) and transfected into Cal51 cells using FuGENE HD (Promega Corporation, Wisconsin, USA) transfection reagent according to manufacturers protocol. The next day, green uorescent protein (GFP)-positive cells were FACS-sorted (MoFlo Astrios, Beckman Coulter, California, USA) as single cells into 96-well plate and clonally propagated. Genomic DNA was sequenced from both sides to map the deletion size (sequencing primers: 50-gattgg acccatcgtttgcg-30, 50-gtcaagggagttggagggaaa-30, 50-gaagaaggagtcgggcagctt-30, 50-gag tggacgactccagttcc-30).

Micropatterned substrates preparation. Fibronectin-coated or activated micropatterned slides with different area, shape or pattern (ref: 109501018; 109500018, or custom containing only 10,000 mm2 or 1,024 mm2 squares) were purchased from CYTOO (CYTOO, Grenoble, France). Laminin coating was performed by incubating activated slides with 20 mg ml 1 laminin (Sigma-Aldrich)

at room temperature (RT) for 2 h. Cadherin coating was performed using a modication of the method described by Czndr et al.47. Briey, substrates were treated with 40 mg ml 1 poly-D-lysine (Merck Millipore, Billerica, USA) for 2 h, incubated with goat anti-human Fc antibody (Jackson Immunoresearch, West

Grove, USA) in 0.2 M, pH 8.5 boric acid at RT for 5 h. Substrates were washed with borate buffer and incubated O/N at 4 C with Recombinant Human E- or N-Cadherin Fc Chimeras (R&D Systems, Minneapolis, USA) in boric acid. Substrates were then washed and equilibrated in DMEM containing 10% fetal bovine serum for 1 h at 37 C. AD-MSCs were seeded at 2 104 cells per cm2

directly on the slides and cultured in standard culture conditions.

To investigate the effect of substrate stiffness on YAP and vinculin expression, cells were seeded onto Fibronectin-coated m-Dish 35 mm, high elastically supported surface 28 and 1.5 kPa (iBIDI, Munich, Germany).

Immunostaining and image analysis. Immunouorescence staining was performed as previously described48. Briey, cells were xed in 4% paraformaldehyde in PBS for 15 min at RT and permeabilized with 0.1% Triton X-100 for 2 min. After incubation with primary antibodies, cells were incubated with the appropriate Alexa uorochrome-conjugated secondary antibodies. F-actin was stained with Alexa Fluor 546 or 647-conjugated phalloidin. Nuclei were counterstained with

40,60-diamidino-2-phenylindole. Samples were embedded in ProLong Gold antifade reagent (Thermo Fisher Scientic) and visualized with Zeiss LSM 780 confocal microscope with 40 (1.3 numerical aperture) and 63 (1.4 numerical aperture)

oil-immersion objective lenses. Z-stacks were acquired with the optimal interval suggested by the software, followed by the application of maximum intensity algorithm.

YAP nucleus/cytoplasm ratio was calculated using an ad hoc developed ImageJ plug-in with the following formula:

PInuc =Anuc

PIcyto =Acyto

where PInuc and PIcyto represent the sum of the intensity values for the pixels in the nuclear and cytoplasmic region respectively, and Anuc and Acyto the area of the corresponding regions.

Focal adhesion quantication was performed by ImageJ software as follows. All the images were acquired at the same magnication and resolution as described in (ref. 49). SUBTRACT BACKGROUND was applied to the channel corresponding to vinculin staining with a SLIDING PARABOLOID option and ROLLING BALL radius of 25 pixels. The images were enhanced by running CLAHE plug in (lock size 19, histogram bins 256, maximum slope 3, no mask and fast) followed

by automatic BRIGHTNESS/CONTRAST and ENHANCE CONTRAST (saturated 0.35). Images were nally binarized using automatic THRESHOLD

command (default settings) and particles analysed (size 0.3015;

circularity 0.000.99).

Primary antibodies used were as follows: rabbit anti-YAP (4912, Cell Signaling Technology, Danvers, USA), mouse anti-YAP/TAZ (sc-101199), rabbit anti-TAZ (sc-48805) (Santa Cruz Biotechnology), rabbit anti-Fibronectin (F3648), mouse anti-Talin (SAB4200041), mouse anti-Vinculin (V9131), mouse anti-Zyxin (Z0377), mouse anti-b-Tubulin (T5076) (Sigma-Aldrich) and mouse anti-Vinculin (ab18058, Abcam, Cambridge, UK). The full list of antibodies and dilutions can be found in Supplementary Data 5.

Isolation of RNA and real-time quantitative PCR analysis. Total RNA was isolated using RNeasy Mini Kit (Qiagen, Hilden, Germany), according to the manufacturers instructions. Complementary DNA was synthesized using the RT2 First Strand Kit (SABiosciences, Frederick, USA). The expression prole of genes involved in different pathways was analysed by the following RT2 Proler PCR Arrays (Qiagen): MSCs (PAHS-082Z), Hippo signalling pathway (PAHS-172Z), human adherens junctions (PAHS-146Z), ECM and adhesion molecules (PAHS-013Z) and FAs (PAHS-145Z). RTPCR was carried out on the LightCycler 480 Real-Time PCR System (Roche, Basel, Switzerland) using the following cycling parameters: 1 cycle at 95 C for 10 min; 45 cycles at 95 C for 15 s and 60 C for 1 min. The internal panel of housekeeping genes set from the manufacturer was used for normalization of expression levels of individual genes and PCR-array data analysed by online resources from the manufacturers website (http://www.sabiosciences.com/pcrarraydataanalysis. php

Web End =http://www.sa http://www.sabiosciences.com/pcrarraydataanalysis. php

Web End =biosciences.com/pcrarraydataanalysis. php ) and statistical R-project software (http://www.r-project.org

Web End =http://www.r-project.org). Genes showing high coefcient of variation among replicas, as well as very low expression in all the experiments (35oCto40) were discarded from the analysis. The results include the heatmaps of quantication cycles (Ct) and the graphs of the mean and s.e.m. values of the fold regulation values obtained by analysing independently two samples per each experimental condition.

Western blotting. Cells were lysed in RIPA buffer (Merck Millipore) with protease and phosphatase inhibitor cocktails (1%, both Sigma-Aldrich) on ice and then centrifuged at 13.000 g for 10 min at 4 C. Supernatants were stored at

80 C. Protein concentration was determined using the BCA method. Forty micrograms of proteins for each sample were loaded in on 10% polyacrylamide gels (Bio-Rad, Marnes-la-Coquette, France), pre-run at 40 V and then at 100 V. The fractionation into soluble, nuclear and cytoskeletal compartments was performed using a ProteoExtract Cytoskeleton Enrichment and Isolation Kit (Merck Milli-pore), according to the manufacturers manual. G/F-actin ratio was quantied using the G-Actin/F-Actin In Vivo Assay Biochem Kit (Cytoskeleton, Denver, USA), according to the product specication. Proteins were transferred to a polyvinylidene diuoride membrane using the Trans-Blot Turbo transfer system (Bio-Rad). Membranes were blocked with 5% BSA in TBST, incubated with diluted primary antibody in 5% BSA in TBST at 4 C O/N and then probed with the proper secondary HRP-linked antibody (Sigma-Aldrich) at RT for 1 h. ChemiDoc imaging system (Bio-Rad) was used to detect chemiluminescence. Band intensities were quantied using Bio-Rad Image Lab software.

Primary antibodies used were as follows: rabbit anti-YAP (4912), rabbit anti-YAP/TAZ (8418), rabbit anti-Integrin a4 (8440), rabbit anti-Integrin a5 (4705), rabbit anti-Integrin aV (4711), rabbit anti-Integrin b1 (9699), rabbit anti-Integrin b3 (13166), rabbit anti-Integrin b4 (4707), rabbit anti-Integrin b-5 (3629), rabbit anti-Colin (3313), rabbit anti-Phospho-Colin (Ser3) (5175), rabbit anti-Ezrin/Radixin/Moesin (3142), rabbit anti-Phospho-Ezrin (Thr567)/Radixin (Thr564)/Moesin (Thr558) (3726), rabbit anti-VASP (3132), rabbit anti-Phospho-VASP (Ser157) (3111), rabbit anti- Phospho-VASP (Ser239) (3114), rabbit anti-FAK (13009), rabbit anti-Phospho-FAK (Tyr397) (8556) (Cell Signaling

10 NATURE COMMUNICATIONS | 8:15321 | DOI: 10.1038/ncomms15321 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15321 ARTICLE

Technology), rabbit anti-TAZ (PA1-46190), mouse anti-GAPDH-HRP conjugate (MA5-15738-HRP) (Thermo Fisher Scientic), mouse anti-Vinculin (ab18058), mouse anti-ROCK-2 (ab56661) (Abcam), mouse anti-Talin (SAB4200041), mouse anti-Zyxin (Z0377), mouse anti-b-Tubulin (T5076) (Sigma-Aldrich) and rabbit anti-Pan-Actin (AAN01, Cytoskeleton). The full list of antibodies and dilutions can be found in Supplementary Data 5. The original blots are shown in Supplementary Fig. 9.

Scanning electron microscopy. Cal51 WT and YAP mutant cells were cultured in bronectin-coated coverslips for two days, xed with a solution of Glutaraldehyde 3% in 100 mM Cacodylate buffer and dehydrated in a series of increasing ethanol concentrations. Samples were mounted on aluminium stubs, sputter-coated with Palladium (JEOL JFC-1300, Tokyo, Japan) and imaged with a Benchtop Scanning Electron Microscope JEOL JCM-6000 (Tokyo, Japan).

Chromatin immunoprecipitation. CAL51 cells were cultured in complete medium for 48 h. Chromatin was immunoprecipitated from three technical replicates using a ChIP-grade anti-YAP antibody (Cell Signaling Technologies) and following the manufacturers instructions (Pierce Agarose ChIP Kit, Thermo Fisher Scientic). Control ChIP was performed by adding mouse immunoglobulins (IgG). The samples were eluted in 30 ml of eluting buffer and stored at -80 C before the analysis.

Library preparation and sequencing. Library preparation was performed and size distribution of each ChIP DNA sample was assessed by running a 1 ml aliquot on Agilent High Sensitivity DNA chip using an Agilent Technologies 2100 Bioanalyzer (Agilent Technologies). The concentration of each DNA sample was determined by using Quant-IT DNA Assay Kit- High Sensitivity and Qubit Fluorometer(Life Technologies). Ten nanograms of puried ChIP DNA were used as starting material for sequencing libraries preparation. Indexed libraries were prepared with TruSeq ChIP Sample Prep Kit (Illumina Inc.). Libraries were sequenced (single read, 1x50 cycles) at a concentration of 10 pM/lane on HiSeq 2500 (Illumina Inc.).

ChIP-Seq data analysis. Data analysis was performed by Genomix4Life S.r.l (Baronissi, Italy). The raw sequence les generated (.fastq) underwent quality control analysis using FASTQC (http://www.bioinformatics.babraham.ac.uk

Web End =http://www.bioinformatics.babraham.ac.uk). Reads were aligned to the human genome (assembly hg19) using bowtie49, allowing up to 1 mismatch and considering uniquely mappable reads. The reads of biological replicates and corresponding input samples were merged for peaks calling as previously described50. ChIP-Seq peaks were identied and analysed using HOMER51 (-F: 2.0, -L: 2.0 and -C: 1.0) with false discovery rateo0.01.

The assignment of YAP peaks to target genes was obtained by the web tool ChIPSeek52. Through this step, it was possible to assign peaks to transcription start site (by default dened from 1 kb to 100 bp), transcription termination site

(by default dened from 100 bp to 1 kb), Exon (Coding), 50-untranslated

region (UTR) Exon, 30 UTR Exon, Intronic or Intergenic. As some annotation overlap, the following order of priority was chosen for the assignment:

Transcription start site (by default dened from 1 kb to 100 bp) Transcription termination site (by default dened from 100 bp to 1 kb) CDS exons 50 UTR exons 30 UTR exons **CpG islands **Repeats Introns Intergenic

Over-represented sequence motifs were dened according to motif descriptors of JASPAR database and computed using PScan-ChIP53. The following parameters were set:

Organism: Homo Sapiens Assembly: hg19 Background: Mixed Descriptors: Jaspar 2016

Nucleotide best occurrence was calculated by weblogo (http://weblogo.berkeley.edu/

Web End =http:// http://weblogo.berkeley.edu/

Web End =weblogo.berkeley.edu/ ) by running Report Best Occurrences analysis on any given transcription factor.

Functional enrichment analysis was performed using IPA (IPA, QIAGEN, Redwood City, http://www.qiagen.com/ingenuity

Web End =www.qiagen.com/ingenuity ) by using the following parameters:

Reference set: Ingenuity Knowledge Base (Genes Only) Relationship to include: Direct and Indirect Filter: Consider only molecules and/or relationship where (species Human)

AND (condence Experimentally Observed).

An unbiased comparison was performed between the single YAP targets obtained in the CAL51 ChIP-Seq analysis and the ones previously reported for

MDA-MB-231 (ref. 27). The mutual targets enrichment for KEGG biological pathways was examined through the Database for Annotation, Visualization and Integrated Discovery 6.8. An a 0,01 for the EASE score P-value was used as a

threshold.

Single- and collective cell invasion assays. Single-cell invasion was assayed by culturing 30,000 CAL51 WT and YAP mutant cells on ECM pre-coated transwell chambers with 8 mm pores and allowed to migrate for 24 h, according to the manufacturers instructions (ECM508 - Chemotaxis Cell Migration Assay, Merck Millipore). Cells in the bottom of the lters were stained with crystal violet, counted in twelve random elds and averaged for the replicate wells.

Growth factor-reduced Matrigel (BD Biosciences, New Jersey, USA) was diluted to 4 mg ml 1 in serum free-cold cell culture media. 10 ml Matrigel containing 1500 wild type or YAP-mutant CAL51 were transferred to each well of a m-Slide (Ibidi, Martinsried, Germany). Matrigel was allowed to solidify at 37 C for 20 min, before adding complete medium. Medium was changed after 3 days.

Youngs modulus mapping by AFM. Standard bioAFM microscope JPK NanoWizard 3 (JPK, Berlin, Germany) was used to perform force mapping procedure. The scanning-by-probe head (maximal visualization range 100-100-15 mm in X-Y-Z axis) of AFM microscope was placed on inverted optical microscope Olympus IX-81, 10 objective was used to nd proper area covered with cells and

to place cantilever in the proper position for force mapping procedure. Plastic Petri dish (TPP, Trasadingen, Switzerland) with either the distilled water for instrument calibration or with cell culture was placed inside the Petri dish heater (JPK) pre-heated to 37 C for 30 min. Non-coated silicon nitride AFM cantilever Hydra 2R-100N (AppNano, California, USA) equipped with pyramidal silicon tip was used for all the experiments. The probe was calibrated before each experiment as described below. Then the laser reection sum was maximized, followed by centering of the laser detector. The AFM probe was introduced in contact with the surface during a standard process of landing. The sensitivity of the AFM setup was determined as a slope of the force-distance curve (FDC) measured by lifting the cantilever with Z-height of 450 nm, time per curve was 3 s. The sensitivity was found in the range 15.07-15.37 nm V 1, cantilever stiffness was calibrated by measurement of its thermal noise and resulted between 17.34 and 19.19 pN nm 1 for different days of experiments.

The bioAFM setting was identical for all the force mapping procedure. SetPoint value was 1.0 nN (relative to baseline value), time per curve 0.5 s, Z-length 15.0 mm, speed of curve recording 30 mm s 1, the FDCs were recorded with data sample rate of 2 kHz. The force mapping procedure was performed as step-by-step recording of

FDCs in the network of 64 64 points on 100 100 mm covering area of single cell

for AD-MSCs or colonies for CAL51. Force mapping process provides a network of FDC (dependency of tip-sample interaction force on tip height above the surface), so called force maps. The absolute value of Youngs modulus can be determined by tting the FDC by Sneddon equation54:

Fd

2E tana p1 u

2

d2

Where F is the measured force, E is Youngs modulus, n is Poisson ration (0.5 for incompressible materials), d is tip-sample separation (obtained by correction of the cantilever height to its bending) and a is half-angle to face of pyramidal tip (reects the tip geometry). Data processing module of the JPK software was used to process the maps of FDCs in a batch mode. Resulted Youngs modulus maps were exported to be post-processed as described in the following chapter.

AFM data analysis and cell mechanics evaluation. The elasticity maps obtained from measurements were pre-processed in open source SPM analysis software Gwyddion to extract geometrically coherent text tables of elasticity values and setpoint positions. Subsequently, data were processed by a Matlab (Mathworks, Natick, Massachussets, USA) script to remove outliers (for example, elasticity values coming from non-contacting points), select the portion of data inherent to the biological sample, observing the values distribution on histogram and calculate the statistical parameters of interest. Since cell-related elasticity value showed non normal distributions (that is, not passing a Lilliefors test), median and inter-quartile range were rstly computed as statistical parameters of confrontation. A ner extraction of the mean elastic modulus was performed with the Matlab distribution tting tool, where data were found to best t an inverse Gaussian distribution. Cell adhesion energy (g) was calculated by using the JKR theory as described in ref. 55, assuming a conservative volume between the oating cell and the fully adherent one, and a zero load radius, through the following equation:

g

8a3E

27pR2

where a is the equivalent zero-loads contact radius (that is, the radius of a circle covering the same area) of an adherent cell, R is the measured contact of a oating cell and E is the elastic modulus measured as described above through AFM force mapping technique. Cell contact area was measured by ImageJ software region of interest tool on adherent CAL51 cells. The same software tools were used to extract the oating radius of the cell strains analysed. A sample number of 24, based on the

NATURE COMMUNICATIONS | 8:15321 | DOI: 10.1038/ncomms15321 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 11

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15321

sample number selected from the elasticity values, was considered selecting data with the minimum distance from their median.

Statistical analysis. Data are presented as means.d. and were calculated using the software package GraphPad Prism v. 6.0. For single cell analysis, a minimum of ten cells per sample was considered. The number of cells in 12 random elds per sample was evaluated in the migration assay. Statistical analyses were performed using Kruskal-Wallis test followed by post hoc multiple comparisons by Dunns test. P40.01 was not considered statistically signicant. AFM data are shown in boxplot format, considering the medianmin/max values calculated for each selected sample. For AD-MSCs data, which did not pass the Lilliesford normality test, statistical signicance was calculated with the Kruskal-Wallis test followed by post hoc Dunns test, and accepted for a Po0.05. For CAL51 data, passing the

Lilliefors normality test but not homoscedasticity F-test, statistical signicance was assessed by an unpaired two-tailed Welchs t-test, or Welchs analysis of variance test followed by GamesHowell test. Sample sizes were based on previously published experiments, in which statistical differences were identied.

Data availability. ChIP-Seq analysis data were submitted to ArrayExpress database (https://www.ebi.ac.uk/arrayexpress/

Web End =https://www.ebi.ac.uk/arrayexpress/) where they can be accessed by the accession number: E-MTAB-5217. The data that support the ndings of this study are available from the corresponding author upon reasonable request.

References

1. Chicurel, M. E., Singer, R. H., Meyer, C. J. & Ingber, D. E. Integrin binding and mechanical tension induce movement of mRNA and ribosomes to focal adhesions. Nature 392, 730733 (1998).

2. Wozniak, M. A., Modzelewska, K., Kwong, L. & Keely, P. J. Focal adhesion regulation of cell behavior. Biochim. Biophys. Acta 1692, 103119 (2004).3. Fenner, J. et al. Macroscopic stiffness of breast tumors predicts metastasis. Sci. Rep. 4, 5512 (2014).

4. Mouw, J. K. et al. Tissue mechanics modulate microRNA-dependent PTEN expression to regulate malignant progression. Nat. Med. 20, 360367 (2014).

5. Wang, K., Degerny, C., Xu, M. & Yang, X.-J. YAP, TAZ, and Yorkie: a conserved family of signal-responsive transcriptional coregulators in animal development and human disease. Biochem. Cell Biol. 87, 7791 (2009).

6. Zhao, B., Li, L., Lei, Q. & Guan, K.-L. The Hippo-YAP pathway in organ size control and tumorigenesis: an updated version. Genes Dev. 24, 862874 (2010).

7. Dupont, S. et al. Role of YAP/TAZ in mechanotransduction. Nature 474, 179183 (2011).

8. Calvo, F. et al. Mechanotransduction and YAP-dependent matrix remodelling is required for the generation and maintenance of cancer-associated broblasts. Nat. Cell Biol. 15, 637646 (2013).

9. Lian, I. et al. The role of YAP transcription coactivator in regulating stem cell self-renewal and differentiation. Genes Dev. 24, 11061118 (2010).

10. Dumbauld, D. W. et al. How vinculin regulates force transmission. Proc. Natl Acad. Sci. USA. 110, 97889793 (2013).

11. Sood, A. K. et al. Biological signicance of focal adhesion kinase in ovarian cancer: role in migration and invasion. Am. J. Pathol. 165, 10871095 (2004).

12. Mo, J.-S., Yu, F.-X., Gong, R., Brown, J. H. & Guan, K.-L. Regulation of the Hippo-YAP pathway by protease-activated receptors (PARs). Genes Dev. 26, 21382143 (2012).

13. Morikawa, Y. et al. Actin cytoskeletal remodeling with protrusion formation is essential for heart regeneration in Hippo-decient mice. Sci. Signal. 8, ra41 (2015).

14. Lucas, E. P. et al. The Hippo pathway polarizes the actin cytoskeleton during collective migration of Drosophila border cells. J. Cell Biol. 201, 875885 (2013).

15. Kim, N.-G. & Gumbiner, B. M. Adhesion to bronectin regulates Hippo signaling via the FAK-Src-PI3K pathway. J. Cell Biol. 210, 503515 (2015).

16. Ma, B. et al. Zyxin-Siah2-Lats2 axis mediates cooperation between Hippo and TGF-b signalling pathways. Nat. Commun. 7, 11123 (2016).

17. Schlegelmilch, K. et al. Yap1 acts downstream of a-catenin to control epidermal proliferation. Cell 144, 782795 (2011).

18. Goldmann, W. H. & Ingber, D. E. Intact vinculin protein is required for control of cell shape, cell mechanics, and rac-dependent lamellipodia formation. Biochem. Biophys. Res. Commun. 290, 749755 (2002).

19. Gumbiner, B. M. & Kim, N.-G. The Hippo-YAP signaling pathway and contact inhibition of growth. J. Cell Sci. 127, 709717 (2014).

20. Azzolin, L. et al. YAP/TAZ incorporation in the b-catenin destruction complex orchestrates the Wnt response. Cell 158, 157170 (2014).

21. Yoshigi, M., Hoffman, L. M., Jensen, C. C., Yost, H. J. & Beckerle, M. C. Mechanical force mobilizes zyxin from focal adhesions to actin laments and regulates cytoskeletal reinforcement. J. Cell Biol. 171, 209215 (2005).

22. Oka, T. & Sudol, M. Nuclear localization and pro-apoptotic signaling of YAP2 require intact PDZ-binding motif. Genes Cells 14, 607615 (2009).

23. Komuro, A., Nagai, M., Navin, N. E. & Sudol, M. WW domain-containing protein YAP associates with ErbB-4 and acts as a co-transcriptional activator for the carboxyl-terminal fragment of ErbB-4 that translocates to the nucleus.J. Biol. Chem. 278, 3333433341 (2003).24. Varelas, X. et al. TAZ controls Smad nucleocytoplasmic shuttling and regulates human embryonic stem-cell self-renewal. Nat. Cell Biol. 10, 837848 (2008).

25. Shao, D. D. et al. KRAS and YAP1 converge to regulate EMT and tumor survival. Cell 158, 171184 (2014).

26. Li, Z. et al. Structural insights into the YAP and TEAD complex. Genes Dev. 24, 235240 (2010).

27. Zanconato, F. et al. Genome-wide association between YAP/TAZ/TEAD and AP-1 at enhancers drives oncogenic growth. Nat. Cell Biol. 17, 12181227 (2015).

28. Mackay, D. J., Esch, F., Furthmayr, H. & Hall, A. Rho- and rac-dependent assembly of focal adhesion complexes and actin laments in permeabilized broblasts: an essential role for ezrin/radixin/moesin proteins. J. Cell Biol. 138, 927938 (1997).

29. Benz, P. M. et al. Differential VASP phosphorylation controls remodeling of the actin cytoskeleton. J. Cell Sci. 122, 39543965 (2009).

30. Burridge, K. & Guilluy, C. Focal adhesions, stress bers and mechanical tension. Exp. Cell Res. 343, 1420 (2016).

31. Grashoff, C. et al. Measuring mechanical tension across vinculin reveals regulation of focal adhesion dynamics. Nature 466, 263266 (2010).

32. Hamadi, A. et al. Regulation of focal adhesion dynamics and disassembly by phosphorylation of FAK at tyrosine 397. J. Cell Sci. 118, 44154425 (2005).

33. Ridley, A. J. et al. Cell migration: integrating signals from front to back. Science 302, 17041709 (2003).

34. Amano, M., Nakayama, M. & Kaibuchi, K. Rho-kinase/ROCK: A key regulator of the cytoskeleton and cell polarity. Cytoskeleton (Hoboken) 67, 545554 (2010).

35. Yang, N. et al. Colin phosphorylation by LIM-kinase 1 and its role in Rac-mediated actin reorganization. Nature 393, 809812 (1998).

36. Nobes, C. D. & Hall, A. Rho, Rac, and Cdc42 GTPases regulate the assembly of multimolecular focal complexes associated with actin stress bers, lamellipodia, and lopodia. Cell 81, 5362 (1995).

37. Piccolo, S., Dupont, S. & Cordenonsi, M. The biology of YAP/TAZ: hippo signaling and beyond. Physiol. Rev. 94, 12871312 (2014).

38. Clyman, R. I., Mauray, F. & Kramer, R. H. Beta 1 and beta 3 integrins have different roles in the adhesion and migration of vascular smooth muscle cells on extracellular matrix. Exp. Cell Res. 200, 272284 (1992).

39. Lin, G. L. et al. Activation of beta 1 but not beta 3 integrin increases cell traction forces. FEBS Lett. 587, 763769 (2013).

40. McBeath, R., Pirone, D. M., Nelson, C. M., Bhadriraju, K. & Chen, C. S. Cell shape, cytoskeletal tension, and RhoA regulate stem cell lineage commitment. Dev. Cell 6, 483495 (2004).

41. Gofn, J. M. et al. Focal adhesion size controls tension-dependent recruitment of alpha-smooth muscle actin to stress bers. J. Cell Biol. 172, 259268 (2006).

42. McClatchey, A. I. & Giovannini, M. Membrane organization and tumorigenesis--the NF2 tumor suppressor, Merlin. Genes Dev. 19, 22652277 (2005).

43. Petrie, R. J., Doyle, A. D. & Yamada, K. M. Random versus directionally persistent cell migration. Nat. Rev. Mol. Cell Biol. 10, 538549 (2009).

44. Hersey, P. et al. A randomized phase 2 study of etaracizumab, a monoclonal antibody against integrin av b3,dacarbazine in patients with stage IV metastatic melanoma. Cancer 116, 15261534 (2010).

45. Moroishi, T., Hansen, C. G. & Guan, K.-L. The emerging roles of YAP and TAZ in cancer. Nat. Rev. Cancer 15, 7379 (2015).

46. Finch-Edmondson, M. L. et al. TAZ protein accumulation is negatively regulated by YAP abundance in mammalian cells. J. Biol. Chem. 290, 2792827938 (2015).

47. Czndr, K. et al. Micropatterned substrates coated with neuronal adhesion molecules for high-content study of synapse formation. Nat. Commun. 4, 2252 (2013).

48. Forte, G. et al. Human cardiac progenitor cell grafts as unrestricted source of supernumerary cardiac cells in healthy murine hearts. Stem Cells 29, 20512061 (2011).

49. Langmead, B., Trapnell, C., Pop, M. & Salzberg, S. L. Ultrafast and memory-efcient alignment of short DNA sequences to the human genome. Genome Biol. 10, R25 (2009).

50. Uuskla-Reimand, L. et al. Topoisomerase II beta interacts with cohesin and CTCF at topological domain borders. Genome Biol. 17, 182 (2016).

51. Heinz, S. et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol. Cell 38, 576589 (2010).

52. Chen, T.-W. et al. ChIPseek, a web-based analysis tool for ChIP data. BMC Genomics 15, 539 (2014).

53. Zambelli, F., Pesole, G. & Pavesi, G. PscanChIP: nding over-represented transcription factor-binding site motifs and their correlations in sequences from ChIP-Seq experiments. Nucleic Acids Res. 41, W535W543 (2013).

12 NATURE COMMUNICATIONS | 8:15321 | DOI: 10.1038/ncomms15321 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15321 ARTICLE

54. Digiuni, S. et al. Single cell wall nonlinear mechanics revealed by a multiscale analysis of AFM force-indentation curves. Biophys. J. 108, 22352248 (2015).

55. Pugno, N. M. Flexible nanovectors. J. Phys. Condens. Matter 20, 474205 (2008).

Acknowledgements

The work was supported by the projects no. LQ1605, LQ1601 and LO1304 from the National Program of Sustainability II or I respectively (MEYS CR), by the project FNUSA-ICRC number CZ.1.05/1.1.00/02.0123 (OP VaVpI) and by CEITEC-Central European Institute of Technology (CZ.1.05/1.1.00/02.0068) from the European Regional Development Fund. Moreover, the work was supported by the European Social Fund and European Regional Development Fund-Project MAGNET (Number CZ.02.1.01/0.0/0.0/ 15_003/0000492). We are grateful to Katerina Melajova for helping with cell culture, transfection and western blotting, Carmen Escobedo-Lucea for critical reading, Federico Tidu for scanning electron microscopy analysis and Enrico Cenni for helping with image analysis. N.M.P. is supported by the European Research Council PoC 2015 Silkene No. 693670, by the European Commission H2020 under the Graphene Flagship Core 1 No. 696656 (WP14 Polymer Nanocomposites) and under the FET Proactive Neurobres No. 732344.

Author contributions

G.N., J.O.D.L.C., J.V., C.M., S.P. and F.M.: performed the experiments. G.F. and G.B.S.: conceived and designed the experiments, and wrote the paper. G.C., M.P., J.P. and Z.S.: performed the AFM analysis. N.M.P.: modelled cell adhesion. Z.M., S.P. and M.H.: performed deep sequencing analysis. A.S.-G.: analysis of gene expression.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Nardone, G. et al. YAP regulates cell mechanics by controlling focal adhesion assembly. Nat. Commun. 8, 15321 doi: 10.1038/ncomms15321 (2017).

Publishers note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional afliations.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2017

NATURE COMMUNICATIONS | 8:15321 | DOI: 10.1038/ncomms15321 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 13