Zebrafish embryonic explants undergo genetically encoded self-assembly
Abstract
Embryonic stem cell cultures are thought to self-organize into embryoid bodies, able to undergo symmetry-breaking, germ layer specification and even morphogenesis. Yet, it is unclear how to reconcile this remarkable self-organization capacity with classical experiments demonstrating key roles for extrinsic biases by maternal factors and/or extraembryonic tissues in embryogenesis. Here, we show that zebrafish embryonic tissue explants, prepared prior to germ layer induction and lacking extraembryonic tissues, can specify all germ layers and form a seemingly complete mesendoderm anlage. Importantly, explant organization requires polarized inheritance of maternal factors from dorsal-marginal regions of the blastoderm. Moreover, induction of endoderm and head-mesoderm, which require peak Nodal-signaling levels, is highly variable in explants, reminiscent of embryos with reduced Nodal signals from the extraembryonic tissues. Together, these data suggest that zebrafish explants do not undergo bona fide self-organization, but rather display features of genetically encoded self-assembly, where intrinsic genetic programs control the emergence of order.
Introduction
Embryonic development entails both progressive specification of cells into distinct cell fates and their respective allocation into different domains within the embryo. Key signaling molecules, also termed morphogens, play critical roles in eliciting both the genetic programs underlying cell fate specification and the morphogenetic cell properties required to shape the developing embryo (Gilmour et al., 2017; Heisenberg and Solnica-Krezel, 2008; Pinheiro and Heisenberg, 2020). While amphibians critically rely on maternal factors for precise spatiotemporal activation of morphogen signaling during early embryonic development, mammalian embryos mostly depend on reciprocal interactions between embryonic and extraembryonic tissues (Mummery and Chuva de Sousa Lopes SM, 2013; Heasman, 2006; Langdon and Mullins, 2011; Stern and Downs, 2012; Tam and Loebel, 2007). Interestingly, zebrafish embryos combine features characteristic of both amphibian and mammalian development. Here, maternal factors initially deposited within the oocyte (Langdon and Mullins, 2011; Marlow, 2020) and signals from extraembryonic tissues, namely the yolk cell and its yolk syncytial layer (YSL), play important roles in germ layer induction and patterning (Carvalho and Heisenberg, 2010; Chen and Kimelman, 2000; Erter et al., 1998; Fan et al., 2007; Gagnon et al., 2018; Hong et al., 2011; Mizuno et al., 1999; Mizuno et al., 1996; Ober and Schulte-Merker, 1999; Schier and Talbot, 2005; van Boxtel et al., 2015; Veil et al., 2018; Xu et al., 2012). Nodal/TGFβ ligands, for instance, key signals involved in mesoderm and endoderm (mesendoderm) specification and patterning, are initially expressed under maternal control and then secreted by the extraembryonic YSL to pattern the embryonic blastoderm (Chen and Schier, 2002; Dougan et al., 2003; Erter et al., 1998; Fan et al., 2007; Fauny et al., 2009; Feldman et al., 1998; Gagnon et al., 2018; Gore et al., 2005; Gritsman et al., 2000; Gritsman et al., 1999; Hong et al., 2011; Rebagliati et al., 1998; Rogers et al., 2017; Rogers and Müller, 2019; Thisse and Thisse, 2015; van Boxtel et al., 2015; Xu et al., 2012).
Pioneering studies have suggested that both human and mouse embryonic stem cell (ES) colonies exhibit a remarkable capacity to self-organize in vitro and recapitulate features associated with early embryogenesis (Baillie-Johnson et al., 2015; Beccari et al., 2018; Doetschman et al., 1985; Shahbazi et al., 2019; Simunovic et al., 2019; Simunovic and Brivanlou, 2017; ten Berge et al., 2008; Turner et al., 2017; Turner et al., 2014; van den Brink et al., 2014; van den Brink et al., 2020; Warmflash et al., 2014). For instance, homogeneous stimulation of human ES colonies cultured on 2-dimensional circular substrates with BMP4 ligands is sufficient to reproducibly generate radially symmetric colonies, where cells at the outmost rim differentiate into extraembryonic tissue, while the remaining cells differentiate progressively into more internal rings composed of the different embryonic germ layers - endoderm, mesoderm and ectoderm (Warmflash et al., 2014). These apparent self-organizing processes are thought to rely on boundary sensing, community effects and the interplay between activators and inhibitors (Chhabra et al., 2019; Etoc et al., 2016; Nemashkalo et al., 2017; Sasai, 2013; Xavier da Silveira Dos Santos and Liberali, 2019). While these pioneering studies suggest that embryonic tissues might have some self-organizing capacity, classical studies have also established the dependency of embryonic development on both maternal factors and extraembryonic tissues, thereby raising key questions as to their specific contributions to robust embryonic development.
Results
To address this question, we established an ex vivo system to culture zebrafish blastoderm tissue explants lacking the yolk cell and its associated YSL. For this, blastoderm explants were prepared from embryos at early cleavage stages (256 cell stage), before YSL formation or germ layer induction, and kept in serum-free medium until stage-matched control embryos had completed gastrulation (bud stage; Figure 1A). We found that these blastoderm explants initially rounded up in culture, but eventually became ‘pear’-shaped by forming an extension (Figure 1B,C, Video 1). Interestingly, this extension began to form when stage-matched control embryos initiated body axis extension at shield stage (Sepich et al., 2005; Solnica-Krezel and Sepich, 2012), pointing at the possibility that these processes might be related. Contrary to intact embryos, however, blastoderm explants formed a large luminal cavity, reminiscent of a blastocoel (Krens et al., 2017), which was typically positioned at the opposite end of the explant extension (Figure 1B, Figure 1—figure supplement 1A,B, Video 2). Importantly, cells within explants appeared to normally complete their rounds of cleavages and divisions following explant preparation, as evidenced by the similar average diameter of explant and stage-matched embryo cells (Figure 1—figure supplement 1C). Moreover, we found that the onset of zygotic genome activation (ZGA), as monitored by the expression of the two zygotic genes tbx16 and chrd (Gagnon et al., 2018; Lee et al., 2013), was similar in stage-matched embryos and blastoderm explants (Figure 1—figure supplement 1D). Finally, comparable morphological changes could be observed when blastoderm explants were prepared from earlier (64 and 128 cell-stage) or later (high stage) embryos (Figure 1—figure supplement 2A–D), suggesting that the exact time point of explant preparation during pre-gastrula stages has no decisive influence on explant morphogenesis.
To analyze how these distinct changes in explant morphology relate to the specification of different progenitor cell types therein, we analyzed the expression of various progenitor cell marker genes by whole mount in situ hybridization. Strikingly, we found distinct expression domains for all tested marker genes outlining ectoderm, neuroectoderm, endoderm and different types of mesoderm progenitor cell populations. While ectoderm (gata2) and neuroectoderm (six3) marker genes were expressed predominantly at the opposite end of the explant extended portion, the region where the cavity was also forming (Figure 1D–G, Figure 1—figure supplement 1A,B,H, Figure 1—figure supplement 3A), all mesendoderm and endoderm marker genes were exclusively found within or adjacent to the extension (Figure 1D–G, Figure 1—figure supplement 1G–J, Figure 1—figure supplement 3A). Specifically, we found that markers of posterior axial mesoderm (ntl and flh) were expressed in the central region of the extension, flanked by paraxial and lateral mesoderm (papc, myoD, tbx6) (Figure 1D–G, Figure 1—figure supplement 1G–I, Figure 1—figure supplement 3A). In contrast, anterior axial mesendoderm markers (hgg and gsc) were found to be expressed at variable positions along the extended region or directly adjacent to it (Figure 1D–G, Figure 1—figure supplement 1G–I, Figure 1—figure supplement 3A). Finally, sporadic endoderm progenitors (sox32) could be found throughout the extension (Figure 1G, Figure 1—figure supplement 1G,H,J, Figure 1—figure supplement 3A). On the explant surface, much like in the intact embryo, we found that cells expressed elevated levels of Keratin 4, a marker of enveloping layer (EVL) cell differentiation (Figure 1—figure supplement 1E,F). Collectively, these findings suggest that the blastoderm explants contain all main germ layer progenitors also found in intact embryos and display a seemingly complete patterning of the mesendoderm anlage. They also indicate that, although explants undergo extensive morphogenesis, mesendoderm progenitors do not clearly internalize and, instead, organize into an extension (Video 3). This is highly reminiscent of the ‘exogastrula’ phenotype observed in embryos with defective mesendoderm internalization and, more recently, in mouse gastruloids (Beccari et al., 2018; Branford and Yost, 2002; Turner et al., 2017; van den Brink et al., 2014).
To ensure that all extraembryonic tissues were efficiently removed during explant preparation and that explants do not form a secondary YSL, we analyzed the expression of marker genes, normally present within the YSL in intact embryos, in blastoderm explants. We found that all of these genes were either not expressed in explants or showed only sporadic expression in small groups of blastoderm cells (Figure 1—figure supplement 4A,B). Moreover, we failed to detect any clearly recognizable syncytium in explants, a hallmark of the YSL in intact embryos (Figure 1—figure supplement 4C). Together, these data strongly support the notion that germ layer specification and mesendoderm patterning in blastoderm explants can occur in the absence of the extraembryonic yolk and YSL.
To understand how mesendoderm tissues are formed within explants, we first analyzed whether Nodal ligands, critical for mesendoderm specification and patterning in vivo (Erter et al., 1998; Feldman et al., 1998; Gritsman et al., 2000; Gritsman et al., 1999; Rebagliati et al., 1998), are also expressed in explants. We observed expression of both Nodal ligands sqt and cyc in blastoderm explants (Figure 2—figure supplement 1A,B), consistent with previous observations in intact embryos (Erter et al., 1998; Fan et al., 2007; Feldman et al., 1998; Rebagliati et al., 1998; van Boxtel et al., 2015). In addition to this, we detected nuclear localization of phosphorylated SMAD2/3 (pSMAD2/3), a well-established readout of Nodal/TGFβ signaling activation (Economou and Hill, 2020; Schier, 2009; van Boxtel et al., 2015; van Boxtel et al., 2018), in cells located close to the wounding site of explants (Figure 2A, Figure 2—figure supplement 1C,D), indicative of local Nodal signaling activation. A comparison of the spatiotemporal dynamics of Nodal signaling establishment between explants and intact embryos further revealed that, while blastoderm explants showed initially fewer and less bright pSMAD2/3 positive nuclei than intact embryos, both parameters continued to increase in explants until stage-matched embryos had reached shield stage (Figure 2A–C, Figure 2—figure supplement 1C–J). This differs from the temporal dynamics of Nodal signaling in intact embryos, where nuclear accumulation of pSMAD2/3 peaks already by 50% epiboly (Figure 2A–C, Figure 2—figure supplement 1C–J). Notably, the wounding site of blastoderm explants, where Nodal signaling was activated corresponds to the former blastoderm margin where explants were removed from the yolk cell (Figure 2D, Figure 2—figure supplement 1E–J). Given that Nodal signaling in intact embryos is initiated in the blastoderm margin (Dubrulle et al., 2015; Harvey and Smith, 2009; Schier, 2009; van Boxtel et al., 2015), this points at the intriguing possibility that Nodal signaling is active in similar cells in embryos and explants. This notion was further substantiated by our observation that cells close to the wounding site in explants display only very limited dispersal during explant maturation (Figure 2—figure supplement 2A–D). Moreover, similar to intact embryos, a gradient of BMP signaling was also detectable in blastoderm explants (Figure 2—figure supplement 2E,F), suggesting that explants might also use a BMP signaling gradient for patterning along their dorsoventral axis.
To test whether and when Nodal signaling is required for mesendoderm induction and morphogenesis in blastoderm explants, we prepared explants from maternal-zygotic (MZ) mutants of the Nodal co-receptor Oep (MZoep), previously shown to be essential for Nodal signaling and proper mesendoderm induction during embryogenesis (Gritsman et al., 1999). We found that blastoderm explants obtained from MZoep mutants failed to extend, indicative of strongly reduced mesendoderm formation (Figure 2E–G). This suggests that, similar to the situation in intact embryos, proper Nodal signaling is critical for mesendoderm induction and morphogenesis in blastoderm explants. To determine when Nodal signaling is required for mesendoderm induction and morphogenesis, we cultured explants in the presence of the selective and reversible Nodal inhibitor SB-505124 (50 µM) (DaCosta Byfield et al., 2004; Fan et al., 2007; Hagos and Dougan, 2007) for different time-periods during explant maturation. We found that blocking Nodal signaling from 256 cell to shield stage was sufficient to completely abolish mesendoderm induction and explant extension (Figure 2E–G, Figure 2—figure supplement 3A–D), while blocking Nodal signaling after shield stage had no strong impact on these processes (Figure 2E–H, Figure 2—figure supplement 3A–D). This is highly reminiscent of the early requirement of Nodal signaling for mesendoderm induction and morphogenesis in intact embryos (Gritsman et al., 1999; Hagos and Dougan, 2007; Pinheiro and Heisenberg, 2020), suggesting that Nodal signaling functions similarly in embryos and blastoderm explants.
Nodal-dependent mesendoderm induction has previously been shown to promote convergence and extension movements and body axis elongation by inducing the non-canonical Wnt/Planar Cell Polarity (Wnt/PCP) signaling module (Djiane et al., 2000; Gray et al., 2011; Tada and Heisenberg, 2012; Wallingford, 2012; Williams and Solnica-Krezel, 2019). To analyze whether explant extension also relies on non-canonical Wnt/PCP signaling, we prepared explants from MZ mutant embryos for different components of the Wnt/PCP signaling module, which display strongly diminished body axis elongation in vivo. Blastoderm explants prepared from MZ mutants for the Wnt/PCP components wnt11, wnt5b and fz7a/b also exhibited reduced or no extension (Figure 3A–D). This suggests that explant extension, similar to body axis elongation in intact embryos, critically depends on non-canonical Wnt/PCP signaling. In line with this, we also found that during explant maturation, the extension length increased, and the width of the non-extended portion decreased (Figure 3—figure supplement 1A–C), and that explant extension was largely independent of cell proliferation (Figure 3—figure supplement 1D–F). This points at the intriguing possibility that explants, like intact embryos, might undergo some type of convergence and extension movements, independent of cell proliferation.
Our data, so far, suggest that mesendoderm specification and morphogenesis in explants relies on the activation of Nodal signaling around the wounding site. In intact embryos, Nodal ligand expression, and thus activation of Nodal signaling, at the blastoderm margin depends both on initial maternal activation of the transcriptional co-activator β-catenin (β-cat) and Nodal signals emanating from the YSL (van Boxtel et al., 2015; Feldman et al., 1998; Gagnon et al., 2018; Hong et al., 2011; Langdon and Mullins, 2011; Marlow, 2020; Xu et al., 2012). Since explants contain, at most, remnants of the YSL (Figure 1—figure supplement 4A–C), we tested whether Nodal signaling activation in this context might be triggered by β-cat. Indeed, analysis of nuclear localization of β-cat in explants shortly after preparation revealed elevated levels of nuclear β-cat in blastoderm cells close to the wounding site (Figure 4A–C). This suggests that in explants, similar to intact embryos, β-cat nuclear accumulation might be responsible for Nodal signaling activation, and consequently, mesendoderm induction and morphogenesis.
Nuclear accumulation of β-cat is thought to be triggered by maternal factors, which are transported to the future dorsal side of the zebrafish embryo upon fertilization (Langdon and Mullins, 2011; Marlow, 2020). To determine whether this β-cat nuclear accumulation at the dorsal blastoderm margin is indeed responsible for Nodal signaling activation and mesendoderm induction in explants, we prepared explants from embryos depleted of maternal dorsal determinants (DD-removed) and analyzed the formation of the Nodal signaling domain therein. To this end, we mechanically removed the vegetal-most part of the zygote cortex, where these factors are thought to localize prior to transport to the future dorsal side of the embryo soon after fertilization (Langdon and Mullins, 2011; Marlow, 2020; Mizuno et al., 1999; Ober and Schulte-Merker, 1999). As expected, in blastoderm explants prepared from DD-removed embryos, a large portion (>50%) of which displayed a fully ventralized phenotype one day after fertilization (1 dpf), nuclear accumulation of β-cat was strongly diminished (Figure 4—figure supplement 1A–D). Interestingly, we found that explants originating from DD-removed embryos showed a clear reduction in Nodal signaling (Figure 4D,E, Video 4) and reduced expression of the pan-mesendodermal marker sebox (Figure 4—figure supplement 1E–G). Consistent with the notion that Nodal-mediated mesendoderm specification was reduced in DD-removed explants, we also found that the proportion of blastoderm explants forming an extension and the degree of explant extension were reduced (Figure 4F–I). Moreover, the extent of mesendoderm induction in DD-removed explants closely correlated with the degree of explant extension at bud stage (Figure 4—figure supplement 1H). This suggests that in blastoderm explants, maternal deposition of dorsal determinants in cells close to the future wounding site is required for Nodal signaling activation and mesendoderm induction, much more than in embryos that retain Nodal signals emanating from the YSL (Figure 4D,E, Figure 4—figure supplement 1E–H, Figure 4—figure supplement 2A–F). To exclude that the effect of β-cat depletion on Nodal signaling is a mere secondary consequence of embryo ventralization (Mizuno et al., 1999; Ober and Schulte-Merker, 1999), we analyzed Nodal signaling in blastoderm explants prepared from embryos ventralized by overexpression of a constitutively active (CA) form of the BMP receptor Alk8 (Payne et al., 2001). CA-Alk8-overexpressing blastoderm explants did not display a strong reduction in the establishment of a Nodal signaling domain (Figure 4—figure supplement 2G,H), supporting the notion that β-cat does not affect Nodal signaling through its effect on dorsoventral embryo patterning.
To further challenge the notion that maternal pre-patterning of blastoderm explants is critical for mesendoderm induction and morphogenesis, we sought to prepare explants lacking all pre-patterning. To this end, we prepared explants from the animal pole (AP) of the blastoderm, which should be devoid of marginal cells and dorsal determinants (Xu et al., 2014). In these explants, no accumulation of β-cat was observed in the nuclei of cells close to the wounding site (Figure 4—figure supplement 1C,D). Moreover, AP-derived explants showed strongly reduced nuclear pSMAD2/3 accumulation, mesendoderm induction and tissue extension (Figure 4D–H, Figure 4—figure supplement 1E–H, Video 4). These findings therefore suggest that maternal pre-patterning by dorsal determinants is needed for proper mesendoderm induction and morphogenesis in blastoderm explants.
To determine whether positional information retained within the blastoderm and subsequently wounding site of blastoderm explants is critical for coordinated mesendoderm induction and patterning, we dissociated blastoderm explants obtained from high stage embryos into their constituent cells and allowed them to reaggregate in culture in the presence of serum-free medium (Figure 4—figure supplement 3A). These reaggregates are thus expected to consist of a random mixture of blastomeres, including some marginal cells that are pre-patterned to express Nodal signals and, thus, induce mesendoderm. We found that blastomere reaggregates underwent progressive compaction and eventually rounded up in culture (Figure 4—figure supplement 3B). Importantly, reaggregated cells were still capable of undergoing cell divisions (Figure 4—figure supplement 3B), as evidenced by their progressive reduction in cell diameter (Figure 4—figure supplement 3C), suggesting that blastomere reaggregates are viable under culture conditions. Interestingly, we found that even when cultured until stage-matched embryos had reached bud stage, blastomere reaggregates failed to form a clearly recognizable extension (Figure 4—figure supplement 3D,E). Consistent with this, we found that expression of the mesodermal marker ntl (Schulte-Merker, 1995; Schulte-Merker et al., 1994) was restricted to a few sparse cells or seemingly randomly scattered cells within the reaggregates (Figure 4—figure supplement 3F–H). Importantly, neither increasing the cell number within reaggregates nor changing the culture media conditions substantially impacted on reaggregate morphogenesis or its ability to undergo coordinated mesoderm specification (Figure 4—figure supplement 3D–H). Collectively, these findings suggest that maternal pre-patterning and positional information within the blastoderm and, subsequently wounding site of blastoderm explants, is required for coordinated mesoderm induction and morphogenesis.
Our analysis of blastoderm explants and reaggregates so far indicates that in isolation from extraembryonic signaling sources, the spatially localized pre-patterning within the blastoderm margin can still drive a large degree of mesendoderm patterning and morphogenesis. This seems to contradict previous claims from embryo analysis that Nodal signaling from the extraembryonic YSL is needed for these processes (Carvalho and Heisenberg, 2010; Fan et al., 2007; Gagnon et al., 2018; Hong et al., 2011; van Boxtel et al., 2015; Veil et al., 2018; Xu et al., 2012). To understand this apparent discrepancy, we directly compared mesendoderm induction in blastoderm explants with intact embryos where Nodal signal production within the YSL is blocked. To block Nodal expression within the YSL, we injected sqt and cyc MOs specifically into the YSL of high stage embryos. Consistent with previous observations (Fan et al., 2007), we found that injection of high amounts of both sqt (2 ng/embryo) and cyc (4 ng/embryo) MOs impaired induction of endoderm (sox32) and head mesoderm (hgg and gsc) (Figure 5A, Figure 5—figure supplement 1A,B), the two mesendoderm progenitor cell types whose induction is thought to require the highest dose of Nodal signaling (Dougan et al., 2003; Gritsman et al., 2000; Schier et al., 1997; Thisse and Thisse, 1999). The induction of notochord (ntl), paraxial (papc) and ventrolateral mesoderm (myoD and tbx6), in contrast, appeared normal in the morphant embryos (Figure 5A, Figure 5—figure supplement 1A,B). This also became apparent by 32 hpf, when most morphant embryos displayed severe cyclopia, indicative of defective induction of anterior mesendoderm structures, but largely normal trunk and tail structures (Figure 5B). Together, these observations suggest that - consistent with previous findings (Fan et al., 2007) - the YSL is needed for robust induction of endoderm and head mesoderm, and that this signaling source might be specifically required to achieve peak levels of Nodal signaling along the blastoderm margin.
To address this possibility, we analyzed the number and intensity of nuclear pSMAD2/3 in cells along the blastoderm margin of embryos injected with control (6 ng/embryo) or sqt (2 ng/embryo) and cyc (4 ng/embryo) MOs into the YSL. We found a strong and persistent reduction of both the number and intensity of nuclear pSMAD2/3 positive cells along the blastoderm margin in sqt/cyc morphants (Figure 5C–G, Figure 5—figure supplement 1C), indicating that YSL-derived Nodal signals are needed to enhance Nodal signaling in the overlying blastoderm. Next, we asked whether injecting lower amounts of sqt/cyc MOs into the YSL would result in a less pronounced reduction in endoderm and head mesoderm specification. YSL-specific injection of lower amounts of sqt (1 or 0.4 ng/embryo) and cyc (2 or 0.8 ng/embryo) MOs led to a dose-dependent reduction in both the number of morphant embryos expressing the ppl marker hgg and the size of the ppl, as determined by the expression area of hgg, in those embryos at bud stage (Figure 6A–C). Likewise, by 32 hpf, the severity of the cyclopia phenotype closely correlated with the amount of sqt/cyc MOs injected within the YSL (Figure 5B, Figure 6—figure supplement 1A). Endoderm induction, visualized by sox32 expression (Kikuchi et al., 2001), was also strongly reduced in morphant embryos even with the lowest amount of YSL-injected sqt/cyc MOs (0.4 ng sqt + 0.8 ng cyc MO; Figure 6—figure supplement 1B,C). Consistent with such dose-dependent reduction in head mesoderm and endoderm specification, we found that Nodal signaling, assayed by pSMAD2/3 nuclear localization, was decreased to a lesser extent in embryos YSL-injected with low compared to high amounts of sqt/cyc MOs (2 ng sqt and 4 ng cyc MO) (Figure 6D–G, Figure 6—figure supplement 1D; compare to Figure 5C–G). Collectively, these findings indicate that robust induction of head mesoderm and endoderm structures in intact embryos critically relies on the right dosage of extraembryonic Nodal ligands.
Finally, to understand how the variable induction of head mesoderm in sqt/cyc MO morphants might relate to our observations in blastoderm explants, we performed a detailed analysis of the induction of these structures within the explants. Notably, we found that in explants the induction of head mesoderm and endoderm was highly variable (Figure 1D, Figure 1—figure supplement 1G,H), and both the size of the ppl and the number of sox32-positive endoderm cells were strongly reduced (Figure 6C, Figure 1—figure supplement 1G,H,J). This is consistent with our observations from the embryo analysis that Nodal signals from the extraembryonic YSL are needed for robust induction of these structures.
Discussion
The distinct ability of embryonic or stem cell aggregates to autonomously generate pattern and form has challenged the classical view that maternal factors and/or extrinsic biases, such as extraembryonic tissues, are essential for proper embryo patterning and morphogenesis (Mummery and Chuva de Sousa Lopes SM, 2013; Langdon and Mullins, 2011; Shahbazi et al., 2019; Simunovic and Brivanlou, 2017; Stern and Downs, 2012; Tam and Loebel, 2007; Turner et al., 2016). By establishing a novel model system to culture zebrafish embryos in the absence of the yolk cell and its associated structures, we were able to show that blastoderm explants can recapitulate a remarkable degree of patterning and morphogenesis reminiscent of intact embryos. Similar observations were also recently reported in a preprint (Trivedi et al., 2019), supporting the notion that zebrafish embryonic tissues can autonomously organize into the structures they are primed to form. Moreover, our findings provide insight into the specific contributions of maternal pre-patterning and extraembryonic tissues for robust development.
Previous studies have provided evidence that maternal factors, initially deposited at the vegetal pole of the oocyte, are required to activate β-cat at the future dorsal side of the embryo and, consequently, initiate Nodal signaling therein (Langdon and Mullins, 2011; Marlow, 2020). Our data support this view and show that this maternal activation of Nodal signaling at the dorsal germ ring margin is required for proper mesoderm induction in blastoderm explants. In particular, our findings that removing dorsal and marginal pre-patterning in AP blastoderm explants or interfering with the spatial distribution of these pre-patterning cues in reaggregates is sufficient to suppress coordinated mesoderm induction and morphogenesis, clearly shows that maternal pre-patterning is indispensable for blastoderm explant organization and morphogenesis. This differs from recent observations in human and mouse embryonic stem cell (hESCs and mESCs, respectively) colonies, where embryoid bodies (EBs) are thought to form by spontaneous symmetry-breaking and self-organization (Turner et al., 2016). Possibly, this difference in self-organizing capacity is due to blastoderm cells displaying reduced stemness when compared to ES cells, a plausible assumption given that blastoderm explants were taken from intact embryos that already underwent some degree of embryo polarization and cell fate restriction (Langdon and Mullins, 2011; Marlow, 2020). Furthermore, we cultured our blastoderm explants in media without serum or other growth factors in order to monitor blastoderm maturation in the absence of any extrinsic signals. In contrast, stem cell colonies or, most recently, also zebrafish explants were cultured in media supplemented with growth factors and/or serum, which might be required for self-organization to emerge in vitro (Baillie-Johnson et al., 2015; Beccari et al., 2018; ten Berge et al., 2008; Turner et al., 2017; Turner et al., 2014; Trivedi et al., 2019; van den Brink et al., 2014; Warmflash et al., 2014).
In zebrafish embryos, extraembryonic Nodal signals from the YSL are thought to be required for mesendoderm induction along the dorsal-ventral axis (Carvalho and Heisenberg, 2010; Fan et al., 2007; Gagnon et al., 2018; Hong et al., 2011; Schier and Talbot, 2005; van Boxtel et al., 2015; Veil et al., 2018; Xu et al., 2012). Our finding that blastoderm explants can not only induce the three germ layers, but also display a seemingly complete patterning of the mesendoderm anlage, seems to argue against such critical function for YSL-derived Nodal signals in mesendoderm patterning. However, closer inspection of the role of the YSL for mesendoderm patterning revealed some striking similarities between embryos with reduced YSL-derived Nodal signals and blastoderm explants: in both cases, induction of head mesoderm and endoderm, two tissues whose induction requires high doses of Nodal signaling (Dougan et al., 2003; Gritsman et al., 2000; Schier et al., 1997; Thisse and Thisse, 1999), is rather variable. This points at the intriguing possibility that Nodal signals from the extraembryonic YSL effectively boost Nodal signaling within the blastoderm margin, needed for robust induction of these two tissues. Yet, mesendoderm patterning defects in embryos with reduced Nodal signals from the YSL are generally more pronounced than in blastoderm explants, suggesting that the YSL is more important for robust induction of head mesoderm and endoderm in the intact embryo compared to explants. One possibility for such embryo-specific requirement of the YSL might be to induce sufficiently high doses of Nodal signaling along the entire blastoderm margin, a requirement less important in explants where the blastoderm margin has collapsed into a single region surrounding the wounding site.
While blastoderm explants were able to recapitulate a remarkable degree of embryo patterning and morphogenesis, there are still several differences between explants and embryos: foremost, mesendoderm progenitors do not undergo internalization in explants, but instead form an ‘exogastrula’, similar to previous observations in mouse gastruloids (Baillie-Johnson et al., 2015; Beccari et al., 2018; Turner et al., 2017; Turner et al., 2014; van den Brink et al., 2014). One possible explanation for this could be that the specific geometry of the embryo is critical for mesendoderm internalization. Alternatively, extraembryonic structures might form a permissive substrate for mesendoderm internalization. Notably, while mesendoderm progenitors fail to internalize in explants, they still form an extension, which - similar to body axis elongation in embryos - depends on Wnt/PCP-signaling (Čapek et al., 2019; Gray et al., 2011; Solnica-Krezel and Sepich, 2012; Tada and Heisenberg, 2012). This suggests that the cellular rearrangements underlying body axis extension in explants and body axis elongation in embryos might be conserved. Moreover, our and previous observations support that cell proliferation does not seem to play a major role in either of these processes (Liu et al., 2017; Quesada-Hernández et al., 2010; Zhang et al., 2008). Yet, whether volumetric growth by cell proliferation might function in later stages of explant maturation, similar to its role, for instance, in tail extension in chick and mouse embryos (Mongera et al., 2019), remains to be determined.
Collectively, our findings suggest that zebrafish blastoderm explants undergo what has been recently termed ‘genetically encoded self-assembly’, where heterogeneities among cells endowed with a similar genetically encoded cell or developmental program drive the organization of those cells, via local interactions, into the structures they were initially primed to form (Turner et al., 2016; Xavier da Silveira Dos Santos and Liberali, 2019). This interpretation implies that both intrinsic heterogeneities (e.g. maternal factors) and extraembryonic signals might not only be required to render embryonic development predictable and robust, but also to drive the emergence of embryo organization. Future work will have to address to what extent stem cell colonies/organoids indeed undergo bona fide self-organization, and whether some of their organization might in fact be due to genetically encoded self-assembly.
Materials and methods
Fish lines and husbandry
Request a detailed protocolZebrafish (Danio rerio) handling was performed as described (Westerfield, 2000). The mutant and transgenic lines used in this study are detailed in the key resources table. Embryos were raised at 25–31°C in Danieau’s solution or E3 medium and staged according to Kimmel et al. (1995).
Blastoderm explant preparation
Request a detailed protocolEmbryos for blastoderm explant preparation were manually dechorionated with watchmaker forceps. Explants were then prepared at 256 cell stage, unless stated otherwise, by removing the entire blastoderm from the yolk cell using forceps, and then washed by gentle pipetting to remove residual yolk granules (schematic representation in Figure 1A). Animal pole explants were prepared by removing only the animal portion of the blastoderm at 256 cell stage using a hair knife instead of forceps. All explants were cultured in Ringer’s solution (116 mM NaCl, 2.9 mM KCl, 1.8 m M CaCl2, 5.0 mM HEPES, pH 7.2) or Danieau’s medium (58 mM NaCl, 0.7 mM KCl, 0.4 mM MgSO4, 0.6 mM Ca(NO3)2, 5 mM HEPES, pH 7.6). Approximately 1 hour (hr) after explant preparation, explants with delayed cell cleavages, as assessed by increased cell size, were removed from the dish. The remaining explants were then staged based on sibling embryos from the same egg lay. To analyze explant morphology, bright-field side view images of blastoderm explants were acquired using a stereo-microscope (Olympus SZX 12) equipped with a QImaging Micropublisher 5.0 camera.
Embryo microinjections
Request a detailed protocolmRNAs were synthesized using the mMessage mMachine Kit (Ambion). Injections at 1 cell stage were carried out as described (Westerfield, 2000) and the plasmids used for mRNA production are detailed in the key resources table. Labelling the interstitial fluid was performed by injecting 2 ng of Dextran Alexa Fluor 647, or Dextran Cascade Blue, between the deep cells of blastoderm explants corresponding to high stage or sphere stage embryos, as described in Krens et al. (2017). For time-lapse imaging of the interstitial fluid, embryos were also injected at 1 cell stage with 80 pg of Membrane-GFP. For high-resolution live imaging of gastrulating blastoderm explants, Tg(gsc::EGFP-CAAX) expressing embryos were injected at 1 cell stage with 50 pg of H2A-chFP to label the cell nuclei. For the clone dispersal analysis, the plasma membrane of all cells was labelled by injecting 80 pg of Membrane-RFP mRNA at 1 cell stage. 6.5 pg of H2A-chFP mRNA was subsequently injected into a single marginal blastomere at 128 cell stage (schematic representation in Figure 2—figure supplement 2A). To constitutively activate BMP signaling, 30 pg of CA-Alk8 mRNA were injected together with 60 pg of Membrane-GFP mRNA at 1 cell stage.
YSL injections were performed right after its formation, between 1 k and high stage. As a control for homogeneous YSL injections, 50 pg of H2B-EGFP mRNA and 0.2% Phenol Red were co-injected into the YSL. To block YSL-derived Nodal signals, 0.4–2 ng of sqt MO (sqt MO: 5′-ATGTCAAATCAAGGTAATAATCCAC-3′ [Feldman and Stemple, 2001]) were co-injected with 0.8–4 ng of cyc MO (cyc MO: 5′-GCGACTCCGAGCGTGTGCATGATG-3′ [Karlen and Rebagliati, 2001]) into the YSL. As a control for MO injections, 6 ng of a standard negative control MO (human β-globin MO: 5′-CCTCTTACCTCAGTTACAATTTATA-3′) were also injected.
Live imaging of blastoderm explants
Request a detailed protocolBright-field time-lapse imaging of explants was performed using a Zeiss LSM 800 upright microscope equipped with a Zeiss Plan-Apochromat 20x/1.0 water immersion objective starting when stage-matched control embryos had reached sphere-early dome stage. High-resolution live imaging of developing blastoderm explants was performed using a LaVision Trim 2-photon microscope equipped with a Zeiss Plan-Apochromat 16x/0.8 water immersion objective and a Ti:Sa laser (Chameleon, Coherent) set to 830 nm and an OPO laser set at 1100 nm. In this case, imaging was started when stage-matched control embryos had reached 50% epiboly-germ ring stage. In both cases, blastoderm explants were oriented in a side view, and the acquired images were then processed using Fiji and/or Imaris.
Analysis of fluid accumulation, YSL formation and EVL differentiation
Request a detailed protocolImaging of the explant interstitial fluid at bud stage was performed using a LaVision Trim 2-photon microscope equipped with a Zeiss Plan-Apochromat 20x/1.0 water immersion objective and a Ti:Sa laser (Chameleon, Coherent) set to 800 nm and an OPO laser set at 1100 nm, or a Zeiss LSM 880 upright confocal microscope equipped with a Zeiss Plan-Apochromat 20x/1.0 water immersion objective. Time-lapse imaging of the explant interstitial fluid was performed using a Zeiss LSM 800 upright microscope equipped with a Zeiss Plan-Apochromat 20x/1.0 water immersion objective, starting when stage-matched control embryos had reached high-oblong stage. In both cases, blastoderm explants were oriented in a side view, and the acquired images were then processed using Fiji and/or Imaris.
Imaging of YSL-like features in blastoderm explants was performed when corresponding stage-matched control embryos reached 50% epiboly-germ ring stage using a LaVision Trim 2-photon microscope as described above. Explants were oriented with the wounding site facing the objective (top view), and the acquired images were then processed using Fiji.
To analyze EVL differentiation in blastoderm explants, both animal-pole oriented embryos and explants were imaged in parallel on a Zeiss LSM 880 upright confocal microscope equipped with a Zeiss Plan-Apochromat 20x/1.0 water immersion objective. To assess whether the entire explant was covered by EVL, explants and embryos were then fixed at bud stage in 4% paraformaldehyde (PFA) at 4°C overnight, washed 5x with PBS and then imaged on a Zeiss LSM 880 upright confocal microscope equipped with a Zeiss Plan-Apochromat 20x/1.0 water immersion objective. The acquired images for both embryos (side view) and explants (side view) were then processed using Fiji and/or Imaris (Bitplane).
Sample preparation for live imaging
Request a detailed protocolDechorionated embryos were mounted in 2% agarose molds on petri dishes and immobilized in 0.7% low melting point agarose. Explants were mounted in smaller agarose molds (800 × 800 µm from Microtissues) and cultured in Ringer’s or Danieau’s solution.
Whole mount in situ hybridization (WMISH)
Request a detailed protocolEmbryos were fixed with 4% PFA overnight at 4°C. Antisense RNA probes were synthesized using SP6, T7 or T3 RNA polymerase from mMessage mMachine kits (ThermoFisher, AM1344) with Roche digoxigenin (DIG)-modified nucleotides, or fluorescein (FITC)-labeled nucleotides, from partial cDNA sequences. WMISHs were performed as previously described (Thisse and Thisse, 2008) with the following modifications: (i) proteinase K treatment was performed only for sox32 in situ hybridization and for double WMISHs; (ii) 5% dextran sulfate was added to the hybridization solution for sox32, mxtx2, tbx16 and bmp2b probes and all double WMISHs; and, (iii) 0.05x SSC was used for more stringent washes for sqt in situ hybridizations. Single WMISHs were performed with DIG-labeled antisense RNA probes, while for double WMISHs, a mixture of DIG- and FITC-labeled RNA probes were added simultaneously during the hybridization step.
For double WMISHs, hybridization and visualization of DIG-labelled probes with the chromogenic NBT/BCIP substrate were performed as previously described (Thisse and Thisse, 2008). When the signal was sufficiently intense, the staining solution was replaced with a solution of 100 mM glycine-HCl pH2.2 for 10 min at room temperature as described in Denker et al. (2008). The samples were then washed several times with fresh PBS + 0.1% (w/v) Tween (PBST) and blocked in blocking buffer (1 × PBST, 2% sheep serum (v/v), 2 mg/ml BSA) for 2 hr at room temperature and incubated overnight at 4°C with Anti-Fluorescein-AP Fab fragments (1:500) in fresh blocking buffer. Afterwards, the samples were washed 6x for 10–20 min at room temperature in PBT and 3x in 0.1 M Tris-HCl pH8.2 + 0.1% (w/v) Triton. The staining solution for red staining was prepared using SIGMAFAST Fast Red TR/Naphthol AS-MX Tablets according to the manufacturer’s protocol and added to the samples. After staining, the samples were cleared for 10 min in 100% EtOH and washed several times in PBT. Blastoderm explants and corresponding stage-matched embryos were always processed in the same tube, and the staining reaction was therefore stopped simultaneously. All WMISHs were imaged on a stereo-microscope (Olympus SZX 12) equipped with a QImaging Micropublisher 5.0 camera.
Whole mount immunofluorescence (WMIF)
Request a detailed protocolα-pSMAD2/3 and α-β-catenin whole mount immunofluorescence was performed as described previously (van Boxtel et al., 2015). In short, dechorionated embryos and explants were fixed in 4% PFA (in PBS) overnight at 4°C, washed in PBS and then directly transferred into 100% MeOH and stored at −20°C for, at least, 2 hr. The samples were then washed 5x in PBSTr (PBS + 1% (w/v) Triton X-100), blocked 2x in blocking solution (PBSTr + 10% goat serum + 1% DMSO) for 2 hr and then incubated overnight at 4°C with α-pSMAD2/3 antibody (1:1000) and/or α-β-catenin antibody (1:500) diluted in blocking solution. Afterwards, the samples were washed 5x in PBSTr for 10 min, washed 3x in PBS + 0.1% (w/v) Triton for 1 hr and then incubated overnight at 4°C in secondary antibody (1:500) diluted in blocking solution. For labelling cell nuclei, the samples were incubated in DAPI (1:1000) diluted in blocking solution for 18 min at room temperature. After secondary antibody incubation, the samples were washed 2x for 5 min in PBSTr and 5x for 10–20 min in PBS + 0.1% Triton.
α-pSMAD1/5 whole mount immunofluorescence was performed as described previously (Pomreinke et al., 2017) with the following modifications: (i) for re-hydration, embryos and explants were directly transferred to PBST; (ii) Anti-Phospho-Smad1/5 (Ser463/465) (1:100) was used to detect pSmad1/5; (iii) secondary antibody incubation (1:500) was performed overnight at 4°C; and (iv) labelling of cell nuclei with DAPI was performed as described above. Explants and corresponding stage-matched control embryos were always processed in the same tube.
Immunostained embryos and blastoderm explants were imaged using an upright confocal microscope (Zeiss LSM 880 or a Zeiss LSM 800 upright microscope equipped with a Zeiss Plan-Apochromat 20x/1.0 water immersion objective) and mounted in 2% agarose molds and subsequently covered with PBS. Embryos and explants were oriented using the Tg(gsc::EGFP-CAAX) line to identify the dorsal side (embryos) or wounding site (explants). Imaging conditions were kept similar between different samples and replicates of the same experiment.
Nodal inhibitor treatment
Request a detailed protocolBlastoderm explants together with stage-matched control embryos were incubated in 0.1% DMSO or in 50 µM of the reversible Nodal inhibitor SB-505124 (DaCosta Byfield et al., 2004; Fan et al., 2007; Hagos and Dougan, 2007; Hagos et al., 2007; Rogers et al., 2017; van Boxtel et al., 2015; Vogt et al., 2011). For inhibitor washout experiments, embryos and explants were washed 5x in Danieau’s solution for 5 min.
Cell division inhibition (HUA treatment)
Request a detailed protocolTo block cell division events during explant extension, blastoderm explants were incubated in 0.1% DMSO or a cocktail of 60 mM Hydroxyurea and 300 μM Aphidicolin (HUA) from shield to bud stage, as determined by stage-matched control embryos. To assess explant morphology, bright-field side view images of bud stage blastoderm explants treated with DMSO or HUA were then acquired using a stereo-microscope (Olympus SZX 12) equipped with a QImaging Micropublisher 5.0 camera.
Removal of dorsal determinants
Request a detailed protocolRemoval of dorsal determinants was performed as described (Mizuno et al., 1999; Ober and Schulte-Merker, 1999). In short, embryos were dechorionated immediately after egg laying in Danieau’s solution and the vegetal-most part of the yolk cell was removed using a hair knife. The whole procedure was carried out within 15 min after the first eggs were laid. Manipulated embryos were then allowed to heal for 10 min before being transferred to fresh Danieau’s medium. Dorsal determinants removed embryos were used as staging controls for the corresponding explants.
Assessment of mesendoderm induction
Request a detailed protocolFor assessing changes in mesoderm induction upon Nodal inhibitor treatment, after dorsal determinants removal and in animal pole explants, corresponding explants were prepared from Tg(sebox::EGFP) transgenic embryos to mark mesendoderm progenitors. At bud stage, these explants together with control explants were fixed in 4% PFA at 4°C overnight and washed 5x in PBS. The samples were then imaged on a Zeiss LSM 800 Upright microscope equipped with a Zeiss Plan-Apochromat 20x/1.0 water immersion objective. Imaging conditions were kept similar between all samples and replicates.
Blastomere reaggregate preparation
Request a detailed protocolFor dissociation, manually dechorionated embryos were transferred to 1 ml of dissociation buffer (calcium-free Ringer’s + 5 mM EGTA for reaggregates cultured in Ringer’s solution and L15 + 5 mM EGTA for reaggregates cultured in L15) at high stage. The embryos were then mechanically dissociated by pipetting and washed 3x with the corresponding culture medium, followed by centrifugation at 200 g for 3 min to remove residual yolk proteins and EGTA. The dissociated cells were then seeded in agarose microwells cast from PDMS molds (Microtissues, diameter 800 µm, depth 800 µm) at a density of one or three dissociated embryo(s)/microwell. After a settlement period of around 30 min, the petri dish with dissociated cells was filled with the respective culture medium (schematic representation in Figure 4—figure supplement 3A). Blastomere reaggregates were analyzed for their morphology and fixed for in situ hybridization, when control embryos had reached bud stage.
To assess whether cells keep dividing in culture - as a proxy for their health status and continued development - small reaggregates (one embryo/microwell), prepared from a mix of embryos with labelled nuclei and cell membranes (1 cell stage injection of 80–100 pg Membrane-GFP and 50 pg of H2A-chFP) and unlabeled embryos, were cultured in Ringer’s solution. These reaggregates were then imaged at a Zeiss LSM 880 upright confocal microscope equipped with a Zeiss Plan-Apochromat 20x/1.0 water immersion objective shortly after reaggregation. The size of 15–30 randomly selected cells on the inside of the reaggregate was measured every 32 min using the Line tool in Fiji.
Analysis of explant morphology
Request a detailed protocolTo quantify the morphology of blastoderm explants, bright-field side view images of explants were acquired as described above. Explants were considered as extended when a clear indentation was observed between the round ‘back’ and extended ‘tip’ of the explant. Explant circularity was then determined by manually outlining the shape of the explant using the Circularity plugin from Fiji. All explants, independently of whether they formed a clearly recognizable extension or not, were pooled for the circularity analysis.
To quantify the relative length of the explant extension, the extension length starting where a clear change in curvature was detectable was measured using the Segmented Line tool from Fiji and then normalized to the total explant length. Only explants with a recognizable extension were included in this analysis. Note that in this analysis, we only included blastoderm explants where the extension could be clearly delineated (except for Figure 2—figure supplement 3D and Figure 4—figure supplement 1H, where not-extended explants were included and their normalized extension length was considered to be 0). To ensure that potential defects in axis extension were indeed due to the different pharmacological and genetic perturbations and not embryo batch variability, only experiments in which at least 70% of control explants showed clearly recognizable extensions were considered for further analysis.
To analyze the change in area during explant extension, side view images of explants at the indicated time points were manually outlined using the freehand tool (schematic representation in Figure 3—figure supplement 1A). The width of the explants was measured at the point where the explant was thickest using the Line tool from Fiji (schematic representation in Figure 3—figure supplement 1B). The length of the extension over time was measured as described above (schematic representation in Figure 3—figure supplement 1C).
Analysis of explant cavity localization
Request a detailed protocolTo quantify the average location of the explant cavity, bright-field side view images of explants were acquired as described above. The segmented line tool in Fiji was used to determine the distances from the tip of the extension to the closest edge of the cavity, and from the opposite edge of the cavity to the other end of the explant. The measured distances were then normalized to the total explant length, as schematically represented in Figure 1—figure supplement 1B. If an explant formed more than one cavity, the location of the biggest lumen was used for this quantification.
Analysis of cell size
Request a detailed protocolCell size in blastoderm explants and stage-matched control embryos was measured using the Line tool in Fiji. For this, explants and embryos were fixed at various developmental stages (high stage, sphere, 30% epiboly, 50% epiboly, germ ring and shield stage) and immunostained for β-catenin to outline cell-cell contacts. Approximately 30 randomly selected deep cells from both explants and embryos were measured for this analysis.
Analysis of gene expression domains
Request a detailed protocolFor this analysis, only blastoderm explants expressing the gene of interest, as determined by in situ hybridization, and where the back-tip axis could be clearly delineated, were analyzed. To quantify the location of a gene expression domain along the back-tip axis of the explant, the distances from the tip of the extension to the start and end of the expression domain were measured and normalized by the total explant length using the Segmented Line tool in Fiji. The normalized length of the explant was then binned in 10 windows and the presence (1) or absence (0) of expression within a certain bin was plotted in a binary fashion and, ultimately, averaged across explants and experimental replicates. For six3-positive explants, which in some cases exhibited more than one expression domain, the measurement was performed for each expression domain separately. For hgg-positive explants, only samples with a single ppl were considered for this analysis.
To compare the size of the hgg expression domains in different experimental conditions, the boundaries of the expression domain, as determined by in situ hybridization, was manually outlined and measured in Fiji. Animal or side views were used respectively for intact embryos and blastoderm explants, and only samples with a single ppl were considered for this analysis.
Analysis of EVL differentiation
Request a detailed protocolAnalysis of early EVL differentiation was performed using Tg(krt4::EGFP-CAAX) embryos, expressing EGFP in differentiating EVL cells. Analysis of EGFP expression in surface cells of both stage-matched embryos and explants , was used as a readout for EVL differentiation and this analysis was started shortly after the first visible accumulation of EGFP in the intact embryo. For quantifying EGFP accumulation, the same number of z planes were SUM projected in explants and embryos using Fiji. A ROI of 400.90 μm2 was drawn in the center of the samples, and EGFP mean intensity was measured every 15 min. For background subtraction, an ROI of 24.91 × 24.91 μm2 was analyzed outside of the sample.
Analysis of pSMAD2/3 nuclear accumulation
Request a detailed protocolFor a quantitative analysis of pSMAD2/3 nuclear accumulation, randomly selected embryos and stage-matched blastoderm explants exhibiting pSMAD2/3 positive nuclei were used. The Imaris Spot detection plugin was used to automatically determine the 3D coordinates of all nuclei based on the DAPI signal. Nuclei belonging to the EVL, undergoing cell division or displayed inhomogeneous and/or barely visible DAPI signal, were manually excluded from the analysis. Moreover, only nuclei within a Z-volume of 100 µm were considered for this analysis in order to avoid depth-related intensity changes. The mean fluorescence intensities of pSMAD2/3 and DAPI as well as the corresponding nuclei 3D coordinates were extracted using the Spot Detection plugin. As a reference point to calculate nuclei distribution within the explants or embryos, the wounding site of the explant, detectable by increased and irregular junctional staining of β-catenin and/or comparatively lower and more irregular nuclei density, or the YSL nuclei in the intact embryo were manually determined using the same plugin and, subsequently, extracted. To determine the distance of each nucleus to the wounding site or YSL, a custom MATLAB script (Source code 1) was used to project all nuclei into a 2D plane along a manually defined axis. For embryos, all nuclei were projected along their Z axis. For explants, an additional reference point was manually provided and used to re-position the wounding site in the center of the image. The geometric distance of all nuclei to the closest YSL/wounding site spot in this 2D plane was then automatically calculated using the aforementioned script.
To determine the number of pSMAD2/3 positive nuclei in blastoderm explants and embryos, the following criteria were used: (i) only nuclei within 150 µm from the reference point (YSL nuclei for embryos, wounding site for explants) were considered in order to avoid imaging artifacts often arising at the sample edges; (ii) following background subtraction, only nuclei with positive nuclear pSMAD2/3 mean intensity were considered for further analysis; (iii) nuclear pSMAD2/3 mean intensity was normalized using DAPI mean intensity to correct depth-related artifacts; and (iv) only nuclei with a normalized pSMAD2/3 to DAPI ratio above 0.1 were considered to be positive. For background subtraction, the mean nuclear pSMAD2/3 intensity between 120–150 µm of the YSL or wounding site was used.
The brightest pSMAD2/3 positive nuclei were defined as the top 30% brightest normalized pSMAD2/3 to DAPI ratio of all positive nuclei, which typically correspond to the first 1–2 cell rows in the intact embryo. The brightest pSMAD2/3 positive nuclei in stage-matched explants and embryos were used for intensity comparison and for calculating the distance to the YSL or wounding site. To plot this distance as cell tiers, we considered the average cell size measured at 50% epiboly, germ ring and shield stage, as described above (see Figure 1—figure supplement 1C).
To classify pSMAD2/3 nuclear accumulation in blastoderm explants and embryos at 50% epiboly as present, strongly reduced or absent, we manually counted the number of bright pSMAD2/3 positive nuclei using Imaris. Wildtype blastoderm explants and embryos had an average of ±22 or ±81 very bright pSMAD2/3 positive nuclei, respectively. A strongly reduced pSMAD2/3 domain was defined as a 50% reduction in the number of bright pSMAD2/3 nuclei as compared to the average number detected at this stage in control explants and embryos. Samples were defined as negative when no pSMAD2/3 nuclei could be detected.
Analysis of cell clone dispersal
Request a detailed protocolTo analyze cell clone dispersal, explants expressing Membrane-RFP in all cells and H2A-chFP in cell clones were imaged on a Zeiss LSM 880 upright confocal microscope equipped with a Zeiss Plan-Apochromat 20x/1.0 water immersion objective. The Imaris Spot detection plugin was used to automatically determine the 3D coordinates of all H2A-chFP-expressing nuclei at the start of acquisition (typically corresponding to late sphere stage) and after 1 and 2 hr. As proxies for cell dispersion, both the distance between all labelled nuclei (averaged over all pairs of cells), and the average distance between each cell and the center of mass of the labelled cell cluster were computed for each time point.
Analysis of nuclear β-catenin positive nuclei
Request a detailed protocolNon-dividing β-catenin positive nuclei were manually selected using the Imaris Spot detection plugin based on their nuclear β-catenin signal. Projection and calculation of distances were performed as described for pSMAD2/3 nuclear immunostainings. To plot the number of β-catenin positive nuclei as a function of the distance to the wounding site, expressed as cell tiers, we measure the diameter of 10 randomly selected cells from each sample using the Line tool using Fiji as described above.
To analyze the clustering of β-catenin positive nuclei, the angular distribution of all positive nuclei in reference to an axis drawn between the wounding site and the brightest β-catenin nuclei was automatically calculated (note that these nuclei were then excluded from the subsequent angle analysis), using the same MATLAB script described for analysis of pSMAD2/3 immunostainings (Source code 1).
Analysis of the sebox::EGFP-expressing domain
Request a detailed protocolTo assess the extent of mesendoderm induction, as assayed by the pan-mesendodermal marker sebox, the same number of z planes were MAX projected in explants of different experimental conditions (Nodal inhibitor treatments, DD-removed and animal pole explants) and corresponding controls using Fiji. To compare the size of the mesendodermal domain, the freehand tool was used to manually outline the sebox::EGFP-positive area. For this analysis, only explants, where the complete sebox::EGFP was detectable, were included. Based on this selection, explants with a domain 50% smaller than the average area in control explants were categorized as reduced mesendoderm induction. Explants, which did not show any obvious EGFP accumulation, were considered to be negative, and the corresponding area of the sebox::EGFP domain was considered to be 0.
Statistics
The number of embryos, blastoderm explants or blastoderm reaggregates analyzed (n) and the number of experimental replicates (N) are indicated in the figure legends. No statistical tests were used to assess sample sizes. No inclusion/exclusion criteria were used and all analyzed samples were included in the analysis. In the graphs, the error bars shown correspond to ± s.d., except in Figure 1E and Figure 1—figure supplement 1I where the error bars denote ± s.e.m. In box plots, the error bars correspond to the minimum and maximum of the dataset distribution. The statistical test used to assess significance is stated in the figure legends and was chosen after the distribution of each group was tested using the Shapiro-Wilk normality test. To compare two groups, we used either a two-tailed Student t test or a two-tailed Mann-Whitney test, depending on whether the dataset shows normal distribution. When using the t test, the variances were assessed using the F-test. To compare more than two groups, we used either ANOVA or Kruskal-Wallis test, depending on whether the dataset shows normal distribution. In these cases, a correction for multiple comparisons was used to increase statistical power. All statistical analyses were performed using GraphPad Prism.
Data availability
All data generated and analyzed in the manuscript are provided as Source data Files. The Custom Script used has been uploaded as Source code 1.
References
-
The yolk syncytial layer in early zebrafish developmentTrends in Cell Biology 20:586–592.https://doi.org/10.1016/j.tcb.2010.06.009
-
The role of the yolk syncytial layer in germ layer patterning in zebrafishDevelopment 127:4681–4689.
-
Role of frizzled 7 in the regulation of convergent extension movements during gastrulation in Xenopus laevisDevelopment 127:3091–3100.
-
The in vitro development of blastocyst-derived embryonic stem cell lines: formation of visceral yolk sac, blood islands and myocardiumJournal of Embryology and Experimental Morphology 87:27–45.
-
Temporal dynamics in the formation and interpretation of nodal and BMP morphogen gradientsCurrent Topics in Developmental Biology 137:363–389.https://doi.org/10.1016/bs.ctdb.2019.10.012
-
The primary role of zebrafish nanog is in extra-embryonic tissueDevelopment 145:dev147793.https://doi.org/10.1242/dev.147793
-
The role of maternal Activin-like signals in zebrafish embryosDevelopmental Biology 309:245–258.https://doi.org/10.1016/j.ydbio.2007.07.010
-
Back and forth between cell fate specification and movement during vertebrate gastrulationCurrent Opinion in Genetics & Development 18:311–316.https://doi.org/10.1016/j.gde.2008.07.011
-
Essential role of MARCKS in cortical actin dynamics during gastrulation movementsJournal of Cell Biology 164:169–174.https://doi.org/10.1083/jcb.200310027
-
The role of ppt/Wnt5 in regulating cell shape and movement during zebrafish gastrulationMechanisms of Development 120:467–476.https://doi.org/10.1016/S0925-4773(03)00004-2
-
Stages of embryonic development of the zebrafishDevelopmental Dynamics 203:253–310.https://doi.org/10.1002/aja.1002030302
-
Molecular regulation of pancreas development in zebrafish, methods in cell biologyMethods in Cell Biology 100:261–280.https://doi.org/10.1016/B978-0-12-384892-5.00010-4
-
Maternal and zygotic control of zebrafish dorsoventral axial patterningAnnual Review of Genetics 45:357–377.https://doi.org/10.1146/annurev-genet-110410-132517
-
Setting up for gastrulation in zebrafishCurrent Topics in Developmental Biology 136:33–83.https://doi.org/10.1016/bs.ctdb.2019.08.002
-
Mechanics of anteroposterior Axis formation in vertebratesAnnual Review of Cell and Developmental Biology 35:259–283.https://doi.org/10.1146/annurev-cellbio-100818-125436
-
Handbook of Stem Cells139–154, Differentiation in early development, Handbook of Stem Cells, Elsevier Inc, 10.1016/B978-0-12-385942-6.00014-7.
-
Functional characterization and genetic mapping of alk8Mechanisms of Development 100:275–289.https://doi.org/10.1016/S0925-4773(00)00541-4
-
Zebrafish gastrulation: putting fate in motionCurrent Topics in Developmental Biology 136:343–375.https://doi.org/10.1016/bs.ctdb.2019.10.009
-
Nodal and BMP dispersal during early zebrafish developmentDevelopmental Biology 447:14–23.https://doi.org/10.1016/j.ydbio.2018.04.002
-
The one-eyed pinhead gene functions in mesoderm and endoderm formation in zebrafish and interacts with no tailDevelopment 124:327–342.
-
Nodal morphogensCold Spring Harbor Perspectives in Biology 1:a003459.https://doi.org/10.1101/cshperspect.a003459
-
Molecular genetics of Axis formation in zebrafishAnnual Review of Genetics 39:561–613.https://doi.org/10.1146/annurev.genet.37.110801.143752
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
No tail (ntl) is the zebrafish homologue of the mouse T (Brachyury) geneDevelopment 120:1009–1015.
-
The zebrafish no tail geneSeminars in Developmental Biology 6:411–415.https://doi.org/10.1016/S1044-5781(06)80005-8
-
Initiation of convergence and extension movements of lateral mesoderm during zebrafish gastrulationDevelopmental Dynamics 234:279–292.https://doi.org/10.1002/dvdy.20507
-
A 3D model of a human epiblast reveals BMP4-driven symmetry breakingNature Cell Biology 21:900–910.https://doi.org/10.1038/s41556-019-0349-7
-
Friction forces position the neural anlageNature Cell Biology 19:306–317.https://doi.org/10.1038/ncb3492
-
Gastrulation: making and shaping germ layersAnnual Review of Cell and Developmental Biology 28:687–717.https://doi.org/10.1146/annurev-cellbio-092910-154043
-
The hypoblast (visceral endoderm): an evo-devo perspectiveDevelopment 139:1059–1069.https://doi.org/10.1242/dev.070730
-
Gene function in mouse embryogenesis: get set for gastrulationNature Reviews. Genetics 8:368–881.https://doi.org/10.1038/nrg2084
-
Antivin, a novel and divergent member of the TGFbeta superfamily, negatively regulates mesoderm inductionDevelopment 126:229–240.
-
Formation of the vertebrate embryo: moving beyond the spemann organizerSeminars in Cell & Developmental Biology 42:94–102.https://doi.org/10.1016/j.semcdb.2015.05.007
-
The specificities of small molecule inhibitors of the tgfß and BMP pathwaysCellular Signalling 23:1831–1842.https://doi.org/10.1016/j.cellsig.2011.06.019
-
Planar cell polarity and the developmental control of cell behavior in vertebrate embryosAnnual Review of Cell and Developmental Biology 28:627–653.https://doi.org/10.1146/annurev-cellbio-092910-154208
-
BookThe Zebrafish Book. A Guide for the Laboratory Use of Zebrafish (Danio Rerio)Eugene, United States: Oregon Press.
-
From single cells to tissue self-organizationThe FEBS Journal 286:1495–1513.https://doi.org/10.1111/febs.14694
-
Nanog-like regulates endoderm formation through the Mxtx2-Nodal pathwayDevelopmental Cell 22:625–638.https://doi.org/10.1016/j.devcel.2012.01.003
Article and author information
Author details
Funding
H2020 European Research Council (MECSPEC742573)
- Carl-Philipp Heisenberg
Austrian Academy of Sciences
- Alexandra Schauer
European Molecular Biology Organization (850-2017)
- Diana Pinheiro
Human Frontier Science Program (LT000429/2018-L2)
- Diana Pinheiro
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Ethics
Animal experimentation: All animal experiments in this study were performed in strict accordance with the guidelines of the Ethics and Animal Welfare Committee (ETK) in Austria. The respective approval number that covers the performed experiments is 66.018/0010-WF/II/3b/2014.
Copyright
© 2020, Schauer 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
-
- 6,410
- views
-
- 903
- downloads
-
- 47
- 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
- Genetics and Genomics
Paternal obesity has been implicated in adult-onset metabolic disease in offspring. However, the molecular mechanisms driving these paternal effects and the developmental processes involved remain poorly understood. One underexplored possibility is the role of paternally induced effects on placenta development and function. To address this, we investigated paternal high-fat diet-induced obesity in relation to sperm histone H3 lysine 4 tri-methylation signatures, the placenta transcriptome, and cellular composition. C57BL6/J male mice were fed either a control or high-fat diet for 10 weeks beginning at 6 weeks of age. Males were timed-mated with control-fed C57BL6/J females to generate pregnancies, followed by collection of sperm, and placentas at embryonic day (E)14.5. Chromatin immunoprecipitation targeting histone H3 lysine 4 tri-methylation (H3K4me3) followed by sequencing (ChIP-seq) was performed on sperm to define obesity-associated changes in enrichment. Paternal obesity corresponded with altered sperm H3K4me3 at promoters of genes involved in metabolism and development. Notably, altered sperm H3K4me3 was also localized at placental enhancers. Bulk RNA-sequencing on placentas revealed paternal obesity-associated sex-specific changes in expression of genes involved in hypoxic processes such as angiogenesis, nutrient transport, and imprinted genes, with a subset of de-regulated genes showing changes in H3K4me3 in sperm at corresponding promoters. Paternal obesity was also linked to impaired placenta development; specifically, a deconvolution analysis revealed altered trophoblast cell lineage specification. These findings implicate paternal obesity effects on placenta development and function as one potential developmental route to offspring metabolic disease.
-
- Developmental Biology
Although the impact of gender-affirming hormone therapy (GAHT) on spermatogenesis in trans women has already been studied, data on its precise effects on the testicular environment is poor. Therefore, this study aimed to characterize, through histological and transcriptomic analysis, the spermatogonial stem cell niche of 106 trans women who underwent standardized GAHT, comprising estrogens and cyproterone acetate. A partial dedifferentiation of Sertoli cells was observed, marked by the co-expression of androgen receptor and anti-Müllerian hormone which mirrors the situation in peripubertal boys. The Leydig cells also exhibited a distribution analogous to peripubertal tissue, accompanied by a reduced insulin-like factor 3 expression. Although most peritubular myoid cells expressed alpha-smooth muscle actin 2, the expression pattern was disturbed. Besides this, fibrosis was particularly evident in the tubular wall and the lumen was collapsing in most participants. A spermatogenic arrest was also observed in all participants. The transcriptomic profile of transgender tissue confirmed a loss of mature characteristics - a partial rejuvenation - of the spermatogonial stem cell niche and, in addition, detected inflammation processes occurring in the samples. The present study shows that GAHT changes the spermatogonial stem cell niche by partially rejuvenating the somatic cells and inducing fibrotic processes. These findings are important to further understand how estrogens and testosterone suppression affect the testis environment, and in the case of orchidectomized testes as medical waste material, their potential use in research.