Dermomyotome-derived endothelial cells migrate to the dorsal aorta to support hematopoietic stem cell emergence
Abstract
Development of the dorsal aorta is a key step in the establishment of the adult blood-forming system, since hematopoietic stem and progenitor cells (HSPCs) arise from ventral aortic endothelium in all vertebrate animals studied. Work in zebrafish has demonstrated that arterial and venous endothelial precursors arise from distinct subsets of lateral plate mesoderm. Here, we profile the transcriptome of the earliest detectable endothelial cells (ECs) during zebrafish embryogenesis to demonstrate that tissue-specific EC programs initiate much earlier than previously appreciated, by the end of gastrulation. Classic studies in the chick embryo showed that paraxial mesoderm generates a subset of somite-derived endothelial cells (SDECs) that incorporate into the dorsal aorta to replace HSPCs as they exit the aorta and enter circulation. We describe a conserved program in the zebrafish, where a rare population of endothelial precursors delaminates from the dermomyotome to incorporate exclusively into the developing dorsal aorta. Although SDECs lack hematopoietic potential, they act as a local niche to support the emergence of HSPCs from neighboring hemogenic endothelium. Thus, at least three subsets of ECs contribute to the developing dorsal aorta: vascular ECs, hemogenic ECs, and SDECs. Taken together, our findings indicate that the distinct spatial origins of endothelial precursors dictate different cellular potentials within the developing dorsal aorta.
Editor's evaluation
In this important study, the authors identify an additional source of dorsal aorta endothelium derived from the somites that are conventionally thought of as defined blocks of skeletal muscle. The authors present convincing data that these "somite-derived endothelial cells" (SDECs) arise from bipotential precursors in the dermamyotome that give rise to endothelium or muscle. The authors conclude that these cells support but do not produce emergent hematopoietic stem cells in the dorsal aorta, findings that will be of interest to stem cell and developmental biologists.
https://doi.org/10.7554/eLife.58300.sa0Introduction
The primitive vascular network, which integrates into all organ systems in the developing organism, arises from endothelial precursors termed angioblasts (Risau and Flamme, 1995). To form a functional vascular network, angioblasts must first differentiate into a variety of distinct arterial and venous endothelial cell (EC) types, including hemogenic, endocardial, and blood brain barrier ECs (Aird, 2007; Herbert and Stainier, 2011). EC differentiation is thought to initiate midway through somitogenesis, during the migration of angioblasts to the embryonic midline, where they coalesce to form the vascular tube (Herbert et al., 2009; Isogai et al., 2003; Jin et al., 2005). The current view is that angioblasts are initially equipotent and undergo successive steps of differentiation that, along with cues from local microenvironments, give rise to specialized subsets of ECs (Atkins et al., 2011; Marcelo et al., 2013). This view, however, does not consider possible differences in EC function due to different developmental origins. Instead, exposure to embryonic signaling cascades mediated via Wnt (Hübner et al., 2017), Hedgehog (Hh) (Gering and Patient, 2005; Vokes and McMahon, 2004; Wilkinson et al., 2012; Williams et al., 2010), Vascular Endothelial Growth Factor (VEGF) (Casie Chetty et al., 2017; Hong et al., 2006; Lawson et al., 2003; Lawson et al., 2002; Wythe et al., 2013), and Notch molecules are thought to differentially instruct equipotent angioblasts to each distinct endothelial cell fate (Fang et al., 2017; Lawson et al., 2001; Siekmann and Lawson, 2007; Zhong et al., 2001).
Lateral plate mesoderm (LPM) is known to be the primary source of ECs across vertebrate phyla (Potente and Mäkinen, 2017). However, recent findings suggest that ECs can arise from distinct mesodermal derivatives, including extraembryonic-derived erythromyeloid progenitors (EMPs) that contribute extensively to the murine kidney vasculature (Plein et al., 2018). Furthermore, lineage tracing studies in zebrafish have demonstrated that ECs supportive of hematopoiesis can derive from endoderm (Nakajima et al., 2023). To better understand the development of EC subsets, we performed single-cell RNA sequencing (scRNA-seq) of ECs purified by flow cytometry over a range of time points during zebrafish embryogenesis. Following the end of gastrulation, the earliest developmental timepoint that nascent ECs can be identified, a variety of distinct molecular signatures were present, including those similar to kidney-specific ECs, brain-specific ECs, and paraxial mesoderm-derived ECs. These data suggest that specific EC fates may be specified much earlier than previously appreciated and that embryonic origins may dictate the development of tissue-specific EC subsets.
Within these datasets, we detected a signature indicative of paraxial mesoderm (PM) origins; we thus further explored this possibility since previous findings suggested that somite-derived ECs (SDECs) are intimately involved with both the development of the dorsal aorta and the emergence of hematopoietic stem and progenitor cells (HSPCs; Nguyen et al., 2014; Pouget et al., 2006). There is evidence in rodents and other amniotes that PM generates a contingent of ECs that contribute to embryonic vasculature (Ambler et al., 2001; Esner et al., 2006; Noden, 1989; Pardanaud et al., 1996; Pouget et al., 2008; Pouget et al., 2006; Wilting et al., 1995; Yvernogeau et al., 2012). Specifically, these ECs arise from a transient somitic compartment known as the dermomyotome (Eichmann et al., 1993; Ema et al., 2006; Pouget et al., 2008; Yvernogeau et al., 2012), where skeletal hypaxial muscle progenitors (skMPs) reside (Tozer et al., 2007). In the chick embryo, SDECs migrate to the dorsal aorta to replace hemogenic endothelium that exits the ventral aorta as emerging HSPCs (Pouget et al., 2006). A recent report in zebrafish has also suggested the existence of an endothelial-producing compartment (termed ‘endotome’) within the central somite that generates SDECs that migrate to the dorsal aorta (Nguyen et al., 2014). Whereas experiments in birds have suggested that SDECs only replace hemogenic endothelium, the experiments in zebrafish by Nguyen et al., 2014 suggested that SDECs may help induce the HSPC program. However, the genetic profile of SDECs and their possible role in HSPC induction remain poorly understood.
Here, using a combination of molecular, genetic, and computational approaches, we characterize a rare population of SDECs in zebrafish that emerges from the trunk dermomyotome. Trunk SDECs migrate and contribute exclusively to rostral regions in the dorsal aorta. Within the somite, EC-fate acquisition occurs sequentially, concomitant with the epithelialization of each somite and the migration of angioblasts toward the midline of the embryo. We show that Wnt signaling is a key regulator of the distribution of ECs within the somite, whereas Notch signaling is necessary for skeletal muscle progenitor cell maintenance. Finally, epistasis experiments indicate that SDECs arise from bipotent precursors within the somite, with skMPs showing competency to become ECs in a meox1- and npas4l (cloche Stainier et al., 1995) – dependent manner. While SDECs integrate into the dorsal aorta, indelible lineage tracing demonstrates that they do not harbor hematopoietic potential. Instead, they appear to act at least in part as a developmental niche to facilitate HSPC emergence. Collectively, our findings indicate that distinct EC subsets are molecularly and functionally distinct as early as the end of gastrulation and describe a cellular mechanism by which the somite regulates the production of hematopoietic precursors.
Results
Molecular differences in endothelial cells underlie cellular diversity within the vasculature
In zebrafish, most ECs originate from the LPM (Jin et al., 2005). However, recent studies have suggested that somites also produce ECs that integrate into the vascular cord (Nguyen et al., 2014), but the role and nature of these ECs remain incompletely defined. To better understand how different ECs are specified in the zebrafish embryo, we utilized single-cell RNA sequencing (scRNA-seq) on a collection of purified ECs from distinct vascular transgenic animals to represent diverse endothelial cohorts. Specifically, we purified cells via fluorescence-activated cell sorting (FACS) from TgBAC(etv2:Kaede)ci6, Tg(fli1:DsRed)um13; Tg(tp1:GFP)um14, and Tg(drl:H2B-dendra) (hereafter termed mixed-vasculature) (Kohli et al., 2013; Mosimann et al., 2015; Parsons et al., 2009; Villefranc et al., 2007) between 22 and 24 hpf, and performed scRNA-seq on each sample. Following quality control (QC), we merged the individual datasets. Altogether, we obtained 1994 single cells that belonged to the endothelial lineage and an additional 5446 erythroid and myeloid cells, which we omitted from our analyses as our focus was on the EC subsets. Following unsupervised clustering of single-cell transcriptomes, we identified erythroid cells, myeloid cells, and eight distinct endothelial cell clusters (Figure 1A and B), which we named based on known marker genes within established tissue lineages (Figure 1—source data 1 and Figure 1—figure supplement 1). First, we identified a sizeable endothelial cluster, general endothelium (GE), that co-expressed canonical endothelial genes and a putative hemogenic endothelium (HE) cluster that co-expressed endothelial and hematopoietic genes, including etv2, dab2, and stab2. We identified a smaller EC cluster corresponding to pre-HSCs based on expression of genes such as cmyb, cebpa, and gfi1aa. We found two clusters, brain vascular endothelial cells 1 and 2 (BVECs-I and BVECs-II), that co-expressed brain and neuronal-associated genes, including tncb, elavl3, and fmoda. Another cluster likely represents endocardial EC (EEC) based upon co-expression of canonical heart genes, including gata5, hand2, and spock3. A kidney vascular EC (KVECs) cluster co-expressed kidney-associated genes, such as pax2a, ap1m2, and myh9a. Finally, we identified an EC cluster that contained the expression of paraxial mesoderm (PM) signature genes, including igf2a, fbn2b, and fn1a. This PM signature suggested that this population may represent somite-derived endothelial cells (SDECs). We used differential expression analysis among clusters to identify distinct gene programs which are enriched within each specific subset (Figure 1B and Figure 1—source data 1).
To test our clustering predictions and verify the somitic origin of the SDEC cluster, we compared its genetic signature to annotated genes expressed within the somite compartment based on the AmiGO annotation database (Consortium, 2019; Figure 1—source data 2). Our comparative analysis resulted in 32 genes that were commonly expressed within the SDEC fraction and the AmiGO annotation database (Figure 1C’ and C”). Interestingly, several of these 32 genes were enriched within the SDEC cluster (Figure 1C” white boxed genes and 1D). We repeated our analysis on the BVEC-I and KVEC clusters. To validate their cluster identity, we compared their transcriptomes to the annotated brain and kidney genes on AmiGO. Our comparative analysis resulted in 32 and 63 common genes, respectively (Figure 1—figure supplement 2A–C” and Figure 1—figure supplement 3A–C”). Furthermore, as in the SDEC cluster, our analyses confirmed that many of these genes were enriched in both clusters (Figure 1—figure supplement 2D and Figure 1—figure supplement 3D). Thus, our scRNA-seq analysis approach identified various tissue-specific EC subsets in the 22–24 hpf embryo, including HE, pre-HSCs, and a rare population of SDECs.
Identified EC populations can be traced back as early as the tailbud stage
After identifying the discrete genetic signatures of each endothelial cell cluster at 22–24 hpf, we queried if these signatures could also be identified at earlier stages. We followed the same protocol to purify and scRNA-seq samples collected at different developmental stages, namely drl:H2B-dendra+ embryos at tailbud (10 hpf) and 12 somite stages (ss) (15 hpf), as well as etv2:Kaede+ embryos at 15 ss (16.5 hpf) (Figure 2A–D). Next, we crossed-referenced the transcriptomes of all EC subsets from the mixed vasculature 22–24 hpf samples (Figures 1A and 2A) with those of cells purified from earlier tailbud, 12-, and 15-ss samples and clustered them. We identified EC clusters with distinct transcriptomes as early as the tailbud stage, suggesting that EC specification started by the end of gastrulation and the initiation of somitogenesis. While the EC clusters were still heterogeneous prior to the 12 ss stage, we could identify most EC subsets, including HE, pre-HSC, KVECs, BVECs, and SDECs, after the 15 ss when their transcriptomes became more defined (Figure 2B). Together, these data suggest a gradual EC specification process that begins as early as the tailbud stage (Figure 2B–D).
To explore how early tissue-specific EC segregation occurs, we compared the changes in expression patterns of EC populations from early TgBAC(etv2:Kaede)ci6 15 ss and later 22 hpf scRNA-seq data. Specifically, we focused on the changes in expression in the 32 overlapping genes between the scRNA-seq transcriptome data and the AmiGo annotated genes (Figure 1C’ and Figure 1—source data 2). From the 32 SDEC genes compared, 16 genes were upregulated in the etv2:Kaede+ 22 hpf sample compared to the 15 ss sample (Figure 2E and F), whereas 16 genes were downregulated (Figure 2G and H). GO function analysis indicated a role for upregulated genes in arterial EC differentiation, whereas downregulated genes were involved in somitogenesis (Figure 2E and G, respectively). The expression and downregulation of somitic genes within etv2+ ECs between 15 ss and 22 hpf highlight their somitic origin and the loss of myogenic cell fate of these cells as they emerge from the somite compartment.
We repeated our analysis on the BVEC-I and KVEC clusters to validate this observation by comparing samples from etv2:Kaede+ 15 ss and 22 hpf. In the BVEC-I cluster, we found 17 upregulated genes that indicate epithelium development, whereas seven genes indicating brain development were downregulated (Figure 2—figure supplement 1A–D). In the KVEC cluster, we found 18 upregulated genes indicative of epithelial cell differentiation in the kidney and 14 downregulated genes indicating kidney development (Figure 2—figure supplement 1E–H). As with SDECs, both BVEC-I and KVEC cells showed evidence of their commitment to a brain and kidney EC fate as early as 15 ss.
Rare endothelial cells emerge from trunk somites located above the dorsal aorta
Our transcriptome analysis identified that SDECs arise as a distinct EC population as early as the tailbud stage. Next, we aimed to characterize the development of SDECs and their potential to contribute to DA formation at later stages. To label somitic cells and follow their trajectories, we crossed a photoconvertible Tg(actb2:nls-Eos) animal with a Tg(fli1:eGFP)y1 vasculature reporter line (Cruz et al., 2015; Isogai et al., 2003). Resulting embryos ubiquitously express a photoconvertible nuclear GFP that can be instantly converted to RFP by UV exposure. To determine which somites generate SDECs and to quantify how many SDECs are made by each somite pair, we collected embryos at different somitogenesis stages ranging from 4 ss to 18 ss and converted single somite pairs using directed UV exposure via confocal microscopy. We avoided converting any LPM-derived EC on the lateral or ventral sides of the somite (Figure 3A and A’). Embryos were allowed to develop to 32–36 hpf, then imaged and examined for the presence of RFP+ SDECs within the dorsal aorta (Figure 3B). We observed that the trunk somites (numbers 5–18), located above the yolk tube extension generated the most SDECs. In contrast, photoconverted somites anterior to the 4th somite rarely gave rise to SDECs (Figure 3C and D). In sum, most SDECs were generated by somites 5–18. Each somite pair contributed between 0–6 SDECs to the DA (Figure 3D and Figure 3—source data 1). Interestingly, we occasionally noticed SDECs within the roof of the DA, suggesting a migratory path that passes through the roof of the DA before moving to the floor of the DA, consistent with previous reports (Zhao et al., 2022). On rare occasions, we also detected SDECs in the posterior caudal vein or intersomitic vessels, consistent with previous findings (Nguyen et al., 2014).
Somites give rise to etv2+ endothelial cells concomitant with somite epithelialization
Next, we aimed to further characterize the development of SDECs by following etv2:eGFP+ cells and their trajectories from the somite. As the SDEC subset became more distinct between 12–15 ss, we focused our analysis on that developmental time frame. To do so, we used confocal time-lapse imaging of Tg(etv2.1:eGFP)zf372 embryos (Veldman and Lin, 2012) injected with mOrange2:CAAX mRNA. These embryos have early ECs marked by GFP with cell boundaries demarcated by mOrange2. As previously described (Veldman and Lin, 2012), beginning at the 10 ss, we observed a line of etv2:eGFP+ along the most medial part of the LPM (Figure 4A). We did not detect any etv2:eGFP+ in the somites at this stage. Initiating at the 12 ss, we noted etv2:eGFP+ in the lateral lip of the somitic compartment (Figure 4B). After the onset of etv2:eGFP expression, within the next 3 hr, somite-derived etv2:eGFP+ rounded up and delaminated from the somite, then integrated into the cohort of LPM-derived etv2:eGFP+ that migrated toward the midline (Figure 4C–E). Time-lapse imaging of Tg(etv2.1:eGFP)zf372; Tg(phlbd1:Gal4-mCherry) embryos resulted in similar observations that support the emergence of ECs from the somites as early as 12 ss (Figure 4—video 1 and Figure 4—video 2).
In mice and chick, SDECs emerge from the same region as skeletal muscle progenitor cells in the hypaxial dermomyotome compartment (Pouget et al., 2006; Tozer et al., 2007). To examine the spatial origins of SDECs in the zebrafish, we performed double fluorescent in situ hybridization (FISH) for the endothelial marker etv2 and the skeletal muscle progenitor marker pax3a (Relaix et al., 2005). At 12 ss, we observed colocalization of etv2 and pax3a within rare cells in the somite (Figure 4F). However, FISH for etv2 and myod, a marker of differentiated muscle cell types (reviewed in Hernández-Hernández et al., 2017), did not show colocalization (Figure 4G). These results suggest that SDECs emerge within the somite from precursor cells shared with the muscle lineage. We observed that somitic etv2+ were localized specifically within the dermomyotome region, the muscle progenitor cell compartment of the somite. As this was previously shown in other organisms, we suggest that a conserved mechanism of SDEC generation is shared amongst vertebrates (Mayeuf-Louchart et al., 2016; Mayeuf-Louchart et al., 2014; Pouget et al., 2006).
The dermomyotome contains progenitors with muscle and endothelial potential
Since our results above suggested the presence of bipotent progenitors with competence for muscle and endothelial cell differentiation, we sought to determine if blocking skeletal muscle differentiation may lead to enhanced SDEC generation by knocking down mesenchyme homeobox 1 (meox1), an essential regulator of muscle cell formation (Mankoo et al., 2003). To knock down meox1, we used a meox1 translation blocking morpholino (MO) that can block translation of a functioning Meox1 protein but does not affect the transcription of a viable meox1 mRNA; thus, it can still be detected by an RNA probe. Using FISH to label meox1 and etv2 expression, we observed ectopic expression of etv2 in morphant animals and an increased number of SDECs (Figure 5A–F). Like normal SDECs, ectopic etv2:eGFP+ emerged at 12 ss. However, in contrast to wild type, meox1 morphants showed SDEC generation as late as 23 hpf in Tg(phldb1:mCherry); Tg(etv2.1:eGFP)zf372 embryos (Distel et al., 2009; Kobayashi et al., 2014; Figure 5G and H). Thus, meox1 loss of function leads to enhanced and prolonged production of SDECs.
Previous work in mice identified the Notch signaling pathway as a positive regulator of muscle fate in the somite (Mayeuf-Louchart et al., 2014). We tested the requirement for Notch signaling in zebrafish SDEC generation by knocking down the essential notch regulator mindbomb (mib) (Itoh et al., 2003) in Tg(etv2.1:eGFP)zf372 animals and examining SDEC formation between 12–14 ss. While loss of meox1 alone led to the ectopic formation of a few SDECs (Figure 5I and J), knockdown of both mib and meox1 led to a profound expansion of etv2:eGFP+ within the somite (Figure 5K).
To examine the role of Notch signaling more precisely in these cells, we examined notch3fh332/fh332 animals (referred to as notch3-/-) since notch3 is the primary notch receptor in the somites at this stage (Kim et al., 2014). We injected meox1 MO into notch3-/- embryos and performed FISH for meox1 and etv2 (Figure 5L–N). Using this strategy, we observed that in notch3-/- homozygous mutant embryos, the formation of SDECs was not impaired (Figure 5L and M). However, following the combined loss of notch3 and meox1, we observed the ectopic formation of etv2+; meox1+ double positive cells by FISH (Figure 5N and Figure 5—figure supplement 1).
To determine which cell populations were explicitly affected by the loss of notch3, we compared the expression levels of markers for muscle progenitors (pax3a, pax7b), differentiated muscle cells (myod, myog), and differentiated endothelial cells (etv2, fli1a) by qRT-PCR in notch3 mutant and heterozygous embryos at 24 hpf (Figure 5O). In notch3-/- animals, we observed a decrease in the expression of muscle progenitor markers (pax3a, pax7b) concomitant with an increase in the expression of muscle differentiation markers (myod, myog) and endothelial differentiation markers (etv2, fli1a) (Figure 5O). Furthermore, we observed premature expression of the muscle differentiation marker myoHII (Beier et al., 2011; Salucci et al., 2015; Sjöblom et al., 2008) by antibody staining at 48 hpf in notch3-/- animals (Figure 5P and Q). Together, these results suggest that Notch signaling is dispensable for the specification of SDECs but is required for pax3+ and pax7+ skeletal muscle progenitor maintenance, as previously described (Relaix et al., 2005; Zhang et al., 2021). Moreover, the absence of notch3 leads to premature differentiation of muscle and SDEC fates (Figure 5R).
npas4l is required for the formation of SDECs npas4l (cloche) is regarded as the most upstream gene required for blood and endothelial cell specification (Reischauer et al., 2016; Stainier et al., 1995). Therefore, we sought to determine if npas4l is also required for SDEC development. WISH in 12 ss embryos showed a complete absence of etv2 expression along the embryonic A-P axis in Clos9/s9 mutants (referred to as cloche) compared to controls (Figure 6A–B), showing that npas4l function is required for both PLM-derived ECs and PM-derived SDECs. Knockdown of meox1 to increase SDEC differentiation from bipotent somitic progenitors was insufficient to rescue SDEC generation, indicating npas4l function is necessary to generate the SDEC lineage. In addition, suppressing muscle specification in the absence of npas4l was insufficient to induce SDEC fate (Figure 6C–D). Similarly, in Tg(etv2.1:eGFP)zf372; cloche embryos injected with mib MO and meox1 MO, we observed a complete absence of etv2:eGFP+, compared to cloche control siblings (Figure 6—figure supplement 1).
Lastly, we performed qRT-PCR for endothelial and muscle cell genes from 48 hpf cloche mutant and control embryos (Figure 6E). As expected, cloche mutant embryos showed decreased endothelial gene expression (etv2 and fli1a). Interestingly, we also observed a concomitant increase in the expression of skeletal muscle differentiation genes (myod and myogenin) (Figure 6E). Together, these results confirm that npas4l is required for specification of all EC subsets, including SDECs, and suggest a negative regulatory role in muscle differentiation in the shared bipotent progenitors of SDECs and skeletal muscle cells (Figure 6F and G).
Wnt signaling regionalizes the formation of SDECs
Previous work has shown that Wnt signaling is required for the differentiation of muscle progenitors through activation of the required skeletal muscle factor myf5 (reviewed in von Maltzahn et al., 2012). In addition, inhibition of Wnt signaling in early presomitic mesoderm leads to an increase in endothelial cells that can integrate into the zebrafish vasculature (Veldman et al., 2013), suggesting that Wnt signaling may also act later in the somite to balance muscle and endothelial cell production from shared muscle progenitor cells. To investigate whether Wnt signaling is active in meox1+ muscle progenitor cells during SDEC development, we performed FISH for meox1 in the background of Tg(7xTCF-Xla.Siam:GFP)ia4, a destabilized Wnt/TCF reporter line +. At the 12 ss, we observed cells positive for GFP and meox1 (Figure 7A), indicating that Wnt signaling is active while SDEC fate decisions occur. Next, to determine if Wnt inhibition affects SDEC development, we treated Tg(etv2.1:eGFP)zf372 embryos with IWP-L6, a potent inhibitor of Wnt protein secretion (Wang et al., 2013), from 2 to 12 ss. Inhibition of the Wnt signal was confirmed by qRT-PCR for the canonical target gene axin2 (Figure 7B), which led to a decrease in axin2 and an increase in etv2 transcripts. In addition, we observed a slight reduction in meox1 expression; however, it was not statistically significant. By examining serial sections, Wnt inhibition led to the formation of ectopic etv2:eGFP+ in the somites (Figure 7C–D’). Finally, we also determined that Wnt signaling is required for meox1 expression (Figure 7E–F’). Together, these results demonstrate that Wnt signaling compartmentalizes the formation of SDECs within the somite to the most lateral region. In addition, they show that Wnt inhibition results in downregulation of meox1 and the expansion of etv2:eGFP+.
SDECs integrate into the dorsal aorta but do not generate HSPCs
Our results confirmed previous studies in the chick (Pouget et al., 2006) and zebrafish (Nguyen et al., 2014) that ECs emerge from the somites and incorporate into the DA (Figure 3). Next, we were interested in understanding the role of SDECs in the developing zebrafish vasculature and hematopoiesis. To this end, we obtained a PM-specific Gal4 transgenic line, Tg(tbx6:Gal4FF:GFP-nls) (Yabe et al., 2016), which has been shown to recapitulate tbx6 mRNA expression (referred to as tbx6:Gal4). We crossed this line to Tg(UAS-Cre) (Butko et al., 2015) to drive Cre recombinase only within the PM. We then crossed it to a ubiquitously expressed reporter line Tg(actb2:loxP-BFP-loxP-DsRed)sd27 (Kobayashi et al., 2014), which upon genetic recombination, switches from a BFP to DsRed cassette (abbreviated as A2BD for clarity). Upon confocal imaging of the triple transgenic tbx6:Gal4; Tg(UAS-Cre); A2BD animals at 48 hpf, we observed DsRed+ cells within the region of the axial vasculature (Figure 8A–C). To confirm that these cells were SDECs arising from a tbx6+ somitic population, we crossed the tbx6:Gal4; Tg(UAS-Cre) line with a previously published endothelial-specific CFP-to-YFP switch line, TgBAC(kdrl:LOXP-AmCyan-LOXP-ZsYellow) (referred to as kdrl:CSY) (Zhou et al., 2011). Imaging the DA in tbx6:Gal4; Tg(UAS-Cre); kdrl:CSY triple transgenic embryos at 4 dpf showed YFP+ cells corresponding to SDECs to localize preferentially to the trunk region of the dorsal aorta (Figure 8D–G and K). We found this labeling pattern consistent among different embryos over a range of developmental time points.
To complement these results, we performed lineage tracing using an LPM-specific switch line, Tg(drl:CreERT2) (Henninger et al., 2017). Likewise, we crossed it to kdrl:CSY for EC-specific tracing. We incubated Tg(drl:CreERT2); kdrl:CSY embryos with tamoxifen (4-OHT) to induce Cre-based recombination between 8 and 24 hpf. At 4 dpf, confocal imaging showed complementary results to our PM-specific tbx6+ lineage tracing experiments. We observed that drl-derived YFP+ ECs contributed to all regions of the vasculature, as would be expected from an LPM source (Figure 8H–J and K). However, we observed that the contribution of these cells was more robust within the posterior region of the DA and inverse to the contribution of the PM-specific, tbx6-labeled SDECs, which contributed ECs mainly to the trunk region of the vasculature (Figure 8D–K). Together, these results further demonstrate that SDECs integrate into the DA in zebrafish and appear to contribute preferentially to the trunk portion of the dorsal aorta.
Since HSPCs derive from hemogenic endothelium, specifically within the DA, we examined whether or not SDECs can generate HSPCs. We lineage traced PM-derived cells into the kidney, the adult hematopoietic organ in teleosts. We dissected kidneys of adult tbx6:Gal4; Tg(UAS-Cre); A2BD transgenic animals and observed no contribution of switched DsRed+ cells to any hematopoietic lineage (Figure 8L). A ubiquitous vascular-specific transgenic Cre driver was utilized as a positive control for hemogenic endothelium, Tg(kdrl:Cre)s898 (Bertrand et al., 2010; Figure 8—figure supplement 1). Taken together, these results demonstrate that SDECs integrate into the dorsal aorta but do not generate HSPCs.
SDECs support the emergence of HSPCs
Since SDECs did not directly contribute to HSPCs, we wanted to explore their role by analyzing the differential expression of SDECs versus lateral plate mesoderm-derived ECs (LPMDEC) that populated the DA. We compared the SDEC to HE and Pre-HSC clusters at 22–24 hpf. From SDECs, we identified many genes previously attributed to ‘niche’ functions required to induce aortic hemogenic endothelium (Figure 9A; Charbord et al., 2014). These results suggested that SDECs might be acting in a paracrine manner to support hematopoietic induction. If so, we reasoned that changes in the number of SDECs might result in the altered formation of hematopoietic cells. Increasing the numbers of SDECs via enforced expression of etv2 with a muscle-specific mylz2 promoter (Ju et al., 2003) or knockdown of meox1 led to an increase in runx1+ by WISH (Figure 9B–E). In contrast, decreasing the number of SDECs via over-expression of meox1 mRNA led to the loss of runx1+ and gata2b+hematopoietic cells in the DA (Figure 9F–I). These results suggest that SDECs act in a paracrine manner to support the induction of hematopoietic cells from neighboring hemogenic endothelial cells.
Discussion
Here, we present an extensive analysis of how a bipotent skeletal muscle progenitor population gives rise to a rare contingent of ECs that are functionally distinct from endothelial cells specified in the LPM (Garcia-Martinez and Schoenwolf, 1992; Psychoyos and Stern, 1996; Schoenwolf et al., 1992; Selleck and Stern, 1991). Through single-cell analysis, we identify unique molecular signatures sufficient to identify subsets of ECs, including somite-derived endothelial cells (SDECs), brain endothelial cells (BVECs), kidney endothelial cells (KVECs), and endocardium as early as the tailbud stage. Together, these results indicate that cell fate restriction of EC subsets occurs earlier than previously thought. Our results show that SDECs emerge from the trunk dermomyotome during somitogenesis and contribute to the dorsal aorta of the zebrafish embryo. We show that modulators of muscle differentiation, such as meox1, notch, and wnt, regulate SDEC number and the timeframe of SDEC differentiation. Furthermore, our findings elucidate how SDECs are needed to support the hemogenic program required for the formation of HSPCs.
Avian and murine studies have previously identified a population of endothelial cells that arise from the somites (Ambler et al., 2001; Esner et al., 2006; Kardon et al., 2002; Mayeuf-Louchart et al., 2016; Mayeuf-Louchart et al., 2014; Pardanaud et al., 1996; Pouget et al., 2006; Wilting et al., 1995; Yvernogeau et al., 2012). Previously, the common belief was that in zebrafish, all endothelial and blood cells arise exclusively from the LPM (Childs et al., 2002; Jin et al., 2005; Kohli et al., 2013; Lawson and Weinstein, 2002; Zhang and Rodaway, 2007). Furthermore, angioblasts are thought to be specified as an equipotent population that undergoes progressive cell-fate restriction only upon migration to the midline (Kobayashi et al., 2014; Lawson et al., 2001). To our knowledge, only two studies have shown that paraxial mesoderm can generate endothelial cells in zebrafish (Martin and Kimelman, 2012; Nguyen et al., 2014). The first noted that presomitic mesoderm could generate ECs that incorporate into caudal blood vessels following early inhibition of the Wnt signaling pathway after gastrulation (Martin and Kimelman, 2012). The second study highlighted a rare population of ECs generated in the somite that migrated to the developing dorsal aorta (Nguyen et al., 2014). Each of these studies was based on retrospective analyses, making it difficult to ascertain precisely when and where SDECs arise and the extent of their contribution. In addition, studies by Murayama et al., 2023; Murayama et al., 2015 have shown that a population of stromal cells important in the maintenance of zebrafish HSPCs derive from the sclerotome, the ventrolateral domain of the somite. These do not appear to give rise to ECs but rather to stromal cells that reside adjacent to blood vessels. As we have occasionally observed rare cells outside of the vasculature in our lineage-tracing studies, some of our marked cells may include this stromal cell subset. Further studies are required to ascertain how this population of somite-derived cells compares to SDECs in the trophic support of hematopoietic precursors.
Labeling individual somite pairs using Tg(actb2:nls-Eos), a pan-nuclear photoconvertible transgenic line, allowed us to quantify the contribution of SDECs from trunk somites. Using Tg(etv2.1:eGFP)zf372, with etv2 as one of the earliest markers of endothelial cell fate acquisition, allowed us to visualize precisely when and where SDECs emerge. Lastly, we could permanently and differentially mark ECs from the LPM or PM using an endothelial-specific switch transgene and follow their respective contributions to the vasculature. Together, both temporal and indelible lineage tracing approaches showed that somites 5–18 have potential to generate SDECs. Emergence of SDECs initiated at the 12-somite stage, concomitant with the epithelization of the somite and the migration of LPM-derived ECs to the midline. Notably, contribution of SDECs to adult hematopoietic tissues was not observed, consistent with many previous studies concluding that the adult hematopoietic program is LPM-derived (Henninger et al., 2017; Jin et al., 2007; Murayama et al., 2006). These results are consistent with a previous study by Nguyen et al., 2014 that identified a rare population of EC precursors born in a central location within somite, which they termed the ‘endotome’, that incorporated into the axial vessels but did not appear to generate blood (Nguyen et al., 2014). By contrast, in our study, we observed that somite-derived etv2+ ECs arise from the hypaxial dermomyotome, which is consistent with earlier observations in avian and mouse embryos (Eichmann et al., 1993; Ema et al., 2006; Pouget et al., 2006; Tozer et al., 2007). Notably, the hypaxial dermomyotome comprises bipotent-skeletal muscle progenitors that generate ECs or myoblasts upon local cues. Since these bipotent progenitors are present in the territory where the hypaxial muscle precursors originate (Christ and Ordahl, 1995), we believe a more appropriate term for this somitic compartment would be the ‘myo-endotome’.
To better understand the signaling pathways regulating SDEC emergence, we tested the roles of notch, wnt, and npas4l. We find that Notch signaling is essential for the maintenance of muscle progenitors, as previously shown in mice (Schuster-Gossler et al., 2007). Loss of the Notch signaling components mib or notch3 resulted in a premature depletion of skeletal muscle progenitor markers and a concomitant increase of differentiated muscle and endothelial cell gene expression. Our results thus suggest that Notch signaling is required to maintain muscle progenitors rather than determining cell fate choice between muscle and endothelium. However, further work is needed to investigate whether the reduced expression of skeletal muscle progenitor markers results from progenitor cell depletion or simply reflects lower expression levels within these cells. In addition, Wnt signaling is necessary for the regionalization of muscle and endothelial cells within the somite, in agreement with previous work (Borello et al., 2006; Martin and Kimelman, 2012). We find that during homeostasis, a few etv2:eGFP+ delaminate from the hypaxial dermomyotome. However, combined loss of Notch signaling and meox1 shows that EC competence within the somite encompasses a much larger territory, extending as far as the median dermomyotome. Interestingly, we find that endothelial potential within bipotent skeletal muscle progenitors is npas4l-dependent. These results are surprising since npas4l has been shown to be the earliest gene required for EC specification and was thought to be restricted to the LPM starting at the tailbud stage (Reischauer et al., 2016). Equally surprising, a recent paper demonstrated that ECs populating the zebrafish caudal hematopoietic tissue are endoderm-derived (Nakajima et al., 2023). Similar to SDECs, this novel population is likewise dependent upon npas4l function (Nakajima et al., 2023). Together, these results suggest that npas4l is required to generate endothelial cell fate regardless of tissue origins.
Recent studies have elucidated a role for non-hemogenic ECs in supporting the hemogenic program and in amplifying HSPC number in various systems (Butler et al., 2012; Butler et al., 2010; Butler and Rafii, 2012; Gori et al., 2017; Gori et al., 2015; Guo et al., 2017; Hadland et al., 2015; Kim et al., 2014; Kobayashi et al., 2010; Lis et al., 2017; Raynaud et al., 2013; Sandler et al., 2014). Interestingly, several groups have successfully differentiated endothelial precursors into transplantable HSPCs, a process that requires an endothelial ‘niche’ population (Lis et al., 2017; Sandler et al., 2014). Our work indicates that changes to the number of SDECs via downregulation or overexpression of meox1 resulted in increased or decreased hematopoietic factors, such as runx1 and gata2b, respectively. This finding and evidence that SDECs lack hemogenic potential suggest that SDECs are important as niche support cells to induce HSPC formation. Future work will elucidate the molecular nature of this interaction and will forward efforts to instruct HSPC fate from human pluripotent precursors.
Materials and methods
Zebrafish husbandry
Request a detailed protocolWild-type AB* and transgenic TgBAC(etv2:Kaede)ci6 (Kohli et al., 2013) referred to as etv2:Kaede, Tg(fli1:DsRed)um13 (Villefranc et al., 2007) referred to as fli1:DsRed, Tg(tp1:GFP)um14 (Parsons et al., 2009) referred to as tp1:GFP, Tg(drl:H2B-dendra) (Mosimann et al., 2015) referred to as drl:H2B-dendra, Tg(actb2:nls-Eos) (Cruz et al., 2015) referred to as actb2:nls-Eos, Tg(fli1:eGFP)y1 (Isogai et al., 2003) referred to as fli1:eGFP, Tg(etv2.1:eGFP)zf372 (Veldman and Lin, 2012) referred to as etv2:eGFP, Tg(phlbd1:Gal4-mCherry) (Distel et al., 2009) referred to as phlbd1:mCherry, notch3fh332/fh332 (Alunni et al., 2013) referred to as notch3-/-, Clos9/s9 (Reischauer et al., 2016) referred to as npas4l-/- (cloche), Tg(7xTCF-Xla.Siam:GFP)ia4 (Moro et al., 2012) referred to as TCF:d2GFP, Tg(tbx6:Gal4FF:GFP-nls) (Yabe et al., 2016) referred to as tbx6:Gal4; Tg(UAS-Cre)* (Butko et al., 2015) referred to as UAS-Cre, Tg(actb2:loxP-BFP-loxP-DsRed)sd27 (Kobayashi et al., 2014) referred to as A2BD, TgBAC(kdrl:LOXP-AmCyan-LOXP-ZsYellow) (Zhou et al., 2011) referred to as kdrl:CSY, Tg(drl:CreERT2) (Henninger et al., 2017) referred to as drl:CreERT2, and Tg(kdrl:Cre)s898 (Bertrand et al., 2010) referred to as kdrl:Cre. Zebrafish embryos and adult fish were raised in a circulating aquarium system (Aquaneering) at 28 °C and maintained in accordance with UCSD Institutional Animal Care and Use Committee (IACUC) guidelines (Protocol number: S04168), the Association for Assessment and Accreditation of Laboratory Animal Care International (AAALAC) (Accreditation number: 000503), and the Public Health Service (PHS) Policy on Humane Care and Use of Laboratory Animals (Assurance number: A3033-1 or D16-00020). All experiments were performed under approved methods of anesthesia, and every effort was made to minimize suffering.
FACS
Request a detailed protocolTgBAC(etv2:Kaede)ci6 embryos were collected at 15 ss and 22 hpf. Tg(drl:H2B-dendra) embryos were collected at tailbud stage, 12 ss, and 22 hpf. Tg(fli1:DsRed)um13; Tg(tp1:GFP)um14 embryos were collected at 24 hpf. Samples were separately dissociated in phosphate-buffered saline (PBS) supplemented with 10% fetal bovine serum (FBS) by homogenizing them with a sterile plastic pestle or pipette. Dissociated cells were filtered through a 35 μm nylon cell strainer (Falcon 2340) and then rinsed with PBS with 10% FBS. Propidium iodide (Sigma) was added (1 μg ml−1) to exclude dead cells and debris. FACS was performed using GFP, Dendra, and DsRed fluorescence gating with a FACS Aria II flow cytometer (Beckton Dickinson).
Single-cell RNA sample preparation
Request a detailed protocolAfter FACS, total cell concentration and viability were ascertained using a TC20 Automated Cell Counter (Bio-Rad). Samples were resuspended in 1XPBS with 10% bovine serum albumin (BSA) at a concentration between 8E2-3E3 cells per ml. Samples were loaded on the 10 X Chromium system and processed as per manufacturer’s instructions (10 X Genomics). Single-cell libraries were prepared per the manufacturer’s instructions using the Single Cell 3’ Reagent Kit v2 (10 X Genomics). Single-cell RNA-seq libraries and barcode amplicons were sequenced on an Illumina HiSeq2500 platform.
Single-cell RNA sequencing analysis
Request a detailed protocolThe Chromium 3’ sequencing libraries were generated using Chromium Single Cell 3’ Chip kit v3 and sequenced with an Illumina HiSeq2500 platform. The Ilumina FASTQ files were used to generate filtered matrices using CellRanger 3.0.0 (10X Genomics) with default parameters and imported into R for exploration and statistical analysis using a Seurat package V4 (La Manno et al., 2018). Counts were normalized according to total expression, multiplied by a scale factor (10,000), and log-transformed. For cell cluster identification and visualization, gene expression values were also scaled according to highly variable genes after controlling for unwanted variation generated by sample identity. Cell clusters were identified based on UMAP of the first 20 principal components of PCA using Seurat’s method, Find Clusters, with an original Louvain algorithm and resolution parameter value 0.5. To find cluster marker genes, Seurat’s method, FindAllMarkers was used. Only genes exhibiting a significant p-value of (<0.022), and a minimal average absolute log2-fold change of 1 between each cluster and the rest of the dataset, were considered differentially expressed. To merge individual 22–24 hpf datasets and remove batch effects, Seurat v3 Integration and Label Transfer standard workflow (Stuart et al., 2019) was used. To project samples of embryonic stage <22 hpf onto the reference 22–24 hpf dataset, Seurat’s MapQuery function was used.
UV-driven photoconversions
Request a detailed protocolTg(actb2:nls-Eos); Tg(fli1:eGFP)y1 embryos were collected and staged at developmental stages ranging from 4 to 18 ss in a glass bottom Petri dish that was precast using a 2% low melting agarose (Fisher Bioreagents, BP160-500) and a custom-made stamp (Idylle-lab) to created semi-spherical wells. Embryos were then staged and photoconverted with newly emerging somites directly at the bottom of the dish using a Leica SP8 confocal microscope. Recently developed somite pairs were selected by setting a region of interest, then photoconverted by applying maximum UV exposure (405 nm; 100% laser power) using a 20 x objective lens and high digital zoom (>X4) for 90 s. To validate the quality and specificity of the somite photoconversion, images were taken of each channel, eGFP (488 nm; 10% laser power) and mCherry (564 nm; 25% laser power). Photoconverted zebrafish embryos were then transferred to a petri dish with 1-phenyl 2-thiourea (PTU) media and incubated at 28.5 °C in the dark. At 32 hpf, embryos were embedded in a glass bottom petri dish covered by 1% agarose and imaged on a Leica SP8 confocal microscope using resonant scanning. Images of the dorsal aorta were taken in a sequential stack manner for eGFP (488 nm; 10% laser power) and mCherry (564 nm; 35% laser power) using a 20 x objective lens.
Quantification of SDECs emerging from somites
Request a detailed protocolMaximum intensity projection (MIP) images were taken for each converted embryo to identify the converted somite pair. Then, SDECS were identified and quantified in the dorsal aorta by examining individual z-stacks and visualizing colocalization of fli1+; actb2:nlsEosRFP converted cells in Tg(actb2:nls-Eos); Tg(fli1:eGFP)y1 embryos. Particularly in fli1:eGFP- embryos, we identified SDECs by observing actb2:nlsEosRed cells on the floor of the DA based on the brightfield channel. Statistical analysis of the number of SDECs per somite was performed using Prism 9 (GraphPad).
Statistical analysis and graphic display
Request a detailed protocolGraphic data representation and statistical analysis were performed using Prism 9 (Graphpad Version 9.5.1). In Figure 3, All graphs with single dot scatter plots indicate the mean as column height, and standard error of the mean (SEM) is represented as error bars. In Figures 5—7, the statistics of the qPCR were performed using unpaired, two-tailed t-tests (n=3), and standard error of the mean (SEM) is represented as error bars.
Microscopy and timelapse videos
Request a detailed protocolEmbryos were embedded and sectioned as described by Kobayashi et al., 2014. Live transgenic embryos and flat-mount or whole-mount two-color double FISH samples were imaged using confocal microscopy (SP5 or SP8, Leica). Since the emission spectra of CFP and YFP overlap, images were captured in two separate sequences to filter overlapping signals by limiting the PMT detection. For CFP a 476 nm laser was used, and the PMT detector was set to collect signals ranging between 480 nm and 505 nm. For YFP, a 514 nm laser was used, and the PMT detector was set to collect signals ranging between 520 nm and 570 nm.
Microinjections of morpholinos, mRNA, and plasmids
Request a detailed protocolAntisense morpholinos (MOs; Gene Tools, LLC) were diluted as 1- or 3 mM stock in H2O. meox1-MO (CTGGCTGACTGTTCCATACTGAAGA) and mib-MO (GCAGCCTCACCTGTAGGCGCACTGT) were injected at the 1–2 cell stage of development. mOrange-CAAX mRNA was synthesized from linearized mOrange-CAAX with the mMessege mMaching kit (Ambion). 100 pg mOrange-CAAX mRNA was injected into one- to two-cell stage embryos. A total of 25 pg of mylz2:etv2 construct was coinjected with 50 pg of transposase mRNA. All microinjections were performed with the indicated RNA or MO concentration in a 1 nl using a PM 1000 cell microinjector (MDI).
WISH
Request a detailed protocolWhole-mount single or double enzymatic in situ hybridization was performed on embryos fixed overnight with 4% paraformaldehyde (PFA) in PBS. Fixed embryos were washed briefly in PBS and transferred into methanol for storage at −20 °C. Embryos were rehydrated stepwise through methanol in PBS +0.1% Tween 20 (PBT). Rehydrated embryo samples were then incubated with 10 μg ml−1 proteinase K in PBT for 5 min for 5–10 somite stage (12–15 hpf) embryos and 15 min for 24–36 hpf embryos. After proteinase K treatment, samples were washed in PBT and refixed in 4% PFA for 20 min at room temperature. After washes in two changes of PBT, embryos were prehybridized at 65 °C for 1 hr in hybridization buffer (50% formamide, 5 x saline-sodium citrate (SSC), 500 μg ml−1 torula tRNA, 50 μg ml−1 heparin, 0.1% Tween 20, 9 mM citric acid (pH 6.5)). Samples were then hybridized overnight in hybridization buffer including digoxigenin (DIG)- or fluorescein-labeled RNA probe. After hybridization, experimental samples were washed stepwise at 65 °C for 15 min each in hybridization buffer in 2×SSC mix (75%, 50%, 25%), followed by two washes with 0.2×SSC for 30 min each at 65 °C. Further washes were performed at room temperature for 5 min each with 0.2×SSC in PBT (75%, 50%, 25%). Samples were incubated in PBT with 2% heat-inactivated goat serum and 2 mg ml−1 BSA (block solution) for 1 hr and then incubated overnight at 4 °C in block solution with diluted DIG-antibodies (1:5,000) conjugated with alkaline phosphatase (AP) (Roche). To visualize WISH signal, samples were washed three times in AP reaction buffer (100 mM Tris, pH 9.5, 50 mM MgCl2, 100 mM NaCl, and 0.1% Tween 20) for 5 min each and then incubated in the AP reaction buffer with NBT/BCIP substrate (Roche).
For two-color double FISH, embryos were blocked in maleic acid buffer (MAB; 150 mM maleic acid, 100 mM NaCl, pH 7.5) with 2% Roche blocking reagent (MABB) for 1 hr at room temperature, after hybridizing at 65 °C with probes as described above. Embryos were incubated overnight at 4 °C in MABB with anti-fluorescein POD. (Roche) at a 1:500 dilution. After four washes in MABB for 20 min, each followed by washes in PBS at room temperature, embryo samples were incubated in TSA Plus Fluorescein Solution (Perkin Elmer) for 1 hr. Embryos were washed for 10 min each in methanol in PBS (25%, 50%, 75%, 100%). Embryos were incubated in 1% H2O2 in methanol for 30 min at room temperature and washed for 10 min each in methanol in PBS (75%, 50%, 25%) and 10 min in PBS. After blocking for 1 hr in MABB, embryos were incubated overnight at 4 °C in MABB with anti-DIG POD. (Roche) at a 1:1000 dilution. As described above, samples were washed and incubated in TSA Plus CY3 solution (Perkin Elmer). Embryos were washed thrice for 10 min each in PBT and refixed in 4% PFA after the complete staining. Antisense RNA probes for the following genes were prepared using probes containing DIG or fluorescein-labeled UTP: etv2 (Clements et al., 2011), pax3a (Minchin et al., 2013), myod (Minchin et al., 2013), meox1 (Nguyen et al., 2014), runx1 (Burns et al., 2005), and gata2b (Butko et al., 2015).
Quantitative real-time RT-PCR
Request a detailed protocolTotal RNA was collected from whole embryos (~20 embryos) using TRIzol reagent (Ambion, Life Technologies) and isolated from notch3fh332/fh332 mutant embryos or control siblings at 48 hpf; cloche mutant embryos or control siblings at 48 hpf; or IWP2 treated and control embryos at 15 hpf with the RNeasy kit (Qiagen). cDNA was generated from total RNA with iScript cDNA synthesis kit (Bio-Rad). The following primers were used for cDNA quantification: ef1α (forward, 5′- GAGAAGTTCGAGAAGGAAGC –3′; reverse, 5′- CGTAGTATTTGCTGGTCTCG –3′), etv2 (forward, 5′- CGAGGTTCTGGTAGGTTTGAG –3′; reverse, 5′- GCACAAAGGTCATGTTCTCAC –3’), fli1a (forward, 5′- CGTCAAGCGAGAGTATGACC –3′; reverse, 5′- AGTTCATCTGAGACGCTTCG –3’), myod1 (forward, 5′- GAAGACGGAACAGCTATGAC –3′; reverse, 5′- GGAGTCTCTGTGGAAATTCG –3’), myf5 (forward, 5′- CCAGACAGTCCAAACAACAGACC –3′; reverse, 5′- TGAGCAAGCAGTGTGAGTAAGCG –3’), pax3a (forward, 5′- ATTCCTTGGAGGTCTCTACG –3′; reverse, 5′- CTACTATCTTGTGGCGGATG –3’), pax7b (forward, 5′- CAGTATTGACGGCATTCTGGGAG –3′; reverse, 5′- TCTCTGCTTTCTCTTGAGCGGC –3’), myf5 (forward, 5′- GAATAGCTACAACTTTGACG –3′; reverse, 5′- GTAAACTGGTCTGTTGTTTG –3’), myog (forward, 5′- GTGGACAGCATAACGGGAACAG –3′; reverse, 5′-TCTGAAGGTAACGGTGAGTCGG-3’).
Immunofluorescence
Request a detailed protocolWhole-mount immunofluorescence staining for myoHII was performed on notch3fh332/fh332 embryos at 48 hpf using a MF-20 (DSHB) antibody. Immunofluorescence was performed as previously described in Alexander et al., 1998.
Pharmacological treatment
Request a detailed protocolIWR (Tocris) was dissolved in DMSO at a concentration of 10 mM. Zebrafish AB* embryos were incubated in 10 μM IWR solution from 11 to 15 hpf, then fixed with 4% PFA or snap frozen for RNA collection.
Analysis of whole kidney marrow cells
Request a detailed protocolWhole kidney marrow cells (WKM) were prepared as previously described in Kobayashi et al., 2008. Briefly, adult tbx6:Gal4; Tg(UAS-Cre); A2BD or kdrl:Cre; A2BD transgenic animals were anesthetized and sacrificed following IACUC regulations. The kidney of individual transgenic animals was dissected then minced with scissors. Hematopoietic cells were mechanically dissociated by trituration in PBS supplemented with 2% of FBS and filtered through 40 µm nylon mesh. Prior to flow cytometry analysis, samples were stained with a viability dye (5 nM of Sytox Red, ThermoFisher) to exclude dead cells from the analysis. Flow cytometry was performed on a FACS Aria II (BD Bioscience). Collected data were analyzed using FlowJo software.
Plasmid construction
Request a detailed protocolPlasmid expression constructs were generated by standard means using PCR from cDNA libraries generated from zebrafish larvae at 24 hpf. These were cloned into pCS2 +downstream of an SP6 promoter. Previously described mylz2 promoter (Ju et al., 2003) was cloned upstream of etv2 CDS flanked by Tol2 sites.
Replicates
Request a detailed protocolAll experiments assessing phenotype and expression patterns were replicated in at least two independent experiments. Embryos were collected from independent crosses, and experimental processing (staining or injections) was carried out on independent occasions. Exceptions to this include the analysis of the kidney hematopoietic cells by flow from tbx6:Gal4; Tg(UAS-Cre); A2BD and the imaging of tbx6:Gal4; Tg(UAS-Cre); kdrl-CSY where many embryos were processed, and the corresponding n are reported in the associated figure legends. For the analysis of qRT-PCR, three independent experiments were performed per condition. One sequencing run for the single-cell sequencing experiments was performed per timepoint from a pool of at least 100 embryos within the corresponding transgenic background.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. scRNA-seq data are available in ArrayExpress under accession number E-MTAB-13196.
-
Dryad Digital RepositorySingle Cell Sequencing Data for Endothelial Cell Types.https://doi.org/10.6075/J0GB22J0
-
ArrayExpressID E-MTAB-13196. Single-cell RNA-seq (10x Chromium) of FACS-sorted drl:H2B-Dendra2, etv2:Kaede, fli1:DsRed, and tp1:eGFP transgenic zebrafish embryos.
References
-
Phenotypic heterogeneity of the endothelium: II. Representative vascular bedsCirculation Research 100:174–190.https://doi.org/10.1161/01.RES.0000255690.03436.ae
-
Endothelial differentiation: molecular mechanisms of specification and heterogeneityArteriosclerosis, Thrombosis, and Vascular Biology 31:1476–1484.https://doi.org/10.1161/ATVBAHA.111.228999
-
Myogenic differentiation of mesenchymal stem cells co-cultured with primary myoblastsCell Biology International 35:397–406.https://doi.org/10.1042/CBI20100417
-
Hematopoietic stem cell fate is established by the Notch-Runx pathwayGenes & Development 19:2331–2342.https://doi.org/10.1101/gad.1337005
-
Generation of a vascular niche for studying stem cell homeostasisMethods in Molecular Biology 904:221–233.https://doi.org/10.1007/978-1-61779-943-3_18
-
Vegf signaling promotes vascular endothelial differentiation by modulating etv2 expressionDevelopmental Biology 424:147–161.https://doi.org/10.1016/j.ydbio.2017.03.005
-
Patterning of angiogenesis in the zebrafish embryoDevelopment 129:973–982.https://doi.org/10.1242/dev.129.4.973
-
Early stages of chick somite developmentAnatomy and Embryology 191:381–396.https://doi.org/10.1007/BF00304424
-
The Gene ontology resource: 20 years and still going strongNucleic Acids Research 47:D330–D338.https://doi.org/10.1093/nar/gky1055
-
Positional control of mesoderm movement and fate during avian gastrulation and neurulationDevelopmental Dynamics 193:249–256.https://doi.org/10.1002/aja.1001930305
-
Vascular niche promotes hematopoietic multipotent progenitor formation from pluripotent stem cellsThe Journal of Clinical Investigation 125:1243–1254.https://doi.org/10.1172/JCI79328
-
Endothelial cells promote expansion of long-term engrafting marrow hematopoietic stem and progenitor cells in primatesStem Cells Translational Medicine 6:864–876.https://doi.org/10.5966/sctm.2016-0240
-
Endothelial jagged-2 sustains hematopoietic stem and progenitor reconstitution after myelosuppressionThe Journal of Clinical Investigation 127:4242–4256.https://doi.org/10.1172/JCI92309
-
Endothelium and NOTCH specify and amplify aorta-gonad-mesonephros-derived hematopoietic stem cellsThe Journal of Clinical Investigation 125:2032–2045.https://doi.org/10.1172/JCI80137
-
Molecular control of endothelial cell behaviour during blood vessel morphogenesisNature Reviews. Molecular Cell Biology 12:551–564.https://doi.org/10.1038/nrm3176
-
The myogenic regulatory factors, determinants of muscle development, cell identity and regenerationSeminars in Cell & Developmental Biology 72:10–18.https://doi.org/10.1016/j.semcdb.2017.11.010
-
Angiogenic network formation in the developing vertebrate trunkDevelopment 130:5281–5290.https://doi.org/10.1242/dev.00733
-
In vivo imaging of embryonic vascular development using transgenic zebrafishDevelopmental Biology 248:307–318.https://doi.org/10.1006/dbio.2002.0711
-
Regulation of endothelial cell differentiation and specificationCirculation Research 112:1272–1287.https://doi.org/10.1161/CIRCRESAHA.113.300506
-
Chamber identity programs drive early functional partitioning of the heartNature Communications 6:8146.https://doi.org/10.1038/ncomms9146
-
Embryonic origins and assembly of blood vesselsThe American Review of Respiratory Disease 140:1097–1103.https://doi.org/10.1164/ajrccm/140.4.1097
-
Notch-responsive cells initiate the secondary transition in larval zebrafish pancreasMechanisms of Development 126:898–912.https://doi.org/10.1016/j.mod.2009.07.002
-
Vascular heterogeneity and specialization in development and diseaseNature Reviews. Molecular Cell Biology 18:477–494.https://doi.org/10.1038/nrm.2017.36
-
Sclerotomal origin of vascular smooth muscle cells and pericytes in the embryoDevelopmental Biology 315:437–447.https://doi.org/10.1016/j.ydbio.2007.12.045
-
VasculogenesisAnnual Review of Cell and Developmental Biology 11:73–91.https://doi.org/10.1146/annurev.cb.11.110195.000445
-
Mesoderm movement and fate during avian gastrulation and neurulationDevelopmental Dynamics 193:235–248.https://doi.org/10.1002/aja.1001930304
-
Notch signalling and the regulation of angiogenesisCell Adhesion & Migration 1:104–106.https://doi.org/10.4161/cam.1.2.4488
-
Alpha-actinin structure and regulationCellular and Molecular Life Sciences 65:2688–2701.https://doi.org/10.1007/s00018-008-8080-8
-
Etsrp/Etv2 is directly regulated by Foxc1a/b in the Zebrafish angioblastCirculation Research 110:220–229.https://doi.org/10.1161/CIRCRESAHA.111.251298
-
Gateway compatible vectors for analysis of gene function in the zebrafishDevelopmental Dynamics 236:3077–3087.https://doi.org/10.1002/dvdy.21354
-
Hedgehog signaling: iguana debuts as a nuclear gatekeeperCurrent Biology 14:R668–R670.https://doi.org/10.1016/j.cub.2004.08.015
-
Wnt signaling in myogenesisTrends in Cell Biology 22:602–609.https://doi.org/10.1016/j.tcb.2012.07.008
-
The development of highly potent inhibitors for porcupineJournal of Medicinal Chemistry 56:2700–2704.https://doi.org/10.1021/jm400159c
-
Angiogenic potential of the avian somiteDevelopmental Dynamics 202:165–171.https://doi.org/10.1002/aja.1002020208
-
ETS factors regulate Vegf-dependent arterial specificationDevelopmental Cell 26:45–58.https://doi.org/10.1016/j.devcel.2013.06.007
Article and author information
Author details
Funding
National Institute of Diabetes and Digestive and Kidney Diseases (R01-DK074482)
- David Traver
- Pankaj Sahai-Hernandez
NIH Office of the Director (TR01-OD026219)
- David Traver
National Institutes of Health (T32-HL086344)
- Pankaj Sahai-Hernandez
California Institute for Regenerative Medicine (EDUC4-12804)
- Shai Eyal
American Heart Association (19POST34380328)
- Ondrej Svoboda
Ministry of Education of the Czech Republic (RVO: 68378050-KAV-NPUI)
- Ondrej Svoboda
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank members of the Traver laboratory for helpful discussions and Stephanie Grainger for discussions and careful editing of the manuscript. We thank the UCSD Institute for Genomic Medicine sequencing core for supporting the scRNA-seq sample preparation and sequencing. We would like to thank Dr. Peng Guo from the Nikon Imaging Center and Jennifer Santini from the Microscopy core at UCSD for their support. We want to thank Anthony Santella from the Bao laboratory at Memorial Sloan Kettering Cancer Center and Kurt Weiss from the Huisken laboratory at the Morgridge Institute for Research for their support and expert advice. Funding This work was supported by the National Institute of Health [R01-DK074482 to DT and PS-H, TR01-OD026219 to DT]; National Institute of Health [T32-HL086344 to PS-H]; California Institute of Regenerative Medicine [EDUC4-12804 to SE]; American Heart Association [19POST34380328 to OS]; and Ministry of Education of the Czech Republic [RVO: 68378050-KAV-NPUI to OS].
Ethics
Zebrafish embryos and adult fish were raised in a circulating aquarium system (Aquaneering) at 28 °C and maintained by UCSD Institutional Animal Care and Use Committee (IACUC) guidelines (Protocol number: S04168), the Association for Assessment and Accreditation of Laboratory Animal Care International (AAALAC) (Accreditation number: 000503), and the Public Health Service (PHS) Policy on Humane Care and Use of Laboratory Animals (Assurance number: A3033-1 or D16-00020). All experiments were performed under approved methods of anesthesia, and every effort was made to minimize suffering.
Copyright
© 2023, Sahai-Hernandez, Pouget, Eyal 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
-
- 865
- views
-
- 119
- downloads
-
- 10
- 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
Cilia defects lead to scoliosis in zebrafish, but the underlying pathogenic mechanisms are poorly understood and may diverge depending on the mutated gene. Here, we dissected the mechanisms of scoliosis onset in a zebrafish mutant for the rpgrip1l gene encoding a ciliary transition zone protein. rpgrip1l mutant fish developed scoliosis with near-total penetrance but asynchronous onset in juveniles. Taking advantage of this asynchrony, we found that curvature onset was preceded by ventricle dilations and was concomitant to the perturbation of Reissner fiber polymerization and to the loss of multiciliated tufts around the subcommissural organ. Rescue experiments showed that Rpgrip1l was exclusively required in foxj1a-expressing cells to prevent axis curvature. Genetic interactions investigations ruled out Urp1/2 levels as a main driver of scoliosis in rpgrip1 mutants. Transcriptomic and proteomic studies identified neuroinflammation associated with increased Annexin levels as a potential mechanism of scoliosis development in rpgrip1l juveniles. Investigating the cell types associated with annexin2 over-expression, we uncovered astrogliosis, arising in glial cells surrounding the diencephalic and rhombencephalic ventricles just before scoliosis onset and increasing with time in severity. Anti-inflammatory drug treatment reduced scoliosis penetrance and severity and this correlated with reduced astrogliosis and macrophage/microglia enrichment around the diencephalic ventricle. Mutation of the cep290 gene encoding another transition zone protein also associated astrogliosis with scoliosis. Thus, we propose astrogliosis induced by perturbed ventricular homeostasis and associated with immune cell activation as a novel pathogenic mechanism of zebrafish scoliosis caused by cilia dysfunction.
-
- Developmental Biology
- Stem Cells and Regenerative Medicine
Stimulation of pancreatic beta cell regeneration could be a therapeutic lead to treat diabetes. Unlike humans, the zebrafish can efficiently regenerate beta cells, notably from ductal pancreatic progenitors. To gain insight into the molecular pathways involved in this process, we established the transcriptomic profile of the ductal cells after beta cell ablation in the adult zebrafish. These data highlighted the protein phosphatase calcineurin (CaN) as a new potential modulator of beta cell regeneration. We showed that CaN overexpression abolished the regenerative response, leading to glycemia dysregulation. On the opposite, CaN inhibition increased ductal cell proliferation and subsequent beta cell regeneration. Interestingly, the enhanced proliferation of the progenitors was paradoxically coupled with their exhaustion. This suggests that the proliferating progenitors are next entering in differentiation. CaN appears as a guardian which prevents an excessive progenitor proliferation to preserve the pool of progenitors. Altogether, our findings reveal CaN as a key player in the balance between proliferation and differentiation to enable a proper beta cell regeneration.