Cleavage modification did not alter blastomere fates during bryozoan evolution
BMC Biology volume 15, Article number: 33 (2017)
Stereotypic cleavage patterns play a crucial role in cell fate determination by precisely positioning early embryonic blastomeres. Although misplaced cell divisions can alter blastomere fates and cause embryonic defects, cleavage patterns have been modified several times during animal evolution. However, it remains unclear how evolutionary changes in cleavage impact the specification of blastomere fates. Here, we analyze the transition from spiral cleavage – a stereotypic pattern remarkably conserved in many protostomes – to a biradial cleavage pattern, which occurred during the evolution of bryozoans.
Using 3D-live imaging time-lapse microscopy (4D-microscopy), we characterize the cell lineage, MAPK signaling, and the expression of 16 developmental genes in the bryozoan Membranipora membranacea. We found that the molecular identity and the fates of early bryozoan blastomeres are similar to the putative homologous blastomeres in spiral-cleaving embryos.
Our work suggests that bryozoans have retained traits of spiral development, such as the early embryonic fate map, despite the evolution of a novel cleavage geometry. These findings provide additional support that stereotypic cleavage patterns can be modified during evolution without major changes to the molecular identity and fate of embryonic blastomeres.
Cleavage is the sequence of cell divisions that turns a zygote into a multicellular embryo, and plays an essential role in the specification of cell fates before the onset of gastrulation. A cleavage pattern can be variable, where the blastomere positions are not predictable (e.g., mouse), or stereotypic (e.g., ascidian), where the embryonic cell divisions form a precise, identifiable three-dimensional pattern of blastomeres . There is evidence that different types of cleavage can dictate different underlying mechanisms of cell fate specification ; however, it is still largely unknown how a stereotypic cleavage pattern affects the evolution of animal morphology . Cleavage patterns are highly diverse, they can even differ between closely related species [3, 4], or remain conserved in different animal lineages over long evolutionary periods . A notable example of the latter is known as spiral cleavage, and is a rich framework to investigate the relation between development and evolution.
Spiral cleavage occurs in molluscs, annelids, nemerteans, and polyclad flatworms [6–16]. In these groups, the fertilized eggs divide through a highly stereotypic cleavage pattern where blastomeres at the 4-cell stage cleave with the mitotic spindles oblique to the animal–vegetal axis, alternating direction (clockwise and counterclockwise) at each division cycle, termed the spiral cleavage pattern [5, 6, 17–19]. This determinate developmental mode allowed for the identification of homologous blastomeres across taxa and unprecedented detail in the comparison of animal embryogenesis, further revealing that spiral-cleaving embryos not only have the same cleavage pattern, but homologous blastomeres between groups have a similar fate in the larval and adult tissues [5, 18]. The study of spiral cleavage thus revealed that, in contrast to late developmental stages, early development can remain conserved for extended evolutionary periods, shaping our current understanding about the relation between ontogeny and phylogeny [20–24].
Even though spiral cleavage has been modified in a multitude of ways throughout evolution with changes in blastomere sizes and cell fate specification [5, 17–19], the cleavage pattern itself remained fairly conserved. Known cases where the spiral cleavage pattern was lost is usually associated with drastic developmental changes, such as the transition to a syncytial blastoderm in cephalopods , or the evolution of extra-embryonic yolk cells in platyhelminthes . However, the recent improvements in the resolution of protostome relationships revealed that the spiral cleavage pattern has been drastically modified or lost in even more groups than previously thought .
Spiral cleavage is a synapomorphy for the Spiralia (Lophotrochozoa sensu lato, after ), a major protostome clade containing all spiral-cleaving groups . However, not all spiralians (i.e., animals belonging to the clade Spiralia) display a spiral cleavage pattern during embryogenesis. Recent spiralian phylogenies [27–32] indicate that clades that do not exhibit oblique cell divisions, such as the bryozoans , brachiopods , gastrotrichs , and rotifers , must have modified or lost the ancestral spiral cleavage pattern during evolution [5, 37] (Fig. 1). For this reason, such groups are essential to understand how cleavage patterns and blastomere fates evolve and can uniquely reveal which developmental traits, if any, remained conserved in the evolutionary transition from spiral to a derived cleavage geometry.
In the current work, we investigate the development of a group that lost the spiral cleavage pattern during evolution – the bryozoans. These sessile colonial invertebrates occur in oceans worldwide and have fairly diverse reproductive strategies and larval stages [38, 39], but none of the species investigated so far display a spiral arrangement of embryonic blastomeres . Bryozoans display a unique stereotypic cleavage pattern with a biradial arrangement of the blastomeres that is widely conserved within the group.
Previous studies of bryozoan embryology [40–45] suggest that the animal-most blastomeres give rise to the apical disc and aboral epithelium of the larva, the vegetal-most derivatives of the animal blastomeres form the ciliated band, and the vegetal blastomeres produce the oral epithelium and endomesoderm [38, 39, 46]. This coarse fate map appears to be overall similar to that of spiral-cleaving embryos [47, 48]. However, cleavage patterns have only been systematically followed until the 64-cell stage [43, 45] and, as of today, there is no detailed cell lineage or fate map of a bryozoan larva. Several basic developmental questions remain unsolved. For example, the relation between the embryonic animal–vegetal axis and the larval body axes is unclear , and the fate of the blastopore remains to be confirmed [39, 41, 44, 49]. Finally, the fate of internalized cells has not been traced  and the source of mesoderm remains an especially contentious topic .
In this study, we investigate the embryogenesis of the cosmopolitan gymnolaemate species Membranipora membranacea (Linnaeus, 1767) to understand the evolutionary transition from a spiral to a biradial cleavage pattern. We take advantage of the vast cell lineage data available for spiralians and the growing literature on spiralian gene expression, to compare the molecular identity and fate of embryonic blastomeres between the bryozoan and other spiral-cleaving embryos with cellular resolution. We were able to identify the embryonic source of most larval tissues of M. membranacea based on 4D microscopy recordings, and to combine this cell lineage data with the activity of the MAPK pathway and expression of several conserved developmental markers, generating a detailed overview of the blastomere identities and fates in the bryozoan. The comparison to a typical spiral development reveals that the early blastomeres of M. membranacea share similar molecular identities and fates with other spiral-cleaving embryos, despite the contrasting cleavage pattern. Given the phylogenetic position of bryozoans, we suggest these coincident developmental traits were inherited from a spiral-cleaving ancestor during the evolutionary transition from spiral to biradial cleavage. The findings support the hypothesis that stereotypic cleavage patterns can be modified during evolution without major changes to blastomere gene expression and fates. Our study highlights the power of the comparative approach to address fundamental questions of development and evolution, such as the relation between cleavage patterns and fate maps.
General development and data overview
Colonies of M. membranacea spawn fertilized discoidal eggs into the water column . The released eggs undergo activation, quickly become spherical (Fig. 2a), and initiate cleavage at around 2 hours post activation (hpa) with a discernible accumulation of yolk at the vegetal pole (Fig. 2b). Throughout development, the embryo maintains close contact with the fertilization envelope via abundant cytoplasmic extensions (Fig. 2a, i). The yolky cells at the vegetal pole are internalized during gastrulation (Fig. 2c–e, j–m) and, by the mid gastrula stage (16 hpa), the primordia of the apical organ (apical disc) and of the ciliated band (corona) are visible (Fig. 2e). The vegetal plate invaginates and the embryo elongates along the animal–vegetal axis forming a late gastrula at 24 hpa with clearly defined larval structures (i.e., apical organ, shell, gut, and corona) (Fig. 2f, g). At this point the fertilization envelope opens at the animal and vegetal ends and the embryo begins to swim by ciliary beating (Fig. 2f). The internal cavity (vestibule) widens in the anteroposterior axis resulting in the typical laterally compressed, triangular shaped and shelled feeding larva of gymnolaemate bryozoans – the cyphonautes (Fig. 2h, i) .
In this study, we recorded four individual embryos from the 2-cell stage at 2 hpa (Fig. 2b) until the late gastrula stage at 24 hpa (Fig. 2f, g). Due to the opaqueness of the embryo, it was necessary to trace animal and vegetal clone populations in different individuals, namely “wild type 1” (wt1) and wt2, respectively (Additional file 1: Video S1). To evaluate the potential variability between individual embryos, we recorded two additional animal pole views (embryos wt3 and wt4, Additional file 2: Video S2). We were able to trace a greater number of individual cells until the mid gastrula stage at 16 hpa (Additional file 3: Figure S1), when the primordia of most larval structures were clearly defined and the fate of the cells could be determined (Fig. 2e).
Additional file 1: Video S1. 4D recordings of M. membranacea embryos wt1 and wt2. Embryo wt1 was recorded from the animal pole (left) and wt2 was recorded from the vegetal pole (right). We focused the different focal planes of each time point into a single image (focus stacking) to reveal the overall embryo morphology. The time-lapse stack-focused images were animated at 25 frames per second (1000× acceleration) and the developmental stages of both embryos were synchronized by accelerating the video of wt2 by 10% in relation to wt1. Duration = 21.8 h (from ~3 hpa to ~25 hpa). (MP4 19251 kb)
Additional file 2: Video S2. 4D recordings of M. membranacea embryos wt3 and wt4. Embryos wt3 (left) and wt4 (right) were recorded from the animal pole. We focused different focal planes of each time point into a single image (focus stacking) to reveal the overall embryo morphology. The time-lapse stack-focused images were animated at 25 frames per second (1000× acceleration) and the developmental stages of both embryos were synchronized by accelerating the video of wt3 by 5% in relation to wt4. Duration = 23.0 h (from ~3 hpa to ~25 hpa). (MP4 19353 kb)
Overall, the data we collected suggests the cell lineage of individual M. membranacea embryos is highly stereotypic and exhibits small variation in the timing of cell divisions (see below). However, due to our limited sample size, we cannot fully account for the cell fate variability that might exist in bryozoan development, particularly at later stages. The results we report below thus reflect the consensus data between the four M. membranacea embryos tracked in this study.
Cleavage pattern and embryonic axes
The cleavage of M. membranacea is biradial as previously described for gymnolaemates (Fig. 3) [38, 39, 48, 49, 52]. At 15 °C, the first cell division occurs between 1 and 2 hpa and produces two equal blastomeres with a meridional cleavage furrow. The second division is also meridional and perpendicular to the first, resulting in four blastomeres of equal sizes around 3 hpa (Fig. 2j; Fig. 3; Fig. 4, 4-cell). We labeled the blastomere that gives rise to the posterior structures of the larval body as “D” (see Fig. 5 for fate map overview and “Methods” for nomenclature details). In most embryos, the cell sister of the D blastomere gives rise to the right side of the embryo . At 4 hpa, an equatorial third division gives rise to four animal blastomeres with lower yolk content (1a–1d), and four equally sized vegetal blastomeres with a greater amount of yolk displaced towards the center of the embryo (1A–1D) (Fig. 2c; Fig. 3; Fig. 4, 8-cell). During the next division at 5.2 hpa, each animal blastomere divides meridionally, parallel to the plane of the first cleavage, forming a 16-cell stage embryo that clearly differs from the canonical spiral cleavage pattern (Fig. 2k; Fig. 3; Fig. 4, 16-cell). Since these blastomeres occupy the same position along the animal–vegetal axis, and thus cannot be objectively labeled with superscript 1 or 2, they received the subscript i or e to indicate their internal or external position in relation to the central axis of the embryo (see “Methods” for nomenclature details). The vegetal blastomeres cleave in the same manner, but slightly after. At the 16-cell stage (5.2 hpa), yolk-rich cells (2A–2D) lie inner to the outer vegetal cells of the second quartet (2a–2d) and the embryo is clearly biradial.
During subsequent stages, the eight animal blastomeres of M. membranacea act as octets, dividing synchronously (Fig. 3). The first octet (animal pole cells 1q, four inner and four outer cells) divides equatorially making a brief 24-cell stage and the octets 1q1 and 1q2 (6.5 hpa). This division is shortly followed by an unequal cleavage originating the third quartet (3a–3d) from the four inner vegetal blastomeres at 6.7 hpa (Fig. 2l; Fig. 3; Fig. 4, 28-cell). Outer vegetal cells of the second quartet (2q) divide parallel to the second division at 7.5 hpa, resulting in 12 outer vegetal cells (3a–3d, 2aR–2dR, 2aL–2dL) that surround four large blastomeres in the vegetal plate (3A–3D) at the 32-cell stage. At 8 hpa, the top animal octet (1q1) divides, forming a 40-cell embryo (Fig. 3; Fig. 4, 40-cell). Finally, the vegetal most animal octet (1q2) divides meridionally at 8.6–9 hpa forming an equatorial row of cells above the vegetal blastomeres (Fig. 3; Fig. 4, 48-cell).
Cell lineage variability
We found little variation between the cell lineages of the four embryos. That is, a particular cell in one embryo has the same lineage history, occupies the same relative position and divides roughly at the same time as the respective cell in a different embryo. A direct comparison between the four M. membranacea embryos reveals that the lineages overlap well, exhibiting only small variations in the timing of cell divisions (Additional file 4: Figure S2A). We quantified this variability by plotting the time of birth of individual cells and calculating the magnitude of variation across embryos (Additional file 4: Figure S2A). The timing is fairly consistent until 9 hpa and homologous cells in different embryos divide less than 20 min apart from each other (Additional file 4: Figure S2B). We also found that embryos wt1, wt3, and wt4 have similar timing, with cell divisions occurring within 10 min of each other (Additional file 4: Figure S2C). These data suggest the development of M. membranacea is highly stereotypical with consistently timed cell divisions between individuals.
Within a single embryo, the cell divisions between the correspondent blastomeres of each quartet are mostly synchronous up to the 64-cell stage at 11 hpa (Fig. 4). At this point, we observe the first significant asynchronies in the cell divisions of a quartet, both occurring in the posterior D quadrant. The cell 1de 11 divides 2 h later than its quartet correspondents 1ae 11, 1be 11, and 1ce 11 (Fig. 4 and Additional file 5: Video S3), while the cell 1di 12 divides approximately 1 h before its partners 1ai 12, 1bi 12, and 1ci 12 (Additional file 6: Figure S3). We also observe, in the four embryos, a 3.5 h delay in the division of 3D, relative to the divisions of 3A–3C. These D quadrant asynchrony events occur with surprising consistency between the different bryozoan embryos at least until the stages analyzed in this study. We could detect a few cases of variability in the timing of divisions, but overall our data indicates the development of M. membranacea varies little between individuals.
Additional file 5: Video S3. Delayed D quadrant cell division in M. membranacea. Development of embryo wt1 (left) and wt3 (right) annotated with cell models from Simi BioCell . The lineage of the quartet 1qe 11 is highlighted in color from the 4-cell stage (3 hpa). The cell 1de 11 (red) divides 2 h later than its quartet correspondents 1ae 11 (green), 1be 11 (blue), and 1ce 11 (purple). This delayed division is one of the first morphological manifestations of the break in biradial symmetry and the establishment of the anteroposterior axis of the embryo. (MP4 19251 kb)
Finally, at a similar time point, we observe the first difference in the orientation of the cleavage plane between quartet cells. While 1de 12 divides equatorially, 1ae 12, 1be 12, and 1ce 12 divide meridionally. The asynchrony in the D quadrant and shift in cleavage orientation are the first morphological events that mark the break in the biradial symmetry of the bryozoan embryo.
Cellular origin of larval tissues
The larval body of M. membranacea develops from the four quadrants in a symmetrical manner, each lineage contributing almost equally to the structures on their respective sides: D = posterior, C = right, B = anterior, and A = left (Fig. 5a and Additional file 7: Video S4).
Additional file 7: Video S4. Cells tracked in the embryo wt1 color-coded by quadrant of origin. Video is the same as Additional file 1: Video S1 but annotated with a maximum projection of the cell models from MaMuT . Cells out of the focus plane are also shown. (MP4 19251 kb)
Progeny of the first quartet of animal blastomeres (1a–1d) gives rise to animal ectodermal structures such as the apical organ, the aboral epithelia and the corona (Fig. 5a). The apical organ is derived from derivatives of the apical-most cells 1a1, 1c1, and 1d1 (Figs. 5b and 6a). Cells 1a and 1c form the lateral and anterior most portion of the apical organ while the posterior cell 1d contributes not only to the posterior portion, but also to the tissues at the base of the apical organ (Fig. 6a). Thus, the cell 1b is the only blastomere of the first animal quartet that does not contribute to the apical organ. Epithelial cells between the apical organ and the corona are mostly derived from the octets 1q11 and 1q12. Outer coronal cells originate from 1q12 and 1q2 while inner coronal cells (turned inwards after the invagination of the vegetal plate) are derived from 1q2 (Fig. 6a). For a detailed overview of cell fates see Additional file 8: Figure S4.
The vegetal blastomeres 1A–1D form the epithelium of the vestibule, the oral/anal ectoderm, as well as the cells internalized during gastrulation, which originate the endoderm and mesoderm of the cyphonautes larva (Fig. 5a and Additional file 9: Video S5). The cellular arrangement at the vegetal plate in a 32-cell embryo (7.5 hpa) consists of 12 outer cells (3a–3d, 2aR/L–2dR/L) and four large inner blastomeres (3A–3D) (Fig. 6b and Additional file 10: Figure S5B). Here, we define gastrulation as the internalization of these four vegetal cells. It occurs by delamination and epiboly in two rounds of division of the outer vegetal 12-tets, which divide radially, pushing the four larger blastomeres internally and outlining a blastopore (Fig. 6b and Additional file 10: Figure S5B–F). At the 90-cell stage (9.7 hpa), 12 cells define the blastopore lip, but this number gets reduced to 8 cells after the next division (Fig. 6b and Additional file 10: Figure S5F–I). From the 12 vegetal cells, one does not divide (3c2) and continues to line the right side of the blastopore lip (Fig. 6b and Additional file 10: Figure S5E). Cells at the vertices of the blastopore at the 90-cell stage at 9.7 hpa (2aL2, 2bR2, 2cL2, and 2dR2) get pushed away from the blastopore lip, which now consists of 8 cells (Fig. 6b and Additional file 10: Figure S5I). Blastomeres not forming the blastoporal lip also undergo the same round of two radial divisions except for 3c1, the sister of 3c2. The derivatives of these 12 vegetal outer blastomeres form the whole ectoderm that invaginates and develops into the epithelia of the vestibule and preoral funnel. Thus, in the course of the invagination of the vegetal plate, and of the animal–vegetal elongation of the embryo, the blastopore in M. membranacea becomes the larval mouth.
Additional file 9: Video S5. Cells tracked in the vegetal ectoderm of embryo wt2 color-coded by generation. Video is the same as Additional file 1: Video S1 but the embryo is rotated 180°. The angle of the embryo was also slightly adjusted in the BigDataViewer  to align the vegetal surface and only a few slices of the z-axis are shown. Cell models were exported from MaMuT . Video shows the internalization of the large vegetal blastomeres (black) and the rounds of cell divisions of the vegetal 12-tets (blue). See Additional file 10: Figure S5 for a detailed view of the events in the vegetal ectoderm. (MP4 19148 kb)
During epiboly (9.7 hpa, around 90 cells), three of the internalized large blastomeres (3A–3C) undergo a round of unequal division forming the basal cells 4a–4c and the apical cells 1A–1C (Fig. 6c). The division of the 3D cell occurs with a 3.5-h delay in comparison to the other blastomeres (13.1 hpa). This round of division sets apart the endoderm (4A–4D) from the mesodermal tissues (4a–4c) of M. membranacea cyphonautes larva. The cell 4d is also formed, but we could not resolve its fate. The cells 4a and 4c divide twice anteroposteriorly, forming a pair of lateral rows of mesodermal cells (Fig. 6c). The most anterior cells (4aA and 4cA) form a bilateral pair of muscle cells extending from the corona to the apical organ (Fig. 7). Interestingly, one anterolateral cell (4aA1) migrates from the corona level until the apical organ during animal–vegetal elongation of the embryo (Additional file 11: Video S6). At the frontal portion of the larva, the cell 4b divides forming a column of cells stacking from the corona until the apical organ; the identity or role of these cells is unknown (Fig. 6c). Blastomeres 4A and 4C undergo anteroposterior divisions while 4B divides meridionally at 15.2 hpa lining up with the blastoporal opening and forming the endodermal tissues of the cyphonautes larva (Figs. 6c and 7).
Additional file 11: Video S6. Migration of a mesodermal cell in the embryo wt1. The anterolateral cell 4aA1 (dashed circle) migrates from the vegetal region near the corona to the animal region near the apical organ, during animal–vegetal elongation of the embryo. This was the only cell migration event we were able to track. The video was exported using MaMuT’s “Extract track stack” action. (MP4 11059 kb)
Previous work revealed that the MAPK signaling pathway might establish the position of the dorsal organizer in molluscan embryos . So far, all investigated molluscs show the asymmetric activation of MAPK in the 3D blastomere [54–57]. Using an antibody against the activated form of MAPK, we found that, in the bryozoan M. membranacea, the first detectable MAPK activity occurs in the 3D vegetal blastomere at the 28-cell stage (6.7 hpa) (Fig. 8b). MAPK activity persists in the 3D cell from the 28-cell to the 90-cell stage (9.7 hpa) and fades prior to the 3D division around 90-cell stage (Fig. 8b–f). MAPK activity is not continued in the progeny of 3D, 4D, or 4d (Fig. 8g–h) and was not detected in later embryonic stages.
The inhibition of the MAPK pathway in molluscs causes defects in the dorsoventral patterning [54–57] while, in annelids, MAPK-inhibited embryos have disorganized muscle and nerve tracts and overall shortened morphology [58–60]. We used the MEK inhibitor U0126 to investigate the role of the MAPK pathway in the development of M. membranacea at 10 °C.
We investigated the effects of different U0126 concentrations (1, 10, 25 μM) on the development of M. membranacea when applied at 2 hpa in the 2-cell stage (Additional file 12: Figure S6A). We found the severity of the phenotype correlates with the concentration of the inhibitor, where the higher concentrations of 10 and 25 μM result in the complete disruption of the normal morphology (Additional file 12: Figure S6A). These embryos show no identifiable larval structures, such as a differentiated apical organ or musculature, have a lower number of cells, and are shorter compared to control samples (Additional file 12: Figure S6B).
The proportion of embryos exhibiting a severe phenotype decreases when the treatment begins at later developmental stages (from 4–8 hpa, 8-cell), even though these time points precede the observed period of MAPK activity of M. membranacea (Additional file 13: Figure S7). Embryos treated from 10 hpa onwards show progressively milder phenotypes (Additional file 13: Figure S7). In treatments beginning at 10–16 hpa (16- to 90-cell), the larval structures, such as apical organ, ciliated band, and gut, are formed but the embryos are shorter and delayed in development in comparison to control embryos, while 18–24 hpa (early/mid gastrula) samples have almost normal morphology (Additional file 13: Figure S7).
Finally, to identify the developmental defects caused by the MEK inhibitor, we recorded M. membranacea embryos treated with 10 μM U0126 under the 4D microscope. We found the earliest abnormality associated with U0126-treated embryos is a misguided fourth cleavage (8–16 cell stage) in individuals exhibiting the severe phenotype, while embryos with milder phenotypes develop slower when compared to wild type, but do not show any obvious cleavage abnormalities (Additional file 14: Figure S8).
In order to complement the cell lineage data, we cloned the bryozoan orthologs of 16 widely conserved patterning genes that have known and largely consistent developmental roles in diverse metazoans, namely the anterior markers six3/6, dlx, otx, pax6, and nk2.1; the foregut genes foxa and gsc; the germline marker nanos, the posterior and hindgut genes bra, cdx, evx, and wnt1; the endodermal marker gata456; and the mesodermal genes twist, foxc, and foxf. We analyzed the expression of these genes during M. membranacea development (Figs. 9, 10 and 11; Additional file 15: Figure S9), and used them as markers to further reveal the identities of the different blastomeres in the bryozoan embryo.
We first detected transcripts of six3/6 – a transcription factor associated to anterior ectodermal patterning in cnidarians , hemichordates , and diverse protostomes  – during early M. membranacea gastrulation in one outer lateral vegetal plate cell (2cR2), one anterior endomesodermal cell (4B), and in five cells of the apical disc (Figs. 9 and 11a). Expression of six3/6 clears from 2cR2 and 4B, but persists in the inner cells of the forming apical organ, a central neural region occupied by serotonergic-positive cells in other cyphonautes larvae . We detected dlx transcripts, a gene involved in neurogenesis and proximodistal patterning of flies and vertebrates , in the eight animal pole cells (1q) of the 16-cell stage (5.2 hpa), broadly in the apical disc during gastrulation and elongation and, finally, localized to the whole apical organ in the late gastrula (Fig. 9).
The gene otx is involved in anterior ectodermal patterning of diverse metazoans [62, 63, 66–71] and endomesoderm specification of deuterostomes [72–74]. In M. membranacea, otx is expressed in all blastomeres between the 2- and 8-cell stages and gets restricted to the apical octet of the 16-cell stage at 5.2 hpa (Fig. 9). At the 32-cell stage (7.5 hpa), otx transcripts localize to the 1q2 octet and during gastrulation there are three rows of cells expressing otx with a posterior gap (Fig. 11b). During mid gastrula (16 hpa), two cells in the apical organ express transcripts of otx (Fig. 11b). In the late gastrula (24 hpa), otx is expressed in the corona and vestibule epithelium (Fig. 11b). Expression of pax6 is first detected during gastrulation, in bilateral patches of the animal ectoderm, and remains as a thin line of expression encircling the embryo above the corona (Fig. 9). The gene nk2.1 is involved in the patterning of the neural plate in vertebrates  and is expressed in anterior and ventral territories, including the apical/neural plate and anterior endoderm in cephalochordates , hemichordates , echinoderms , and annelids . Transcripts of nk2.1 are present in the progeny of the vegetal cells 2b and 3b in the early gastrula stage at 9.7 hpa (Fig. 9). These cells occupy an anterior vegetal position abutting the anterior blastopore lip until the edge of the vegetal plate. After the invagination of the vegetal plate, nk2.1-positive cells are lining the anterior portion of the preoral funnel, next to the mouth.
Expression of foxa is related to endoderm specification and commonly associated with the blastopore lip and foregut tissues in echinoderms  and annelids [79, 80]. At the 16-cell stage (5.2 hpa), we detected faint expression of foxa in the outer vegetal blastomeres and in 10 (out of 12) cells surrounding the four large blastomeres at the 32-cell stage (2q and 3q, except posterior cells 2dL and 2dR) (Figs. 9 and 11c). Expression of foxa persisted in the daughter cells of the next division forming two rows of cells around the blastopore with a gap at the posterior end (Fig. 11c). With the invagination of the vegetal plate, this region occupies an anterior/lateral position in the vestibule wall, surrounding the mouth region of the late gastrula (Figs. 9 and 11d). We only found transcripts of gsc at the early gastrula stage (9.7 hpa) in two anterior and a bilateral pair of cells at the vegetal plate (Figs. 9 and 11e). In the late gastrula (24 hpa), gsc is expressed in bilateral domains of the vestibule wall, which fuse anteriorly.
The widely conserved germline marker nanos [81, 82] is expressed at 7.5 hpa in two posterior cells of the vegetal plate at the 32-cell stage (2dL and 3d) (Fig. 9). In subsequent stages, nanos continues to be restricted to two cells at the posterior portion of the vegetal plate, localizing to the internal sac region of the cyphonautes larva (Fig. 11f).
The posterior/hindgut and mesodermal markers that we tested only initiate expression during gastrulation. The gene bra can have multiple roles, but it is generally related to mesoderm and posterior/hindgut patterning in several metazoans . Expression of M. membranacea bra at 9.7 hpa in the early gastrula occurs at the vegetal plate in a posterior band of cells near the blastopore lip (Figs. 10 and 11g). It localizes to 6–8 cells at the posterior end of the mid gastrula and a broad portion of the posterior and lateral vestibule ectoderm (Fig. 11g). M. membranacea bra expression domain reaches the posterior portion of the preoral funnel as well as the future hindgut area of the larva (Fig. 10). A single posterior vegetal plate cell (2dR2) and its daughter cells (2dR21 and 2dR22) express the posterior/hindgut markers cdx and evx at the early gastrula (Figs. 10 and 11h, j). During gastrulation, cdx and evx continue to be expressed at the posterior edge of the vegetal plate (Figs. 10 and 11j, l) and localize to the posterior vestibule ectoderm (hindgut) of the late gastrula (Figs. 10 and 11l). At this stage, evx is also found in the posterior region of the gut (Fig. 11l, m). We also detected a transient evx expression in the two internalized blastomeres 4a and 4c of the early gastrula. Finally, wnt1 is expressed in a row of 3–5 cells (including 2dL2, 2dR2, and 3d2) posterior to the blastopore during gastrulation (Figs. 10 and 11n). At the late gastrula (24 hpa), wnt1 is detected at the posterior-most vestibule ectoderm, positioned between the corona and hindgut (Figs. 10 and 11n).
Expression of twist, a central regulator in mesoderm differentiation in several metazoans , occurs in a narrow time window in the early gastrula of M. membranacea. We detected a colorimetric signal in bilateral internalized cells of the early gastrula – possibly 4a, 4c, or derivatives – as well as at the anterior end of the late gastrula (Figs. 10 and 11o). Transcripts of foxc, commonly expressed in anterior and posterior mesodermal domains in flies , annelids , and brachiopods , are present in one unidentified posterior vegetal plate cell of the early gastrula and two similarly positioned cells during mid gastrulation (Fig. 10). In the late gastrula, foxc expression is located in the internal sac area. The gene foxf is a transcription factor involved in mesoderm patterning and expressed mainly in visceral and anterior territories in flies [88, 89], cephalochordates [86, 90], and brachiopods . In M. membranacea it is expressed in the mesodermal cell 4b in the early and mid gastrula stages (Fig. 10). This cell and its descendants divide subsequently from basal to apical, forming a distinct frontal row of mesodermal cells expressing foxf at the anterior portion of the late gastrula.
We found two copies of the endomesodermal marker gata456  in the transcriptome of M. membranacea. While the gene gata456a is not expressed at detectable levels in any of the analyzed stages, gata456b is strongly expressed in endodermal tissues of the bryozoan. The expression of gata456b initiates early, in the vegetal 3D blastomere at the 32-cell stage (7.5 hpa) (Figs. 10 and 11p). The expression expands to adjacent lateral blastomeres 4A and 4C in the early gastrula, and in subsequent stages gata456b continues to be expressed in the endodermal tissues forming the gut of the cyphonautes larva (Figs. 10 and 11p).
The phylogenetic position of bryozoans provides a valuable opportunity to investigate the evolution of developmental traits. Even though the kinship of Bryozoa remains inconclusive – the group is more closely related either to Entoprocta and Cycliophora  or to Phoronida and Brachiopoda , or both  – most phylogenetic analyses place the bryozoans nested within the Spiralia [27, 29–32]. This indicates that the ancestral cleavage pattern of Spiralia – spiral cleavage – must have been modified in the bryozoan lineage during evolution . Here, we examined the similarities and differences between the embryogeneses of the bryozoan M. membranacea and those of spiral-cleaving embryos, by integrating cell lineage and molecular data, and provide a hypothesis for the evolution of bryozoan development from a spiral-cleaving ancestor.
Specification of the D quadrant
One critical event of animal embryogenesis is the establishment of dorsoventral polarity. In spiral-cleaving embryos, this event is tied to the specification of the D quadrant during development . In species where the first two embryonic cell divisions are unequal, the D quadrant is determined early by the asymmetric distribution of maternal cytoplasmic determinants, while in spiral-cleaving species that form equal-sized blastomeres at the 4-cell stage, the D quadrant is specified around the 24- to 32-cell stage by inductive interactions mediated by cell contacts between micromeres and macromeres [93–97]. In the current work, we found evidence that the specification of the D quadrant in the equal, biradial-cleaving bryozoan M. membranacea resembles that of equal, spiral-cleaving molluscs in the timing of specification, pattern of MAPK activation, and asynchrony of the D quadrant cell divisions post-specification.
In equal-cleaving molluscs, the specification of the D quadrant correlates with the activation of the MAPK pathway in the 3D macromere only [55, 57]. In M. membranacea, whose equal-sized blastomeres at the 4-cell stage give rise to perfectly symmetrical embryonic quadrants, that are indistinguishable from each other until gastrulation, the earliest molecular asymmetry we could detect is the activation of the MAPK pathway in a single vegetal blastomere that produces the posterior portion of the larval body. As in equal-cleaving molluscs, MAPK is activated in the bryozoan 3D blastomere on the fifth round of cell divisions, suggesting the D quadrant of M. membranacea is specified as early as the 28-cell stage. This might indicate that bryozoans and equal-cleaving molluscs undergo similar developmental mechanisms of D quadrant specification (but see below). Interestingly, most equal-cleaving spiralians studied so far exhibit a single MAPK-activated blastomere during early development, while unequal-cleaving species show diverse patterns of activation (see Additional file 16: Table S1), thus suggesting that this pattern of MAPK activity is a common feature of equal-cleaving embryogenesis independent of its cleavage geometry.
Blocking the MAPK pathway during mollusc embryogenesis results in radialized larvae that lack muscles, shell, and foot, suggesting that MAPK activation might signal for the specification of the D quadrant [54–57]. We tested if the MAPK pathway could have a similar developmental role in the bryozoan by using the MEK inhibitor U0126 . We found that blocking the MAPK pathway early in development severely disrupts the normal development of M. membranacea, a phenotype analogous to the radialized larvae of molluscs. However, later treatments do not result in axial defects, even if the drug is applied before the activation of the MAPK in the 3D blastomere. Thus, the occurrence of the severe phenotype does not correlate with the period of MAPK activation in M. membranacea, in contrast to the snail Crepidula fornicata, where the embryos become radialized in all treatments before and during MAPK activation – but are not disrupted if the drug is applied after this critical period .
The fact that blocking the 3D MAPK activation does not lead to axial defects in the bryozoan suggests that the MAPK pathway might not have a role in specifying the D quadrant in M. membranacea. However, our experimental dataset does not exclude alternative explanations. One possibility is that the U0126 concentration that we used for the timed experiments (10 μM) is not sufficient to completely inactivate MEK, and the remaining MAPK (ERK) – although undetected by immunohistochemistry – would still transduce the signal and form larvae without axial defects. For instance, similar mildly abnormal phenotypes were also observed in molluscs treated with the same U0126 concentration, which can indicate the partial inactivation of the MAPK signaling [57, 99]. In addition, we cannot exclude the possibility that the early disruption we observe in M. membranacea is due to an undetected period of MAPK activation before the 3D activation, or even due to off-target effects of the inhibitor. For these reasons, although we present preliminary evidence that MAPK inhibition alters M. membranacea development, the developmental role of the MAPK pathway in bryozoans remains unclear and dependent upon future work using higher U0126 concentrations and additional experimental conditions.
Once the D quadrant has been determined, it typically shows asynchronous cell divisions in relation to the other quadrants of spiral-cleaving embryos . For instance, the 3D macromere in the mollusc Patella vulgata  and the 1d derivatives of Ilyanassa obsoleta [101, 102] undergo a late division. Our analyses of M. membranacea cell lineage indicate similar asynchronous cell divisions in the D quadrant, which include the 3D blastomere and 1d derivatives of the bryozoan. Therefore, the specification of the D quadrant seems to be correlated with subsequent changes in the cell cycle timing in both M. membranacea and spiral-cleaving embryos.
Overall, M. membranacea exhibits a similar pattern and timing of MAPK activation, as well as equivalent asynchronous cell divisions in the D quadrant, when compared to equal-cleaving molluscs. Given the phylogenetic position of bryozoans, these similarities might suggest that some of the underlying traits of spiral-cleaving embryos were maintained during the evolutionary transition from spiral to biradial cleavage. The comparison also reveals that equal cleavage might be associated with a single D quadrant MAPK-activated blastomere in spiralian development. Nevertheless, the MAPK pathway is still poorly sampled in spiralians, and other spiral and non-spiral-cleaving groups, such as phoronids, nemerteans, polyclads, rotifers, and gastrotrichs, need to be investigated to properly understand the roles and the evolution of MAPK signaling in spiralian development.
Comparative spiralian fate maps
The stereotypic nature of spiral cleavage supports the identification of putative homologous blastomeres between different spiralian lineages, and therefore enables the unprecedented comparison of blastomere fates across clades [6, 9, 12, 103, 104]. The comparative study of spiral cleavage has revealed that homologous blastomeres share mostly-similar fates in various clades [5, 18, 19, 48, 105]. The cleavage of M. membranacea clearly differs from the spiral cleavage pattern, which complicates the identification of homologous blastomeres between the bryozoan and a spiral-cleaving embryo. However, we established a common developmental feature to base our comparative cell lineage and gene expression analyses.
In both spiral and bryozoan embryogenesis, the vegetal blastomeres sequentially give rise to quartets of daughter cells, while remaining at the vegetal-most portion of the embryo until being internalized during gastrulation. We thus compare the quartets of M. membranacea to the quartets of spiral-cleaving embryos in terms of gene expression and fate in the larval tissues. We find the quartets have a similar molecular identity and contribute to the same set of structures in the larvae of bryozoan and spiral-cleaving groups, and that the subset of blastomeres that gives rise to these structures partially overlaps (Fig. 12). This indicates that bryozoans might share a common embryonic patterning of early blastomere fates with other spiralians, and that, in the current phylogenetic scenario, such developmental trait has remained conserved despite the drastic modification in the cleavage pattern from spiral to biradial.
First quartet: apical organ and ciliated band
The first quartet of micromeres in spiral-cleaving embryos contributes to the apical organ, the ciliated band, and all the ectoderm in between [6, 9, 11, 106–114]. In the M. membranacea, the first quartet of animal blastomeres also gives rise to these ectodermal structures of the cyphonautes larva. This suggests that, in both the bryozoan and spiral-cleaving embryos, the third cleavage demarcates a split in the embryonic fate map, in which the first quartet of animal blastomeres only gives rise to the ectodermal structures placed towards the animal pole, while the progeny of the vegetal blastomeres (i.e., the second, third, and fourth quartets) produces a different set of larval structures (see the next sections). When we compare in more detail the specific fates of the descendants of the first quartet, we find that some blastomeres contributing to the apical organ or ciliated band of the cyphonautes larva indeed contribute to the respective structures of spiral-cleaving larvae – but that this similarity is not complete, and different subsets of blastomeres contribute to the apical organ and ciliated band.
The apical organ, for example, is usually formed by the progeny of the apical-most 1q1 micromeres in groups with spiral cleavage [48, 105]. While the apical organ of M. membranacea larva is also derived from the apical-most subset of 1qi 1 (=1q11), descendants of 1qe 1 (=1q21) also contribute to the structure. We find a similar situation when comparing the embryonic origin of the corona (i.e., the ciliated band of the cyphonautes larva) with the prototroch – a ciliated band considered to be an ancestral trait for the larval stages of trochozoan spiralians . The prototroch of annelids and molluscs is formed by 1q (accessory and primary trochoblasts) and 2a–c (secondary trochoblasts) descendants [106, 113, 116]. Our data reveals that the corona of M. membranacea is formed by blastomeres equivalent to the accessory and primary trochoblasts of the prototroch , but unlike spiral-cleaving embryos, the second quartet does not contribute to the ciliated band of the bryozoan larva (Additional file 17: Figure S10). Thus, despite being located in a vegetal domain of the bryozoan embryo, such as the telotroch of some spiral-cleaving larvae, the corona shares a common embryonic origin with the prototroch.
In general, we find that equivalent early blastomeres of the bryozoan and spiral-cleaving embryos contribute to similar larval structures, but that the fate of the progeny of these early blastomeres only partially overlap between the bryozoan and spiral-cleaving embryos. These observations suggest that, during the evolution of the bryozoans, shifts in the blastomere fates occurred in late embryogenesis while the early embryonic patterning, presumably inherited from a spiral cleavage ancestor, might have remained conserved.
The bryozoan fate map indicates that the apical organ, outer ectoderm, and corona of the cyphonautes larva have a similar embryonic origin as the apical organ, pretrochal elements, and prototroch of spiral-cleaving embryos, respectively. They all derive from the first quartet blastomeres. The fate map similarity is paralleled by molecular data, since the genes expressed in this region of M. membranacea have an equivalent spatial arrangement in other spiralian embryos. six3/6 and dlx are expressed at the animal end, while otx is expressed in the vegetal-most progeny of the first quartet. Other spiralian embryos display a similar arrangement of these transcripts [63, 68, 69, 71, 117–122] (Additional file 18: Table S2). Therefore, the region derived from the first quartet of M. membranacea match the pretroch region of spiral-cleaving embryos.
In M. membranacea the blastomeres that form the apical organ express six3/6 and dlx from the 16-cell stage, suggesting that these genes might be involved in the establishment of the embryonic animal–vegetal identities, and possibly in the molecular patterning of the cyphonautes apical organ. The expression of otx in the bryozoan is associated with the ciliated band of the cyphonautes larva, similar to other spiralians where otx is expressed near or in the larval ciliated band [68, 69, 123]. The gene is an interesting example because it provides the opportunity to integrate cell lineage and gene expression data between the bryozoan and spiral-cleaving embryos. As explained above, the ciliated band of trochophore larvae, the prototroch, is formed by the contribution of first quartet and second quartet blastomeres [106, 116], while the corona of M. membranacea derives solely from first quartet blastomeres, which are putatively homologous to the primary trochoblasts of the prototroch. In the mollusc Patella vulgata, otx is expressed in all prototroch cells . Interestingly, the second quartet blastomeres of M. membranacea – the set of blastomeres that form the secondary trochoblasts in the prototroch – also express otx, as observed in P. vulgata, even though these cells do not contribute to the corona of the cyphonautes larva. This observation suggests that presumptive homologous blastomeres between the bryozoan and other spiralians might still share a similar molecular identity, even though they do not form similar tissues.
Overall, our work reveals that the first quartet of M. membranacea embryo and the first quartet of spiral-cleaving embryos give rise to a similar set of larval structures, and give rise to a larval body region with similar molecular profile. Thus, the outer ectodermal region of the cyphonautes larva corresponds, in developmental terms, to the head region of other spiralians.
Second and third quartet: larval mouth
The second and third blastomere quartets of spiral-cleaving embryos contribute to a diverse set of ectodermal structures, such as the foregut, ciliated bands, neurons, the mollusc shell gland and foot, the annelid trunk and nerve cords, as well as ecto-mesodermal muscle cells [107–114, 124, 125]. In M. membranacea, these blastomeres form the whole vegetal ectoderm that gives rise to the vestibule epithelium, including the preoral funnel and posterior ectoderm of the cyphonautes larva. In most spiralians, the second and third quartets are the blastomeres surrounding the blastopore – the orifice formed at the site of endomesoderm internalization , whose developmental fate has been a significant trait for the discussions about metazoan evolution [127–129]. Nevertheless, the fate of the blastopore in bryozoans remained open to discussion [39, 41, 44, 49].
Even though in most gymnolaemate bryozoans the blastopore closes after gastrulation [41, 43], or in some cases, an orifice is not formed at all , an ultrastructural study in M. membranacea revealed that its blastopore remains open until the larval stage . Our cell lineage data indicate that cells at the blastopore lip give rise to the preoral funnel of M. membranacea, and that the endodermal cells lining the blastopore form the anterior portion of the larval gut. We also found the foregut marker foxa is expressed around most of the blastopore lip, except for a couple of posterior rows, and that foxa expression persists around the future larval mouth, indicating that most cells associated with the blastopore of M. membranacea have a foregut molecular identity. Thus, independent ultrastructural, molecular, and cell lineage data provides robust evidence for a persistent blastopore and the protostomic development of M. membranacea, as previously suggested .
We found that the vegetal ectoderm – the cells derived from the second and third quartet – exhibits an anteroposterior polarity, as revealed by the differential expression of molecular markers. The anterior/foregut markers nk2.1, foxa, and gsc are expressed in a region opposed to posterior/hindgut markers bra, cdx, evx, and wnt1, which are generally restricted to the D quadrant. Transcripts of nk2.1 are restricted to the B quadrant in a comparable position, in relation to the cyphonautes anteroposterior axis, to the anterior/ventral expression found in other bilaterians [62, 71, 122, 129, 130]. In a similar fashion, transcripts of wnt1, a gene commonly expressed at the posterior end of bilaterians [131–134], occur at the posterior region of the vegetal ectoderm of M. membranacea (see also Additional file 18: Table S2 for a comparison of spiralian-specific gene expression patterns). This suggests that at least some molecular aspects of the bilaterian axial patterning has remained conserved in the cyphonautes larva.
In some cases, the transcripts of M. membranacea are not only located at similar body regions, but in the putative homologous blastomeres of spiral-cleaving embryos. An example is the expression of foxa between the bryozoan and the annelid Hydroides elegans . In both, foxa is expressed in the second quartet blastomeres early in development and in the cells that surround the blastopore during gastrulation, with a peculiar posterior gap . Another comparable cellular expression is the gene bra, expressed in the second and third quartet progeny at the posterior lip of the blastopore of the molluscs Patella vulgata  and Haliotis asinina . Therefore, the M. membranacea data indicates that the molecular identity of the blastomeres remained conserved to a certain extent within spiralians, despite the modified cleavage geometry and vegetal placement of the second and third quartet in the bryozoan embryo.
Fourth quartet: muscle and mesenchymal cells
The embryonic source of mesoderm in bryozoans has been a contentious topic . Classical works suggest that mesodermal cells derive from endodermal blastomeres, but could not demonstrate the embryonic origin with cellular resolution [40–43, 45, 135]. However, recent ultrastructural data in M. membranacea suggests an ectodermal origin for the bryozoan mesoderm, from the delamination of an ectodermal cell during gastrulation . Our cell lineage data indicate that the first mesodermal cells of M. membranacea derive from the fourth quartet (4a–4c). The lateral cells 4aA and 4cA form the anterior muscles of the cyphonautes larva, while the progeny of 4b1 gives rise to a stack of mesenchymal cells that express the anterior mesoderm marker foxf. We did not observe the delamination of an anterior ectodermal cell as described by Gruhl , but cannot discard the existence of other cells contributing to the mesoderm of M. membranacea. Our work corroborates previous classical studies of bryozoan embryology by revealing that the mesoderm of the bryozoan M. membranacea originates from multiple fourth-quartet blastomeres.
The source of mesodermal tissues in spiral-cleaving embryos has been extensively studied and discussed [5, 9, 60, 108, 109, 112–114, 124, 125, 136–141]. There are generally two sources, an anterior mesoderm derived from the third quartet blastomeres (known as ectomesoderm) and a posterior mesoderm derived from the 4d blastomere. While the blastomeres contributing to the anterior mesoderm (usually 3a and 3b) are more variable [113, 141], most spiral-cleaving species have the 4d as the sole endomesodermal contributor. However, there are exceptions. In the annelid Capitella teleta, for example, mesoderm formation shifted to the 3c and 3d blastomeres [10, 114], even though the 4d blastomere still produces the germline . This reveals how a particular contribution to later tissues can be decoupled and geared to different blastomeres in spiral-cleaving embryos . We find that the source of mesodermal tissues of the bryozoan M. membranacea differs from other spiralians because (1) the third quartet does not contribute to the anterior mesoderm and (2) multiple blastomeres of the fourth quartet give rise to mesodermal tissues. In addition, the blastomeres 4a–4c often give rise to endodermal tissues in spiral-cleaving embryos [108, 113]. Therefore, our data suggests that the specification of the anterior mesoderm in the bryozoan might have been shifted in time and in position from the third to the fourth quartet.
Although only a subset of the genes we analyzed is expressed in the mesoderm of M. membranacea (six3/6, evx, twist, foxf), the patterns indicate that the bryozoan mesoderm is already regionalized at early gastrulation. For instance, we found that 4d expresses evx and cdx, genes commonly associated to posterior mesodermal and hindgut fates in spiralians [120, 121, 129, 143–149]. Even so, we could not resolve the fate of the 4d blastomere in M. membranacea although the expression data indicates that the 4d blastomere might contribute to the posterior mesoderm and hindgut of the bryozoan. The lateral mesodermal cells 4a/4c and derivatives express evx, and possibly twist, in the same manner as the expression of evx orthologs in the annelid Capitella teleta during early development . Finally, at the anterior mesoderm, we find the expression of foxf, also observed in the brachiopod Terebratalia transversa . Overall, these molecular data reveal that M. membranacea mesoderm is regionalized and that at least some of the expression patterns are conserved with other spiralians.
The 4d cell and its descendants also form the germline and are known to express nanos in spiral-cleaving embryos . Germ cells have not been identified during the embryogenesis of any bryozoan and were only found in zooids after metamorphosis . The expression of nanos in M. membranacea differs from the pattern of spiral-cleaving embryos, since nanos is expressed in two posterior cells of the second and third quartet. These blastomeres divide repeatedly, but nanos expression is always retained in two cells that become part of the larval internal sac – the structure that persists during metamorphosis giving rise to the outer case of the zooid . Thus, we hypothesize the nanos-positive cells might be stem cells contributing to the differentiation of the internal sac, but further analysis in competent larvae and metamorphosed juveniles are needed to clarify the fate and molecular identity of these cells.
4Q blastomeres: gut
The 4Q blastomeres are the largest and yolkier cells of the gastrulating bryozoan embryo and, in M. membranacea, they produce the cyphonautes larval gut. Similarly, the macromeres of spiral-cleaving embryos generally give rise to endodermal tissues [6, 9, 11, 103, 108, 110–114, 153], even though in the polyclad Hoploplana inquilina these cells break up and the gut is likely derived from 4d [16, 109]. Nevertheless, the endodermal fate of these large vegetal cells appears to be a common feature between M. membranacea and spiral-cleaving embryos.
The expression of the endodermal marker gata456b corroborates the molecular identity of the 4Q blastomeres in the bryozoan. In spiralians, gata456 expression is mainly associated to endodermal and, in some cases, mesodermal tissues [80, 87, 121, 129, 154–156]. The expression of gata456b in M. membranacea is clearly associated with the larval gut, suggesting that the molecular patterning of endodermal structures in the bryozoan is conserved with other spiralian groups, and possibly to other bilaterians . One exception is the bryozoan Bugula neritina, where gata456 is expressed at the apical organ . However, this might be related to the fact that the coronate larva of B. neritina does not have a gut, and requires further investigation. The early expression of gata456b in M. membranacea does differ slightly from other spiralians, since neither the 4B blastomere nor the 4a–4c blastomeres show gata456 transcripts [80, 154, 156]. Nevertheless, the common fate and gata456 expression between M. membranacea 4Q blastomeres (4A, 4C, and 4D) and the macromeres of other spiralians suggest that the bryozoan shares similar molecular and developmental traits for the endoderm patterning with spiral-cleaving embryos.
A modified spiral cleavage
Our investigation shows a series of similarities between the embryonic development of the bryozoan M. membranacea and the embryogenesis of annelids, molluscs, nemerteans, and polyclads. The vegetal blastomeres sequentially give rise to quartets of daughter cells, the first asynchronous cell divisions occur in the posterior quadrant, the quadrant identities can be identified at the 32-cell stage, the MAPK activity resembles that of equal-cleaving molluscs, and several genes are expressed in equivalent blastomeres or embryonic regions (Fig. 12). In addition, the early blastomeres of M. membranacea and of spiral-cleaving embryos have similar fates in the larval tissues (Fig. 12). That is, the first animal blastomeres form the whole region from the apical organ to the corona – equivalent to the pretrochal region – and the ciliated band itself, the second and third quartets contribute to the oral ectoderm, the fourth quartet gives rise to the mesoderm of the larva, and the four large vegetal blastomeres are internalized and become endoderm (Fig. 12). Since the phylogeny of Spiralia indicates spiral cleavage is ancestral and bryozoans are nested within the clade [27, 29–32], we interpret these developmental similarities as inherited traits from an ancestral spiral-cleaving embryogenesis.
In this context, during the evolution of gymnolaemate bryozoans the ancestral spiral cleavage pattern, characterized by the alternating oblique cell divisions, was modified to biradial cell divisions. While the cleavage pattern changed and the anterior mesoderm was reallocated to the fourth quartet, some aspects of the development have remained conserved, such as the D quadrant specification, MAPK activity, and overall fate map of early blastomeres.
Spiral cleavage has been modified not only in bryozoans, but in different spiralian branches as well, such as flatworms, molluscs, and brachiopods . In most of these cases, the embryonic development has changed to such extent that no traces of spiral cleavage are found, and the ancestral cleavage geometry can only be inferred by the phylogenetic position of the clade (e.g., the discoidal cleavage of cephalopods ). Remnants of spiral cleavage usually consist of oblique mitotic spindles, as in the flatworm Macrostomum lignano, which displays a typical spiral cleavage pattern until the third cleavage, when the embryonic development becomes considerably modified . The bryozoan M. membranacea differs from these previously known cases because we can recognize shared cell lineage and developmental traits with spiral-cleaving embryos that are not the cleavage geometry itself. The evolutionary mechanisms involved in the transition from spiral to biradial cleavage remain unclear, but changes in the orientation of the mitotic spindle have a genetic basis and are beginning to be uncovered using molecular and computational approaches [160–162]. Our data suggests that the quartet-divisions, cell fates, and other traits commonly associated with a spiral cleavage program were maintained in the bryozoan development despite the evolutionary modification to a biradial cleavage pattern.
Evolution of cleavage patterns
The fate of the early embryonic blastomeres is thought to be causally related to the cleavage pattern during development [2, 163]. In such case, a change in the cleavage geometry would lead to a change in the cell fates. However, the bryozoan cell lineage illustrates a case where the cleavage pattern and blastomere fates are not evolutionary coupled. We find that the fates of the early blastomeres are similar between M. membranacea and spiral-cleaving embryos despite the modified biradial cleavage pattern. The relative positioning of the second and third quartets even differs between the bryozoan and a typical spiral-cleaving embryo (Fig. 12), but these blastomeres still contribute to similar tissues, suggesting the early cell fate determination remained relatively conserved during bryozoan evolution. The bryozoan cell lineage illustrates how a widely conserved determinate cleavage pattern – spiral cleavage – can evolve without major changes in other developmental traits, such as the blastomere fates and molecular identity.
We found several developmental genes expressed in a similar spatial arrangement between bryozoans and other spiralians, as revealed by MAPK (3D blastomere), otx (vegetal-most blastomeres of the first quartet), and bra and foxa (second and third quartets). This molecular map is similar not only to the typical spiral-cleaving embryos, but also to brachiopod embryos [87, 129], whose embryos have a much greater number of cells and no stereotypic cleavage pattern . Thus, a single cell in the bryozoan embryo expressing gata456 might be homologous to a whole region of gata456 expression in the brachiopod embryo , as suggested by Hejnol . This reinforces the hypothesis that cell fate determination is not tied to a particular cleavage pattern, but depends on the underlying molecular framework established early in development . This is the case for nematodes, whose cleavage patterns diverged drastically between groups without a corresponding change in the resulting phenotype . However, a clearer parallel case to the spiral-to-biradial evolution is the transition from the spiral cleavage of polychaete annelids to the derived cleavage of clitellate annelids; despite the differences in the cleavage pattern, the clitellate fate map does not deviate significantly from the annelid ground plan . Overall, our findings support the hypothesis that, in evolutionary terms, the causal ontogenetic connection between cleavage pattern and blastomere fates, if any, can be broken . In the case of the bryozoan M. membranacea, the molecular identity and fate of the early blastomeres might have been maintained, despite the modification in the geometry of cell divisions. Further comparative cell lineage studies with other non-spiral spiralian lineages, such as gastrotrichs and rotifers, will be crucial to establish the ancestral traits of spiralian development and to better comprehend the relation between cleavage patterns and cell fates during evolution.
The embryonic development of M. membranacea provides a unique comparative standpoint to the typical spiral cleavage pattern. It reveals that spiral cleavage is not an all-or-nothing character and has been extensively modified in the diverse spiralian lineages. In particular, we suggest that the cleavage geometry of the bryozoan embryo evolved independently from other spiralian developmental traits. Therefore, modifying spiral cleavage does not require drastic developmental changes such as the ones found in cephalopods or parasitic flatworms. More generally, our data suggests that determinate cleavage patterns can be modified without major changes to the molecular identity of blastomeres and cell fates, which challenges the idea that the cleavage pattern is evolutionarily coupled to the specification of cell fates. In this perspective, the evolutionary conservation of cell fates in spiral-cleaving clades might be a consequence of a conserved underlying molecular patterning, overlaid by a determinate cleavage pattern. Overall, our work highlights the importance of comparative data to better understand the evolution of spiralian development.
Collection, spawning, and cultures
We collected M. membranacea in the fjord waters of Hjellestadosen (60°15′23.9" N 5°14′20.1" E) in Bergen, Norway, between May and September. We handpicked kelp blades with ripe bryozoan colonies from floating boat docks, and maintained the fronds in tanks with flowing sea water at 10 °C. To induce spawning, we cut a portion of the kelp blade with mature colonies, usually the ones with more opaque whitish/pinkish zooids, and transferred it to a glass bowl with sea water sterilized with UV-light and filtered through a 0.2 μm mesh (UVFSW). The bowl was placed under a stereomiscroscope with direct light and a digital thermometer to monitor the water temperature. Ripe colonies began to spawn in around 5 minutes or more, or usually when the temperature reached 15 °C. Once the temperature rose to 16 °C, the bowl was cooled down on ice with no direct light.
A spawning colony was sequentially transferred to new bowls with UVFSW at 10 °C to distribute the vast amounts of eggs. For each bowl with eggs, we added a solution of ethylenediaminetetraacetic acid (EDTA) to a final volume of 0.1 mM (usually ~20 μL of 0.5 M EDTA) to induce egg activation . The bowl was then placed in a incubator at 15 °C for 30–60 min. Activated eggs were concentrated by swirling, distributed to smaller glass bowls with UVFSW, and washed twice to remove the EDTA. We adjusted the number of eggs per bowl so that eggs sitting on the bottom of the dish no longer touched each other. We kept the cultures at 15 °C. One colony could be re-used for spawning multiple times. M. membranacea colonies maintained in the flowing tanks remained viable to developmental studies for a week.
4D recordings and cell tracing
We pipetted embryos at the 2-cell stage to a glass slide coated with poly-L-lysine. We mounted the embryos under a cover slip with supporting clay feet, completed the volume with UVFSW, and sealed the cover slip with vaseline. The slide was put under an automated 4D microscopy system (Caenotec, email@example.com)  with a cooling ring around the objective to keep the temperature at 15 °C. We recorded whole-embryo stacks (40–60 optical slices) every 40 s under differential interference contrast for four individuals (wt1–wt4). Development was recorded for approximately 24 h, until the embryos became ciliated and swam away from the field of view (Additional file 1: Video S1 and Additional file 2: Video S2).
We loaded the image sequence data into the tracking software Simi BioCell (Simi®)  and manually traced individual cells. We programmatically analyzed the cell lineage data using simi.py, a Python library we wrote for this purpose and available at https://github.com/nelas/simi.py. To conduct further analyses and visualization of the data, we used simi.py to convert the data from Simi BioCell into the format of MaMuT (Massive Multi-view Tracker) , a cell lineage tracker integrated with the BigDataViewer  for visualization of large image datasets on the Fiji open source platform .
The source cell lineage files, the processed and converted cell lineages, and the data files used for all analyses have been deposited in a data repository and are available for download .
Cell lineage nomenclature
We annotated the individual cells of M. membranacea using the spiral cleavage nomenclature [6, 9, 11] with modifications to better describe the unique features of the bryozoan cell lineage. The blastomeres at the 4-cell stage were labeled as A, B, C, and D according to their fate, as identified by the video recordings. For example, the blastomere giving rise to the dorsal/posterior region of the cyphonautes larva was assigned to the D quadrant. The quartets derived from the four large blastomeres were labeled as in spiral cleavage (1q–4q). At the 8-cell stage, the animal blastomeres were named 1a–1d and the vegetal blastomeres 1A–1D. Where applicable, the progeny of a blastomere received the standard superscripts of spiral cleavage, with 1 for the apical daughter cell and 2 for the basal daughter cell. We have adapted the nomenclature to accommodate the first cleavage of the animal blastomeres, where cells occupy the same position in the animal–vegetal axis at the 16-cell stage. The four central cells received the subscript i (for internal), and the four outer cells received the subscript e (for external). As an example, the cell 1ai 1 is the apical progeny of the internal A quadrant animal blastomere. If cells divided along the anteroposterior axis, they received the A or P superscript. In addition, we labeled the progeny of meridional divisions with the superscript R, for the cell to the right, and L, for the cell to the left when viewed from the animal pole.
We fixed representative developmental stages for antibody staining in 4% formaldehyde for 1 h at room temperature, washed the embryos in PTw (1x PBS + 0.1% Tween-20) and stored them in PTw at 4 °C. For in situ hybridization, we fixed the samples in a solution of 4% formaldehyde/0.2% glutaraldehyde solution to avoid tissue damage during the protocol. After 1 h fixation at room temperature, we washed the embryos in PTw, dehydrated them through a methanol series, and kept the samples in 100% methanol at –20 °C.
Gene cloning and in situ hybridization
We assembled the Illumina RNA-seq reads from M. membranacea (NCBI SRA Project: SRX1121923) with Trinity  and used known genes to identify putative orthologs in the transcriptome. We performed PCR using gene-specific primer pairs on cDNA synthesized with the SMARTer RACE cDNA Amplification kit (Clontech). Primers were designed with Primer3 . Gene sequences and the corresponding primer pairs were deposited in the GenBank (NCBI) with the accession numbers KY565381–KY565397, and are also available in the data repository . We synthesized antisense DIG-labeled riboprobes with MEGAscript kit (Ambion) and performed colorimetric in situ hybridization according to an established protocol . We observed no significant variability in the gene expression patterns between individual embryos.
Orthology was assigned by aligning amino acid sequences of M. membranacea against annotated genes from diverse metazoans using MAFFT 7.271 , retaining only informative portions of the alignment with GBlocks 0.91b with relaxed parameters  and running a Maximum Likelihood phylogenetic analysis with RAxML 8.2.4  using automatic model recognition and rapid bootstrap. We manually verified the alignments using UGENE . Resulting trees from the maximum likelihood analysis were rendered into cladograms using the ETE Toolkit  (Additional file 19: Figure S11). Gene orthology runs and source files are available in the data repository .
Immunohistochemistry and MAPK antibody
We permeabilized the embryos with several washes in PTx (1x PBS + 0.2% Triton X-100) for 2 h and blocked with two washes of 1 h in PTx + 0.1% bovine serum albumin (BSA), succeeded by 1 h incubation in PTx + 5% normal goat serum. Samples were incubated with the primary antibody for the MAPK diphosphorylated ERK-1&2 (Sigma M9692-200UL) diluted 1:200, and stored overnight at 4 °C on a nutator. We removed the MAPK antibody with three 5 min and four 30 min washes in PTx + 0.1% BSA, blocked in PTx + 5% normal goat serum for 1 h, and incubated nutating overnight at 4 °C with the secondary antibody Anti-Mouse-POD conjugate (Jackson) diluted 1:250. We removed the secondary antibody with three 5 min followed by several washes in PTx + 0.1% BSA for 2 h.
To amplify and develop the signal, we incubated the embryos for 3–5 min with the provided reagent solution and fluorochrome from TSA reagent kit Cy5 (Perkin Elmer). We stopped the reaction with two washes in a detergent solution (50% formamide, 2x SSC, 1% SDS) at 60 °C to reduce background, followed by PTw washes at room temperature. We stained nuclei by incubating permeabilized embryos in DAPI 1:500, Sytox Green 1:1000, or Propidium Iodide 1:500 for 2 h. Nuclei staining was combined with f-actin staining by the addition of BODIPY FL phallacidin 5 U/mL previously evaporated to remove methanol.
We repeated the immunostaining several times under different conditions to optimize the signal-to-noise ratio for different developmental stages to identify the developmental sequence of MAPK activation and using at least 100 embryos per well to account for individual variability.
To verify the identity of the single MAPK-activated blastomere – the embryo is still biradial at this stage – we carefully cross-checked cell lineage, in situ hybridization, and immunohistochemistry data. The relative timing of division between the 3Q blastomeres is highly consistent between the four M. membranacea embryos analyzed in this study. The blastomeres 3A and 3C divide first (in synchrony), the 3B divides next and, around 3.5 h later, the division of the 3D blastomere finally occurs. By analyzing MAPK immunostainings of closely timed developmental stages, we were able to stage the embryos and fully resolve identity of the 3Q blastomeres. In addition, we performed immunostaining after the in situ hybridization of nk2.1 and verified that the MAPK-activated blastomere is opposite to the nk2.1 territory (restricted to the B quadrant) (see ).
We blocked the MAPK pathway in M. membranacea by soaking the embryos with the MEK inhibitor U0126 (Promega) – a compound that inhibits the activation of the MAPK (ERK 1&2) by inhibiting the kinase activity of MAP Kinase Kinase (MAPKK or MEK 1/2) . U0126 has been extensively used to study the MAPK role in spiralian development [54–60]. We resuspended U0126 in dimethyl sulfoxide (DMSO) to a stock concentration of 10 mM, aliquoted to smaller volumes and stored the samples at –20 °C for no more than a week. We diluted the stock solution to working concentrations using UVFSW and incubated the treated embryos. Controls were incubated in UVFSW and in DMSO diluted in UVFSW 1:400, i.e., the maximum amount of DMSO used in treated wells. The MAPK experiments were conducted at 10 °C usually in two setups, a fine-picked where only 2-cell embryos were collected to assure we were seeing the inhibitor effects in healthy embryos, and a course-picked sample where a larger amount of embryos was picked.
Microscopy and image processing
We mounted in situ embryos in 70% glycerol in PTw. Embryos from antibody staining were mounted in 97% 2,2’-thiodiethanol [180, 181], 80% glycerol in PBS, or SlowFade® Gold Antifade (ThermoFisher). We imaged the samples with a Zeiss AxioCam HRc mounted on a Zeiss Axioscope A1, using differential interference contrast technique (Nomarski) for in situ hybridizations and a fluorescent lamp for the MAPK antibody staining. We used a Confocal Leica TCS SP5 to image fluorescent samples. Colorimetric in situ hybridizations were also scanned under the confocal using reflection microscopy . We processed all resulting confocal stacks in Fiji . When necessary, we adjusted the distribution of intensity levels to improve contrast with Fiji or GIMP. We created vector graphics and assembled the figure plates using Inkscape.
To convert the time-lapse image stacks from the 4D microscopy system to video, we blended the in-focus information from different focal planes (focus stacking) using the enfuse program of the Enblend/Enfuse software for each timepoint, and animated the resulting image sequence to 25 frames per second using FFmpeg. We exported the cell models generated by Simi BioCell and MaMuT and overlaid in the videos to show the position of the tracked cells in the embryo. The code for the image/video processing steps are available in our previous study .
days post activation
ethylened iami netetraacetic acid
hours post activation
phosphate-bu ffered saline
filtered sea water sterilized with UV light
Davidson EH. How embryos work: a comparative view of diverse modes of cell fate specification. Development. 1990;108:365–89.
Wray GA. The evolution of cell lineage in echinoderms. Am Zool. 1994;34:353–63. doi:10.1093/icb/34.3.353.
Raff RA. Direct-developing sea urchins and the evolutionary reorganization of early development. Bioessays. 1992;14:211–8. doi:10.1002/bies.950140403.
Schulze J, Schierenberg E. Evolution of embryonic development in nematodes. Evodevo. 2011;2:18. doi:10.1186/2041-9139-2-18.
Hejnol A. A twist in time—the evolution of spiral cleavage in the light of animal phylogeny. Integr Comp Biol. 2010;50:695–706. doi:10.1093/icb/icq103.
Wilson EB. The cell-lineage of Nereis. A contribution to the cytogeny of the annelid body. J Morphol. 1892;6:361–466. doi:10.1002/jmor.1050060301.
Lillie FR. The embryology of the Unionidae. A study in cell-lineage. J Morphol. 1895;10:1–100. doi:10.1002/jmor.1050100102.
Child CM. A preliminary account of the cleavage of Arenicola cristata, with remarks on the mosaic theory. Zool Bull. 1897;1:71–94. doi:10.2307/1535469.
Conklin EG. The embryology of Crepidula, A contribution to the cell lineage and early development of some marine gasteropods. J Morphol. 1897;13:1–226. doi:10.1002/jmor.1050130102.
Eisig H. Zur Entwicklungsgeschichte der Capitelliden. Mittheilungen Aus der Zoologischen Station Zu Neapel. 1898;13:1–292.
Child CM. The early development of Arenicola and Sternaspis. Wilhelm Roux Arch Entwickl Mech Org. 1900;9:587–723. doi:10.1007/BF02156195.
Treadwell AL. The cytogeny of Podarke obscura Verrill. J Morphol. 1901;17:399–486. doi:10.1002/jmor.1050170304.
Wilson EB. Experiments on cleavage and localization in the nemertine-egg. Wilhelm Roux Arch Entwickl Mech Org. 1903;16:411–60. doi:10.1007/BF02152027.
Nelson JA. The early development of Dinophilus: a study in cell-lineage. Proc Acad Natl Sci Phila. 1904;56:687–737. doi:10.2307/4062961.
Zeleny C. Experiments on the localization of developmental factors in the nemertine egg. J Exp Zool. 1904;1:293–329. doi:10.1002/jez.1400010205.
Surface FM. The early development of a polyclad, Planocera Inquilina Wh. Proc Acad Natl Sci Phila. 1907;59:514–59. doi:10.2307/4063205.
Costello DP, Henley C. Spiralian development: a perspective. Am Zool. 1976;16:277–91. doi:10.1093/icb/16.3.277.
Henry J, Martindale MQ. Conservation and innovation in spiralian development. Hydrobiologia. 1999;402:255–65. doi:10.1093/icb/16.3.277.
Lambert JD. Developmental patterns in spiralian embryos. Curr Biol. 2010;20:R72–7. doi:10.1016/j.cub.2009.11.041.
Maienschein J. Cell lineage, ancestral reminiscence, and the biogenetic law. J Hist Biol. 1978;11:129–58. doi:10.1007/BF00127773.
Dohle W. Zur Frage der Homologie ontogenetischer Muster. Zool Beitr (NF). 1989;32:355–89.
Dohle W. Differences in cell pattern formation in early embryology and their bearing on evolutionary changes in morphology. Geobios Mem Spec. 1989;22(Supplement 2):145–55. doi:10.1016/S0016-6995(89)80015-4.
Guralnick R. A recapitulation of the rise and fall of the cell lineage research program: the evolutionary-developmental relationship of cleavage to homology, body plans and life history. J Hist Biol. 2002;35:537–67. doi:10.1023/A:1021119112943.
Scholtz G. Homology and ontogeny: pattern and process in comparative developmental biology. Theory Biosci. 2005;124:121–43. doi:10.1007/BF02814480.
Wadeson PH, Crawford K. Formation of the blastoderm and yolk syncytial layer in early squid development. Biol Bull. 2003;205:179–80.
Martín-Durán JM, Egger B. Developmental diversity in free-living flatworms. Evodevo. 2012;3:7. doi:10.1186/2041-9139-3-7.
Kocot KM. On 20 years of Lophotrochozoa. Org Divers Evol. 2016;16:329–43. doi:10.1007/s13127-015-0261-3.
Dunn CW, Giribet G, Edgecombe GD, Hejnol A. Animal phylogeny and its evolutionary implications. Annu Rev Ecol Evol Syst. 2014;45:371–95. doi:10.1146/annurev-ecolsys-120213-091627.
Nesnidal MP, Helmkampf M, Meyer A, Witek A, Bruchhaus I, Ebersberger I, et al. New phylogenomic data support the monophyly of Lophophorata and an Ectoproct-Phoronid clade and indicate that Polyzoa and Kryptrochozoa are caused by systematic bias. BMC Evol Biol. 2013;13:253. doi:10.1186/1471-2148-13-253.
Struck TH, Wey-Fabrizius AR, Golombek A, Hering L, Weigert A, Bleidorn C, et al. Platyzoan paraphyly based on phylogenomic data supports a noncoelomate ancestry of Spiralia. Mol Biol Evol. 2014;31:1833–49. doi:10.1093/molbev/msu143.
Laumer CE, Bekkouche N, Kerbl A, Goetz F, Neves RC, Sørensen MV, et al. Spiralian phylogeny informs the evolution of microscopic lineages. Curr Biol. 2015;25:2000–6. doi:10.1016/j.cub.2015.06.068.
Kocot KM, Struck TH, Merkel J, Waits DS, Todt C, Brannock PM, et al. Phylogenomics of Lophotrochozoa with consideration of systematic error. Syst Biol. 2017;66(2):256–82. doi:10.1093/sysbio/syw079.
Santagata S. Ectoprocta. In: Evolutionary Developmental Biology of Invertebrates 2. Vienna: Springer; 2015. p. 247–62. doi:10.1007/978-3-7091-1871-9\_11.
Santagata S. Brachiopoda. In: Evolutionary Developmental Biology of Invertebrates 2. Vienna: Springer; 2015. p. 263–77. doi:10.1007/978-3-7091-1871-9\_12.
Hejnol A. Gastrotricha. In: Evolutionary Developmental Biology of Invertebrates 2. Vienna: Springer; 2015. p. 13–9. doi:10.1007/978-3-7091-1871-9\_2.
Hejnol A. Gnathifera. In: Evolutionary Developmental Biology of Invertebrates 2. Vienna: Springer; 2015. p. 1–12. doi:10.1007/978-3-7091-1871-9\_1.
Dunn CW, Hejnol A, Matus DQ, Pang K, Browne WE, Smith SA, et al. Broad phylogenomic sampling improves resolution of the animal tree of life. Nature. 2008;452:745–9. doi:10.1038/nature06614.
Reed CG. Bryozoa. In: Reproduction of Marine Invertebrates. Giese AC, Pearse JS, Pearse VB, editors. Pacific Grove: The Boxwood Press; 1991. pp. 85–245.
Zimmer RL. Phoronids, brachiopods, and bryozoans, the lophophorates. In: Gilbert SF, Raunio AM, editors. Embryology: Constructing the organism. Sunderland, MA: Sinauer Associates, Inc.; 1997. p. 279–305.
Barrois J. Recherches sur l’embryologie des bryozoaires. Trav Stn Zool Wimereux. 1877;1:1–305.
Prouho H. Contribution a l’histoire des Bryozoaires. Arch Zool Exp Gen. 1892;2:557–656.
Calvet L. Contributions à l’histoire naturelle des bryozoaires ectoproctes marins. Travaux de l’Institut de Zoologie de l’Université de Montpellier. Montpellier: Coulet et Fils, Editeurs; 1900.
Pace RM. On the early stages in the development of Flustrella hispida (Fabricius), and on the existence of a “yolk nucleus” in the egg of this form. Q J Microsc Sci. 1906;s2-50:435–78.
Marcus E. Bryozoarios Marinhos Brasileiros II. Bol Fac Phil Sc Letr Univ S Paulo Zoologia. 1938;2:3–137.
Corrêa DD. A embriologia de Bugula flabellata (J. V. Thompson) (Bryozoa Ectoprocta). Bol Fac Fil Ciênc e Letr Univ São Paulo - Zoologia. 1948;13:7–71.
Hyman LH. The lophophorate Coelomates—Phylum ectoprocta. In: The invertebrates: Smaller coelomate groups. New York: McGraw-Hill Book Company, Inc; 1959. p. 275–501.
Nielsen C. Phylum Ectoprocta. In: Animal Evolution: Interrelationships of the Living Phyla. 2nd ed. New York: Oxford University Press; 2001. p. 244–63.
Nielsen C. Trochophora larvae: Cell-lineages, ciliary bands and body regions. 2. Other groups and general discussion. J Exp Zool B Mol Dev Evol. 2005;304:401–47. doi:10.1002/jez.b.21050.
Gruhl A. Ultrastructure of mesoderm formation and development in Membranipora membranacea (Bryozoa: Gymnolaemata). Zoomorphology. 2009;129:45–60. doi:10.1007/s00435-009-0099-3.
Temkin MH. Gamete spawning and fertilization in the gymnolaemate bryozoan Membranipora membranacea. Biol Bull. 1994;187:143–55. doi:10.2307/1542237.
Stricker SA, Reed CG, Zimmer RL. The cyphonautes larva of the marine bryozoan Membranipora membranacea. I. General morphology, body wall, and gut. Can J Zool. 1988;66:368–83. doi:10.1139/z88-054.
Reed CG. Phylum Bryozoa. In: Strathmann MF, editor. Reproduction and Development of Marine Invertebrates of the Northern Pacific Coast: Data and Methods for the Study of Eggs, Embryos, and Larvae. Seattle: University of Washington Press; 1987. p. 494–510.
Martín-Durán JM, Vellutini BC, Hejnol A. Embryonic chirality and the evolution of spiralian left-right asymmetries. Philos Trans R Soc Lond B Biol Sci. 2016;371(1710). doi:10.1098/rstb.2015.0411.
Lambert JD, Nagy LM. MAPK signaling by the D quadrant embryonic organizer of the mollusc Ilyanassa obsoleta. Development. 2001;128:45–56.
Lambert JD, Nagy LM. The MAPK cascade in equally cleaving spiralian embryos. Dev Biol. 2003;263:231–41. doi:10.1016/j.ydbio.2003.07.006.
Henry JJ, Perry KJ. MAPK activation and the specification of the D quadrant in the gastropod mollusc, Crepidula fornicata. Dev Biol. 2007;313:181–95. doi:10.1016/j.ydbio.2007.10.019.
Koop D, Richards GS, Wanninger A, Gunter HM, Degnan BM. The role of MAPK signaling in patterning and establishing axial symmetry in the gastropod Haliotis asinina. Dev Biol. 2007;311:200–12. doi:10.1016/j.ydbio.2007.08.035.
Amiel AR, Henry JQ, Seaver EC. An organizing activity is required for head patterning and cell fate specification in the polychaete annelid Capitella teleta: New insights into cell-cell signaling in Lophotrochozoa. Dev Biol. 2013;379:107–22. doi:10.1016/j.ydbio.2013.04.011.
Pfeifer K, Schaub C, Domsch K, Dorresteijn A, Wolfstetter G. Maternal inheritance of twist and analysis of MAPK activation in embryos of the polychaete annelid Platynereis dumerilii. PLoS One. 2014;9:e96702. doi:10.1371/journal.pone.0096702.
Kozin VV, Filimonova DA, Kupriashova EE, Kostyuchenko RP. Mesoderm patterning and morphogenesis in the polychaete Alitta virens (Spiralia, Annelida): Expression of mesodermal markers Twist, Mox, Evx and functional role for MAP kinase signaling. Mech Dev. 2016;140:1–11. doi:10.1016/j.mod.2016.03.003.
Sinigaglia C, Busengdal H, Leclère L, Technau U, Rentzsch F. The bilaterian head patterning gene six3/6 controls aboral domain development in a cnidarian. PLoS Biol. 2013;11:e1001488. doi:10.1371/journal.pbio.1001488.
Lowe CJ, Wu M, Salic A, Evans L, Lander E, Stange-Thomann N, et al. Anteroposterior patterning in hemichordates and the origins of the chordate nervous system. Cell. 2003;113:853–65. doi:10.1016/S0092-8674(03)00469-0.
Steinmetz PRH, Urbach R, Posnien N, Eriksson J, Kostyuchenko RP, Brena C, et al. Six3 demarcates the anterior-most developing brain region in bilaterian animals. Evodevo. 2010;1:14. doi:10.1186/2041-9139-1-14.
Nielsen C, Worsaae K. Structure and occurrence of cyphonautes larvae (Bryozoa, Ectoprocta). J Morphol. 2010;271:1094–109. doi:10.1002/jmor.10856.
Panganiban G, Rubenstein JLR. Developmental functions of the Distal-less/Dlx homeobox genes. Development. 2002;129:4371–86.
Boncinelli E, Gulisano M, Broccoli V. Emx and Otx homeobox genes in the developing mouse brain. J Neurobiol. 1993;24:1356–66. doi:10.1002/neu.480241008.
Umesono Y, Watanabe K, Agata K. Distinct structural domains in the planarian brain defined by the expression of evolutionarily conserved homeobox genes. Dev Genes Evol. 1999;209:31–9.
Arendt D, Technau U, Wittbrodt J. Evolution of the bilaterian larval foregut. Nature. 2001;409:81–5. doi:10.1038/35051075.
Steinmetz PRH, Kostyuchenko RP, Fischer A, Arendt D. The segmental pattern of otx, gbx, and Hox genes in the annelid Platynereis dumerilii. Evol Dev. 2011;13:72–9. doi:10.1111/j.1525-142X.2010.00457.x.
Martín-Durán JM, Janssen R, Wennberg S, Budd GE, Hejnol A. Deuterostomic development in the protostome Priapulus caudatus. Curr Biol. 2012;22:2161–6. doi:10.1016/j.cub.2012.09.037.
Marlow H, Tosches MA, Tomer R, Steinmetz PR, Lauri A, Larsson T, et al. Larval body patterning and apical organs are conserved in animal evolution. BMC Biol. 2014;12:7. doi:10.1186/1741-7007-12-7.
Harada Y, Okai N, Taguchi S, Tagawa K, Humphreys T, Satoh N. Developmental expression of the hemichordate otx ortholog. Mech Dev. 2000;91:337–9. doi:10.1016/S0925-4773(99)00279-8.
Mitsunaga-Nakatsubo K, Akasaka K, Sakamoto N, Takata K, Matsumura Y, Kitajima T, et al. Differential expression of sea urchin Otx isoform (hpOtxE and HpOtxL) mRNAs during early development. Int J Dev Biol. 2003;42:645–51.
Hinman VF, Nguyen AT, Davidson EH. Expression and function of a starfish Otx ortholog, Am Otx: a conserved role for Otx proteins in endoderm development that predates divergence of the eleutherozoa. Mech Dev. 2003;120:1165–76. doi:10.1016/j.mod.2003.08.002.
Shimamura K, Hartigan DJ, Martinez S, Puelles L, Rubenstein JL. Longitudinal organization of the anterior neural plate and neural tube. Development. 1995;121:3923–33.
Venkatesh TV, Holland ND, Holland LZ, Su M-T, Bodmer R. Sequence and developmental expression of amphioxus AmphiNk2–1: Insights into the evolutionary origin of the vertebrate thyroid gland and forebrain. Dev Genes Evol. 1999;209:254–9. doi: 10.1007/s004270050250 .
Takacs CM, Amore G, Oliveri P, Poustka AJ, Wang D, Burke RD, et al. Expression of an NK2 homeodomain gene in the apical ectoderm defines a new territory in the early sea urchin embryo. Dev Biol. 2004;269:152–64. doi:10.1016/j.ydbio.2004.01.023.
Oliveri P, Walton KD, Davidson EH, McClay DR. Repression of mesodermal fate by foxa, a key endoderm regulator of the sea urchin embryo. Development. 2006;133:4173–81. doi:10.1242/dev.02577.
Arenas-Mena C. Embryonic expression of HeFoxA1 and HeFoxA2 in an indirectly developing polychaete. Dev Genes Evol. 2006;216:727–36. doi:10.1007/s00427-006-0099-y.
Boyle MJ, Seaver EC. Expression of FoxA and GATA transcription factors correlates with regionalized gut development in two lophotrochozoan marine worms: Chaetopterus (Annelida) and Themiste lageniformis (Sipuncula). Evodevo. 2010;1:2. doi:10.1186/2041-9139-1-2.
Extavour CGM, Akam M. Mechanisms of germ cell specification across the metazoans: Epigenesis and preformation. Development. 2003;130:5869–84. doi:10.1242/dev.00804.
Juliano CE, Swartz SZ, Wessel GM. A conserved germline multipotency program. Development. 2010;137:4113–26. doi:10.1242/dev.047969.
Technau U. Brachyury, the blastopore and the evolution of the mesoderm. Bioessays. 2001;23:788–94. doi:10.1002/bies.1114.
Technau U, Scholz CB. Origin and evolution of endoderm and mesoderm. Int J Dev Biol. 2003;47:531–9.
Häcker U, Kaufmann E, Hartmann C, Jürgens G, Knöchel W, Jäckle H. The Drosophila fork head domain protein crocodile is required for the establishment of head structures. EMBO J. 1995;14:5306–17.
Shimeld SM, Boyle MJ, Brunet T, Luke GN, Seaver EC. Clustered Fox genes in lophotrochozoans and the evolution of the bilaterian Fox gene cluster. Dev Biol. 2010;340:234–48. doi:10.1016/j.ydbio.2010.01.015.
Passamaneck YJ, Hejnol A, Martindale MQ. Mesodermal gene expression during the embryonic and larval development of the articulate brachiopod Terebratalia transversa. Evodevo. 2015;6:10. doi:10.1186/s13227-015-0004-8.
Zaffran S, Küchler A, Lee HH, Frasch M. Biniou (FoxF), a central component in a regulatory network controlling visceral mesoderm development and midgut morphogenesis in Drosophila. Genes Dev. 2001;15:2900–15. doi:10.1101/gad.917101.
Pérez Sánchez C, Casas-Tintó S, Sánchez L, Rey-Campos J, Granadino B. DmFoxF, a novel Drosophila fork head factor expressed in visceral mesoderm. Mech Dev. 2002;111:163–6. doi:10.1016/S0925-4773(01)00603-7.
Mazet F, Amemiya CT, Shimeld SM. An ancient Fox gene cluster in bilaterian animals. Curr Biol. 2006;16:R314–6. doi:10.1016/j.cub.2006.03.088.
Patient RK, McGhee JD. The GATA family (vertebrates and invertebrates). Curr Opin Genet Dev. 2002;12:416–22. doi:10.1016/S0959-437X(02)00319-2.
Hejnol A, Obst M, Stamatakis A, Ott M, Rouse GW, Edgecombe GD, et al. Assessing the root of bilaterian animals with scalable phylogenomic methods. Proc Biol Sci. 2009;276:4261–70. doi:10.1098/rspb.2009.0896.
Freeman G, Lundelius JW. Evolutionary implications of the mode of D quadrant specification in coelomates with spiral cleavage. J Evol Biol. 1992;5:205–47. doi:10.1046/j.1420-9101.1992.5020205.x.
van den Biggelaar JA. Development of dorsoventral polarity and mesentoblast determination in Patella vulgata. J Morphol. 1977;154:157–86. doi:10.1002/jmor.1051540111.
Martindale MQ, Doe CQ, Morrill JB. The role of animal-vegetal interaction with respect to the determination of dorsoventral polarity in the equal-cleaving spiralian, Lymnaea palustris. Wilhelm Roux Arch Entwickl Mech Org. 1985;194:281–95. doi:10.1007/BF01152174.
Martindale MQ. The “organizing” role of the D quadrant in an equal-cleaving spiralian, Lymnaea stagnalis as studied by UV laser deletion of macromeres at intervals between third and fourth quartet formation. Int J Invertebr Reprod Dev. 1986;9:229–42. doi:10.1080/01688170.1986.10510198.
Gonzales EE, van der Zee M, Dictus WJAG, van den Biggelaar J. Brefeldin A or monensin inhibits the 3D organizer in gastropod, polyplacophoran, and scaphopod molluscs. Dev Genes Evol. 2006;217:105–18. doi:10.1007/s00427-006-0118-z.
Favata MF, Horiuchi KY, Manos EJ, Daulerio AJ, Stradley DA, Feeser WS, et al. Identification of a novel inhibitor of mitogen-activated protein kinase kinase. J Biol Chem. 1998;273:18623–32. doi:10.1074/jbc.273.29.18623.
Lartillot N, Lespinet O, Vervoort M, Adoutte A. Expression pattern of Brachyury in the mollusc Patella vulgata suggests a conserved role in the establishment of the AP axis in Bilateria. Development. 2002;129:1411–21.
Guralnick RP, Lindberg DR. Reconnecting cell and animal lineages: What do cell lineages tell us about the evolution and development of Spiralia? Evolution. 2001;55:1501–19. doi:10.1111/j.0014-3820.2001.tb00671.x.
Clement AC. Experimental studies on germinal localization in Ilyanassa. I. The role of the polar lobe in determination of the cleavage pattern and its influence in later development. J Exp Zool. 1952;121:593–625. doi:10.1002/jez.1401210310.
Goulding MQ. Cell lineage of the Ilyanassa embryo: Evolutionary acceleration of regional differentiation during early development. PLoS One. 2009;4:e5506. doi: 10.1371/journal.pone.0005506 .
Mead AD. The early development of marine annelids. J Morphol. 1897;13:227–326. doi:10.1002/jmor.1050130202.
Wilson EB. Considerations on cell-lineage and ancestral reminiscence, based on a re-examination of some points in the early development of annelids and polyclades. Science. 1898;7:225–6. doi:10.2307/1624852.
Nielsen C. Trochophora larvae: cell-lineages, ciliary bands, and body regions. 1. Annelida and Mollusca. J Exp Zool B Mol Dev Evol. 2004;302:35–68. doi:10.1002/jez.b.20001.
Damen P, Dictus WJ. Cell lineage of the prototroch of Patella vulgata (Gastropoda, Mollusca). Dev Biol. 1994;162:364–83. doi:10.1006/dbio.1994.1094.
Dictus WJ, Damen P. Cell-lineage and clonal-contribution map of the trochophore larva of Patella vulgata (Mollusca). Mech Dev. 1997;62:213–26. doi:10.1016/S0925-4773(97)00666-7.
Henry JJ, Martindale MQ. Conservation of the spiralian developmental program: Cell lineage of the nemertean, Cerebratulus lacteus. Dev Biol. 1998;201:253–69. doi:10.1006/dbio.1998.8966.
Boyer BC, Henry JJ, Martindale MQ. The cell lineage of a polyclad turbellarian embryo reveals close similarity to coelomate spiralians. Dev Biol. 1998;204:111–23. doi:10.1006/dbio.1998.9084.
Maslakova SA, Martindale MQ, Norenburg JL. Fundamental properties of the spiralian developmental program are displayed by the basal nemertean Carinoma tremaphoros (Palaeonemertea, Nemertea). Dev Biol. 2004;267:342–60. doi:10.1016/j.ydbio.2003.10.022.
Henry JQ, Okusu A, Martindale MQ. The cell lineage of the polyplacophoran, Chaetopleura apiculata: variation in the spiralian program and implications for molluscan evolution. Dev Biol. 2004;272:145–60. doi:10.1016/j.ydbio.2004.04.027.
Ackermann C, Dorresteijn A, Fischer A. Clonal domains in postlarval Platynereis dumerilii (Annelida: Polychaeta). J Morphol. 2005;266:258–80. doi:10.1002/jmor.10375.
Hejnol A, Martindale MQ, Henry JQ. High-resolution fate map of the snail Crepidula fornicata: The origins of ciliary bands, nervous system, and muscular elements. Dev Biol. 2007;305:63–76. doi:10.1016/j.ydbio.2007.01.044.
Meyer NP, Boyle MJ, Martindale MQ, Seaver EC. A comprehensive fate map by intracellular injection of identified blastomeres in the marine polychaete Capitella teleta. Evodevo. 2010;1:8. doi:10.1186/2041-9139-1-8.
Rouse GW. Trochophore concepts: ciliary bands and the evolution of larvae in spiralian Metazoa. Biol J Linn Soc Lond. 1999;66(4):411–64. doi:10.1111/j.1095-8312.1999.tb01920.x.
Henry JQ, Hejnol A, Perry KJ, Martindale MQ. Homology of ciliary bands in Spiralian Trochophores. Integr Comp Biol. 2007;47:865–71. doi:10.1093/icb/icm035.
Pineda D, Saló E. Planarian gtsix3, a member of the six/so gene family, is expressed in brain branches but not in eye cells. Mech Dev. 2002;119 Suppl 1:S167–71. doi:10.1016/S0925-4773(03)00111-4.
Santagata S, Resh C, Hejnol A, Martindale MQ, Passamaneck YJ. Development of the larval anterior neurogenic domains of Terebratalia transversa (Brachiopoda) provides insights into the diversification of larval apical organs and the spiralian nervous system. Evodevo. 2012;3:3. doi:10.1186/2041-9139-3-3.
Hiebert LS, Maslakova SA. Hox genes pattern the anterior-posterior axis of the juvenile but not the larva in a maximally indirect developing invertebrate, Micrura alaskensis (Nemertea). BMC Biol. 2015;13:23. doi:10.1186/s12915-015-0133-5.
Hiebert LS, Maslakova SA. Expression of Hox, Cdx, and Six3/6 genes in the hoplonemertean Pantinonemertes californiensis offers insight into the evolution of maximally indirect development in the phylum Nemertea. Evodevo. 2015;6:26. doi:10.1186/s13227-015-0021-7.
Martín-Durán JM, Vellutini BC, Hejnol A. Evolution and development of the adelphophagic, intracapsular Schmidt’s larva of the nemertean Lineus ruber. Evodevo. 2015;6:28. doi:10.1186/s13227-015-0023-5.
Perry KJ, Lyons DC, Truchado-Garcia M, Fischer AHL, Helfrich LW, Johansson KB, et al. Deployment of regulatory genes during gastrulation and germ layer specification in a model spiralian mollusc Crepidula. Dev Dyn. 2015;244(10):1215–48. doi:10.1002/dvdy.24308.
Nederbragt AJ, te Welscher P, van den Driesche S, van Loon AE, Dictus WJAG. Novel and conserved roles for orthodenticle/otx and orthopedia/otp orthologs in the gastropod mollusc Patella vulgata. Dev Genes Evol. 2002;212:330–7. doi:10.1007/s00427-002-0246-z.
Render J. Cell fate maps in the Ilyanassa obsoleta embryo beyond the third division. Dev Biol. 1997;189:301–10. doi:10.1006/dbio.1997.8654.
Lyons D, Perry K, Henry J. Spiralian gastrulation: germ layer formation, morphogenesis, and fate of the blastopore in the slipper snail Crepidula fornicata. Evodevo. 2015;6:24. doi:10.1186/s13227-015-0019-1.
Lankester ER. Notes on the embryology and classification of the animal kingdom: comprising a revision of speculations relative to the origin and significance of the germ-layers. Q J Microsc Sci. 1877;17:399–454.
Martindale MQ, Hejnol A. A developmental perspective: changes in the position of the blastopore during bilaterian evolution. Dev Cell. 2009;17:162–74. doi:10.1016/j.devcel.2009.07.024.
Hejnol A, Martín-Durán JM. Getting to the bottom of anal evolution. Zool Anz. 2015;256:61–74. doi:10.1016/j.jcz.2015.02.006.
Martín-Durán JM, Passamaneck YJ, Martindale MQ, Hejnol A. The developmental basis for the recurrent evolution of deuterostomy and protostomy. Nat Ecol Evol. 2016;1:0005. doi:10.1038/s41559-016-0005.
Tessmar-Raible K, Raible F, Christodoulou F, Guy K, Rembold M, Hausen H, et al. Conserved sensory-neurosecretory cell types in annelid and fish forebrain: insights into hypothalamus evolution. Cell. 2007;129:1389–400. doi:10.1016/j.cell.2007.04.041.
Holland LZ, Holland NN, Schubert M. Developmental expression of AmphiWnt1, an amphioxus gene in the Wnt1/wingless subfamily. Dev Genes Evol. 2000;210:522–4. doi:10.1007/s004270050342.
Prud’homme B, de Rosa R, Arendt D, Julien J-F, Pajaziti R, Dorresteijn AWC, et al. Arthropod-like expression patterns of engrailed and wingless in the annelid Platynereis dumerilii suggest a role in segment formation. Curr Biol. 2003;13:1876–81. doi:10.1016/j.cub.2003.10.006.
Bolognesi R, Beermann A, Farzana L, Wittkopp N, Lutz R, Balavoine G, et al. Tribolium Wnts: evidence for a larger repertoire in insects with overlapping expression patterns that suggest multiple redundant functions in embryogenesis. Dev Genes Evol. 2008;218:193–202. doi:10.1007/s00427-007-0170-3.
Pani AM, Mullarkey EE, Aronowicz J, Assimacopoulos S, Grove EA, Lowe CJ. Ancient deuterostome origins of vertebrate brain signalling centres. Nature. 2012;483:289–94. doi:10.1038/nature10838.
d’Hondt J-L. Sur l’évolution des quatre macromères du pôle végétatif chez les embryons de Bryozoaires Eurystomes. Cah Biol Mar. 1983;24:177–85.
Boyer BC, Henry JQ, Martindale MQ. Dual origins of mesoderm in a basal spiralian: Cell lineage analyses in the polyclad turbellarian Hoploplana inquilina. Dev Biol. 1996;179:329–38. doi:10.1006/dbio.1996.0264.
Lartillot N, Le Gouar M, Adoutte A. Expression patterns of fork head and goosecoid homologues in the mollusc Patella vulgata supports the ancestry of the anterior mesendoderm across Bilateria. Dev Genes Evol. 2002;212:551–61. doi:10.1007/s00427-002-0274-8.
Lambert JD. Mesoderm in spiralians: the organizer and the 4d cell. J Exp Zool B Mol Dev Evol. 2008;310:15–23. doi:10.1002/jez.b.21176.
Gline SE, Nakamoto A, Cho S-J, Chi C, Weisblat DA. Lineage analysis of micromere 4d, a super-phylotypic cell for Lophotrochozoa, in the leech Helobdella and the sludgeworm Tubifex. Dev Biol. 2011;353:120–33. doi:10.1016/j.ydbio.2011.01.031.
Lyons DC, Perry KJ, Lesoway MP, Henry JQ. Cleavage pattern and fate map of the mesentoblast, 4d, in the gastropod Crepidula: a hallmark of spiralian development. Evodevo. 2012;3:21. doi:10.1186/2041-9139-3-21.
Lyons DC, Henry JQ. Ins and outs of Spiralian gastrulation. Int J Dev Biol. 2014;58:413–28. doi:10.1387/ijdb.140151dl.
Dill KK, Seaver EC. Vasa and nanos are coexpressed in somatic and germ line tissue from early embryonic cleavage stages through adulthood in the polychaete Capitella sp. I. Dev Genes Evol. 2008;218:453–63. doi:10.1007/s00427-008-0236-x.
Le Gouar M, Lartillot N, Adoutte A, Vervoort M. The expression of a caudal homologue in a mollusc, Patella vulgata. Gene Expr Patterns. 2003;3:35–7. doi:10.1016/S1567-133X(02)00091-1.
de Rosa R, Prud’homme B, Balavoine G. Caudal and even-skipped in the annelid Platynereis dumerilii and the ancestry of posterior growth. Evol Dev. 2005;7:574–87. doi:10.1111/j.1525-142X.2005.05061.x.
Fröbius AC, Seaver EC. ParaHox gene expression in the polychaete annelid Capitella sp. I. Dev Genes Evol. 2006;216:81–8. doi:10.1007/s00427-005-0049-0.
Kulakova MA, Cook CE, Andreeva TF. ParaHox gene expression in larval and postlarval development of the polychaete Nereis virens (Annelida, Lophotrochozoa). BMC Dev Biol. 2008;8:61. doi:10.1186/1471-213X-8-61.
Hui JHL, Raible F, Korchagina N, Dray N, Samain S, Magdelenat G, et al. Features of the ancestral bilaterian inferred from Platynereis dumerilii ParaHox genes. BMC Biol. 2009;7:43–3. doi:10.1186/1741-7007-7-43.
Samadi L, Steiner G. Conservation of ParaHox genes’ function in patterning of the digestive tract of the marine gastropod Gibbula varia. BMC Dev Biol. 2010;10:74. doi:10.1186/1471-213X-10-74.
Fritsch M, Wollesen T, Wanninger A. Hox and ParaHox gene expression in early body plan patterning of polyplacophoran mollusks. J Exp Zool B Mol Dev Evol. 2016;326:89–104. doi:10.1002/jez.b.22671.
Seaver EC, Yamaguchi E, Richards GS, Meyer NP. Expression of the pair-rule gene homologs runt, Pax3/7, even-skipped-1 and even-skipped-2 during larval and juvenile development of the polychaete annelid Capitella teleta does not support a role in segmentation. Evodevo. 2012;3:8. doi:10.1186/2041-9139-3-8.
Rebscher N. Establishing the germline in spiralian embyos. Int J Dev Biol. 2014;58:403–11. doi:10.1387/ijdb.140125nr.
Stricker SA. Metamorphosis of the marine bryozoan Membranipora membranacea: An ultrastructural study of rapid morphogenetic movements. J Morphol. 1988;196:53–72. doi:10.1002/jmor.1051960106.
Henry JJ, Collin R, Perry KJ. The slipper snail, Crepidula: an emerging lophotrochozoan model system. Biol Bull. 2010;218:211–29. doi:10.1086/BBLv218n3p211.
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. doi:10.1111/j.1525-142X.2007.00216.x.
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. doi:10.1111/j.1525-142X.2006.00136.x.
Wong KS-Y, Arenas-Mena C. Expression of GATA and POU transcription factors during the development of the planktotrophic trochophore of the polychaete serpulid Hydroides elegans. Evol Dev. 2016;18:254–66. doi:10.1111/ede.12196.
Martín-Durán JM, Hejnol A. The study of Priapulus caudatus reveals conserved molecular patterning underlying different gut morphogenesis in the Ecdysozoa. BMC Biol. 2015;13:29. doi:10.1186/s12915-015-0139-z.
Fuchs J, Martindale MQ, Hejnol A. Gene expression in bryozoan larvae suggest a fundamental importance of pre-patterned blastemic cells in the bryozoan life-cycle. Evodevo. 2011;2:13–3. doi:10.1186/2041-9139-2-13.
Willems M, Egger B, Wolff C, Mouton S, Houthoofd W, Fonderie P, et al. Embryonic origins of hull cells in the flatworm Macrostomum lignano through cell lineage analysis: Developmental and phylogenetic implications. Dev Genes Evol. 2009;219:409–17. doi:10.1007/s00427-009-0304-x.
Kuroda R. A twisting story: how a single gene twists a snail? Mechanogenetics. Q Rev Biophys. 2015;48:445–52. doi:10.1017/S0033583515000098.
Davison A, McDowell GS, Holden JM, Johnson HF, Koutsovoulos GD, Liu MM, et al. Formin is associated with left-right asymmetry in the pond snail and the frog. Curr Biol. 2016;26:654–60. doi:10.1016/j.cub.2015.12.071.
Brun-Usan M, Marín-Riera M, Grande C, Truchado-Garcia M, Salazar-Ciudad I. A set of simple cell processes are sufficient to model spiral cleavage. Development. 2017;144(1):54–62. doi:10.1242/dev.140285.
Valentine JW. Cleavage patterns and the topology of the metazoan tree of life. Proc Natl Acad Sci U S A. 1997;94:8001–5.
Long JA, Stricker SA. Brachiopoda. In: Reproduction of Marine Invertebrates: Echinoderms and Lophophorates. Giese AC, Pearse JS, Pearse VB, editors. Pacific Grove: The Boxwood Press; 1991. pp. 47–84.
Henry JJ, Klueg KM, Raff RA. Evolutionary dissociation between cleavage, cell lineage and embryonic axes in sea urchin embryos. Development. 1992;114:931–8.
Kuo D-H. The polychaete-to-clitellate transition: An EvoDevo perspective. Dev Biol. 2017. doi:10.1016/j.ydbio.2017.01.016.
Hejnol A, Schnabel R. What a couple of dimensions can do for you: Comparative developmental studies using 4D microscopy—examples from tardigrade development. Integr Comp Biol. 2006;46:151–61. doi:10.1093/icb/icj012.
Schnabel R, Hutter H, Moerman D, Schnabel H. Assessing normal embryogenesis in Caenorhabditis elegans using a 4D microscope: variability of development and regional specification. Dev Biol. 1997;184:234–65. doi:10.1006/dbio.1997.8509.
Wolff C, Tinevez J-Y, Pietzsch T, Stamataki E, Harich B, Preibisch S, et al. Reconstruction of cell lineages and behaviors underlying arthropod limb outgrowth with multi-view light-sheet imaging and tracking. bioRxiv. 2017;112623. doi:10.1101/112623.
Pietzsch T, Saalfeld S, Preibisch S, Tomancak P. BigDataViewer: visualization and processing for large image data sets. Nat Methods. 2015;12:481–3. doi:10.1038/nmeth.3392.
Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: An open-source platform for biological-image analysis. Nat Methods. 2012;9:676–82. doi:10.1038/nmeth.2019.
Vellutini BC, Martín-Durán JM, Hejnol A. Supplementary files for the bryozoan cell lineage project. Figshare. 2017. doi:10.6084/m9.figshare.4563574.v2.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52. doi:10.1038/nbt.1883.
Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, et al. Primer3–new capabilities and interfaces. Nucleic Acids Res. 2012;40:e115–5. doi:10.1093/nar/gks596.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol Biol Evol. 2013;30:772–80. doi:10.1093/molbev/mst010.
Talavera G, Castresana J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol. 2007;56:564–77. doi:10.1080/10635150701472164.
Stamatakis A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3. doi:10.1093/bioinformatics/btu033.
Okonechnikov K, Golosova O, Fursov M, UGENE team. Unipro UGENE: A unified bioinformatics toolkit. Bioinformatics. 2012;28:1166–7. doi:10.1093/bioinformatics/bts091.
Huerta-Cepas J, Dopazo J, Gabaldón T. ETE: A python environment for tree exploration. BMC Bioinformatics. 2010;11:24. doi:10.1186/1471-2105-11-24.
Staudt T, Lang MC, Medda R, Engelhardt J, Hell SW. 2,2′-Thiodiethanol: a new water soluble mounting medium for high resolution optical microscopy. Microsc Res Tech. 2007;70:1–9. doi:10.1002/jemt.20396.
Asadulina A, Panzera A, Verasztó C, Liebig C, Jékely G. Whole-body gene expression pattern registration in Platynereis larvae. Evodevo. 2012;3:27. doi:10.1186/2041-9139-3-27.
Jékely G, Arendt D. Cellular resolution expression profiling using confocal detection of NBT/BCIP precipitate by reflection microscopy. Biotechniques. 2007;42:751–5. doi:10.2144/000112462.
Riedl RJ. Gnathostomulida from america. Science. 1969;163:445–52. doi:10.1126/science.163.3866.445.
Pennerstorfer M, Scholtz G. Early cleavage in Phoronis muelleri (Phoronida) displays spiral features. Evol Dev. 2012;14:484–500. doi:10.1111/ede.12002.
Marcus E. Bryozoarios Marinhos Brasileiros III. Bol Fac Phil, Sc Letr Univ S Paulo, Zoologia. 1939;3:113–299.
Merkel J, Wollesen T, Lieb B, Wanninger A. Spiral cleavage and early embryology of a loxosomatid entoproct and the usefulness of spiralian apical cross patterns for phylogenetic inferences. BMC Dev Biol. 2012;12:11. doi:10.1186/1471-213X-12-11.
Boell LA, Bucher G. Whole-mount in situ hybridization in the rotifer Brachionus plicatilis representing a basal branch of lophotrochozoans. Dev Genes Evol. 2008;218:445–51. doi:10.1007/s00427-008-0234-z.
Smith JM, Cridge AG, Dearden PK. Germ cell specification and ovary structure in the rotifer Brachionus plicatilis. Evodevo. 2010;1:5. doi:10.1186/2041-9139-1-5.
Lapan SW, Reddien PW. dlx and sp6-9 Control optic cup regeneration in a prototypic eye. PLoS Genet. 2011;7:e1002226. doi:10.1371/journal.pgen.1002226.
Callaerts P, Munoz-Marmol AM, Glardon S, Castillo E, Sun H, Li W-H, et al. Isolation and expression of a pax-6 gene in the regenerating and intact planarian Dugesia(G)tigrina. Proc Natl Acad Sci. 1999;96:558–63. doi:10.1073/pnas.96.2.558.
Garcia-Fernàndez J, Baguñà J, Saló E. Genomic organization and expression of the planarian homeobox genes Dth-1 and Dth-2. Development. 1993;118:241–53.
Martín-Durán JM, Amaya E, Romero R. Germ layer specification and axial patterning in the embryonic development of the freshwater planarian Schmidtea polychroa. Dev Biol. 2010;340:145–58. doi:10.1016/j.ydbio.2010.01.018.
Martín-Durán JM, Romero R. Evolutionary implications of morphogenesis and molecular patterning of the blind gut in the planarian Schmidtea polychroa. Dev Biol. 2011;352:164–76. doi:10.1016/j.ydbio.2011.01.032.
Sato K, Shibata N, Orii H, Amikura R, Sakurai T, Agata K, et al. Identification and origin of the germline stem cells as revealed by the expression of nanos-related gene in planarians. Dev Growth Differ. 2006;48:615–28. doi:10.1111/j.1440-169X.2006.00897.x.
Handberg-Thorsager M, Saló E. The planarian nanos-like gene smednos is expressed in germline and eye precursor cells during development and regeneration. Dev Genes Evol. 2007;217:403–11. doi:10.1007/s00427-007-0146-3.
Lee SE, Jacobs DK. Expression of Distal-less in molluscan eggs, embryos, and larvae. Evol Dev. 1999;1:172–9.
Tomarev SI, Callaerts P, Kos L, Zinovieva R, Halder G, Gehring W, et al. Squid pax-6 and eye development. Proc Natl Acad Sci U S A. 1997;94:2421–6.
Kranz AM, Tollenaere A, Norris BJ, Degnan BM, Degnan SM. Identifying the germline in an equally cleaving mollusc: Vasa and Nanos expression during embryonic and larval development of the vetigastropod Haliotis asinina. J Exp Zool B Mol Dev Evol. 2010;314:267–79. doi:10.1002/jez.b.21336.
Rabinowitz JS, Chan XY, Kingsley EP, Duan Y, Lambert JD. Nanos is required in somatic blast cell lineages in the posterior of a mollusk embryo. Curr Biol. 2008;18:331–6. doi:10.1016/j.cub.2008.01.055.
McDougall C, Korchagina N, Tobin JL, Ferrier DE. Annelid Distal-less/Dlx duplications reveal varied post-duplication fates. BMC Evol Biol. 2011;11:241. doi:10.1186/1471-2148-11-241.
Boyle MJ, Yamaguchi E, Seaver EC. Molecular conservation of metazoan gut formation: Evidence from expression of endomesoderm genes in Capitella teleta (Annelida). Evodevo. 2014;5:39. doi:10.1186/2041-9139-5-39.
Arendt D, Tessmar K, de Campos-Baptista M-IM, Dorresteijn A, Wittbrodt J. Development of pigment-cup eyes in the polychaete Platynereis dumerilii and evolutionary conservation of larval eyes in Bilateria. Development. 2002;129:1143–54.
Denes AS, Jékely G, Steinmetz PRH, Raible F, Snyman H, Prud’homme B, et al. Molecular architecture of annelid nerve cord supports common origin of nervous system centralization in Bilateria. Cell. 2007;129:277–88. doi:10.1016/j.cell.2007.02.040.
Quigley IK, Xie X, Shankland M. Hau-Pax6A expression in the central nervous system of the leech embryo. Dev Genes Evol. 2007;217:459–68. doi:10.1007/s00427-007-0156-1.
Rebscher N, Zelada-González F, Banisch TU, Raible F, Arendt D. Vasa unveils a common origin of germ cells and of somatic stem cells from the posterior growth zone in the polychaete Platynereis dumerilii. Dev Biol. 2007;306:599–611. doi:10.1016/j.ydbio.2007.03.521.
Kang D, Pilon M, Weisblat DA. Maternal and zygotic expression of a nanos-class gene in the leech Helobdella robusta: primordial germ cells arise from segmental mesoderm. Dev Biol. 2002;245:28–41. doi:10.1006/dbio.2002.0615.
Passamaneck YJ, Furchheim N, Hejnol A, Martindale MQ, Lüter C. Ciliary photoreceptors in the cerebral eyes of a protostome larva. Evodevo. 2011;2:6. doi:10.1186/2041-9139-2-6.
Vellutini BC, Hejnol A. Expression of segment polarity genes in brachiopods supports a non-segmental ancestral role of engrailed for bilaterians. Sci Rep. 2016;6:32387. doi:10.1038/srep32387.
Petersen CP, Reddien PW. Smed-catenin-1 is required for anteroposterior blastema polarity in planarian regeneration. Science. 2008;319:327–30. doi:10.1126/science.1149943.
Adell T, Salò E, Boutros M, Bartscherer K. Smed-Evi/Wntless is required for beta-catenin-dependent and -independent processes during planarian regeneration. Development. 2009;136:905–10. doi:10.1242/dev.033761.
Nederbragt AJ, Lespinet O, van Wageningen S, van Loon AE, Adoutte A, Dictus WJAG. A lophotrochozoan twist gene is expressed in the ectomesoderm of the gastropod mollusk Patella vulgata. Evol Dev. 2002;4:334–43. doi:10.1046/j.1525-142X.2002.02020.x.
Song MH, Huang FZ, Chang GY, Weisblat DA. Expression and function of an even-skipped homolog in the leech Helobdella robusta. Development. 2002;129:3681–92.
Seaver EC, Kaneshige LM. Expression of ‘segmentation’ genes during larval and juvenile development in the polychaetes Capitella sp. I and H. elegans. Dev Biol. 2006;289:179–94. doi:10.1016/j.ydbio.2005.10.025.
Dill KK, Thamm K, Seaver EC. Characterization of twist and snail gene expression during mesoderm and nervous system development in the polychaete annelid Capitella sp. I. Dev Genes Evol. 2007;217:435–47. doi:10.1007/s00427-007-0153-4.
Pfeifer K, Schaub C, Wolfstetter G, Dorresteijn A. Identification and characterization of a twist ortholog in the polychaete annelid Platynereis dumerilii reveals mesodermal expression of Pdu-twist. Dev Genes Evol. 2013;223:319–28. doi:10.1007/s00427-013-0448-6.
Altenburger A, Martinez P, Wanninger A. Homeobox gene expression in Brachiopoda: the role of Not and Cdx in bodyplan patterning, neurogenesis, and germ layer specification. Gene Expr Patterns. 2011;11:427–36. doi:10.1016/j.gep.2011.07.001.
We thank the S9 members for fruitful feedback and discussions, in particular Anette Elde, Jonas Bengtsen, and Anlaug Boddington for helping with the bryozoan collections, and Sabrina Schiemann and Ralf Schnabel for the assistance with the 4D microscopy and Simi BioCell particulars. We also thank three anonymous reviewers for the thoughtful comments and constructive feedback.
The study was funded by the core budget of the Sars Centre and by The European Research Council Community’s Framework Program Horizon 2020 (2014–2020) ERC grant agreement 648861 to AH. The funding bodies had no role in the design, collection, analyses, or interpretation of data.
Availability of data and materials
AH and BCV designed the study. BCV collected the samples, performed the 4D recordings of bryozoans, and traced the cell lineages. JMMD and BCV did gene cloning and in situ hybridization. BCV analyzed the data and wrote a draft manuscript. AH, JMMD and BCV revised and contributed to the text. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Overview of the cell tracking data of four wild type embryos of M. membranacea. (A) Raw cell lineages tracked in Simi BioCell . wt1 is a recording of the animal pole providing most of the data for the aboral epithelium, wt2 is a vegetal pole view providing detailed information for the vegetal ectoderm, and wt3 and wt4 are additional recordings of the animal pole. Development time measured in hours post activation (hpa) is shown in the Y axis. (B) Number of tracked cells per embryo showing the proportion of cells by quadrant. wt1 is the most-complete cell lineage and the embryo with best coverage of the B quadrant. (C) Relative density of tracked cells per time for each embryo. The plot is complementary to the raw lineages and show that embryos wt1 and wt2 were tracked for a longer period than wt3 and wt4. The peaks indicate the moment that the maximum number of tracked cells was reached and when cell births begin to decrease, which is an indicator for the increase in the number of untracked cells. The data and code for generating the plots (B) and (C) are available at . (PNG 2314 kb)
Cell lineage variability in M. membranacea. (A) Overlap between the cell lineages of embryos wt1 to wt4 (left) and the birth time of individual cells up to 11 hpa (right). Only cells tracked in the four embryos are shown. The black horizontal lines indicate the mean birth time of a cell between the embryos, and vertical black lines the standard error. The embryo wt2 lags behind the other three embryos, but the variability in the timing of cell divisions is low. (B) Standard deviation for the birth time of a cell between different embryos by time of development. The timing of cell divisions between homologous cells does not surpass 20 min until 9 hpa. After that, the variability increases. (C) Same as (B) but only for embryos tracked from the animal pole (wt1, wt3, and wt4). Timing variability also increases over time, but the range of variation is contained within 20 min even after 9 hpa, values that are in the same order of magnitude of the variability observed in C. elegans . The data and code for generating the plots are available at . (PNG 420 kb)
Quartet synchrony in animal pole embryos of M. membranacea. Color gradient represents the time a cell took to divide after the first cell of its quartet had divided. Each column represents one embryo that has been tracked from the animal pole (wt1, wt3, and wt4). Quartets 1qi 12 and 1qe 11 show a consistent asynchronous event between the three embryos. Raw data, code to generate the plot, and additional comparative plots including all embryos are available at . (PNG 52 kb)
Detailed fate map of M. membranacea. The data reflect the consensus between embryos wt1, wt2, wt3, and wt4. Quadrant A, B, C, and D give rise to the left, anterior, right, and posterior regions of the embryo. Question marks indicate cells whose fate could not be determined. (PNG 972 kb)
Sequence of cell divisions in the vegetal ectoderm of wt2. Images from Additional file 9: Video S5, but oriented with the D quadrant to the bottom and mirrored due to the reverse chirality of this embryo (see ). (A) 28-cell stage (6.7 hpa) showing the quartets 2q (purple), 3q (blue), and 3Q (black). (B) The second quartet divides forming the founders of the vegetal 12-tet (blue). (C) The 3Q quartet begins to be internalized. White dashed line demarcates the blastoporal lip. (D) Second generation of outer (blue) and central (yellow) 12-tets. (E) Eleven of the 12 cells have divided (yellow) and surround the 4Q cells (black). The cell 3c divides later than the others. This delay was only observed in embryo wt2 and might indicate a developmental variability between embryos (see synchrony plots in ). (F) After 3c divides, the blastopore (white dashed line) is demarcated by 12 cells from the same generation. (G) The cells lining the blastopore begin to divide forming central (red) cells. Some of the outer (blue) cells also begin to divide (green). (H) The blastopore is now demarcated by seven cells third generation (red) and 3c2 from the second generation (yellow). The blastopore is now narrower. (I) The cells at the vertices of the second generation (yellow) divide forming cells at the vertices of the embryo that are not part of the blastoporal lip. (PNG 2522 kb)
Development of M. membranacea under different concentrations of the MEK inhibitor U0126. (A) Maximum intensity projection of a confocal stack for the most representative phenotype of each U0126 treatment. Ratio in the lower left corner shows the number of embryos scored for the shown phenotype versus the total number of embryos in the treatment. Phenotypes in an additional seawater-only treatment (no DMSO) were indistinguishable from DMSO control and had a ratio of 31/37 (not shown). (B) Measurements for the number of nuclei (y axis), embryo width (point size), and embryo height (color scale) for the confocal scans of (A). Each colored point represents one embryo, the black rhombus and error bars shows the mean number of nuclei with standard error. Number of embryos scanned and measured per treatment: DMSO = 3, 1 μM = 3, 10 μM = 5, 25 μM = 4. In all treatments, U0126 was added at 3 hpa (2-cell stage), the embryos developed at 10 °C and were fixed at 44 hpa. Scale bars = 20 μm. Raw data and code to generate the plot are available at . (PNG 829 kb)
M. membranacea embryos treated with the MEK inhibitor U0126 (10 μM) from different developmental stages. All treatments developed at 10 °C and were fixed at 72 hpa. Arrowhead indicates when the U0126 treatments began for each experimental condition, represented by the horizontal colored lines. Representative phenotypes are shown for the 4, 8, and 18 hpa treatments. We scored 100 embryos under light microscopy for each treatment to obtain the ratio of severe/mild phenotypes. Treatments showing the severe phenotype are shown in orange, with the percentage of severe phenotypes indicated at the right end. Treatments without severe phenotypes were colored in yellow. (PNG 865 kb)
4D recordings of M. membranacea embryos treated with the MEK inhibitor U0126. Each row corresponds to the timeline of an individual embryo. All recordings were synchronized by the timing of the second cleavage (4-cell stage = 0 h). Time scale shows the number of hours after 4-cell stage. The exact developmental time is shown on panels that do not correspond to the time shown in the main scale (top right corner). Frame number of each panel is shown on the bottom left corner (a frame was captured every 40 s). The orange rectangle indicates the cleavage abnormality observed in embryos that exhibit the severe phenotype. (PNG 6748 kb)
Gene expression throughout M. membranacea cell lineage. (A) Various developmental stages illustrating the gene expression patterns in the animal ectoderm, vegetal ectoderm, and endomesoderm with cellular resolution. (B) Cell lineage diagrams indicating the lineages where the above genes are expressed. Vivid colors indicate gene expression while more transparent branches indicate absence of expression for each particular gene analyzed. (PNG 1047 kb)
Cell lineage comparison between the larval ciliated bands of M. membranacea and Patella vulgata . (PNG 139 kb)
Orthology assignment for the bryozoan genes used in this study. (A) six3/6, (B) dlx and evx, (C) otx and gsc, (D) pax6, (E) nk2.1, (F) foxa, foxc, and foxf, (G) nanos, (H) bra, (I) cdx, (J) twist, (K) gata456, and (L) wnt1. Cladograms show branch support values and bryozoan orthologs in red. (PDF 3134 kb)
About this article
Cite this article
Vellutini, B.C., Martín-Durán, J.M. & Hejnol, A. Cleavage modification did not alter blastomere fates during bryozoan evolution. BMC Biol 15, 33 (2017). https://doi.org/10.1186/s12915-017-0371-9