Skip to main content
  • Research article
  • Open access
  • Published:

Evolutionary origin of gastrulation: insights from sponge development

Abstract

Background

The evolutionary origin of gastrulation—defined as a morphogenetic event that leads to the establishment of germ layers—remains a vexing question. Central to this debate is the evolutionary relationship between the cell layers of sponges (poriferans) and eumetazoan germ layers. Despite considerable attention, it remains unclear whether sponge cell layers undergo progressive fate determination akin to eumetazoan primary germ layer formation during gastrulation.

Results

Here we show by cell-labelling experiments in the demosponge Amphimedon queenslandica that the cell layers established during embryogenesis have no relationship to the cell layers of the juvenile. In addition, juvenile epithelial cells can transdifferentiate into a range of cell types and move between cell layers. Despite the apparent lack of cell layer and fate determination and stability in this sponge, the transcription factor GATA, a highly conserved eumetazoan endomesodermal marker, is expressed consistently in the inner layer of A. queenslandica larvae and juveniles.

Conclusions

Our results are compatible with sponge cell layers not undergoing progressive fate determination and thus not being homologous to eumetazoan germ layers. Nonetheless, the expression of GATA in the sponge inner cell layer suggests a shared ancestry with the eumetazoan endomesoderm, and that the ancestral role of GATA in specifying internalised cells may antedate the origin of germ layers. Together, these results support germ layers and gastrulation evolving early in eumetazoan evolution from pre-existing developmental programs used for the simple patterning of cells in the first multicellular animals.

Background

It has been long argued and widely accepted that early metazoan evolution included the progressive addition and elaboration of cell layers, with the first animals being little more than a small sphere of epithelial cells. The addition of a second inner layer is thought to have followed, resulting in the so-called ‘gastraea’ that possessed an ectoderm and endoderm. A middle mesodermal layer is proposed to have evolved last [1, 2]. It has been noted that this progression is reflected in contemporary early branching animals, with sponges and placozoans being interpreted as having one, two or no germ layers, cnidarians and ctenophores possessing two or possibly three germ layers and bilaterians having three or possibly four [2, 3]. The generation of two or more germ layers by gastrulation and the consistency of the fate of these layers unites the Eumetazoa (though exceptions exist; for example, vertebrate tail-bud mesoderm [2]). Further, comparative developmental genetic evidence supports the homology of the endomesoderm [4] and ectodermal neurogenesis [5] across eumetazoans. However, it remains controversial whether the cell layers in sponges are homologous to eumetazoan germ layers and whether sponges undergo gastrulation.

Sponges consist of four classes - Demospongiae, Calcarea, Homoscleromorpha and Hexactinellida [6, 7] - of simple animals that share a common body organisation. As juveniles and adults, they are sessile, tri-layered animals that lack true muscles and nerves. Sponges possess an internal network of canals and ciliated choanocyte chambers lined with epithelial cells, primarily endopinacocytes and choanocytes, and are separated from the external environment by another epithelial layer, the exopinacoderm (Figure 1; hexactinellids are constructed as a syncytium). Sandwiched between these epithelial layers is the collagenous mesohyl, which is enriched with multiple types of amoebocytes, including the pluripotent stem cell type, the archeocyte (Figure 1). Water currents created by ciliary beating in choanocyte chambers enable most of the physiological requirements of sponges, including feeding, respiration, excretion and reproduction. This juvenile/adult body plan is the outcome of the dramatic reorganisation of the radially-symmetrical, bi- or trilayered larva at metamorphosis [8].

Figure 1
figure 1

Evolution of germ layers and gastrulation. Animal multicellularity, with the larval and juvenile/adult body plans of extant metazoans possessing multiple cell layers, is depicted above. The arrow represents metamorphosis from larval to juvenile/adult forms. Cell layers are coloured: orange, inner layer; green, outer layer; and blue, middle layer. In the case of cnidarians and bilaterians, these correspond to endoderm, ectoderm and mesoderm (bilaterians only), respectively. The phylogenetic relationship of these clades and the sister group to metazoans, the choanoflagellates, is shown below, along with the evolutionary origin of multicellularity. The origin of germ layers and gastrulation is debatable, occurring either before or after the divergence of sponges and eumetazoans.

Sponges utilise recognisable gastrulation-like morphogenetic movements during embryogenesis (for example, delamination, ingression, egression and invagination) and metamorphosis (for example, epithelial-mesenchyme transitions, EMT) (summarised in [8, 9]). In marine haplosclerid demosponges, for instance, morphological evidence suggests that early cleavage leads to the formation of a stereoblastula, which then delaminates to form a multilayered parenchymella larva (for example, Amphimedon queenslandica[10]). In homoscleromorphs (for example, Oscarella spp.), multipolar egression of a stereoblastula generates a single-layered cinctoblastula larva [11], while in a hexactinillid (for example, Oopsacas minuta), cellular delamination of a coeloblastula gives rise to the trichimella larva [12, 13]. Interestingly, the pattern of embryogenesis can be polymorphic. In a halisarcid demosponge (for example, Halisarca dujardini), multipolar ingression and/or invagination of a coeloblastula or a stereoblastula forms a multilayered parenchymella or disphaerula (‘a hollow sphere within a sphere’) larva [14, 15]. Epithelial-mesenchymal transitions are common at metamorphosis in diverse homoscleromorphs, calcareans and demosponges (reviewed in [8]).

Despite the occurrence of gastrulation-like morphogenetic processes in sponges, it is still unclear how sponge development relates to eumetazoan gastrulation, which is defined as a morphogenetic event that leads to the establishment of germ layers. Some authors consider cell movements that generate multi-layered embryos during sponge embryogenesis gastrulation (for example, [12, 16]), although some argue that gastrulation occurs during sponge metamorphosis when an adult internal layer develops by EMT and transdifferentiation of larval ciliated epithelial cells (for example, [17]). Yet it remains unclear whether sponge cell layers display the defining characteristic of eumetazoan gastrulation. This uncertainty led others to question the existence of germ layers, and hence gastrulation, in sponges (for example, [18]).

