Ancestral and novel roles of Pax family genes in mollusks
© The Author(s). 2017
Received: 16 November 2016
Accepted: 18 February 2017
Published: 16 March 2017
Pax genes are transcription factors with significant roles in cell fate specification and tissue differentiation during animal ontogeny. Most information on their tempo-spatial mode of expression is available from well-studied model organisms where the Pax-subfamilies Pax2/5/8, Pax6, and Paxα/β are mainly involved in the development of the central nervous system (CNS), the eyes, and other sensory organs. In certain taxa, Pax2/5/8 seems to be additionally involved in the development of excretion organs. Data on expression patterns in lophotrochozoans, and in particular in mollusks, are very scarce for all the above-mentioned Pax-subfamilies, which hampers reconstruction of their putative ancestral roles in bilaterian animals. Thus, we studied the developmental expression of Pax2/5/8, Pax6, and the lophotrochozoan-specific Paxβ in the worm-shaped mollusk Wirenia argentea, a member of Aplacophora that together with Polyplacophora forms the Aculifera, the proposed sister taxon to all primarily single-shelled mollusks (Conchifera).
All investigated Pax genes are expressed in the developing cerebral ganglia and in the ventral nerve cords, but not in the lateral nerve cords of the tetraneural nervous system. Additionally, Pax2/5/8 is expressed in epidermal spicule-secreting or associated cells of the larval trunk and in the region of the developing protonephridia. We found no indication for an involvement of the investigated Pax genes in the development of larval or adult sensory organs of Wirenia argentea.
Pax2/5/8 seems to have a conserved role in the development of the CNS, whereas expression in the spicule-secreting tissues of aplacophorans and polyplacophorans suggests co-option in aculiferan skeletogenesis. The Pax6 expression pattern in Aculifera largely resembles the common bilaterian expression during CNS development. All data available on Paxβ expression argue for a common role in lophotrochozoan neurogenesis.
The morphological diversity of the phylum Mollusca is represented by eight recent clades, including the well-known gastropods, bivalves, and cephalopods. This bodyplan plasticity renders mollusks an ideal target for evolutionary and developmental studies. Although numerous aspects of intra-molluscan relationships are still controversially discussed, recent phylogenomic analyses show a basal dichotomy comprising the monophyletic Aculifera (including the eight-shelled Polyplacophora and the shell-less, spicule-bearing aplacophoran clades Neomeniomorpha or Solenogastres and Chaetodermomorpha or Caudofoveata) and the primarily single-shelled Conchifera (Monoplacophora, Gastropoda, Cephalopoda, Scaphopoda, Bivalvia) [1–3]. The vermiform Neomeniomorpha represents one of the least investigated molluscan taxa, although recent comparative studies have demonstrated its importance to reconstruct ancestral aculiferan traits [4, 5].
Adult neomeniomorphs lack eyes but do exhibit several sensory organs, including a vestibular (atrial) sense organ and a dorsoterminal sense organ that is often considered a homolog of the osphradium of other molluscan clades [6, 14–16] (Fig. 1b, c, d). Additionally, some neomeniomorphs exhibit a pedal commissural sac that was suggested to serve as a geosensory organ . Similar to many other spiralians, neomeniomorph larvae exhibit an apical organ with a ciliary tuft [13, 18]. Interestingly, and in contrast to the aplacophoran clades, representatives of the Polyplacophora exhibit ventrally positioned posttrochal larval eyes, which are lost some time after metamorphosis.
The highly conserved paired box (Pax) genes encode for a family of metazoan transcription factors, which are crucial for cell fate specification and tissue differentiation and are known to be involved in the development of excretory organs, myogenesis, neurogenesis, biomineralisation processes (skeletogenesis), and the development of visual and geosensory systems in numerous bilaterians (e.g., [19–21]). Although neomeniomorph larvae lack a number of sensory organs such as eyes and statocysts, they possess Pax genes whose orthologs have a conserved function in neurogenesis and sensory organ development. This poses the question as to where these genes are expressed during the ontogeny of neomeniomorphs.
Pax2/5/8 is known to be involved in the development of excretory systems of certain annelids, onychophorans, vertebrates, as well as in the formation of auditory/geosensory systems and in establishing the midbrain/hindbrain boundary of vertebrates as well as the deutocerebrum/tritocerebrum boundary of ecdysozoans [22–24]. The developmental expression and putative function of Pax2/5/8 has also been studied in a few lophotrochozoan taxa [25–28]. In mollusks, Pax2/5/8 is expressed during development of multimodal sensory systems of gastropods , polyplacophorans, and cephalopods . Furthermore, Pax2/5/8 is expressed in the mantle of gastropods, polyplacophorans, bivalves, and cephalopods, which is known to be rich in sensory structures [25, 27]. In the polychaete annelid Platynereis dumerilii Pax2/5/8 expression was found during development of photoreceptor cells of regenerating segments , while in the leech Helobdella austinensis Pax2/5/8 expression is confined to the neuroectoderm of the developing ventral nerve cords and to the developing nephridia . A recent study of Pax6 and Pax2/5/8 in two brachiopods suggested that these genes have a putative role as regulators of the segment polarity gene engrailed, although further details remain vague . Pax6 is generally involved in the development of the bilaterian central nervous system (CNS) and is a key player in eye gene regulatory networks of most bilaterians, including cephalopod and polyplacophoran mollusks [30–32]. Paxβ represents the lophotrochozoan-specific ortholog of the recently discovered Paxα/β subfamily and has hitherto been exclusively studied in the leech Helobdella austinensis [33–35]. Late embryos of H. austinensis show broad, segmentally restricted mesodermal expression as well as expression in the CNS and eyes during organogenesis . So far, the expression pattern of Paxα (the deuterostome and ecdysozoan Paxβ ortholog) is exclusively known from the onychophoran Euperipatoides rowelli .
In order to further assess putative ancestral versus novel roles of these important developmental regulators, we here provide the first detailed study of selected Pax family genes in a hitherto largely neglected but evolutionarily highly important molluscan clade, the Neomeniomorpha.
Adult and developmental stages of the neomeniomorph Wirenia argentea Odhner, 1921 were collected, maintained, and reared from January to May 2012, November 2012 to February 2013, and from November to December 2013, respectively, following Redl et al.  with the following modifications during the last season: The sediment sample was kept in 20μm-filtered and UV-sterilized sea water with a salinity of 35‰ (FSSW) precooled to 4°C. Every four days the adult specimens were transferred to clean plastic jars with fresh FSSW, which increased egg laying productivity. The freshly laid eggs were transferred into clean plastic jars with FSSW and kept under the same conditions as the adults. After hatching Wirenia develops via the so-called pericalymma or test cell larva. Herein, age of the larvae is given in days post hatching (dph). We fixed four morphologically distinguishable stages: freshly hatched test cell larva (0–1 dph), early test cell larva (6–7 dph), mid-stage test cell larva (10–11 dph), and late test cell larva (14–15 dph), whereby each stage encompasses a developmental time range of approximately 24 h.
RNA extraction and fixation of animals
All molecular biological procedures, ranging from RNA extraction until the end of the in situ hybridization protocol, were conducted using RNase-free (or diethylpyrocarbonate (DEPC)-treated) water. Total RNA was extracted from larvae using the Qiagen RNeasy Mini Kit (Qiagen, Venlo, The Netherlands) with the QIAshredder homogenizer (Qiagen). Prior to RNA extraction the larval material was either shock frozen on dry ice or conserved in RNAlater. Additional RNA was extracted from adult specimens. In order to separate animal tissue from the gut content and to avoid contamination, all adult specimens were starved for at least two weeks prior to RNA extraction. Adult RNA was extracted using TRI reagent (Sigma-Aldrich) according to the manufacturer’s instructions with the optional centrifugation step after homogenization and the following modifications: The animals were shock frozen with dry ice immediately before homogenization and RNA precipitation was performed with a 1:1 mixture of isopropanol and a high salt precipitation solution containing 0.8mol/l trisodium citrate dihydrate and 1.2mol/l sodium chloride. RNA pellets were dissolved in water and stored at −80°C.
Advanced larval stages were relaxed prior to fixation for 20 to 30 min at 4°C by adding a 3.2% magnesium chloride solution. Developmental stages were fixed for in situ hybridization with 4% paraformaldehyde (PFA) in 0.1M MOPS buffer (with 0.5M/l NaCl, 2mM/l MgSO4, and 1mM/l EGTA added) for 45 min at room temperature (RT). Fixed larvae were stepped into precooled (4°C) or prechilled (−20°C) 75% EtOH, washed three times for a total period of 15–30 min in precooled (4°C) or for 45 min to 1.5 h in prechilled (−20°C) 75% EtOH, and stored in fresh 75% EtOH at −20 °C.
RNAseq and transcriptome assembly
Larval and adult RNA samples were pooled and sequenced by Illumina technology (Eurofins, Ebersberg, Germany). Sequencing and bioinformatic processing of the resulting paired-end libraries as well as the transcriptome assembly were performed as described in Redl et al. .
Gene identification and orthology assignment
We identified candidate genes in our transcriptome database by reciprocal BLAST  searches using bilaterian Pax gene sequences from NCBI GenBank as queries. Nucleotide sequences of the best-fitting contigs were translated into amino acid sequences using Geneious, Version 6.1.6 (Biomatters, Auckland, New Zealand).
GenBank accession numbers of genes used in the phylogenetic analysis
GenBank accession number
Gene cloning and probe synthesis
Primers used for PCR
Fragment length (in bases)
Whole-mount in situ hybridization
Samples stored in 75% EtOH were stepped into 4% PFA in 1x Roti-Stock phosphate buffered saline (pH = 7.4; PBS; Carl Roth) with 0.05M/l EGTA (PPE) and decalcified in PPE for 1h at RT. Subsequently, samples were washed six times for 5 min each in PBS with 0.1% Tween 20 (PBT; Carl Roth) at RT and then heated to 37°C in a water bath during the last washing step. Proteinase K treatment was done with a solution of 10μg/ml Proteinase K (Roche Life Science) in PBT for 10 min at 37°C without agitation. The specimens were then washed twice for 5 min each at RT in PBT, twice for 5 min each in 1% triethanolamine (TEA) in PBT, four times for 5 min each in an instantly made mixture of 0.3% acetic anhydride and 1% TEA in PBT, and again twice for 5 min each in PBT. Afterwards, the samples were postfixed in 4% PFA in PBS for 45 min at RT and washed five times for 5 min each at RT in PBT. Samples were subsequently stepped into the hybridization buffer (HB) consisting of 50% formamide with 0.075M/l trisodium citrate, 0.75M/l sodium chloride, 5mM/l EDTA, 50μg/ml heparin sodium salt (Sigma-Aldrich), 1x Denhardt’s Solution (Carl Roth), 100μg/ml RNA from torula yeast, Type VI (Sigma-Aldrich), and 5% dextran sulfate sodium salt from Leuconostoc spp. (Sigma-Aldrich). The samples were then transferred into fresh HB, heated in a water bath to 56°C (in case of the Pax6 probe) or 60°C (in case of Pax2/5/8 and Paxβ), respectively, and prehybridized at these specific temperatures for 15–20 h. All RNA probes were diluted in HB to a final concentration of 1–2μg/ml, denatured for 10 min at 85°C, and applied to the samples. Hybridization was conducted for 24–26 h at the specific temperatures of 56°C and 60°C, respectively. Then, the samples were kept at hybridization temperature and washed three times for 20 min each in pre-warmed 50% formamide with 0.06M/l trisodium citrate, 0.6M/l sodium chloride, and 0.1% Tween 20, twice for 20 min each in 50% formamide with 0.03M/l trisodium citrate, 0.3M/l sodium chloride, and 0.1% Tween 20, and three times for 15 min each in 50% formamide with 0.015M/l trisodium citrate, 0.15M/l sodium chloride, and 0.1% Tween 20. The samples were then placed at RT to cool down. Afterwards, they were washed thrice for 20 min each at RT in 0.015M trisodium citrate solution with 0.15M/l sodium chloride and 0.1% Tween 20. Subsequently, the samples were washed thrice for 5 min each at RT in 0.1M maleic acid buffer (MAB) with 0.15M/l sodium chloride and 0.1% Tween 20 (pH = 7.5). Blocking of unspecific binding sites was done for 3 h at RT with a 2% solution of Blocking Reagent (Roche Life Science) in MAB. Anti-Digoxigenin (DIG)-AP, Fab fragments (Roche Life Science, Ref. 11093274910) were applied in a 1:2,500 or 1:5,000 dilution in 2% block solution for 13–16 h at 4°C. After incubation, the samples were rinsed eight times for 20 min each at RT in PBT, twice for 5 min each without agitation in 0.1M Tris buffer with 0.1M/l sodium chloride (pH = 9.5; AP buffer) and 0.1% Tween 20, and twice for 10 min each without agitation in AP buffer with 50mM/l magnesium chloride and 0.1% Tween 20. Finally, all samples were transferred into staining buffer (AP buffer with 50mM/l magnesium chloride, 7.5% polyvinyl alcohol, and 20μl/ml NBT/BCIP stock solution (Roche Life Science)). Color reaction took place at 4°C for 3–8 h (depending on probe, probe concentration, and DIG concentration) and was stopped by washing twice in 0.1M glycine in PBT (pH = 2.2) followed by additional two washes in PBT, 5 min each. Stained specimens were then fixed for 12–24 h at 4°C in 4% PFA in PBS, subsequently rinsed four times for at least 5 min each in PBT, and finally stored in PBT at 4°C.
Mounting and clearing
Larval stages were stepped into deionized water and washed four times for 5 min each in deionized water. Larvae were then stepped into EtOH and washed thrice for 5 min each in 100% EtOH. Next, the larvae were transferred into a 1:1 mixture of benzyl benzoate and benzyl alcohol, and subsequently mounted in this medium on microscope slides. Approximately 250 specimens were processed and investigated in total and 81 (25 with Pax2/5/8 expression, 31 with Pax6 expression, and 25 with Paxβ expression) were scanned with a confocal microscope.
Microscopy, 3D rendering, and image processing
Specimens were analyzed and light micrographs were taken on an Olympus BX53 microscope equipped with an Olympus DP73 camera and the software cellSens Standard, Version 1.11 (Olympus Corporation, Shinjuku, Tokyo, Japan). Confocal laser scanning microscopy was conducted using a Leica DMI6000 CFS microscope equipped with a Leica TCS SP5 II scanning system (Leica Microsystems, Wetzlar, Germany) and software LAS AF, Version 2.6.0 or 2.6.3. The autofluorescent signal was scanned in fluorescence mode using a 405 nm laser and gene expression signal was scanned in reflection mode using a 633nm laser (see [44–47]). The obtained confocal image stacks were processed and used to prepare 3D reconstructions of the gene expression signal as well as 3D renderings of larval tissue with Imaris x64, Version 7.3.1 (Bitplane AG, Zurich, Switzerland). The same software was used to create video files of larval stages including gene expression patterns. Generated images were finally processed with Adobe Photoshop CS6 Extended, Version 13.0.1 x64, and Adobe Photoshop CC 2015 (Adobe Systems). Figures and schematic drawings were generated with Adobe Illustrator CS5, Version 15.0.0 and Adobe Illustrator CC 2015, Version 1.0 (Adobe Systems).
Larval morphology of the neomeniomorph Wirenia argentea
Additional file 1: Video of an autofluorescence scan and subsequent 3D rendering of an early test cell larva (6–7 days post hatching) of Wirenia argentea (Neomeniomorpha). (AVI 17546 kb)
Larval development of Wirenia is characterized by the outgrowth of a posterior trunk from the base of the pseudo-blastopore, resulting in a mushroom-like appearance of later stages. The first external indication of this outgrowing trunk is visible in the early test cell larva, in which the major part of the developing trunk is still masked by the outer apical cap (Fig. 3). This mode of formation results in an epidermal fold enclosing the “peri-imaginal space” (“Peri-Imaginalraum” sensu Salvini-Plawen ), which is lined by the epidermal cell layer of the outgrowing trunk on both sides (Fig. 3d, e). The early test cell larva already exhibits a ventral mouth opening and the developing foregut but lacks the two posteriormost test cells of the freshly hatched test cell larva. The typical mushroom-like appearance of later stages is recognizable in the mid-stage test cell larva, where the trunk is elongated. Late test cell larvae are already in the settlement phase and the apical cap and the covering test cells start to degenerate. The length of this phase varies individually and can last for several days . At this stage, the posterior trunk is already thickened, elongated, and dorsally and laterally covered by spicules. The ventral trunk exhibits a ventrolongitudinal ciliary band, which marks the developing creeping sole.
Gene orthologs and phylogenetic analysis
Additional file 2: Pax2/5/8 expression in a freshly hatched test cell larva. (AVI 25082 kb)
Additional file 3: Pax2/5/8 expression in an early test cell larva. (AVI 29128 kb)
Mid-stage test cell larvae express War-Pax2/5/8 (Additional file 4) in cells of the anlagen of the cerebral ganglia. This cerebral expression can be subdivided in two pretrochal groups of War-Pax2/5/8-expressing cells. The most anterior group lies close to the cerebral commissure and adjacent and anterior to the epithelial cells of the developing foregut (#1 in Fig. 6k, l, n, o). The second group of pretrochal War-Pax2/5/8-expressing cells is located more posteriorly and lies adjacent to the trochoblasts (#2 in Fig. 6k-o). A third bilateral group of War-Pax2/5/8-expressing cells is located posttrochally and ventrally but slightly anterior to the invagination of the peri-imaginal space and the outgrowing trunk (#4 in Fig. 6k-o). Two further bilateral expression domains are located on both sides of the ventral mouth opening, which lies at the base of the invagination of the peri-imaginal space (#5 in Fig. 6o). Scattered War-Pax2/5/8 expression is also detected in cells on both sides of and adjacent to the developing creeping sole, i.e., the region where the paired ventral nerve cords develop (#6 in Fig. 6l, n, o). As in early test cell larvae, mid-stage test cell larvae still express War-Pax2/5/8 in the epidermal cell layer of the outgrowing trunk, which also lines the peri-imaginal space (Fig. 6k-o). A ventral gap of this epidermal expression is still present in cells of the developing creeping sole. As in early test cell larvae, a bilateral group of War-Pax2/5/8-expressing cells is located ventrally at the posterior pole of the trunk (#3 in Fig. 6l, n, o).
Additional file 4: Pax2/5/8 expression in a mid-stage test cell larva. (AVI 25505 kb)
The pretrochal expression of War-Pax2/5/8 in late test cell larvae (Additional file 5) is restricted to a marginal expression in the cerebral ganglia, with a more distinct expression close to the cerebral commissure (#1 in Fig. 6p, q, s). Similar to the expression pattern in mid-stage test cell larvae another bilateral expression domain of War-Pax2/5/8 is located posttrochally and rather ventrally of the mid-horizontal section plane, and anterior to the invagination of the peri-imaginal space and the outgrowing trunk (#4 in Fig. 6p-s). As in mid-stage test cell larvae, War-Pax2/5/8 is expressed in cells on both sides of and adjacent to the developing creeping sole, i.e., in the region where the paired ventral nerve cords are located (#6 in Fig. 6q, s). Furthermore, War-Pax2/5/8 is still expressed in the epidermal cell layer of the outgrowing trunk (Fig. 6p-t). Thus, the ectodermal expression of War-Pax2/5/8 in the epidermal layer of the outgrowing trunk is present throughout all developmental stages investigated (Fig. 6f, k, p). The same applies to the characteristic ventral gap of expression in the region of the developing creeping sole.
Additional file 5: Pax2/5/8 expression in a late test cell larva. (AVI 29202 kb)
Additional file 6: Pax6 expression in a freshly hatched test cell larva. (AVI 16204 kb)
Early test cell larvae exhibit strong expression of War-Pax6 (Additional file 7) in cells associated with the developing cerebral ganglia (#1 in Fig. 8e-h). Furthermore, strong ectodermal expression is present in two ventral stripes of epidermal cells, which line the peri-imaginal space on both sides and therefore show a hook-like appearance (#2 in Fig. 8e, f, g i, h). The outgrowing trunk exhibits three dorsal bilateral expression domains, each most likely composed of a pair of single War-Pax6-expressing epidermal cell (#3, #4 and #5 in Fig. 8f, g ii, h). A further paired bilateral expression domain (#6) is formed by a single cell belonging to the outer dorsolateral epidermal lining of the peri-imaginal space (Fig. 8f, g ii, h).
Additional file 7: Pax6 expression in an early test cell larva. (AVI 20713 kb)
No dorsal expression of War-Pax6 was observed in mid-stage test cell larvae. As in early test cell larvae, War-Pax6 is strongly expressed in the anlagen of the cerebral ganglia (#1 in Fig. 8i-l; Additional file 8). Likewise, the hook-like expression domains of War-Pax6 in the epidermal lining of the peri-imaginal space are constantly present in early and mid-stage test cell larvae (#2´ in Fig. 8i-k). Two longitudinal expression domains of War-Pax6 (#2´´) are separated by a small gap and are situated in direct posterior extension of the hook-like domains of #2´ (Fig. 8i-l). War-Pax6-expressing cells of #2´´ are epidermal cells flanking the epidermal cells of the developing creeping sole (Fig. 8i, k).
Additional file 8: Pax6 expression in a mid-stage test cell larva. (AVI 20261 kb)
The expression of War-Pax6 is generally weaker in late test cell larvae, but is still present in the cerebral ganglia (#1 in Fig. 8m-p; Additional file 9). The epidermal layer with spicule-secreting cells is now restricted to the elongated trunk, which causes the longitudinal appearance of the formerly hook-like expression domain of #2´ in this late stage (Fig. 8m, n, p). War-Pax6-expressing cells of #2´´ are also present in late test cell larvae and, together with the cells of #2´, form two longitudinal stripes of weak expression on both sides of the developing creeping sole, i.e., in the area where the ventral nerve cords develop (Fig. 8m, n). Separated from #2´, additional War-Pax6-expressing cells flank the ventral mouth opening at the base of the invagination of the peri-imaginal space (#7 in Fig. 8m, n, o).
Additional file 9: Pax6 expression in a late test cell larva. (AVI 19856 kb)
Additional file 10: Paxβ expression in a freshly hatched test cell larva. (AVI 30190 kb)
As is the case for War-Pax6, early test cell larvae express War-Paxβ (Additional file 11) in cells of the developing cerebral ganglia (#1) as well as in two ventral stripes of epidermal cells, which line the peri-imaginal space on both sides and therefore show a hook-like appearance (#2 in Fig. 9e-h). Furthermore, War-Paxβ is expressed in two rather dorsally located spots at the posterior pole of the outgrowing trunk (#3 in Fig. 9f). Each spot is most likely formed by a single War-Paxβ-expressing cell.
Additional file 11: Paxβ expression in an early test cell larva. (AVI 18566 kb)
In mid-stage test cell larvae War-Paxβ-expressing cell groups are located in similar expression domains as described for the early test cell larvae and thus closely resemble the War-Pax6 expression pattern of mid-stage test cell larvae (Additional file 12). Strong War-Paxβ expression is present in the anlagen of the cerebral ganglia (#1) and in two ventral stripes of ectodermal cells adjacent to the cells of the developing creeping sole (#2 in Fig. 9i-l). As is the case for the War-Pax6 expression, the stripes of War-Paxβ-expressing ectodermal cells of group 2 also line the fold of the peri-imaginal space, which causes a hook-like appearance at the anterior pole of the expression domain.
Additional file 12: Paxβ expression in a mid-stage test cell larva. (AVI 16661 kb)
Expression of War-Paxβ in late test cell larvae (Additional file 1) is comparatively low. The most intense expression is still present in the developing cerebral ganglia (#1 in Fig. 9m, n, p), which is congruent with the expression of War-Pax6 in late test cell larvae. A further analogy to the expression pattern of War-Pax6 is the location of two spots of War-Paxβ-expressing cells flanking the mouth opening at the base of the invagination of the peri-imaginal space (#4 in Fig. 9m, n, o). The ventral stripe-like War-Paxβ expression on both sides adjacent to the creeping sole (#2; similar to the younger stages) is hardly detectable in late test cell larvae and, if visible, is mostly limited to a faint expression in the ventral part of the trunk (#2 in Fig. 9m, n, p).
Additional file 13: Paxβ expression in a late test cell larva. (AVI 11934 kb)
Putative ancestral and co-opted roles of Pax2/5/8 in neomeniomorph mollusks
Pax2/5/8 is an ortholog of the Drosophila Sv/Spa and the vertebrate Pax2, Pax5, and Pax8 genes. These play a crucial role in the development and regionalization of the highly centralized and complex brains of vertebrates and insects as well as cephalopods, where they were most likely independently recruited into similar functions [24, 27]. Based on its expression in sensory structures of chordates (e.g., [23, 51]), insects , and mollusks , Pax2/5/8 was proposed to have a conserved role in the formation of structures responsible for balance and geotaxis in eumetazoans .
For Mollusca, our data on an aplacophoran supplement existing Pax2/5/8 expression data on polyplacophorans, gastropods, bivalves, and cephalopods [25, 27, 53]. Although Pax2/5/8 was found to be expressed in the adult brain of the gastropod Haliotis , there is no evidence that this gene is involved in its development . However, developmental expression of Pax2/5/8 in Haliotis was found in the statocysts, eyes, foot, and in putative chemosensory organs of the pallial chamber . The expression of Pax2/5/8 in putative (precursors of) sensory cells was also confirmed for a polyplacophoran, where it is probably expressed in the developing esthetes (shell eyes) and the larval ampullary system (cf. ). Though, it was proposed that Pax2/5/8 expression probably predates sensory cell development in multimodal sensory systems during molluscan ontogeny . Although neomeniomorphs do possess distinct sensory organs, such as the atrial sense organ or the dorsoterminal sense organ, we did not find any evidence for Pax2/5/8 expression in cells that could be unambiguously assigned to a distinct sensory structure. However, Wirenia exhibits Pax2/5/8 expression in cells dedicated to the developing cerebral ganglia and the ventral nerve cords, which is in stark contrast to the other investigated molluscan representatives with a considerable simple (i.e., little concentrated) CNS, such as bivalves, gastropods, and polyplacophorans, but is in accordance with Pax2/5/8 expression in the highly centralized brain of cephalopods. Since data for only one species per class-level taxon are currently available, more comparative analyses are required to further interpret this finding. Nevertheless, and in contrast to War-Pax6 and War-Paxβ, War-Pax2/5/8 is not equally expressed in the region of the developing cerebral ganglia, but instead shows increased expression in cells adjacent to the apical organ and the cerebral commissure (#1 in Fig. 6 f-t), as well as in a slightly more posterior location where the future pedal ganglia develop (#2 in Fig. 6 f-o). Although neomeniomorphs do not exhibit a highly centralized CNS, this specific Pax2/5/8 expression could hint towards a similar function in CNS regionalization comparable to its proposed function in cephalopods (cf. ).
Aside from its putative function in patterning neural or sensory structures, Pax2/5/8 is also expressed during development of excretory organs, such as the nephridia of onychophorans and annelids [26, 35] and the kidneys of vertebrates [55, 56]. Interestingly, Wirenia is the only investigated mollusk, which shows Pax2/5/8 expression in the region where the paired protonephridia of the larva develop (#4 in Fig. 6 k-t, see also ). Whether this hints towards a conserved function of Pax2/5/8 in bilaterian nephrogenesis or whether this represents a coopted function of the aplacophoran lineage remains unclear until more data on various clades become available.
Remarkably, Pax2/5/8 is expressed in epidermal spicule-secreting or associated cells of neomeniomorph as well as polyplacophoran mollusks . Although it was not explicitly mentioned for Polyplacophora, the respective light micrographs show a distinct expression signal in the spicule-secreting perinotum (or girdle) of Acanthochitona crinita (cf. Fig. 6 g in ). Not much is known about spicule secretion and its molecular background in aculiferan mollusks, but a recent study on the evolution of molluscan secretomes and biomineralization processes proposed that independent co-option of gene families are important driving forces acting on molluscan biomineralization . Since the secretion of calcareous spicules most likely constitutes an aculiferan autapomorphy and Pax2/5/8 appears not to be involved in molluscan shell development [25, 27], it seems likely that Pax2/5/8 expression in the spicule-secreting tissues of Polyplacophora and Neomeniomorpha evolved along the line leading towards the Aculifera after their split from their conchiferan allies.
Conserved Pax6 expression in aculiferan molluscs
Pax6 is generally expressed in developing nervous systems and in various photosensory organs, regardless of their overall morphology in numerous metazoans, which underlines its important role as regulatory gene in eye developmental networks (see, e.g., [58–66]). Prior to our study, molluscan Pax6 expression had only been studied in developmental stages of cephalopods [30, 31, 67, 68] and a polyplacophoran , as well as in an adult gastropod . The latter study showed that neural expression of Pax6 in adults of the gastropod Haliotis asinina is in the cerebral and pleuropedal ganglia as well as in sensory structures such as the eyes, tentacles, and gills . Cephalopod developmental Pax6 expression in the squids Loligo opalescens and Euprymna scolopes is confined to their brain lobes, eyes, and olfactory organ [30, 67]. The polyplacophoran Leptochiton asellus likewise shows Pax6 expression in the developing nervous system as well as in the ventrally and posttrochally positioned larval eyes . Since the latter study largely focused on eye development, we analyzed Pax6 expression during neurogenesis of another polyplacophoran representative, Acanthochitona crinita (unpublished data). Our data confirm Pax6 expression in the region of the developing ventral nerve cords but not in the lateral nerve cords of the tetraneural nervous system of this polyplacophoran species. Interestingly, this finding is congruent with Pax6 expression in later developmental stages of the neomeniomorph Wirenia argentea, which argues for a conserved expression of this gene during neurogenesis of the ventral nerve cords of protostomes (cf., e.g., data on various crustaceans, insects, onychophorans, planarians, annelids, and polyplacophorans) [32, 62, 63, 66, 69–71]. Moreover, a conserved bilaterian mode of expression of Pax6 is also present in the larval episphere of W. argentea, where it correlates with the development of the cerebral ganglia.
Freshly hatched test cell larvae of Wirenia express Pax6 in two patches at the ventral base of the pseudo-blastopore, i.e., in the region where a pair of neurite bundles emerges from the posterior neurogenic domain (cf. ). Cells of the apical organ, the only larval sensory organ known from neomeniomorphs, do not exhibit Pax6 expression. The most complex Pax6 expression pattern occurs in early test cell larvae of Wirenia (#3, #4, #5, and #6 in Fig. 8f, g ii, h). Since this expression pattern is neither complemented by coexpression of Pax2/5/8 or Paxβ nor supported by neurotransmitter distribution of 5-HT (serotonin) and FMRF-amide, the neurogenic nature of these Pax6 expressing cells is questionable. Remarkably, late test cell larvae of Wirenia exhibit two patches of Pax6-expressing cells at the anterior margin of the larval hyposphere and lateral to the ventral mouth opening, which is congruent with two patches of Paxβ-expressing cells (compare #7 in Fig. 8m, n, o and #4 in Fig. 9m, n, o). By comparing the polyplacophoran trochophore larva with the neomeniomorphan test cell larva we found that this specific location of the Pax6- and Paxβ-expressing cells correlates with the location of the polyplacophoran larval eyes and its associated Pax6 expression. Although, as far as currently known, neomeniomorph larvae and adults lack distinct photosensory organs (eyes), this correlation and the proposed close relationship of Polyplacophora and Neomeniomorpha calls for further morphological investigations combined with expression studies of suitable candidate genes that should focus on putative (rudimentary?) photoreceptors in this region.
The putative role of Paxα/β during bilaterian development
The Paxα/β group is the most recently discovered Pax subfamily and has tentatively been named “Pax?” . Initially, sequences of the “Pax?” subfamily have been shown to be present in Ctenophora, Deuterostomia (Echinodermata, Hemichordata), Nematoda, Annelida, and Platyhelminthes . Another study revealed a lophotrochozoan-specific Paxβ clade , which could also be assigned to the “Pax?” assemblage. Recently, a Paxα clade was described as supposed sister group to the lophotrochozoan-specific Paxβ clade . The Paxα clade comprises sequences of Panarthropoda as well as Deuterostomia, whereas the previously included ctenophore PaxA/B sequences received less support. Therefore, “Pax?” is considered a distinct bilaterian Pax subfamily termed “Paxα/β“ . The ancestral role of the Paxα/β subfamily in bilaterians remains unknown. However, the Paxβ expression data of the single investigated lophotrochozoan so far, the leech Helobdella, reveal expression of two Paxβ homologs, Hsp-Paxβ1 and Hsp-Paxβ2, in the developing CNS and, in case of Paxβ1, additionally at the site of eye formation . A similar situation was found in the onychophoran Euperipatoides rowelli, where Ero-Paxα is expressed in the CNS, more precisely in the lateral brain tissue associated with each eye of the adult animal . These findings are now complemented by our data on the expression of War-Paxβ in the region of the developing cerebral ganglia and the ventral nerve cords of the neomeniomorph mollusk W. argentea. Although the role of Paxα/β is obviously not exclusively assigned to developmental processes, the common expression in the CNS of all species hitherto investigated argues for an ancestral role in CNS specification at least in protostomes. Since data of Paxα expression in Deuterostomia are still lacking, the ancestral role of the Paxα/β subfamily for the entire Bilateria remains obscure. However, both investigated lophotrochozoans, W. argentea and Helobdella austinensis, show distinct Paxβ expression in the CNS of their anterior (“head”) region as well as in their ventral nerve cords (present study, ). The results for W. argentea are in accordance with the findings in another aculiferan mollusk, the polyplacophoran Acanthochitona crinita, where it is expressed in the lateral region of the cerebral commissure and in the ventral nerve cords (unpublished data). These congruent expression patterns in H. austinensis, A. crinita, and W. argentea argue for a conserved expression of Paxβ during neurogenesis of the ventral nerve cords of annelids and (aculiferan) mollusks.
Expression of Pax2/5/8, Pax6, and Paxβ and homology of lophotrochozoan ventral nerve cords
Our data provide first insights into the molecular background of aplacophoran neurogenesis on gene expression level and thus add significant insight into the putative roles of the herein investigated Pax genes in molluscan development. The Pax6 expression pattern in the aculiferan clades Neomeniomorpha and Polyplacophora largely resembles the common bilaterian expression during CNS development. The expression of Paxβ in the CNS of the neomeniomorph Wirenia argentea and the leech Helobdella austinensis argues for an ancestral role in lophotrochozoan neurogenesis, in particular during formation of the cerebral ganglia and the ventral nerve cords. Furthermore, we found indication for a conserved role of Pax2/5/8 in CNS development in Bilateria, while its expression in the spicule-secreting or associated cells in both neomeniomorphs and polyplacophorans suggests a novel function of this gene in aculiferan skeletogenesis. None of the investigated Pax genes are involved in the development of the lateral (visceral) nerve cords in either Neomeniomorpha or Polyplacophora, suggesting a different molecular background of this tetraneural subset on the one hand and homology of the ventral nerve cords within Spiralia (or Lophotrochozoa) on the other.
We thank Henrik Glenner (Department of Biology, University of Bergen, Norway) for providing boat time, laboratory space, and logistic support. We thank the crew of the RV Hans Brattström (University of Bergen) for assistance with collection of animals. MS thanks Martin Fritsch for providing valuable advice on confocal laser scanning microscopy and the IMARIS software. André Luiz de Oliveira gratefully acknowledges the financial support of CAPES and the Brazilian program Science without Borders (CsF – project 6090/13-3). This study was supported by a grant of the Austrian Science Foundation (FWF) to AW (grant number P24276-B22).
Availability of data and materials
All sequence data used in this data are available through GenBank.
MS and AW designed the study. MS designed the Wirenia argentea primers, amplified and cloned the gene fragments, synthesized the probes, performed in situ hybridization experiments, analyzed the data, created 3D reconstructions and illustrations, and drafted the manuscript. TW contributed the Pax gene sequences for the polyplacophoran Acanthochitona crinita. MS, ER, and TW extracted the RNA and synthesized the cDNA. ALO performed the transcriptome assembly and Pax gene phylogeny. MS, ER, CT, and TW reared the study material in Bergen, Norway. CT coordinated and supervised research in Bergen. AW supervised the project and contributed to data interpretation and writing of the manuscript. All authors provided input and read and approved the final version of the manuscript.
The authors declare that they have no competing interest.
Consent for publication
The entire study was carried out adhering to national and international laws. The species used in this study is neither protected nor endangered and work on it not restricted in any way.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Kocot KM, Cannon JT, Todt C, Citarella MR, Kohn AB, Meyer A, et al. Phylogenomics reveals deep molluscan relationships. Nature. 2011;477:452–6.View ArticlePubMedPubMed CentralGoogle Scholar
- 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–7.View ArticlePubMedGoogle Scholar
- Smith SA, Wilson NG, Goetz FE, Feehery C, Andrade SCS, Rouse GW, et al. Corrigendum: Resolving the evolutionary relationships of molluscs with phylogenomic tools. Nature. 2013;493:708.View ArticleGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Scherholz M, Redl E, Wollesen T, Todt C, Wanninger A. From complex to simple: myogenesis in an aplacophoran mollusk reveals key traits in aculiferan evolution. BMC Evol Biol. 2015;15:201.View ArticlePubMedPubMed CentralGoogle Scholar
- Todt C, Büchinger T, Wanninger A. The nervous system of the basal mollusk Wirenia argentea (Solenogastres): a study employing immunocytochemical and 3D reconstruction techniques. Mar Biol Res. 2008;4:290–303.View ArticleGoogle Scholar
- Wanninger A. Shaping the Things to Come: Ontogeny of Lophotrochozoan Neuromuscular Systems and the Tetraneuralia Concept. Biol Bull. 2009;216:293–306.View ArticlePubMedGoogle Scholar
- Sigwart JD, Sumner-Rooney LH. 18: Mollusca: Caudovofeata, Monoplacophora, Polyplacophora, Scaphopoda, and Solenogastres. In: Schmidt-Rhaesa A, Hartzsch S, Purschke G, editors. Structure and Evolution of Invertebrate Nervous Systems. Oxford: Oxford University Press; 2016. p. 172–89.Google Scholar
- Voronezhskaya EE, Croll RP. 20: Mollusca: Gastropoda. In: Schmidt-Rhaesa A, Hartzsch S, Purschke G, editors. Structure and Evolution of Invertebrate Nervous Systems. Oxford: Oxford University Press; 2016. p. 196–221.Google Scholar
- Wanninger A. 19: Mollusca: Bivalvia. In: Schmidt-Rhaesa A, Hartzsch S, Purschke G, editors. Structure and Evolution of Invertebrate Nervous Systems. Oxford: Oxford University Press; 2016. p. 190–5.Google Scholar
- Okusu A. Embryogenesis and Development of Epimenia babai (Mollusca, Neomeniomorpha). Biol Bull. 2002;203:87–103.View ArticlePubMedGoogle Scholar
- Todt C, Wanninger A. Of tests, trochs, shells, and spicules: Development of the basal mollusk Wirenia argentea (Solenogastres) and its bearing on the evolution of trochozoan larval key features. Front Zool. 2010;7:6.View ArticlePubMedPubMed CentralGoogle Scholar
- Redl E, Scherholz M, Todt C, Wollesen T, Wanninger A. Development of the nervous system in Solenogastres (Mollusca) reveals putative ancestral spiralian features. EvoDevo. 2014;5:48.View ArticlePubMedPubMed CentralGoogle Scholar
- Salvini-Plawen Lv. Early evolution and the primitive groups. In: Trueman ER, Clarke MR, editors. The Mollusca: Evolution (Volume 10). London: Academic; 1985. p. 59–150.Google Scholar
- Todt C. Salvini-Plawen Lv. Ultrastructure and histochemistry of the foregut in Wirenia argentea and Genitoconia rosea (Mollusca, Solenogastres). Zoomorphology. 2004;123:65–80.View ArticleGoogle Scholar
- Lindberg DR, Sigwart JD. What is the molluscan osphradium? A reconsideration of homology. Zool Anz. 2015;256:14–21.View ArticleGoogle Scholar
- Haszprunar G. Feinmorphologische Untersuchungen an Sinnesstrukturen ursprünglicher Solenogastres. Zool Anz. 1986;217:345–62.Google Scholar
- Wanninger A, Wollesen T. Lophotrochozoa (Spiralia): Mollusca. In: Wanninger A, editor. Evolutionary Developmental Biology of Invertebrates 2. Vienna: Springer Vienna; 2016. p. 103–53.Google Scholar
- LeClair E, Bonfiglio L, Tuan RS. Expression of the paired-box genes Pax-1 and Pax-9 in limb skeleton development. Dev Dynam. 1999;214:101–15.View ArticleGoogle Scholar
- Thompson JA, Ziman MR. Pax genes during neural development and their potential role in neuroregeneration. Progr Neurobiol. 2011;95:334–51.View ArticleGoogle Scholar
- Blake JA, Ziman MR. Pax genes: regulators of lineage specification and progenitor cell maintenance. Development. 2014;141:737–51.View ArticlePubMedGoogle Scholar
- Pfeffer PL, Gerster T, Lun K, Brand M, Busslinger M. Characterization of three novel members of the zebrafish Pax2/5/8 family: dependency of Pax5 and Pax8 expression on the Pax2.1 (noi) function. Development. 1998;125:3063–74.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Lichtneckert R, Reichert H. Insights into the urbilaterian brain: conserved genetic patterning mechanisms in insect and vertebrate brain development. Heredity. 2005;94:465–77.View ArticlePubMedGoogle Scholar
- 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–78.View ArticlePubMedGoogle Scholar
- Lang B, Shankland M. Expression of Pax-2/5/8 in Helobdella sp. Austin: Honor´s thesis; University of Texas, USA; 2010. http://hdl.handle.net/2152/13422. Accessed 30 Aug 2016.Google Scholar
- Wollesen T, Rodríguez Monje SV, Todt C, Degnan BM, Wanninger A. Ancestral role of Pax2/5/8 in molluscan brain and multimodal sensory system development. BMC Evol Biol. 2015;15:231.View ArticlePubMedPubMed CentralGoogle Scholar
- Vellutini BC, Hejnol A. Expression of segment polarity genes in brachiopods supports a non-segmental ancestral role of engrailed for bilaterians. Sci Rep. 2016;6:32387.View ArticlePubMedPubMed CentralGoogle Scholar
- Backfisch B, Rajan VBV, 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.View ArticlePubMedGoogle Scholar
- Hartmann B, Lee PN, Kang YY, Tomarev S, de Couet HG, Callaerts P. Pax6 in the sepiolid squid Euprymna scolopes: evidence for a role in eye, sensory organ and brain development. Mech Dev. 2003;120:177–83.View ArticlePubMedGoogle Scholar
- Navet S, Andouche A, Baratte S, Bonnaud L. Shh and Pax6 have unconventional expression patterns in embryonic morphogenesis in Sepia officinalis (Cephalopoda). Gene Expr Patterns. 2009;9:461–7.View ArticlePubMedGoogle Scholar
- Vöcking O, Kourtesis I, Hausen H. Posterior eyespots in larval chitons have a molecular identity similar to anterior cerebral eyes in other bilaterians. EvoDevo. 2015;6:40.View ArticlePubMedPubMed CentralGoogle Scholar
- Schmerer M, Savage RM, Shankland M. Paxβ: a novel family of lophotrochozoan Pax genes. Evol Dev. 2009;11:689–96.View ArticlePubMedGoogle Scholar
- Hill A, Boll W, Ries C, Warner L, Osswalt M, Hill M, et al. Origin of Pax and Six gene families in sponges: Single PaxB and Six1/2 orthologs in Chalinula loosanoffi. Dev Biol. 2010;343:106–23.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Redl E, Scherholz M, Wollesen T, Todt C, Wanninger A. Cell proliferation pattern and twist expression in an aplacophoran mollusk argue against segmented ancestry of Mollusca. J Exp Zool B Mol Dev Evol. 2016; 326:422–36.Google Scholar
- Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinf. 2009;10:1.View ArticleGoogle Scholar
- Katoh K, Stanley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Hall TA. BioEdit: A user-friendly biological sequence alignment program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.Google Scholar
- Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Darriba D, Taboada GL, Doallo R, Posada D. ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics. 2011;27:1164–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Jones DT, Taylor WR, Thornton JM. The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992;8:275–82.PubMedGoogle Scholar
- Rambaut A. FigTree-v1.4.2. 2014. http://treebioedacuk/software/figtree/. Accessed 10 Mar 2016.Google Scholar
- Jékely G, Arendt D. Cellular resolution expression profiling using confocal detection of NBT/BCIP precipitate by reflection microscopy. Biotechniques. 2007;42:751–5.View ArticlePubMedGoogle Scholar
- Altenburger A, Martinez P, Wanninger A. Homeobox gene expression in Brachiopoda: The role of Not and Cdx in bodyplan patterning, neurogenesis, and germ layer specification. Gene Expr Patterns. 2011;11:427–36.View ArticlePubMedGoogle Scholar
- Fritsch M, Wollesen T, de Oliveira AL, Wanninger A. Unexpected co-linearity of Hox gene expression in an aculiferan mollusk. BMC Evol Biol. 2015;15:151.View ArticlePubMedPubMed CentralGoogle Scholar
- Fritsch M, Wollesen T, Wanninger A. Hox and ParaHox Gene Expression in Early Body Plan Patterning of Polyplacophoran Mollusks. J Exp Zool (Mol Dev Evol). 2016;326:89–104.View ArticleGoogle Scholar
- Salvini-Plawen Lv. Was ist eine Trochophora? – Eine Analyse der Larventypen mariner Protostomier. Zool Jb Anat. 1980;103:389–423.Google Scholar
- Friedrich M. Evo-Devo gene toolkit update: at least seven Pax transcription factor subfamilies in the last common ancestor of bilaterian animals. Evol Dev. 2015;17:255–7.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- McCauley DW, Bronner-Fraser M. Conservation of Pax gene expression in ectodermal placodes of the lamprey. Gene. 2002;287:129–39.View ArticlePubMedGoogle Scholar
- Kavaler J, Fu W, Duan H, Noll M, Posakony JW. An essential role for the Drosophila Pax2 homolog in the differentiation of adult sensory organs. Development. 1999;126:2261–72.PubMedGoogle Scholar
- 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–57.View ArticleGoogle Scholar
- Haszprunar H, 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.View ArticlePubMedGoogle Scholar
- Dahl E, Koseki H, Balling R. Pax genes and organogenesis. Bioessays. 1997;19:755–65.View ArticlePubMedGoogle Scholar
- Bouchard M, Souabni A, Mandler M, Neubüser A, Busslinger M. Nephric lineage specification by Pax2 and Pax8. Genes Dev. 2002;16:2958–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Kocot KM, Aguilera F, McDougall C, Jackson DJ, Degnan BM. Sea shell diversity and rapidly evolving secretomes: insights into the evolution of biomineralization. Front Zool. 2016;13:23.View ArticlePubMedPubMed CentralGoogle Scholar
- Callaerts P, Halder G, Gehring WJ. Pax-6 in Development and Evolution. Annu Rev Neurosci. 1997;20:483–532.View ArticlePubMedGoogle Scholar
- Walther C, Gruss P. Pax-6, a murine paired box gene, is expressed in the developing CNS. Development. 1991;113:1435–49.PubMedGoogle Scholar
- Arendt D, Tessmar K, de Campos-Baptista M-I M, Dorresteijn A, Wittbrodt J. Development of pigment-cup eyes in the polychaete Platynereis dumerilii and evolutionary conservation of larval eyes in Bilateria. Development. 2002;129:1143–54.PubMedGoogle Scholar
- Callaerts P, Clements J, Francis C, Hens K. Pax6 and eye development in Arthropoda. Arthropod Struct Dev. 2006;35:379–91.View ArticlePubMedGoogle Scholar
- Denes AS, Jékely G, Steinmetz PRH, Raible F, Snyman H, Prud’homme B, et al. Molecular Architecture of Annelid Nerve Cord Supports Common Origin of Nervous System Centralization in Bilateria. Cell. 2007;129:277–88.View ArticlePubMedGoogle Scholar
- Quigley IK, Xie X, Shankland M. Hau-Pax6A expression in the central nervous system of the leech embryo. Dev Genes Evol. 2007;217:459–68.View ArticlePubMedGoogle Scholar
- Clements J, Hensa K, Merugu S, Dichtl B, de Couet HG, Callaerts P. Mutational analysis of the eyeless gene and phenotypic rescue reveal that an Intact Eyeless Protein is Necessary for normal Eye and Brain Development in Drosophila. Dev Biol. 2009;334:503–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Passamaneck YJ, Furchheim N, Hejnol A, Martindale MQ, Lüter C. Ciliary photoreceptors in the cerebral eyes of a protostome larva. EvoDevo. 2011;2:6.View ArticlePubMedPubMed CentralGoogle Scholar
- Eriksson BJ, Samadi L, Schmid A. The expression pattern of the genes engrailed, pax6, otd and six3 with special respect to head and eye development in Euperipatoides kanangrensis Reid 1996 (Onychophora: Peripatopsidae). Dev Genes Evol. 2013;223:237–46.View ArticlePubMedPubMed CentralGoogle Scholar
- Tomarev IS, Callaerts P, Kos L, Zinovieva R, Halder G, Gehring WJ, et al. Squid Pax-6 and eye development. Proc Natl Acad Sci U S A. 1997;94:2421–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Yoshida M-a, Yura K, Ogura A. Cephalopod eye evolution was modulated by the acquisition of Pax-6 splicing variants. Sci Rep. 2014;4:4256.PubMedPubMed CentralGoogle Scholar
- Kammermeier L, Leemans R, Hirth F, Flister S, Wenger U, Walldorf U, et al. Differential expression and function of the Drosophila Pax6 genes eyeless and twin of eyeless in embryonic central nervous system development. Mech Dev. 2001;103:71–8.View ArticlePubMedGoogle Scholar
- Blackburn DC, Conley KW, Plachetzki DC, Kempler K, Battelle B-A, Brown NL. Isolation and expression of Pax6 and atonal homologues in the American Horseshoe Crab, Limulus polyphemus. Dev Dyn. 2008;237:2209–19.View ArticlePubMedPubMed CentralGoogle Scholar
- Martín-Durán JM, Monjo F, Romero R. Morphological and molecular development of the eyes during embryogenesis of the freshwater planarian Schmidtea polychroa. Dev Genes Evol. 2012;222:45–54.View ArticlePubMedGoogle Scholar