www.nature.com/scientificreports OPEN Expression of segment polarity genes in brachiopods supports a non-segmental ancestral role of received: 10 February 2016 engrailed for bilaterians Accepted: 09 August 2016 Published: 26 August 2016 Bruno C. Vellutini & Andreas Hejnol The diverse and complex developmental mechanisms of segmentation have been more thoroughly studied in arthropods, vertebrates and annelids—distantly related animals considered to be segmented. Far less is known about the role of “segmentation genes” in organisms that lack a segmented body. Here we investigate the expression of the arthropod segment polarity genes engrailed, wnt1 and hedgehog in the development of brachiopods—marine invertebrates without a subdivided trunk but closely related to the segmented annelids. We found that a stripe of engrailed expression demarcates the ectodermal boundary that delimits the anterior region of Terebratalia transversa and Novocrania anomala embryos. In T. transversa, this engrailed domain is abutted by a stripe of wnt1 expression in a pattern similar to the parasegment boundaries of insects—except for the expression of hedgehog, which is restricted to endodermal tissues of the brachiopod embryos. We found that pax6 and pax2/5/8, putative regulators of engrailed, also demarcate the anterior boundary in the two species, indicating these genes might be involved in the anterior patterning of brachiopod larvae. In a comparative phylogenetic context, these findings suggest that bilaterians might share an ancestral, non-segmental domain of engrailed expression during early embryogenesis. Annelids (e.g. earthworms), arthropods (e.g. insects) and vertebrates (e.g. humans) are bilaterally symmet- ric animals that have a major part of their adult body organized into segments—sophisticated morphological units repeated along the anteroposterior axis1. The developmental mechanisms of segmentation are diverse and complex, but arthropods and vertebrates do share some molecular similarities in the patterning of their body segments2–7. These findings stimulated a debate about the evolution of segmentation that resulted in two conflict- ing evolutionary hypotheses. Either the similarities are the outcome of common descent, and thus support the homology of bilaterian body segments2,8–10, or they represent evolutionary convergences11–13. These hypotheses, however, are based on the examination of distantly related groups. Annelids, arthropods and vertebrates belong to distinct branches of the Bilateria, Spiralia, Ecdysozoa and Deuterostomia, respectively (Fig. 1a). This sole comparison over long evolutionary distance can be misleading to reconstruct the evolution of segmentation, because the ancestral conditions within more closely related taxa remain unknown. Thus, dis- tinguishing homology from convergence requires a robust phylogenetic framework and dense taxonomic sam- pling11,14. To better understand the evolution of segmentation, it is necessary to investigate the related groups of annelids, arthropods and vertebrates. Most of these lineages do not have body segments, but can display a variety of serially repeated structures, such as neural ganglia, excretory organs, coeloms, and others1,12,15. Since segmen- tation is not an “all-or-nothing” trait15,16, further data on the morphology and gene networks in these lineages can help to elucidate the evolution of the developmental mechanisms of segmentation. A crucial mechanism of segmentation is the ability to establish tissue boundaries during animal develop- ment17. One of the best studied developmental boundaries are the parasegmental boundaries of the fruit fly Drosophila melanogaster. Fly segmentation occurs by the simultaneous division of the trunk into fourteen parasegments, developmental regions that lie out of register with the larval segments18. At the molecular level, the parasegment boundary is established by the abutting domains of engrailed (en) and wingless (wnt1)19, and is maintained by a positive feedback loop mediated by hedgehog (hh)20. Therefore, en, wnt1 and hh are commonly Sars International Centre for Marine Molecular Biology, University of Bergen, Thormøhlensgate 55, 5006 Bergen, Norway. Correspondence and requests for materials should be addressed to A.H. (email: [email protected]) Scientific RepoRts | 6:32387 | DOI: 10.1038/srep32387 1 www.nature.com/scientificreports/ Figure 1. Phylogenetic relationships between bilaterally symmetric animals and the morphological features of larval brachiopods. (a) Phylogenetic tree based on recent data29,60,61,74 highlighting the three major bilaterian clades, Deuterostomia, Ecdysozoa and Spiralia. Groups traditionally regarded as truly segmented are marked in bold. Brachiopoda is highlighted in orange. (b) Ectodermal and mesodermal larval traits of the brachiopods T. transversa and N. anomala, as representatives of the Rhynchonelliformea and Craniiformea47. referred to as segment polarity genes21–23, and their expression and function is conserved among the body seg- ments of other arthropods7,24–26. The expression of segment polarity genes in annelids is more variable and, in general, do not suggest a role in the patterning of body segments5,27. So far, the only exception is the annelid Platynereis dumerilii, where the expression of en, wnt1 and hh correlates with the segmental boundaries10,28. Thus, to better understand the ancestral developmental functions of the segment polarity genes in the spiralians, it is important to investigate gene expression in groups more closely related to annelids. Brachiopods are bivalved marine invertebrates with close relations to other lophophorates, nemerteans, anne- lids and molluscs29 (Fig. 1a). Even though adult brachiopods do not have body segments, the larval stages show Scientific RepoRts | 6:32387 | DOI: 10.1038/srep32387 2 www.nature.com/scientificreports/ Figure 2. Developmental stages of the trilobed larva of T. transversa (a–e) and the bilobed larva of N. anomala (f–j). Each panel shows a ventral view (top) and lateral view (bottom) of a maximum intensity projection of embryos stained with DAPI. Anterior is top in all panels and ventral is to the right in all lateral views. White arrowheads mark the apical/mantle boundary and black arrowheads mark the mantle/pedicle boundary. al: apical lobe, ml: mantle lobe, pl: pedicle lobe. Scale bars = 20 μ m. putative repeated structures in ectodermal and mesodermal tissues. Externally, brachiopod larvae exhibit two to four lobes along the anteroposterior axis that are divided by transverse ectodermal boundaries. These lobes were once homologized to annelid segments30, but the observation that the ectodermal boundaries do not involve the underlying mesoderm weakened this hypothesis31. The mesoderm morphology is variable between brachiopod species, it can be partitioned into two or four pairs of coelomic sacs32,33 or not subdivided. The discovery of a brachiopod larvae with serially arranged coelomic sacs33 has revived the idea that brachiopods had a segmented ancestor34,35. Therefore, these putative segmented structures, and the closer phylogenetic position to annelids, place brachiopods as an interesting group to test the involvement of the segment polarity genes in other develop- mental boundaries, and to better comprehend the evolution of segmentation mechanisms in protostomes. To investigate the developmental and molecular features of the boundaries in brachiopod embryos, we stud- ied the trilobed larva of Terebratalia transversa36, and the bilobed larva of Novocrania anomala33,37, species that belong to distinct brachiopod lineages (Fig. 1b). In these two species, we analyzed the expression of the seg- ment polarity genes en, wnt1, and the core components of the Hedgehog signaling pathway to test whether their expression correlate with the development of the ectodermal and mesodermal boundaries of T. transversa and N. anomala larvae. Furthermore, we examined upstream control genes of en and discovered similarities between the molecular profile of a brachiopod larval boundary and the embryonic patterning of anterior boundaries in deu- terostomes, such as the hemichordate collar/trunk boundary and the vertebrate fore/midbrain boundary, further suggesting a non-segmental ancestral role of en for bilaterians. Results Ectodermal and mesodermal boundaries of larval brachiopods. In T. transversa, gastrulation occurs by invagination and results in a radially symmetric gastrula that elongates in the anteroposterior axis (Fig. 2a–c). A transverse ectodermal furrow (indentation of the epidermis) appears above the midline of the bilateral gastrula stage, extending from the blastopore lips to the dorsal side (Fig. 2c), with epidermal cells posterior to the inden- tation being more elongated (Supplementary Fig. S1). The ectodermal furrow and cell shape differences at the bilateral gastrula are the first morphological manifestation of the boundary between the future apical and mantle Scientific RepoRts | 6:32387 | DOI: 10.1038/srep32387 3 www.nature.com/scientificreports/ lobes in T. transversa (apical/mantle boundary). At subsequent developmental stages, this furrow deepens and clearly divides the apical lobe from the remainder of the embryo (Fig. 2d). T. transversa exhibits a third lobe at the posterior end, the pedicle lobe (Fig. 2e). The boundary between the mantle and pedicle lobes (mantle/pedicle boundary) is identified by the narrowing of the posterior portion of the embryo (Fig. 2d) and by the subsequent outgrowth of the mantle lobe (Fig. 2e). The morphology of the man- tle/pedicle boundary differs from the apical/mantle boundary because it is demarcated by an ectodermal fold, rather than by an indentation furrow. For this reason, the ectodermal boundaries of the trilobed brachiopod larva are not repeated, but unique structures of the larval body. In N. anomala, the apical and mantle lobes are also demarcated by an ectodermal furrow, as seen in T. transversa (Fig. 2f–i), but its larva does not form a pedicle lobe (Fig. 2i). Thus, the morphology of the apical/mantle boundary is conserved between both species of brachiopods. The mesoderm of T. transversa has prominent projections associated to the chaetae sacs in the mantle lobe, but there are no partitions (Supplementary Fig. S2). In contrast, the mesoderm of N. anomala has lateral constrictions that individualize four pairs of coelomic sacs, but the posterior pouches are also tightly associated to the three pairs of dorsal chaetae sacs, while remaining interconnected more medially in the ventral side (Supplementary Fig. S2). Expression of the segment polarity genes hh, en and wnt1. To compare the ectodermal and mes- odermal boundaries of brachiopods with the segment boundaries of arthropods and annelids, we analyze the expression of the segment polarity genes en, wnt1 and hh during the embryonic development of T. transversa and N. anomala. Expression of hh localizes to the blastoporal lip and invaginating endomesoderm during early gas- trulation, and is restricted to the endoderm in later stages (Figs 3a and 4a). In N. anomala, we detected an addi- tional transverse ventral domain of hh near the animal pole that disappears in the bilobed stage (Fig. 4a). Since the Hedgehog ligand can signal across embryonic layers, we further analyzed the expression of the Hedgehog recep- tors ptc and smo and the transcription factor gli. Transcripts are expressed in the mesoderm of T. transversa and N. anomala, as well as in the ectodermal apical and mantle territories (Supplementary Figs S3 and S4). Altogether, the expression of hh and downstream pathway components does not correlate spatially or temporally with the development of the larval ectodermal boundaries of T. transversa and N. anomala. We first detect transcripts of en in T. transversa in the radial gastrula stage. The expression forms an almost-complete ring in the midregion of the embryo and, progressively fades from the anterior and posterior end (Supplementary Fig. S5), leaving a pair of lateral domains of en in the asymmetric gastrula stage (Fig. 3b). These domains extend ventrally and dorsally without reaching the blastopore on the ventral side of the bilateral gastrula stage (Fig. 3b). Transcripts of en border the apical/mantle furrow in the subsequent stage and fade in the trilobed larvae (Fig. 3b). In N. anomala a single dorsal ectodermal stripe of en expression occurs in the radial gastrula (Supplementary Fig. S5) and a second stripe emerges at the posterior end of the asymmetric gastrula (Fig. 4e). Similar to T. transversa, en stripes extend ventrally without reaching the embryo midline and fade in the developed larva (Fig. 4e). Therefore, expression of en is consistent between the two species, with the apical/mantle furrow forming immediately anterior to a pair of en stripes. In addition to the striped domains demarcating the apical/mantle boundary, we identified two more posterior en territories on the ventral and dorsal regions of larval brachiopods, localized in the anterior-most region of the pedicle lobe. One is located on the ventral side of T. transversa and consists of a pair of en ectodermal bands associated with the blastopore in the bilateral gastrula stage (Fig. 3b). The other is located on the opposite side, at the center of the dorsal surface, as a triangular-shaped domain of en expression (Supplementary Fig. S5). A corre- spondent domain also occurs on the dorsal side of N. anomala and localizes to the shell rudiment in the bilobed larva stage (Supplementary Fig. S5). Mesodermal expression of en initiates later than the ectodermal domains. In T. transversa, transcripts of en form a pair of bands in the pedicle mesoderm that are connected to the ectodermal domains, while the mesoder- mal expression in N. anomala localizes directly inner to the lateral ectodermal domains of en (Supplementary Fig. S6). Expression of en is located to the posterior portion of the second and third pair of coelomic sacs, but not in the first or fourth (Supplementary Fig. S6). Overall, the mesodermal expression of en occurs in contiguity with the preceding ectodermal domains in both species. Expression of wnt1 in the asymmetric gastrula of T. transversa is associated with the posterior portion of the blastopore (Fig. 3c) and shows no spatial correlation to the lateral en domains (Fig. 5a). At the bilateral gastrula, a transverse pair of stripes of wnt1 expression appears in the apical lobe of T. transversa, bordering the apical/man- tle furrow anteriorly (Fig. 3c). The expression abuts en at the apical/mantle boundary without overlap (Fig. 5b). As the apical/mantle furrow deepens in the bilobed larva, wnt1 and en transcripts resolve into well-defined non-overlapping stripes in T. transversa (Fig. 5c). The stripes demarcate precisely the morphological furrow at the apical/mantle boundary of T. transversa with wnt1 positioned anterior and en positioned posterior to the fur- row (Fig. 5c). In contrast, the posterior domain of wnt1 develops a tight expression overlap with the pedicle lobe domain of en, on the ventral side of T. transversa (Fig. 5b,c). Such precise coexpression pattern also occurs on the dorsal side, when wnt1 transcripts are first detected within the triangle-shaped territory of en (Fig. 5c). N. anomala does not exhibit wnt1 domains in the apical lobe. The transcripts are exclusively at the posterior end, associated with the blastopore (Fig. 4c). We did not detect any coexpression sites for en and wnt1 in N. anom- ala. Thus, the expression of wnt1 at the apical/mantle boundary is not consistent between brachiopods, and only T. transversa has wnt1 stripes demarcating the morphology of the furrow. Expression of the putative engrailed regulators pax6, pax2/5/8 and fgf8/17/18. The only gene consistently related to the apical/mantle boundary in both brachiopod larvae is en. Its early expression, however, does not encircle the circumference of the embryo (Fig. 3b), suggesting that there must be upstream factors positioning the apical/mantle boundary. Regulators of en known from D. melanogaster segmentation, such as even-skipped or sloppy paired38,39, do not correlate to the en domains at the apical/mantle boundary of T. Scientific RepoRts | 6:32387 | DOI: 10.1038/srep32387 4 www.nature.com/scientificreports/ Figure 3. Whole mount in situ hybridization of hh, en and wnt1 orthologs in representative developmental stages of the brachiopod T. transversa. Anterior is top in all panels and ventral is to the right in all lateral views. (a) Expression of hh. (b) Expression of en. (c) Expression of wnt1. Line drawings represent expression at the bilateral gastrula stage. White arrowheads mark the apical/mantle boundary and black arrowheads mark the mantle/pedicle boundary. Double arrowheads mark the ectodermal expression of en. Asterisk marks area of unspecific staining due to the rudiment of the larval shell (b, trilobed). cs: chaetae sacs. Scale bars = 20 μ m. transversa40. Therefore, we selected further candidate genes related to the patterning of the vertebrate neuroec- toderm, where en plays a crucial role in the establishment of boundaries41. In the mid/hindbrain boundary, en is upregulated by pax2 in a positive feedback loop, represses and is repressed by pax6, and, together with fgf8, is essential to maintain the boundary between the diencephalon and mesencephalon42–45. Thus, we analyzed the expression of the correspondent brachiopod orthologs pax6, pax2/5/8 and fgf8/17/18 in relation to the expression of en and the development of the apical/mantle boundary. Transcripts of pax6 form a pair of broad anterior domains that occupy the entire anterior half of T. transversa radial gastrula, and localize to the apical lobe in subsequent stages (Fig. 6a). In N. anomala, pax6 expression ini- tiates in the radial gastrula as an apical ring and also localizes to the apical lobe in later stages (Fig. 7a). In both species the posterior limit of pax6 expression borders the apical/mantle furrow. T. transversa pax2/5/8 expression begins as a posterior half-ring in the radial gastrula (Fig. 6b) that overlaps with the pax6 territory at the future apical/mantle boundary (Fig. 8a). The transcripts cover the whole exten- sion of the mantle ectoderm in the asymmetric gastrula, except for the chaetae sac primordia and blastoporal lips (Fig. 6b), and contain the lateral domains of en (Fig. 8d). The anterior limit of pax2/5/8 expression crosses the apical/mantle furrow (Fig. 6b) and overlaps with pax6 expression at the posteriormost region of the apical lobe (Fig. 8b,c). We corroborated this anterior limit of pax2/5/8 transcripts by the combined analysis with en expression (Fig. 8e,f); see Supplementary Fig. S7 for additional data. In N. anomala, pax2/5/8 is expressed in a Scientific RepoRts | 6:32387 | DOI: 10.1038/srep32387 5 www.nature.com/scientificreports/ Figure 4. Whole mount in situ hybridization of hh, en and wnt1 orthologs in representative developmental stages of the brachiopod N. anomala. Anterior is top in all panels and ventral is to the right in all lateral views. (a) Expression of hh. (b) Expression of en. (c) Expression of wnt1. Line drawings represent expression at the bilateral gastrula stage. White arrowheads mark the apical/mantle boundary. Scale bars = 20 μ m. posterodorsal territory in the mantle lobe, limited by the apical lobe (Fig. 7b). Different than T. transversa, we observe mesodermal expression of pax2/5/8 in two pairs of domains between the chaetae sacs of N. anomala (Fig. 7b). In summary, pax6 and pax2/5/8 are expressed in complementary domains, with a narrow overlap, that encircle the whole embryo before the lateral patches of en form transverse stripes. Thus, the apical/mantle furrow in both species is demarcated by the posterior limit of pax6 expression with abutting domains of en in the mantle lobe. Transcripts of fgf8/17/18 form a pair of transverse ventral bands in the apical lobe of T. transversa asymmetric gastrula (Fig. 6c). These bands do not border the stripes of en posterior to the furrow (Supplementary Fig. S7). In addition, developing chaetae sacs express fgf8/17/18 (Fig. 6c) and are cleared from en and pax2/5/8 expression (Supplementary Fig. S7). In N. anomala, fgf8/17/18 expression initiates in the asymmetric gastrula with a sim- ilar anterior ventral band and a posterior domain encircling blastopore (Fig. 7c). This posterior demarcates the anterior-most region of the mantle lobe at the bilateral gastrula but fades in the bilobed larva (Fig. 7c). Similarly to T. transversa, chaetae primordia express fgf8/17/18 (Fig. 7c). Thus, the expression of fgf8/17/18 is mainly asso- ciated with the chaetae sacs in both species and does not show a clear correlation with ectodermal boundaries. Over-activation of the Wnt pathway in T. transversa. The pervading role of Wnt signaling in the axial patterning of metazoans46, led us to investigate the role of the canonical Wnt pathway in the development of T. transversa. We found that the over-activation of Wnt signaling prevents the formation of the apical/mantle furrow (Supplementary Fig. S8). The posterior larval structures are expanded, the mantle lobe is not formed, and the Scientific RepoRts | 6:32387 | DOI: 10.1038/srep32387 6 www.nature.com/scientificreports/ Figure 5. Whole mount double fluorescent in situ hybridization with en (vermillion) and wnt1 (green) in T. transversa. (a) Asymmetric gastrula. Same specimens of Fig. 2b. (b) Bilateral gastrula. Same specimens of Fig. 2c. (c) Bilobed larva. Ventral view is the same specimen of Fig. 2d. Ventral and lateral views are maximum intensity projections. Optical sections show additional counter staining with DAPI (gray). Merge (yellow) highlights the areas of en and wnt1 co-expression. Scale bars = 20 μ m. expression of en, wnt1 and pax6 shifts more anteriorly (Supplementary Fig. S9). Early treatments entirely abolish the anterior domains of en and wnt1, but do not suppress the expression of pax6 (Supplementary Fig. S9). Discussion The evolutionary lineages of T. transversa (Rhynchonelliformea) and N. anomala (Craniiformea) diverged at least 500 million years ago47. It has been hypothesized that planktotrophy is ancestral for the Brachiopoda48 and that larvae of rhynchonelliforms and craniiforms might have evolved a lecithotrophic mode of development inde- pendently49,50. However, there are no extant planktotrophic larvae in the rhynchonelliforms, and for this reason, it is not possible to infer if the planktotrophy found in the larvae of extant linguliforms—the sister group of the craniiforms—is ancestral or derived. In our study, we compared species from the rhynchonelliforms and cranii- forms to determine which traits are conserved or derived. We found that the trilobed larva of T. transversa and the bilobed larva of N. anomala share an apical/mantle boundary with a conserved ectodermal furrow morphology. A similar ectodermal furrow also delimits the apical lobe of the planktotrophic larva of Lingula anatina51, suggest- ing the apical/mantle boundary is an ancestral feature of brachiopod embryogenesis. We found that T. transversa and N. anomala larvae show surprisingly consistent patterns of gene expression. Most expression domains in a species have a reciprocal, similarly positioned territory in the other species. In the trilobed larva of T. transversa, for example, en is expressed in a pair of anterior domains that correlate closely with the apical/mantle boundary, and in a pair of posterior domains bordering the mantle/pedicle boundary. The bilobed larva of N. anomala not only exhibits the correspondent en domains lining its apical/mantle boundary, Scientific RepoRts | 6:32387 | DOI: 10.1038/srep32387 7 www.nature.com/scientificreports/ Figure 6. Whole mount in situ hybridization of pax6, pax2/5/8 and fgf8/17/18 orthologs in representative developmental stages of the brachiopod T. transversa. (a) Expression of pax6. (b) Expression of pax2/5/8. (c) Expression of fgf8/17/18. Anterior is top in all panels and ventral is to the right in all lateral views. Line drawings represent expression at the bilateral gastrula stage. White arrowheads mark the apical/mantle boundary and black arrowheads mark the mantle/pedicle boundary. Scale bars = 20 μ m. Scientific RepoRts | 6:32387 | DOI: 10.1038/srep32387 8 www.nature.com/scientificreports/ Figure 7. Whole mount in situ hybridization of pax6, pax2/5/8 and fgf8/17/18 orthologs in representative developmental stages of the brachiopod N. anomala. (a) Expression of pax6. (b) Expression of pax2/5/8. (c) Expression of fgf8/17/18. Anterior is top in all panels and ventral is to the right in all lateral views. Line drawings represent expression at the bilateral gastrula stage. White arrowheads mark the apical/mantle boundary. Scale bars = 20 μ m. Scientific RepoRts | 6:32387 | DOI: 10.1038/srep32387 9 www.nature.com/scientificreports/ Figure 8. Whole mount double fluorescent in situ hybridization of pax6 + pax2/5/8 (a–c) and en + pax2/5/8 (d–f) in the brachiopod T. transversa. (a,d) Asymmetric gastrula in lateral view. (b,e) Bilateral gastrula in ventral view. (c,f) Bilobed larva in ventral view. Anterior is top in all panels and ventral is to the right in all lateral views. Side panels show the relation between gene expression and the apical/mantle boundary (striped yellow line). White arrowheads mark the apical/mantle boundary and black arrowheads mark the mantle/ pedicle boundary. Scale bars = 20 μ m. but also the reciprocal posterior pair. We do not interpret these domains in N. anomala as a vestigial expression of en. First, in known cases of vestigial expression, an incipient aspect of the morphology is still present52,53. But there is no morphological evidence of a differentiated pedicle lobe in N. anomala. Second, the evolutionary origin of the pedicle lobe remains unresolved and it is unclear if it is an ancestral feature of brachiopod larvae. Nonetheless, the consistent expression patterns between T. transversa and N. anomala allow us to suggest that the presence of two pairs of lateral en domains is the ancestral condition for brachiopods, independent of the trilobed or bilobed larval morphology. In order to better understand the evolution of en expression across bilaterians, we compared the brachiopod data with the available data from other animals (Fig. 9). We found that most domains correlate to some kind of epithelial boundary such as segment boundaries of annelids, shell borders in molluscs, arthropod segment boundaries, seam cells in nematodes, collar/trunk boundary in hemichordates, brain boundaries in vertebrates and somite boundaries in cephalochordates (Supplementary Table S1). A closer look reveals that at its earliest instance of embryonic expression, most groups have an ectodermal expression of en in laterodorsal domains (Fig. 9 and Supplementary Table S1). In annelids and brachiopods, a pair of lateral en domains are located at the head/trunk and apical/mantle boundary, respectively. In D. melanogaster the first stripe of en occurs at the cephalic furrow54 and in hemichordates en is expressed at the collar/trunk boundary55. This comparative data suggests that en expression is mainly ectodermal and might have been associated to the interface between an anterior (e.g. head) and a posterior region of the embryo (e.g. trunk), a division shared by most protostomes and deuterostomes. In addition, our data indicates that the early patterning of the apical/mantle boundary of brachiopods is char- acterized by the overlapping expression of pax6 and pax2/5/8 (Fig. 10a). Interestingly, the complementary pat- terns of pax6 and pax2/5/8 expression during early embryogenesis occur in the ectoderm of hemichordates55, cephalochordates56,57 and in the neuroectoderm of vertebrates43,44 (Fig. 10b), suggesting that deuterostomes share a similar arrangement of gene products along the anteroposterior axis58. No nervous structures have been described at the apical/mantle boundary of larval brachiopods, so far40,59. Thus, the similarity to the brachiopod pattern could indicate that such genes were also involved in the axial patterning of the last common protos- tome ancestor. Or simply, that they have been independently coopted to these structures. Unfortunately, data on the expression of pax6 and pax2/5/8, specially during early embryogenesis, are still scarce outside arthropods (Supplementary Table S1). Thus, it remains to be determined if the expression patterns of these genes in brachio- pods are ancestral to the protostomes. Overall, the spatial arrangement and temporal expression patterns of our candidate genes suggest the apical/ mantle boundary of brachiopods is not patterned by a segment polarity mechanism (Supplementary Fig. S10). Nevertheless, we show that abutting stripes of en and wnt1 are not exclusive for the segment and parasegment boundaries of P. dumerilii28 and arthropods7,19,24–26, respectively—it also occurs in a non-segmental boundary of a larval brachiopod. Since the expression of segment polarity genes in annelids is variable and does not support a role in segment formation5, and recent phylogenies suggest the protostome ancestor did not have a segmented trunk29,60, the segment polarity genes might be an arthropod-specific feature. Therefore, we interpret the adjacent Scientific RepoRts | 6:32387 | DOI: 10.1038/srep32387 10