Here, we show by cell-labelling experiments in the demosponge A. queenslandica that the cell layers established during embryogenesis lack fixed identities and do not directly give rise to the cell layers of the juvenile. In addition, inner epithelial choanocytes in the juvenile can alter their fate by transdifferentiating into a range of cell types, including cells of the outer epithelium. Thus, sponge cell layers lack fate determination and stability. However, we find that expression of the sponge orthologue to the highly conserved eumetazoan endomesodermal marker GATA[4] is consistently restricted to the inner layer of larvae and juveniles, suggesting a deep shared ancestry between the sponge internal cell layer and the eumetazoan endomesoderm.

Results and discussion

Germ layers by definition exhibit restricted cell fate [2, 3]. Thus, to address whether sponge embryonic and early juvenile cell layers are germ layers, we employed cell-labelling techniques in the demosponge A. queenslandica to examine the fate of larval epithelial and internal cells at metamorphosis, as well as the stability and fate of differentiated choanocytes in newly established juvenile tissues at the so-called ‘rhagon’ stage (staging described in Additional file 1: Figure S1).

The swimming A. queenslandica larva metamorphoses into a feeding juvenile with choanocyte chambers within three days of attaching its anterior end onto a natural substratum (crustose coralline algae) at 25°C. Within hours of settling, the three-layered, bullet-shaped larva transforms into an encrusting mat. A previous cell-labelling study in A. queenslandica indicated that external larval cells transdifferentiated into internal choanocytes in juveniles [10]. However, it was not determined whether all cell types found on the larval exterior are transdifferentiating at metamorphosis, and whether the larval epithelial cells differentiate directly into choanocytes or indirectly via an intermediate cell type. In other sponges, the larval epithelial layer has been reported to shed entirely [19], be phagocytised by archeocytes [2022], differentiate into choanocytes through a non-ciliated amoebocyte intermediate [2325], or directly differentiate into choanocytes without loss of cilia [26] (reviewed in [8]). Using the terminal deoxynucleotidyl transferase dUTP nick end-labelling (TUNEL) assay for detecting apoptotic DNA fragmentation, we demonstrate that a large proportion of A. queenslandica larval epithelial cells undergo programmed cell death shortly after the onset of metamorphosis [see Additional file 2: Figure S2]. TUNEL-labelled nuclear fragments can be observed readily in the metamorphosing animal, often as phagocytised bodies within archeocytes [see Additional file 2: Figure S2].

To study the fate of non-apoptosing larval epithelial cells, the outer layer was labelled with the plasma-membrane dye CM-DiI. We assume that the cell lineage will be indicated by retention of this label in their plasma membrane, while the plasma membrane of the cells that phagocytised apoptotic fragments of larval epithelial origin will not be labelled. An overnight incubation of free-swimming A. queenslandica larvae in seawater containing 10 μM CM-DiI strongly labelled the plasma membrane of ciliated epithelial cell types, the columnar epithelial cell and the flask cell (Figure 2A). During metamorphosis, these labelled cells internalised, resorbing their cilia and differentiating into archeocytes, which migrated throughout the animal [see Additional file 3: Figure S3]. Within two to three days, a number of juvenile cell types were labelled, including choanocytes and exopinacocytes (Figure 2B, C), suggesting that these differentiated from pluripotent archeocytes, which were originally derived from the larval epithelium. Therefore, at metamorphosis ciliated epidermal cells of the larvae can either undergo apoptosis or transdifferentiate, via an archeocyte intermediate, into a range of juvenile cell types that are found in the inner, outer or middle layer of the juvenile.

Figure 2
figure 2

Cell labelling and lineage-tracing through metamorphosis in Amphimedon queenslandica. A, D: Distributions of labelled cells in free-swimming larvae. Anterior (swimming axis) is to the bottom. B, C, E, F: Descendants of the labelled cells in juveniles. Nuclei are stained with DAPI. In E and F, the juveniles are labelled with an antibody against tyrosinated tubulin (tyrTub). A: Longitudinal confocal sections through the centre of a free-swimming larva incubated with CM-DiI, showing strong labelling in ciliated epidermal cell types, the columnar epithelial cell (co) and the flask cell (fc) (inset), with little labelling in inner cell mass (icm). B: Choanocytes in chambers (ch). In some cases, a subset of choanocytes is CM-DiI-labelled in a single choanocyte chamber (arrowhead in inset), suggesting that multiple precursor cells can be involved in development of a single chamber. C: Labelled exopinacocytes (arrowheads). D: A confocal longitudinal section through the centre of a free-swimming larva pulse-labelled with EdU. Note that the labelled cells localise in the inner cell mass (icm) and are likely to be proliferating archeocytes with characteristic large nucleoli (nu) [see Additional file 3: Figure S3]. E: An EdU-positive choanocyte in the chamber (arrowhead). F: An EdU-positive exopinacocyte (arrowhead). Other abbreviations: ep, outer layer epithelium; ci, cilium; ex external environment; mh, mesohyl. Scale bars: 100 μm (A, D), 10 μm (B, C, E, F, inset in D), 5 μm (inset in A). DAPI, 4',6-diamidino-2-phenylindole; EdU, 5-ethynyl-2′-deoxyuridine.

Using the thymidine analog, 5-ethynyl-2′-deoxyuridine (EdU), to label the nuclei of proliferating cells, we found that only amoeboid cells localised to the inner cell mass, identified here as archeocytes [see Additional file 4: Figure S4], and no other cell types, continue to proliferate in the swimming larva (Figure 2D). Taking advantage of this observation we followed the fate of these EdU-labelled cells during metamorphosis. We found that these archeocytes of the larval inner cell mass origin also gave rise to choanocytes and exopinacocytes during metamorphosis (Figure 2E, F). Together, these observations indicate that a range of larval cell types can generate the diversity of cell types present in the juvenile and that there is no correspondence between the cell layers established during embryogenesis and those produced at metamorphosis.

