The vertebrate embryo undergoes a series of dramatic morphological changes as the body extends to form the complete anterior-posterior axis during the somite-forming stages. The molecular mechanisms regulating these complex processes are still largely unknown. We show that the Hippo pathway transcriptional coactivators Yap1 and Wwtr1 are specifically localized to the presumptive epidermis and notochord, and play a critical and unexpected role in posterior body extension by regulating Fibronectin assembly underneath the presumptive epidermis and surrounding the notochord. We further find that Yap1 and Wwtr1, also via Fibronectin, have an essential role in the epidermal morphogenesis necessary to form the initial dorsal and ventral fins, a process previously thought to involve bending of an epithelial sheet, but which we now show involves concerted active cell movement. Our results reveal how the Hippo pathway transcriptional program, localized to two specific tissues, acts to control essential morphological events in the vertebrate embryo.https://doi.org/10.7554/eLife.31065.001
A key step in the development of the vertebrate embryonic body is the change from the roughly spherical-shaped embryo at the end of gastrulation to the elongated body present at the end of the somite-forming stages when the initial anterior-posterior body plan is fully established (reviewed in Bénazéraf and Pourquié, 2013; Henrique et al., 2015; Kimelman, 2016; Wilson et al., 2009). This process involves complex morphogenetic changes as cells migrate out of a bipotential neuromesodermal progenitor population located at the most posterior end of the embryo and into the presomitic mesoderm and posterior neural tube (Bénazéraf et al., 2010; Goto et al., 2017; Kanki and Ho, 1997; Lawton et al., 2013; reviewed in McMillen and Holley, 2015; Serwane et al., 2017; Steventon et al., 2016).
The factors and mechanisms regulating these morphogenetic events are still largely unknown, but are an active area of research. The mesodermal transcription factors Tbx6/16 and Mesogenin play an essential role in this process as mutants in these genes fail to form the posterior body axis normally (Chalamalasetty et al., 2014; Chapman and Papaioannou, 1998; Fior et al., 2012; Kimmel et al., 1989; Manning and Kimelman, 2015; Nowotschin et al., 2012; Yabe and Takada, 2012; Yoon and Wold, 2000), although the relative importance of each factor is species specific (Kimelman, 2016; Nowotschin et al., 2012). The signaling factors Wnt and Fgf also play essential roles by controlling various aspects of cell movement at the posterior end of the embryo (Bénazéraf et al., 2010; Goto et al., 2017; Lawton et al., 2013; Steventon et al., 2016).
An additional factor that was shown to have a major role in overall early vertebrate body morphogenesis is the Hippo pathway transcription factor Yap1, a key mechanotransduction sensor whose activity is regulated both through protein degradation and nucleocytoplasmic shuttling (Piccolo et al., 2014). In this study (Porazinski et al., 2015), a medaka yap1 mutant (hirame) was shown to have a severely distorted body form that was described as being a flattened embryo, although the specific role of Yap1 in forming the posterior body during somitogenesis was not examined. Yap1 was thus proposed to regulate tissue tension throughout the medaka embryo, thereby opposing the forces of gravity. Based on additional studies using 3D spheroids of human retinal epithelial cells, Yap1 was proposed to act by regulating the transcription of one or more Rho GTPase activating protein genes (arhgaps), although these genes were not shown to be functional targets of Yap1 in the medaka embryo (Porazinski et al., 2015).
Using zebrafish embryos, we show here that simultaneous loss of Yap1 and its paralog Wwtr1 (previously called Taz, Hilman and Gat, 2011; Piccolo et al., 2014), causes a severe defect in the elongation of the embryonic body during the mid-late somitogenesis stages. Surprisingly, we find that Yap1 and Wwtr1 are not ubiquitously expressed throughout the embryo but are specifically localized to the presumptive epidermis and notochord. Analysis of transcriptional changes in yap1-/-;wwtr1-/- double mutants reveals downregulation of a spectrum of genes expressed within the presumptive epidermis and notochord, but with little or no effect on the expression of arhgaps. Examination of the presumptive epidermis in yap1;wwtr1 double mutants also revealed profound defects, leading to a new understanding of the basis of the initial steps in the formation of the median fin fold, the precursor to the dorsal and ventral unpaired fins. Importantly, we discovered that the yap1;wwtr1 double mutants have aberrant Fibronectin (Fn) deposition, and show that interfering with Fn function results in embryos that appear similar to the yap1;wwtr1 double mutants, with inhibition of posterior elongation and a failure of presumptive epidermal cells to move into the nascent fin fold. We propose that Yap1 and Wwtr1 regulate transcription programs specifically within presumptive epidermal and notochord cells, allowing them to adhere to the surrounding tissues and form the extracellular matrix that is necessary for both posterior body elongation and presumptive epidermal cell migration to produce the median fin fold.
We generated new mutant alleles of yap1 and wwtr1 using CRISPR (Lai et al., 2017). yap1bns19 contains a 41 bp deletion in exon 1, which encodes part of the TEAD binding domain, whereas wwtr1bns35 contains a 29 bp insertion in exon 2, which encodes part of the WW domain (Figure 1A). Both frameshift mutations are predicted to lead to truncated proteins (p.Ile39Argfs*72 and p.Pro145Glnfs*15, respectively). yap1;wwtr1 double homozygous mutants (yap1-/-;wwtr1-/-) exhibit a severe defect in forming the posterior body, first clearly visible at the 15–16 somite stage (Figure 1B and Video 1). Although not studied, the same defect was observed previously in zebrafish double mutant embryos with different yap1 and wwtr1 mutations (Figure 1G in Miesfeld et al., 2015; Supplemental Figure 6I in Nakajima et al., 2017), demonstrating that the defects we saw were not due to the particular alleles we had identified. Prior to the 15–16 somite stage we did not observe any differences among the embryos during gastrulation or early somitogenesis, indicating that the role of zygotic Yap1/Wwtr1 becomes essential as the tail begins to form. yap1;wwtr1 double mutants formed normally proportioned somites, although the somites never obtained the chevron shape seen in wild-type or sibling embryos (Figure 1C). The body was not flattened by gravity (Figure 1E,F) as described for medaka yap1 mutants (Porazinski et al., 2015). We also observed that among the embryos that arose from crosses of yap1+/-;wwtr1+/- adults, a proportion showed less severe defects after the embryonic stages, including milder posterior defects (Figure 1D). When these embryos were genotyped (n = 15) they were all found to be yap1-/-;wwtr1+/-, showing that while one copy of yap1 alone is sufficient for posterior body formation, one copy of wwtr1 is not sufficient, in the context of the absence of its paralogue.
During the somitogenesis stages, separate groups have reported somewhat different distributions of yap1 mRNA in early zebrafish embryos (He et al., 2015; Jiang et al., 2009; Loh et al., 2014) and our data not shown). In medaka, only the temporal expression of yap1 using PCR analysis was described (Porazinski et al., 2015). The expression of wwtr1 mRNA in neither system has been reported during the somitogenesis stages. Since the protein distribution of these factors, whose levels and cellular localization are controlled post-translationally (Piccolo et al., 2014), have not been analyzed during somitogenesis in any vertebrate embryo, we examined the expression of Yap1 and Wwtr1 in zebrafish at mid-somitogenesis, when the mutant phenotype first becomes clear, in the posterior of the embryo. Surprisingly, Yap1 was strongly expressed in the presumptive epidermis at the 18-somite stage, as well as along the axis (Figure 2A). In both cases, the staining was nuclear, which is indicative of active Yap1 (Piccolo et al., 2014). The presumptive epidermal expression of Yap1 was confirmed with fluorescent in situ hybridization (FISH) using probes for the previously characterized presumptive epidermal markers keratin4 and keratin8 (Figure 2—figure supplement 1). The expression of Yap1 protein in the notochord was ascertained using FISH with a probe for the notochord marker no tail, also called ta1 (Schulte-Merker et al., 1994, Figure 2B). Importantly, yap1 mutants did not exhibit any Yap1 antibody staining, demonstrating that the antibody is specific to Yap1 (Figure 2C). Similarly, Wwtr1 was also localized to the presumptive epidermis and notochord (Figure 2—figure supplement 2), in keeping with the overlapping role of these two factors (Zanconato et al., 2015). Based on these findings we also examined the expression of Yap1 in medaka to determine if zebrafish and medaka were significantly different. As in zebrafish, medaka Yap1 is also predominantly expressed in the presumptive epidermis and notochord in the posterior (Figure 2—figure supplement 2). Altogether, these data suggest that Yap1 and Wwtr1 are not acting ubiquitously throughout the embryo but function in localized domains of the posterior body to control morphogenesis.
In order to understand precisely how Yap1 and Wwtr1 are regulating posterior morphogenesis, the tailbud and presomitic mesoderm up to the third newest somite (S-III) were manually dissected at the 16–18 somite stage from embryos collected from a cross of yap1+/-;wwtr1+/- adults, and which had been separated by phenotype into siblings and yap1;wwtr1 double mutant embryos, and then subjected to RNA-seq analysis (Supplementary file 1). As shown in Table 1, most of the top hits have not been previously analyzed, with only four having been localized to the presumptive epidermis or notochord during the somitogenesis stages (data from ZFIN.org). We therefore produced in situ hybridization probes to the remaining genes shown in Table 1 and found that their expression all localized to the presumptive epidermis or notochord at the 18-somite stage (Table 1). For example, ecrg4b is expressed in the presumptive epidermis at this stage, but is not detected in yap1;wwtr1 double mutants, confirming the RNA-seq results (Figure 3A,B). Similarly, wu:fc23c09, a zebrafish podocan family member, is expressed throughout the notochord in sibling embryos but is mostly missing from the posterior notochord in yap1;wwtr1 double mutants (Figure 3—figure supplement 1A,B), which is within the region that was isolated for the RNA-seq analysis.
We also used Gene Ontology (GO) analysis to further analyze our RNA-seq results. In accord with the primary targets of Yap1/Wwtr1 being expressed in the presumptive epidermis, categorization by Biological Process revealed genes expressed in either the epidermis or the fin, which is primarily an epidermal structure (Figure 4A). Notochord did not show up as statistically significant in this analysis both because only a small number of Yap1/Wwtr1 targets are expressed in the notochord (Table 1) and because several genes expressed in the notochord in zebrafish embryos are not identified as being expressed in the notochord in the GO database. Interestingly, the major targets identified using Kegg Pathway analysis comprised genes involved in the cytoskeleton, or in the interactions of cells with their environment (Figure 4B). Finally, we examined genes previously identified as Hippo pathway transcriptional targets in large-scale screens in other systems (Zhang et al., 2009; Zhao et al., 2008). The top 10 targets are shown in Table 2, and all except one of these are expressed in the presumptive epidermis or notochord. The one exception is nuak1a, which is weakly expressed in the ventral mesoderm and thus appears to be another tissue regulated by Yap1/Wwtr1, albeit for a very limited set of target genes. Because the genes expressed in Table 2 have more modest fold-changes, we also validated that their expression decreased using quantitative PCR (Figure 4—figure supplement 1).
To further understand the regulation of these genes by Yap1 and Wwtr1, we made use of two transgenic lines that allow temporal regulation of Yap1/Wwtr1, Tg(hsp70:RFP-DNyap), which expresses a dominant-negative form of Yap1 under heat shock control, and Tg(hsp70:RFP-CAyap), which expresses a constitutively active Yap1 under heat shock control (Mateus et al., 2015). These lines are referred to here as HS:DN-yap and HS:CA-yap, respectively. Embryos were heat shocked at the 4-somite stage and then fixed at 18 somites. Inhibition of Yap1/Wwtr1 function with HS:DN-yap caused a downregulation of ecrg4b expression, although not as strongly as observed in the yap1;wwtr1 double mutants (Figure 3C,D). Intriguingly, enhancement of Yap1/Wwtr1 levels with HS:CA-yap caused an upregulation of ecrg4b, indicating that Yap1/Wwtr1 are not working to activate transcription to maximum levels (Figure 3E,F). In contrast, while HS:DN-yap caused an inhibition of wu:fc23c09 expression in the posterior notochord similar to that observed in the mutants, HS:CA-yap did not upregulate wu:fc23c09, indicating that Yap1/Wwtr1 already functions at maximal levels in this tissue, or at least for this target (Figure 3—figure supplement 1C,E). Interestingly, the increased expression of ecrg4b caused by CA-yap was restricted to the presumptive epidermis (Figure 3G,H), and the absence of wu:fc23c09 expression caused by DN-yap was restricted to the posterior notochord (Figure 3—figure supplement 1D). These observations indicate that Yap1/Wwtr1 must work with tissue-specific factors to control gene expression.
Since the analysis of a yap1 mutant in medaka proposed that arhgaps were key targets of Yap1 (Porazinski et al., 2015), we examined arhgaps expression in our RNA-seq data. Of the many arhgaps found in zebrafish, only two weakly expressed arhgaps (based on the RNA-seq data), arhgap12a and arhgap27, were observed to change between siblings and mutants, with arhgap12a down 1.8-fold in the mutants and arhgap27 down 1.9-fold (Supplementary file 1). However, we also found that many presumptive epidermal genes were decreased approximately 2-fold in the mutants, including all of the abundant keratins. For example, krt4 was down 2.4-fold, krt8 was down 2.7-fold, and krt18 was down 1.8-fold. We also examined the expression of these two arhgap genes by in situ hybridization to determine where they were expressed and to look for possible changes in expression in the mutants. We found that both arhgaps are expressed in the presumptive epidermis, and found no obvious spatial difference in their expression between siblings and mutants (Figure 3—figure supplement 2). One arhgap gene is expressed in the posterior notochord, arhgap17a (data from ZFIN.org), but this gene showed no change between siblings and mutants in our RNA-seq data. Our data suggest that arhgaps are not the key targets of Yap1/Wwtr1 in regulating posterior body morphogenesis in zebrafish (see Discussion), and instead identify a number of genes whose expression is highly regulated by these transcription factors in either the presumptive epidermis or notochord, which are the same tissues in which we found Yap1 and Wwtr1 expression.
The localization of Yap1 and Wwtr1 to the presumptive epidermis and notochord, together with the changes in gene expression in these tissues, indicated that defects in one or both tissues were the source of the posterior morphogenetic defects. Since the notochord cells undergo active morphogenetic movements during somitogenesis in many species (Jiang and Smith, 2007; Keller et al., 2000; Lee and Anderson, 2008; McMillen and Holley, 2015), it seemed possibly to be the key target tissue for Yap1/Wwtr1 regulated posterior morphogenesis. A zebrafish mutant that lacks all notochord formation, noto1/floating head, was previously identified in large-scale genetic screens, and shown to have a shortened body axis at 24 hpf (30 somite stage, Odenthal et al., 1996; Stemple et al., 1996; Talbot et al., 1995). Since posterior morphogenesis at earlier stages had not been characterized, we obtained the original noto1n1 mutant line (Talbot et al., 1995) and examined the formation of the posterior body. At the 16-somite stage there was a minor but recognizable difference between noto mutants and their siblings, and this difference became more pronounced as embryos developed to the 24-somite stage, paralleling the time course of alterations in the yap1;wwtr1 double mutants (Figure 4—figure supplement 2). However, since the defects in noto mutants were much less severe than those in yap1;wwtr1 double mutants, our results suggested that alterations in the notochord are only a small part of the overall yap1;wwtr1 double mutant phenotype.
To determine whether the presumptive epidermis was altered in yap1;wwtr1 double mutants we first examined the expression of cytokeratins using two antibodies, panKr1-8, which recognizes a broad spectrum of cytokeratins, and 79.14, which recognizes a more limited subset (Conrad et al., 1998). We found that PanKr1-8 labeled the entire presumptive epidermis, whereas 79.14 only labeled the outermost cells (Figure 5A,C). With both antibodies we observed that both the dorsal and ventral midlines at the 24-somite stage were 1–2 cell layers thick in mutants, whereas in siblings they were several layers thick (Figure 5A–D).
Because the cytokeratin antibodies do not work well at earlier stages, we turned to an antibody to the transcription factor Tp63, which marks the formation of the basal layer of the nascent presumptive epidermis from mid-gastrula stages onward (Bakkers et al., 2002; Lee and Kimelman, 2002) in order to examine the presumptive epidermis during mid-somitogenesis stages. From the 18-somite stage to the 24-somite stage (a period of 3 hr), the embryo undergoes a rapid change in the thickness of the presumptive epidermis along the dorsal and ventral midlines to form the nascent median fin fold (MFF) that ultimately gives rise to the dorsal and ventral unpaired fins (Figure 6A), as was previously observed in histological sections (Abe et al., 2007; Dane and Tucker, 1985). Importantly, the increase in thickness is due to an increase in the number of Tp63-positive cells, while the number of Tp63-negative presumptive epidermal cells does not noticeably change. At the 18-somite stage, yap1;wwtr1 double mutants appeared unaffected, with a single layer of Tp63-positive cells (Figure 6B, and data not shown). During the next three hours, however, while the siblings increased the thickness of the presumptive epidermis at the midline (Figure 6C), the double mutants were essentially unchanged from the 18-somite stage (Figure 6D). Examination of epidermal cell proliferation using an antibody to phospho-Histone H3 and of cell death using acridine orange did not show noticeable differences between siblings and yap1;wwtr1 double mutants (data not shown), suggesting that the failure to increase the epidermal thickness at the midline in the double mutants is due to morphogenetic defects.
Previous studies of fin fold formation focusing on the dorsal side of the embryo and using cross sections of temporally staged embryos suggested that the MFF forms from shape changes of generally static cells (Abe et al., 2007; Dane and Tucker, 1985, Figure 7A), whereas analysis using in situ hybridization with the MFF markers dlx5a and dlx6a suggested movement of the whole presumptive epidermal sheet toward the midline as the underlying neural tissue rolled up to form the neural rod (Heude et al., 2014). Because previous studies had used only fixed tissue, we used live imaging to understand the cellular dynamics in both wild-type and mutant embryos as the presumptive epidermis changes shape. We also imaged the ventral fin fold as during these early stages it undergoes much more dramatic changes in overall morphology than the dorsal fin fold, although we generally observed the same results on both sides. Embryos were injected with mRNA that encodes both a membrane GFP and a nuclear version of the photoconvertible protein Kikume (NLS-Kikume/EGFP-CAAX), which was photoconverted right before imaging in order to be able to follow membranes and nuclei in separate fluorescent channels, and then live imaged on a spinning disc confocal microscope between the 20- and 26-somite stages. Surprisingly, we found that in wild-type embryos the presumptive epidermal cells are highly active (Video 2). The cells were constantly changing shape as well as their interactions with neighbors, with new cells moving in at the basal side of the fin fold (Figure 7B, arrowheads). By following the newly emerging cells backwards in time through different Z-planes, we could observe cells moving in from both lateral sides toward the midline. Once cells arrived at the midline they tended to stay there (Figure 8A, black lines). In yap1;wwtr1 double mutants, the cells were still moving and changing shape (Video 3), although not as dramatically as in wild-type. New cells showed up at the midline over time, but other cells departed such that there was no net accumulation of new cells at the midline (Figure 7C), in contrast to what was observed in the siblings. Following the cells over time through different Z-planes, we could watch individual cells sliding across the midline and only occasionally staying in place (Figure 8A, red lines). We also observed that mutant cells moved faster than wild-type cells (Figure 8B), and remained in an elongated shape along the body axis (length) whereas wild-type cells extended along the apical-basal axis (height, Figure 8C). We conclude that formation of the MFF is normally a very active process with new presumptive epidermal cells constantly moving into the midline and staying there to build up the presumptive epidermal thickening, in contrast to the bending of a presumptive epidermal sheet envisioned from the static images. In yap1;wwtr1 double mutants, the presumptive epidermal cells are highly motile, but are unable to coalesce at the midline, resulting in the failure to form the nascent MFF.
We wondered if there might be a common cause that connected the posterior body and presumptive epidermal morphogenetic defects, and if it was related to the observed function of Yap1/Wwtr1 in the presumptive epidermis and notochord. A series of elegant studies has demonstrated that Integrin-Fibronectin interactions play an essential role in body elongation during the somitogenesis stages as well as forming the boundaries between somites (Dray et al., 2013; Jülich et al., 2005; Koshida et al., 2005; McMillen and Holley, 2015). Intriguingly, Fn is present ubiquitously but it assembles specifically at the border between somites, at the interface between the somites and presumptive epidermis, and at the border between the somites and notochord, due to the activation of Integrins within the somitic and presomitic mesoderm at these borders (Dray et al., 2013; Jülich et al., 2015). Within the rest of the mesoderm, Cadherins act to prevent integrin activation, leading to Fn assembly only at the borders (Jülich et al., 2005). Importantly, the Fn underlying the presumptive epidermis and surrounding the notochord is proposed to provide the ECM necessary for body elongation by mechanically coupling the paraxial mesoderm to the presumptive epidermis and notochord (reviewed in McMillen and Holley, 2015). We therefore examined Fn distribution in sibling and mutant embryos.
In sibling embryos, Fn was clearly seen in the intersomitic furrows (Figure 9A), and in a continuous line underlying the presumptive epidermis and surrounding the notochord at the embryonic midline (Figure 9G). At the lateral edges of the somites, where they contact the presumptive epidermis, rosettes of Fn fibers were observed (Figure 9C,E). The distribution of Fn in medaka embryos at the same stage was the same (Figure 9—figure supplement 1). yap1;wwtr1 double mutants had normal Fn distribution between the somites (Figure 9B), whereas the Fn distribution at the lateral edges was very abnormal, with gaps in the rosettes as well as regions with enhanced Fn distribution (Figure 9D,F). Similarly, Fn at the midline also showed pronounced gaps in the region underlying the presumptive epidermis (Figure 9H, arrows). Expression of Fn around the most newly formed notochord at the posterior end of the embryo appeared mostly unaffected in yap1;wwtr1 double mutants, but it disappeared in the region where new somites form at the anterior end of the presomitic mesoderm (Figure 9H, arrowheads), indicating that the notochord initially assembles Fn but is unable to maintain it. Importantly, the defects in Fn only occurred at the borders with the tissues that express Yap1/Wwtr1 (the notochord and presumptive epidermis), and not within the mesoderm where these factors are not expressed.
Embryos not expressing Fn are severely deformed, forming almost no recognizable body (Jülich et al., 2005). To achieve a partial knockdown of Fn function, we produced a dominant-negative Fn1b (Fn∆C), since fn1b is the predominant fn expressed in the tailbuds in our RNA-seq data, based on a construct used to inhibit Fn in Xenopus embryos (Rozario et al., 2009; Schwarzbauer and DeSimone, 2011). Injection of 400 pg fn∆C mRNA into 1 cell embryos caused gastrulation defects in approximately half of the embryos. Embryos that completed gastrulation showed various degrees of body extension defects, with some embryo having phenotypes similar to those seen in yap1;wwtr1 double mutants (Figure 10A–C). Embryos with an elongation defect had somites that did not become chevron-shaped (Figure 10D,E), and presumptive epidermal cells that did not form the multi-layered MFF (Figure 10F,G) as in yap1;wwtr1 double mutants. Interestingly, in these embryos somite borders still formed despite the importance of Fn for somite boundaries (Jülich et al., 2005; Koshida et al., 2005), indicating that posterior body and presumptive epidermal morphogenesis is more sensitive to disruption of Fn than is somitogenesis.
In order to understand the defects caused by disrupting Fn in more detail, we live imaged embryos expressing Fn∆C (Video 4). In general we saw the same trends as with the yap1;wwtr1 double mutants, although the effects were not as strong, potentially due to the variable distribution of protein from injected mRNA. As was observed in the fixed tissue (Figure 10G), the presumptive epidermis in the Fn∆C-expressing embryos was a single layer thick (Figure 10—figure supplement 1A). As with the yap1;wwtr1 double mutant cells, the cells expressing Fn∆C often remained at the midline for short intervals of time (Figure 10—figure supplement 1B). The Fn∆C-expressing cells showed variable rates of movement, with some cells moving faster than wild-type cells and none moving slower (Figure 8B). Moreover, Fn∆C-expressing cells did not adopt the wild-type morphology over time, and instead remained in an elongated shape along the body axis as was observed in the yap1;wwtr1 double mutant cells (Figure 8C). Thus, disruption of Fn causes similar effects to those see in yap1;wwtr1 double mutants, both in terms of overall phenotype as well as in the behavior of presumptive epidermal cells.
Since Yap1/Wwtr1 are mechanosensors and regulated by extracellular signals (Piccolo et al., 2014), we asked whether alterations in Fn might also regulate their activity. Fn∆C expressing embryos exhibiting a yap1;wwtr1 double mutant-like phenotype (Figure 10B,C) were analyzed using the anti-Yap1 antibody. Expression of Yap1 was the same as in wild-type (Figure 2A), with strong nuclear expression in the epidermis and notochord (data not shown). Thus, the alterations in Fn seen in the yap1;wwtr1 double mutants are not in turn regulating Yap1.
The disruptions in Fn surrounding the notochord and under the presumptive epidermis suggested that these tissues might have reduced adhesion to their neighboring tissues. An undulating notochord is indicative of reduced mechanical coupling between the notochord and paraxial mesoderm (Dray et al., 2013), and we also observed an undulating notochord in yap1;wwtr1 double mutant embryos (Figures 1F and 9H). We also observed that the presumptive epidermis seemed less firmly attached to the underlying tissue when we microdissected tailbuds for RNA-seq analysis. To examine this phenotype more directly, we injected embryos with mRNA encoding a membrane-targeted GFP and then live imaged embryos at the 16- and 20-somite stages. Whereas the presumptive epidermis remained attached to the underlying tissues in sibling embryos (Figure 11A,C), in yap1;wwtr1 double mutant embryos the presumptive epidermis began to separate from the somites at the 16-somite stage, and this phenotype was very pronounced at the 20-somite stage (Figure 11B,D). Taken together, our results support a model in which Yap1 and Wwtr1 regulate gene expression in the notochord and presumptive epidermis in order to promote cell adhesion and fibronectin deposition, which is necessary for both posterior body elongation and presumptive epidermal morphogenesis. We elaborate on this model below.
We show here that zebrafish yap1;wwtr1 double mutants have a severe defect in body elongation beginning at approximately the 15–16 somite stage. Although this process was not studied in medaka, the yap1/hirame mutants also lack an extended tail (Figure 1 in Porazinski et al., 2015), demonstrating that this is a conserved process. We also found that yap1;wwtr1 double mutants fail to form the nascent MFF, which was not mentioned for yap1/hirame, most likely because the flattening of the medaka embryo precluded conclusions about the state of the MFF. What we did not find similar was the flattening of the body in zebrafish yap1;wwtr1 double mutants (Figure 1B,E,F), which was the basis for the proposed key role of Yap in opposing the forces of gravity (Porazinski et al., 2015). Why these two species differ in this regard is not clear at this point but possibly due to the fact that medaka has a hard chorion whereas zebrafish has a very soft chorion, and thus medaka may have evolved to specifically resist the pressure of the embryo against the chorion.
We find that Yap1 and Wwtr1 are not present ubiquitously throughout the embryo but instead are both specifically localized to the presumptive epidermis and notochord. Intriguingly, we also find that medaka Yap1 is similarly localized, which was not previously known. The localization of the zebrafish proteins is in strong accord with our RNA-seq data, which shows major changes in gene expression in presumptive epidermal and notochordal genes. While these two tissues are very different in function and lineage, we see alterations in Fn assembly, which underlies the presumptive epidermis and overlies the notochord. Importantly, we do not see defects in Fn assembly between the somites, which supports the finding that Yap1 and Wwtr1 do not function within the mesoderm, and also shows that there is not a general defect in Fn assembly. Defects in Fn assembly were also seen in yap1/hirame mutant embryos between the lens and retina (Porazinski et al., 2015). Surprisingly, altered Fn in yap1/hirame mutants was also reported throughout the neural tube whereas we do not find any Fn assembled in the neural tube in zebrafish, as previously reported (Dray et al., 2013), or in medaka, which is more in keeping with the normal assembly of extracellular matrix between tissues. The importance of Fn in zebrafish is shown with the use of a dominant-negative mutant, Fn∆C, that phenocopies both the extension defect and presumptive epidermal cell migration defects seen in yap1;wwtr1 double mutants. These results, taken together with our evidence of decreased adhesion of the notochord and presumptive epidermis to the surrounding tissues, supports a model in which Yap1 and Wwtr1, functioning redundantly within the presumptive epidermis and notochord, work to promote posterior body elongation and presumptive epidermal morphogenesis necessary for MFF formation (Figure 12A).
The defects in yap1;wwtr1 double mutants begin at the 15–16 somite stage, which is the time when the tailbud everts from the body, and the body begins to achieve a straight conformation along the anterior-posterior axis rather than being curved around the yolk. Indeed, yap1;wwtr1 double mutants remain in a curved conformation (Figure 1B). Part of this defect can be attributed to defects in the notochord as we found that noto mutants also begin to show a phenotype at the 15–16 somite stage, but the defect is much less severe than that seen in yap1;wwtr1 double mutants. These results point to the presumptive epidermis, the other major site of Yap1/Wwtr1 expression, as also being involved in formation of the posterior body. The key factor linking the presumptive epidermis and notochord to body extension is the ECM, and specifically Fibronectin (reviewed in McMillen and Holley, 2015). Fn couples the underlying mesoderm to the notochord and presumptive epidermis, which allows forces from convergence-extension within the paraxial mesoderm to allow the body to extend. The notochord also undergoes convergence-extension, and this helps the body extend as long as the notochord is coupled to the mesoderm (Glickman et al., 2003; Keller et al., 1992; Steventon et al., 2016; Weliky et al., 1991; Wilson and Keller, 1991). Finally, the spinal cord, which is connected to the presumptive epidermis, and the notochord undergo volumetric growth as the body extends, further assisting in the extension of the body (Steventon et al., 2016). In noto mutants, the notochord component is removed but the presumptive epidermis remains. When this process is disrupted such as in the case of our dominant-negative Fn injections, or when Integrin function is knocked down (Dray et al., 2013), posterior extension is blocked. However, unlike in the Integrin knockdown experiments, yap1;wwtr1 double mutants form normal somites because the intersomitic Fn deposition is unaffected. Intriguingly, the observation that the dose of dominant-negative Fn we used produced embryos with defective body extension but normal somites argues that body extension is more sensitive to disruption than the formation of somite boundaries. Possibly because the Fn underlying the presumptive epidermis and surrounding the notochord is under much more force from the convergence-extension movements of the presomitic mesoderm that extend the body (Steventon et al., 2016) than is the Fn that separates the somites, this pool of Fn may be much more sensitive to even minor perturbations.
We also found decreases in notochordal and presumptive epidermal cell adhesion in yap1;wwtr1 double mutants. In living embryos we see a separation of the prospective epidermis from the underlying somitic tissue, which becomes progressively more pronounced as the overall phenotype of the embryos increases. In addition, the undulating notochord we observed in the mutants is indicative of decreased mechanical coupling between the notochord and paraxial mesoderm as is also seen in Integrin knockdown embryos (Dray et al., 2013). Since we see no changes in fibronectin expression levels and little or no changes in integrin levels in our RNA-seq analysis (Supplementary file 1), our results suggest that there are defects within the presumptive epidermis and notochord that results in a failure to assemble Fn. While recent results have shown how Cadherins suppress Fibronectin assembly between adjacent cells within the paraxial mesoderm by keeping Integrins in an inactive conformation (Jülich et al., 2015), the factors that regulate the activation of Integrin within the early embryo are not known. Our RNA-seq data of genes down-regulated in yap1;wwtr1 double mutants provides a valuable dataset of genes worth further investigation, most of which have not been analyzed even at the level of in situ hybridization let alone functionally studied.
Our analysis of the presumptive epidermis in yap1;wwtr1 double mutants has provided new insight into the mechanisms that form the nascent MFF from presumptive epidermal cells. We find that the presumptive epidermal cells are highly active and constantly changing shape when imaged in live embryos. They actively move toward the dorsal and ventral midlines as individual cells and then change their overall shape, with an increase along the apical-basal axis. At later stages, as laminin is used to stabilize the nascent MFF (Webb et al., 2007) by forming connections between apposing cells (Dane and Tucker, 1985), we expect that cell motility is reduced allowing the unpaired fins to extend away from the body. In yap1;wwtr1 double mutants, the presumptive epidermal cells are still highly motile but they transit through the midline rather than staying there, and they remain in an elongated shape along the body axis with no increase in the apical-basal direction, all resulting in a failure to form the MFF. Intriguingly, we can largely recapitulate this effect by inhibiting Fn assembly using dominant-negative Fn injections. In these embryos, the body fails to elongate, the presumptive fin fold does not form, and the cells transit through the midline and remain in an elongated shape. Interestingly, inhibition of cell movement with the myosin inhibitor Blebbistatin did not result in defects in Fn assembly under the epidermis (or around the notochord), demonstrating that the Yap1/Wwtr1-dependent defects in Fn assembly cause migration defects, rather than defective migration of cells leading indirectly to altered Fn assembly (data not shown). These results suggest that in addition to its role in posterior elongation, the Fn matrix provides a substrate for the presumptive epidermal cells to adhere to, which allows them to converge from the right and left sides of the embryo and push against each other to initiate MFF formation (Figure 12B). We call this process concerted migration, since the fin fold forms from the action of the presumptive epidermal cells migrating dorsally or ventrally and pushing against each other, creating the force to form the fin fold, and increasing the cells’ height. These processes then allow the basal surfaces of juxtaposed MFF cells to adhere to each other, with the resulting formation of cross-linking fibers to stabilize the structure (Dane and Tucker, 1985). Thus, Fn does double duty in extending the body and creating the fins that the elongated body uses to swim.
Transcriptional analysis of the genes regulated in the tailbud by Yap1/Wwtr1 identified many genes with large changes in expression, including 8 of our top 12 hits which had not previously been analyzed by in situ hybridization in the early embryo. Importantly, all these genes were expressed in the presumptive epidermis or notochord, exactly the tissues in which we found Yap1 and Wwtr1 expressed. While many of the genes further down the list have not been characterized, those that have are also expressed in the notochord or presumptive epidermis. An outlier at position 36 on our list with a 4.4-fold change is myosin light chain (mylpfa). However, in the tailbud mylpfa is expressed specifically in the adaxial cells (von Hofsten et al., 2008), the small population of future slow muscle cells that originate directly next to the notochord. Since adaxial cell formation depends on interactions with the notochord (Jackson and Ingham, 2013), the decreased mylpfa expression may be a secondary effect of alterations to the notochord.
Our results suggest that the arhgaps are unlikely to be the key genes regulated by Yap1/Wwtr1 in zebrafish. Two weakly expressed presumptive epidermal arhgaps, arhgap27 and arhgap12a, are found far down on our list of dysregulated genes ordered by fold-change (positions 254 and 300, respectively), among a long list of genes that change approximately 2-fold. Moreover, the only notochord expressed arhgap, arhgap17a, is not changed in its expression levels in yap1;wwtr1 double mutant tails. In human 3D spheroids formed using a human retina pigmented epithelial cell line, ARHGAP18 was identified as one of 40 genes down-regulated when YAP levels were reduced, and reduction of ARHGAP18 by siRNA caused a similar phenotype to siRNA reduction of YAP (Porazinski et al., 2015). In medaka yap1 mutants, however, arhgap18 was reduced just 24% from control levels, and inactivation of arhgap18 did not produce a recognizable phenotype (Porazinski et al., 2015). In zebrafish, arhgap18 is not expressed in the tailbud, nor is its closest relative, arhgap28 (37% amino acid identity). Similarly, ARHGAP29 was found to be down-regulated 50% in a human gastric cancer cell line after YAP1 knock out, but the levels of arhgap29 were unchanged in the tailbud of zebrafish yap1;wwtr1 double mutants. These results show that the Yap1 transcriptional program is context dependent and the results found in cell culture can not be inferred to be true in the embryo (Qiao et al., 2017).
Our data also show that Yap1/Wwtr1 is mostly regulating a completely different suite of genes in the presumptive epidermis than in the notochord. Of the known genes expressed in both tissues, the highest scoring gene is laminin α5 (position 86), which is decreased 3-fold in yap1;wwtr1 double mutants. However, we previously identified a mutant in this gene and saw no effect on posterior body formation, and no effect on fin-fold formation prior to 30 hpf (Webb et al., 2007). Our results suggest therefore that Yap1 and Wwtr1 are responsible for broad transcriptional regulatory networks within both the presumptive epidermis and notochord that allow these cells to adhere to the surrounding tissue and assemble the extracellular matrix.
Our results also reveal that Yap1/Wwtr1 do not activate a generic transcriptional program in the early embryo, and instead work in concert with tissue specific factors. This was clearly shown in the case of ecrg4b which is expressed in the epidermis but not in the notochord, whereas Yap1 and Wwtr1 are active in both tissues. Moreover, when a constitutively active Yap was ubiquitously expressed, ecrg4b expression increased only in the epidermis but was not activated in other tissues. Finally, well known Hippo pathway transcriptional targets such as cyclin D1 and aurora kinase A were not down-regulated in the yap1;wwtr1 double mutants according to our RNA-seq data, also supporting the idea that there is not a generic pathway activated by these factors. Instead, our data argue that Yap1 and Wwtr1 enhance the expression of particular genes in concert with tissue-specific transcription factors, but they themselves do not provide tissue specificity.
Intriguingly, a newly published paper has shown that YAP1, which is expressed in the murine presomitic mesoderm (PSM), has an essential role in repressing oscillations of the somite clock (Hubaud et al., 2017). Since this repression is a key part of the regulatory system (Hubaud et al., 2017), and since Yap1 (or Wwtr1) is not expressed in the fish PSM in either zebrafish or medaka, it raises the question of what factor might be playing a role in fish similar to that identified for YAP1 in mouse. As the key targets of YAP1 in the murine PSM are identified, it will be important to compare their regulation in fish.
The yap1bns19 and wwtr1bns35 CRISPR mutant lines are described elsewhere and genotyped as described (Lai et al., 2017). Double mutant yap1-/-;wwtr1-/- embryos were identified by phenotype whereas the sibling embryos had a wild-type phenotype. The noton1 line (Talbot et al., 1995) was obtained from the Zebrafish International Resource Center and genotyped using the provided instructions. The transgenic lines Tg(hsp70:RFP-DNyap)zf621 and Tg(hsp70:RFP-CAyap)zf622 (Mateus et al., 2015) were obtained from Dr. Ken Poss (Duke). These transgenic fish were genotyped by PCR with the primers Yap 394 F (ACTACACCTGCGACTTCAAGA) and Yap 394 R (CCAGCAGTACCTGCATCTGTA), which produce a 394 bp fragment. Embryos were heat shocked at the 4-somite stage at 40 C for 30 min in a circulating water bath. Fixed medaka embryos were kindly provided by Siddharth Ramakrishnan (University of Puget Sound).
A construct was produced in the vector CS2 that has a nuclear localized version of the photoconvertible protein Kikume, followed by the viral 2A sequence, followed by a membrane localized EGFP (NLS-Kikume/EGFP-CAAX). This was used to make mRNA using the mMessage Machine SP6 kit (Ambion). 15 pg of mRNA were injected at the 1 cell stage, and imaged on a spinning disc confocal microscope (3I) using a 40x water immersion lens. Embryos were mounted laterally in a glass-bottomed dish and immobilized with tricaine. Images were acquired every 5 min with Z-axis intervals of 1 μm for 5 hr. Image processing was done using the Fiji plugin for ImageJ.
A fragment of the Fn1b coding region from the start codon to amino acid 630 was amplified with the primers ggcggcggatccgAAATGACCCGTGAGTCAGTAAAGAG and ggcggcatcgATGGGAATTGGGCTGATTTCCAG, and cloned into the BamH1 and ClaI sites of the CS2 vector to make mRNA using the mMessage Machine SP6 kit. 1 cell embryos were injected with 100 or 150 pg of mRNA.
Alkaline phosphatase in situ hybridizations were performed using digoxygenin-labeled probes. The primers used to produce probes for this paper not previously published are shown in Supplementary file 3, and all were cloned into the Bluescript vector (Stratagene). Fluorescent in situ hybridizations were as described (Lauter et al., 2011), except that a primary antibody (anti-Yap or anti-Wwtr1) was used after the first tyramide step, followed by the use of a fluorescent secondary antibody. For immunofluorescence, embryos were fixed in 4% paraformaldehyde (PFA) 2 hr at room temperature or overnight at 4 C, except for the keratin antibodies, which were fixed with Dent’s fixative for 3 hr at room temperature. The anti-Tp63 antibody was purchased from GeneTex (Irvine, CA, USA; GTX124660) and used at a 1:800 dilution according to the manufacturer’s protocol. The anti-Fibronectin antibody was purchased from Sigma (F3648) and used at a 1:100 dilution following a protocol kindly provided by Scott Holley. The anti-Yap (4912) and Wwtr1 (D24E4) antibodies were purchased from Cell Signaling Technology and both were used at a 1:100 dilution with standard immunofluorescence conditions. Note that the Wwtr1 antibody is described as a Yap/Taz antibody, but in zebrafish it is Wwtr1-specific (Figure 2—figure supplement 2; Lai et al., 2017; Miesfeld et al., 2015). The keratin antibodies, pankr1-8 (Progen Biotechnik, Heidelberg, Germany; 61006) and 79.14 (Developmental Studies Hybridoma Bank) were both used at a 1:10 dilution with standard immunofluorescence conditions. For phalloidin staining, embryos were fixed in 4% PFA, digested for 10 mins with 10 ug/ml proteinase K, then re-fixed in PFA and incubated in Alex Fluro 568 phalloidin (Invitrogen) overnight at 4° C. All whole mount in situ images were performed on a minimum of 12 embryos from replicate experiments and confocal images were from a minimum of 10 samples from replicate experiments.
Two replicate samples of twenty-four tails of 16–18 somite stage embryos from a cross of yap1+/-;wwtr1+/- adults were collected on seven separate occasions. Embryos from a cross of yap1+/-;wwtr1+/- adults were separated into mutants or siblings solely by phenotype. We collected the tails in a manner that randomized the collection order of the samples (mutant, siblings) and replicates (replicate 1, 2 and 3; a third one as backup). The collection order was not identical from one collection to the next. Tails are tissues of the whole tail from the tailbud to the third newest somite (S-III). RNA from two replicates of mutant and sibling samples were isolated using the miRNeasy micro Kit (Qiagen) and treated with in-column DNase digestion. Quality checks on RNA for library preparation were performed with a BioAnalyzer 2100 (Agilent Technologies) and LabChip Gx Touch 24 (Perkin-Elmer). A total of 300 ng of total RNA was used as input for Truseq Stranded mRNA library preparation. Sequencing was performed on the NextSeq 500 instrument (Illumina) using v1 chemistry, obtaining a minimum of 10–15M reads (75 bp paired-end) per library. Quality checks were performed on raw reads with FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc). Adapter sequences and poor quality base calls were trimmed from the raw reads with Trimmomatic v0.33 (parameters: LEADING:3 TRAILING:3 SLIDINGWINDOW:5:20 CROP:500 MINLEN:15)(Bolger et al., 2014). The processed reads were then aligned to Zv10 version of the zebrafish genome with STAR 2.4.0a (parameters: --outFilterScoreMinOverLread 0 --outFilterMatchNminOverLread 0 --outFilterMatchNmin 30 --outFilterMultimapNmax 999 --alignEndsProtrude 10 ConcordantPair) (Dobin et al., 2013). The resulting bam alignment file was passed through the Cufflinks suite v2.2.1 with standard parameters (Trapnell et al., 2012) to calculate the Fragments Per Kilobase of transcript per Million mapped reads (FPKM) for each gene as well as the statistics on the fold-changes. The summarized data can be found in Supplementary file 1 and the complete dataset is deposited on GEO with the accession number GSE102606.
Enrichment of genes for Gene Ontology (GO) terms and Kegg Pathways were performed using the Database for Annotation, Variation and Integrated Discovery (DAVID) web tool (Huang et al., 2009). Genes that are significantly downregulated (adjusted p-value<0.05, and sibling FPKM >8500) were defined as ‘Gene List’ in the input, and compared against all zebrafish genes. GO terms and Kegg Pathways were significantly enriched when the corrected P-value (Benjamini) <0.05.
We identified Yap1/Wwtr1 target genes from published datasets (Zhang et al., 2009; Zhao et al., 2008). The list of genes induced by YAP1 overexpression in MCF10A cells are found in Supplementary file 3 of Zhang et al., 2009 and Supplementary file 3 of Zhao et al. (2008). These two lists were overlapped to obtain a list of Yap1/Wwtr1 target genes. We then identified which of these genes are significantly downregulated in the tail buds of yap1;wwtr1 double mutants.
Tail bud samples were collected in the same manner as in the RNA-seq experiment. A total of 4 biological replicates (pool of 6 tail buds) from siblings and mutants were collected. Total RNA was isolated with the Direct-zol kit (Zymogen, R2061) following manufacturer’s instructions. RNA-cDNA reaction was performed with the Maxima First Strand cDNA Synthesis Kit for RT-qPCR (Thermo Scientific, K1671). A total of 80–130 ng of RNA per sample was used in the cDNA synthesis reaction. Primers for the qPCR reactions can be found in Supplementary file 4. The DyNAmo Flash SYBR Green qPCR Kit (Thermo Scientific, F415) was used for the qPCR reaction. The amplification cycle was: 95°C for 15 s; 65°C for 15 s; 72°C for 5 s. All genes of interest Cq values were normalized to rpl13 Cq values (ΔCq). ΔCqsibling – ΔCqmutant is the log2(Fold-change) for each gene of interest. The one-tailed two-sample t-test was performed to determine whether the ΔCq values are significantly different between siblings and mutants. P-values were corrected for multiple hypotheses testing with the FDR (false-discovery rate) method.
Formation and segmentation of the vertebrate body axisAnnual Review of Cell and Developmental Biology 29:1–26.https://doi.org/10.1146/annurev-cellbio-101011-155703
Trimmomatic: a flexible trimmer for Illumina sequence dataBioinformatics 30:2114–2120.https://doi.org/10.1093/bioinformatics/btu170
Biochemical identification and tissue-specific expression patterns of keratins in the zebrafish Danio rerioCell and Tissue Research 293:195–205.https://doi.org/10.1007/s004410051112
Yes-associated protein (Yap) is necessary for ciliogenesis and morphogenesis during pronephros development in zebrafish (Danio Rerio)International Journal of Biological Sciences 11:935–947.https://doi.org/10.7150/ijbs.11346
Neuromesodermal progenitors and the making of the spinal cordDevelopment 142:2864–2875.https://doi.org/10.1242/dev.119768
The evolutionary history of YAP and the hippo/YAP pathwayMolecular Biology and Evolution 28:2403–2417.https://doi.org/10.1093/molbev/msr065
Control of muscle fibre-type diversity during embryonic development: the zebrafish paradigmMechanisms of Development 130:447–457.https://doi.org/10.1016/j.mod.2013.06.001
yap is required for the development of brain, eyes, and neural crest in zebrafishBiochemical and Biophysical Research Communications 384:114–119.https://doi.org/10.1016/j.bbrc.2009.04.070
Mechanisms of convergence and extension by cell intercalationPhilosophical Transactions of the Royal Society B: Biological Sciences 355:897–922.https://doi.org/10.1098/rstb.2000.0626
Tales of tails (and trunks): forming the posterior body in vertebrate embryosCurrent Topics in Developmental Biology 116:517–536.https://doi.org/10.1016/bs.ctdb.2015.12.008
Multiple roles for Hippo pathway effector Wwtr1 in cardiac wal maturaion in zebrafishDevelopment.
Regulated tissue fluidity steers zebrafish body elongationDevelopment 140:573–582.https://doi.org/10.1242/dev.090381
The tissue mechanics of vertebrate body elongation and segmentationCurrent Opinion in Genetics & Development 32:106–111.https://doi.org/10.1016/j.gde.2015.02.005
Yap and Taz regulate retinal pigment epithelial cell fateDevelopment 142:3021–3032.https://doi.org/10.1242/dev.119008
The biology of YAP/TAZ: hippo signaling and beyondPhysiological Reviews 94:1287–1312.https://doi.org/10.1152/physrev.00005.2014
Fibronectins, their fibrillogenesis, and in vivo functionsCold Spring Harbor Perspectives in Biology 3:a005041.https://doi.org/10.1101/cshperspect.a005041
Laminin alpha5 is essential for the formation of the zebrafish finsDevelopmental Biology 311:369–382.https://doi.org/10.1016/j.ydbio.2007.08.034
Stem cells, signals and vertebrate body axis extensionDevelopment 136:1591–1604.https://doi.org/10.1242/dev.021246
Genome-wide association between YAP/TAZ/TEAD and AP-1 at enhancers drives oncogenic growthNature Cell Biology 17:1218–1227.https://doi.org/10.1038/ncb3216
TEAD transcription factors mediate the function of TAZ in cell growth and epithelial-mesenchymal transitionJournal of Biological Chemistry 284:13355–13362.https://doi.org/10.1074/jbc.M900843200
TEAD mediates YAP-dependent gene induction and growth controlGenes & Development 22:1962–1971.https://doi.org/10.1101/gad.1664408
Marianne BronnerReviewing Editor; California Institute of Technology, United States
In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.
Thank you for choosing to send your work, "Regulation of posterior body and ectodermal morphogenesis in zebrafish by localized Yap1 and Wwtr1", for consideration at eLife. Your article has been reviewed by three peer reviewers, and the evaluation has been overseen by a Senior/Reviewing Editor.
While all of the reviewers are enthusiastic about the study, they also agree that major revisions are necessary to bring this to publishable form. In particular, additional data are needed to support the conclusion that the defects are due to fibronectin misregulation. Genetic mosaic experiments are needed to determine whether the fin defect is due to a migration defect. Also, more detailed presentation of the expression profiling data is in order. Included below are the detailed comments of all three reviewers for your reference.
The authors present the analysis of double mutant zebrafish embryos for yap and wwrt1. This represents one of the first substantial analyses of the role of Yap in the embryonic development of zebrafish and represents an important contribution to the field.
Despite opting for a genome wide analysis of the genes that are mis-regulated upon mutation in yap-/-;wwtr1-/-mutants, only a very limited analysis of this dataset is performed in the manuscript. At the very least it would be of interest to look into the top genes that are upregulated as well as downregulated. Although, it is not clear from Table 1 whether the 'fold change' that is referred to is positive or negative. In addition, there are a number of other common targets for yap that have been found in other contexts (such as CTGF), it would be useful to see whether these are also modulated in the context of the zebrafish. I think that these analyses are necessary to support the conclusions made about notochord and ectodermal markers being the predominant sets of genes that we mis-regulated in the mutants.
In order to conclude about the role of the notochord phenotype in generating defects in posterior body development, it is necessary to include a more detailed comparison to the yap-/-;wwtr1-/- mutants. While the text states that the defects paralleled those of the yap-/-;wwtr1-/- mutant, it is not clear what aspects of the phenotype are the same or different. For example, the timing of somitogenesis and overall somite number, as well at the timing of axial truncation could easily be determined for both mutants and displayed alongside one another. As it stands, I am not sure that anything substantial is gained in this comparison. Furthermore, the fact that defects in yap-/-;wwtr1-/-mutants precede those in embryos where the notochord is absent suggest that there is an additional function of yap/wwtr1 at earlier stages that precede both its role in the notochord and in fin formation. The authors could make use of their heat-shock inducible lines in order to temporally separate these events.
Similarly, a more detailed spatiotemporal analysis of the onset of disruptions in the Fn deficient embryos in comparison to the yap-/-;wwtr1-/-mutants is required in order to conclude that this is a phenocopy.
The authors claim on the basis of their live imaging experiments that the cells forming the ventral fin fold are 'highly migratory' and that these behaviours are moderately reduced in the mutants. It is not clear from the one control movie shown in the paper whether the cells are migratory, as it seems that only very little cell rearrangement is occurring. While the cells are undergoing some cell shape changes, this is not quantified or correlated with the formation of the fin fold in time and space. In order to attribute specific cell behaviours to fin fold formation, it is necessary to quantify these behaviours and compare them to ectodermal cells in regions that do not form a fin, or at stages preceding fin formation. These same parameters can be then compared to the mutant situation. At the very least the authors should check cell velocity, directionality and alterations in cell shape in order to support their conclusions. Also, the number of videos that have been analysed in each condition should be detailed.
The study by Kimelman et al., addresses the timely question whether and how Yap signaling is involved in vertebrate body elongation. In recent years, Yap has been implicated in mechanotransduction. While a number of recent studies have examined the role of mechanics in body elongation, none have address the obvious question of how Yap signaling may operate during elongation. The results are surprising in that Yap signaling is required for normal zebrafish body elongation but that Yap's function is restricted to the notochord and epidermis. Moreover, Yap does not appear to function by regulated arhgaps, as suggested by studies in organoids, but rather perturbs the Fibronectin matrix. The current study is thorough, and includes the generation of zebrafish knockouts for Yap, its paralog wwtr1, and characterization of double mutants. Transcriptional profiling of dissected tailbuds is used to characterize changes in gene regulation. Heat shock dominant negative and constitutively active Yap transgenes provide corroborating evidence. Overall, the study is very interesting, the manuscript is clearly written and, in general, the conclusions are supported by the data presented. However, there are a few conclusions that need to be better substantiated, as detailed below.
Expression of the dn-Fn construct approximates the yap1;wwtr1 double mutant phenotype, but the authors should perform timelapse analysis of the dn-Fn fin fold morphogenesis to confirm that the phenotype is consistent at the level of cell behavior. The phenotypes could otherwise arise from slightly different mechanisms.
The proposition that adhesion of the ectoderm is reduced in the yap1;wwtr1 double mutant needs to be strengthened. The authors are correct that fixation/over-fixation may obscure this phenotype. It could be revealing to image live yap1;wwtr1 double mutants expressing a membrane GFP/ and nuclear RFP. In transverse sections, one would predict that there would be evident separation of the ectoderm from the dorsal somite.
The major question that comes to mind, particularly with respect to the ectoderm phenotype, is whether the phenotype is due to failure to create a normal Fibronectin matrix, which is then affecting cell migration, or whether the cells migrate abnormally and this affects the morphology of the Fibronectin matrix? Genetic mosaics could potentially address this question.
The authors confirm previous findings from Medaka indicating that Yap and Wwtr1 are essential for posterior axis elongation. Whereas, previously it has been suggested that axis elongation in Yap/Wwtr1 loss of function is caused by an inability to combat gravity these authors suggest a role for Yap and Wwtr1 in ectodermal and notochord morphogenesis. They show that Yap and Wwrt1 are expressed specifically in tailbud ectoderm and notochord. Moreover, through RNAseq analysis they find that YAP and Wwrt1 control the expression of exrg4b and wu:fc23c09 in the tailbud. In contrast, Yap and Wwrt1 do not control arhgap12a/27 expression as previously reported. The authors then investigate ectodermal morphogenesis in mutants and observe a loss of ectodermal thickening associated with dorsal fin formation. Through live imaging the authors note that, in WT embryos, ectodermal cells appear to migrate toward the midline. In contrast, mutant ectodermal cells do not accumulate in the midline despite behaving dynamically. As Fibronectin is an important factor in axis elongation and tissue adhesion the authors analyze FN deposition and find that FN is not maintained in newly forming somites where yap and wwrt1 is expressed. Finally, they indicate that epithelial integrity is, therefore, lost in mutants.
This article is appropriate for publication in eLife with major revisions.
Fin fold formation is thought to occur either through shape changes in static cells or movement of the whole epidermal sheet toward the midline. The authors attempt to clarify this controversy through live imaging. The authors find cells that enter the midline over time and determine that this occurs through cell migration. However, these data do not exclude that convergent extension like movements could be sufficient to move cells toward the midline. Here, the appearance of cells in the midline could be mistaken for cell migration when their appearance is actually the consequence of a cell rearrangement driven, perhaps, through t1 transitions. Their conclusions are therefore overstated and should be experimentally addressed for publication.
Although their data are highly correlative with a role for FN downstream of Yap and wwrt1, these data do not confirm this interaction. An epistasis experiment is required to determine whether FN is downstream of Yap and Wwrt1.https://doi.org/10.7554/eLife.31065.036
- David Kimelman
- Natalie L Smith
- Jason Kuan Han Lai
- Didier YR Stainier
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
We thank Young Kwon for suggesting (and providing) the Yap1 antibody, Scott Holley for sharing his Fibronectin staining protocol, Siddharth Ramakrishnan for providing medaka embryos, as well as Michelle Collins and Stefan Günther for help preparing the RNA-seq samples.
Animal experimentation: This study was performed in strict accordance with institutional (UW and MPG) and national ethical and animal welfare guidelines. All of the animals were handled according to approved institutional animal care protocols (Permission No. B2/1068 for DS and IACUC protocol 2387-02 for DK).
- Marianne Bronner, Reviewing Editor, California Institute of Technology, United States
© 2017, Kimelman et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.