Abstract
In vertebrates, the four transcription factors Sox2, c-Myc, Pou5f1 and Klf4 are involved in the differentiation of several tissues during vertebrate embryogenesis; moreover, they are normally co-expressed in embryonic stem cells and play roles in pluripotency, self-renewal, and maintenance of the undifferentiated state in adult cells. The in vitro forced co-expression of these factors, named Yamanaka factors (YFs), induces pluripotency in human or mouse fibroblasts. Botryllus schlosseri is a colonial tunicate undergoing continuous stem cell-mediated asexual development, providing a valuable model system for the study of pluripotency in the closest living relatives of vertebrates. In this study, we identified B. schlosseri orthologs of human Sox2 and c-Myc genes, as well as the closest homologs of the vertebrate-specific Pou5f1 gene, through an in-depth evolutionary analysis of the YF gene families in tunicates and other deuterostomes. Then, we studied the expression of these genes during the asexual cycle of B. schlosseri using in situ hybridization in order to investigate their possible involvement in tissue differentiation and in pluripotency maintenance. Our results show a shared spatio-temporal expression pattern consistent with the reported functions of these genes in invertebrate and vertebrate embryogenesis. Moreover, Myc, SoxB1 and Pou3 were expressed in candidate stem cells residing in their niches, while Pou2 was found expressed exclusively in the immature previtellogenic oocytes, both in gonads and circulating in the colonial vascular system. Our data suggest that Myc, SoxB1 and Pou3 may be individually involved in the differentiation of the same territories seen in other chordates, and that, together, they may play a role in stemness even in this colonial ascidian.
Introduction
Stem cells (SCs) are undifferentiated cells capable of self-renewal and with the ability to differentiate into specialised cells. Given their pivotal role in animal development and regeneration, understanding the molecular mechanisms controlling their potency and evolution is of critical importance. Studies in SC biology are mainly based on mammals or other vertebrates, despite their limited regenerative capacities (Zhao et al., 2016). Tunicates are the sister group of vertebrates (Delsuc et al., 2006), and colonial tunicates are the only chordates able to regenerate entire organisms from circulating SCs (Rinkevich et al., 1995; Sunanaga et al., 2006; Voskoboynik et al., 2007; Brown et al., 2009; Kassmer et al., 2019; Manni et al., 2019; Scelzo et al., 2019). As such, this group offers unique opportunities to investigate the similarities and differences in SC biology between vertebrate and invertebrate chordates.
Among colonial tunicates, the ascidian Botryllus schlosseri is recognized as a model species for the study of SCs and asexual reproduction (Manni et al., 2007; Voskoboynik and Weissman, 2015; Manni et al., 2019). This species reproduces both sexually and asexually. In a colony, each individual (zooid) is hermaphroditic and produces swimming larvae via sexual reproduction. Larvae metamorphose into sessile individuals that reproduce asexually by budding (blastogenesis), forming clonal zooids grouped in star-shaped systems and embedded in a common tunic (Figures 1A,B). In each colony, three generations coexist: adult zooids, their primary buds (1Bs), and secondary buds (2Bs). In laboratory conditions (18°C), each week all adult zooids are resorbed during a process called takeover and replaced by their 1Bs; simultaneously, 2Bs become 1Bs and produce a new generation of 2Bs. Each new 2B originates from a thickening of the peribranchial epithelium of a 1B, which further develops into a vesicle undergoing organogenesis. Organogenesis nears completion when the 2B becomes a 1B where tissue differentiation occurs. B. schlosseri SCs are able to freely circulate in the colony haemocoel where, over generations, they contribute to tissue differentiation in buds (Kowarsky et al., 2021). Candidate SCs (cSCs) were identified in B. schlosseri as hemoblast, i.e., round cells, with a diameter of about 5 μm and a high nucleus-cytoplasm ratio (Rosental et al., 2018; Kowarsky et al., 2021). Two cSC transient niches have also been proposed in adult zooids (Figures 1E–H): the haemocoel, ventral to the anterior endostyle (the gland on the pharynx floor which produces the mucus net for filtration) (Voskoboynik et al., 2008; Rosental et al., 2018), and the “cell islands'', i.e., groups of haemocytes in the ventral body wall (Rinkevich et al., 2013). Cell islands are numerous in adult individuals belonging to early-cycle colonies and decrease in number toward the late-cycle. cSCs reside and proliferate in these transient niches and migrate before takeover through the haemocoel to relocate in newly formed niches of the following generations (Rinkevich et al., 2013).
FIGURE 1
The four transcription factors c-Myc, Sox2, Pou5f1 (also named Oct3/4) and Klf4 are of interest given their role in tissue differentiation and in both regulating and inducing pluripotency in vertebrates. Sox2 and Pou5f1 are involved in maintaining pluripotency of embryonic stem cells (Niwa et al., 2000; Avilion et al., 2003; Boyer et al., 2005), while c-Myc and Klf4 play a key role in self-renewal and maintenance of the undifferentiated state (Cartwright et al., 2005; Li et al., 2005). Remarkably, the forced co-expression of these four transcription factors, also named Yamanaka Factors (YFs), is sufficient to reprogram murine and human fibroblasts into induced pluripotent stem (iPS) cells (Takahashi and Yamanaka, 2006; Takahashi et al., 2007).
c-Myc belongs to the MYC gene family, characterized by a C-terminal basic-helix-loop-helix-leucine zipper (bHLHZ) domain mediating protein dimerization and DNA binding, and two N-terminal trans-activation motifs (Atchley and Fitch, 1995; James and Eisenman, 2002). Vertebrates possess multiple Myc genes (c-Myc, N-Myc, L-Myc), with c-Myc involved in cell division and growth control, and considered a central hub for pluripotent SC regulatory networks (Chappell and Dalton, 2013; Fagnocchi and Zippo, 2017). In contrast, two genes are found in a cnidarian with high regenerative ability (Hartl et al., 2010), and only a single homolog is present in Drosophila and most other invertebrates (Gallant, 2006; McFerrin and Atchley, 2011). Despite differences in gene number, the expression of the single Myc gene in the undifferentiated cells of invertebrates suggests that its function is similar to that of c-Myc in vertebrates (Kawamura et al., 2008; Hartl et al., 2010; Fujiwara et al., 2011).
Sox2 is a member of the ancient SOX gene family that has undergone frequent gene duplications and is therefore divided into ten subgroups (A-J) based on the similarity of the DNA-binding domain (Wegner, 1999; Bowles et al., 2000). In vertebrates, Sox2, Sox1 and Sox3 comprises the SoxB1 group, while the close SoxB2 group consists of Sox14 and Sox21. Invertebrates have only a single member each of the SoxB1 and SoxB2 groups (Bowles et al., 2000; Heenan et al., 2016). The role of SoxB1 in pluripotency appears to be conserved among invertebrates (Onal et al., 2012; Niwa et al., 2016), where it is fundamental for neural development and regeneration (Mizuseki et al., 1998; Buescher et al., 2002; Bylund et al., 2003; Graham et al., 2003; Miya and Nishida, 2003; Meulemans and Bronner-Fraser, 2007; Marlow et al., 2009; Ross et al., 2018) and is expressed in multiple adult self-renewing epithelia (Arnold et al., 2011).
Pou5f1 belongs to the metazoan POU gene family, characterized by a modular DNA-binding domain (POU) consisting of a N-terminal POU-specific domain (POUs) and a C-terminal homeodomain (POUh) separated by a linker of variable length. Based upon similarity across the entire POU domain, Pou genes are grouped into six classes (Rosenfeld, 1991), that have evolved through multiple lineage-specific events of gene duplications and losses (Frankenberg and Renfree, 2013; Gold et al., 2014; Onichtchouk, 2016). Pou5f1 belongs to the vertebrate-specific POU5 class (also including Pou5f3 and the primate/rodent-specific Pou5f2), hypothesised to be arisen from a Pou3-like ancestor together with the bilaterian-specific POU2 class (Gold et al., 2014).
Klf4 is a member of the Kruppel like factors (Klf) family, characterized by the presence of three C-terminal highly conserved zinc finger motifs and a high diversification of the associated transactivation/repression domains (Vangapandu and Ai, 2009; Presnell et al., 2015). Depending on their structure and/or function, the Klf members are partitioned into different groups whose number varies across studies (McConnell and Yang, 2010; Pei and Grishin, 2015; Presnell et al., 2015; Pollak et al., 2018). Pei and Grishin (2015) found a minimum of 18 Klf genes in vertebrates and classified them into seven groups (A-G), with invertebrate deuterostomes possessing a single member of each group and the loss of KlfE in the ascidian Ciona intestinalis. Together with Klf1/2/17, the human Klf4 belongs to the group KlfA.
In this study we investigate the presence and the evolution of the YF genes in ascidians, and analyze when and where the orthologs of the vertebrate YF genes are expressed in the blastogenetic cycle of B. schlosseri by using in situ hybridization (ISH) experiments. Our results show patterns of expression consistent with the reported function of these genes in invertebrate and vertebrate development and differentiation. Moreover, the expression of Myc, SoxB1 and Pou3 is retrieved in cSCs residing in cSC niches, suggesting these transcription factors may play a role in the maintenance of undifferentiated state, and highlighting the necessity to better investigate whether similar molecular mechanisms control pluripotency in vertebrates and tunicates.
Methods
Identification and Evolution of Yamanaka Factors in Chordates
The in silico identification of ascidian YF genes and the evolutionary analyses of the YF gene families in chordates were performed on the species listed in Supplementary Table 1, as detailed in Supplementary Data Sheet 1.
B. schlosseri Colonies Rearing
B. schlosseri colonies were collected in the Venice Lagoon and reared at the Department of Biology (University of Padova) in standard laboratory condition (Sabbadin, 1955; Sabbadin, 1960). Colonies were observed daily and fixed at 6 developmental phases as defined by the staging method developed by Sabbadin (1955) (Figure 1A).
RNA Extraction, Yamanaka Factors cDNA Sequencing, and in Situ Hybridization
In order to validate the B. schlosseri predictions of SoxB1, Myc, Pou2 and Pou3, and to develop probes for ISH experiments, total RNA was extracted from B. schlosseri single colonies according to the protocol of Campagna et al. (2016), followed by cDNA synthesis, cloning and sequencing. Whole-mount ISHs were carried out as described in Franchi and Ballarin (2014) using at least three samples per stage. Detailed protocols are in Supplementary Data Sheet 1.
Histology and Transmission Electron Microscopy
Samples have been treated according to the classical protocols for histology (Hematoxylin-Eosin) and Transmission Electron Microscopy (TEM); details on the protocols are reported in Supplementary Data Sheet 1.
Results
Annotation of the YF Gene Families in Ascidians
A total of 47 coding sequences (CDS) of the YF families were re-annotated or annotated ex novo in the reference genomes of the 11 analysed ascidians (see genes with “_p” suffix in Supplementary Tables 2–5, Fasta sequences in Supplementary Data Sheet 2–5, and details in Supplementary Data Sheet 1). Original misannotations, causing partial or anomalous CDS identification, were mainly due to: failure in exon identification, erroneous exon boundaries, incorrect assignment of exons to genes, assembly anomalies, or sequencing errors.
A Single Ancestral Myc Gene Is Expressed Extensively During Bud Development
A single Myc gene was present in B. schlosseri and in all the other ascidians and invertebrate deuterostomes analyzed, and our evolutionary tree clearly shows that this single gene gave rise to the three vertebrate paralogs (L-Myc, N-Myc and C-Myc) through two duplication rounds (Figure 2A). Noteworthy, the evolutionary history of each paralog corresponds exactly with the species tree. Indeed the ascidian sequences cluster according to taxonomic classification in orders and families with high statistical support (Supplementary Table 1 and bootstrap values in Figure 2A).
FIGURE 2
The Myc gene structure is highly conserved in ascidians and consists of two coding exons separated by a very long intron likely responsible for the incorrect Myc annotation in all ascidians but Ciona robusta (range: 10.1–23.8 kb in all analysed ascidians except for the 3–4.6 kb of Molgulidae; see Supplementary Table 2 and Supplementary Data Sheet 1). Remarkably, in all ascidians the second exon encodes for both the bHLHZ domain and a long C-terminal region (mean length and standard deviation: 118 ± 31 amino acids), which are absent in all other analyzed vertebrates and invertebrates genes.
In the blastogenetic cycle, ISHs demonstrated that Myc was expressed in all the tissues across all developmental stages except in the nervous system rudiment both in 2Bs, from the early stages of development, and in 1Bs (Figures 1C,D, 2). Myc expression was weak as buds progressed in development, disappearing in adult tissues, with the exception of cells showing a morphology comparable to cSCs (round cells, with a diameter of 4.92 ± 0.85 µm and a high nucleus-cytoplasm ratio) (Figures 2H,I, 3A,B, Supplementary Data Sheet 6), aggregated in the ventral region in the proximity of the endostyle, and in cell islands in all the phases. These cells were the only labelled hemocytes (Ballarin and Cima, 2005): immunocytes (Figures 2C,G,H) exhibited non-specific labelling, whereas phagocytes (Figure 2H) and vacuolated cells (Figures 2C,D) were never labelled.
FIGURE 3
SoxB1 Is Amply Expressed During Bud Development but Persists in Adult Zooids Only in Candidate Stem Cells
In order to unambiguously identify the ascidian orthologs of the human Sox2 gene, we reconstructed the evolutionary history of both the SoxB1 and SoxB2 groups in deuterostomes. Our tree clearly delineates the two groups (Figure 3C) and shows a general agreement between the gene tree and the ascidian species tree (see taxonomy classification in Supplementary Table 1), although not all nodes have statistically significant bootstrap values (Figure 3C). Both SoxB1 and SoxB2 have conserved gene structure in ascidians (with 2 and 4 coding exons, respectively), and were re-annotated in B. schlosseri and in several other ascidians (Supplementary Table 3 and Supplementary Data Sheet 1).
SoxB1 expression was absent in early stages of blastogenesis until the late 2B stages (stages 4-5–6). Here, SoxB1 transcript was found in the developing branchial epithelium, oesophagus, heart and neural complex (Figures 1C,D, 3). In 1Bs, the gene was additionally expressed in stigmata and endostyle. In adult zooids in early-, mid-, and late-cycle, the expression was restricted to the stigmata and to a few cSCs (round cells, with a diameter of 5.3 ± 1.06 µm and a high nucleus-cytoplasm ratio; Figure 3A; Supplementary Data Sheet 6) aggregated in the region ventral to the endostyle and in cell islands. Among hemocytes, these cells were the only positive ones during the blastogenetic cycle (Figure 3L). However, cSCs expressing SoxB1 were less numerous than those expressing the other YFs, probably in relationship with their functional state (Figure 3B, Supplementary Data Sheet 6).
Pou2 Expression Is Female Gonad-Restricted Whereas Pou3 has a Wider Expression Pattern
Due to its complexity, the evolutionary history of the POU family was reconstructed including the POU2, POU3, POU4 and POU5 classes, with vertebrate representatives of the POU6 class as outgroups (Figure 4 and Supplementary Image 1). In the reconstructed tree all basal nodes remain unresolved and a polytomous POU3 “bush” can be observed (Figure 4A). Increasing the number of replicates does not improve the support of these nodes, indeed bootstopping tests did not converge even after 1,000 trees. However, five highly supported monophyletic clades are identified (see nodes with dots in Figure 4A), corresponding to: 1) the POU2 class; 2) the POU4 class; 3) the vertebrate-specific POU5 class; 4) all analysed vertebrate Pou3f1 genes; and 5) four fast-evolving ascidian sequences (POU-X clade in Figure 4A). Thus our tree confirms the absence of POU5 in ascidians and allows the unambiguous identification of the ascidian orthologs of POU2 and POU4 classes, remarkably present in all 11 analysed ascidian genomes (Figure 4A). Moreover, the topology of both the POU2 and POU4 ascidian clades corresponds to the expected species tree based on the ascidian classification, and a nearly conserved intron/exon structure is observed within each of these clades (see column “Ascidian Classification” and “N° of coding exons” in Supplementary Table 4).
FIGURE 4
Since a published phylogenetic reconstruction performed on all metazoans and all POU classes hypothesized a Pou3-like gene as the last common ancestor of the POU5/3/2 polytomous clade (Gold et al., 2014), we focused our expression pattern analyses on Pou2 and Pou3 (see Discussion).
Two Pou2 genes, named Pou2a and Pou2b, have been identified in the B. schlosseri genome (Bs_g12092_Pou2a and Bs_g13218_Pou2b in Figure 4A). They differ only for the first and last exon while the remaining intron/exon structure is identical (10 exons, with 16 over a total of 18 intron/exon boundaries falling in exactly the same protein position). Thus, the encoded proteins are nearly identical (97% amino acid identity over 255 alignable sites) but have short and completely different N- and C-terminal ends, and different total lengths (359 for Pou2a vs. 293 amino acids for Pou2b). Since the C-terminal region of Pou2a is more similar than Pou2b to the homologous regions of the closest species B. leachii, only Pou2a was experimentally validated (Supplementary Data Sheet 1) and used in ISH experiments. Indeed, we hypothesize that Pou2b is the result of errors in genome assembly as the entire mapping regions and even the intron sequences of Pou2a and Pou2b are almost identical (94% nucleotide identity with 2% gap over 4.6 kb). However, a very recent gene duplication cannot be ruled out.
The POU-X clade only includes representatives of the families Styelidae and Pyuridae and was putatively assigned to the POU3 class based on two considerations: the single-exon gene structure, typical of all other Pou3 of deuterostomes here analysed (Supplementary Table 4), and the presence of elements known to impair the tree resolution (i.e., the very long branches of POU-X sequences associated with the intrinsic polytomy of the POU3 class). Therefore, the Bs_g52290 sequence of POU-X clade (Figure 4A) was considered a valid B. schlosseri Pou3 and used in subsequent ISH experiments.
Pou2 was found expressed in the female germline and, more specifically, in oocytes in previtellogenesis in 1Bs, 2Bs, adult zooids, and circulating in the haemocoel (Figures 1C,D, 4A–D). Moreover, its transcript was retrieved in cells of early developing 2B. Pou3 was found expressed in early-cycle of blastogenic development (Figures 1C,D, 4E–I). Subsequently its expression was restricted to the germline and the rudiments of the neural complex, heart, and digestive system. In the following phases of bud development, the expression persisted at high levels only in the female and male germlines, while lower levels were found in the nervous system. In adult zooids, the presence of Pou3 transcript persisted in germline and in cSCs (round cells, with a diameter of 5.2 ± 1.08 µm and a high nucleus-cytoplasm ratio; Figure 3A, Supplementary Data Sheet 6) found near the endostyle and in cell islands in all the phases. Moreover, it could be retrieved in the endostyle, with an antero-posterior pattern of respectively high and low expression levels. As reported for Myc and SoxB1, cSCs were the only hemocytes expressing Pou3 during the blastogenetic cycle (Figure 4H).
The Klf Gene Family
The KLF proteins are characterized by a high N-terminal sequence variability that makes their alignment unreliable except in the three zinc finger motifs. The reconstructed evolutionary tree, based on the 93 reliably alignable sites of the zinc finger domain and including the KlfA/B/C/D/F/G groups, is almost completely unresolved (Supplementary Image 2) and shows a mildly supported KlfA clade (75% bootstrap) inside a wide polytomous KlfA/B/C/D clade. Only a partial KlfA sequence, corresponding only to the Zn finger domain, was found and annotated by us in an earlier assembly version of the B. schlosseri genome (see “Notes” in Supplementary Table 5), suggesting that the KlfA absence in the current genome version was due to assembly problems. The very high conservation among groups of the available KlfA region (i.e., the Zn finger domain) prevented us from designing target-specific ISH probes able to discriminate between the different Klf groups, therefore no gene expression analyses were performed on this gene.
Yamanaka Factors Are Individually Expressed in Several Organ Rudiments of the Early Bud and, Except Pou2, in Candidate Adult Stem Cells in Proposed Niches
An overview of the spatio-temporal expression pattern of Myc, SoxB1, Pou2 and Pou3 during the blastogenic cycle is summarised in Figures 1C,D. The endostyle, neural complex, heart, gut, undifferentiated germline, and female germline all expressed three factors during the same developmental stages: Myc, SoxB1 and Pou3. The same factors were expressed during the entire cycle in some cSCs, showing comparable dimensions among samples (Figure 3A), in both the SC niches of adult individuals and in hemocele. Almost 70% of adult cSCs in niches expressed Myc and Pou3 (Figures 3A,B, Supplementary Data Sheet 6), suggesting that these two factors should be coexpressed in at least 20% of cSCs. However, SoxB1 is expressed in a lower percentage of cells (Figures 3A,B, Supplementary Data Sheet 6). Of note is that only the female gonad additionally expressed Pou2. Myc and Pou (both Pou2 and Pou3) genes were expressed in early 2Bs. Myc and SoxB1 were recognised during the development of the branchial and peribranchial epithelia, whereas SoxB1 and Pou3 were present during brain formation. Finally, the testis expressed only Myc, and only during development, whereas the adult stigmata maintain the expression of SoxB1.
Discussion
All YF gene families have a complex evolutionary history characterized by several gene duplications, some of which also occurred in the ancestor of vertebrates. Our evolutionary analyses show that a one-to-many relationship exists between ascidian and human orthologs of some YF genes, which means that the ascidian gene has more than one ortholog in humans. Indeed, the ascidian SoxB1 is the ortholog not only of the human/vertebrate YF Sox2, but also of Sox1 and Sox3. Likewise, the ascidian Myc gene is the ortholog not only of the human/vertebrate YF c-Myc, but also of N-Myc and L-Myc. This implies that the functional role and the expression pattern of these ascidian genes could have been inherited by one or many human/vertebrate orthologs. With that in mind, we have analysed the expression pattern of the YF genes in B. schlosseri, a highly regenerative model species belonging to the taxon evolutionarily closest to vertebrates.
Our results show that Myc is expressed in all developing tissues of buds and in cSCs in niches (Figures 1C,D), suggesting a crucial role in undifferentiated territories. Noteworthy, Myc expression in cSCs in their niches is a new finding, not observed in previous studies on Myc expression in two other colonial ascidians. Indeed, in Botryllus primigenus, Myc was proposed as a marker of undifferentiated, but not proliferative, somatic and germline SCs, as it was expressed by multipotent epithelial cells in buds, undifferentiated coelomic cells, and developing gonads (Kawamura et al., 2008). In Polyandrocarpa misakiensis, Myc expression was observed in the atrial epithelium in the organ-forming region of developing bud, and in fibroblast-like cells, which participate in the organogenesis together with the epithelial cells (Fujiwara et al., 2011). Moreover, the noteworthy identification of an ascidian-specific C-terminal tail, located just downstream the bHLHZ domain, incentivises further investigation on the role and the possible additional functions performed by this protein region.
As summarized in Figure 1C, SoxB1 is widely expressed during B. schlosseri asexual development. It is found in the developing nervous system, as has also been reported for Sox2 and other SoxB1 group genes in many metazoans (Wood and Episkopou, 1999; Penzel et al., 1997; Rex et al., 1997; Meulemans and Bronner-Fraser, 2007; Ferrero et al., 2014; Ross et al., 2018). In the developing nervous system we also detected the expression of Pou3 (Figure 1C), known to be involved in neural system development in vertebrates and other invertebrates, (Candiani et al., 2002; Cole and Arnone, 2009; Wollesen et al., 2014; Garner et al., 2016; Wei et al., 2016; Prünster et al., 2019). The expression, in samples at the same developmental stages, of SoxB1 and Pou3 in the developing nervous system, observed here and in previous studies (Iwafuchi-Doi et al., 2011; Hui et al., 2015), suggests a co-regulation of these two factors during neurogenesis. We also retrieved SoxB1 expression in the developing oesophagus, stomach and proximal intestine, suggesting a role in the differentiation and maintenance of these epithelia. This result corresponds with the Sox2 expression in the adult stomach tissues of mammals (Arnold et al., 2011). Moreover, SoxB1 was found expressed in the endostyle of 1B, and in cSCs of the endostyle niche and cell islands, indicating potential involvement in stemness. Supporting this hypothesis, Niwa et al. (2016) demonstrated that C. robusta and other invertebrate SoxB1 proteins are able to replace the pluripotency-associated function of mouse Sox2. Thus the role of SoxB1 in the control of pluripotency in B. schlosseri is plausible.
Our evolutionary analyses demonstrate the absence of the POU5 class in ascidians, indicating that POU5 arose after the tunicate-vertebrate divergence. Since there is no true ortholog of Pou5f1 to be investigated in B. schlosseri, the ancestor of the POU5 class could instead be investigated in its place. Indeed, it remains unclear whether the role of vertebrate Pou5f1 in SC pluripotency was inherited from the ancestral paralog (see discussion in Gold et al., 2014). Our POU evolutionary tree shows a basal polytomy that prevents identifying the POU5 ancestor (Figure 4A). However, our tree is in agreement with a previous wide evolutionary analysis of the entire POU family (including POU1) in all Metazoa (Gold et al., 2014). In that study, a polytomous clade including POU5, POU3 and POU2 classes was identified, and the POU4 class was found as the paralog that diverged after the POU6+POU1 basal clade. Relying on the ancestral state reconstruction and on the taxonomic distribution of POU2 and POU5, Gold et al. (2014) hypothesized that a Pou3-like gene was the last common ancestor of the POU5/3/2 polytomous cluster, in which POU2 and POU5 classes represent bilaterian-specific and vertebrate-specific duplications, respectively. On this basis, and given the goal of our study, we investigated the expression pattern in B. schlosseri of both Pou3, as representative of the Pou3-like ancestor of Pou5f1, and Pou2, as the other closest paralog of Pou5f1.
Remarkably, while POU2 and POU4 are present in all analysed ascidians, Pou3 was found only in ascidians of two Stolidobranchia families (i.e., Styelidae and Pyuridae; Supplementary Table 1 and Figure 4A) and seems to be absent even in the two widely studied Ciona species with highly curated genomic sequence annotations. This leads us to hypothesize a Pou3 loss in most ascidian groups, although the long branch/fast evolution of the ascidian Pou3 clade (Figure 4A) could also explain the failure to identify these genes. As summarized in Figure 1D, we found that B. schlosseri Pou3 is expressed in the endostyle, even in the adult phase, especially in the anterior-most part, close to the previously identified SC niche (Rinkevich et al., 2013). Moreover, Pou3 is expressed in the few pluripotent cells of the early 2B and cSCs in their niches. This expression pattern highlights a likely role of B. schlosseri Pou3 in undifferentiated cells, as also recently evidenced in Botrylloides diegensis/leachii (Viard et al., 2019; Kassmer et al., 2020). Overall, our and previous results support a role in SC pluripotency for the Pou3-like ancestor that gave rise to the vertebrate Pou5f1, here represented by the Pou3 of ascidians. Finally, Pou3 was found expressed in the differentiating heart, gut and gonads, and in the developing nervous system (see above).
As for Pou2, it is expressed exclusively in the immature previtellogenic oocytes, both in gonads and circulating in the colonial vascular system. Remarkably, this expression pattern significantly differs from those of vertebrates and invertebrates, where Pou2 is expressed in many more territories. For example, Pou2 is involved in neurogenesis in cephalopods (together with Pou3) (Wollesen et al., 2014), in the establishment of the dorsal/ventral axis in echinoderms (Range and Lepage, 2011) and in immunity in vertebrates (Strubin et al., 1995; Karnowski et al., 2012). Notably, this is the first report of Pou2 expression in tunicates.
Overall, these results point to an expression pattern of the vertebrate Pou5f1 more similar to that of B. schlosseri Pou3 rather than Pou2, thus supporting the Gold’s and coll. (2014) hypothesis of a Pou3-like gene as ancestor of Pou5f1.
To summarize, we identified in the model species B. schlosseri the orthologs of vertebrate YFs through deep evolutionary analyses. Then, we found that Myc, SoxB1 and Pou3 are individually expressed in several differentiating territories, consistently with the reported function of these genes in chordate embryogenesis. Moreover, we found their expression in cells with features of cSCs exhibiting the same morphological and morphometric features. These cells are located close to the endostyle and in cell islands, and according to our analyses, should coexpress at least Pou3 and Myc. Further studies, such as a deep analysis of niche dynamic and double in situ hybridization experiments, are needed to understand the actual co-expression of these factors in cSCs, which could in turn suggest a shared ancestral pluripotency program in chordates.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository and accession numbers can be found below: GenBank, OL828248-OL828252.
Author contributions
Conception and design: LM, CG; Mariculture, observation and sample collection: FG, MM, VV; Gene annotation and evolutionary analyses: CG, MS, CA; Cloning and sequencing: MS, MM, FG; ISH: MM, VV, MS, FG, LM; Manuscript writing: CG, MS, LM, VV. All authors participated in the data discussion and manuscript finalization.
Funding
This study was supported by MIUR, Italy (PRIN-2015, grant number: 2015NSFHXF) to LM and CG; by the Molecular Biodiversity Laboratory of the Italian node of Lifewatch (CNR, Consiglio Nazionale delle Ricerche) to CG; by the University of Padova, grant BIRD213252 to LM.
Acknowledgments
The authors thank Ayelet Voskoboynik and Mark Kowarsky for providing contigs of an earlier assembly of the B. schlosseri genome and for discussion on gene annotation; Timothy Frawley for critical review and edits; and a reviewer for valuable suggestions.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2022.782722/full#supplementary-material
References
1
ArnoldK.SarkarA.YramM. A.PoloJ. M.BronsonR.SenguptaS.et al (2011). Sox2+ Adult Stem and Progenitor Cells Are Important for Tissue Regeneration and Survival of Mice. Cell Stem Cell9, 317–329. 10.1016/j.stem.2011.09.001
2
AtchleyW. R.FitchW. M. (1995). Myc and Max: Molecular Evolution of a Family of Proto-Oncogene Products and Their Dimerization Partner. Proc. Natl. Acad. Sci.92, 10217–10221. 10.1073/pnas.92.22.10217
3
AvilionA. A.NicolisS. K.PevnyL. H.PerezL.VivianN.Lovell-BadgeR. (2003). Multipotent Cell Lineages in Early Mouse Development Depend on SOX2 Function. Genes Dev.17, 126–140. 10.1101/gad.224503
4
BallarinL.CimaF. (2005). Cytochemical Properties of Botryllus Schlosseri Haemocytes: Indications for Morpho-Functional Characterisation. Eur. J. Histochem.49, 255–264.
5
BowlesJ.SchepersG.KoopmanP. (2000). Phylogeny of the SOX Family of Developmental Transcription Factors Based on Sequence and Structural Indicators. Dev. Biol.227, 239–255. 10.1006/dbio.2000.9883
6
BoyerL. A.LeeT. I.ColeM. F.JohnstoneS. E.LevineS. S.ZuckerJ. P.et al (2005). Core Transcriptional Regulatory Circuitry in Human Embryonic Stem Cells. Cell122, 947–956. 10.1016/j.cell.2005.08.020
7
BrownF. D.TiozzoS.RouxM. M.IshizukaK.SwallaB. J.De TomasoA. W. (2009). Early Lineage Specification of Long-Lived Germline Precursors in the Colonial Ascidian Botryllus Schlosseri. Development136, 3485–3494. 10.1242/dev.037754
8
BrozovicM.DantecC.DardaillonJ.DaugaD.FaureE.GinesteM.et al (2018). ANISEED 2017: Extending the Integrated Ascidian Database to the Exploration and Evolutionary Comparison of Genome-Scale Datasets. Nucleic Acids Res.46, D718–D725. 10.1093/nar/gkx1108
9
BuescherM.HingF. S.ChiaW. (2002). Formation of Neuroblasts in the Embryonic central Nervous System of Drosophila melanogaster Is Controlled by SoxNeuro. Development129, 4193–4203. 10.1242/dev.129.18.4193
10
BylundM.AnderssonE.NovitchB. G.MuhrJ. (2003). Vertebrate Neurogenesis Is Counteracted by Sox1-3 Activity. Nat. Neurosci.6, 1162–1168. 10.1038/nn1131
11
CampagnaD.GaspariniF.FranchiN.VituloN.BallinF.ManniL.et al (2016). Transcriptome Dynamics in the Asexual Cycle of the Chordate Botryllus Schlosseri. BMC Genomics17, 275. 10.1186/s12864-016-2598-1
12
CandianiS.CastagnolaP.OliveriD.PestarinoM. (2002). Cloning and Developmental Expression of AmphiBrn1/2/4, a POU III Gene in Amphioxus. Mech. Dev.116, 231–234. 10.1016/s0925-4773(02)00146-6
13
CartwrightP.McLeanC.SheppardA.RivettD.JonesK.DaltonS. (2005). LIF/STAT3 Controls ES Cell Self-Renewal and Pluripotency by a Myc-dependent Mechanism. Development132, 885–896. 10.1242/dev.01670
14
ChappellJ.DaltonS. (2013). Roles for MYC in the Establishment and Maintenance of Pluripotency. Cold Spring Harb. Perspect. Med.3, a014381. 10.1101/cshperspect.a014381
15
ColeA. G.ArnoneM. I. (2009). Fluorescent In Situ Hybridization Reveals Multiple Expression Domains for SpBrn1/2/4 and Identifies a Unique Ectodermal Cell Type that Co-expresses the ParaHox Gene SpLox. Gene Expr. Patterns9, 324–328. 10.1016/j.gep.2009.02.005
16
DelsucF.BrinkmannH.ChourroutD.PhilippeH. (2006). Tunicates and Not Cephalochordates Are the Closest Living Relatives of Vertebrates. Nature439, 965–968. 10.1038/nature04336
17
FagnocchiL.ZippoA. (2017). Multiple Roles of MYC in Integrating Regulatory Networks of Pluripotent Stem Cells. Front. Cel Dev. Biol.5, 7. 10.3389/fcell.2017.00007
18
FerreroE.FischerB.RussellS. (2014). SoxNeuro Orchestrates central Nervous System Specification and Differentiation in Drosophila and Is Only Partially Redundant with Dichaete. Genome Biol.15, R74. 10.1186/gb-2014-15-5-r74
19
FranchiN.BallarinL. (2014). Preliminary Characterization of Complement in a Colonial Tunicate: C3, Bf and Inhibition of C3 Opsonic Activity by Compstatin. Dev. Comp. Immunol.46 (2), 430–438. 10.1016/j.dci.2014.05.014
20
FrankenbergS.RenfreeM. B. (2013). On the Origin of POU5F1. BMC Biol.11, 56. 10.1186/1741-7007-11-56
21
FujiwaraS.IsozakiT.MoriK.KawamuraK. (2011). Expression and Function of Myc during Asexual Reproduction of the Budding Ascidian Polyandrocarpa Misakiensis. Dev. Growth Differ.53, 1004–1014. 10.1111/j.1440-169X.2011.01312.x
22
GallantP. (2006). Myc/Max/Mad in Invertebrates: the Evolution of the Max Network. Curr. Top. Microbiol. Immunol.302, 235–253. 10.1007/3-540-32952-8_9
23
GarnerS.ZyskI.ByrneG.KramerM.MollerD.TaylorV.et al (2016). Neurogenesis in Sea Urchin Embryos and the Diversity of Deuterostome Neurogenic Mechanisms. Development143, 286–297. 10.1242/dev.124503
24
GoldD. A.GatesR. D.JacobsD. K. (2014). The Early Expansion and Evolutionary Dynamics of POU Class Genes. Mol. Biol. Evol.31, 3136–3147. 10.1093/molbev/msu243
25
GrahamV.KhudyakovJ.EllisP.PevnyL. (2003). SOX2 Functions to Maintain Neural Progenitor Identity. Neuron39, 749–765. 10.1016/s0896-6273(03)00497-5
26
GaspariniF.ManniL.CimaF.ZanioloG.BurighelP.CaicciF.et al (2015). Sexual and Asexual Reproduction in the Colonial Ascidian Botryllus schlosseri. Genesis53, 105–20. 10.1002/dvg.22802
27
HartlM.MitterstillerA.-M.ValovkaT.BreukerK.HobmayerB.BisterK. (2010). Stem Cell-specific Activation of an Ancestral Myc Protooncogene with Conserved Basic Functions in the Early Metazoan Hydra. Proc. Natl. Acad. Sci. USA107, 4051–4056. 10.1073/pnas.0911060107
28
HeenanP.ZondagL.WilsonM. J. (2016). Evolution of the Sox Gene Family within the Chordate Phylum. Gene575, 385–392. 10.1016/j.gene.2015.09.013
29
HuiS. P.NagT. C.GhoshS. (2015). Characterization of Proliferating Neural Progenitors after Spinal Cord Injury in Adult Zebrafish. PLoS ONE10, e0143595. 10.1371/journal.pone.0143595
30
Iwafuchi-DoiM.YoshidaY.OnichtchoukD.LeichsenringM.DrieverW.TakemotoT.et al (2011). The Pou5f1/Pou3f-dependent but SoxB-independent Regulation of Conserved Enhancer N2 Initiates Sox2 Expression during Epiblast to Neural Plate Stages in Vertebrates. Dev. Biol.352, 354–366. 10.1016/j.ydbio.2010.12.027
31
JamesL.EisenmanR. N. (2002). Myc and Mad bHLHZ Domains Possess Identical DNA-Binding Specificities but Only Partially Overlapping Functions In Vivo. Proc. Natl. Acad. Sci. USA99, 10429–10434. 10.1073/pnas.162369299
32
KassmerS. H.LangenbacherA. D.De TomasoA. W. (2020). Integrin-alpha-6+ Candidate Stem Cells Are Responsible for Whole Body Regeneration in the Invertebrate Chordate Botrylloides Diegensis. Nat. Commun.11, 4435. 10.1038/s41467-020-18288-w
33
KassmerS. H.LangenbacherA.De TomasoA. W. (2019). Primordial Blasts, a Population of Blood Borne Stem Cells Responsible for Whole Body Regeneration in a Basal Chordate. BioRxiv. 10.1101/647578
34
KawamuraK.TachibanaM.SunanagaT. (2008). Cell Proliferation Dynamics of Somatic and Germline Tissues during Zooidal Life Span in the Colonial Tunicate Botryllus Primigenus. Dev. Dyn.237, 1812–1825. 10.1002/dvdy.21592
35
KowarskyM.AnselmiC.HottaK.BurighelP.ZanioloG.CaicciF.et al (2021). Sexual and Asexual Development: Two Distinct Programs Producing the Same Tunicate. Cell Rep34, 108681. 10.1016/j.celrep.2020.108681
36
LiF.WangY.ZellerK. I.PotterJ. J.WonseyD. R.O’DonnellK. A.et al (2005). Myc Stimulates Nuclearly Encoded Mitochondrial Genes and Mitochondrial Biogenesis. Mol. Cel. Biol.25, 6225–6234. 10.1128/MCB.25.14.6225-6234.2005
37
ManniL.GaspariniF.HottaK.IshizukaK. J.RicciL.TiozzoS.et al (2014). Ontology for the Asexual Development and Anatomy of the Colonial Chordate Botryllus schlosseri. PLoS One9, e96434. 10.1371/journal.pone.0096434
38
ManniL.AnselmiC.CimaF.GaspariniF.VoskoboynikA.MartiniM.et al (2019). Sixty Years of Experimental Studies on the Blastogenesis of the Colonial Tunicate Botryllus Schlosseri. Dev. Biol.448, 293–308. 10.1016/j.ydbio.2018.09.009
39
ManniL.ZanioloG.CimaF.BurighelP.BallarinL. (2007). Botryllus Schlosseri: a Model Ascidian for the Study of Asexual Reproduction. Dev. Dyn.236, 335–352. 10.1002/dvdy.21037
40
MarlowH. Q.SrivastavaM.MatusD. Q.RokhsarD.MartindaleM. Q. (2009). Anatomy and Development of the Nervous System of Nematostella vectensis, an Anthozoan Cnidarian. Dev. Neurobiol.69, 235–254. 10.1002/dneu.20698
41
McConnellB. B.YangV. W. (2010). Mammalian Krüppel-like Factors in Health and Diseases. Physiol. Rev.90, 1337–1381. 10.1152/physrev.00058.2009
42
McFerrinL. G.AtchleyW. R. (2011). Evolution of the Max and Mlx Networks in Animals. Genome Biol. Evol.3, 915–937. 10.1093/gbe/evr082
43
MeulemansD.Bronner-FraserM. (2007). Insights from Amphioxus into the Evolution of Vertebrate Cartilage. PLoS ONE2, e787. 10.1371/journal.pone.0000787
44
MiyaT.NishidaH. (2003). Expression Pattern and Transcriptional Control of SoxB1 in Embryos of the Ascidian Halocynthia Roretzi. Zool. Sci.20, 59–67. 10.2108/zsj.20.59
45
MizusekiK.KishiM.MatsuiM.NakanishiS.SasaiY. (1998). Xenopus Zic-Related-1 and Sox-2, Two Factors Induced by Chordin, Have Distinct Activities in the Initiation of Neural Induction. Development125, 579–587. 10.1242/dev.125.4.579
46
NiwaH.MiyazakiJ.SmithA. G. (2000). Quantitative Expression of Oct-3/4 Defines Differentiation, Dedifferentiation or Self-Renewal of ES Cells. Nat. Genet.24, 372–376. 10.1038/74199
47
NiwaH.NakamuraA.UrataM.Shirae-KurabayashiM.KurakuS.RussellS.et al (2016). The Evolutionally-Conserved Function of Group B1 Sox Family Members Confers the Unique Role of Sox2 in Mouse ES Cells. BMC Evol. Biol.16, 173. 10.1186/s12862-016-0755-4
48
OnalP.GrünD.AdamidiC.RybakA.SolanaJ.MastrobuoniG.et al (2012). Gene Expression of Pluripotency Determinants Is Conserved between Mammalian and Planarian Stem Cells. EMBO J.31, 2755–2769. 10.1038/emboj.2012.110
49
OnichtchoukD. (2016). Evolution and Functions of Oct4 Homologs in Non-mammalian Vertebrates. Biochim. Biophys. Acta1859, 770–779. 10.1016/j.bbagrm.2016.03.013
50
PattengaleN. D.AlipourM.Bininda-EmondsO. R. P.MoretB. M. E.StamatakisA. (2010). How many Bootstrap Replicates Are Necessary?J. Comput. Biol.17, 337–354. 10.1089/cmb.2009.0179
51
PeiJ.GrishinN. V. (2015). C2H2 Zinc finger Proteins of the SP/KLF, Wilms Tumor, EGR, Huckebein, and Klumpfuss Families in Metazoans and beyond. Gene573, 91–99. 10.1016/j.gene.2015.07.031
52
PenzelR.OschwaldR.ChenY.TackeL.GrunzH. (1997). Characterization and Early Embryonic Expression of a Neural Specific Transcription Factor xSOX3 in Xenopus laevis. Int. J. Dev. Biol.41, 667–677.
53
PollakN. M.HoffmanM.GoldbergI. J.DrosatosK. (2018). Krüppel-like Factors: Crippling and Un-crippling Metabolic Pathways. JACC Basic Transl. Sci.3, 132–156. 10.1016/j.jacbts.2017.09.001
54
PresnellJ. S.SchnitzlerC. E.BrowneW. E. (2015). KLF/SP Transcription Factor Family Evolution: Expansion, Diversification, and Innovation in Eukaryotes. Genome Biol. Evol.7, 2289–2309. 10.1093/gbe/evv141
55
PrünsterM. M.RicciL.BrownF. D.TiozzoS. (2019). De Novo neurogenesis in a Budding Chordate: Co-option of Larval Anteroposterior Patterning Genes in a Transitory Neurogenic Organ. Dev. Biol.448, 342–352. 10.1016/j.ydbio.2018.10.009
56
RangeR.LepageT. (2011). Maternal Oct1/2 Is Required for Nodal and Vg1/Univin Expression during Dorsal-Ventral axis Specification in the Sea Urchin Embryo. Dev. Biol.357, 440–449. 10.1016/j.ydbio.2011.07.005
57
RexM.OrmeA.UwanoghoD.TointonK.WigmoreP. M.SharpeP. T.et al (1997). Dynamic Expression of Chicken Sox2 and Sox3 Genes in Ectoderm Induced to Form Neural Tissue. Dev. Dyn.209, 323–332. 10.1002/(SICI)1097-0177(199707)209:3<323::AID-AJA7>3.0.CO;2-K
58
RinkevichB.ShlembergZ.FishelsonL. (1995). Whole-body Protochordate Regeneration from Totipotent Blood Cells. Proc. Natl. Acad. Sci. USA92, 7695–7699. 10.1073/pnas.92.17.7695
59
RinkevichY.VoskoboynikA.RosnerA.RabinowitzC.PazG.OrenM.et al (2013). Repeated, Long-Term Cycling of Putative Stem Cells between Niches in a Basal Chordate. Dev. Cel24, 76–88. 10.1016/j.devcel.2012.11.010
60
RosenfeldM. G. (1991). POU-domain Transcription Factors: Pou-Er-Ful Developmental Regulators. Genes Dev.5, 897–907. 10.1101/gad.5.6.897
61
RosentalB.KowarskyM.SeitaJ.CoreyD. M.IshizukaJ. K.PalmeriK. J.et al (2018). Complex Mammalian-Like Haematopoietic System Found in a Colonial Chordate. Nature564, 425–429. 10.1038/s41586-018-0783-x
62
RossK. G.MolinaroA. M.RomeroC.DockterB.CableK. L.GonzalezK.et al (2018). Soxb1 Activity Regulates Sensory Neuron Regeneration, Maintenance, and Function in Planarians. Dev. Cel47, 331–347. 10.1016/j.devcel.2018.10.014
63
SabbadinA. (1960). Nuove Ricerche Sull’ Inversione Sperimentale Situs Viscerum Botryllus Schlosseri Sabbadin. Estratto dall’Archivio di Oceanografia e Limnologia11.
64
SabbadinA. (1955). Osservazioni sullo sviluppo, l’accrescimento e la riproduzione di Botryllus schlosseri (Pallas), in condizioni di laboratorio. Bolletino di zoologia22, 243–263. 10.1080/11250005509439204
65
ScelzoM.AliéA.PagnottaS.LejeuneC.HenryP.GillettaL.et al (2019). Novel Budding Mode in Polyandrocarpa Zorritensis: a Model for Comparative Studies on Asexual Development and Whole Body Regeneration. Evodevo10, 7. 10.1186/s13227-019-0121-x
66
SunanagaT.SaitoY.KawamuraK. (2006). Postembryonic Epigenesis of Vasa-Positive Germ Cells from Aggregated Hemoblasts in the Colonial Ascidian, Botryllus Primigenus. Dev. Growth Differ.48, 87–100. 10.1111/j.1440-169X.2006.00849.x
67
TakahashiK.TanabeK.OhnukiM.NaritaM.IchisakaT.TomodaK.et al (2007). Induction of Pluripotent Stem Cells from Adult Human Fibroblasts by Defined Factors. Cell131, 861–872. 10.1016/j.cell.2007.11.019
68
TakahashiK.YamanakaS. (2006). Induction of Pluripotent Stem Cells from Mouse Embryonic and Adult Fibroblast Cultures by Defined Factors. Cell126, 663–676. 10.1016/j.cell.2006.07.024
69
VangapanduH. V.AiW. (2009). Kruppel Like Factor 4 (KLF4): A Transcription Factor with Diverse Context-Dependent Functions. Gene Ther. Mol. Biol.13a, 194–204.
70
ViardF.RobyC.TuronX.BouchemousseS.BishopJ. (2019). Cryptic Diversity and Database Errors Challenge Non-indigenous Species Surveys: An Illustration with Botrylloides Spp. In the English Channel and Mediterranean Sea. Front. Mar. Sci.6. 10.3389/fmars.2019.00615
71
VoskoboynikA.Simon-BlecherN.SoenY.RinkevichB.De TomasoA. W.IshizukaK. J.et al (2007). Striving for Normality: Whole Body Regeneration through a Series of Abnormal Generations. FASEB J.21, 1335–1344. 10.1096/fj.06-7337com
72
VoskoboynikA.SoenY.RinkevichY.RosnerA.UenoH.ReshefR.et al (2008). Identification of the Endostyle as a Stem Cell Niche in a Colonial Chordate. Cell Stem Cell3, 456–464. 10.1016/j.stem.2008.07.023
73
VoskoboynikA.WeissmanI. L. (2015). Botryllus Schlosseri, an Emerging Model for the Study of Aging, Stem Cells, and Mechanisms of Regeneration. Invertebr. Reprod. Dev.59, 33–38. 10.1080/07924259.2014.944673
74
WegnerM. (1999). From Head to Toes: the Multiple Facets of Sox Proteins. Nucleic Acids Res.27, 1409–1420. 10.1093/nar/27.6.1409
75
WeiZ.AngererL. M.AngererR. C. (2016). Neurogenic Gene Regulatory Pathways in the Sea Urchin Embryo. Development143, 298–305. 10.1242/dev.125989
76
WollesenT.McDougallC.DegnanB. M.WanningerA. (2014). POU Genes Are Expressed during the Formation of Individual Ganglia of the Cephalopod central Nervous System. Evodevo5, 41. 10.1186/2041-9139-5-41
77
WoodH. B.EpiskopouV. (1999). Comparative Expression of the Mouse Sox1, Sox2 and Sox3 Genes from Pre-gastrulation to Early Somite Stages. Mech. Dev.86, 197–201. 10.1016/S0925-4773(99)00116-1
78
ZhaoA.QinH.FuX. (2016). What Determines the Regenerative Capacity in Animals?Bioscience66, 735–746. 10.1093/biosci/biw079
Summary
Keywords
ascidians, asexual development, Botryllus schlosseri, cMyc, Pou, Klf4, Sox2, Yamanaka factors
Citation
Vanni V, Salonna M, Gasparini F, Martini M, Anselmi C, Gissi C and Manni L (2022) Yamanaka Factors in the Budding Tunicate Botryllus schlosseri Show a Shared Spatio-Temporal Expression Pattern in Chordates. Front. Cell Dev. Biol. 10:782722. doi: 10.3389/fcell.2022.782722
Received
24 September 2021
Accepted
09 February 2022
Published
07 March 2022
Volume
10 - 2022
Edited by
Alberto Stolfi, Georgia Institute of Technology, United States
Reviewed by
Federico David Brown, University of São Paulo, Brazil
Filomena Ristoratore, Zoological Station Anton Dohrn, Italy
Updates
Copyright
© 2022 Vanni, Salonna, Gasparini, Martini, Anselmi, Gissi and Manni.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lucia Manni, lucia.manni@unipd.it; Carmela Gissi, carmela.gissi@uniba.it
†These authors have contributed equally to this work and share first authorship
‡These authors have contributed equally to this work and share last authorship
This article was submitted to Morphogenesis and Patterning, a section of the journal Frontiers in Cell and Developmental Biology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.