Sponge metamorphosis has been regarded by some as gastrulation (summarised in [8, 17]). As observed here and in previous studies, morphogenetic processes operational during sponge metamorphosis occur also during eumetazoan gastrulation (for example, EMT; reviewed in [8, 27]). Therefore, we sought to determine if cell layers established at the end of A. queenslandica metamorphosis are developmentally determined. To test this, we labelled choanocytes that form internal choanocyte chambers and a small number of archeocytes surrounding choanocyte chambers by incubating juveniles with CM-DiI (Figure 3A, B). Within one week, the majority of the labelled choanocytes had undergone EMT and dedifferentiated into archeocytes (Figure 3C, D; Additional file 5: Figure S5), of which some had migrated away from the chambers and differentiated into other cell types including exopinacocytes (Figure 3E) and sclerocytes (Figure 3F). Thus, inner epithelial choanocytes are capable of transdifferentiating into outer epithelial cells via an archeocyte intermediate, demonstrating the lack of cell fate determination.

Figure 3
figure 3

Choanocyte labelling and lineage-tracing in juvenile Amphimedon queenslandica. A, B: Choanocytes in chambers (ch) in juveniles are labelled with CM-DiI; note strong plasma membrane labelling in choanocytes (inset). C-F: Descendants of juveniles whose chambers were labelled with CM-DiI seven days earlier (as in A). Multiple other juvenile cell types are now labelled and spread through the sponge body. Note in (D) the limited number of CM-DiI-labelled choanocytes remaining in the choanocyte chambers (ch) seven days after labelling and the presence of labelled archeocytes with large nuclei (nu; inset), consistent with transdifferentiation of choanocytes into archeocytes. Note the presence of nuclear fragments (fr) presumably resulting from phagocytosis that has occurred during de-differentiation [see Additional file 5: Figure S5]. Arrowheads in E and F show a labelled exopinacocyte and a sclerocyte, respectively, indicating these cell types can also be produced from de-differentiating choanocytes. Nuclei are stained with DAPI. Scale bars: 100 μm (A, C), 10 μm (B, D-F). DAPI, 4',6-diamidino-2-phenylindole.

The lack of fate determination in cell layers in A. queenslandica is consistent with this demosponge lacking germ layers as classically defined. In light of this finding, genes shown to be instrumental in eumetazoan ectoderm or endomesoderm formation may not be expected to be present in demosponges. Indeed, of the genes previously implicated in the differentiation of the endomesoderm across Cnidaria and Bilateria, namely twist, snail, forkhead and GATA[4], only GATA was identified in the genome of the demosponge A. queenslandica[28] [see Additional file 6: Figure S6]. Also, amongst the 51 transcription factor- and signalling molecule-encoding genes shown to be ‘endomesodermally’ expressed in the cnidarian Nematostella vectensis[29], only nine orthologues are present in the A. queenslandica genome [see Additional file 7: Table S1].

A. queenslandica GATA transcripts are localised to a substantial fraction of internal cells of late-stage embryos, free-swimming larvae (Figure 4B-D) and juveniles (Figure 4F-I), but not detectable in outer layer epithelial cells at any developmental stage, similar to predominantly internal (that is, endomesodermal) expression of bilaterian GATA456[3032] and cnidarian GATA[4]. Thus, although this study argues against the existence of true germ layers in demosponges, the striking parallel in GATA expression patterns suggests a common ancestry between the sponge inner cell layer and the eumetazoan endomesoderm.

Figure 4
figure 4

GATA mRNA expression during development in A. queenslandica . Late blastulae (A; sensu Leys and Degnan [10]), mid-stage embryos with a pigment spot (ps) (B), late-stage embryos with a pigment ring (pr) (C), a free-swimming larva (D), a settlement stage postlarva (E), and a juvenile (F-I) were labelled with the antisense riboprobe AqGATA. The juvenile in I is labelled with DAPI and an anti-tyrTub antibody. In B-D, the specimens are viewed from the lateral side with the posterior pigmented structures (ps and pr) placed up in each panel. In E and F, the specimens are viewed from the top of postlarva. Medial optical sections exposing the internally localised subepithelial cells are shown. Insets in B and C show GATA expression in individual subepithelial cells (arrowheads), and an inset in D shows the demarcation between the GATA-negative epidermal layer (ep) and the strongly GATA-positive subepithelial domain (su). An inset in E shows the downregulation of GATA expression in the subepithelial domain (su) relative to the outer layer epithelium at settlement. Arrowheads in F indicate the internal spicule-forming tent-pole-like structures in the juvenile that strongly express GATA. In G, an inset shows a GATA-expressing choanocyte with a cilium (ci) and a basal nucleus (nu). An arrowhead in H shows a cluster of GATA-expressing spicule (sp)-forming cells at the apex of the tent-pole-like structures shown in F; note in I that their cell bodies (arrowhead) are situated basal to the pinacoderm (pi). Scale bars: 100 μm (A-F), 50 μm (H, I, and insets in B-E), 10 μm (G). DAPI, 4',6-diamidino-2-phenylindole.

GATA genes are present in other early branching taxa, including placozoans, ctenophores and the other classes of sponges [3335], although A. queenslandica is the only species for which there currently are published expression patterns. Understanding the evolution of germ layers will be assisted by the analysis of the evolution and expression of genes conserved in bilaterian mesoderm specification and formation, along with further analysis of developmental processes in these taxa. For example, recent analysis of the genome of the ctenophore Mnemiopsis leidyi failed to reveal many ‘mesodermal’ genes despite it appearing to have genuine endomesoderm [34], and it is known that cell movements from one embryonic layer to another occur during metamorphosis in some medusozoan cnidarians [36, 37]. These approaches, along with resolving the exact order in which these early-branching metazoans diverged from the lineage leading to bilaterians (that is, whether sponges or ctenophores are the sister group to all other animals and whether sponges are monophyletic [34, 38]), provide an opportunity to reconstruct the evolution of germ layers and gastrulation.

Conclusions

Comparison of demosponge and eumetazoan development and body plans minimally provides insights into their last common ancestor, which appears to have existed some 800 million years ago [39]. Although the possibility that their last common ancestor possessed germ layers and these were lost in the demosponge lineage cannot be excluded, the lack of developmental commitment of larval and juvenile cell types in A. queenslandica and homology of internal layers across a demosponge and eumetazoans are consistent with eumetazoan germ layers and gastrulation evolving after the divergence of demosponge and eumetazoan lineages. Thus, we propose that ancestrally metazoans consisted of multiple cell layers whose fate was labile. The evolution of the gene regulatory network endowing progressive determination of cell layers resulted in the emergence of primary germ layers and mechanisms to segregate these layers in eumetazoans (that is, gastrulation). The origin of the majority of genes that have conserved roles in eumetazoan gastrulation and germ layer determination evolved after the divergence of demosponge and eumetazoan lineages (for example, Fkh, Twist, Snail), giving credence to this supposition. The conserved expression of GATA in cells located inside the body of all metazoans, along with conserved differential expression of Wnt and TGF-β ligands along demosponge and eumetazoan body axes [40], suggests that these genes were essential in providing critical positional information to cells in the first multicellular metazoans. These, along with eumetazoan-specific genes, later became essential components of the specification and determination of eumetazoan germ layers.

