- Open Access
Gene expression in bryozoan larvae suggest a fundamental importance of pre-patterned blastemic cells in the bryozoan life-cycle
EvoDevo volume 2, Article number: 13 (2011)
Bryozoa is a clade of aquatic protostomes. The bryozoan life cycle typically comprises a larval stage, which metamorphoses into a sessile adult that proliferates by asexual budding to form colonies. The homology of bryozoan larvae with other protostome larvae is enigmatic. Bryozoan larvae exhibit blastemic tissues that contribute to build the adult during morphogenesis. However, it remains unclear if the cells of these tissues are pre-determined according to their future fate or if the cells are undifferentiated, pluripotent stem cells. Gene expression studies can help to identify molecular patterning of larval and adult tissues and enlighten the evolution of bryozoan life cycle stages.
We investigated the spatial expression of 13 developmental genes in the larval stage of the gymnolaemate bryozoan Bugula neritina. We found most genes expressed in discrete regions in larval blastemic tissues that form definitive components of the adult body plan. Only two of the 13 genes, BnTropomyosin and BnFoxAB, were exclusively expressed in larval tissues that are discarded during metamorphosis.
Our results suggest that the larval blastemas in Bugula are pre-patterned according to their future fate in the adult. The gene expression patterns indicate that some of the bryozoan blastemas can be interpreted to correspond to homologous adult tissues of other animals. This study challenges an earlier proposed view that metazoan larvae share homologous undifferentiated "set-aside cells", and instead points to an independent origin of the bryozoan larval stage with respect to other lophotrochozoans.
Bryozoa (Ectoprocta) is a monophyletic group of sessile, colonial invertebrates and includes over 6,000 species in aquatic habitats worldwide . Bryozoan life history, reproduction and anatomy are so fundamentally different from other metazoan groups (for example they lack typical circulatory structures or nephridia and the nervous systems of larvae and adults are unique), that traditional morphological investigations and the fossil record failed to clarify their evolutionary history. Bryozoa belong to Lophotrochozoa , but their phylogenetic position within the group is still ambiguous [3–6]. Recently, a few molecular studies indicated a close relationship of bryozoans with the clade Entoprocta + Cycliophora, but with low support [7, 8]. Within Bryozoa, three major clades are recognized, Gymnolaemata (Eurystomata), Stenolaemata (Cyclostomata), and Phylactolaemata, but the phylogenetic interrelationships of these groups remain controversial [9–11].
Bryozoans have indirect development and their life cycle includes a sexually produced larval stage as well as asexual reproduction by cloning to give rise to colonial adult forms . The most species rich bryozoan clade with over 5,000 species, the Gymnolaemata, has evolved a fascinating diversity of reproductive mechanisms and larval forms. Less than 20 species release their eggs directly into the surrounding water, where they develop into the well known cyphonautes larvae, planktotrophic (feeding) larvae with characteristic shells. Even fewer species produce shelled larvae with a non-functioning gut called pseudocyphonautes. The great majority of gymnolaemates, however, have evolved brood protection and their embryos develop into lecitotrophic (non-feeding) "coronate" larvae. Coronate larvae usually lack both a shell and a gut, but traces of a non-functioning gut were observed in some species, which was interpreted as a vestigial gut [13, 14]. Despite the above outlined differences among gymnolaemate larvae, they usually share a set of morphological characters including transitory larval structures such as larval muscles and nerves (such as an apical organ), a glandulo-sensory organ (pyriform organ), internal cells, as well as blastemic cells that give rise to definitive adult tissues . During a drastic metamorphosis, the transitory larval structures are discarded and only a single species is known to retain its larval gut in the adult . Typically, all adult structures (for example gut, nervous system) are built de novo from the blastemic cells, which are found in various positions in the larvae and can give rise to different adult structures during metamorphosis in individual gymnolaemate species [16, 17]. However, some authors propose that only ectodermal and mesodermal cells are involved in adult body plan formation from the larva, and this process mirrors the asexual budding process, which involves proliferation of the polypide (gut and lophophore) from the pluripotent body wall [12, 14, 18, 19].
A comprehensive cell fate map does not exist for any species of bryozoan, so it is unknown if the blastemic cells in the larvae are undifferentiated, pluri- or multipotent set-aside/stem cells, or if these cells are determined for their future fate in the adult. To help to clarify the question to what extent the adult tissues are already determined in the larval stage, we investigated the spatial expression of 13 developmental genes in the newly released, coronate larvae of the gymnolaemate bryozoan Bugula neritina. This species is one of the better-studied bryozoans with respect to larval morphology and metamorphosis, although many details are still lacking. This initial set of genes was chosen according to the specific expression of genes in certain metazoan germ layers or organs as well as their putative conserved functions in the development among taxa, and represents a foundation for future molecular studies. Amongst the investigated genes, Tropomyosin is known to be a general marker of bilaterian musculature; the genes Hox4, SoxB2, SoxE and FoxB have functions in neural development; the genes FoxA, GATA456 and Cdx (Caudal) have been shown to be involved in gut patterning among taxa; expression of GATA123 and BAMBI was observed in the ectoderm of metazoans, and Wnt genes have been shown to be involved in multiple events including axial patterning (see Discussion). Genes engaged in bilaterian gut formation were investigated to determine if there is any molecular indication for a larval gut being present as anlage of the adult gut in B. neritina. We discuss the gene expression patterns of the bryozoan larval stage in the light of available data for other metazoans.
The larval anatomy of Bugula neritina has been previously studied and relevant data are reviewed here together with our own observations (Figure 1). The larvae have a roundish body shape (length, 260 μm; width, 190 μm) and the apical disc, which leads in swimming direction, marks the apical pole (Figure 1A, B). The apical disc comprises a central neural area (apical organ), which is surrounded by an outer epidermal blastema and an inner mesodermal blastema (Figure 1E) [14, 20–22]. The apical organ is detected in semi-thin sections (not shown). The apical disc is encircled by infolded ectodermal cells, the pallial epithelium (Figure 1A, E). Opposite of the apical disc lies the large internal sac, which constitutes nearly half of the larval interior (Figure 1B). The internal sac is an ectodermal invagination used for attachment to the substrate during metamorphosis (see below) and comprises the neck, wall, and roof region . The wall region is extensive and is folded upon itself (Figure 1D, E). The roof cells are elongated and contain basal inclusions of electron transparent vesicles . The neck cells are recognized by their dark inclusions (Figure 1D, E). Here, we define the region of the internal sac opening as "abapical" instead of the formerly used "oral", since coronate larvae lack a mouth. The anterior side of the larva is defined by the ciliated cleft and the pyriform organ, a typical gymnolaemate larval organ (Figure 1C, D). It is a glandular organ, used for crawling, sensing, and transitorily attaching the larva to the substrate prior to metamorphosis . The pyriform organ extends between the mesodermal blastema and the internal sac (Figure 1B, E). The larva possesses additional nerves, muscles, and yolky cells (Figure 1D) [12, 20, 23]. Some of the internal (yolky) cells are probably endodermal and mesodermal ; however, neither their origin nor their contribution to the adult body plan is resolved. Most of the larval surface is covered by several hundred coronal cells, which constitute the larval swimming organ (Figure 1A, D) .
The terminology that describes the blastemic cell layers in bryozoans is currently confusing and it will be useful to change the nomenclature for the blastemic cells in the future, when comprehensive cell lineage studies become established, which trace the embryological origin of these cells. Herein, we define the cells of the internal sac, the pallial epithelium, the epidermal blastema, and the mesodermal blastema as blastemic cells.
All data outlined in this section and in Figure 2 were gathered from previous studies. Bugula neritina is a brooder and embryos are nourished by a special placenta-like system . The free swimming period of the larva is short, probably between one and 30 hours before settlement . Initial attachment of the larva to the substrate and first morphogenetic movements take only a few minutes and the metamorphosis to the first feeding adult is accomplished within several hours (Figure 2A-I) [18, 21]. The blastemic tissue of the internal sac aids in attachment of the larva to the substrate. The neck region of the internal sac contains dark inclusions aiding in initial attachment before being discarded, while the roof region contains secretory material and forms the permanent attachment disc of the first individual of the colony (ancestrula) [16, 19]. The wall region of the internal sac subsequently forms the epidermis of the free surface of the ancestrula . The epidermal blastema forms the digestive tract and the lophophore of the adult bryozoan . The mesodermal blastema contributes the lophophoral coelomic lining and the splanchnic peritoneum, while the origin of the somatic peritoneum is unclear . The pallial epithelium forms the epithelium of the tentacle sheath of the adult. Transitory larval structures such as the apical organ, the pyriform organ, and the corona are discarded at metamorphosis.
Gene expression patterns in the bryozoan larval stage
Here we describe the spatial expression of 13 developmental genes in the larval stage of Bugula neritina according to gene families.
In the larvae of Bugula neritina, we found BnTropomyosin expression at the upper edge of the apical disc, along the pallial epithelium, surrounding the larva in an equatorial position, as well as alongside the ciliated cleft and internally around the glandular area of the pyriform organ (Figure 3A-F). Additional expression is found in the abapical region below the internal sac (Figure 3A, B). Overall, BnTropomyosin expression is strongest in the coronal cells at the larval surface as well as in cells just below the coronal cells, suggesting the presence of myoepithelial cells in these areas (Figure 3D).
We observed BnHox4 expression in the wall and roof regions of the internal sac in the larvae (Figures 4A-C and 5A, B). Expression is also found internally in the glandular cell complex of the pyriform organ, as well as in the epidermal and mesodermal blastemas (Figure 4B).
SoxB2 and SoxE
BnSoxB2 expression is similar to BnHox4 expression (Figures 4D-F and 5C). In the bryozoan larva, BnSoxE is expressed in a subset of the BnSoxB2 expressing cells. The expression of BnSoxE is within two sickle shaped areas in the mesodermal blastema, which have their greatest extension in posterio-lateral position (Figures 4G-I and 5D). Additional faint expression is observed in the wall region of the internal sac (not shown).
Expression of BnBAMBI is found in a defined area of a trapezoid shape on the posterior side of the bryozoan larva, including cells of the epidermal and mesodermal blastemas, as well as some internal cells. In addition, cells positioned in the abapical wall region show expression, which gives a horseshoe-shaped appearance (Figure 4J-L).
FoxA, FoxB and FoxAB
BnFoxA is expressed in a continuous ring in the epidermal and mesodermal blastemas, as well as in the wall and roof of the internal sac (Figures 6A-C and 7A, B). BnFoxB is co-expressed with BnFoxA in the wall cells of the internal sac, but also shows expression in the circular pallial epithelium (Figure 6D-F). BnFoxAB differs from the expression patterns of BnFoxA and BnFoxAB and is found in epidermal cells at the posterior side of the larvae along the edges of the ciliated groove, and the expression continues from anterior to posterior on the abapical side (Figures 6G-I and 7C, D).
Similar to BnFoxA, BnCdx is expressed as a broad ring in the apical disc in cells of the epidermal and mesodermal blastemas, the wall region of the internal sac, and weak expression in the roof cells (Figure 6J-L).
GATA123 and GATA456
The BnGATA factor 123 shows expression in a ring of cells in the pallial epithelium (Figures 8A-C and 9A, B). Expression is also observed in a few cells of the epidermal blastema and beneath the mesodermal blastema in the glandular area. Weak expression is present in the neck region of the internal sac (Figure 8A, B). In contrast, BnGATA456 is expressed in a single spot on the posterior side of the apical organ in the epidermal blastema of the Bugula neritina larvae (Figure 8D-F).
Wnt1 and Wnt4
BnWnt1 is expressed in cells of the neck region of the internal sac, appearing as a ring at the abapical pole. Further expression is found in cells of the internal sac apical to the neck cells and in the pallial epithelium (Figures 8G-I and 9C, D). BnWnt4 is co-expressed with BnWnt1 in some cells apical of the neck cells (Figure 8J-L). Additional expression is found in the abfrontal half of the pallial epithelium (not shown).
The lecitotrophic larva of Bugula neritina contains only a few truly larval organs, including the apical sense organ and associated nerves, the swimming organ (corona), the glandulo-sensory organ (pyriform organ), internal yolk cells, and muscles. These transitory larval structures are all likely to be discarded at metamorphosis . Blastemic tissues that are known to form the adult during metamorphosis are the epidermal blastema, the mesodermal blastema, the pallial epithelium, and the internal sac. Of the 13 developmental genes that we investigated, 11 genes are expressed in discreet and highly reproducible regions of one or more of the blastemic tissues in the coronate larva that form the adult during metamorphosis. The results indicate that the blastemic cells are probably molecularly pre-determined in the bryozoan larva. Some of the 11 genes show additional expression in the pyriform organ, which is used by the larval stage and resorbed at metamorphosis. Only BnTropomyosin and BnFoxAB are solely expressed in transitory larval structures, the musculature and the pyriform organ (ciliated cleft), respectively (Figures 3, 10). In the following sections, we discuss the gene expression patterns of the bryozoan larval stage in respect to available data of other metazoans.
Genes involved in neural development
Hox4, SoxB, SoxE, BAMBI and FoxB
In the animals investigated so far, Hox genes have been shown to play a crucial role in body plan regionalization along the anterior-posterior axis and are to a great extent expressed in ectodermal and mesodermal derivates. In Acoela, expression of a central Hox gene is found in ectodermal cells in early developmental stages which later give rise to putative neural precursor cells . In several polychaetes, Hox4 is expressed in larval ectoderm and developing neural structures [27–29]. Accordingly, in the mollusk Haliotis, Hox4 expression is in presumptive neuroectoderm and developing ganglia in the trochophore stage and later in the mantle, suggesting partial co-option of the gene for shell formation . Rather similar to Hox4, SoxB genes probably have a conserved role in nervous system development in cnidarians and spiralians [31–36]. In the bryozoan Bugula, BnHox4 and BnSoxB2 are widely co-expressed in both transitory larval structures as well as blastemic cells (Figure 10A, C). Ectodermal expression of Hox4 and SoxB2 in Bugula is similar to their expression in other metazoans, where they are involved in neuroectodermal patterning. Prior to our study, SoxE orthologs have only been investigated in deuterostomes, and in vertebrates they change the fate of neural stem cells into glial stem cells [37–39]. In the bryozoan, BnSoxE is expressed in a small subset of the BnSoxB2 expressing, blastemic cells (Figure 10E). BAMBI is an inhibitor of TGF-β signaling and expression was earlier observed in the beetle Tribolium, where it is co-expressed with BMP/Dpp in the dorsal ectoderm , and in deuterostomes, where BAMBI is co-expressed with Bmp2/4 orthologs [41–44]. BnBAMBI is expressed in a subset of the BnSoxB2 expressing cells, which form parts of the adult bryozoan (Figure 10G). FoxB orthologs play a role in neural development in the cnidarian Nematostella, the ecdysozoan Drosophila, and in chordates [33, 45–48]. In Bugula, BnFoxB is exclusively expressed in future ectodermal tissues (Figure 10K).
The above considerations reveal that the bryozoan orthologs of Hox4, SoxB2, SoxE, BAMBI, and FoxB are to a great extent expressed in limited domains of the blastemas, and the expression patterns parallel that of other animal taxa to some degree. It is feasible that some of the investigated orthologs also play a role in nervous system patterning in bryozoans, and additional studies of gene function will be helpful to evaluate conserved and novel gene functions in the bryozoan life cycle.
Genes involved in gut development
FoxA, Cdx and GATA456
The genes FoxA, Cdx, and GATA456 have previously been shown to be involved in bilaterian gut development. FoxA (Forkhead, HNF3) is important for the development of components of the digestive tract in protostomes as well as in deuterostomes [49–54]. In an acoel, FoxA is expressed in the endoderm surrounding the mouth, suggesting an ancestral role of FoxA in the endoderm and a later co-option in oral ectoderm in bilaterians . In annelids, FoxA is expressed during gut formation [56–58], and in the mollusk Patella, FoxA expression is in the endoderm and in the developing foregut . In the bryozoan larva, BnFoxA is expressed in different blastemas that give rise to adult tissues during metamorphosis (Figure 10I). Expression of BnFoxA in the epidermal blastema, which is supposed to form the adult digestive tract, suggests a function in gut development similar to other metazoans. Cdx/Caudal orthologs are involved in metazoan hindgut formation [60, 61]. In addition, Cdx orthologs are expressed in the brains of both an acoel and an annelid [26, 62]. In the mollusk Patella, Cdx expression is observed in the posterior neuroectoderm and mesodermal cells . In the bryozoan larva, BnCdx is co-expressed with BnFoxA (Figure 10O), and it seems likely that BnCdx is partially involved in adult gut formation as in other lophotrochozoans. However, compared to other taxa, neither FoxA nor Cdx show positional pre-patterning in the bryozoan. The gene GATA456 is involved in endodermal specification in annelids [57, 58]. In the bryozoan larva, BnGATA456 is co-expressed with BnCdx and BnFoxA in a confined domain in the epidermal blastema (Figure 10S).
Our results lead to the conclusion that the genes FoxA, Cdx, and GATA456 are probably involved in the formation of the adult digestive tract in bryozoans and other metazoans and that the epidermal blastema is pre-patterned according to its future fate. In this study, we find no sign of a vestigial larval gut in Bugula neritina, corroborating former histological investigations .
Genes involved in ectoderm specification
The GATA123 factor appears to be involved in ectodermal lineage specification in annelids [57, 58, 64], similar to what has been described for Drosophila. We observed BnGATA123 expression to the greatest extent in ectodermal larval structures as well as in blastemic tissues that form the ectoderm of the adult (Figure 10Q), consistent with the ectodermal expression of the gene in other taxa.
Genes for axial patterning
Wnt signaling might be involved in the development of the primary body axis in most Metazoa . However, from this suggested ancestral role, Wnt diversified and regulates many processes in metazoan taxa. A common feature of Wnt expression is a staggered arrangement along the anterior-posterior axis with partly overlapping domains in, for example, Nematostella, leech, or Capitella[67, 68]. In Bugula, Wnt1 and Wnt4 expression is in specific, partly overlapping subsets of the future adult ectoderm, suggesting that the internal sac is regionalized and its cells might be differentially involved in adult body wall patterning (see also Figure 10U, W).
Gene expression in transitory larval tissues
Tropomyosin and FoxAB
Most of the 11 genes mentioned above are largely expressed in subsets of progenitors of adult tissues, with only BnTropomyosin and BnFoxAB expressed exclusively in transitory larval tissues. Tropomyosin is a general marker of bilaterian musculature. The F-actin component of the larval musculature of Bugula species was investigated in previous studies and revealed prominent longitudinal and radial central muscles positioned between the internal sac and the apical organ, as well as body wall muscles [69, 70]. The BnTropomyosin expression appears to conform with body wall muscles, and corresponds closely with a previous ultrastructural investigation of the Bugula neritina larva, which showed (1) the presence of myoepithelial cells between coronal cells, (2) a collarette of myoepithelial cells that join the infolded pallial epithelium, (3) a pair of myoepithelial cells extending along the lateral sides of the ciliated cleft and (4) a collarette joining the oral margin of the corona  (compare with Figure 3). The lack of Tropomyosin expression in blastemic tissues indicates that the adult musculature differentiates after metamorphosis.
Our study showed that BnFoxAB expression is confined to epidermal cells along the ciliated cleft and along a part of the abapical side of the bryozoan larva (Figure 10M). In a brachiopod larva, FoxAB expression was observed in cells of the apical pole instead (Hejnol et al., unpublished data). Since no such expression is found in the bryozoan larva and comparable gene expression data of other metazoans are currently lacking, our conclusions must be provisional, but recruitment of FoxAB to larval structures may also be found in other taxa.
Pluripotent stem cells versus pre-patterned adult tissues
It was earlier proposed that "homologous set-aside cells" exist in "homologous larvae of protostomes and deuterostomes" . These set-aside cells were supposed to have a rather unlimited division capacity and remain pluripotent and undifferentiated until late embryogenesis. This idea is rooted in the hypothesis that early metazoans were similar to modern larvae, and that the bilaterian adult stage evolved by the innovation of set-aside cells, distinct from the larval cells . However, this hypothesis would presume a similar gene regulatory system in all modern larvae. Here, we have shown that the blastemic cells ("set-aside cells") in the bryozoan larval stage express several metazoan developmental genes. The result indicates that the fate of the blastemic cells is already determined in the bryozoan larval stage and the blastemic cells of bryozoans are probably not pluripotent stem cells or set-aside cells. Our study rather suggests similarities between developmental gene expression in the bryozoan blastemas and adult tissues of other metazoans. In conclusion, this study does not support homology of the bryozoan larval stage and other lophotrochozoan larvae. Instead, the gene expression patterns presented herein indicate that planktonic larvae might have secondarily evolved in bryozoans. The latter hypothesis needs to be evaluated by further studies of gene expression and gene function during the development of bryozoans as well as other lophotrochozoans.
We have shown here the gene expression of 13 developmental genes (Tropomyosin, Hox4, SoxB2, SoxE, BAMBI, Cdx, FoxA, FoxB, FoxAB, GATA123, GATA456, Wnt1, and Wnt4) in the coronate larva of the bryozoan Bugula neritina. Eleven of the 13 genes are expressed in blastemic cells, which are precursors of adult tissues. Probably only a few of the investigated genes have their main function in the larval stage, as for example, BnTropomyosin, which is expressed in larval muscle cells. This study reveals that the blastemic cells in the bryozoan larval stage are most probably pre-patterned according to their future fate in the adult and are not pluripotent, undifferentiated set-aside cells as previously stated . Our results contradict the idea that larval and adult bodies are different entities, but that there is a natural transition from the embryo to the adult with an intermediate larval stage. A comparison of expression patterns among metazoans reveals molecular similarities between bryozoan blastemas and adult tissues of other metazoans. Hence, our study does not indicate homology of the bryozoan larva with other lophotrochozoan larvae, but instead suggests conserved patterns of developmental gene expression amongst lophotrochozoan and metazoan adults. This study adds important data to the fundamental discussion about the evolution of metazoan larval stages and should trigger the interest in investigating gene expression in the "set-aside cells" of the larval stage of other lophotrochozoans as, for example nemerteans. Cell-lineage studies of bryozoans and additional gene expression studies during bryozoan development will also contribute to our understanding of the evolution of metazoan life cycles.
Collection of bryozoans
Colonies of Bugula neritina were collected from submerged hard substrates in shallow water depth (0 to 3 m) in harbors of Honolulu ("La Mariana" and Kewalo Basin) and Pearl Harbor, Oahu, Hawaii in May and June 2009. The colonies were kept in flowing seawater tables in the dark at the Kewalo Marine Laboratory (Hawaii) for a minimum of one day. By exposing colonies to pointed light sources, larvae were released from spawning colonies. They were immediately collected from the water surface under dissecting microscopes and prepared for further investigation.
RNA isolation and cDNA synthesis
Larvae were fixed and stored in RNAlater at 4°C. Larval mRNA was obtained using DynaBeads mRNA DIRECT Kit (Invitrogen, Carlsbad, CA, USA) according to the supplier and stored at -80°C. Complementary DNA (cDNA) synthesis was achieved using the Advantage RT-for-PCR Kit protocol (Clontech Laboratories, Mountain View, CA, USA) following the supplier's instructions. cDNA was stored at -20°C.
The sequences for the genes BnTropomyosin, BnBAMBI, BnSOXB2, BnSoxE BnHox4, BnFoxB, BnFoxAB, WNT1 were gained from a public EST library of Bugula neritin a . Fragments for BnCdx, BnFoxA, BnWNT4 and both GATA factors were gained using degenerate primers with larval cDNA as template. Sequences of the genes from the EST library were amplified using gene specific primers and degenerate fragments were extended using rapid amplification of cDNA ends (RACE) with a SMART RACE cDNA amplification kit (Clontech Laboratories, Mountain View, CA, USA). All fragments were cloned into pGEM-T Easy vectors (Promega Corporation, Madison, WI, USA), transformed into E.coli, and clones sequenced at Macrogen, Inc. (Seoul, South Korea). Fragments obtained from B. neritina were used for probe synthesis in in situ hybridization reactions, which are described below. Primer sequences are listed in Additional file 1 Table S1. Genes were deposited at NCBI GenBank (see below).
Gene accession numbers
BnBAMBI [GenBank: HQ914790]; BnCdx [GenBank: HQ914791]; BnFoxA [GenBank: HQ914792]; BnFoxAB [GenBank: HQ914793]; BnFoxB [GenBank: HQ914794]; BnGATA123 [GenBank: HQ914795]; BnGATA456 [GenBank: HQ914796]; BnGsc [GenBank: HQ914797]; BnHox4 [GenBank: HQ914798]; BnSoxB2 [GenBank: HQ914799]; BnSoxE [GenBank: HQ914800]; BnTropomyosin [GenBank: HQ914801]; BnWnt1 [GenBank: HQ914802]; BnWnt4 [GenBank: HQ914803]; BnWnt8 [GenBank: HQ914804].
Gene orthology assignment
Gene orthologies for all genes (except for BnBAMBI and BnTropomyosin of which the orthology was detected by alignments) were determined by phylogenetic analyses using PhyML 3.0 . Alignments were conducted using MUSCLE  and corrected by hand. ProtTest  was used to determine the best-fitting model. A total of 1,000 to 3,000 bootstraps were calculated respectively, see Additional data files.
In situ hybridization protocol
Larvae of B. neritina were relaxed in 7.14% MgCl2, prefixed in glutaraldehyde fixative (0.3% glutaraldehyde, 3.7% formaldehyde in seawater) for two minutes, fixed in 3.7% formaldehyde for one hour at 4°C, washed five times in PTw (1x PBS = phosphate buffered saline + 0.1% Tween 20) and once in distilled water, which was replaced three times by 100% methanol. The larvae were finally stored in the latter at -20°C. All further steps follow the protocol of Hejnol et al.. Following the development of the in situ probes, larvae were then transferred into 70% glycerol and whole-mounted on glass slides and expression patterns were imaged with a Nikon DXM1200 digital camera mounted on a Nikon Eclipse E1000 microscope. For detailed examination of the expression patterns, larvae were also prepared for histology.
Histology of larva after in situ hybridization
For the histological preparations, in situ hybridized larvae stored in 70% glycerol were washed three times in PBS and dehydrated in an ascending ethanol series with a final step of dehydration in 100% propylene oxide. The larvae were transferred into a 1:1 mixture of 100% propylene oxide and Low Viscosity Resin (LVR) (Agar Scientific, Stansted, UK) over night for infiltration. The animals were embedded in 100% LVR and semithin serial sections (2 μm) were performed on a Leica RM2255 microtome using a Diatome Histo Jumbo Diamond Knife (Diatome, Hatfield, PA, USA). Sections were stained with a 1% solution of basic fuchsine (p-Rosanilin) in 70% ethanol or toluidine blue and embedded in LVR on slides. Imaging was performed with the Nikon setup described in the previous section.
expressed sequence tag
messenger ribonucleic acid
phosphate buffered saline
polymerase chain reaction
low viscosity resin
rapid amplification of cDNA ends
Gordon DP, Taylor PD, Bigey FP: Phylum Bryozoa: moss animals, sea mats, lace corals. Kingdom Animalia: Radiata, Lophotrochozoa, Deuterostomia. New Zealand Inventory of Biodiversity. Edited by: Gordon DP. 2009, Christchurch, New Zealand: Canterbury University Press, 1: 271-279.
Halanych KM, Bacheller JD, Aguinaldo AM, Liva SM, Hillis DM, Lake JA: Evidence from 18S ribosomal DNA that the lophophorates are protostome animals. Science. 1995, 267: 1641-1643. 10.1126/science.7886451.
Dunn CW, Hejnol A, Matus DQ, Pang K, Browne WE, Smith SA, Seaver E, Rouse GW, Obst M, Edgecomb GD, Sørensen MV, Haddock SHD, Schmidt-Rhaesa A, Okusu A, Møbjerg Kristensen R, Wheeler WC, Martindale MQ, Giribet G: Broad phylogenomic sampling improves resolution of the animal tree of life. Nature. 2008, 452: 745-749. 10.1038/nature06614.
Paps J, Baguñà J, Riutort M: Bilaterian phylogeny: a broad sampling of 13 nuclear genes provides a new lophotrochozoan phylogeny and supports a paraphyletic basal Acoelomorpha. Mol Biol Evol. 2009, 26: 2397-2406. 10.1093/molbev/msp150.
Helmkampf M, Bruchhaus I, Hausdorf B: Phylogenomic analyses of lophophorates (brachiopods, phoronids and bryozoans) confirm the Lophotrochozoa concept. Proc R Soc B: Biol Sci. 2008, 275: 1927-1933. 10.1098/rspb.2008.0372.
Helmkampf M, Bruchhaus I, Hausdorf B: Multigene analysis of lophophorate and chaetognath phylogenetic relationships. Mol Phylogenet Evol. 2008, 46: 206-214. 10.1016/j.ympev.2007.09.004.
Hejnol A, Obst M, Stamatakis A, Ott M, Rouse GW, Edgecombe GD, Martinez P, Bagunà J, Bailly X, Jondelius U, Wiens M, Müller WEG, Seaver E, Wheeler WC, Martindale MQ, Giribet G, Dunn CW: Assessing the root of bilaterian animals with scalable phylogenomic methods. Proc R Soc B: Biol Sci. 2009, 276: 4261-4270. 10.1098/rspb.2009.0896.
Paps J, Baguñà J, Riutort M: Lophotrochozoa internal phylogeny: new insights from an up-to-date analysis of nuclear ribosomal genes. Proc R Soc B: Biol Sci. 2009, 276: 1245-1254.
Fuchs J, Obst M, Sundberg P: The first comprehensive molecular phylogeny of Bryozoa (Ectoprocta) based on combined analysis of nuclear and mitochondrial genes. Mol Phylogenet Evol. 2009, 52: 225-233. 10.1016/j.ympev.2009.01.021.
Tsyganov-Bodounov A, Hayward PJ, Porter JS, Skibinski DOF: Bayesian phylogenetics of Bryozoa. Mol Phylogenet Evol. 2009, 52: 904-910. 10.1016/j.ympev.2009.05.010.
Waeschenbach A, Cox CJ, Littlewood DTJ, Porter JS, Taylor PD: First molecular estimate of cyclostome bryozoan phylogeny confirms extensive homoplasy among skeletal characters used in traditional taxonomy. Mol Phylogenet Evol. 2009, 52: 241-251. 10.1016/j.ympev.2009.02.002.
Reed CG: Bryozoa. Marine Invertebrates: Echinoderms and Lophophorates. Edited by: Giese AC, Pearse JS, Pearse VB. 1991, Pacific Grove, CA: Boxwood Press, 6: 86-222.
Nielsen C, Worsaae K: Structure and occurence of cyphonautes larvae (Bryozoa, Ectoprocta). J Morphol. 2010, 271: 1094-1109. 10.1002/jmor.10856.
Zimmer RL, Woollacott RM: Structure and classification of gymnolaemate larvae. Biology of Bryozoans. Edited by: Woollacott R, Zimmer RL. 1977, New York: Academic Press, 57-87.
Wood TS: Development and metamorphosis of cyphonautes larvae in the freshwater ctenostome bryozoan, Hislopia malayensis Annadale, 1916. Bryozoan Studies 2007. 2008, Martinsville, VA: Virginia Museum of Natural History, 339-346.
Reed CG, Cloney RA: The settlement and metamorphosis of the marine bryozoan Bowerbankia gracilis (Ctenostomata: Vesicularioidea). Zoomorphology. 1982, 101: 103-132. 10.1007/BF00312018.
Stricker SA: Settlement and metamorphosis of the marine bryozoan Membranipora membranacea. Bull Mar Sci. 1989, 45: 387-405.
Woollacott RM, Zimmer RL: Metamorphosis of cellularioid bryozoans. Proceedings of the Symposium on Settlement and Metamorphosis of Marine Invertebrate Larvae: American Zoological Society Meeting, Toronto, Ontario, Canada. Edited by: Chia FS, Rice ME. 1978, New York: Elsevier, 49-63.
Zimmer RL, Woollacott RM: Metamorphosis, ancestrulae, and coloniality in bryozoan life cycles. Biology of Bryozoans. Edited by: Woollacott R, Zimmer RL. 1977, New York: Academic Press, 91-139.
Calvet L: Contribution a l'histoire naturelle des Bryozoaires ectoproctes marins. 1900, Coulet et fils, Masson et cie: Montpellier, Paris
Reed CG, Woollacott RM: Mechanisms of rapid morphogenetic movements in the metamorphosis of the bryozoan Bugula neritina (Cheilostomata, Cellularioidea). 1. Attachment to the substratum. J Morphol. 1982, 172: 335-348. 10.1002/jmor.1051720308.
Woollacott RM, Zimmer RL: Attachment and metamorphosis of the cheilo-ctenostome bryozoan Bugula neritina (Linné). J Morphol. 1971, 134: 351-382. 10.1002/jmor.1051340307.
Gruhl A: Ultrastructure of mesoderm formation and development in Membranipora membranacea (Bryozoa: Gymnolaemata). Zoomorphology. 2010, 129: 45-60. 10.1007/s00435-009-0099-3.
Wendt DE: Energetics of larval swimming and metamorphosis in four species of Bugula (Bryozoa). Biol Bull. 2000, 198: 346-356. 10.2307/1542690.
Barrois J: Embryogenie de Bryozoaires; essai d'une théorie générale du développement basée sur l´étude de la metamorphose. J Anat Physiol. 1882, 18: 124-157.
Hejnol A, Martindale M: Coordinated spatial and temporal expression of Hox genes during embryogenesis in the acoel Convolutriloba longifissura. BMC Biol. 2009, 7: 65-10.1186/1741-7007-7-65.
Fröbius AC, Matus DQ, Seaver EC: Genomic organization and expression demonstrate spatial and temporal Hox gene colinearity in the lophotrochozoan Capitella sp. I. PLoS One. 2008, 3: 1-17.
Irvine SQ, Martindale MQ: Expression patterns of anterior Hox genes in the polychaete Chaetopterus: Correlation with morphological boundaries. Dev Biol. 2000, 217: 333-351. 10.1006/dbio.1999.9541.
Kulakova M, Bakalenko N, Novikova E, Cook CE, Eliseeva E, Steinmetz PRH, Kostyuchenko RP, Dondua A, Arendt D, Akam M, Andreeva T: Hox gene expression in larval development of the polychaetes Nereis virens and Platynereis dumerilii (Annelida, Lophotrochozoa). Dev Genes Evol. 2007, 217: 39-54. 10.1007/s00427-006-0119-y.
Hinman VF, O'Brien EK, Richards GS, Degnan BM: Expression of anterior Hox genes during larval development of the gastropod Haliotis asinina. Evol Dev. 2003, 5: 508-521. 10.1046/j.1525-142X.2003.03056.x.
Kerner P, Simionato E, Le Gouar M, Vervoort M: Orthologs of key vertebrate neural genes are expressed during neurogenesis in the annelid Platynereis dumerilii. Evol Dev. 2009, 11: 513-524. 10.1111/j.1525-142X.2009.00359.x.
Le Gouar M, Guillou A, Vervoort M: Expression of a SoxB and a Wnt2/13 gene during the development of the mollusc Patella vulgata. Dev Genes Evol. 2004, 214: 250-256. 10.1007/s00427-004-0399-z.
Magie CR, Pang K, Martindale MQ: Genomic inventory and expression of Sox and Fox genes in the cnidarian Nematostella vectensis. Dev Genes Evol. 2005, 215: 618-630. 10.1007/s00427-005-0022-y.
O'Brien EK, Degnan BM: Expression of POU, Sox, and Pax genes in the brain ganglia of the tropical abalone Haliotis asinina. Mar Biotechnol. 2000, 2: 545-557. 10.1007/s101260000039.
Phochanukul N, Russell S: No backbone but lots of Sox: invertebrate Sox genes. Int J Biochem Cell Biol. 2010, 42: 453-464. 10.1016/j.biocel.2009.06.013.
Taguchi S, Tagawa K, Humphreys T, Satoh N: Group B Sox genes that contribute to specification of the vertebrate brain are expressed in the apical organ and ciliary bands of hemichordate larvae. Zool Sci. 2002, 19: 57-66. 10.2108/zsj.19.57.
Stolt CC, Lommes P, Friedrich RP, Wegner M: Transcription factors Sox8 and Sox10 perform non-equivalent roles during oligodendrocyte development despite functional redundancy. Development. 2004, 131: 2349-2358. 10.1242/dev.01114.
Stolt CC, Lommes P, Sock E, Chaboissier MC, Schedl A, Wegner M: The Sox9 transcription factor determines glial fate choice in the developing spinal cord. Genes Dev. 2003, 17: 1677-1689. 10.1101/gad.259003.
Wegner M, Stolt CC: From stem cells to neurons and glia: a Soxist's view of neural development. Trends Neurosci. 2005, 28: 583-588. 10.1016/j.tins.2005.08.008.
Van der Zee M, da Fonseca RN, Roth S: TGFß signaling in Tribolium: vertebrate-like components in a beetle. Dev Genes Evol. 2008, 218: 203-213. 10.1007/s00427-007-0179-7.
Grotewold L, Plum M, Dildrop R, Peters T, Rüther U: Bambi is co-expressed with Bmp-4 during mouse embryogenesis. Mech Dev. 2001, 100: 327-330. 10.1016/S0925-4773(00)00524-4.
Lowe CJ, Terasaki M, Wu M, Freeman RM, Runft L, Kwan K, Haigo S, Aronowicz J, Lander E, Gruber C, Smith M, Kirschner M, Gerhart J: Dorsoventral patterning in hemichordates: insights into early chordate evolution. PLoS Biol. 2006, 4: 291-10.1371/journal.pbio.0040291.
Onichtchouk D, Chen YG, Dosch R, Gawantka V, Delius H, Massagué J, Niehrs C: Silencing the TGF-ß signalling by the pseudoreceptor BAMBI. Nature. 1999, 401: 480-485. 10.1038/46794.
Yu JK, Satou Y, Holland ND, Shin IT, Kohara Y, Satoh N, Bronner-Fraser M, Holland LZ: Axial patterning in cephalochordates and the evolution of the organizer. Nature. 2007, 445: 613-617. 10.1038/nature05472.
Gamse JT, Sive H: Early anteroposterior division of the presumptive neurectoderm in Xenopus. Mech Dev. 2001, 104: 21-36. 10.1016/S0925-4773(01)00358-6.
Grinblat Y, Gamse J, Patel M, Sive H: Determination of the zebrafish forebrain: induction and patterning. Development. 1998, 125: 4403-4416.
Häcker U, Grossniklaus U, Gehring WJ, Jäckle H: Developmentally regulated Drosophila gene family encoding the fork head domain. Proc Natl Acad Sci USA. 1992, 89: 8754-8758. 10.1073/pnas.89.18.8754.
Mazet F, Shimeld SM: The evolution of chordate neural segmentation. Dev Biol. 2002, 251: 258-270. 10.1006/dbio.2002.0831.
Gaudet J, Mango SE: Regulation of organogenesis by the Caenorhabditis elegans FoxA protein PHA-4. Science. 2002, 295: 821-825. 10.1126/science.1065175.
Kormish JD, Gaudet J, McGhee JD: Development of the C. elegans digestive tract. Curr Opin Genet Dev. 2010, 20: 346-354. 10.1016/j.gde.2010.04.012.
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-158. 10.1016/j.ydbio.2010.01.018.
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-4181. 10.1242/dev.02577.
Taguchi S, Tagawa K, Humphreys T, Nishino A, Satoh N, Harada Y: Characterization of a hemichordate fork head/ HNF-3 gene expression. Dev Genes Evol. 2000, 210: 11-17. 10.1007/PL00008181.
Weigel D, Jürgens G, Kuttner F, Seifert E, Jäckle H: The homeotic gene fork head encodes a nuclear protein and is expressed in terminal regions of the Drosophila embryo. Cell. 1989, 57: 645-658. 10.1016/0092-8674(89)90133-5.
Hejnol A, Martindale M: Acoel development indicates the independent evolution of the bilaterian mouth and anus. Nature. 2008, 456: 382-386. 10.1038/nature07309.
Arenas-Mena C: Embryonic expression of HeFoxA1 and HeFoxA2 in an indirectly developing polychaete. Dev Genes Evol. 2006, 216: 727-736. 10.1007/s00427-006-0099-y.
Boyle MJ, Seaver EC: Developmental expression of foxA and gata genes during gut formation in the polychaete annelid, Capitella sp. I. Evol Dev. 2008, 10: 89-105. 10.1111/j.1525-142X.2007.00216.x.
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-10.1186/2041-9139-1-2.
Lartillot N, Le Gouar ML, 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-561. 10.1007/s00427-002-0274-8.
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-584. 10.1111/j.1525-142X.2005.05061.x.
Wu LH, Lengyel JA: Role of caudal in hindgut specification and gastrulation suggests homology between Drosophila amnioproctodeal invagiation and vertebrate blastopore. Development. 1998, 125: 2433-2442.
Fröbius AC, Seaver EC: ParaHox gene expression in the polychaete annelid Capitella sp. I. Dev Genes Evol. 2006, 216: 81-88. 10.1007/s00427-005-0049-0.
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-37. 10.1016/S1567-133X(02)00091-1.
Gillis WJ, Bowerman BA, Schneider SQ: Ectoderm- and endomesoderm-specific GATA transcription factors in the marine annelid Platynereis dumerilii. Evol Dev. 2007, 9: 39-50. 10.1111/j.1525-142X.2006.00136.x.
Heitzler P, Haenlin M, Ramain P, Calleja M, Simpson P: A genetic analysis of pannier, a gene necessary for viability of dorsal tissues and bristle positioning in Drosophila. Genetics. 1996, 143: 1271-1286.
Petersen CP, Reddien PW: Wnt Signaling and the polarity of the primary body axis. Cell. 2009, 139: 1056-1068. 10.1016/j.cell.2009.11.035.
Cho S, Vallès Y, Giani VC, Seaver EC, Weisblat DA: Evolutionary dynamics of the wnt gene family: a lophotrochozoan perspective. Mol Biol Evol. 2010, 27: 1645-1658. 10.1093/molbev/msq052.
Kusserow A, Pang K, Sturm C, Hrouda M, Lentfer J, Schmidt HA, Technau U, von Haeseler A, Hobmayer B, Martindale MQ, Holstein TW: Unexpected complexity of the Wnt gene family in a sea anemone. Nature. 2005, 433: 156-160. 10.1038/nature03158.
Gruhl A: Muscular systems in gymnolaemate bryozoan larvae (Bryozoa: Gymnolaemata). Zoomorphology. 2008, 127: 143-159. 10.1007/s00435-008-0059-3.
Santagata S: The morphology and evolutionary significance of the ciliary fields and musculature among marine bryozoan larvae. J Morphol. 2008, 269: 349-364. 10.1002/jmor.10592.
Peterson KJ, Cameron RA, Davidson EH: Set-aside cells in maximal indirect development: evolutionary and developmental significance. BioEssays. 1997, 19: 623-631. 10.1002/bies.950190713.
Raff RA: Origins of the other metazoan body plans: the evolution of larval forms. Philos Trans R Soc Lond B: Biol Sci. 2008, 363: 1473-1479. 10.1098/rstb.2007.2237.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1197. 10.1093/nar/gkh340.
Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-2105. 10.1093/bioinformatics/bti263.
Reed CG, Woollacott RM: Mechanisms of rapid morphogenetic movements in the metamorphosis of the bryozoan Bugula neritina (Cheilostomata, Cellularioidea): II. The role of dynamic assemblages of microfilaments in the pallial epithelium. J Morphol. 1983, 177: 127-143. 10.1002/jmor.1051770202.
The authors thank the personnel of the M. Q. Martindale lab and Michiel Boekhout for help in the molecular laboratory, as well as the support staff of the Kewalo Marine Laboratory. We are grateful to Brian Nedved and Elizabeth Perotti for help with species collection and to Christiane Todt for help with the histological preparations. Furthermore, we wish to thank three anonymous reviewers for their valuable comments on the manuscript. The work was made possible with Swedish travel grants from the Adlerbertska Stipendium Stiftelse and the Helge Ax:son Johnsons Stiftelse to JF. The research was funded by NSF AToL grant (EF05-31558), NASA Exobiology and the DFG (HE5183/2-2).
The authors declare that they have no competing interests.
JF and AH performed the molecular laboratory work and JF conducted histology and imaging. AH coordinated the project and performed the phylogenetic analyses. JF wrote the draft manuscript and AH and MQM participated in manuscript preparation. All authors read and approved the final version of the manuscript.
Electronic supplementary material
Additional file 1: Primer sequences and gene orthology analyses. Degenerate primer sequences (Table S1) and trees of the orthology analyses of the Bugula neritina genes BnSoxB2, BnSoxE, BnFoxA, BnFoxB, BnFoxAB, BnHox4, BnCdx, BnGATA123, BnGATA456, BnWnt1, BnWnt4, and BnWnt8 (Figures S1 - S5). (PDF 818 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.