Wnt11 acts on dermomyotome cells to guide epaxial myotome morphogenesis
Abstract
The dorsal axial muscles, or epaxial muscles, are a fundamental structure covering the spinal cord and vertebrae, as well as mobilizing the vertebrate trunk. To date, mechanisms underlying the morphogenetic process shaping the epaxial myotome are largely unknown. To address this, we used the medaka zic1/zic4-enhancer mutant Double anal fin (Da), which exhibits ventralized dorsal trunk structures resulting in impaired epaxial myotome morphology and incomplete coverage over the neural tube. In wild type, dorsal dermomyotome (DM) cells reduce their proliferative activity after somitogenesis. Subsequently, a subset of DM cells, which does not differentiate into the myotome population, begins to form unique large protrusions extending dorsally to guide the epaxial myotome dorsally. In Da, by contrast, DM cells maintain the high proliferative activity and mainly form small protrusions. By combining RNA- and ChIP-sequencing analyses, we revealed direct targets of Zic1, which are specifically expressed in dorsal somites and involved in various aspects of development, such as cell migration, extracellular matrix organization, and cell-cell communication. Among these, we identified wnt11 as a crucial factor regulating both cell proliferation and protrusive activity of DM cells. We propose that dorsal extension of the epaxial myotome is guided by a non-myogenic subpopulation of DM cells and that wnt11 empowers the DM cells to drive the coverage of the neural tube by the epaxial myotome.
Editor's evaluation
This study addresses an interesting and underexplored question in developmental biology, specifically cell migration and muscle development. It builds upon prior analysis of the medaka Double anal fin (Da) mutants by using detailed bioinformatic and time-lapse analysis to explain dorsal somite extension. The authors show that dorsal muscle morphogenesis is actively guided by dorsal dermomyotome cells, rather than being passively shaped by physical constraints alone. Looking downstream of Da, they show that Wnt signaling is central to dorsal extension of the epaxial myotome and propose that similar functions may shape the dorsal musculature across vertebrates.
https://doi.org/10.7554/eLife.71845.sa0Introduction
Active locomotion, which is powered by skeletal muscles in vertebrates, is critical for animals to survive. Vertebrate skeletal muscles consist of axial muscles (head, trunk, and tail muscles) and appendicular muscles (limb muscles). Axial muscles first arose in the chordates to stabilize and enable side-to-side movement of the body axis. In jawed vertebrates, the subdivision of axial muscles into epaxial (dorsal) and hypaxial (ventral) muscles led to an increased range of movement: dorsoventral undulation in fish and lateral movements in terrestrial vertebrates (Goodrich, 2004; Romer and Parsons, 1986; Fetcho, 1987; Sefton and Kardon, 2019). Among these muscles, epaxial muscles are characterized by their unique anatomical structure, which extends dorsally and surrounds the vertebrae. This morphology also ensures mechanical support and protection of the vertebrae and the spinal cord inside. While we have a detailed understanding of how the myotome, precursors of epaxial and hypaxial muscles, differentiates from a somitic compartment called the dermomyotome (DM) (Kalcheim et al., 1999; Gros et al., 2009; Hollway and Currie, 2005), we only begin to understand the cellular and molecular mechanisms of the subsequent morphogenetic processes generating the epaxial muscles. Previous studies in rats and mice suggested that myocytes of the epaxial myotome do not actively migrate dorsally but are guided by external forces (Deries et al., 2010; Deries et al., 2012). However, what exerts such forces to drive extension of epaxial myotome is still unclear.
Fish have been extensively utilized to study myotome development thanks to the transparency of their bodies throughout embryonic development (Wolff et al., 2003; Stellabotte et al., 2007; Nguyen et al., 2017; Ganassi, 2018). Additionally, their epaxial trunk muscles have a simple structure consisting of only one anatomical unit (reviewed in Sefton and Kardon, 2019). Like other vertebrates, fish myotomes, on either side of the neural tube, extend dorsally after somite differentiation and eventually cover the neural tube by the end of embryonic development (Figure 1A). The spontaneous medaka (Oryzias latipes) mutant Double anal fin (Da) displays a particular epaxial myotome morphology, in which the dorsal ends of the left and right epaxial myotome fail to extend sufficiently and thus do not cover the neural tube at the end of embryonic development. Previous studies demonstrated that the dorsal trunk region of the Da mutant is transformed into the ventral one, including not only the myotome but also the body shape, skeletal elements, pigmentation, and fin morphology (Figure 1B and D; Ishikawa, 1990; Ohtsuka et al., 2004). Given the unique morphological features, the medaka Da mutant is an excellent model to study the morphogenesis of epaxial myotome. Genetic analysis of the Da mutant revealed that this phenotype is due to a dramatic reduction of the expression of the transcription factors zic1 and zic4 in the dorsal somites (Figure 1C and E), and identified zic1/zic4 as master regulators of trunk dorsalization (Ohtsuka et al., 2004; Kawanishi et al., 2013). The down-regulation of zic1/zic4 specifically in the dorsal somites is caused by the insertion of a large transposon, disrupting the dorsal somite enhancer of zic1/zic4 (Moriyama et al., 2012; Inoue et al., 2017). While the function and downstream targets of Zic1 and Zic4 have been studied in nervous system development and somitogenesis (Li, 2006; Pan et al., 2011; Himeda et al., 2013; Hong and Saint-Jeannet, 2017; Aruga and Millen, 2018), the molecular mechanism of how these Zic genes control dorsal trunk morphogenesis has not been investigated so far.
Here, we describe the morphogenetic process of the formation of the epaxial myotome of the back, which we termed ‘dorsal somite extension’ (Figure 1A). By in vivo time-lapse imaging we uncovered its cellular dynamics; during dorsal somite extension, dorsal DM cells reduce their proliferative activity and subsequently form unique large protrusions extending dorsally, guiding the epaxial myotome towards the top of the neural tube. We also found that these DM cells form a subpopulation that gives rise to non-myotomal cell lineages during embryonic development. In the Da mutant, by contrast, DM cells keep their high proliferative activity and mainly form small protrusions. Mechanistically, we identified a Zic1 downstream-target gene, wnt11 (former wnt11r Postlethwait et al., 2019), as a crucial factor for dorsal somite extension, and demonstrated that Wnt11 regulates cellular behavior of dorsal DM cells by promoting protrusion formation and negatively regulating proliferation. We thus propose an unprecedented process of epaxial myotome morphogenesis driven by a non-myogenic population of DM cells during embryogenesis.
Results
Our previous study showed that zic1 and zic4 expression starts at embryonic stages and persists throughout life (Kawanishi et al., 2013). Phenotypic analysis of homozygous adult Da mutants implies long-term participation of Zic-downstream genes in the formation of dorsal musculatures, which eventually affects the external appearance of the fish adult trunk. Here, we examined the initial phase of this long-term dorsalization process. In the following of the study, we will focus on zic1, since zic1 and zic4 are expressed in an identical fashion with overlapping functions in trunk dorsalization of medaka, and zic4 is expressed slightly weaker than zic1 (Moriyama et al., 2012; Kawanishi et al., 2013).
The dorsal myotome of the Da mutant fails to cover the neural tube
In wild-type (Wt) medaka, the dorsal ends of the myotomes first came in contact at 7 days post fertilization (dpf, stage 37) and formed the tight, thick myotome layer covering the neural tube at the end of embryonic development (9 dpf, stage 39) (Figure 1F–G, Figure 1—figure supplement 1A-F). In the ventralized Da mutant, however, the dorsal ends of the myotomes did not extend sufficiently and failed to cover the neural tube at the end of embryonic development (Figure 1H–I, Figure 1—figure supplement 1G-L). The ends of the ventralized dorsal myotome in the mutant displayed a round shape (not pointed as found in Wt) which resembled the morphology of the ventral myotome.
We wondered if there are other morphological differences between Wt dorsal myotome and the ventralized dorsal myotome of the Da mutant. Indeed, the cross-sectional area in the Da mutant was significantly larger compared to Wt (Figure 1L). Possible explanations for a larger cross-sectional area in the Da mutant could be a larger myofiber diameter or a higher number of myofibers, which make up the myotome. We measured the diameter of dorsal myotome muscle fibers in Wt and Da mutant embryos but could not observe a difference, suggesting that dorsal myotome of the Da mutant has a higher number of myotomal cells (Figure 1—figure supplement 1M).
Proliferative activity of the dorsal DM cells is enhanced in the Da mutant
In fish, as in other vertebrates, the DM gives rise to muscle precursor cells, which ultimately differentiate into myofibers. In medaka the DM is a one cell-thick, Pax3/7-positive cell layer encompassing the myotome (Figure 2A–B’’; Hollway et al., 2007; Abe et al., 2019). A high proliferative activity of the dorsal DM could explain a larger cross-sectional area of the dorsal myotome of the Da mutant. To test this, we performed immunohistochemistry against the mitotic marker phospho-histone H3 (pH3) and the DM marker Pax3/7 on Wt and Da embryos (Figure 2C–E). In both Wt and Da embryos, pH3-positive cells were randomly distributed in the dorsal DM without obvious bias (Figure 2C–D). At the 12-somite stage (12 ss, 1.7 dpf stage 23), when zic1 expression in the somites becomes restricted to the dorsal region (Kawanishi et al., 2013), the number of pH3-positive DM cells per dorsal somite was not significantly different between Wt and Da. Remarkably, from 16 ss (1.8 dpf, stage 24) onwards, the number of pH3-positive cells became reduced in the Wt, whereas in Da, no such reduction was observed (Figure 2E). At 35 ss (3.4 dpf, stage 30), pH3-positive cells increased both in the Wt and Da, but the mutant DM cells were more proliferative (Figure 2E). Immunohistochemistry against another proliferation marker PCNA confirmed these findings (Figure 2—figure supplement 1A-C). The number of pH3-positive cells in the ventral DM was not significantly different in Wt and Da embryos at 22 ss (2.75 dpf, stage 26) (Figure 2—figure supplement 2A). These results suggest that zic1 reduces proliferative activity of the DM, which becomes evident following the confinement of its expression to the dorsal somite region.
Wt dorsal DM cells form numerous large, motile protrusions at the onset of dorsal somite extension
The epaxial myotome, on either side of the neural tube, extends dorsally to cover the neural tube by the end of embryonic development. In the Da mutant, the myotome is unable to cover the neural tube despite increased dorsal myotome growth at the hatching stage (Figures 1L and 2). This suggest that additionally to physical extension an active process might support the dorsal movement of somites. To examine the behavior of zic1-positive cells underlying this dorsal somite extension, we performed in vivo time-lapse imaging of dorsal somites using the transgenic line Tg(zic1:GFP,zic4:DsRed), which expresses GFP under the control of the zic1 promoter and enhancers to visualize the dorsal somitic cells (Kawanishi et al., 2013) (hereafter called Tg(zic1:GFP) since the DsRed fluorescence was negligible in the following analyses). Intriguingly, around 22 ss and onwards, cells at the tip of the dorsal somites started to form numerous large protrusions extending dorsally towards the top of the neural tube (Figure 3A, Video 1). We defined the beginning of protrusion formation as the onset of dorsal somite extension. Close-up views of the time-lapse images (Figure 3A, Video 1) showed that these protrusions were motile and dynamically formed new branches at their dorsal tips (Figure 3J, Figure 3—figure supplement 1D). Immunohistochemistry revealed that the protrusion-forming cells belong to the DM (Figure 3—figure supplement 1A-A'').
To characterize the protrusions, we classified them according to their length into small ( < 8 μm, Figure 3C arrowheads) and large ( ≥ 8 μm, Figure 3C arrow) protrusions (Figure 3—figure supplement 1B-C). Based on their shape, we reasoned that the small protrusions correspond to lamellipodia (Figure 3C, arrowheads), while the large protrusions appeared more complex. To investigate the nature of large protrusions, we injected Actin-Chromobody GFP mRNA to visualize the actin skeleton. The large protrusions were found to exhibit a complex architecture consisting of a lamellipodia-like core structure (Figure 3D, bracket) with additional multiple bundles of filopodia (protrusions with linear arranged actin filaments) branching out from their dorsal tips (Figure 3D, arrows, summarized in Figure 3E).
Interestingly, time-lapse in vivo imaging of Tg(zic1:GFP);Da showed that in the Da background, protrusions started to form later (Figure 3B, Video 2) and the number of large protrusions and protrusions in total was significantly lower than in Wt (Figure 3I–I’’). In addition, protrusions were transient and mostly failed to form new branches at their dorsal tips (Figure 3K, Figure 3—figure supplement 1E). While no difference in the actin skeleton of small protrusions could be observed, filament bundles branching out from large protrusions of Da DM cells contained fewer and shorter filopodia compared to Wt (Figure 3G, arrowheads, summarized in Figure 3H). These results indicate that the protrusive activity, especially the ability to form large protrusions, is significantly reduced in the Da mutant. The large protrusions of the dorsal DM cells continuously appeared at later stages of dorsal somite extension, too (Figure 3—figure supplement 2).
To investigate the role of large protrusions during the onset of dorsal somite extension, we inhibited filopodia formation using ML141 (Figure 3L and M). ML141 specifically inhibits Cdc42/Rac1 GTPases, which are critical for filopodia formation (Hong et al., 2013; Fantin, 2015). Intriguingly, dorsal somites of embryos treated with ML141 extended significantly less dorsally compared to control embryos treated with DMSO (Figure 3N). Taken together, our data suggest that the unique large protrusions of the dorsal DM are involved in guiding the epaxial myotome dorsally and zic1 might promote this function.
DM cells delaminate and accumulate between opposing somites during the late phase of dorsal somite extension
We continued to trace the behavior of the DM tip cells until the somites reach the top of the neural tube. Indeed, in vivo imaging of Tg(zic1:GFP) revealed that dorsal DM cells continued to form protrusions, and additionally, some of them delaminated to become Zic1-positive mesenchymal cells accumulating in the space between the dorsal ends of the left and the right somites (Figure 4A, star-shaped cells, arrowheads) from 4.5 dpf (stage 33) onwards. This is consistent with previous observations of strongly zic1 expressing mesenchymal cells in Wt at late embryonic stages (Ohtsuka et al., 2004). As dorsal somite extension proceeded, the number of these mesenchymal cells increased, filling the space between the two myotomes (Figure 4A–C, arrowheads, Videos 3–4, arrowheads indicate representative mesenchymal cells). These mesenchymal cells formed protrusions towards neighboring mesenchymal cells and DM cells at the tip of somites, creating a dense cellular network between the dorsal ends of the somites. Mosaic cell-labelling demonstrated that the mesenchymal cells originated from the DM (Figure 4—figure supplement 1A-F, arrowheads). Interestingly, while mesenchymal cells dynamically formed protrusions, they showed no extensive migratory behavior and were rather stationary (Videos 3–4, Figure 4—figure supplement 2). This could suggest that these protrusions fulfill a non-migratory function. When the opposing somites came in contact with each other at 8 dpf (stage 38, 1 day before hatching), the mesenchymal cells tended to attach to the nearest DM cells at the tip of the somite, bridging the gap between the left and right DM cells (Figure 4I, Figure 4—figure supplement 1M-M'''').
In Da mutants, mesenchymal DM cells were also detected in the space between the two myotomes, but the timing of their appearance was delayed, that is between 5 dpf (stage 34) and 5.5 dpf (stage 35) (Figure 4D–F) (4.5 dpf in Wt). Additionally, Da mesenchymal cells exhibited a rounder morphology and formed significantly fewer protrusions compared to Wt while the number of mesenchymal cells was not largely affected (Figure 4G and H, Figure 4—figure supplement 1G-L).
To examine the function of mesenchymal DM cells, we ablated these cells between the 10th somites of 5.5 dpf Tg(zic1:GFP) embryos using an UV laser (Figure 4J–K’, Figure 4—figure supplement 3A-B'). Intriguingly, the distance between the left and the right tips of the 10th somites increased after ablating the mesenchymal DM cells while the neighboring 9th somites continued to shorten the gap (Figure 4L). This suggests that mesenchymal DM cells hold the left and the right somite together to promote the dorsal somite extension at this late phase of myotome development. After 16 hr post ablation, dorsal extension of the 10th somites eventually resumed as the mesenchymal DM cells regenerated at the ablation site (Figure 4L).
Collectively, dorsal DM cells and the mesenchymal cells derived from them seem to actively participate in the entire process of dorsal somite extension, from its onset to neural tube coverage at the end.
Finally, we examined the fate of the dorsal DM cells that derive the mesenchymal cells. We employed a photoconversion technique mediated by a photoconvertible protein KikGR to specifically label a single dorsal DM cell exhibiting protrusions at the tip of a 10th and 20th somites during dorsal somite extension (Figure 4M, arrowhead) and tracked them until the hatching stage. We found that the labeled dorsal DM cells eventually differentiated into cells surrounding blood vessels (likely mural cells, n = 8/9 and 4/8 for 10th and 20th somites, respectively; Figure 4N and N’, arrowheads) and mesenchymal cells in the dorsal fin fold (n = 4/8 for 20th somites; Figure 4O and O’, arrowheads); however, we did not observe any labeled axial muscles (n = 0/9 and 0/8 for 10th and 20th somites, respectively). This suggests that the dorsal DM cells at the tip of somites, guiding myotome extension, do not become myotomal cells themselves at least until the end of embryonic development (hatching stage, 9 dpf).
Zic1 regulates the expression of dorsal-specific genes during somite differentiation
We then addressed the molecular machinery controlling the dorsal somite extension investigated above. Since somite extension is impaired in the zic1-enhancer mutant Da (Moriyama et al., 2012; Kawanishi et al., 2013), we reasoned that downstream genes of Zic1 are regulators of this process.
First, we identified genes which are specifically expressed in the dorsal somites. For this, we collected somites of the transgenic line Tg(zic1:GFP) by removing other tissues after pancreatin treatment, FACS sorted them into the dorsal (GFP+) and ventral (GFP-) somitic cells including the DM cells, and performed RNA-seq and ATAC-seq on both cell populations (Figure 5A, Figure 5—figure supplement 1A). The RNA-seq identified 1,418 differentially expressed genes. Among them 694 genes showed higher expression in the dorsal somites (termed hereafter dorsal-high genes), and 724 genes showed higher expression in the ventral somites (termed hereafter dorsal-low genes) (Figure 5C). We confirmed that zic1 and zic4 were found among the dorsal-high genes (Figure 5C, Supplementary file 1).
Next, we identified potential direct Zic1-target genes by investigating Zic1 binding sites. Since there were no suitable antibodies available to perform ChIP-seq against medaka Zic1, we created a transgenic line expressing a Myc tagged zic1 in the Da background under the control of zic1 promoter and enhancers (Tg(zic1:zic1-Myc,zic4:DsRed);Da, called Tg(zic1:zic1-Myc);Da hereafter). This transgene (zic1:zic1-Myc) efficiently rescued the ventralized phenotype (dorsal and anal fin shape, pigmentation pattern, and body shape) of the Da mutant (Figure 5—figure supplement 1B-C), verifying the full functionality of the tagged Zic1 protein. Somites of the transgenic line Tg(zic1:zic1-Myc);Da were dissected and subjected to ChIP-seq using antibodies against Myc (Figure 5B). Since somites contain a low number of cells, we applied ChIPmentation (Schmidl et al., 2015) to identify genome-wide Zic1 binding sites. From two biological replicates, 5247 reliable ChIP peaks were identified, and we confirmed that the enriched DNA motif among these peaks was showing high similarity with previously identified binding motifs of ZIC family proteins (Figure 5—figure supplement 2A-B). Then, we associated each Zic1 peak to the nearest gene within 50 kb and identified 3232 genes as Zic1 target genes. By comparing the ChIPmentation with RNA-seq data, we found that Zic1 target genes are overrepresented in dorsal-high and dorsal-low genes. While 47% of the dorsal-high genes and 27% of the dorsal-low genes were found to be Zic1 target genes, only 12% of genes which showed no differential expression in dorsal or ventral somites were potential Zic1 downstream target genes (Figure 5D, Supplementary file 1). This suggests that Zic1 can function as transcriptional activator and repressor of versatile genes, but the former role seems to be dominant.
wnt11 is a direct downstream target of Zic1 and down-regulated in the dorsal somites of the Da mutant
To identify potential regulators of dorsal somite extension, we further investigated the differentially expressed Zic1 target genes. Gene Ontology (GO) analysis indicated that both dorsal-high and dorsal-low gene groups, regardless of whether they are Zic1 targets or non-Zic1 targets, were significantly enriched in development related GO terms (Figure 6A, Figure 6—figure supplement 1A, Supplementary files 2 and 3). This indicates that Zic1 regulates a number of developmental genes both directly and indirectly. These results are consistent with the fact that Zic1 regulates various dorsal-specific morphologies of somite-derivatives (Kawanishi et al., 2013).
In the dorsal-high Zic1 targets, GO terms related to cell migration showed higher enrichment (e.g. ‘chemotaxis’ (p = 5.62E-10), ‘locomotion’ (6.65E-15), ‘ameboidal-type cell migration’ (p = 9.32E-10)) than dorsal-low genes or non-Zic1 target genes (Supplementary file 3). Interestingly, Wnt signaling pathway genes (e.g.: axin2, wnt11, sp5, lrp5, fzd10, prickle1a) and semaphorin-plexin signaling pathway genes (e.g.: sema3a, sema3c, sema3g, plxna1, plxna2, plxnb2, plxnb3, nrp2) were included in these gene groups. Indeed, the GO terms ‘Wnt signaling pathway’ and ‘Semaphorin-plexin signaling pathway’ themselves were also significantly enriched (p = 4.45E-7 and 1.94E-8, respectively) in dorsal-high Zic1 target genes (Supplementary file 2). This suggests that Zic1 directly regulates Wnt pathway, semaphorin and plexin genes, possibly to regulate cell movement in the dorsal somites. We also noticed that the terms ‘extracellular matrix organization’ (p = 7.74E-7; e.g.: adamts20, fbln1), ‘cell communication’ (p = 1.95E-5; e.g.: efna5, epha3) are enriched, suggesting that these genes also affect dorsal somite cell behavior.
Wnt signaling pathway components also ranked high among the dorsal-high Zic1 target genes by pathway enrichment analysis (Figure 6B, Figure 6—figure supplement 1B). Among them, we focused on wnt11 for further analyses due to the following reasons: First, wnt11 was one of the most differentially expressed genes in dorsal somites (Figure 5C, Supplementary file 1). Second, previous studies implicated Wnt11 in protrusion formation and cell migration (Ulrich et al., 2003; De Calisto et al., 2005; Garriock and Krieg, 2007; Matthews et al., 2008). This is particularly interesting since DM cells also exhibit protrusions and migration activity during dorsal somite extension, which are defective in Da mutants.
At the wnt11 locus, peaks of the Zic1-ChIP overlapped with intergenic open chromatin regions downstream of wnt11. These sites were more accessible in dorsal somites than in ventral somites, suggesting that Zic1 regulates wnt11 directly (Figure 6C). Additionally, in situ hybridizations against wnt11 performed on Wt and Da embryos indicated that wnt11 expression is significantly reduced at the dorsal tip of Da dorsal somites including the DM (Figure 6—figure supplement 1C-D). Already when the zic1 expression becomes restricted in the dorsal somites (12 ss) and with proceeding development when zic1 expression gets further restricted to the most dorsal part of the somites and mesenchymal DM cells (Ohtsuka et al., 2004), the expression of wnt11 in the dorsal somites of the Da mutant was reduced compared to Wt (Figure 6D–F, Figure 6—figure supplement 1C-F').
Taken together, we identified wnt11 as promising downstream target gene of Zic1, which could be a novel somite dorsalization factor and play a role in dorsal somite extension.
Wnt11 regulates cell behavior of dorsal DM
Our previous RNA-seq dataset revealed that the expression of wnt11 starts before gastrulation and increases as development proceeds (Figure 7—figure supplement 1A; Nakamura et al., 2021). Injection of single guide RNAs (sgRNAs) against wnt11 to knock out the gene revealed that the resultant genetically mosaic F0 embryos displayed delayed epiboly movement and subsequent morphological defects, including shorter body axis and impaired trunk development (Figure 7—figure supplement 2), implying the difficulty of assessment of the wnt11 function at later stages using the wnt11 mutant. To circumvent the problem, we took two different approaches to knock-down wnt11 during dorsal somite extension. Firstly, we used a wnt11 anti-sense morpholino (wnt11 MO) and determined a concentration which resulted in maximal knock-down effects during dorsal somite extension with minimal gastrulation phenotypes (see Materials and methods). Secondly, we performed temporally controlled knock-down of wnt11 after gastrulation using photo-cleavable Photo-Morpholinos (PhotoMOs) (Tallafuss, 2012; Figure 7A). Strikingly, in wnt11 Photo-Morphants with severe phenotypes the dorsal myotome failed to cover the neural tube (n = 7, Figure 7B–C), a phenotype similar to the Da mutant myotome. Additionally, the myofibers of the Photo-Morphants at 9 dpf were shorter and less organized. This phenotype is consistent with a previously reported role of Wnt11 during early myogenesis, where it regulates the elongation and orientation of myoblasts (Gros et al., 2009).
Next, we investigated the proliferative activity of dorsal DM cells in wnt11 morphants and found that knock-down of wnt11 (either by conventional MO or by PhotoMO) induced significantly more pH3-positive cells per somite, compared to control embryos at 22 ss (Figure 7D). These findings are similar to the observations previously made in the dorsal DM of Da mutants, and suggest that Wnt11 is negatively regulating proliferation of dorsal DM cells.
We then explored whether the onset of dorsal somite extension in wnt11 morphants is similarly impaired as in the Da by time-lapse in vivo imaging. Remarkably, we observed a protrusion formation behavior of wnt11 morphant DM cells similar to Da DM cells, namely delayed onset of protrusion formation and the formation of fewer, shorter protrusions (Figure 7E, Video 5). Quantification of protrusions indicated that wnt11 morphants have significantly fewer large protrusions and protrusions in total, compared to control embryos (Figure 7F–H’’).
Overall, knock-down of the Zic1 target gene wnt11 recapitulated the phenotype of Da DM cells (Figure 3G–G’’), suggesting the crucial role of Wnt11 in regulating protrusion formation of DM cells.
To further confirm the importance of Wnt11 during dorsal somite extension, we performed rescue experiments in Da embryos. At 18 ss (2.1 dpf, stage 25), we injected a mix of human recombinant Wnt11 (hrWnt11) protein (or BSA in the control group) and Dextran Rhodamine onto the top of the 10th somite of Tg(zic1:GFP);Da embryos (Figure 8A–B’’). Strikingly, the dorsal DM of Da embryos injected with Wnt11 formed significantly more large protrusions and more protrusions in total compared to Da embryos injected with BSA only (Figure 8C–C’’). This indicates that the ventralized Da DM protrusion phenotype can be partially rescued by Wnt11 protein injections and further emphasizes the importance of Wnt11 in somite dorsalization and during dorsal somite extension.
Wnt11 acts through the Wnt/Ca2+ signaling pathway at the onset of somite extension
Finally, we investigated through which signaling pathway the non-canonical Wnt11 acts. In Xenopus, Wnt11 acts through the Wnt/Ca2+ pathway regulating the migration of cells from the dorsal somite and the neural crest into the dorsal fin fold (Garriock and Krieg, 2007). To examine whether this signaling pathway also plays a role during dorsal somite extension, we inhibited the Wnt/Ca2+ signaling pathway using KN-93 in Tg(zic1:GFP) embryos from 4 ss (1.3 dpf, stage 20)–22 ss (Figure 9A). KN-93 specifically inhibits CaMKII, a component of the Wnt/Ca2+ pathway, by binding to it and preventing its phosphorylation, thus keeping it in an inactivated state (Sumi et al., 1991; Tombes et al., 1995; Wu and Cline, 1998; Garriock and Krieg, 2007; Rothschild et al., 2013). Embryos treated with KN-93 showed a higher number of pH3-positive dorsal DM cells per somite, compared to embryos in the control group (Figure 9B). Furthermore, DM cells at the tip of the dorsal somite of embryos treated with KN-93 formed significantly fewer large and fewer protrusions in total, compared to embryos in the control group (Figure 9C–C’’), although the effect was less significant compared to the wnt11 morphants (Figure 7).
From these results, we suggested that Wnt11 potentially acts through the Wnt/Ca2+ signaling pathway at the onset of somite extension. Wnt/Ca2+ signaling pathway is known to regulate actin polymerization (Choi and Han, 2002; Kohn and Moon, 2005) which could explain the dynamic protrusive activity of the DM cells (Figure 3), emphasizing the importance of Wnt11 during epaxial myotome morphogenesis.
Discussion
Here, we elucidated the key developmental process underlying the epaxial myotome morphogenesis to cover the neural tube in the teleost, medaka. DM cells have been known to serve as a progenitor pool for myotomal and dermal cells (Ben-Yair and Kalcheim, 2005; Hollway et al., 2007; Stellabotte and Devoto, 2007). In our study, we showed that DM cells at the tip of the dorsal somite form unique large protrusions, guiding the myotome towards the top of the neural tube. Furthermore, these DM cells do not give rise to myotomal cells during embryonic development. We found that zic1 and its direct downstream target wnt11 are crucial factors for this morphogenetic process. We thus revealed a novel role of non-myogenic DM cells during epaxial myotome morphogenesis. Furthermore, to our best knowledge, this work demonstrates for the first time that the neural tube coverage by myotomes is driven by active cell movement, rather than passively by muscle growth. Consistent with this, in the absence of Zic1 activity, the myotome is unable to cover the neural tube at the hatching stage, despite increased dorsal myotome growth (Figures 1L and 2), supporting the involvement of the active cellular mechanism in this process. The model we propose for dorsal somite extension is shown in Figure 10.
In the present study, we described the characteristic behavior of dorsal DM cells and propose their guiding role for the epaxial myotome moving towards the top of the neural tube. During dorsal somite extension, DM cells at the tip of the dorsal somites form large, motile protrusions, which contain multiple bundles of filopodia dynamically branching out from the tip of the protrusions. These long protrusions could be beneficial for invading the restricted open space between the neural tube and the ectoderm. Additionally, previous studies showed that in migrating mesenchymal cells, the formation of lamellipodia is associated with higher migratory speed, whereas filopodia play an exploratory role and are associated with high directionality (Leithner et al., 2016; Innocenti, 2018). Large protrusions of wild-type DM cells consist of a lamellipodia-like core and multiple filopodia, which could account for fast dorsal somite extension with high accuracy. The detailed structure and function of the large protrusions needs to be further investigated in the future studies.
In zebrafish, hypaxial muscle precursors delaminate from the ventral tip of the DM and migrate collectively in a cell stream to the prospective pectoral fin bud region while the cells at the leading edge form long filopodia (Haines et al., 2004; Talbot et al., 2019). Thus, DM cells, regardless of whether dorsal or ventral, seem to have the potential to form protrusions and to exhibit migratory behavior. However, regulation of migratory behavior could be different as the migration of hypaxial muscle precursors is mediated by the receptor tyrosine kinase Met (Haines et al., 2004) which is not expressed in zebrafish dorsal DM cells.
At later stages of dorsal somite extension, we observed that DM cells delaminate from the tip of the dorsal somites and progressively occupy the space between the opposing dorsal somites. These mesenchymal DM cells actively form protrusions towards neighboring mesenchymal DM cells and DM cells of the somites, thus forming a dense cellular network on top of the neural tube. This could provide a communication platform for the opposing somites to meet at the right position, exactly on top of the neural tube. Indeed, ablation of mesenchymal DM cells resulted in transient opening of the gap between the somites, suggesting that these cells at least exert a tension to hold the positions of the dorsal tips of the left and right somite together during neural tube coverage. Interestingly, our lineage tracing experiment revealed that the dorsal DM cells at the tip of somites differentiated into cells of non-myotomal lineages including mural cells of blood vessels and mesenchymal cells in the dorsal fin fold. These observations are consistent with previous studies reporting their somitic origins in zebrafish and Xenopus (Garriock and Krieg, 2007; Lee et al., 2013; Ando et al., 2016). However, our data also showed that these DM cells do not appear to contribute to the myotome at the hatching stage. We thus propose that the neural tube coverage is driven by a specific population of dorsal DM cells, which do not contribute to myotomal cells during embryonic development.
The ability of the motile DM cells to guide the underlying myotome dorsally assures complete coverage of the neural tube by the epaxial myotomes. Considering the fact that DM cells are in close cell-cell contact with the underlying myotome, it is plausible that they are biomechanically coupled with myotomal cells, facilitating dorsal myotome extension. A similar mechanism has been investigated in mouse neural tube formation. During mouse neural tube closure, the surface ectoderm, overlying the neuroepithelium, forms cell protrusions towards the midline. Disruption of membrane ruffles, a form of lamellipodia, in the surface ectoderm results in incomplete fusion of the subjacent neuroepithelium (Rolo et al., 2016).
Our data suggest that in medaka dorsal somites Wnt11 exerts its effect through promotion of protrusion formation and down-regulation of cell proliferation in the dorsal DM. Regarding the protrusion formation, previous studies of migrating mesendodermal cells in the zebrafish mutant silberblick and neural crest cells in Xenopus embryos consistently reported that Wnt11 is involved in the oriented elongation and stabilization of protrusions (Ulrich et al., 2003; De Calisto et al., 2005; Matthews et al., 2008). Similarly, myocardial cells in the heart of Wnt11 mutant mice also form fewer protrusions (Zhou et al., 2007). Thus, regulation of cell protrusions is a conserved function of Wnt11, observed in various developmental processes. However, the relationship between Wnt11 and cell proliferation could be context-dependent, as it negatively regulates proliferation in mouse neonatal hearts (Touma, 2017), while it promotes cell proliferation in mouse intestinal epithelial cell culture (Ouko et al., 2004). Hence, the function of Wnt11 in medaka somites is unique in that it promotes tissue elongation by regulating a balance between proliferative and migrative activity.
By inhibiting the Wnt/Ca2+ signaling pathway we showed that Wnt11 probably acts through this non-canonical Wnt signaling pathway during dorsal somite extension. Likewise, previous studies in Xenopus have shown that Wnt11 acts through this pathway during cell migration into the dorsal fin fold. There it regulates epithelial–mesenchymal transition in a distinct dorsal somite cell population which, together with a population of neural crest cells, contribute to the mesenchyme of the dorsal fin fold (Garriock and Krieg, 2007). Furthermore, during convergent extension in vertebrate gastrulation, Wnt/Ca2+ signaling pathway can regulate cell adhesion by promoting actin polymerization (Choi and Han, 2002; Kohn and Moon, 2005). Whether Wnt11 functions cell-autonomously or not is context dependent (De Calisto et al., 2005; Garriock and Krieg, 2007; Gros et al., 2009). In our study, administration of exogenous Wnt11 protein rescued the protrusive behavior of the surrounding DM cells, suggesting that Wnt11 can function cell non-autonomously during dorsal somite extension.
Importantly, our RNA-seq and ChIPmentation analyses revealed that Zic1 has diverse downstream target genes including various developmental genes. This suggests that the dorsal myotome is established via pleiotropic actions of Zic1; Wnt11 may not be a sole factor for dorsal somite extension, although it was suggested to be crucial in the present study. The semaphorin-plexin pathway may play a role, since previous studies suggested, besides its implication in axon guidance and neural cell migration, a role in non-neural cell migration (Alto and Terman, 2017). Furthermore, cell migration and tissue deformation are often linked with the extracellular matrix (ECM), which provides guiding or restraining cues influencing cell movements. Previous studies showed that during mouse epaxial myotome development, ECM composition changes dynamically, which is tightly accompanied by epaxial muscle morphogenesis; while the laminin content in the ECM decreases, increasing tenascin and stable fibronectin contents potentially promotes the alignment of myofibers and their final organization (Deries et al., 2012). Furthermore, recent studies suggested that cells can remodel the surrounding ECM and thereby increase their own motility. For example, during zebrafish gastrulation, the metalloproteinase mmp14 is expressed by migrating endoderm cells and degrades laminin and fibronectin, components of the ECM (Hu et al., 2018). In this context, of particular importance is our identification of Adamts20, encoding a proteoglycanase, as a dorsal-high Zic1-target gene in the somites (Supplementary file 1). Since Adamts20 is known to play a pivotal role in embryonic melanoblast migration by remodeling the dermal ECM (Rao et al., 2003; Silver et al., 2008), Adamnts20 could thus facilitate migration of Wnt11-expressing DM cells through remodeling of the ECM. Further functional analysis of dorsal-high Zic1-target genes identified in somites will provide useful insight into the molecular network driving dorsal somite extension.
Finally, in vertebrates and especially in fish, body shape and muscle morphology are closely linked, since a majority of the body mass consist of muscular tissue. Previous studies in fish populations have shown that speciation and adaptation to a specific aquatic habitat are associated with a changes in body depth, a measurement of the trunk dorsoventral axis (Tobler, 2008; Elmer et al., 2010; Weese et al., 2012; Fruciano, 2016). In this context, the external appearance of the adult Da mutant is intriguing, in that it exhibits a teardrop body shape (typical for fish swimming in the middle layer like tuna), instead of a dorsally flattened one (typical for surface swimming fish like medaka). Since our study shows that the activity of wnt11 could influence the body shape by regulating cell proliferation and the behavior of muscle progenitor cells, Wnt11 could be one of the crucial factors in evolution and diversity of body shape in fish. Furthermore, the expression of zic1 (Ohtsuka et al., 2004; Sun Rhodes and Merzdorf, 2006; Houtmeyers et al., 2013) and wnt11 (Marcelle et al., 1997; Olivera-Martinez et al., 2004; Garriock et al., 2005; Matsui et al., 2005) are strongly conserved among vertebrates, and we hypothesize that Wnt11-mediated morphogenesis of the somites represents an evolutionarily conserved mechanism that acts across vertebrates.
Materials and methods
Animals and transgenic lines
Request a detailed protocolFish were raised and maintained under standard conditions. All experimental procedures and animal care were performed according to the animal ethics committee of the University of Tokyo. Sex was randomly assigned to experimental groups. Medaka d-rR stain was used as wild type, the Da mutant used in this study was previously described (Ohtsuka et al., 2004). The pre-existing transgenic line Tg(zic1:GFP,zic4:DsRed) (Kawanishi et al., 2013) was used, and the transgenic line Tg(zic1:GFP,zic4:DsRed);Da was created by crossing Da mutants with Tg(zic1:GFP,zic4:DsRed). The transgenic line Tg(zic1:zic1-Myc,zic4:DsRed);Da was generated by modifying the BAC used to generate Tg(zic1:GFP,zic4:DsRed) (Kawanishi et al., 2013), by replacing the ORF of GFP with the ORF of zic1 containing a sequence for a Myc-tag fused to its C-terminus. To establish the transgenic line, the BAC(zic1:zic1-Myc,zic4:DsRed) was co-injected with I-SceI Meganuclease (NEB) into one-cell stage Da embryos, as previously described (Thermes, 2002).
Visualization of actin skeleton of protrusions
View detailed protocolThe AC-TagGFP2 sequence from the Actin-Chromobody plasmid (TagGFP2) (Chromotek) was cloned into the pMTB vector for mRNA generation. To investigate the actin skeleton of protrusions, cells of Wt and Da embryos were mosaically labeled using Actin-Chromobody-GFP (AC-GFP) mRNA. Embryos were injected at one-cell stage with 152 ng/µl membrane-mCherry mRNA. At four-cell stage, one cell was injected with 184 ng/µl AC-GFP mRNA. Embryos were raised until the desired stage and in vivo imaging was performed.
Lineage tracing of mesenchymal DM cells
Request a detailed protocolMosaic labeling of cells was achieved by co-injecting 20 ng/µl pMTB-memb-mTagBFP2 plasmid with Tol2 mRNA into one-cell stage embryos of the transgenic line Tg(zic1:GFP,zic4:DsRed). Embryos were raised until 5 dpf and labeled cells were continuously observed until 8 dpf. Observation of cells was performed according to “In vivo imaging and in vivo time-lapse imaging”.
Photoconversion of Dorsal DM Cells
Wild type embryos were injected with 200 ng/µl kikGR mRNA at one-cell stage. Single dorsal DM cells exhibiting protrusions at the dorsal tip of 10th and 20th somites at 3.5 dpf (st. 29) were selectively photoconverted with a 405 nm laser under an LSM 710 confocal microscope and imaged again at 9 dpf (hatching stage). During the experiment, 8 embryos were labeled at both 10th and 20th somites, and 1 embryo was labeled only at the 10th somite.
Morpholino injection
Request a detailed protocolwnt11 MO was injected at one-cell stage. We used 12.5 µM wnt11 MO after finding that injection of 25 µM wnt11 MO gives rise to gastrulation delay, leading to smaller head and trunk phenotypes at the onset of dorsal somite extension. Injection of 12.5 µM wnt11 MO mostly did not hamper gastrulation. We excluded any MO-injected embryos showing gastrulation defects in the analyses at later timepoints.
Photo-Morpholino injection
Request a detailed protocolwnt11 Sense-Photo-Morpholino (PhotoMO) and wnt11 antisense Morpholino were annealed in a ratio 2:1. Microinjection of 25 µM of the annealed oligonucleotides was performed at one-cell stage, and embryos were raised until 4 ss in the dark. Photo-cleavage was performed using the 10 x objective and the DAPI filter of a Keyence BZ-9000 Biorevo microscope (Keyence). Embryos were mounted, dorsal side facing up, in 1% Methylcellulose in a glass bottom dish (Wako) and illuminated for 30 min. After Photo-cleavage, embryos were dechorionated and raised until the desired stage for subsequent analysis.
Generation of sgRNAs
Request a detailed protocolThe online tool CCTop (Stemmer et al., 2015) was used to design sgRNAs. sgRNAs were generated as previously described by Hwang et al., 2013. Oligonucleotides were annealed and ligated into DR274 vector (Addgene), which was previously linearized using BsaI-HF (New England Biolabs). Template for in vivo transcription of sgRNA was amplified by PCR. In vivo transcription was performed using the HiScribe T7 Quick High Yield RNA Synthesis Kit (New England Biolabs) and RNA was purified using the RNeasy Mini kit (Qiagen). Quantity and quality of sgRNA were analyzed using a NanoDrop and agarose gel electrophoresis. sgRNA cleaving activity was confirmed by an in vivo assay.
sgRNA description | Target site [PAM] |
---|---|
wnt11_1 | TTCTGTCGGGAGACTTTA[TGG] |
wnt11_2 | GCTGGATGTTTAACGGAG[TGG] |
wnt11 sgRNA injection
Request a detailed protocolMicroinjections were performed using 25 ng/µl wnt11_1 sgRNA, 25 ng/µl wnt11_2 sgRNA and 100 ng/µl Cas9-NLS (New England Biolabs) at one-cell stage. To analyze phenotypes, embryos were mounted in 1% methylcellulose and imaged using a Leica M165 FC stereo microscope.
Injection of human recombinant Wnt11 protein into Da mutant somite
Request a detailed protocolTo immobilize embryos, embryos from the Tg(zic1:GFP,zic4:DsRed);Da transgenic line were injected with 25 ng/µl α-bungarotoxin mRNA at one-cell stage. At 18 ss embryos were mounted in 1% low melting agarose in 1 x Yamamoto’s Ringer Solution and oriented with the dorsal side facing upwards. Embryos were injected on top of the 10th somite with a mix containing Dextran Rhodamine (Thermo Fisher) and 1.7 ng hrWnt11 protein (R&D Systems) or BSA (Sigma-Aldrich) and raised to 24 ss, followed by in vivo imaging and analysis.
Laser ablation of mesenchymal DM cells
Request a detailed protocolTg(zic1:GFP,zic4:DsRed) embryos were immobilized by injecting 25 ng/µl α-bungarotoxin mRNA at one-cell stage and raised until 5.5 dpf. Z-stacks of the 9th (control somite) and 10th somite were imaged before ablation using the 40 x objective of a Zeiss LSM 710 confocal microscope. Laser ablation protocol was adapted from Morsch et al., 2017 and Volpe et al., 2020. Prior to laser ablation the 40 x objective and additional the 2.2 x zoom were used to focus only on the ablation site, the mesenchymal DM cells between the 10th somites. The ablation site was chosen to prevent damage of dorsal somites and the neural tube. Laser ablation was perforemd using a 405 nm laser at 80% power for 10 min on ‘continous’ mode. The 9th and 10th somite were subject to in vivo imaging directly after laser ablation and every 8 h until 40 hpa.
KN-93 treatment
Request a detailed protocolDechorionated 4 ss embryos were treated with 30 μM KN-93 (Wako) or DMSO (Sigma-Aldrich) in 1 x Yamamoto’s Ringer Solution at 28 °C, in the dark until 22 ss.
ML141 treatment
Request a detailed protocolEmbryos at 4 ss were dechorionated and treated with 500 μM ML141 (Sigma-Aldrich) or DMSO in 1 x Yamamoto’s Ringer Solution at 28 °C, in the dark until 24 ss.
In vivo imaging and in vivo time-lapse Imaging
Request a detailed protocolTo immobilize embryos for in vivo imaging, embryos were injected at one-cell stage with 25 ng/µl α-bungarotoxin mRNA (Swinburne et al., 2015; Lischik et al., 2019). Embryos were mounted in 1% Low melting agarose (Sigma-Aldrich) in 1 x Yamamoto’s Ringer Solution in a glass-bottomed petri dish (IWAKI) and oriented with the dorsal side facing down. Imaging was performed using a Zeiss LSM 710 confocal microscope system (Zeiss) equipped with an inverted stand and a Zeiss AXIO Observer Z1 and a T-PMT detector. The embryos were positioned with the 5th or 10th somite in the center and images were acquired using a 40 x water objective. For the in vivo time-lapse imaging, the 10th or 15th somite was positioned in the center, z-stacks were imaged in a 600 s interval for 10–15 hr. Image analysis was performed in Fiji using the “Image Stabilizer” Plugin, the FFT Bandpass filter and the “Draw_arrows” Plugin to draw customized arrows (Li, 2008; Daetwyler et al., 2020).
Whole mount in situ hybridization
Request a detailed protocolWhole mount in situ hybridization was performed as previously described with the following modifications (Takashima et al., 2007). Embryos were fixed in 4% PFA/1.5 x PTW at 4 °C, overnight. Hybridization was performed at 65 °C, overnight. Samples were treated with alkaline-phosphatide anti-DIG-AP Fab fragments (1:2000, Roche). Signals were developed using 4-nitro blue tetrazolium chloride (NBT, Roche) and 5-bromo-4-chloro-3-indolyl phosphate (BCIP, Roche).
Whole-mount immunohistochemistry
Request a detailed protocolEmbryos were fixed in 4% PFA/PBS for 2 hr at room temperature or at 4 °C overnight. Samples were permeabilized with 0.5% TritonX-100 (Wako) in 1 x PBS for 1–2 hr and blocked in blocking solution (2% BSA (Sigma-Aldrich), 1% DMSO, 0.2% TritonX-100 in 1 x PBS) for 2–4 hr at room temperature. Samples were incubated with respective primary antibody diluted in blocking solution at 4 °C, overnight. After an additional 4 hr blocking step, samples were incubated with respective secondary antibodies diluted in blocking solution, at 4 °C, overnight. Samples were stored in 1 x PBS at 4 °C until imaging.
Vibratome sectioning
Request a detailed protocolSamples were mounted in 4% agarose in 1 x PBS. 40 μm or 200 µm sections were obtained by a Vibratome (Leica, Vibratome). The sections were mounted on a glass slide (Matsunami) in 60% glycerol (Merck, Wako) and stored at 4 °C until imaging. Images were acquired using the 40 x water objective of a Zeiss LSM 710 confocal microscope.
Histological sections
Request a detailed protocolDechorionated embryos were fixed in Bouin’s solution overnight, followed by a gradual dehydration using ethanol. Samples were embedded in Technovit 7,100 (Heraeus Kulzer) and sectioned into 5–6 µm thick sections. Sections were stained with hematoxylin (Wako) and imaged using the 1.6 x objective of a Leica M165 FC fluorescent stereo microscope.
Image processing and statistical analysis
Request a detailed protocolImage processing was performed with the image processing software Fiji. The 3D-recreation of in vivo imaging date was created using FluoRender (Wan et al., 2009). Measurements of morphological features (distance between myotome tips, area of cross-section of dorsal somites, diameter of myofibers, distance between dorsal somite tip and the tip of neural tube, somite height) was performed by averaging the analysis of the feature from three consecutive Y planes. RStudio was used for the statistical analysis and representation of the data. In bar plots mean and error limits, defined by the standard deviation, are indicated. In box plots median first and third quantiles are indicated. Statistical significance was determined by un-paired t-tests, a p-value p < 0.05 was considered as significant. In the figure legends sample size (n) and number of individuals used in the experiment are stated. Sample sizes were not predetermined using statistical methods, but the sample sizes used are similar to those generally used in the field. To compare experimental groups, the allocation was performed randomly, without blinding.
RT-PCR of cDNA generated from embryonic tails
Request a detailed protocolTo investigate the gene expression in tails of embryos, tails were dissected anterior from the first somite. Ten tails were pooled together and RNA was isolated using Isogen (Nippon Gene). RNA was purified using the RNeasy Mini kit (Qiagen) and reverse transcribed to cDNA using the Super Script III Kit (Invitrogen). RT-PCR was performed using the Thunderbird Sybr qPCR Mix (Toyobo) following manufacturer’s instructions and run in the Agilent Mx3000P qPCR System (Agilent). Normalization of relative quantities was performed against gapdh expression, followed by analysis with excel and RStudio.
Isolation of dorsal and ventral somite cells for ATAC-Seq and RNA-Seq
Request a detailed protocolYolk and head region were removed from 22 ss Tg(zic1:GFP,zic4:DsRed) embryos, and remaining trunk-tail pieces were incubated with 10 mg/ml pancreatin (Wako) at room temperature for 5–10 min to loosen the adhesion between tissues while maintaining the integrity within each tissue. Epidermis and intestinal tissues were removed, and then somites were isolated from the neural tube. Collected somites were subsequently dissociated into individual cells in 0.5% (w/v) Trypsin (Nacalai Tesque) at 37 °C for 10 min, and the dissociation was stopped by adding the same volume of 15% (v/v) FBS / Leiboviz’s L-15 (Life Technologies). Dissociated cells were washed with PBS, and sorted into GFP-positive (dorsal) and negative (ventral) cells using FACSAria III (BD Biosciences). Dead cells were detected by Propidium iodide (Life Technologies) and removed.
RNA-Seq
Request a detailed protocolTotal RNA was extracted from sorted somite cells using RNeasy Mini kit (Qiagen). mRNA was enriched by poly-A capture and mRNA-seq libraries were generated using KAPA Stranded mRNA-seq Kit (KAPA Biosystems). Libraries were generated from two biological replicates, and sequenced using the Illumina HiSeq 1,500 platform.
RNA-Seq data processing
Request a detailed protocolThe sequenced reads were pre-processed to remove low-quality bases and adapter derived sequences using Trimmomatic v0.32 (Bolger et al., 2014), and aligned to the medaka reference genome (HdrR, ASM223467v1) by STAR (Dobin et al., 2013). Reads with mapping quality (MAPQ) larger than or equal to 20 were used for the further analyses.
ATAC-Seq
Request a detailed protocolATAC-seq was performed as previously described (Buenrostro et al., 2013) with some modifications. Approximately 4000 sorted somite cells were used for each experiment. After washing with PBS, cells were resuspended in 500 μl cold lysis buffer (10 mM Tris-HCl pH 7.4, 10 mM NaCl, 3 mM MgCl2, 0.1% Igepal CA-630), centrifuged for 10 min at 500 g, supernatant was removed. Tagmentation reaction was performed as described previously (Buenrostro et al., 2013) with Nextera Sample Preparation Kit (Illumina). After tagmented DNA was purified using MinElute kit (Qiagen), two sequential PCRs were performed to enrich small DNA fragments. First, a 9-cycle PCR was performed using indexed primers from Nextera Index Kit (Illumina) and KAPA HiFi HotStart ReadyMix (KAPA Biosystems), amplified DNA was size selected for a size less than 500 bp using AMPure XP beads (Beckman Coulter). A second 7-cycle PCR was performed using the same primers as for the first PCR. PCR product was purified by AMPure XP beads. Libraries were generated from two biological replicates, and sequenced using the Illumina HiSeq 1,500 platform.
ChIPmentation
Request a detailed protocolYolk and head region were removed from 22 ss Tg(zic1:zic1-Myc,zic4:DsRed);Da embryos, followed by an incubation with 10 mg/ml pancreatin (Wako) at room temperature for 5–10 min. Epidermis and intestinal tissues were removed and somites were isolated from the neural tube. ChIP was performed as previously described with the following modifications (Nakamura et al., 2014). Isolated somites were fixed with 1% formaldehyde for 8 min at room temperature then quenched by adding glycine (200 mM final) and incubating on ice for 5 min. After washing with PBS, cell pellets were stored at –80 ℃. Approximately 1.8 × 106 cells were thawed on ice, suspended in lysis buffer (50 mM Tris-HCl (pH 8.0), 10 mM EDTA, 1% SDS, 20 mM Na-butyrate, complete protease inhibitors, 1 mM PMSF) and sonicated 10 times using a Sonifier (Branson) at power 5. Chromatin lysates were collected by centrifugation and diluted 10-fold with RIPA ChIP buffer (10 mM Tris-HCl (pH 8.0), 140 mM NaCl,1 mM EDTA, 0.5 mM EGTA, 1% Triton X-100, 0.1% SDS, 0.2% sodium deoxycholate, 20 mM Na-butyrate, complete protease inhibitors, 1 mM PMSF) followed by an incubation with antibody/protein A Dynabeads (Invitrogen) complex at 4 °C, overnight, while rotating. Immunoprecipitated samples were washed three times with RIPA buffer (10 mM Tris-HCl (pH 8.0), 140 mM NaCl,1 mM EDTA, 0.5 mM EGTA, 1% Triton X-100, 0.1% SDS, 0.2% sodium deoxycholate) and once with TE buffer. After the washing steps 150 μl of Tris-HCl was added.
Library preparation for ChIPmentation was performed as previously described (Schmidl et al., 2015) with the following modifications. Twenty-four μl of Tagmentation reaction mix (10 mM Tris-HCl pH8.0, 5 mM MgCl2, 10%(v/v) N,N-dimethyl formamide) and 1 μl of Tagment DNA Enzyme from Nextera Sample Preparation Kit (Illumina) were added to the DNA-beads complex and incubated for 70 s at 37 °C. A total of 150 μl ice-cold RIPA buffer was added and incubated for 5 min on ice. The DNA-beads complex was washed with RIPA buffer and TE buffer, suspended in 50 μl lysis buffer and 3 μl of 5 M NaCl, and incubated at 65 °C, overnight. The sample was incubated for 2 hr with 2 μl of 20 mg/ml ProteinasK (Roche), and subjected to AMPure XP beads (Beckman Coulter) purification. The library was amplified by 18-cycle PCR using indexed primers from Nextera Index Kit (Illumina) and KAPA HiFi HotStart ReadyMix (KAPA Biosystems). For the input chromatin, tagmentation reaction was performed after DNA purification. Libraries were generated from two biological replicates, and sequenced using the Illumina HiSeq 1,500 platform.
ChIPmentation and ATAC-Seq data processing
Request a detailed protocolThe sequenced reads were pre-processed to remove low-quality bases and adapter derived sequences using Trimmomatic v0.32 (Bolger et al., 2014) and aligned to the medaka reference genome (HdrR, ASM223467v1) by BWA (Li and Durbin, 2009). Reads with mapping quality (MAPQ) larger than or equal to 20 were used for the further analyses. MACS2 (version 2.1.1.20160309) (Zhang et al., 2008) was used to call peaks and generate signals per million reads tracks using following options; ChIPmentation: -g 600000000 -B --SPMR --keep-dup 2, ATAC-seq: --nomodel --extsize 200 --shift –100 g 600000000 -q 0.01 -B --SPMR. The input chromatin was used as control of ChIPmentation.
For ChIPmentation, peak regions called by two biological replicates were used as reliable peaks.
Motif analyses of Zic1 ChIPmentation peaks
Request a detailed protocolMotifs enriched at reliable ChIPmentation peaks were analyzed by findMotifsGenome command of HOMER (Heinz et al., 2010) using default parameters.
Identification of Zic1 target genes
Request a detailed protocolDifferentially expressed genes were identified using DESeq2 (padj <0.01). Each reliable ChIP peak was associated to the nearest TSS, and the gene was defined as Zic1-target gene if the distance between the peak and the TSS was closer than 50 kb.
Gene ontology and pathway analyses
Request a detailed protocolThe gene ontology enrichment analyses and pathway enrichment analyses were performed using the Gene Ontology Resource (Ashburner et al., 2000; Gene Ontology Consortium, 2021; Gene Ontology Consortium, 2021).
Data availability
Sequencing data generated in this study have been submitted to the DDBJ BioProject database under accession number PRJDB11712.
-
DDBJ BioProjectID PRJDB11712. The function of transcription factor Zic1 in dorsalization of the medaka somite.
References
-
Semaphorins and their Signaling MechanismsMethods in Molecular Biology 22:11–14.https://doi.org/10.1016/s0896-6273(00)80672-8
-
Clarification of mural cell coverage of vascular endothelial cells by live imaging of zebrafishDevelopment (Cambridge, England) 143:1328–1339.https://doi.org/10.1242/dev.132654
-
ZIC1 Function in Normal Cerebellar Development and Human Developmental PathologyAdvances in Experimental Medicine and Biology 1046:249–268.https://doi.org/10.1007/978-981-10-7311-3_13
-
Gene Ontology: tool for the unification of biologyNature Genetics 25:25–29.https://doi.org/10.1038/75556
-
Lineage analysis of the avian dermomyotome sheet reveals the existence of single cells with both dermal and muscle progenitor fatesDevelopment (Cambridge, England) 132:689–701.https://doi.org/10.1242/dev.01617
-
Trimmomatic: A flexible trimmer for Illumina sequence dataBioinformatics (Oxford, England) 30:2114–2120.https://doi.org/10.1093/bioinformatics/btu170
-
Pax group III genes and the evolution of insect pair-rule patterningDevelopment (Cambridge, England) 128:3445–3458.https://doi.org/10.1242/dev.128.18.3445
-
Essential role of non-canonical Wnt signalling in neural crest migrationDevelopment (Cambridge, England) 132:2587–2597.https://doi.org/10.1242/dev.01857
-
Developmental fate of the mammalian myotomeDevelopmental Dynamics 239:2898–2910.https://doi.org/10.1002/dvdy.22425
-
Extracellular matrix remodeling accompanies axial muscle development and morphogenesis in the mouseDevelopmental Dynamics 241:350–364.https://doi.org/10.1002/dvdy.23703
-
STAR: ultrafast universal RNA-seq alignerBioinformatics (Oxford, England) 29:15–21.https://doi.org/10.1093/bioinformatics/bts635
-
Local variation and parallel evolution: morphological and genetic diversity across a species complex of neotropical crater lake cichlid fishesPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 365:1763–1782.https://doi.org/10.1098/rstb.2009.0271
-
Myogenin promotes myocyte fusion to balance fibre number and sizeNature Communications 9:4232.https://doi.org/10.1038/s41467-018-06583-6
-
The Gene Ontology resource: enriching a GOld mineNucleic Acids Research 49:D325–D334.https://doi.org/10.1093/nar/gkaa1113
-
A two-Step mechanism for myotome formation in chickDevelopmental Cell 6:875–882.https://doi.org/10.1016/j.devcel.2004.05.006
-
Met and Hgf signaling controls hypaxial muscle and lateral line development in the zebrafishDevelopment (Cambridge, England) 131:4857–4869.https://doi.org/10.1242/dev.01374
-
Vertebrate myotome developmentBirth Defects Research. Part C, Embryo Today 75:172–179.https://doi.org/10.1002/bdrc.20046
-
Characterization of a Cdc42 protein inhibitor and its use as a molecular probeThe Journal of Biological Chemistry 288:8531–8543.https://doi.org/10.1074/jbc.M112.435941
-
The ZIC gene family encodes multi-functional proteins essential for patterning and morphogenesisCellular and Molecular Life Sciences 70:3791–3811.https://doi.org/10.1007/s00018-013-1285-5
-
Efficient genome editing in zebrafish using a CRISPR-Cas systemNature Biotechnology 31:227–229.https://doi.org/10.1038/nbt.2501
-
Myotome formation: A multistage processCell and Tissue Research 296:161–173.https://doi.org/10.1007/s004410051277
-
Modular development of the teleost trunk along the dorsoventral axis and zic1/zic4 as selector genes in the dorsal moduleDevelopment (Cambridge, England) 140:1486–1496.https://doi.org/10.1242/dev.088567
-
An exclusively mesodermal origin of fin mesenchyme demonstrates that zebrafish trunk neural crest does not generate ectomesenchymeDevelopment (Cambridge, England) 140:2923–2932.https://doi.org/10.1242/dev.093534
-
The Xfeb gene is directly upregulated by Zic1 during early neural developmentDevelopmental Dynamics 235:2817–2827.https://doi.org/10.1002/dvdy.20896
-
Fast and accurate short read alignment with Burrows-Wheeler transformBioinformatics (Oxford, England) 25:1754–1760.https://doi.org/10.1093/bioinformatics/btp324
-
Coordinate actions of BMPs, Wnts, Shh and noggin mediate patterning of the dorsal somiteDevelopment (Cambridge, England) 124:3955–3963.https://doi.org/10.1242/dev.124.20.3955
-
Wnt11r is required for cranial neural crest migrationDevelopmental Dynamics 237:3404–3409.https://doi.org/10.1002/dvdy.21758
-
Triggering Cell Stress and Death Using Conventional UV Laser Confocal MicroscopyJournal of Visualized Experiments 120:54983.https://doi.org/10.3791/54983
-
Large hypomethylated domains serve as strong repressive machinery for key developmental genes in vertebratesDevelopment (Cambridge, England) 141:2568–2580.https://doi.org/10.1242/dev.108548
-
CTCF looping is established during gastrulation in medaka embryosGenome Research 31:968–980.https://doi.org/10.1101/gr.269951.120
-
Molecular mechanisms controlling dorsal dermis generation from the somitic dermomyotomeThe International Journal of Developmental Biology 48:93–101.https://doi.org/10.1387/ijdb.15272374
-
Wnt11 signaling promotes proliferation, transformation, and migration of IEC6 intestinal epithelial cellsThe Journal of Biological Chemistry 279:26707–26715.https://doi.org/10.1074/jbc.M402877200
-
A role for Zic1 and Zic2 in Myf5 regulation and somite myogenesisDevelopmental Biology 351:120–127.https://doi.org/10.1016/j.ydbio.2010.12.037
-
A defect in A novel ADAMTS family member is the cause of the belted white-spotting mutationDevelopment (Cambridge, England) 130:4665–4672.https://doi.org/10.1242/dev.00668
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
Connecting muscle development, birth defects, and evolution: An essential role for muscle connective tissueCurrent Topics in Developmental Biology 132:137–176.https://doi.org/10.1016/bs.ctdb.2018.12.004
-
Dynamic somite cell rearrangements lead to distinct waves of myotome growthDevelopment (Cambridge, England) 134:1253–1257.https://doi.org/10.1242/dev.000067
-
The newly synthesized selective Ca2+/calmodulin dependent protein kinase II inhibitor KN-93 reduces dopamine contents in PC12h cellsBiochemical and Biophysical Research Communications 181:968–975.https://doi.org/10.1016/0006-291x(91)92031-e
-
The zic1 gene is expressed in chick somites but not in migratory neural crestGene Expression Patterns 6:539–545.https://doi.org/10.1016/j.modgep.2005.10.006
-
Phenotypic analysis of a novel chordin mutant in medakaDevelopmental Dynamics 236:2298–2310.https://doi.org/10.1002/dvdy.21245
-
Muscle precursor cell movements in zebrafish are dynamic and require Six family genesDevelopment (Cambridge, England) 146:10.https://doi.org/10.1242/dev.171421
-
Turning gene function ON and OFF using sense and antisense photo-morpholinos in zebrafishDevelopment (Cambridge, England) 139:1691–1699.https://doi.org/10.1242/dev.072702
-
I-SceI meganuclease mediates highly efficient transgenesis in fishMechanisms of Development 118:91–98.https://doi.org/10.1016/S0925-4773(02)00218-6
-
Toxic hydrogen sulfide and dark caves: Phenotypic and genetic divergence across two abiotic environmental gradients in Poecilia mexicanaEvolution; International Journal of Organic Evolution 62:2643–2659.https://doi.org/10.1111/j.1558-5646.2008.00466.x
-
G1 cell cycle arrest and apoptosis are induced in NIH 3T3 cells by KN-93, an inhibitor of CaMK-II (the multifunctional Ca2+/CaM kinaseCell Growth & Differentiation 6:1063–1070.
-
Slb/Wnt11 controls hypoblast cell migration and morphogenesis at the onset of zebrafish gastrulationDevelopment (Cambridge, England) 130:5375–5384.https://doi.org/10.1242/dev.00758
-
Confocal Microscope-Based Laser Ablation and Regeneration Assay in Zebrafish Interneuromast CellsJournal of Visualized Experiments 159:60966.https://doi.org/10.3791/60966
-
An interactive visualization tool for multi-channel confocal microscopy data in neurobiology researchIEEE Transactions on Visualization and Computer Graphics 15:1489–1496.https://doi.org/10.1109/TVCG.2009.118
-
Stabilization of dendritic arbor structure in vivo by CaMKIIScience (New York, N.Y.) 279:222–226.https://doi.org/10.1126/science.279.5348.222
-
Model-based analysis of ChIP-Seq (MACSGenome Biology 9:R137.https://doi.org/10.1186/gb-2008-9-9-r137
Article and author information
Author details
Funding
Japan Society for the Promotion of Science (JP15H05859)
- Hiroyuki Takeda
Japan Society for the Promotion of Science (JP19K23741)
- Toru Kawanishi
Japan Society for the Promotion of Science (JP18K14620)
- Ryohei Nakamura
Japan Science and Technology Agency (JPMJCR13W3)
- Hiroyuki Takeda
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank the members of the Takeda laboratory for constructive feedback and discussions on the project. We are greatful for Y Yamagichi and M Funato for fish husbandry. This work was supported by Japan Society for the Promotion of Science (JSPS) KAKENHI Grant Numbers JP15H05859 (HT), JP19K23741 (TK) and JP18K14620 (RN) and Japan Science and Technology Agency CREST Grant Number JPMJCR13W3 (HT).
Ethics
All experimental procedures and animal care were performed according to the animal ethics committee of the University of Tokyo (Approval No. 20-02).
Copyright
© 2022, Heilig 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.
Metrics
-
- 1,175
- views
-
- 136
- downloads
-
- 7
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
Organogenesis requires the proper production of diverse cell types and their positioning/migration. However, the coordination of these processes during development remains poorly understood. The gonad in C. elegans exhibits a mirror-symmetric structure guided by the migration of distal tip cells (DTCs), which result from asymmetric divisions of somatic gonadal precursors (SGPs; Z1 and Z4). We found that the polarity of Z1 and Z4, which possess mirror-symmetric orientation, is controlled by the redundant functions of the LIN-17/Frizzled receptor and three Wnt proteins (CWN-1, CWN-2, and EGL-20) with distinct functions. In lin-17 mutants, CWN-2 promotes normal polarity in both Z1 and Z4, while CWN-1 promotes reverse and normal polarity in Z1 and Z4, respectively. In contrast, EGL-20 inhibits the polarization of both Z1 and Z4. In lin-17 egl-20 cwn-2 triple mutants with a polarity reversal of Z1, DTCs from Z1 frequently miss-migrate to the posterior side. Our further analysis demonstrates that the mis-positioning of DTCs in the gonad due to the polarity reversal of Z1 leads to mis-migration. Similar mis-migration was also observed in cki-1(RNAi) animals producing ectopic DTCs. These results highlight the role of Wnt signaling in coordinating the production and migration of DTCs to establish a mirror-symmetric organ.
-
- Cell Biology
- Developmental Biology
Chondrocyte columns, which are a hallmark of growth plate architecture, play a central role in bone elongation. Columns are formed by clonal expansion following rotation of the division plane, resulting in a stack of cells oriented parallel to the growth direction. In this work, we analyzed hundreds of Confetti multicolor clones in growth plates of mouse embryos using a pipeline comprising 3D imaging and algorithms for morphometric analysis. Surprisingly, analysis of the elevation angles between neighboring pairs of cells revealed that most cells did not display the typical stacking pattern associated with column formation, implying incomplete rotation of the division plane. Morphological analysis revealed that although embryonic clones were elongated, they formed clusters oriented perpendicular to the growth direction. Analysis of growth plates of postnatal mice revealed both complex columns, composed of ordered and disordered cell stacks, and small, disorganized clusters located in the outer edges. Finally, correlation between the temporal dynamics of the ratios between clusters and columns and between bone elongation and expansion suggests that clusters may promote expansion, whereas columns support elongation. Overall, our findings support the idea that modulations of division plane rotation of proliferating chondrocytes determines the formation of either clusters or columns, a multifunctional design that regulates morphogenesis throughout pre- and postnatal bone growth. Broadly, this work provides a new understanding of the cellular mechanisms underlying growth plate activity and bone elongation during development.