Methods

Animal collection

Adult A. queenslandica were collected at the Heron Island Research Station, Queensland, Australia. Larvae were collected upon release from adults and metamorphosing individuals were obtained by allowing larvae to settle and metamorphose on a plastic petri dish as previously described [41].

Cell lineage tracing with CM-DiI

The lipophilic tracer CM-DiI (Molecular Probes, C7000, Mulgrave, Victoria, AUS) was used to label ciliated epithelial cells in larvae and choanocytes in the juveniles. Larvae or juveniles were incubated in seawater with 10 μM CM-DiI. To differentially label the plasma membrane of ciliated epithelial cells, larvae were incubated for 12 hours. Incubation for one hour was sufficient to differentially label the majority of choanocytes in juveniles. Little labelling occurred in internal cells in larvae or juveniles. The labelled specimens were rinsed in fresh seawater several times and then allowed to develop further. Specimens were fixed as described below at various stages of development. Nuclei were labelled with the fluorescent dye 4',6-diamidino-2-phenylindole (DAPI; 1:1,000, Molecular Probes), and the specimens were mounted in ProlongGold antifade reagent (Molecular Probes).

EdU pulse-and-chase experiments

Larvae or juveniles were incubated in seawater containing 200 μM of the thymidine analogue, EdU (Click-iT EdU AlexaFluor 488 cell proliferation kit, C10337, Molecular Probes), for 16 to 18 hours to label nuclei of proliferating larval archeocytes, and for three hours to label S-phase nuclei in the juveniles. Following washes in fresh seawater, the juveniles were immediately fixed. The larvae with labelled archeocytes were allowed to metamorphose into the juvenile stage, which were then fixed. Fixed specimens underwent immunohistochemistry as described below. Following the immunohistochemistry procedure, fluorescent labelling of incorporated EdU was conducted according to the manufacture’s recommendations prior to DAPI labelling.

Immunohistochemistry, TUNEL and confocal microscopy

Larvae, postlarvae and juveniles were fixed as previously described [42]. Fixed specimens were rehydrated and washed in PBST (phosphate-buffered saline + 0.1% Tween20) and blocked with 1% bovine serum albumin (BSA) in PBST. Primary antibodies were incubated with the specimens in PBST. We used antibodies against phosphorylated histone H3 (phospho S10) (rabbit, 1:500, Abcam ab5176, Cambridge, MA, USA) and tyrosinated-∂-tubulin (mouse, 1:500, Sigma T9028, Castle Hill, NSW, AUS) overnight at 4°C. Following washes in PBST and blocking in 1% BSA in PBST, secondary antibodies were incubated with the specimens; we used AlexaFluor 568 (anti-rabbit or -mouse. 1:200, Molecular Probes) and AlexaFluor 647 (anti-rabbit or -mouse, 1:200, Molecular Probes). Filamentous actin was labelled using AlexaFluor 488-conjugated phallacidin (1:25, Molecular Probes). Specimens were incubated with phallacidin together with secondary antibodies. After washes in PBST, the TUNEL assay was performed by using an In Situ Cell Death Detection Kit (TMR red cat no. 1684795, Roche, Indianapolis, IN, USA) following the manufacturer’s protocol. Nuclei were labelled with DAPI and the specimens were mounted in ProlongGold antifade reagent as described above. Images were recorded using a Zeiss LSM 510 META Confocal Microscope and confocal stacks were viewed using ImageJ.

Gene isolation and whole-mount in situhybridization

AqGATA gene fragments were amplified with gene specific primers, by using complimentary DNA from mixed developmental stages as a PCR template. Gene specific primers were as follows: F1, ATGGAGAAAGCAGACGACGCTATGC; F2, CTACCAAGACAGTTTTGTGG; R1, CTACATCAATTGCTGTGGTTGCATGG and R2, ACTTTGACTTCTTTGACTCG. Primers F2 and R1 were used to generate a 776 base template for a riboprobe synthesis. Riboprobe synthesis and in situ hybridization were conducted as described previously [42]. An antisense digoxigenin-labeled riboprobe was hybridised at a final concentration of 1 ng/μl.

Transmission electron microscopy

Free-swimming larvae were processed for transmission electron microscopy (TEM) by using the high-pressure freezing (HPF) technique. Procedures followed the previously established HPF protocol for shark retinal tissues [43] with the following modifications: 20% BSA in artificial seawater was used as a cryoprotectant instead of the teleost ringer/0.7% agarose. In addition, 1% OsO4 was used in the freeze substitution cocktail. Longitudinal ultrathin sections were prepared and analysed.

