Ancestral role of Pax2/5/8 in molluscan brain and multimodal sensory system development
BMC Evolutionary Biology volume 15, Article number: 231 (2015)
Mollusks represent the largest lophotrochozoan phylum and exhibit highly diverse body plans. Previous studies have demonstrated that transcription factors such as Pax genes play important roles during their development. Accordingly, in ecdysozoan and vertebrate model organisms, orthologs of Pax2/5/8 are among others involved in the formation of the midbrain/hindbrain boundary, the auditory/geosensory organ systems, and the excretory system.
Pax2/5/8 expression was investigated by in situ hybridization during the development of representatives of the two major molluscan subclades, Aculifera and Conchifera.
Compared to the investigated polyplacophoran and bivalve species that lack larval statocysts as geosensory organs and elaborate central nervous systems (CNS), cephalopods possess highly centralized brains and statocysts. Pax2/5/8 is expressed in regions where sensory cells develop subsequently during ontogenesis. Expression domains include esthetes and the ampullary system in polyplacophorans as well as the eyes of cephalopods. No Pax2/5/8 expression was observed in the less centralized CNS of bivalve, polyplacophoran, and gastropod embryos, thus arguing for a loss of Pax2/5/8 involvement in CNS development in these lineages. In contrast, Pax2/5/8 is expressed among others in brain lobes along the trajectory of the esophagus that divides the cephalopod brain.
Our results, along with those on Otx- and Hox-gene expression, demonstrate that the cephalopod condition is similar to that in mouse and fruit fly, with Otx being expressed in the anterior-most brain region (except for the vertical lobe) and a Pax2/5/8 expression domain separating the Otx-domain from a Hox-gene-expressing posterior brain region. Thus, Pax2/5/8 appears to have been recruited independently into regionalization of non-homologous complex brains of organisms as different as squid, fruit fly, and mouse. In addition, Pax2/5/8 is expressed in multimodal sensory systems in mollusks such as the esthetes and the ampullary system of polyplacophorans as well as the eyes of cephalopods. Pax2/5/8-expressing cells are present in regions where the future sensory cells such as the polyplacophoran esthetes are situated and hence Pax2/5/8 expression probably predates sensory cell development during ontogeny. In mollusks, Pax2/5/8 is only expressed in derivatives of the ectoderm and hence an ancestral role in molluscan ectoderm differentiation is inferred.
A coordinated tempo-spatial expression of transcription factors is required for cell type specification and the differentiation of the three germ layers into distinct organ systems during bilaterian ontogeny . Studies on ecdysozoan and deuterostome model organisms have elucidated the multifold roles of transcription factors during development; however, a major gap in knowledge does exist for the Lophotrochozoa, one of the three bilaterian superphyla . For inferences concerning the evolutionary conservation versus plasticity of transcription factor expression patterns and their putative roles, a broad comparative approach is imperative. Among the lophotrochozoans, Mollusca is the most speciose phylum and exhibits extremely diverse body plans and adaptative capacities ranging from worm-like spicule-bearing aplacophorans and sessile bivalved clams to highly motile squids (Fig. 1). No less diverse are the ontogenies, with direct and indirect developers and larval forms ranging from trochophore and pericalymma larvae to veligers . For the detection of settlement cues, predator avoidance, location of prey, as well as navigation in a complex environment, the larvae are equipped with sensory organs. The most prominent larval sensory organ is the apical organ that is present in representatives of all molluscan classes investigated so far, except for cephalopods [4–6]. It is equipped with flask-shaped (ampullary) cells and in some mollusks with additional peripheral cells. A so-called ampullary system has been described in polyplacophoran larvae and appears to constitute an apomorphy of the class . It comprises four pairs of sensory cells that are located dorsolaterally and ventrolaterally in the episphere of the trochophore. Cells are of the ampullary sensory cell type, exhibit serotonin-like and FMRFamide-like immunoreactivity and are innervated by the cerebral commissure . Another larval sensory organs, the so-called post-anal organ, only appears to exist in protobranch bivalves . It is composed of sensory cells with two types of cilia and a cavity in proximity to the anus. Coleoid cephalopods do not possess larvae as direct developers but a plethora of sensory organs such as complex camera eyes. These have been subject of a recent study suggesting that a vast number of genes expressed in the squid camera-type eye is not found in the transcriptomes of the nautiloid pinhole eye . Sensory organs as well as other epidermal epithelia or the nervous system are ectodermal derivatives . In mollusks, sensory organs have been described by histological and immunochemical techniques, however, the molecular underpinnings of their development and evolution remain largely obscure (but see ).
Members of the paired box (Pax) gene family encode transcription factors that are involved in a variety of processes such as the development of excretory organs, the differentiation of the musculature, neuroectoderm specification, and the development of auditory/geosensory systems [10, 11]. In bilaterians, nine orthologous groups of Pax proteins are known that can be distinguished on the basis of the presence or absence of highly conserved structural domains . Members of the Pax2/5/8 group exhibit a conserved N-terminal paired domain, an octapeptide, a partial homeodomain, and a C-terminal transactivation domain . The N-terminal domain and the homeodomain serve as DNA binding sites, the octapeptide is the site for activation or repression of the gene, and the C-terminal domain is the protein/protein interaction domain . In vertebrates and ecdysozoans, Pax2/5/8 is involved among others in the genesis of the auditory/geosensory system as well as in the formation of the midbrain/hindbrain boundary [13, 14]. For lophotrochozoans, so far only one study has been published on the expression domains and the putative function of Pax2/5/8 during development. In this study, Pax2/5/8 expression was found in the statocysts, i.e. the geosensory organ, of the gastropod veliger larva .
In order to elucidate putative roles of Pax2/5/8 during molluscan development, we investigated the expression of its orthologs in representatives of both major molluscan lineages, the Aculifera and the Conchifera [15, 16] (Fig. 1). As a polyplacophoran, Acanthochitona crinita has been shown to exhibit putative ancestral features of the Aculifera , while the conchiferan Nucula tumidula is a protobranch that holds a basal position among the bivalves and possesses a unique pericalymma-type larva [4, 18]. A. crinita and N. tumidula do not possess statocysts in their larvae and exhibit a less centralized CNS. In contrast, the conchiferan cephalopod Idiosepius notoides has a highly centralized CNS and statocysts. This study will test whether Pax2/5/8 is expressed in ectodermal, endodermal, or mesodermal derivatives during the ontogeny of all three molluscan species and it will shed light on the shared Pax2/5/8 expression among bilaterian representatives.
Pax2/5/8 gene orthologs and phylogenetic analysis
The multiple sequence alignment of Pax2/5/8 orthologs shows the highly conserved N-terminal paired domain, the octapeptide domain, and the lysine-arginine rich region (Fig. 2b). In particular the octapeptide is less conserved among the bilaterian Pax2/5/8 orthologs. The C-terminal partial homeodomain was not included in the alignment since it is even less well conserved. The phylogenetic analysis included proteins of all Pax families and demonstrates that all three molluscan Pax2/5/8 amino acid sequences cluster with their bilaterian orthologs (Fig. 3).
Pax2/5/8 expression in the polyplacophoran Acanthochitona crinita
The gastrula stage is reached between 3–5 h post fertilization (hpf) at 20 °C. During this stage Acr-Pax2/5/8 is expressed in ectodermal cells of the apical region (Fig. 4a–c). Early trochophore larvae hatch at approximately 12 hpf. They are lecithotrophic and are slightly ovoid in shape, with an apical organ and apical tuft in the episphere (Fig. 4d–h). The apical organ comprises three serotonin-like immunoreactive flask-shaped cells and faint FMRFamide-like immunoreactivity (Fig. 4d, e). Two other groups comprising four FMRFamide-like immunoreactive cells each are located in the episphere close to the prototroch (Fig. 4e). The episphere is separated from the hyposphere by two rows of trochoblasts that give rise to the prototroch. Acr-Pax2/5/8 is expressed in different domains, here referred to as “groups”, located mainly in the larval episphere (Fig. 4e–l). Several Acr-Pax2/5/8-expressing cells are located in the ectoderm of the dorsal episphere adjacent to both rows of trochoblasts (group 1 in Figs. 4e, j, k, l and 5). In addition, two other groups are located bilaterally and anterior to group 1 (group 2 in Figs. 4f, j, l and 5). Each group 2 comprises approximately four ectodermal cells that are located subepidermally and project dendrites through the epidermis (group 2 in Figs. 4f, j and 5). Two single Acr-Pax2/5/8-expressing cells are located bilaterally, immediately adjacent to the two rows of trochoblasts and more ventrally to group 2 in the episphere (group 3 in Figs. 4g, i and 5). In the center of the ventral episphere, a group of 5–10 ectodermal Acr-Pax2/5/8-expressing cells is present (group 4 in Figs. 4h, k, l and 5). The only Acr-Pax2/5/8 expression domain that is located in the hyposphere of the early trochophore is situated in the region of the nascent shell fields in the dorsal hyposphere (arrowheads in Figs. 4f, g, j, k and 5). The apical organ does not express Acr-Pax2/5/8 (red dashed circle in Fig. 4f, j, l). In further developed trochophore larvae (35 hpf), Acr-Pax2/5/8 expression is less pronounced in the episphere, while Acr-Pax2/5/8 expression in the hyposphere gains in intensity (Fig. 5; arrowheads in Fig. 6a, c). This and other Acr-Pax2/5/8-expressing cell groups are located in similar expression domains as described for earlier developmental stages (Fig. 5; cf. arrowheads in Fig. 4j, k with arrowheads in Fig. 6a, c). Acr-Pax2/5/8-expressing cells are present in the expression domains of groups 1, 2, and 4 of previous stages (c.f. Figs. 5 and 6a, c, d). Furthermore, serotonin-like immunoreactive cell somata are located in corresponding regions, such as the lateroventral ampullary cells corresponding to Acr-Pax2/5/8-expressing group 2, or the additional sensory cells that are located in the dorsal episphere that correspond to group 1 (Fig. 6b). In the expression domains of Acr-Pax2/5/8-expressing groups 3 and 4, FMRFamide-like immunoreactive cell somata are located (not shown). In further developed trochophore larvae (35 hpf), the apical organ does not exhibit Acr-Pax2/5/8 expression, however, very close to it, Acr-Pax2/5/8-expressing cells are located (Fig. 6c). In metamorphic competent trochophore larvae (65 hpf), the Acr-Pax2/5/8 expression pattern is similar to the one of earlier trochophores (35 hpf) (Figs. 5 and 6e–k). Strong Acr-Pax2/5/8 expression is present in cells of the nascent shell fields (Figs. 5 and 6e). These Acr-Pax2/5/8-expressing cells are located in the epidermis and project dendrites to the periphery (inset in Fig. 6j). In addition, Acr-Pax2/5/8-expressing cells are numerous in groups 1–4 during earlier development (Figs. 5 and 6f–j). Interestingly, all these Acr-Pax2/5/8 expression domains house FMRFamide-like-immunoreactive (usually flask-shaped) cells (Fig. 6j).
An additional Acr-Pax2/5/8 expression domain is present close to the lateral borders of the foot (Fig. 5, group 5 in Fig. 6h, i). Both groups 5 are interconnected via a thin band of Acr-Pax2/5/8-expressing cells that runs along the abapical side of the mouth (Figs. 5 and 6h, i). Few Acr-Pax2/5/8-expressing cell are present in the posterior mantle fold (Figs. 5 and 6k), while no Acr-Pax2/5/8-expressing cells are located in the region of the apical organ (dashed circle in Fig. 6k).
Settled animals (130 hpf) commence to reduce their episphere and the foot becomes more pronounced. Metamorphosis is largely completed with the shedding of the prototroch, although the eighth shell plate is not formed until probably several weeks after settlement. Early settled animals solely possess few faintly stained Acr-Pax2/5/8-expressing cells in the region of the shell fields (Fig. 6l) and only non-specific staining was observed in the perinotum that surrounds the plates (Fig. 6l).
Pax2/5/8 expression in the bivalve Nucula tumidula
Embryos of N. tumidula reach the 2–4 cell stage by 12 hpf at 6.5 °C and the zygotes gastrulate by approximately 24 hpf. During gastrulation, Ntu-Pax2/5/8 is expressed in approximately five ectodermal cells close to the blastopore (Fig. 7a). Early pericalymma larvae swim by approximately 48 hpf, are roundish, and possess 5–7 Ntu-Pax2/5/8–expressing cells on the future dorsal side below the test-cells (Fig. 7b–d). Slightly further developed pericalymmae are characterized by an apical organ with a ciliary tuft, while mouth and anus are located on the abapical side (Fig. 5). They are covered with large, multiciliary calymma (test) cells. Subsequently, three rows of test cells bear three bands of longer cilia, while shorter cilia still cover the area around the apical tuft as well as the abapical side around mouth and anus. These larvae express Ntu-Pax2/5/8 in the area of the dorsal mantle (Fig. 7e) and the expression domain also includes the shell gland that develops as an invagination after gastrulation on the future dorsal side (Fig. 7f, g). During subsequent development the apical-abapical axis elongates and the foot develops on the ventral side. Meanwhile, the shell field on the dorsal side has secreted two shell valves. Further developed larvae express Ntu-Pax2/5/8 in more ventral and apical mantle regions (Fig. 7h). Late pericalymma larvae express Ntu-Pax2/5/8 along the dorsal mantle region as well as in ventral domains of the mantle (Figs. 5 and 7i–k). Larvae that are metamorphic competent are more ovoid in shape than earlier larvae and express Ntu-Pax2/5/8 mostly in dorsal domains but also in ventral ones (summarized in Figs. 5 and 7l). Lecithotrophic larvae swim until 15–20 dpf in the water column and finally settle. During the settlement process that may last for several hours, larvae appear to probe the substrate and metamorphosis commences with shedding and ingestion of the test cells and the apical organ (not shown). Larvae express Ntu-Pax2/5/8 mainly bilaterally in two postero-ventral domains of the central mantle (Figs. 5 and 8a, b). Faint expression is still visible on the dorsal side close to the hinge of the developing valves (Figs. 5 and 8a, b). In individuals that have lost their test cells, Ntu-Pax2/5/8 is expressed in the posteroventral region of the mantle and along the mantle margin (Figs. 5 and 8c, d). Early juveniles express Ntu-Pax2/5/8 in lateral portions of the ventral mantle (Figs. 5 and 8e, f).
Pax2/5/8 expression in the cephalopod Idiosepius notoides
As a typical cephalopod, I. notoides is a direct developer with discoidal cleavage. Hatching of the juvenile occurs within 9–10 days after oviposition at 25 °C. In the present study, we adhere to the common designation of embryological orientation in cephalopods, i.e., the mantle apex is considered dorsal and the cephalic region including the brachial crown that surrounds the mouth is considered ventral. The funnel is located posteriorly, while the opposite side including the supraesophageal mass is considered anterior.
Anlagen of the CNS, the arms, the funnel, and the mantle are discernable from stage 18 onward (for a common staging system, see ). Early stage 19 individuals of the pygmy squid I. notoides express Ino-Pax2/5/8 in the mantle, the funnel, and all ten arms (Figs. 5 and 9a–c). Faint Ino-Pax2/5/8 expression is present in the developing gills of stage 19 individuals (Fig. 5; arrowheads in Fig. 9a). During subsequent development, Ino-Pax2/5/8 expression is located in the same domains as in stage 19 individuals (Fig. 9d–g). Ino-Pax2/5/8 is additionally expressed in the tissue that interconnects the arms (Fig. 9g) and in nervous tissue of the optic ganglia (arrowheads in Fig. 9e). Stage 24 individuals strongly express Ino-Pax2/5/8 in all arms, the funnel, eyes, and the epidermis of the cephalopedal region (Figs. 5 and 10a, b). Additionally, Ino-Pax2/5/8-expressing cells are present in the region of the interbasal lobes that both belong to the posterior basal lobes of the brain (Fig. 5; arrowheads in Fig. 10a). During stages 24–25, Ino-Pax2/5/8 is also expressed in the eyes (Figs. 5 and 10b, c). Stage 25 individuals still express Ino-Pax2/5/8 in the gills and in large areas of the epidermis of the cephalopedal region (Fig. 10c, d). Ino-Pax2/5/8 expression is only located in the proximal portions of all arms and the proximal portion of the funnel (Fig. 10c, d). Ino-Pax2/5/8 expression domains in the CNS include the anterior basal and the interbasal lobes that both belong to the supraesophageal mass (Fig. 10d). The vertical lobe complex as highest integration center does not express Ino-Pax2/5/8 (Fig. 10d). Ino-Pax2/5/8 is also faintly expressed in the posterior subesophageal mass (Fig. 10d). The posterior perikaryal layer of the middle subesophageal mass also expresses Ino-Pax2/5/8 (Fig. 10d). In late prehatching developmental stages such as stage 27, Ino-Pax2/5/8 is expressed in the gills (Fig. 10e), in the posterior perikaryal layer of the middle subesophageal mass, and around the buccal mass (Fig. 10f). Ino-Pax2/5/8 expression in the interbasal lobes persists in stage 27 individuals (Fig. 10g, h). Faint Ino-Pax2/5/8 expression is still present at the base of the funnel (Fig. 10h).
Evolution of Pax2/5/8 proteins
Our phylogenetic analysis confirms the identities of all three molluscan Pax2/5/8 proteins presented herein. The sequences of Acanthochitona crinita, Idiosepius notoides, and Nucula tumidula exhibit a paired domain, an octapeptide, the lysine-arginine rich region, as well as the partial homeodomain (not shown) diagnostic for Pax2/5/8.
Pax2/5/8 is expressed in the development of multimodal sensory systems of mollusks
Among bilaterians, Pax2/5/8 expression is particularly well-investigated in the Deuterostomia and Ecdysozoa [13, 20]. Studies on various deuterostomes revealed that orthologs of Pax2/5/8 play a role during development of auditory/ geosensory organs [21–23]. Interestingly, the phylogenetically distantly related gastropod mollusk Haliotis asinina also expresses Pax2/5/8 in its geosensory organ, the statocysts . This present study demonstrated that homeobox genes may be recruited into the development of organs with similar function but of different evolutionary origin (i.e., in non-homologous systems). It also poses the question as to how Pax2/5/8 is recruited during the ontogeny of organisms without auditory and equilibrial sense.
In our study, we included the chiton Acanthochitona crinita and the bivalve Nucula tumidula, both lophotrochozoan representatives without larval statocysts, in order to assess the role of Pax2/5/8 during their development. Pax2/5/8 is expressed in several groups of ectodermal cells in the episphere and the hyposphere of the polyplacophoran trochophore. Although A. crinita and N. tumidula possess larval apical (sensory) organs, no Pax2/5/8 expression was detected in components of this organ. In A. crinita cells of groups 2–3 are located in a similar location as cells of the larval ampullary sensory system that is an autapomorphy of Polyplacophora (Figs. 5 and 6 in present study; ). Cells of the ampullary system are flask-shaped and probably chemosensory rather than mechanosensory . Pax2/5/8 expression and FMRFamide-like immunoreactivity may be co-localized in ectodermal cells in the dorsal episphere adjacent to the prototroch (group 2, c.f. Figs. 5 and 6j). Further Pax2/5/8-expressing ectodermal cells (group 4) are located in the ventral episphere and might correspond to the four serotonin-like immunoreactive flask-shaped cells (c.f. group 4 in Figs. 4h, 5, 6b). The immunochemical and molecular phenotype, their flask-shaped somata, and their subepithelial location argue for a sensory role of the above-mentioned cells . This might also be true for the Pax2/5/8-expressing groups 5 that are located on the ventral side along the borders of the foot of advanced trochophore larvae (Fig. 5). Both groups are connected to each other via a small band along the abapical side of the mouth/ apical side of the foot (Fig. 5). This expression domain also includes a region where the “Schwabe organ”, a putative chemosensory organ of a polyplacophoran ingroup, is located . Although the chitonid A. crinita has not been studied with respect to the presence of this organ, it is likely that it is not present in A. crinita since other congeneric species of the Chitonida lack this sensory structure. Pax2/5/8-expressing cells are also located in the posterior mantle fold that exhibits the chemosensory osphradia of adult chitons (arrowheads in Fig. 6k; see also . In both described expression domains, however, no Pax2/5/8 expression is found in postmetamorphic chitons.
A prominent group of Pax2/5/8-expressing cells is present in the shell fields of A. crinita (arrowheads in Figs. 4j, k, 5 and 6a). The majority of these cells is located in central regions of each shell field and exhibits processes that penetrate the epidermis and the nascent calcareous layers of the shell plates (inset in Fig. 6j). These processes exhibit FMRFamide-like immunoreactivity in earlier developmental stages than the associated cell somata (c.f. Figs. 4j, 5 and 6a, j). Judging by their location and morphology, these Pax2/5/8-expressing and FMRFamide-like immunoreactive cells are probably esthetes. Esthetes are primarily photoreceptive organs that constitute a polyplacophoran autoapomorphy and penetrate the surface of the shell plates of adult specimens [25–27]. During ontogeny, esthetes may first be identified by the narrow holes of their canals in shell plates of metamorphosed specimens .
In the pygmy squid, Pax2/5/8 is expressed in various organs that are known to be equipped with sensory cells. During a short developmental time frame (stages 24–25), Ino-Pax2/5/8 is expressed in the eyes (Figs. 5 and 10b, c), resembling the situation of the adult gastropod Haliotis asinina . The brachial crown, i.e. all five arm pairs, expresses Pax2/5/8 in epidermal cells until stage 24 (Fig. 10a, b). In stage 25 individuals, Pax2/5/8 expression is restricted to the proximal part of each arm/tentacle. Expression fades in subsequent developmental stages.
Cephalopod arms and tentacles house a plethora of sensory cells. Although cell proliferation experiments during cephalopod arm growth are still lacking, it is tempting to speculate that Pax2/5/8 is expressed in ectodermal sensory cells in regions of high proliferation activity, such as the proximal portion of the brachial crown in stage 25 specimens. The funnel shows a similar pattern as the brachial crown. During early development the entire funnel expresses Pax2/5/8, while only the proximal portion expresses Pax2/5/8 in subsequent stages until it fades entirely in late prehatching stages. Here, Pax2/5/8 appears to be expressed in the epidermis as it has also been reported for vertebrates . Along the lateral sides of the foot the future gills develop and Pax2/5/8 expression predates their development.
Expression of Pax2/5/8 in the molluscan mantle
Pax2/5/8 is expressed in different regions of the developing mantle of the aculiferan polyplacophoran Acanthochitona crinita, the basal bivalve Nucula tumidula, and the cephalopod Idiosepius notoides. This corroborates data on the veliger larva of Haliotis asinina, where Pax2/5/8 is expressed in the anterior region of the mantle that is known to be richly equipped with sensory cells . In contrast to the esthetes as distinct expression domain of A. crinita, Pax2/5/8 is rather broadly expressed in the mantle of N. tumidula and I. notoides. During earlier development of N. tumidula, Pax2/5/8 expression is, however, first restricted to the dorsal mantle epithelium, including the shell gland. In contrast, the shell gland of H. asinina does not exhibit Pax2/5/8 expression (present study and ). Shortly before metamorphosis, Pax2/5/8 expression shifts to the ventral side of the mantle and disappears in the dorsal region after metamorphosis. This may be explained by the onset of sensory cell development on the ventral side of the mantle that is exposed to water currents in the postmetamorphic, settled animal. Other ectodermal expression domains comprise the cephalopedal epidermis of stage 24 individuals of I. notoides (Fig. 10a).
Pax2/5/8 plays a role in cephalopod brain regionalization
Studies on mouse and Drosophila embryos revealed that several gene orthologs play a role in brain regionalization. The anterior-most brain regions of both organisms, i.e. the murine forebrain/midbrain as well as the fruit fly’s protocerebrum/deuterocerebrum, express Otd/Otx (Fig. 11; ). The adjacent more posterior brain regions, that are patterned in a spatially collinear fashion by anterior Hox gene orthologs, are separated from the anterior-most brain regions by a Pax2/5/8 expression domain (“midbrain/hindbrain boundary”; Fig. 11; ). Interestingly, a similar molecular fingerprint is manifested in the developing brain of the cephalopod Idiosepius notoides (present study) and other cephalopods investigated. Here, Pax2/5/8 expression is restricted to the cerebral ganglia during early development. In subsequent developmental stages transcripts are restricted to the anterior basal and the interbasal lobes that are located slightly anterior to the esophagus, thus dividing the cephalopod brain into a supraesophageal and a subesophageal mass (Fig. 11). Anterior portions of the posterior subesophageal mass that are located close to the esophagus also exhibit faint Pax2/5/8 expression. Notably, brain lobes that are located more anteriorly, such as the median basal lobe or the vertical lobe complex, do not express Pax2/5/8 in I. notoides. In late-stage prehatching individuals, Pax2/5/8 is only expressed in the interbasal lobes as well as in the posterior perikaryal layer of the middle subesophageal mass.
Otx is expressed in the supraesophageal mass of the cuttlefish Sepia officinalis and the pygmy squid I. notoides (Fig. 11; ; Wollesen, unpublished data). Compared to Pax2/5/8 it is, however, expressed in more anterior regions of the supraesophageal mass; i.e., in components of the vertical lobe complex (i.e. the superior frontal and subpedunculate lobe) and the median basal lobe of the supraesophageal mass. In contrast to Pax2/5/8 expression, no Otx expression was found in the subesophageal mass (present study, , Wollesen, unpublished data). In cephalopods, Hox genes are expressed in the CNS with exception of the cerebral ganglia/supraesophageal mass . Although no collinear spatial expression is reported from cephalopods, the Otx- and Pax2/5/8 expression patterns adhere to the overall pattern reported for Drosophila and mouse . Data on other bilaterians suggest a Six3 expression domain in the anterior-most region of the cephalopod brain, i.e. the vertical lobe complex (white box in Fig. 11; ). Surprisingly, Pax2/5/8 is not expressed in the developing nervous system of the polyplacophorans, bivalves, or gastropods investigated so far (present study; ). In contrast to cephalopods, the fruit fly, or mouse none of the latter mollusks possess a highly centralized brain.
Pax2/5/8 expression is restricted to the ectoderm in mollusks
Among lophotrochozoans, Pax2/5/8 expression during development has only been studied in mollusks (present study; ; but see  on the role of Pax2/5/8 in regenerating tissue of the adult polychaete Platynereis dumerilii). In all mollusks investigated so far, Pax2/5/8 is exclusively expressed in ectodermal domains such as the epidermis, the shell gland, the mantle, sensory cells, as well as the nervous system (present study; ). Considering that Pax2/5/8 is consistently expressed in ectodermal epithelia in all metazoans investigated so far, with solely the Onychophora and the Chordata expressing Pax2/5/8 in mesodermal domains such as the nephridia, the most parsimonious explanation is an independent recruitment of Pax2/5/8 in mesodermal epithelia development in these clades (Fig. 12). In nematodes and cnidarians the Pax2/5/8 orthologs are involved in the formation of endodermal structures such as the hindgut or the endoderm at the aboral end of the pharynx (Fig. 12). Since the anterior-most and posterior-most portions of the digestive tract are located at the endodermal-ectodermal boundary, the Pax2/5/8 expression domains in the nematode hindgut and the aboral end of the pharynx of the cnidarian Nematostella vectensis might be an evolutionary reminiscent of ancestral ectodermal expression domains in the ancestors of these taxa (see Fig. 12).
The present study shows that Pax2/5/8 is expressed in multimodal sensory systems in mollusks such as the esthetes and the ampullary system of polyplacophorans, the eyes of squids, but not in the larval apical organ of A. crinita and N. tumidula. Pax2/5/8 expression probably predates sensory cell development during ontogenesis, since Pax2/5/8-expressing cells are present in regions where the future sensory cells are situated. Compared to other bilaterians, Pax2/5/8 is not expressed in the development of the less centralized nervous system of bivalves, gastropods, and polyplacophorans, and hence it most likely lost its role in brain development in these mollusks. Interestingly, Pax2/5/8 is expressed along the trajectory of the esophagus that divides the cephalopod brain into a supraesophageal and a subesophageal mass. Together with Otx and Hox genes Pax2/5/8 might have been recruited into brain regionalization, thus representing an extreme case of convergent evolution of gene function to the situation found in vertebrates (mouse) and insects (fruit fly). Since Pax2/5/8 is largely expressed in ectodermal domains throughout the Bilateria, its ancestral role was most likely in the differentiation of this outer most germ layer of bilaterians.
Collection and culture of animals
Adults of the polyplacophoran Acanthochitona crinita were collected in the intertidal zone close to the Station Biologique Roscoff in Roscoff, France, during the summers of 2013 and 2014. Adults were kept at 20 °C in running seawater. Prior to spawning, adults were isolated and maintained separately in glass dishes. Spawning was induced by exposing adults to sunlight for 3–4 h or thrice each for 15 min to alternating water temperatures, i.e. 25–30 °C and 10–15 °C. Released sperm and oocytes were separated and the latter rinsed several times in Millipore-filtered seawater (MFSW) and fertilized immediately. Early cleavage stages were kept at 20 °C in running seawater, while post-gastrulation stages were kept in MFSW with the addition of antibiotics against bacterial and fungal growth (50 mg streptomycin sulfate (Sigma-Aldrich) and 60 mg penicillin (Sigma-Aldrich) per liter MFSW). Water was changed every other day and settlement of metamorphic competent larvae was induced by the addition of substrate on which adult animals had been found.
Adults of the protobranch bivalve Nucula tumidula were collected from sediment that was sampled with a hyperbenthic sled at 180–220 m depth on muddy seafloor in Hauglandsosen (Bergen, Norway) during the winters 2012 and 2013. Adult individuals were kept at 6.5 °C in MFSW (UV-treated) at the marine living animal facilities of the Department of Biology, University of Bergen. N. tumidula is dioecious and spawning of gametes was induced by keeping adults in seawater of alternating temperatures. Accordingly, adults were exposed thrice to water of 10–15 and 2 °C for 10 min each. Males that released sperm were separated immediately to avoid polyspermy and the oocytes were rinsed several times in MFSW. Oocytes were fertilized and developmental stages were cultured at 6.5 °C in MFSW in glass bowls with water changes every other day.
Adults of the pygmy squid Idiosepius notoides were dip-netted in the sea grass beds of Moreton Bay, Queensland, Australia. Adults were kept in closed aquaria facilities at the School of Biological Sciences of the University of Queensland in Brisbane and fed with various crustaceans. Fertilization is internal and females attached egg clutches to sea grass or to the glass surface of the aquaria. Embryos were cultured and staged as described previously .
RNA extraction and fixation of animals
Several hundred individuals of early cleavage stages, larvae, metamorphic competent individuals to early juveniles of Acanthochitona crinita and Nucula tumidula were collected and stored in RNAlater (Lifetechnologies, Vienna, Austria) at −20 to −80 °C. RNA was extracted with a RNA extraction kit (Qiagen, Roermond, Netherlands) and stored at −80 °C. For I. notoides, egg jelly and chorion of embryos were removed and RNA from approximately 300 specimens including freshly laid zygotes (stage 1) to hatchlings (stage 30) was extracted using TriReagent according to the manufacturer’s instructions (Astral Scientific Pty. Ltd., Caringbah, Australia; see also ). Several hundred individuals of representative developmental stages of all three species were fixed for in situ hybridization experiments as previously described [36, 37]. Specimens were either stored in 75 % EtOH or in 100 % Methanol at −20 °C.
RNAseq and transcriptome assembly
Pooled total RNA of all above-mentioned developmental stages was used for the preparation of amplified short insert cDNA libraries (150–250 bp insert size) for Nucula tumidula and Acanthochitona crinita (Eurofins, Ebersberg, Germany). Both libraries (Kit version TruSeq SBS Kit v3) were sequenced together with another two bar-coded libraries in 1 channel of HiSeq 2000 with Illumina chemistry v3.0. Sequences were demultiplexed according to the 6 bp index code with 0 mismatch allowed. In both cases a PhiX library was added before sequencing to estimate the error rate of the sequences. Sequencing resulted in a total amount of 8.160 Mbp for N. tumidula and 7.147 Mbp for A. crinita. Paired-end reads with an average read length of 100 bp were obtained. These reads were subsequently filtered (rRNA removal) and adapter and low quality sequences were trimmed, normalized, and assembled de novo into contigs with the assembler Trinity . The transcriptomes of A. crinita and N. tumidula comprised 166,556 contigs and 224,633 contigs, respectively. The sequencing strategy for developmental stages of I. notoides was described previously . Briefly, RNA of developmental stages was sequenced by 454 technology by Eurofins (Ebersberg, Germany). After filtering, adapter and low quality read trimming, the reads were normalized and assembled de novo by Eurofins.
Alignment and phylogenetic analysis
Amino acid sequences of bilaterian Pax2/5/8 orthologs retrieved from NCBI were used in blastp searches against the assembled transcriptomes of A. crinita, N. tumidula, and I. notoides. The amino acid sequences with the highest sequence similarities of these blastp hits were aligned in ClustalX 2.1 (for accession numbers see Table 1) and edited manually in Geneious Pro 5.5.6 (Biomatters, Auckland, New Zealand, www.geneious.com). One alignment including only metazoan Pax2/5/8 orthologs was performed in order to highlight the Pax2/5/8-specific domains (Fig. 2b), while another alignment comprising amino acid sequences of various Pax genes was done to show that Acr-Pax2/5/8, Ntu-Pax2/5/8, and Ino-Pax2/5/8 cluster with their metazoan orthologs (Fig. 2a). The conserved motifs of the paired domains of the latter alignment were used to construct a maximum-likelihood consensus tree (100 bootstrap replicates) with the program Phylip v.3.695 .
Molecular isolation of Pax2/5/8 sequence orthologs
RNA pooled from different developmental stages of A. crinita, N. tumidula, and I. notoides, respectively, was used for first-strand cDNA synthesis by reverse transcription using the First strand cDNA Synthesis Kit for rt-PCR (Roche Diagnostics GmbH, Mannheim, Germany). Gene-specific primers were designed from identified Pax2/5/8 orthologs of A. crinita, N. tumidula, and I. notoides and transcripts were amplified via standard PCR. PCR products were size-fractioned by gel electrophoresis and gel bands of the expected length were excised and cleaned up using a QIAquick Gel Extraction Kit (Qiagen). Subsequently, cleaned-up products were cloned by insertion into pGEM-T Easy Vectors (Promega, Mannheim, Germany). Plasmid minipreps were grown overnight, cleaned-up, and sent for sequencing. All Pax2/5/8-gene orthologs were identified using the BLASTx algorithm screening the database of the National Center for Biotechnology Information (NCBI). All sequences and phylogenetic data that have been obtained in this study have been deposited in appropriate data bases.
Probe syntheses and whole-mount in situ hybridization
From the miniprepped plasmids the probe template was amplified via standard PCR using M13 forward and reverse primers, and in vitro transcription reactions were performed with these templates, digoxigenin-UTP (DIG RNA Labeling Kit, Roche Diagnostics GmbH), and SP6/ T7 polymerase (Roche Diagnostics GmbH) for the syntheses of antisense riboprobes, according to the manufacturer’s instructions. Whole-mount in situ hybridization experiments were carried out as described previously [36, 37]. Briefly, specimens were rehydrated into PBT (PBS + 0.1 % Tween-20). They were treated with Proteinase-K (50-60 μg/ml for Acanthochitona crinita, 10 μg/ml for Nucula tumidula, and 25 μg/ml for Idiosepius notoides) in PBT at 37 °C for 15 min and prehybridized in hybridization buffer for 4 h or overnight at 65 °C for A. crinita and I. notoides, and at 56 °C for N. tumidula. Hybridization with a probe concentration of 0.5 to 1 μg/ml was carried out overnight at the same temperatures as those for prehybridization. For I. notoides, a minimum of 20 individuals per stage was investigated and for A. crinita and N. tumidula approximately 50 individuals per stage. In addition, negative controls were carried out with sense probes for all genes and developmental stages. The majority of whole-mount preparations were cleared in a 3:1 solution of benzyl-benzoate and benzyl alcohol, mounted on objective slides, and analyzed. The results were documented with an Olympus BX53 microscope (Olympus, Hamburg, Germany) and polyplacophoran and bivalve developmental stages were additionally scanned with a Leica confocal SP5 II microscope (Leica Microsystems,Wetzlar, Germany) using bright-field and autofluorescence scans as well as the reflection mode . If necessary, images were processed with Adobe Photoshop 9.0.2 software (Adobe Systems, San Jose, USA) to adjust contrast and brightness. Sketch drawings were done with Adobe Illustrator CS5 software (Adobe Systems).
Specimens stored in 75 % EtOH or 100 % MetOH at −20 °C were rehydrated and washed several times in PBS. Samples were rinsed thrice for 20 min each in PBS with 2 % Triton X-100 (PBT) to increase tissue permeability and were then blocked for 4 h in PBT with 3 % normal swine serum (NSS; Jackson ImmunoResearch, West Grove, USA) at RT. Specimens were incubated for 12 to 24 h at RT in a cocktail of primary antibodies directed against serotonin (raised in rabbit, polyclonal, Immunostar, Hudson, USA) or FMRFamide-related peptides (raised in rabbit, polyclonal; Biotrend, Cologne, Germany), and a primary antibody against acetylated α-tubulin (raised in mouse, monoclonal; Sigma Aldrich), all diluted 1:800 in PBT with 2 % NGS. After rinsing specimens five times for 20 min each in PBS at RT, secondary fluorochrome–coupled antibodies were applied in a 1:1.000 dilution over night at RT. These cocktails contained PBT with 1 % NGS and Alexa Fluor 488 (anti-rabbit; Invitrogen) and Alexa Fluor 633 (anti-mouse, Invitrogen), as well as 1 % Hoechst 33342 ((2'-[4-ethoxyphenyl]-5-[4-methyl-1-piperazinyl]-2,5'-bi-1H-benzimidazole trihydrochloride trihydrate); Thermo Scientific) for cell nuclei staining. Samples were rinsed in PBS five times for 20 min each at RT and subsequently mounted on glass slides in Fluoromount G (Southern Biotech, Birmingham, Alabama, USA). Glass slides were stored in the dark at 4 °C until the mounting medium had hardened. Specimens were investigated and results analyzed as described above. Isosurfaces of FMRFamide-like reactive and serotonin-like immunoreactive elements were created by surface rendering algorithms using the 3D reconstruction imaging software IMARIS (Bitplane, Zurich, Switzerland). Specificity of the FMRFamide-like and serotonin-like antibodies was tested by omission of the primary antibodies and rendered no signal.
Statement of ethical approval
Animals were collected, anesthetized, and fixed according to internationally recognized standards (University of Queensland Animal Welfare Permit No. 158/09 “The cultivation of Idiosepius (pygmy squid) for studies in developmental biology” to BMD). Field work permission was granted by the ethics committee for every collection site of I. notoides used in this study.
Schep AN, Adryan B. A comparative analysis of transcription factor expression during Metazoan embryonic development. Plos One. 2013;8:6.
Gazave E, Guillou A, Balavoine G. History of a prolific family: the Hes/Hey-related genes of the annelid Platynereis. EvoDevo. 2014;5:29.
Haszprunar G, Wanninger A. Molluscs. Curr Biol. 2012;22:13.
Zardus JD, Morse MP. Embryogenesis, morphology and ultrastructure of the pericalymma larva of Acila castrensis (Bivalvia: Protobranchia: Nuculoida). Invert Biol. 1998;117:221–44.
Haszprunar G, Friedrich S, Wanninger A, Ruthensteiner B. Fine structure and immunocytochemistry of a new chemosensory system in the chiton larva (Mollusca: Polyplacophora). J Morphol. 2002;251:210–8.
Nielsen C. Animal Evolution. Interrelationships of the living phyla. 3rd ed. Oxford: Oxford University Press; 2012.
Yoshida MA, Ogura A, Ikeo K, Shigeno S, Moritaki T, Winters GC, et al. Molecular evidence for convergence and parallelism in evolution of complex brains of cephalopod molluscs: Insights from visual systems. Integr Comp Biol. 2015. doi:10.1093/icb/icv049.
Fioroni P. Allgemeine und vergleichende Embryologie der Tiere. 2nd ed. Berlin: Springer; 1992.
O’Brien EK, Degnan BM. Expression of Pax258 in the gastropod statocyst: insights into the antiquity of metazoan geosensory organs. Evol Dev. 2003;5:572–8.
Thompson JA, Ziman M. Pax genes during neural development and their potential role in neuroregeneration. Progr Neurobiol. 2011;95:334–51.
Blake JA, Ziman MR. Pax genes: regulators of lineage specification and progenitor cell maintenance. Development. 2014;141:737–51.
Balczarek KA, Lai ZC, Kumar S. Evolution and functional diversification of the paired box (Pax) DNA-binding domains. Mol Biol Evol. 1997;14:829–42.
Heller N, Brändli AW. Xenopus Pax-2/5/8 orthologues: Novel insights into Pax gene evolution and identification of Pax-8 as the earliest marker for otic and pronephric cell lineages. Dev Genet. 1999;24:208–19.
Lichtneckert R, Reichert H. Insights into the urbilaterian brain: conserved genetic patterning mechanisms in insect and vertebrate brain development. Heredity. 2005;94:465–77.
Smith SA, Wilson NG, Goetz FE, Feehery C, Andrade SCS, Rouse GW, et al. Resolving the evolutionary relationships of molluscs with phylogenomic tools. Nature. 2011;480:364–U114.
Kocot KM, Cannon JT, Todt C, Citarella MR, Kohn AB, Meyer A, et al. Phylogenomics reveals deep molluscan relationships. Nature. 2011;477:452–7.
Scherholz M, Redl E, Wollesen T, Todt C, Wanninger A. Aplacophoran mollusks evolved from ancestors with polyplacophoran-like features. Curr Biol. 2013;23:2130–4.
González VL, Andrade SCS, Bieler R, Collins TM, Dunn CW, Mikkelsen PM, et al. A phylogenetic backbone for Bivalvia: an RNA-seq approach. Proc Roy Soc B. 2015;282:42332.
Yamamoto M. Normal embryonic stages of the pygmy cuttlefish, Idiosepius pygmaeus paradoxus Ortmann. Zool Sci. 1988;5:989–98.
Dahl E, Koseki H, Balling R. Pax genes and organogenesis. Bioessays. 1997;19:755–65.
Wada H, Saiga H, Satoh N, Holland PWH. Tripartite organization of the ancestral chordate brain and the antiquity of placodes: insights from ascidian Pax-2/5/8, Hox and Otx genes. Development. 1998;125:1113–22.
Riley BB, Chiang MY, Farmer L, Heck R. The delta gene of zebrafish mediates lateral inhibition of hair cells in the inner ear and is regulated by pax2.1. Development. 1999;126:5669–78.
McCauley DW, Bronner-Fraser M. Conservation of Pax gene expression in ectodermal placodes of the lamprey. Gene. 2002;287:129–39.
Sigwart JD, Sumner-Rooney LH, Schwabe E, Heß M, Brennan GP, Schrödl M. A new sensory organ in “primitive” molluscs (Polyplacophora: Lepidopleurida), and its context in the nervous system of chitons. Front Zool. 2014;11:7.
Eernisse DJ, Reynolds PD. Polyplacophora. Microscopic Anatomy of the Invertebrates, Mollusca I, vol. 5. 1994. p. 55–110.
Reindl S, Salvenmoser W, Haszprunar G. Fine structural and immunocytochemical studies on the eyeless aesthetes of Leptochiton algesirensis, with comparison to Leptochiton cancellatus (Mollusca, Polyplacophora). J Submicrosc Cytol Pathol. 1997;29:135–51.
Speiser DI, Eernisse DJ, Johnsen S. A chiton uses aragonite lenses to form images. Curr Biol. 2011;21:665–70.
Kniprath E. Ontogenetic plate and plate field development in two chitons, Middendorffia and Ischnochiton. Wilhelm Roux’s Arch. 1980;189:97–106.
O’Brien EK, Degnan BM. Expression of POU, Sox and Pax genes in the brain of the tropical abalone, Haliotis asinina. Mar Biotech. 2000;2:545–57.
Yolder BK. Ciliary function in mammalian development. Oxford; Academic Press. 2008.
Buresi A, Baratte S, Da Silva C, Bonnaud L. Orthodenticle/otx ortholog expression in the anterior brain and eyes of Sepia officinalis (Mollusca, Cephalopoda). Gene Expr Patterns. 2012;12:109–16.
Lee PN, Callaerts P, de Couet HG, Martindale MQ. Cephalopod Hox genes and the origin of morphological novelties. Nature. 2003;424:1061–5.
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.
Backfisch B, Veedin Rajan VB, Fischer RM, Lohs C, Arboleda E, Tessmar-Raible K, et al. Stable transgenesis in the marine annelid Platynereis dumerilii sheds new light on photoreceptor evolution. Proc Natl Acad Sci U S A. 2013;110:193–8.
Wollesen T, Cummins SF, Degnan BM, Wanninger A. FMRFamide gene and peptide expression during central nervous system development of the cephalopod mollusk, Idiosepius notoides. Evol Dev. 2010;12:113–30.
Wollesen T, McDougall C, Degnan BM, Wanninger A. POU genes are expressed during the formation of individual ganglia of the cephalopod central nervous system. EvoDevo. 2014;5:41.
Fritsch M, Wollesen T, de Oliveira AL, Wanninger A. Unexpected co-linearity of Hox gene expression in an aculiferan mollusc. BMC Evol Biol. 2015;15:151.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptomes assembly from RNA-seq data without a reference genome. Nat Biotechnol. 2011;15:644–52.
Felsenstein J. PHYLIP - Phylogeny Inference Package (Version 3.2). Cladistics. 1989;5:164–6.
Jekely G, Arendt D. Cellular resolution expression profiling using confocal detection of NBT/BCIP precipitate by reflection microscopy. BIOTECH. 2007;42:751–5.
Vinther J, Sperling EA, Briggs DG, Peterson KJ. A molecular palaeobiological hypothesis for the origin of aplacophoran molluscs and their derivation from chiton-like ancestors. Proc R Soc B. 2011;279:1259–68.
Matus DQ, Pang K, Daly M, Martindale MQ. Expression of Pax gene family members in the anthozoan cnidarian, Nematostella vectensis. Evol Dev. 2007;9:25–38.
Pani AM, Mullarkey EE, Aronowicz J, Assimacopoulos S, Grove EA, Lowe CJ. Ancient deuterostome origins of vertebrate brain signalling centres. Nature. 2012;483:289.
Röttinger E, Saudemont A, Duboc V, Besnardeau L, McClay D, Lepage T. FGF signals guide migration of mesenchymal cells, control skeletal morphogenesis and regulate gastrulation during sea urchin development. Development. 2008;135:353–65.
Bassham S, Canestro C, Postlethwait PH. Evolution of developmental roles of Pax2/5/8 paralogs after independent duplication in urochordate and vertebrate lineages. BMC Biol. 2008;6:35.
Kozmik Z, Holland ND, Kalousova A, Paces J, Schubert M, Holland LZ. Characterization of an amphioxus paired box gene, AmphiPax2/5/8: developmental expression patterns in optic support cells, nephridium, thyroid-like structures and pharyngeal gill slits, but not in the midbrain-hindbrain boundary region. Development. 1999;126:1295–304.
Chamberlin HM, Palmer RE, Newman AP, Sternberg PW, Baillie DL, Thomas JH. The Pax gene egl-38 mediates developmental patterning in Caenorhabditis elegans. Development. 1997;124:3919–28.
Franke FA, Schumann I, Hering L, Mayer G. Phylogenetic analysis and expression patterns of Pax genes in the onychophoran Euperipatoides rowelli reveal a novel bilaterian Pax subfamily. Evol Dev. 2015;17:3–20.
Czerny T, Bouchard M, Kozmik Z, Busslinger M. The characterization of novel Pax genes of the sea urchin and Drosophila reveal an ancient evolutionary origin of the Pax2/5/8 subfamily. Mechan Dev. 1997;67:179–92.
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:7188.
Thomas Rattei und Thomas Eder are thanked for advice with transcriptome assembly. Dan Jackson (Göttingen), Martin Fritsch (Vienna), and Oliver Vöcking (Bergen) are thanked for technical advice. Emanuel Redl (Vienna) is thanked for collecting some of the postmetamorphic N. tumidula. TW thanks the staff of the research vessel Neomys (SBR in Roscoff) for collecting adult scaphopods and Henrik Glenner (Department of Biology, University of Bergen) for providing lab space and boat time. We are grateful to the crew of the RV Hans Brattström (University of Bergen) for assistance with the collection of adults of N. tumidula. TW was supported by ASSEMBLE (Association of European Marine Biological Laboratories) during his stay in at the SBR in Roscoff (France) and thanks the Faculty of Life Sciences, University of Vienna, for generous support. This work was support by an Australia Research Council grant to BMD and a grant by the Austrian Science Fund (FWF; grant number P24276-B22) on molluscan EvoDevo to AW.
The authors declare that they have no competing interests.
TW designed the project with AW. TW and SVR carried out the experiments and TW analyzed all data. TW drafted the manuscript and AW contributed to writing of the manuscript. SVR and TW performed the phylogenetic analysis. CT and BMD provided facilities in Bergen and Brisbane to rear Nucula tumidula and Idiosepius notoides, respectively, and helped to collect adults and rear developmental stages. All authors read and approved the final version of the manuscript.
About this article
Cite this article
Wollesen, T., Rodríguez Monje, S.V., Todt, C. et al. Ancestral role of Pax2/5/8 in molluscan brain and multimodal sensory system development. BMC Evol Biol 15, 231 (2015). https://doi.org/10.1186/s12862-015-0505-z
- Sensory cell