Ectodermal Wnt signaling, cell fate determination, and polarity of the skate gill arch skeleton
Abstract
The gill skeleton of cartilaginous fishes (sharks, skates, rays, and holocephalans) exhibits a striking anterior–posterior polarity, with a series of fine appendages called branchial rays projecting from the posterior margin of the gill arch cartilages. We previously demonstrated in the skate (Leucoraja erinacea) that branchial rays derive from a posterior domain of pharyngeal arch mesenchyme that is responsive to Sonic hedgehog (Shh) signaling from a distal gill arch epithelial ridge (GAER) signaling centre. However, how branchial ray progenitors are specified exclusively within posterior gill arch mesenchyme is not known. Here, we show that genes encoding several Wnt ligands are expressed in the ectoderm immediately adjacent to the skate GAER, and that these Wnt signals are transduced largely in the anterior arch environment. Using pharmacological manipulation, we show that inhibition of Wnt signalling results in an anterior expansion of Shh signal transduction in developing skate gill arches, and in the formation of ectopic anterior branchial ray cartilages. Our findings demonstrate that ectodermal Wnt signalling contributes to gill arch skeletal polarity in skate by restricting Shh signal transduction and chondrogenesis to the posterior arch environment and highlights the importance of signalling interactions at embryonic tissue boundaries for cell fate determination in vertebrate pharyngeal arches.
Editor's evaluation
Building on previous work dissecting the polarisation of a gill arch epithelial ridge (GAER) of a cartilaginous fish (skate), this paper uses RNA–Sequencing, multiplexed HCR, cell fate mapping, and pharmacological manipulation to uncover the contribution of the Wnt signalling pathway to the patterning of the gill arch, highlighting the role tissues boundaries in signalling and patterning embryonic tissues.
https://doi.org/10.7554/eLife.79964.sa0Introduction
The pharyngeal arches of vertebrates are a series of paired columns of tissue that form on either side of the embryonic head. Pharyngeal arches form as iterative outpockets of foregut endoderm contact the overlying surface ectoderm, and this meeting of endoderm and ectoderm generates columns lined laterally by ectoderm, medially by endoderm, and containing a core of mesoderm and neural crest-derived mesenchyme (Graham and Smith, 2001). In fishes, endodermal outpockets fuse with the overlying surface ectoderm, giving rise to the gill slits and the respiratory surfaces of the gills (Gillis and Tidswell, 2017). In amniotes, endodermal outpockets give rise to glandular tissues, such as the tonsils, parathyroid, and ultimobranchial glands (Grevellec and Tucker, 2010). The largely neural crest-derived mesenchyme of the pharyngeal arches gives rise to the pharyngeal skeleton – i.e., the skeleton of the jaws and gills in fishes, and of the jaw, auditory ossicles, and larynx in amniotes (Jiang et al., 2002; Kague et al., 2012; Sleight and Gillis, 2020; Couly and Le Douarin, 1990) – with mesenchymal derivatives of each arch receiving patterning and polarity information via signals from adjacent epithelia (Veitch et al., 1999; Gillis et al., 2009a, Couly et al., 2002; Brito et al., 2006). Pharyngeal arches are named according to their ancestral skeletal derivatives: the 1st pharyngeal arch is termed the mandibular arch, the 2nd pharyngeal arch is the hyoid arch, while the caudal pharyngeal arches are collectively termed the gill arches. In fishes, the gill arches give rise to a series of skeletal ‘arches’ that support the lamellae of the gills.
All jawed vertebrates belong to one of two lineages: cartilaginous fishes (sharks, skates, rays, and holocephalans) or bony vertebrates (ray- and lobe-fined fishes, with the latter including tetrapods). While the gill arch skeleton of both cartilaginous and bony fishes ancestrally consisted proximally of two principal gill arch cartilages (the epi- and ceratobranchials), the gill arch skeleton of cartilaginous fishes additionally includes a distal series of fine cartilaginous rods called branchial rays (Figure 1A–C). These rays reflect the clear anteroposterior polarity of the gill arch skeleton of cartilaginous fishes, originating along the posterior margin of the epi- and ceratobranchial cartilage, and curving posteriorly as they project into the interbranchial septum of each arch (Gillis et al., 2009b). Elasmobranch cartilaginous fishes (sharks, skates, and rays) possess five sets of branchial rays, associated with their hyoid and first four gill arches. Holocephalans, on the other hand, possess a single set of branchial rays supporting their hyoid arch-derived operculum (Gillis et al., 2011).
The branchial rays of cartilaginous fishes develop under the influence of the GAER (gill arch epithelial ridge): a Sonic hedgehog (Shh)-expressing signaling centre located within the posterior–distal epithelium of the gill arches (Gillis et al., 2009b; Gillis and Hall, 2016). As the gill arches undergo a prolonged phase of lateral expansion, Shh signaling from the GAER is asymmetrically transduced within the posterior arch environment, as evidenced by posterior expression of Ptc2 (Gillis and Hall, 2016) – a Shh co-receptor and transcriptional readout of Shh signaling (Pearse et al., 2001). This posterior transduction of Shh signaling, in turn, appears to underlie the anteroposterior polarity of the gill arch skeleton: in skate, branchial rays derive exclusively from posterior–distal (Shh-responsive/GAER-adjacent) gill arch mesenchyme, application of exogenous Shh protein within gill arches is sufficient to induce ectopic branchial ray formation, and targeted or systemic inhibition of hedgehog signaling using cyclopamine results in branchial ray deletion (Gillis et al., 2009b; Gillis and Hall, 2016). However, whether other signaling mechanisms function alongside or in conjunction with GAER Shh signaling to establish and maintain gill arch skeletal polarity remains unexplored.
Here, we show in the little skate (L. erinacea) that the ectoderm immediately adjacent to the GAER expresses several genes encoding Wnt ligands. These Wnt signals are transduced in the anterior arch environment in a pattern that is broadly complementary to the posterior transduction of GAER Shh signaling. Inhibition of Wnt signaling in developing skate embryos results in an anterior expansion of Shh signaling transduction, and in the formation of ectopic branchial rays in the anterior gill arch. We propose that Wnt signaling from the pharyngeal arch ectoderm contributes to the maintenance of anterior–posterior polarity of the skate gill arch skeleton by repressing Shh signaling and chondrogenesis to the posterior gill arch territory.
Results
Wnt gene expression in the skate GAER and GAER-adjacent epithelium
In developing skate gill arches, Shh is expressed in the GAER (Figure 2A), and branchial rays develop from a GAER-responsive domain of posterior arch mesenchyme (Gillis and Hall, 2016). To further explore the transcriptional environment of the GAER and to discover additional gene expression features that may contribute to gill arch polarity, we conducted a comparative transcriptomic and differential gene expression analysis of GAER and non-GAER regions of the first gill arch of the embryonic skate. Briefly, we manually dissected the (1) GAER and (2) non-GAER (control) regions from the first gill arch of stage (S)26 skate embryos (n=5; Figure 2B–Bi), and we performed RNA extraction, library preparation, and RNAseq analysis on these samples. For the GAER region of the first gill arch, we included the whole distal tip of the arch in our dissection to ensure that the Shh-expressing GAER and adjacent Shh-responsive (i.e. Ptc2+) tissues were captured for this analysis (Figure 2C–D). Following de novo transcriptome assembly, we tested for differential gene expression between GAER and non-GAER tissues (after first testing for differential expression of Shh in the GAER region, to confirm dissection accuracy). We used a false discovery rate of 5% and no log fold-change (LogFC) cut-off to generate a list of 401 genes that were upregulated in the GAER region. We sorted this list for genes encoding signaling pathway components and transcription factors with a log-fold change >2 (Figure 2E). Differentially expressed genes were subjected to functional enrichment analysis, and this analysis indicated enrichment of Wnt signaling pathway components in the GAER region (Figure 2—figure supplement 1 for enrichment analysis). Differentially expressed genes in the GAER region included those encoding the Wnt ligands Wnt2b, Wnt3, Wnt4, Wnt7b, and Wnt9b, the transmembrane Wnt inhibitors Kremen1 (Mao et al., 2002) and APCDD1 (Shimomura et al., 2010), and the secreted Wnt antagonist Notum (Zhang et al., 2015). We, therefore, chose to further investigate the involvement of the Wnt pathway in GAER signaling by spatially validating the expression of these genes, using the most highly expressed isoforms for mRNA in situ hybridization by chain reaction (HCR) probe design.
To characterize the expression of genes encoding Wnt ligands relative to Shh expression in the GAER, we used multiplexed ISH by HCR to visualize transcript localization in sections of S26 skate gill arches. We found that while Wnt7b was co-expressed with Shh in the cells of the GAER (Figure 3A), Wnt2b, Wnt3, Wnt4, and Wnt9b (Figure 3B–E) were all expressed predominantly in the epithelium immediately adjacent to the GAER. We next examined the expression of the Wnt signaling downstream target genes (Figure 4) using Shh or Wnt7b as a marker for the GAER. We found that Apcdd1 was expressed in the GAER and anterior–distal mesenchyme (Figure 4A), Notum in the anterior–distal mesenchyme (Figure 4B), and Kremen1 in the anterior (GAER-adjacent) epithelium (Figure 4C). Finally, broad transcription of Axin2 throughout distal gill arch tissues (Figure 4D) indicates that Wnt signaling within this territory is occurring through the canonical/β-catenin pathway (Lustig et al., 2002). While these findings do not allow us to attribute expression of specific transcriptional readouts to signaling by particular Wnt ligands, our spatial expression data nevertheless indicate that some Wnt signals emanating from the GAER or GAER-adjacent epithelium are transduced preferentially within the anterior arch environment, in a pattern that is largely complimentary to the posterior epithelial and mesenchymal transduction of Shh signals from the GAER.
Wnt and Shh expression at the ectoderm–endoderm interface of skate gill arches
The epithelium of vertebrate pharyngeal arches derives from both ectoderm and endoderm, but the specific germ layer origin of the GAER of cartilaginous fishes (endoderm vs. ectoderm) is not known. In the skate, Shh is initially expressed broadly in the anterior endodermal epithelium of each pharyngeal pouch, which subsequently gives rise to the posterior endodermal epithelium of each pharyngeal arch (Gillis and Tidswell, 2017). To test whether the GAER derives from this initially broad endodermal Shh expression domain, we labeled the pharyngeal endoderm of early skate embryos by microinjecting the lipophilic dye CM-DiI into the pharyngeal cavity at stage S18. At this stage, the pharyngeal pouches have not yet fused with the overlying surface ectoderm, allowing for specific CM-DiI labeling of the pharyngeal endodermal epithelium (Figure 5A and B). Injected embryos were then grown to S25–S29, by which time the GAER of the hyoid and gill arches is detectable by expression of Shh (Figure 5C).
Of 29 embryos analyzed, all retained some CM-DiI label within the endodermal lining of the pharyngeal arches. In 8/29 embryos, we recovered CM-DiI labeling up to and including the Shh-expressing cells of the GAER in one or more pharyngeal arches (Figure 5C–C’’), indicating the endodermal origin of these cells. The remaining 21 embryos did not have sufficient labeling of pharyngeal epithelia to inform the origin of the GAER (likely due to label dilution during the long incubation time between injection and analysis). Importantly, however, no CM-DiI labeling was ever observed in epithelial cells immediately anterior to the GAER, suggesting that this epithelium is of ectodermal origin. Together, these observations indicate that the GAER is of endodermal origin and arises at the ectoderm–endoderm interface during pharyngeal arch development and that the Wnt-expressing epithelium immediately anterior to the GAER derives from the ectoderm.
Wnt signaling suppresses GAER Shh signal transduction in the anterior gill arch
To explore the function of Wnt signaling during skate pharyngeal arch development, we inhibited canonical Wnt signaling in skate embryos using IWR1, a small molecule tankyrase inhibitor that antagonizes Wnt signaling by stabilizing the Axin/β-catenin destruction complex (Lu et al., 2009; Kulak et al., 2015). Briefly, skate embryos were maintained in a bath of 50 μM IWR1 or vehicle-only control (DMSO) in seawater from S25 for 72 hr to assess gene expression changes in response to drug treatment (Figure 4—figure supplement 1), or for 6–8 weeks (until S31/32) to assess effects of Wnt signaling inhibition on gill arch skeletal patterning, with replacement of drug- or control-seawater every 48 hr. IWR1-treated embryos showed a marked reduction in the expression of Axin2 throughout their gill arches (Figure 4—figure supplement 1A–B), as well as a reduction in the expression of Apcdd1 in the GAER and anterior–distal gill arch mesenchyme (Figure 4—figure supplement 1C–D) and a complete loss of Notum expression in anterior–distal gill arch mesenchyme (Figure 4—figure supplement 1E–F). These findings indicate that our drug treatment effectively suppressed canonical Wnt signaling from the GAER and GAER-adjacent epithelium, including those signals transduced asymmetrically in the anterior gill arch environment.
Shh is a pro-chondrogenic signal from the GAER in skate gill arches (Gillis et al., 2009b), and this signal is transduced asymmetrically with the posterior gill arch environment (Gillis and Hall, 2016, see above Figure 2D). To test whether anteriorly transduced Wnt signaling may act to restrict Shh signal transduction to the posterior gill arch, we examined the expression of Shh and Ptc2 in skate embryos 72 hr after the onset of IWR1 treatment (Figure 6). While we observed no noticeable change in the expression of Shh or in the spatial distribution of Ptc2 transcripts in gill arch mesenchyme with IWR1 treatment, we did note a qualitative increase in Ptc2 expression within the normal mesenchymal domain and a clear anterior expansion of Ptc2 expression within the epithelium adjacent to the GAER (Figure 6A–D). To quantify this epithelial expansion, we compared the mean number of Ptc2+ nuclei anterior to the Shh-expressing GAER (using DAPI staining and counted manually) across two or three sections at equivalent positions in DMSO control (n=3) or IWR1 (n=4) treated embryos and found a significant anterior expansion of Ptc2 expression within the epithelium adjacent to the GAER in IWR1 treated embryos (p=0.036; Figure 6E–G). This points to a role for GAER-adjacent Wnt signaling in restricting the transduction of pro-chondrogenic Shh signaling in the anterior gill arch territory in skate.
Wnt signaling inhibits branchial ray development in the anterior gill arch environment
Finally, to test for an effect of Wnt inhibition on the gill arch skeleton, we reared skate embryos in DMSO control seawater or IWR1 from S25–S31/32 (Figure 7). IWR1-treated embryos had several distinct phenotypic changes, including smaller size, a loss of gill filaments and disrupted fin morphology (the latter consistent with previous findings of Nakamura et al., 2015; Figure 7A–B). We specifically investigated the impact of IWR1 treatment on gill arch skeletal morphology. The skate hyoid and gill arch skeleton include branchial rays that articulate exclusively with the posterior margin of the gill arch cartilages (Gillis et al., 2009a – see also Figure 1). These rays also invariably develop on the posterior side of the interbranchial muscle plate – a histologically and molecularly distinct derivative of the paraxial mesodermal core that extends down the middle of each pharyngeal arch (Daniel, 1934; Graham, 2008; Hirschberger et al., 2021). Skate embryos reared in DMSO control seawater showed no gill arch skeletal defects (n=4 embryos; Figure 7Ai), while embryos reared in IWR1 seawater possessed conspicuous ectopic branchial rays in the anterior gill arch territory (n=7/8 embryos; Figure 7Bi). These ectopic branchial rays occurred along the entire dorsoventral axis of the arch and ranged from short cartilage rods located in the anterior–distal gill arch to nearly full-length branchial rays that extended much of the way toward the gill arch cartilages. In section, branchial rays of skate embryos reared in DMSO control seawater always developed, as per normal, on the posterior side of the mesodermally-derived interbranchial muscle plate of each gill arch: this was evident from histochemical staining and from HCR for the chondrocyte markers Sox9 and Col2a1 (Figure 8A–B). Conversely, in skate embryos reared in IWR1 seawater, ectopic branchial rays were consistently located anterior to the interbranchial muscle plate (Figure 8C–D). The induction of ectopic branchial rays with IWR1 treatment suggests that Wnt signaling contributes to the maintenance of anteroposterior polarity of the skate gill arch skeleton by repressing chondrogenesis in the anterior gill arch and may do so by restricting transduction of GAER Shh signaling to the posterior gill arch territory.
Discussion
The molecular mechanisms governing spatial regulation of cell fate decisions are crucial to the establishment and maintenance of anatomical polarity within developing tissues and organs. Here, we show that the GAER signaling centre that forms on the hyoid and gill arches of cartilaginous fishes is of endodermal origin, arising precisely at the boundary of surface ectoderm and pharyngeal endoderm, and that pro-chondrogenic Shh signals from the GAER are transduced preferentially within the posterior gill arch environment. We further show that the ectoderm immediately anterior and adjacent to the GAER is a signaling centre that expresses several Wnt family genes and that Wnt signals from this centre are transduced preferentially in the anterior arch environment. Pharmacological inhibition of canonical Wnt signaling in skate embryos leads to an expansion of Shh signal transduction within the anterior territory of the gill arch, and to the formation of ectopic anterior branchial ray cartilage. These findings suggest that Wnt–Shh signaling antagonism at the pharyngeal ectodermal–endodermal tissue boundary is a key molecular regulator of cell fate determination and anatomical polarity with the developing gill arch skeleton of the skate (Figure 9).
Previous analyses have reported the expression of several Wnt genes in developing craniofacial tissues of bony vertebrates (Summerhurst et al., 2008), including the expression of WNT2B and WNT9B in the pharyngeal arch ectoderm of the chick (Geetha-Loganathan et al., 2009), Wnt9b in the first and second pharyngeal arches in zebrafish (Jezewski et al., 2008), Wnt4a in zebrafish pharyngeal ectoderm (Choe et al., 2013), and Wnt9a and Wnt5b in zebrafish oral epithelium (Dougherty et al., 2013; Rochard et al., 2016). While β-catenin-independent (i.e. non-canonical) Wnt signaling plays a crucial role in epithelial remodeling during the early formation of pharyngeal endodermal pouches in zebrafish (Choe et al., 2013) and in maturation and morphogenesis of the zebrafish oropharyngeal skeleton (Rochard et al., 2016; Ling et al., 2017), canonical Wnt signaling has been implicated in patterning of the pharyngeal arch-derived skeleton in multiple taxa. Blocking canonical Wnt signaling causes disrupted facial cartilage formation, including reductions of Meckel’s and ceratohyal cartilages in zebrafish (Alexander et al., 2014) and loss of ectodermal canonical Wnt signaling in mouse causes loss of facial bones with hypoplasia of all facial prominences (Reid et al., 2011). Furthermore, in mouse and human, cleft palate phenotypes may arise as a consequence of mutations in Wnt3 and Wnt9b (Niemann et al., 2004; Juriloff et al., 2006), and perturbation of Wnt/β-catenin signaling has been implicated in CATSHL syndrome, a human developmental disorder involving craniofacial bone malformation and mispatterning of the pharyngeal arches (Sun et al., 2020). Our finding of a role for canonical Wnt signaling in patterning the gill arch skeleton of skate may, therefore, reflect broadly conserved roles for this pathway in pharyngeal skeletal development across jawed vertebrates, though with the skate’s possession of branchial rays offering a unique anatomical readout of anterior–posterior polarity defects that arise within the pharyngeal arch-derived skeleton in response to perturbations.
Our findings are also consistent with previously reported complex and context-dependent roles for Wnt signaling in the regulation of cell fate determination events during vertebrate skeletogenesis. In some instances, Wnt signaling functions to promote skeletogenesis: for example, in chick, misexpression of WNT5A/5B promotes early chondrogenesis of limb bud mesenchyme in vitro and delays the terminal differentiation of growth plate chondrocytes in vivo (Church et al., 2002), and in mouse, canonical WNT signaling from cranial ectoderm induces specification of osteoblast progenitors of cranial dermal bone within underlying mesenchyme (Goodnough et al., 2016). However, in other contexts, Wnt signaling functions to inhibit chondrogenic differentiation or homeostasis. In chick, signaling through WNT3, 4, 7 A, 14 or FZ7 have all been shown to inhibit chondrogenesis in vitro or in vivo (Rudnicki and Brown, 1997), while in mouse conditionally induced haploinsufficiency of the gene encoding the β-catenin degradation complex component APC (i.e. activation of canonical Wnt signaling) results in loss of resting zone chondrocytes and their clonal progeny in developing growth plates (Hallett et al., 2021). In the limb bud, genes encoding WNT ligands are expressed in the ectoderm (Kengaku et al., 1998; Parr et al., 1993; Geetha-Loganathan et al., 2005; Summerhurst et al., 2008), where they function synergistically with fibroblast growth factor signaling from the apical ectodermal ridge to inhibit chondrogenesis of limb bud mesenchyme and to promote soft connective tissue fates (ten Berge et al., 2008). Downstream effectors of Wnt signaling are expressed out-of-phase with Sox9+ digit progenitors in the mouse limb bud, with Wnt signaling functioning to repress digit chondrogenesis as part of a Turing-like mechanism that underlies periodic patterning of the distal limb endoskeleton (Raspopovic et al., 2014). Our findings, therefore, demonstrate a striking parallel between developing skate gill arches and tetrapod limb buds, with Wnt signals emanating from ridge-adjacent epithelia functioning to inhibit chondrogenic differentiation, thereby ensuring differentiation of cartilaginous appendages in a spatially controlled manner.
We have previously shown that Shh signaling from the GAER is required for branchial ray chondrogenesis in skate (Gillis and Hall, 2016), and that the application of exogenous Shh protein to the anterior gill arch is sufficient to induce ectopic branchial rays (Gillis et al., 2009b). However, it is unclear whether Shh from the GAER is a direct inducer of chondrogenesis in the posterior arch mesenchyme, or rather induces chondrogenesis indirectly, via a secondary signal (i.e. emanating from the Shh-responsive epithelium or mesenchyme). We have observed that inhibition of Wnt signaling in skate embryos has no effect on the expression of Shh in the GAER, but rather results in an anterior expansion of Shh signal transduction in the gill arch ectoderm, which in turn, correlates with ectopic chondrogenesis in the anterior gill arch. This observation indicates that the pro-chondrogenic influence of GAER Shh signaling may be indirect, occurring via a secondary Shh-depending epithelial signal whose expression is typically posteriorly restricted by Wnt signaling from GAER-adjacent ectoderm. The Shh and Wnt signaling pathways share multiple downstream components (reviewed in Ding and Wang, 2017). For example, components of the canonical Wnt signaling pathway can positively regulate Shh signaling via GSK3β, by phosphorylating the downstream component SUFU and promoting the release of SUFU from Gli (Takenaka et al., 2007), and by β-catenin affecting Gli1 transcriptional activity via TCF/LEF (Maeda et al., 2006). It has also been found that Gli3 may function as a downstream effector of the Wnt pathway and that Wnt signaling represses Shh activity in the dorsal neural tube through the regulation of Gli3 expression (Alvarez-Medina et al., 2008). Finally, the Shh signaling pathway may regulate Wnt signaling through Gli1 and Gli2, with these factors positively regulating the expression of the secreted Wnt inhibitor frizzled-related protein-1 (sFRP-1) (He et al., 2006). These mechanisms, or others, could account for the apparent cross-regulation between Wnt and Shh signaling that we have observed during the growth and patterning of skate pharyngeal arches.
Finally, it remains to be determined whether antagonistic ectodermal Wnt and endodermal Shh signaling is a derived mechanism of cell fate determination and pharyngeal skeletal patterning in cartilaginous fishes or a conserved feature of the pharyngeal arches of jawed vertebrates. In bony vertebrates, a Shh-expressing posterior endodermal margin (PEM) promotes the proliferative expansion of the second (hyoid) pharyngeal arch (Wall and Hogan, 1995; Richardson et al., 2012; Rees and Gillis, 2022). In bony fishes, the hyoid arch gives rise to an operculum that is supported by a series of dermal bones, and this operculum functions as a protective cover for the gills that arise from the posterior pharyngeal arches, whereas in amniotes, the hyoid arch expands caudally and ultimately fuses with the cardiac eminence to enclose the cervical sinus and close the surface of the neck (Richardson et al., 2012). Amniotes lack homologs of the opercular series of bony fishes or the branchial rays of cartilaginous fishes, respectively, but fossil evidence points to the operculate condition of bony fishes as ancestral for jawed vertebrates (Dearden et al., 2019). Whether juxtaposed Wnt- and Shh-expressing epithelia influence fate determination (e.g. dermal bone vs. connective tissue) and patterning within the opercular series of bony fishes is not known. Further comparative work between bony and cartilaginous fishes with disparate pharyngeal skeletal architectures and tissue compositions will permit the inference of ancestral functions and interactions of epithelial signals during vertebrate pharyngeal arch skeletal development.
Materials and methods
Embryo harvesting
Request a detailed protocolSkate (L. erinacea) embryos were obtained from the Marine Resources Center at the Marine Biological Laboratory in Woods Hole, MA, U.S.A., were reared as described in Gillis et al., 2012 and staged according to Ballard et al., 1993 and Maxwell et al., 2008. Skate embryos were euthanized with an overdose of MS-222 (1 g/L in seawater), and all embryos were fixed in 4% paraformaldehyde (Thermofisher) as per Gillis et al., 2012.
Sectioning and histochemical staining
Request a detailed protocolEmbryos were embedded in paraffin and sectioned at 7 μm as described in O’Neill et al., 2007. Sectioned embryos were stained with Masson’s Trichrome as per Witten and Hall, 2003 or Mayer’s Hematoxylin and Eosin by clearing in two rinses of Histosol (National Diagnostics), rehydration through serial dilutions of ethanol, staining for 15 min in Mayer’s Haematoxylin (Sigma), rinsing in running tap water for 20 min, rinsing briefly in 95% ethanol, and staining in 0.1% w/v Eosin Y (Sigma) in 95% ethanol for 2 min. Slides were then washed briefly in 95% ethanol, washed briefly with 100% ethanol, cleared with Histosol, and mounted with permount (Sigma).
Lineage tracing
Request a detailed protocolEndodermal lineage tracing was performed by microinjection of CellTracker CM-DiI (1,1'-dioctadecyl-3,3,3'3'-tetramethylindocarbocyanine perchlorate; Thermofisher) into the pharyngeal cavity of the stage (S)18 skate embryos with a pulled glass needle. CM-DiI was prepared as per Gillis et al., 2012. Labeled embryos were then grown to S25–29, euthanized, fixed, embedded, and sectioned as described above, and CM-DiI localization was assessed relative to the GAER based on mRNA in situ hybridization for Shh (see below).
mRNA in situ hybridization
Request a detailed protocolChromogenic mRNA in situ hybridization (ISH) was performed on paraffin sections and in whole mount as described in O’Neill et al., 2007, with modifications according to Gillis et al., 2012. An ISH probe against skate Shh (GenBank EF100667) was generated by in vitro transcription using standard methods. Following ISH, whole mounts were rinsed, post-fixed in 4% paraformaldehyde, and graded into 75% glycerol for imaging, while slides were rinsed in PBS, post-fixed in 4% paraformaldehyde, and coverslipped with DAPI-Fluoromount G (SouthernBiotech). Third-generation mRNA ISH by chain reaction (HCR) was performed as per the Choi et al., 2018 protocol for formaldehyde-fixed, paraffin-embedded sections, with modifications as per Criswell and Gillis, 2020. Probes, buffers, and hairpins were purchased from Molecular Instruments (Los Angeles, California, USA). HCR probe set lot numbers are as follows: skate Shh (Lot PRA753), Ptc2 (Lot PRA754), Wnt2b (Lot PRE300), Wnt3 (Lot PRG814), Wnt4 (Lot PRE301), Wnt7b (Lot PRE302), Wnt9b (Lot PRE303), Notum (Lot PRG817), Kremen1 (Lot PRG816), Axin2 (Lot PRG818), Apcdd1 (Lot PRG815), Col2a1 (Lot PRB574), and Sox9 (Lot PRB571). All mRNA ISH experiments were replicated, at minimum, in biological triplicate.
Pharmacological manipulations
Request a detailed protocolFor systemic inhibition of canonical Wnt signaling, experimental and control skate embryos were reared at 15 °C in Petri dishes containing 50 μM IWR1 (Selleck chemicals) in artificial seawater or artificial seawater containing an equivalent volume of vehicle (DMSO) only, respectively. IWR1 has been used previously to inhibit canonical Wnt signaling in skate embryos by Nakamura et al., 2015 and in shark embryos by Thiery et al., 2022. For these experiments, IWR1 was diluted to a working concentration from a 25 mM stock solution in DMSO. Skate embryos were reared in IWR1- or vehicle-containing seawater from S25 until S31/32, with drug or control seawater changes every 48 hr. Once embryos reached the desired stage, they were euthanized, fixed, and processed for whole mount skeletal preparation according to the protocol of Gillis et al., 2009a or paraffin histology as described above.
RNA-Seq analysis of differentially expressed genes in the GAER region
Request a detailed protocolGAER or non-GAER (‘control’) regions of gill arch 1 were manually dissected from skate embryos at S26 using tungsten needles and flash frozen in lysis buffer using liquid nitrogen. In all cases, GAER and control sample pairs were collected from the same arch within the same embryos, or from opposite arches within the same embryo. RNA was extracted from each sample using the RNAaqueous-Micro Kit. cDNA was synthesized from extracted RNA according to the Smart-seq2 protocol (Picelli et al., 2014) and libraries were prepared using the Nextera XT kit (Illumina). Prior to sequencing, barcoded libraries were pooled in equal proportions and checked for quality, insert size, and quantity using Qubit 2.0, and Agilent2100. Sequencing of the 20 libraries generated was conducted by Novogene on an Illumina Hi-Seq-XTen generating 150 bp paired-end reads. The Cambridge Service for Data-Driven Discovery (CSD3) high performance computer was used for cleaning, normalization, quality assessment, assembly, abundance quantification, and annotation, as per Hirschberger et al., 2021. Raw reads were cleaned to remove adapter contamination using Trim Galore v0.4.0 (parameters: --paired -q 20 –nextera), for quality (Phred score 20) and minimum read length (149 bp) using ea-utils tool fastq-mcf (parameters: -q 20 l 100) (Aronesty, 2011). Prior to assembly, all libraries were normalized using Trinity v2.6.6 in silico read normalization (Haas et al., 2013) (parameters: --JM 100 G --max_cov 75). Trinity v2.6.6 was used to assemble the normalized reads (parameters: --seqType fq --max_memory 250 G --CPU 50 min_contig_length 450) (Nishimura et al., 2017; Langmead et al., 2009; Marçais and Kingsford, 2011). Basic Local Alignment Search Tool (BLAST) was used to compare the assembly to the Uniprot protein database with an e-value cut-off of e-20 to search for sequence similarity (Altschul et al., 1990). Transcript abundance was estimated using Salmon, an alignment-free quantification method, via Trinity v2.6.6 utilities (Haas et al., 2013; Patro et al., 2017). A matrix was built for each of the libraries of gene-level abundance estimates (raw counts) using the Trinity abundance_estimates_to_matrix.pl script. MDS principal component analysis was used to check for batch effects and outliers. Individual embryo origin was a batch effect and so was accordingly accounted for in the additive general linear model. Differentially expressed genes between control and GAER tissue samples were identified using the edgeR_3.28.0 package in R version 3.6.2 (Robinson et al., 2010) using a generalized linear model likelihood ratio test (Chen, 2008) using a false discovery rate of 5% and no log fold-change (LogFC) cut-off. See supplementary information for assembly statistics, differential expression analysis, and all R scripts used.
Significantly upregulated genes in the S26 GAER tissue (logFC >2, FDR <0.05) were extracted and annotated using blastx sequence similarity searches against the Swissprot/Uniprot database (taxid: Gallus gallus). Annotation IDs were inputted into the STRING (v11.5) webserver and tested for predicted protein-protein interactions. In the protein-protein interaction analysis, line thickness was set to show confidence, and a medium confidence threshold was set. The interaction network was then tested for functional enrichment against a series of databases (Gene Ontology Biological Process, KEGG Pathways, WikiPathways, Pfam Domains, and InterPro Domains) using the default whole genome as the background. Results from the functional enrichment were downloaded and Rstudio (GGplot2) was used to generate dot plots to visualize the most significantly enriched terms.
Differentially upregulated genes with a log-fold change over two were manually screened for genes involved in known signaling pathways and transcription factors to generate a list for candidate gene validation by ISH. Wnt family members were highly upregulated in the GAER region so were validated as candidate markers for this region through third-generation ISH. Transcripts for probes in skate were based on the most highly expressed isoforms.
Imaging and figure preparation
Request a detailed protocolImages were taken on a Zeiss Axioscope.A1 compound microscope with a Zeiss colibri 7 fluorescence LED light source using a Zeiss Axiocam 305 color or 503 mono camera and ZenPro software or a Leica M165FC stereomicroscope, a Leica DFC7000 T camera and LAS X software. All figures were assembled using Fiji and Adobe creative cloud, with some images flipped or colors inverted for clarity and consistency where needed.
Data availability
All data and R scripts for analysis have been deposited on Figshare: 10.6084/m9.figshare.19615779. RNA sequencing data are available at NCBI-SRA under BioProject ID: PRJNA825354. Biosample accessions: SAMN27512544, SAMN27512545, SAMN27512546, SAMN27512547, SAMN27512548, SAMN27512549, SAMN27512550, SAMN27512551, SAMN27512552, SAMN27512553, SAMN27512554, SAMN27512555, SAMN27512556, SAMN27512557, SAMN27512558, SAMN27512559, SAMN27512560, SAMN27512561, SAMN27512562, SAMN27512563. Probes were purchased from Molecular Instruments (Los Angeles, California, USA). This included the following: for skate Shh (Lot PRA753), Ptc2 (Lot PRA754), Wnt2b (Lot PRE300), Wnt3 (Lot PRG814), Wnt4 (Lot PRE301), Wnt7b (Lot PRE302), Wnt9b (Lot PRE303), Notum (Lot PRG817), Kremen1 (Lot PRG816), Axin2 (Lot PRG818), Apcdd1 (Lot PRG815), Col2a1 (Lot PRB574) and Sox9 (Lot PRB571).
-
NCBI BioProjectID PRJNA825354. Ectodermal Wnt signaling, cell fate determination and polarity of the skate gill arch skeleton.
-
figshareScripts and data on Wnt signalling, cell fate determination and anteroposterior polarity of the skate gill arch skeleton.https://doi.org/10.6084/m9.figshare.19615779.v2
References
-
Basic local alignment search toolJournal of Molecular Biology 215:403–410.https://doi.org/10.1016/S0022-2836(05)80360-2
-
A series of normal stages for development ofscyliorhinus canicula, the lesser spotted dogfish(chondrichthyes: scyliorhinidae)Journal of Experimental Zoology 267:318–336.https://doi.org/10.1002/jez.1402670309
-
Wnt-dependent epithelial transitions drive pharyngeal pouch formationDevelopmental Cell 24:296–309.https://doi.org/10.1016/j.devcel.2012.12.003
-
Wnt regulation of chondrocyte differentiationJournal of Cell Science 115:4809–4818.https://doi.org/10.1242/jcs.00152
-
BookThe Elasmobranch FishesBerkeley,Calif: University of California press.https://doi.org/10.5962/bhl.title.10180
-
Antagonism between hedgehog and Wnt signaling pathways regulates tumorigenicityOncology Letters 14:6327–6333.https://doi.org/10.3892/ol.2017.7030
-
Ectodermal wnt-6 promotes myf5-dependent avian limb myogenesisDevelopmental Biology 288:221–233.https://doi.org/10.1016/j.ydbio.2005.09.035
-
Expression of Wnt signalling pathway genes during chicken craniofacial developmentDevelopmental Dynamics 238:1150–1165.https://doi.org/10.1002/dvdy.21934
-
The origin of vertebrate gillsCurrent Biology 27:729–732.https://doi.org/10.1016/j.cub.2017.01.022
-
Deconstructing the pharyngeal metamereJournal of Experimental Zoology. Part B, Molecular and Developmental Evolution 310:336–344.https://doi.org/10.1002/jez.b.21182
-
The pharyngeal pouches and clefts: development, evolution, structure and derivativesSeminars in Cell & Developmental Biology 21:325–332.https://doi.org/10.1016/j.semcdb.2010.01.022
-
Suppressing Wnt signaling by the hedgehog pathway through sFRP-1The Journal of Biological Chemistry 281:35598–35602.https://doi.org/10.1074/jbc.C600200200
-
Conserved and unique transcriptional features of pharyngeal arches in the skate (leucoraja erinacea) and evolution of the jawMolecular Biology and Evolution 38:4187–4204.https://doi.org/10.1093/molbev/msab123
-
Tissue origins and interactions in the mammalian skull vaultDevelopmental Biology 241:106–116.https://doi.org/10.1006/dbio.2001.0487
-
Wnt9B is the mutated gene involved in multifactorial nonsyndromic cleft lip with or without cleft palate in A/wysn mice, as confirmed by a genetic complementation testBirth Defects Research. Part A, Clinical and Molecular Teratology 76:574–579.https://doi.org/10.1002/bdra.20302
-
Disruption of Wnt/β-catenin signaling and telomeric shortening are inextricable consequences of tankyrase inhibition in human cellsMolecular and Cellular Biology 35:2425–2435.https://doi.org/10.1128/MCB.00392-15
-
Structure-Activity relationship studies of small-molecule inhibitors of Wnt responseBioorganic & Medicinal Chemistry Letters 19:3825–3827.https://doi.org/10.1016/j.bmcl.2009.04.040
-
Negative feedback loop of Wnt signaling through upregulation of conductin/axin2 in colorectal and liver tumorsMolecular and Cellular Biology 22:1184–1193.https://doi.org/10.1128/MCB.22.4.1184-1193.2002
-
Enhancement of GLI1-transcriptional activity by beta-catenin in human cancer cellsOncology Reports 16:91–96.
-
Homozygous WNT3 mutation causes tetra-amelia in a large consanguineous familyAmerican Journal of Human Genetics 74:558–563.https://doi.org/10.1086/382196
-
Full-Length RNA-seq from single cells using smart-seq2Nature Protocols 9:171–181.https://doi.org/10.1038/nprot.2014.006
-
Ectodermal Wnt/β-catenin signaling shapes the mouse faceDevelopmental Biology 349:261–269.https://doi.org/10.1016/j.ydbio.2010.11.012
-
The presence of an embryonic opercular flap in amniotesProceedings. Biological Sciences 279:224–229.https://doi.org/10.1098/rspb.2011.0740
-
Inhibition of chondrogenesis by Wnt gene expression in vivo and in vitroDevelopmental Biology 185:104–118.https://doi.org/10.1006/dbio.1997.8536
-
Gsk3Beta positively regulates hedgehog signaling through Sufu in mammalian cellsBiochemical and Biophysical Research Communications 353:501–508.https://doi.org/10.1016/j.bbrc.2006.12.058
-
Pharyngeal arch patterning in the absence of neural crestCurrent Biology 9:1481–1484.https://doi.org/10.1016/s0960-9822(00)80118-9
Article and author information
Author details
Funding
Wellcome Trust (214953/Z/18/Z)
- Jenaid M Rees
Royal Society (UF130182)
- J Andrew Gillis
Royal Society (URF\R\191007)
- J Andrew Gillis
Royal Society (RG140377)
- J Andrew Gillis
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication. For the purpose of Open Access, the authors have applied a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission.
Acknowledgements
With thanks to Dr Kate Criswell and Dr Christine Hirschberger for advice, and to the University of Cambridge Wellcome PhD. Programme in Developmental Mechanisms. The authors were funded by a Wellcome PhD studentship (214953/Z/18/Z) to JMR, and by a Royal Society University Research Fellowship (UF130182 and URF\R\191007) and Royal Society Research Grant (RG140377) to JAG.
Copyright
© 2023, Rees 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
-
- 867
- views
-
- 96
- downloads
-
- 3
- 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
-
- Cell Biology
- Developmental Biology
How cells regulate the size of their organelles remains a fundamental question in cell biology. Cilia, with their simple structure and surface localization, provide an ideal model for investigating organelle size control. However, most studies on cilia length regulation are primarily performed on several single-celled organisms. In contrast, the mechanism of length regulation in cilia across diverse cell types within multicellular organisms remains a mystery. Similar to humans, zebrafish contain diverse types of cilia with variable lengths. Taking advantage of the transparency of zebrafish embryos, we conducted a comprehensive investigation into intraflagellar transport (IFT), an essential process for ciliogenesis. By generating a transgenic line carrying Ift88-GFP transgene, we observed IFT in multiple types of cilia with varying lengths. Remarkably, cilia exhibited variable IFT speeds in different cell types, with longer cilia exhibiting faster IFT speeds. This increased IFT speed in longer cilia is likely not due to changes in common factors that regulate IFT, such as motor selection, BBSome proteins, or tubulin modification. Interestingly, longer cilia in the ear cristae tend to form larger IFT compared to shorter spinal cord cilia. Reducing the size of IFT particles by knocking down Ift88 slowed IFT speed and resulted in the formation of shorter cilia. Our study proposes an intriguing model of cilia length regulation via controlling IFT speed through the modulation of the size of the IFT complex. This discovery may provide further insights into our understanding of how organelle size is regulated in higher vertebrates.
-
- Developmental Biology
- Genetics and Genomics
Repurposing of pleiotropic factors during execution of diverse cellular processes has emerged as a regulatory paradigm. Embryonic development in metazoans is controlled by maternal factors deposited in the egg during oogenesis. Here, we explore maternal role(s) of Caspar (Casp), the Drosophila orthologue of human Fas-associated factor-1 (FAF1) originally implicated in host-defense as a negative regulator of NF-κB signaling. Maternal loss of either Casp or it’s protein partner, transitional endoplasmic reticulum 94 (TER94) leads to partial embryonic lethality correlated with aberrant centrosome behavior, cytoskeletal abnormalities, and defective gastrulation. Although ubiquitously distributed, both proteins are enriched in the primordial germ cells (PGCs), and in keeping with the centrosome problems, mutant embryos display a significant reduction in the PGC count. Moreover, the total number of pole buds is directly proportional to the level of Casp. Consistently, it’s ‘loss’ and ‘gain’ results in respective reduction and increase in the Oskar protein levels, the master determinant of PGC fate. To elucidate this regulatory loop, we analyzed several known components of mid-blastula transition and identify the translational repressor Smaug, a zygotic regulator of germ cell specification, as a potential critical target. We present a detailed structure-function analysis of Casp aimed at understanding its novel involvement during PGC development.