References

  1. Hyman LH: The Invertebrates: Protozoa Through Ctenophora, Volume 1. 1940, McGraw-Hill Book Company, Inc: New York and London

    Google Scholar 

  2. Hall BK: Germ layers and the germ-layer theory revisited - primary and secondary germ layers, neural crest as a fourth germ layer, homology, and demise of the germ-layer theory. Evol Biol. 1998, 30: 121-186.

    Article  Google Scholar 

  3. Martindale MQ: The evolution of metazoan axial properties. Nat Rev Genet. 2005, 6: 917-927. 10.1038/nrg1725.

    Article  CAS  PubMed  Google Scholar 

  4. Martindale MQ, Pang K, Finnerty JR: Investigating the origins of triploblasty: ‘mesodermal’ gene expression in a diploblastic animal, the sea anemone Nematostella vectensis (phylum, Cnidaria; class, Anthozoa). Development. 2004, 131: 2463-2474. 10.1242/dev.01119.

    Article  CAS  PubMed  Google Scholar 

  5. Nakanishi N, Renfer E, Technau U, Rentzsch F: Nervous systems of the sea anemone Nematostella vectensis are generated by ectoderm and endoderm and shaped by distinct mechanisms. Development. 2012, 139: 347-357. 10.1242/dev.071902.

    Article  CAS  PubMed  Google Scholar 

  6. Erpenbeck D, Worheide G: On the molecular phylogeny of sponges (Porifera). Zootaxa. 2007, 1668: 107-126.

    Google Scholar 

  7. Gazave E, Lapebie P, Ereskovsky AV, Vacelet J, Renard E, Cardenas P, Borchiellini C: No longer demospongiae: Homoscleromorpha formal nomination as a fourth class of porifera. Hydrobiologia. 2012, 687: 3-10. 10.1007/s10750-011-0842-x.

    Article  Google Scholar 

  8. Ereskovsky AV: The Comparative Embryology of Sponges. 2010, Dordrecht; New York: Springer

    Book  Google Scholar 

  9. Leys SP: Gastrulation in sponges. Gastrulation, From Cells to Embryo. Edited by: Stern CD. 2004, New York: Cold Spring Harbor Laboratory Press, 23-31.

    Google Scholar 

  10. Leys SP, Degnan BM: Embryogenesis and metamorphosis in a haplosclerid demosponge: gastrulation and transdifferentiation of larval ciliated cells to choanocytes. Invert Biol. 2002, 121: 171-189.

    Article  Google Scholar 

  11. Ereskovsky AV, Boury-Esnault N: Cleavage pattern in Oscarella species (Porifera, Demospongiae, Homoscleromorpha): transmission of maternal cells and symbiotic bacteria. J Nat Hist. 2002, 36: 1761-1775. 10.1080/00222930110069050.

    Article  Google Scholar 

  12. Boury-Esnault N, Efremova S, Bezac C, Vacelet J: Reproduction of a hexactinellid sponge: first description of gastrulation by cellular determination in the Porifera. Invert Reprod Dev. 1999, 35: 187-201. 10.1080/07924259.1999.9652385.

    Article  Google Scholar 

  13. Leys SP, Cheung E, Boury-Esnault N: Embryogenesis in the glass sponge Oopsacas minuta: formation of syncytia by fusion of blastomeres. Integr Comp Biol. 2006, 46: 104-117. 10.1093/icb/icj016.

    Article  CAS  PubMed  Google Scholar 

  14. Gonobobleva E, Ereskovsky A: Polymorphism in free-swimming larvae of Halisarca dujardini (Demospongiae, Halisarcida). Sponge Science in the New Millennium: Papers Contributed to the VI International Sponge Conference, Rapallo (Italy), 29 September-5 October 2002, Volume 68. Edited by: Pansini M, Pronzato R, Bavestrello G, Manconi R. 2004, Genova, Italy: Boll Mus Ist Biol. Universitá di Genova, 349-356.

    Google Scholar 

  15. Ereskovsky AV, Gonobobleva E: New data on embryonic development of Halisarca dujardini Johnson, 1842 (Demospongiae, Halisarcida). Zoosystema. 2000, 22: 355-368.

    Google Scholar 

  16. Leys SP, Eerkes-Medrano D: Gastrulation in calcareous sponges: in search of Haeckel’s gastraea. Integr Comp Biol. 2005, 45: 342-351. 10.1093/icb/45.2.342.

    Article  PubMed  Google Scholar 

  17. Simpson TL: Cell Biology of Sponges. 1984, New York: Springer-Verlag

    Book  Google Scholar 

  18. Ereskovsky AV, Korotkova GP: The reasons of sponge sexual morphogenesis peculiarities. Berliner Geowiss Abh. 1997, 20: 25-

    Google Scholar 

  19. Bergquist PR, Green CR: Ultrastructural-study of settlement and metamorphosis in sponge larvae. Cahiers De Biologie Marine. 1977, 18: 289-302.

    Google Scholar 

  20. Meewis H: Contribution a l’étude de l’embryogenése de chalinulidae: Haliclona limbata. Ann Soc R Zool Belg. 1939, 70: 201-243.

    Google Scholar 

  21. Misevic GN, Burger MM: The molecular basis of species specific cell-cell recognition in marine sponges, and a study on organogenesis during metamorphosis. Prog Clin Biol Res B. 1982, 85: 193-209.

    CAS  Google Scholar 

  22. Misevic GN, Schlup V, Burger MM: Larval metamorphosis of Microciona prolifera: evidence against the reversal of layers. New Perspectives in Sponge Biology. Edited by: Rutzler K. 1990, Washington, DC: Smithsonian Institution Press, 182-187.

    Google Scholar 

  23. Minchin EA: Note on the larva and the postlarval development of Leucosolenia variabilis, H. sp., with remarks on the development of other Asconidae. Proc R Soc Lond. 1896, 60: 42-52. 10.1098/rspl.1896.0013.

    Article  Google Scholar 

  24. Amano S, Hori I: Metamorphosis of coeloblastula performed by multipotential larval flagellated cells in the calcareous sponge Leucosolenia laxa. Biol Bull. 2001, 200: 20-32. 10.2307/1543082.

    Article  CAS  PubMed  Google Scholar 

  25. Amano S, Hori I: Metamorphosis of calcareous sponges. 2. Cell rearrangement and differentiation in metamorphosis. Invert Reprod Dev. 1993, 24: 13-26. 10.1080/07924259.1993.9672327.

    Article  Google Scholar 

  26. Ereskovsky AV, Tokina DB, Bezac C, Boury-Esnault N: Metamorphosis of cinctoblastula larvae (Homoscleromorpha, Porifera). J Morphol. 2007, 268: 518-528. 10.1002/jmor.10506.

    Article  PubMed  Google Scholar 

  27. Ereskovsky AV, Renard E, Borchiellini C: Cellular and molecular processes leading to embryo formation in sponges: evidences for high conservation of processes throughout animal evolution. Dev Genes Evol. 2013, 223: 5-22. 10.1007/s00427-012-0399-3.

    Article  CAS  PubMed  Google Scholar 

  28. Srivastava M, Simakov O, Chapman J, Fahey B, Gauthier ME, Mitros T, Richards GS, Conaco C, Dacre M, Hellsten U, Larroux C, Putnam NH, Stanke M, Adamska M, Darling A, Degnan SM, Oakley TH, Plachetzki DC, Zhai YF, Adamski M, Calcino A, Cummins SF, Goodstein DM, Harris C, Jackson DJ, Leys SP, Shu SQ, Woodcroft BJ, Vervoort M, Kosik KS, et al: The Amphimedon queenslandica genome and the evolution of animal complexity. Nature. 2010, 466: 720-723. 10.1038/nature09201.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  29. Rottinger E, Dahlin P, Martindale MQ: A framework for the establishment of a cnidarian gene regulatory network for “endomesoderm” specification: the inputs of ss-catenin/TCF signaling. PLoS Genetics. 2012, 8: e1003164-10.1371/journal.pgen.1003164.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  30. Gillis WJ, Bowerman B, Schneider SQ: Ectoderm- and endomesoderm-specific GATA transcription factors in the marine annelid Platynereis dumerilli. Evol Dev. 2007, 9: 39-50. 10.1111/j.1525-142X.2006.00136.x.

    Article  CAS  PubMed  Google Scholar 

  31. Chiodin M, Borve A, Berezikov E, Ladurner P, Martinez P, Hejnol A: Mesodermal gene expression in the acoel Isodiametra pulchra indicates a low number of mesodermal cell types and the endomesodermal origin of the gonads. PLoS One. 2013, 8: e55499-10.1371/journal.pone.0055499.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  32. Boyle MJ, Seaver EC: Developmental expression of foxA and gata genes during gut formation in the polychaete annelid, Capitella sp I. Evol Dev. 2008, 10: 89-105. 10.1111/j.1525-142X.2007.00216.x.

    Article  CAS  PubMed  Google Scholar 

  33. Srivastava M, Begovic E, Chapman J, Putnam NH, Hellsten U, Kawashima T, Kuo A, Mitros T, Salamov A, Carpenter ML, Signorovitch AY, Moreno MA, Kamm K, Grimwood J, Schmutz J, Shapiro H, Grigoriev IV, Buss LW, Schierwater B, Dellaporta SL, Rokhsar DS: The Trichoplax genome and the nature of placozoans. Nature. 2008, 454: 955-960. 10.1038/nature07191.

    Article  CAS  PubMed  Google Scholar 

  34. Ryan JF, Pang K, Schnitzler CE, Nguyen AD, Moreland RT, Simmons DK, Koch BJ, Francis WR, Havlak P, Comparative Sequencing Program NISC, Smith SA, Putnam NH, Haddock SH, Dunn CW, Wolfsberg TG, Mullikin JC, Martindale MQ, Baxevanis AD: The genome of the ctenophore Mnemiopsis leidyi and its implications for cell type evolution. Science. 2013, 342: 1242592-10.1126/science.1242592.

    Article  PubMed Central  PubMed  Google Scholar 

  35. Riesgo A, Farrar N, Windsor PJ, Giribet G, Leys SP: The analysis of eight transcriptomes from all poriferan classes reveals surprising genetic complexity in sponges. Mol Biol Evol. 2014, doi:10.1093/molbev/msu057

    Google Scholar 

  36. Yuan D, Nakanishi N, Jacobs DK, Hartenstein V: Embryonic development and metamorphosis of the scyphozoan Aurelia (Cnidaria, Scyphozoa). Dev Genes Evol. 2008, 218: 525-539. 10.1007/s00427-008-0254-8.

    Article  PubMed  Google Scholar 

  37. Weis V, Buss L: Ultrastructure of metamorphosis in Hydractinia echinata. Postilla. 1987, 199: 1-20.

    Google Scholar 

  38. Sperling EA, Peterson KJ, Pisani D: Phylogenetic-signal dissection of nuclear housekeeping genes supports the paraphyly of sponges and the monophyly of Eumetazoa. Mol Biol Evol. 2009, 26: 2261-2274. 10.1093/molbev/msp148.

    Article  CAS  PubMed  Google Scholar 

  39. Erwin DH, Laflamme M, Tweedt SM, Sperling EA, Pisani D, Peterson KJ: The Cambrian conundrum: early divergence and later ecological success in the early history of animals. Science. 2011, 334: 1091-1097. 10.1126/science.1206375.

    Article  CAS  PubMed  Google Scholar 

  40. Adamska M, Degnan SM, Green KM, Adamski M, Craigie A, Larroux C, Degnan BM: Wnt and TGF-beta expression in the sponge Amphimedon queenslandica and the origin of metazoan embryonic patterning. PLoS One. 2007, 2: e1031-10.1371/journal.pone.0001031.

    Article  PubMed Central  PubMed  Google Scholar 

  41. Leys SP, Larroux C, Gauthier M, Adamska M, Fahey B, Richards GS, Degnan SM, Degnan BM: Isolation of Amphimedon developmental material. CSH protocols. 2008, 2008: pdb.prot5095

    Google Scholar 

  42. Larroux C, Fahey B, Liubicich D, Hinman VF, Gauthier M, Gongora M, Green K, Worheide G, Leys SP, Degnan BM: Developmental expression of transcription factor genes in a demosponge: insights into the origin of metazoan multicellularity. Evol Dev. 2006, 8: 150-173. 10.1111/j.1525-142X.2006.00086.x.

    Article  CAS  PubMed  Google Scholar 

  43. Harahush BK, Green K, Webb R, Hart NS, Collin SP: Optimal preservation of the shark retina for ultrastructural analysis: an assessment of chemical, microwave, and high-pressure freezing fixation techniques. Microsc Res Tech. 2012, 75: 1218-1228. 10.1002/jemt.22052.

    Article  PubMed  Google Scholar 

  44. Woollacott RM: Structure and swimming behavior of the larva of Haliclona tubifera (Porifera, Demospongiae). J Morphol. 1993, 218: 301-321. 10.1002/jmor.1052180306.

    Article  Google Scholar 

  45. Felsenfeld G, Evans T, Reitman M: An erythrocyte-specific DNA-binding factor recognizes a regulatory sequence common to all chicken globin genes. Proc Natl Acad Sci U S A. 1988, 85: 5976-5980. 10.1073/pnas.85.16.5976.

    Article  PubMed Central  PubMed  Google Scholar 

  46. Mnemiopsis Genome Project Portal. [http://research.nhgri.nih.gov/mnemiopsis/]

  47. Basic Local Alignment Search Tool. [http://blast.ncbi.nlm.nih.gov/Blast.cgi]

  48. Acropora digitifera Genome (Ver 1.1). [http://marinegenomics.oist.jp/genomes/viewer?project_id=3&current_assembly_version=oist_v1.1]

  49. Compagen: A comparative genomics platform for early branching Metazoa. [http://compagen.zoologie.uni-kiel.de/]

  50. Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  51. Gascuel O, Guindon S: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.

    Article  PubMed  Google Scholar 

  52. Whelan S, Goldman N: A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach. Mol Biol Evol. 2001, 18: 691-699. 10.1093/oxfordjournals.molbev.a003851.

    Article  CAS  PubMed  Google Scholar 

  53. Simionato E, Ledent V, Richards G, Thomas-Chollier M, Kerner P, Coornaert D, Degnan BM, Vervoort M: Origin and diversification of the basic helix-loop-helix gene family in metazoans: insights from comparative genomics. BMC Evol Biol. 2007, 7: 33-10.1186/1471-2148-7-33.

    Article  PubMed Central  PubMed  Google Scholar 

  54. Larroux C, Luke GN, Koopman P, Rokhsar DS, Shimeld SM, Degnan BM: Genesis and expansion of metazoan transcription factor gene classes. Mol Biol Evol. 2008, 25: 980-996. 10.1093/molbev/msn047.

    Article  CAS  PubMed  Google Scholar 

  55. Larroux C, Fahey B, Degnan SM, Adamski M, Rokhsar DS, Degnan BM: The NK homeobox gene cluster predates the origin of Hox genes. Curr Biol. 2007, 17: 706-710.

    Article  CAS  PubMed  Google Scholar 

  56. Ryan JF, Mazza ME, Pang K, Matus DQ, Baxevanis AD, Martindale MQ, Finnerty JR: Pre-bilaterian origins of the Hox cluster and the Hox code: evidence from the sea anemone, Nematostella vectensis. PLoS One. 2007, 2: e153-10.1371/journal.pone.0000153.

    Article  PubMed Central  PubMed  Google Scholar 

  57. Kerner P, Degnan SM, Marchand L, Degnan BM, Vervoort M: Evolution of RNA-binding proteins in animals: insights from genome-wide analysis in the sponge Amphimedon queenslandica. Mol Biol Evol. 2011, 28: 2289-2303. 10.1093/molbev/msr046.

    Article  CAS  PubMed  Google Scholar 

  58. Adamska M, Larroux C, Adamski M, Green K, Lovas E, Koop D, Richards GS, Zwafink C, Degnan BM: Structure and expression of conserved Wnt pathway components in the demosponge Amphimedon queenslandica. Evol Dev. 2010, 12: 494-518. 10.1111/j.1525-142X.2010.00435.x.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank the staff at the Heron Island Research Station for the use of the facilities, Dr. Kathryn Green at the Centre for Microscopy and Microanalysis at the University of Queensland for preparing EM sections, and Degnan lab members for constructive discussions. This research was supported by an Australian Research Council grant to BMD.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Bernard M Degnan.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

NN and SS carried out the experiments and microscopy. NN and BMD conceived of the study and wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12915_2014_761_MOESM1_ESM.tiff

Additional file 1: Figure S1: Stages of development during metamorphosis in A. queenslandica. A: Free-swimming larva stage; the (parenchymella) larva is egg-shaped and swims with the pigmented ring and associated bundles of long cilia directed posteriorly [10]. The anterior is down. In B-F, metamorphosing postlarvae are viewed from the top. B: Settlement stage. The anterior region of the larva is attached to the substrate, onto which the larva flattens. C: Mat-formation postlarval stage, cells of the metamorphosing postlarva migrate laterally on the substrate to form a mat-like structure. Note that former posterior-ring pigmentation (pi) is still visible but is disappearing. D: Chamber-formation postlarval stage, with developing canals (ca) lined by choanocyte chambers (ch) and endopinacoderm. Note the aquiferous system composed of canals lined by choanocytes and endopinacocytes becomes evident. E: Tent-pole-formation postlarval stage, the exopinacotes covering the outer surface of the metamorphosing postlarva are lifted upwards by formation of tent-pole-like structures consisting of vertically oriented clusters of spicules and associated cells. Arrowheads show the internal tent-pole-like structures, visible here as clustering of cells. F: Juvenile (rhagon) stage with an osculum (os), marking the establishment of the functional aquiferous system. Abbreviations: lc long cilia; pr pigment ring; cc cuboidal cells. Scale bar: 100 μm (A, inset in F), 1 mm (B-F). (TIFF 1 MB)

12915_2014_761_MOESM2_ESM.jpeg

Additional file 2: Figure S2: Apoptosis, loss of epithelial integrity and phagocytosis at early metamorphosis in A. queenslandica. A-C: Free-swimming larva. D-F: Settled postlarva, ≤7 hours post-settlement. G: Settled postlarva, ≤12 hours post-settlement. H: Mat-formation postlarval stage. DAPI is used to label DNA. In A, B, D, E and G, DNA fragmentation is detected by TUNEL assays (TUNEL). In C, F and G, specimens are labelled with an anti-Tyrosinated-Tubulin (tyrTub). A: A confocal longitudinal section through the larva; anterior is down. A boxed region in A is shown in B; note that TUNEL-positive fragmented DNA (arrowheads) localises within the most anterior domain where cuboidal cells reside. C: A confocal section through the outer layer epithelium; note that the epithelium consists of a continuous sheet of cells with apical cilia (ci) and basal, typically anucleolated, nuclei (nu). D: Confocal medial sections viewed from the top of the metamorphosing postlarva. A single section through the outer layer epithelium in the boxed region in D is shown in E. F: A confocal section showing the surface view of the distintegrating outer layer epithelium; note the breakup into smaller, tightly packed cells with anucleolated nuclei (the fragmented epithelium; fr). G: A confocal section of an outer-most region of the metamorphosing postlarva; note the lack of an epithelium, and the presence of phagocytised TUNEL-positive apoptotic bodies (arrowheads). H: A confocal section through an archeocyte in a mat-formation-stage postlarva labelled with DAPI and an anti-tyrTub, in which larval archeocytes and ciliated epidermal cells were simultaneously labelled (with EdU and CM-DiI, respectively) and followed through metamorphosis. Note intracellular CM-DiI-positive cellular fragments containing DAPI-positive nuclear materials (arrowhead), the presumptive apoptotic bodies, in an EdU-positive archeocyte of larval origin. Abbreviations: ep outer layer epithelium; icm inner cell mass. Scale bars: 100 μm (A, D, E), 10 μm (B, C, F-H). (JPEG 398 KB)

12915_2014_761_MOESM3_ESM.jpeg

Additional file 3: Figure S3: Larval ciliated epithelial cells lose cilia and transdifferentiate into choanocytes via an archeocyte intermediate during metamorphosis in A. queenslandica. Confocal sections through a flask cell (fc) at the free-swimming larva stage (A) and a transdifferentiating flask cell at the settlement stage labelled with CM-DiI (B). Note the loss of apical structures including the cilium (ci). C: Confocal sections viewed from the top at settlement, showing that many DiI-labelled cells are internalised (arrowheads). D: Confocal sections viewed from the top during metamorphosis; note that DiI-labelled cells have differentiated into archeocytes with characteristic large nucleoli (nu; inset magnified from the boxed cells). E: Confocal sections of individual archeocytes of the larval epidermal origin undergoing further differentiation at the chamber-formation stage; an arrowhead shows a differentiating choanocyte with a cilium (ci in inset) whose cell body has yet to complete cytokinesis within the mother archeocyte (ar). F: Confocal sections of a multi-nucleated archeocyte-choanocyte intermediate at the chamber-formation stage. Differentiating choanocytes with cilia (ci) are evident (arrowheads). Abbreviations: ap apical non-ciliary process; ba basal process; co columnar epithelial cell. Scale bar: 10 μm (A, B, inset in D, E, F), 100 μm (C, D). (JPEG 290 KB)

12915_2014_761_MOESM4_ESM.jpeg

Additional file 4: Figure S4: EdU is incorporated into proliferating archeocytes in the inner cell mass in A. queenlandica larvae. A, B: Confocal sections of the inner cell mass in a free-swimming larva pulse-labelled with EdU. Nuclei are labelled with DAPI. Arrowheads show EdU-positive nuclei. The boxed region in B is magnified in an inset. Note that the EdU-positive nuclei are 3.5 to 4.0 μm in diameter and contain prominent nucleoli (nu), as observed in larval archeocytes in the inner cell mass in another demosponge Haliclona tubifera[44]. C: A transmission electron micrograph showing prominent nucleoli (nu) in nuclei of archeocytes (ar) in the inner cell mass within an A. queenslandica larva. Scale bar: 10 μm (A, B), 5 μm (C). (JPEG 476 KB)

12915_2014_761_MOESM5_ESM.jpeg

Additional file 5: Figure S5: Epithelial choanocytes can de-differentiate into archeocytes in A. queenslandica. Choanocytes in chambers (ch) in juveniles were labelled with CM-DiI as in Figure 3A, B and descendants of the labelled cells in juveniles one day later are shown. Some choanocytes are seen to extend basal processes (ba) laterally, engulfing neighboring choanocytes (arrowheads in A, B). Some cells appear multinucleated (arrowheads in C), presumably as a result of phagocytic or fusion events. Abbreviations: ch choanocyte chamber; ci cilium. Scale bar: 10 μm. (JPEG 181 KB)

12915_2014_761_MOESM6_ESM.jpeg

Additional file 6: Figure S6: GATA sequence alignment and phylogenetic analyses. A: Partial GATA protein multiple sequence alignment showing the amino acid sites used for phylogenetic analyses. GATA genes encode zinc-finger transcription factors that bind to the (T/A (GATA) A/G) cis-regulatory motif [45] and are characterised by the presence of two zinc-finger domains (N-terminal ‘N-finger’ and C-terminal ‘C-finger’). The alignment encompasses N-finger (underlined in red) and C-finger (underlined in blue) zinc-finger domains each encoding CXNC (boxed) X17CNXC (boxed). B: GATA maximum likelihood phylogeny rooted with GATA-type sequences. Amphimedon queenslandica sequence is highlighted in red. No clear GATA orthologues were identified outside Metazoa, and in the genome of the ctenophore Mnemiopsis leidyi[46]. C: Unrooted GATA maximum likelihood phylogeny. Sequence alignment and phylogenetic analyses were performed on the Geneious platform (v.5.1.7). Related sequences were retrieved via the protein BLAST search using the Amphimedon GATA-like sequence as queries, from GenBank at NCBI [47], Acropora digitifera genome (Version 1.1) portal at OIST [48], Mnemiopsis genome project portal at NIH [46] and Compagen at the University of Kiel (Oscarella sequences; [49]). Sequence IDs/accession numbers are shown with the name of each sequence in B and C. Peptide sequences were aligned with MUSCLE (v3.7) [50] (default settings), and ambiguous regions were manually removed. Phylogenetic trees were reconstructed using the maximum likelihood method implemented in the PhyML program [51]. The WAG substitution model [52] was selected assuming an estimated proportion of invariant sites and four gamma-distributed rate categories to account for rate heterogeneity across sites. The gamma shape parameter was estimated directly from the data. Reliability for internal branches of maximum likelihood trees was assessed using the bootstrapping method (100 bootstrap replicates). Support values are shown at each node except when lower than 50%. The unit of the branch length is the number of substitutions per site. (JPEG 1 MB)

12915_2014_761_MOESM7_ESM.xlsx

Additional file 7: Table S1: A list of the presence or absence of A. queenslandica homologues (orthologues and gene family members) to cnidarian Nematostella vectensis endomesodermal genes (listed in [29]) based on previously published data [28, 5358] and sequence similarity searches in NCBI databases (via BLASTp). (XLSX 11 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( https://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Cite this article

Nakanishi, N., Sogabe, S. & Degnan, B.M. Evolutionary origin of gastrulation: insights from sponge development. BMC Biol 12, 26 (2014). https://doi.org/10.1186/1741-7007-12-26

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1741-7007-12-26

Keywords