- Open Access
Analyses of nervous system patterning genes in the tardigrade Hypsibius exemplaris illuminate the evolution of panarthropod brains
EvoDevo volume 9, Article number: 19 (2018)
Both euarthropods and vertebrates have tripartite brains. Several orthologous genes are expressed in similar regionalized patterns during brain development in both vertebrates and euarthropods. These similarities have been used to support direct homology of the tripartite brains of vertebrates and euarthropods. If the tripartite brains of vertebrates and euarthropods are homologous, then one would expect other taxa to share this structure. More generally, examination of other taxa can help in tracing the evolutionary history of brain structures. Tardigrades are an interesting lineage on which to test this hypothesis because they are closely related to euarthropods, and whether they have a tripartite brain or unipartite brain has recently been a focus of debate.
We tested this hypothesis by analyzing the expression patterns of six3, orthodenticle, pax6, unplugged, and pax2/5/8 during brain development in the tardigrade Hypsibius exemplaris—formerly misidentified as Hypsibius dujardini. These genes were expressed in a staggered anteroposterior order in H. exemplaris, similar to what has been reported for mice and flies. However, only six3, orthodenticle, and pax6 were expressed in the developing brain. Unplugged was expressed broadly throughout the trunk and posterior head, before the appearance of the nervous system. Pax2/5/8 was expressed in the developing central and peripheral nervous system in the trunk.
Our results buttress the conclusion of our previous study of Hox genes—that the brain of tardigrades is only homologous to the protocerebrum of euarthropods. They support a model based on fossil evidence that the last common ancestor of tardigrades and euarthropods possessed a unipartite brain. Our results are inconsistent with the hypothesis that the tripartite brain of euarthropods is directly homologous to the tripartite brain of vertebrates.
How brains evolved is one of the most perplexing questions in biology. Recent debates have centered on how to interpret similarities in brain development of distantly related animals [1,2,3,4,5,6,7,8,9]. During brain development in mice, Otx2 is expressed in the forebrain and midbrain. A paralogous group of paired box genes, Pax2, Pax5, and Pax8, along with the homeobox gene Gbx2, exhibit strong expression near the midbrain–hindbrain boundary. In the hindbrain and more posterior regions of the developing central nervous system, Hox genes are expressed . Orthologs of these genes are expressed in a similar staggered anteroposterior pattern in the tripartite brain of flies . Based on these correspondences, it has been hypothesized that the ancient ancestor of mice and flies—the last common ancestor of Nephrozoa —had a tripartite brain [2, 10, 12,13,14,15].
Similarities extend to specific regions of the brains of protostomes and deuterostomes. Six3 orthologs are expressed in the anteriormost region of the developing brain in representatives of both protostomes and deuterostomes—a region that gives rise to neurosecretory cells in both lineages . Pax6 is expressed in a lateral region of the brain in many protostomes and deuterostomes that have been investigated [17,18,19]. Structural and developmental similarities between the vertebrate pallium and annelid mushroom bodies , and between the vertebrate basal ganglia and the arthropod central complex  have also been identified. These similarities add support to the model of a nephrozoan ancestor with a complex brain [2, 3, 7, 9, 10, 14].
Objections against this view of brain evolution have been raised [4, 5, 8, 22]. Each main lineage of bilaterians includes representatives that lack complex brains, and several bilaterian lineages are characterized by diffuse nervous systems, rather than centralized nervous systems . Cladistic analyses suggest that the nephrozoan ancestor exhibited a nerve net and that centralized nervous systems evolved independently between five and nine times [4, 5], while brains evolved up to four times independently . Intriguingly, orthologs of the genes that pattern the tripartite brains of flies and mice exhibit regionalized expression patterns during development of the hemichordate Saccoglossus kowalevskii, even though hemichordates lack a brain, and instead exhibit a much more diffuse nervous system [24,25,26]. Likewise, genes that pattern tripartite brains show similar regionalized expression patterns in molluscs that exhibit diffuse nervous systems [27,28,29,30,31]. Therefore, genes that control development of the tripartite brains of insects and vertebrates exhibit conserved regionalized expression patterns in animals with diffuse nervous systems—animals that are predicted to have inherited their diffuse nervous systems from the last common ancestor of Nephrozoa [4, 5]. If so, this would indicate that these genes were independently co-opted for the evolution of tripartite brains several times in the nephrozoan lineage [1, 24, 25].
To trace the evolutionary history of complex nervous systems in Nephrozoa, additional taxa from diverse metazoan lineages must be investigated [1, 32, 33]. Studies of tardigrades may help elucidate nervous system evolution. Tardigrada is closely related to Euarthropoda [34,35,36] and, like euarthropods, tardigrades exhibit a segmented centralized nervous system (Fig. 1a) [37,38,39]. Several paired lobes are recognizable in the tardigrade brain, which have been interpreted as homologs of the proto-, deuto-, and tritocerebral brain segments of insects and other euarthropods [40,41,42,43]. By contrast, our recent analysis of Hox genes in the tardigrade Hypsibius exemplaris—formerly misidentified as H. dujardini and renamed H. exemplaris to reflect its burgeoning status as a model system —suggested that the brain of this species is homologous to just the protocerebrum of euarthropods  and therefore exhibits unipartite morphology. Distinguishing between these hypotheses has important implications for our understanding of the evolution of tripartite brain morphology.
To test the homology of the tardigrade brain to the brains of other animals, we investigated the expression patterns of several genes in H. exemplaris that have been implicated in the development of brains in both vertebrates and insects. Our results build on our analysis of Hox genes, demonstrating that tardigrades have a unipartite brain. These results support a model in which euarthropods evolved a tripartite brain after they diverged from tardigrades.
Maintaining H. exemplaris cultures
We used a standard protocol to culture H. exemplaris .
To visualize the developing nervous system, we stained H. exemplaris embryos with a β-tubulin antibody (E7, Developmental Studies Hybridoma Bank). To do this, we modified a method that was successfully implemented to stain the nervous system of hatchling  and adult H. exemplaris specimens  using the β-tubulin antibody, a method used to detect antibody localization in H. exemplaris embryos , and a method for embryonic in situ hybridization for this species . Staged embryos were washed for 1 h in a permeabilization buffer (5 units chitinase (Sigma-Aldrich C6137), 10 mg chymotrypsin (Sigma-Aldrich C4129), 1 ml 50 mM potassium phosphate buffer (pH 6.0)), followed by three 5-minute washes in 0.5X PBTw (0.5X phosphate-buffered saline, 0.1% Tween-20, pH 7.4). Embryos were fixed in 4% formaldehyde/33% heptane in 0.5X PBTw for 30 min at RT. Embryos were then washed five times for 5 min with 0.5X PBTw. Next, embryos were taken through a MeOH dilution series, consisting of 25, 50, 70, and 90% MeOH in 0.5X PBTw, followed by three washes in 100% MeOH. At this stage, embryos were kept in a − 20 °C freezer for at least 20 min. (Embryos can be stored indefinitely at this point.) Embryos were then taken through a reverse dilution series of MeOH and washed three times with 0.5X PBTw. Embryos were cut out of their eggshells with a 25-gauge needle and then washed three times for 10 min and four times for 30 min in 0.2% bovine serum albumin in 0.5X PBTw (BSA). Next, they were washed two times for 30 min in 5% normal goat serum in 0.5X PBTw (NGS). Embryos were then incubated overnight in a 1:100 dilution of the β-tubulin antibody in 5% NGS at 4 °C. The next day, embryos were washed three times for 5 min and four times for 30 min in 0.5X PBTw. This was followed by two 30-minute washes in NGS and an overnight wash in a 1:200 dilution of a goat anti-mouse Cy3-conjugated secondary antibody (Jackson ImmunoResearch) in NGS. The following day, embryos were washed three times for 5 min and six times for 30 min with 0.5X PBTw.
Identifying candidate genes and phylogenetic analyses
We performed reciprocal BLAST searches using human and D. melanogaster sequences as queries to identify candidate genes from H. exemplaris. Initial analyses focused on our draft genome assembly . Three studies identified containment sequences in our draft assembly for H. exemplaris [50,51,52]. Therefore, candidate genes that we identified with this method were then used as queries in BLAST search analyses of three independent transcriptome data sets [52,53,54] and two independent genome assemblies for H. exemplaris [51, 52], to verify that they were true H. exemplaris sequences. For phylogenetic analyses, we aligned sequences using MUSCLE . We trimmed our matrix using Gblocks . We performed both Bayesian and maximum likelihood analyses using the LG model  with an estimated proportion of invariable sites and an estimated gamma shape parameter with four substitution rate categories. We used PhyML for maximum likelihood analyses , with branch supports calculated by 500 bootstrap replicates. We used MrBayes for Bayesian analyses . The Bayesian analyses ran for 450,000 generations after the standard deviation of split frequencies dropped below 0.01, with trees sampled every 100 generations. Posterior probabilities were calculated from 4500 trees from the posterior tree distribution. We identified domains in the predicted protein sequences of candidate H. exemplaris orthologs with CD-Search . We identified intron/exon boundaries by comparing transcriptome sequences to genome sequences using Splign . When we could not find an ortholog of interest in H. exemplaris data sets, we tested for presence of the ortholog in the genome of Ramazzottius varieornatus  using reciprocal BLAST searches.
Genes of interest were amplified with PCR from H. exemplaris embryonic cDNA, or from H. exemplaris genomic DNA in the case of He-pax6. Primers used for this study were as follows: He-elav, 5′-GCATCCAGAACAAGAACATCAAGG-3′, 5′-ACTGGGAAAAGCGAAGTGTCTAGC-3′; He-otd, 5′-GTTCCCGCACCGAGGAAACAG-3′, 5′-CTCTCACGTCCTCCACGCTGA-3′; He-pax2/5/8, 5′-CGTTTTCCTTCAGACTTTCGTCGT-3′, 5′-TCCGATAACTCGTCTCGTTTCCTC-3′; He-pax6, 5′-CGTTTTATTTGCACACAGCGAGATA-3′, 5′-ATCTACCGGATTGCAAAGTTCTGG-3′; He-six3, 5′-ATCTTCACTTGACGCGATTGTGGT-3′, 5′-GTCCTTGCTGTTATCCTCGTCCAT-3′; He-unplugged (unpg), 5′-TTGCGAGAGAAACAAAACTGGATG-3′, 5′-CAAAACAAACGCGCCAAGTG-3′; An-elav outer, 5′-AAACGATGACACAAGACGAAATTA-3′, 5′-CAGCATACCTGTACCTATTCATGG-3′; An-elav inner, 5′-GTTCAGTTGGTGCTATTGAGTCAT-3′, 5′-AGGTAGATAGGATGGTGCTAATGG-3′; An-pax2/5/8 outer, 5′-TATAAACATTTGGTGGAGACGACA-3′, 5′-GAATCAATCAACTTGGAGGAGTTT-3′; An-pax2/5/8 inner, 5′-ATGGTGAGTGCGACTATCATCTTCG-3′, 5′-ATGGTGAGTGCGACTATCATCTTC-3′, 5′-ACATCAGTCGACAGTTACGAGTGT-3′; An-poxm outer, 5′-CGATTTGCATAACGTAACGTACTC-3′, 5′-GAACTGGAAAGAAATGATCGAACT-3′; An-poxm inner, 5′-GCAGACTTTTATTGATGTTGTTCG-3′, 5′-TTTCAAAGTGATTCAAACCAAGAA-3′. Genes were cloned into the pCR™4-TOPO® (Invitrogen) vector following the manufacturers protocol. Minipreps for each clone were made using a QIAprep Spin Miniprep Kit (Qiagen) and following the manufacturers protocol. Miniprep plasmids were sequenced to verify their identity.
In situ hybridization
To make in situ probes, we linearized plasmids with PCR by using the M13 forward and reverse primers that come with the TOPO™ TA Cloning™ Kit for Sequencing (Invitrogen). Linearized template was used in a transcription reaction using either T7 or T3 RNA polymerase (Promega), DIG RNA labeling mix (Roche), and following the manufacturer’s protocols. Synthesized probes were cleaned using an RNAeasy® Mini Kit and following the manufacturer’s instructions. In situ hybridization was performed as previously described . For counter staining the nervous system of in situ hybridization specimens, the β-tubulin antibody was added after we washed out the anti-DIG-AP (Roche) antibody. Instead of washing out the anti-DIG-AP antibody with maleic acid buffer, as we normally do for in situ hybridizations, for counter staining with the β-tubulin antibody, we washed it out with 0.5X PBTw. After these washes, we followed the β-tubulin immunostaining method described above. After completion of the immunostaining procedure, we washed the embryos with maleic acid buffer (pH 9.5) and continued through the development steps of the in situ hybridization method.
After in situ hybridization and/or β-tubulin immunostaining, embryos were mounted in DAPI Fluoromount-G (SouthernBiotech). Fluorescent images were collected on a Zeiss 710 LSM. Maximum projections were produced in ImageJ. DIC images were taken on a Nikon Eclipse 800 microscope. Minimum and maximum displayed pixel values were adjusted in ImageJ and/or Photoshop. In cases where the in situ stain occluded β-tubulin or DAPI signal, we used either the Cyan Hot or glow LUT in ImageJ to better visualize the β-tubulin or DAPI signal.
Identification of an elav ortholog
ELAV/HuD family splicing factors are RNA-binding proteins that are almost exclusively expressed in differentiating and mature neurons . We used elav expression to identify early stages in nervous system development in H. exemplaris embryos. We identified two candidate elav genes in our H. exemplaris nucleotide data. One of these predicted genes was present in all data sets that we analyzed, except for the Yoshida et al.  transcriptome (Additional file 1: Table S1). We refer to this sequence as He-elav. The other candidate elav ortholog was missing in all other data sets that we analyzed except for our H. exemplaris genome assembly (Additional file 2: Table S2) . We were unable to amplify this gene from genomic DNA using nested PCR (Additional file 3: Fig. S1). Furthermore, this predicted gene sequence is located on a small scaffold of only 7049 nucleotides in our genome assembly and it is the only gene predicted to be located on this scaffold. Therefore, we concluded that this sequence most likely represents a contaminant in our genome assembly. We call this sequence anonymous elav (An-ELAV), referring to its likely inclusion as a contaminant from an unknown organism.
We performed phylogenetic analyses on a matrix of 239 amino acid positions from 14 ELAV sequences stemming from 10 species, in addition to He-ELAV and An-ELAV (Additional file 1: Table S1; Additional file 2: Table S2; Additional file 4: FASTA alignments). Both of these sequences were nested within a well-supported clade (470/500 bootstrap support; 1.0 posterior probability) that separated all deuterostome elav genes from all protostome elav genes (Fig. 1a). We interpret this as phylogenetic support for the assignment of both sequences to the elav orthology group. A CD-Search detected an ELAV/HuD family splicing factor domain (Accession TIGR1661; E-value = 6.88e-123) characteristic of ELAV proteins in the He-elav sequence that we identified (Fig. 1b, c). Interestingly, 37 amino acids within the ELAV/HuD family splicing factor domain were not recognized as part of this domain by the CD-Search algorithm (between amino acid positions 250 and 287, Fig. 1c). We also detected the nucleotide sequence that is predicted to give rise to this intervening region in the He-elav sequence that we amplified from embryonic cDNA.
Expression pattern of He-elav and β-tubulin
We used a β-tubulin antibody to visualize the nervous system of H. exemplaris. The central nervous system of H. exemplaris includes a brain that is housed in the head and ventral trunk nervous system . Inner connectives (ic) extend from the inner brain region (ib) to the anteriormost trunk ganglion (ga1; Fig. 2a) . The inner brain region is composed of neuropil . Outer connectives extend between the outer brain region (ob) and the anteriormost trunk ganglion (Fig. 2a) . Neurites of the inner and outer brain regions meet at a thick band of dorsal neuropil . Trunk ganglia adjoin adjacent trunk ganglia by paired connectives (cn) . We detected this basic nervous system architecture with the β-tubulin antibody in 45 h post-laying (hpl) embryos (Additional file 5: Movie 1). The nervous system of our 45 hpl embryos appears to be in the final stage of development, based on a staging series for H. exemplaris .
For in situ hybridizations, we focused on three developmental stages: 24 hpl, when evidence of ectodermal segmentation is apparent, but the nervous system and legs are not; 35 hpl, when developing ganglia are apparent, and leg buds are discernible; and 45 hpl, when legs have reached their final size and the central nervous system exhibits the connectivity that is visible post-embryonically (Additional file 5: Movie 1) . Strong He-elav staining in 45 hpl embryos was restricted to the region where the central nervous system develops (Fig. 2b–d). The He-elav in situ pattern closely matched the pattern of β-tubulin localization in 45 hpl embryos (compare Fig. 2b–d to Additional file 5: Movie 1). At 24 hpl, expression of He-elav staining was noticeably variable between specimens (Fig. 2e). Only one specimen exhibited He-elav staining near the future position of the central nervous system (black arrowhead, Fig. 2e), and He-elav staining was relatively light in this specimen compared to 35 hpl and 45 hpl embryos (Fig. 2f, g). This indicates that specification of neuronal identity begins in embryos at 24 hpl or shortly after.
Identification of a six3 ortholog
The Six family includes three groups of genes—Six1/2/sine oculis, Six3/6/optix, and Six4/5 . We identified three-candidate Six family genes in available H. exemplaris genome and transcriptome assemblies (Additional file 1: Table S1; Additional file 2: Table S2). We performed a phylogenetic analysis using a matrix that included 164 amino acid positions (Additional file 4: FASTA alignments). Each candidate gene from H. exemplaris fell within a different Six family group in our phylogenetic analysis. If we rooted the tree with the candidate gene that appeared to be nested within the Six3/6/optix group, this group was no longer monophyletic. This suggests that this candidate gene represents an H. exemplaris ortholog of the Six3/6/optix group. This gene contains a predicted SIX domain (SIX1_SD) and a homeodomain (Fig. 3b, c), both of which are characteristics of Six family genes . Because the direct ortholog of Six3/optix is referred to as six3 in most animals, we refer to the H. exemplaris ortholog as He-six3.
Expression pattern of He-six3
Strong He-six3 signal was detected during all developmental periods that we investigated (Fig. 4a–c). At 24 hpl, He-six3 signal was detected broadly across the ectoderm in the developing head (Fig. 4d, g). He-six3 signal did not extend to the posterior border of the dorsal head at 24 hpl (Fig. 4d). At 35 hpl, He-six3 transcripts were detected in regions of the head that will give rise to both the outer (ob) and inner (ib) brain regions (Fig. e, h). We also detected transcripts near the developing stomodeal complex in 35 hpl embryos (st; Fig. 4e). At 45 hpl, He-six3 signal was detected broadly across the lateral ectoderm of the head (Fig. 4f). We also detected He-six3 signal in the region of the inner brain where neuropil develops, but not in the dorsal neuropil (dnp; compare Fig. 4j to frames 19–27 of Additional file 5: Movie 1).
Identification of an otd ortholog
We previously identified an otd ortholog in H. exemplaris . In order to better characterize this gene from a phylogenetic perspective, we performed an analysis that included He-otd and candidate H. exemplaris PRD Class homeobox genes from our genome assembly (Additional file 1: Table S1; Additional file 2: Table S2). We followed a previously published naming scheme for PRD Class homeobox families . Our matrix included the 60 amino acid long homeodomains of 74 PRD Class homeobox proteins stemming from several species (Additional file 4: FASTA alignments). In the maximum likelihood tree, He-OTD was nested within a highly supported monophyletic OTX clade (Fig. 5a; 432/500, bootstrap support; 1.00 Bayesian posterior probability; Additional file 6: Fig. S2). The homeodomain was the only conserved domain present in the predicted He-OTD protein (Fig. 5b, c), as is the case for D. melanogaster OTD. We did not detect Pax3/7 prd or Pax eyg subfamily members in any of the H. exemplaris databases that we analyzed. We analyzed the genome assembly of a second tardigrade species—R. varieornatus  using reciprocal BLAST searches, but did not find matches for either of these Pax subfamily members.
Expression pattern of He-otd
At all stages that we analyzed, He-otd signal was strongly localized to the developing head (Fig. 6a–j). At 24 hpl, He-otd signal was detected broadly throughout the ectodermal layer of the head, extending to near the posterior border of the head (Fig. 6d) . At 35 hpl, He-otd signal appeared strongest within the head, including in or near the developing stomodeal complex (st; Fig. 6e). At 45 hpl, we detected strong He-otd signal in the inner brain (ib) and outer brain (ob) neuropil (compare Fig. 6f to frames 11–27 of Additional file 5: Movie 1). He-otd signal colocalized with β-tubulin signal in the inner brain region (Fig. 6j). Strong He-otd signal was also detected in the developing stomodeal complex at 45 hpl (Fig. 6j).
Identification of a pax6 ortholog
PAX6 proteins generally contain a paired box domain (PAX) and a homeodomain [68, 69]. We identified a single H. exemplaris ortholog of pax6 in our analysis of PRD Class homeodomains (Fig. 7a). He-PAX6 was nested within a well-supported monophyletic PAX6 clade in this tree (416/500, bootstrap support; 1.00 Bayesian posterior probability). We also performed phylogenetic analyses on a matrix that included 112 amino acids of the paired box domain (Additional file 4: FASTA alignments). This same gene was also nested within a monophyletic PAX6 clade in our analysis of PAX domains (Fig. 7a; Additional file 1: Table S1; Additional file 2: Table S2). Support for this clade was low in our PAX domain analysis, but we recovered good support for the more inclusive clade PAX4/6 (492/500, bootstrap support; 1.00 Bayesian posterior probability). A Pax domain and a homeodomain were the only conserved domains that we detected in He-PAX6 (Fig. 7b, c). We also detected an ortholog of paxα (Fig. 7a), a gene that is predicted to have been present in the nephrozoan ancestor, but independently lost in the chordate and Drosophila lineages . We did not include PAX EYG proteins in our analysis of PAX domains because the D. melanogaster PAX EYG members—Eyegone and Twin of Eyegone—have highly divergent PAX domains, and we had already established that H. exemplaris was missing a pax eyg ortholog (Fig. 5a). As with our analysis of PRD Class homeodomains (Fig. 5a), we did not detect an H. exemplaris PAX3/7 member in our analysis of PAX domains (Fig. 7a).
Expression pattern of He-pax6
At 24 hpl, the in situ pattern was consistent between embryos, but signal intensity appeared variable (Fig. 8a), suggesting that there are dynamic changes in expression levels of He-pax6 at this development stage. We did not detect He-pax6 transcripts at 35 hpl or 45 hpl (Fig. 8b, c), when nervous system morphology is apparent and He-elav signal is strong (Fig. 2). At 24 hpl, He-pax6 signal was detected in a relatively large lateral region of the outer ectoderm of the developing head (Fig. 8d–f). We did not detect He-pax6 signal in the internal layer of cells in the developing head, which we predict give rise to the foregut (fg; Fig. 8e). In the developing trunk, we detected He-pax6 signal in paired ventromedial domains in the ectodermal cell layer (Fig. 8e, f).
Identification of an unpg ortholog
Gbx/unpg codes for an ANTP class homeodomain protein. In order to identify an ortholog of Gbx/unpg, we looked for ANTP class homeobox genes in available H. exemplaris sequence data, using Ryan et al.  as a guide to choosing human and D. melanogaster ANTP class homeodomain proteins to use in our BLAST searches. We identified 35 candidate ANTP class homeobox genes with this method. We performed phylogenetic analyses on a matrix of 60 amino acids of the ANTP homeodomain from 167 sequences (Additional file 1: Table S1; Additional file 2: Table S2; Additional file 4: FASTA alignments). Of these candidates, one predicted protein sequence was most closely related to a monophyletic group of GBX/UNPG proteins (Fig. 9a; Additional file 7: Fig. S3). Together with the GBX/UNPG clade, this sequence formed a monophyletic group with good support (366/500, bootstrap support; 1.00 Bayesian posterior probability). Although it is not nested within the clade of previously identified GBX/UNPG sequences, it is a reciprocal best blast hit to both Homo sapiens and D. melanogaster GBX/UNPG sequences. Like human GBX1 and GBX2 and D. melanogaster UNPG, the homeodomain is the only conserved domain in the candidate GBX/UNPG sequence (Fig. 9b, c). We refer to this gene as He-unpg, since it is referred to as unpg in onychophorans  and flies .
Expression pattern of He-unpg
We detected strong He-unpg signal at 24 hpl (Fig. 8g). He-unpg signal at 35 hpl and 45 hpl was either weak or absent (Fig. 8h, i), suggesting that it does not play an important role in later stages of nervous system development. At 24 hpl, He-unpg signal was detected broadly throughout the trunk (Fig. 8j–l). He-unpg signal was mostly absent in the developing head, but was detected in a relatively small dorsoposterior domain (Fig. 8j). Signal was detected most strongly between the nuclei of the ectodermal layer and the nuclei of the endomesodermal layer (em), but was still detectible within the endomesodermal layer (Fig. 8k).
Identification of a pax2/5/8 ortholog
We identified two candidate pax2/5/8 orthologs in the initial BLAST search analysis of our genome assembly . One of these candidates was found in all data sets except the transcriptome assemblies from Yoshida et al.  and Levin et al. . We refer to this sequence He-pax2/5/8. The second candidate sequence was not identified in any other database that we analyzed (Additional file 2: Table S2). We were unable to amplify this sequence from genomic DNA using nested PCR (Additional file 3: Fig. S1). Furthermore, this predicted gene sequence is located on a small scaffold of just 2867 nucleotides in our genome assembly. It is the only predicated gene sequence on this scaffold. Therefore, we predict that this sequence represents a contaminant in our genome assembly. We refer to this sequence as anonymous pax2/5/8 (An-pax2/5/8). We detected another PAX sequence in our genome assembly  that was not found in any other database. This predicted sequence was recovered in a highly supported PAX1/9 clade in our phylogenetic analysis (Fig. 7a). We were unable to amplify this gene from genomic DNA using nested PCR (Additional file 3: Fig. S1). Therefore, we refer to this sequence as anonymous poxm (An-poxm).
We analyzed a matrix of Pax genes that included 112 amino acids (Additional file 4: FASTA alignments). He-PAX2/5/8 and An-PAX2/5/8 were both nested within a monophyletic PAX2/5/8 clade in our analysis of PAX domains (Fig. 7a), supporting our assignment of these sequences to the pax2/5/8 subfamily. The PAX domain was the only domain that we detected in He-pax2/5/8 (Fig. 7d, e). Likewise, a PAX domain was the only conserved domain detected in CD-Search analyses of the D. melanogaster and Euperipatoides rowelli (Onychophora) pax2/5/8 orthologs that we used in our phylogenetic analyses.
Expression pattern of He-pax2/5/8
He-pax2/5/8 signal was detected at all developmental stages of our study (Fig. 10a–c). Strong He-pax2/5/8 signal was detected in the developing trunk at all stages investigated (Fig. 10a–i). At 24 hpl, strong He-pax2/5/8 signal was detected between the ectodermal layer and endomesodermal layer and within the endomesodermal layer of the developing trunk (Fig. 10d). At 35 hpl, He-pax2/5/8 signal was detected in the developing leg buds (L1–L2; Fig. 10e). At 45 hpl, He-pax2/5/8 signal was detected in the developing trunk ganglia (ga1–ga4; Fig. 10f) and legs (L4; Fig. 10c; L1–L3; Fig. 10j).
An earlier reconstruction of the nervous system of the ancient nephrozoan ancestor relied on developmental data from just two species—mice and flies . Analyses of additional taxa are required to more confidently resolve the evolution of nervous systems in Nephrozoa [1, 32, 33]. D. melanogaster is part of the highly diverse lineage Euarthropoda. Euarthropoda is part of a larger lineage called Panarthropoda that includes Tardigrada and Onychophora. Like euarthropods, both tardigrades [37, 41, 42, 47, 71,72,73,74] and onychophorans [75,76,77,78,79] exhibit centralized nervous systems and complex brains. Reconstructing the evolution of nervous systems in Panarthropoda requires consideration of nervous system development in representatives of all three panarthropod lineages. According to two phylogenomic studies of Panarthropoda that carefully controlled for long branch attraction artifacts, Tardigrada is the outgroup of a euarthropod + onychophoran lineage [34, 36], which makes the tardigrade lineage, especially important for reconstructing nervous system evolution in Panarthropoda. Here, we compare the results of our study to studies of euarthropods and onychophorans to reconstruct nervous system evolution in Panarthropoda.
Tardigrades have a unipartite brain
To determine whether tardigrades have a tripartite or unipartite brain, it is necessary to determine how the head of tardigrades—the brain housing unit—relates to segments of other panarthropods. Previous studies have compared expression patterns of developmental genes to homologize segments between chelicerates and mandibulates [80,81,82] and between euarthropods and onychophorans [83, 84]. The expression patterns of developmental genes in anterior segments are remarkably similar across the euarthropod/onychophoran clade (Fig. 11a). Six3 is expressed in the anterior part of the first segment [16, 85]. Otd is expressed broadly in a more posterior domain of the first segment and in a more restricted mid-ventral region of more posterior segments [82, 85,86,87,88,89]. Pax6 is expressed in a large lateral domain in the first segment and in more restricted ventral domains in more posterior segments [17, 19, 85, 89,90,91]. The anterior expression border of unpg is in the second segment in D. melanogaster [10, 92] and onychophorans . Pax2/5/8 is expressed segmentally in D. melanogaster  and onychophorans , with the anteriormost expression domain located in the first segment. The anterior expression borders of labial (lab), proboscipedia (pb), and Hox3 are in the third segment, Deformed (Dfd) in the fourth segment, and fushi tarazu (ftz) in the fifth segment in both euarthropods [93, 94] and onychophorans [83, 84]. Results of these studies indicate that—in terms of homology—anterior segments of euarthropods and onychophorans align one to one in anterior–posterior order.
The tripartite brain hypothesis suggests that the tardigrade head is homologous to the first three segments of euarthropods and onychophorans . This hypothesis predicts that genes that are expressed in the first three segments of euarthropods and onychophorans should be expressed in the head of tardigrades. Our analyses of gene expression patterns in H. exemplaris contradict the tripartite hypothesis. The tripartite hypothesis predicts that six3 should only be expressed in the anterior part of the first segment of a three-segment head (Fig. 11a). However, He-six3 is expressed broadly in the head (Fig. 4), rather than being restricted to a small anterior domain. The tripartite hypothesis predicts that broad otd expression should be restricted to the first segment of a three-segment head (Fig. 11a). By contrast, broad expression of He-otd continues to the posterior border of the head (Fig. 6d). The same contradiction applies for pax6 (Fig. 7a, b). The tripartite hypothesis predicts that both labial and Hox3 should be expressed in the tardigrade head (Fig. 11a). In actuality, neither of these genes is expressed in the head of tardigrades . Furthermore, this hypothesis predicts that the anterior expression border of Dfd and ftz should lie in the first and second trunk segment, respectively, in tardigrades (Fig. 11a). However, the anterior expression borders of these genes lie in the third and fourth trunk segments in tardigrades .
The unipartite hypothesis suggests that the tardigrade head is homologous to just the first segment of euarthropods and onychophorans . In contrast to the tripartite brain hypothesis, the unipartite brain hypothesis correctly predicts expression patterns of several developmental genes in H. exemplaris. This hypothesis predicts that six3, otd, and pax6 should be expressed broadly in the head of tardigrades (Fig. 11a), as is the case (Figs. 4, 6, 8a–c), although He-six3 is even more broadly expressed than predicted (see below). Furthermore, this hypothesis correctly predicts the anterior expression domains of Hox3, Dfd, and ftz (Fig. 11a).
While the expression patterns of some genes support the unipartite hypothesis, expression patterns of He-unpg, He-lab, and He-pax2/5/8 are not predicted by either the tripartite hypothesis or the unipartite hypothesis. The tripartite hypothesis predicts that anterior border of unpg expression should lie in the middle of the tardigrade head (Fig. 11a), a position that would correspond to the middle segment of a three-segment head. The unipartite hypothesis predicts that the anterior border of unpg should lie in the first trunk segment of tardigrades (Fig. 11a). However, the anterior border of He-unpg expression is in a posterior region of the developing head (Fig. 8j). Although this result is not predicted by either hypothesis, to our knowledge, in Panarthropoda, unpg expression has only been investigated in D. melanogaster and the onychophoran Euperipatoides kanangrensis. Even between these species, there appears to be variation in the exact position of the anterior border of unpg expression. Therefore, unpg expression should be investigated in additional panarthropod representatives to determine whether the anterior expression boundary of He-unpg is atypical. Concerning He-lab, the tripartite brain hypothesis predicts that the anterior expression boundary of this gene should be in the tardigrade head, in a region that corresponds to the third brain segment of euarthropods (Fig. 11a). The unipartite hypothesis predicts that the anterior expression boundary of lab should lie in the second trunk segment in tardigrades (Fig. 11a). In actuality, the anterior boundary of He-lab expression is in the first trunk segment . If the unipartite hypothesis is correct, as suggested by the expression domains of several other genes (see above), then the anterior expression boundary of lab must have expanded into a more anterior segment in the tardigrade lineage, or retracted into a more posterior segment in the lineage leading to Euarthropoda and Onychophora.
Lastly, the tripartite hypothesis predicts that there should be two or three segmental expression domains of pax2/5/8 in the tardigrade head, based on the expression patterns of this gene in D. melanogaster and onychophorans (Fig. 11a). The unipartite hypothesis predicts that He-pax2/5/8 should exhibit a single segmental expression domain in the head (Fig. 11a). Matching neither hypothesis, strong He-pax2/5/8 expression is restricted to the developing trunk (Fig. 10). In the case of He-pax2/5/8, comparisons to more distantly related animals reveal clues about the composition of the tardigrade brain. In most non-panarthropod animals that have been investigated, the anterior border of pax2/5/8 abuts the posterior border of otx/otd expression . In H. exemplaris, the anteriormost border of strong pax2/5/8 expression closely aligns with the posterior border of strong He-otd expression where the developing ventral nerve cord meets the developing brain (compare Fig. 10c to Fig. 6f). Therefore, concerning the discrepancy between the anteriormost expression domain of pax2/5/8 of H. exemplaris and other panarthropods, H. exemplaris most likely retains the ancestral condition. In sum, the weight of evidence supports the unipartite hypothesis for the composition of the tardigrade brain. Changes in expression domains of developmental genes have previously been implicated in the diversification of animal body plans [95,96,97], raising the possibility that changes in the expression domains of lab, unpg, and pax2/5/8 played a role in the diversification of panarthropod body plans.
Six3 and otd patterns overlap extensively in the developing brain of H. exemplaris
Across Nephrozoa, expression and function of six3 and otd/otx overlap minimally during development of the anterior nervous system—the brain of many animals . Six3 is expressed in the anteriormost median domain of the body axis of nephrozoans [29, 85, 98], where it is thought to play an ancient conserved role in specifying neurosecretory cells . Otx/otd is typically expressed in a more posterolateral part of the developing brain, where it regulates eye development [16, 99,100,101]. In contrast to other animals, He-six3 appeared to be expressed broadly across the developing head, including in lateral regions (Fig. 4), rather than being restricted to an anteromedial domain. It appeared that expression of He-six3 and He-otd broadly overlapped in H. exemplaris, rather than being restricted to nearly non-overlapping expression domains, as seen in most other animals . In fact, we detected expression of both genes in the developing inner brain region and stomodeal complex (Figs. 4e, f, j; 6e, f, j). Furthermore, in the outer parts of the head, strong He-otd signal was restricted to brain neuropil (Fig. 6f, inset), while He-six3 was more broadly expressed (Fig. 4e, f). Comparing our results to a previous study of euarthropods, onychophorans, and annelids , it appears that six3 expression expanded from an anteromedial domain into more posterior and lateral regions of the developing brain in the tardigrade lineage, after this lineage diverged from other animals. Determining where in tardigrade phylogeny expansion of six3 evolved will require studies of additional tardigrade species, including species within Heterotardigrada—the most distant tardigrade relatives of H. exemplaris. Additionally, the exact location of the cell bodies where otd and six3 are expressed in tardigrades would further illuminate the composition of the tardigrade brain.
Pax6 and unpg exhibit unique temporal dynamics
The canonical nervous system patterning genes that we investigated in this study all exhibited strong expression during the earliest stage that we investigated, 24 hpl. Of these genes, only two—He-pax6 and He-ungp—were not expressed strongly at the later stages that we investigated (Fig. 8b, c, h, i). Orthologs of the canonical nervous system patterning genes, including pax6 and unpg orthologs, continue to be expressed during development of onychophorans and euarthropods after appendages and the central nervous system are apparent [10, 16,17,18,19, 70, 82, 85,86,87,88,89,90,91,92], stages that are presumably later than the stages when pax6 and unpg are expressed in H. exemplaris, since legs and the nervous system are not morphologically visible at 24 hpl in this species. While it is likely that pax6 and ungp play roles in setting up the general regionalized pattern of the H. exemplaris nervous system, our results suggest that—unlike in onychophorans and euarthropods—these genes do not play later roles in nervous system patterning, such as specifying or maintaining neural cell types . However, it is possible that these genes are expressed at later stages at a low level that is difficult to detect with our in situ hybridization method, but biologically significant. Testing this possibility will require establishment of new methods of visualizing gene expression in H. exemplaris and functional data for these genes.
The evolution of the panarthropod Pax family
A recent analysis of Pax family genes that included transcriptome data from H. exemplaris and the onychophoran E. rowelli revealed a highly supported gene clade that was referred to as paxα . Based on this study, it appears that the nephrozoan ancestor had at least seven Pax gene subfamilies : (1) Pax alpha, (2) Pax eyg, (3) Pox neuro, (4) Pax4/6/10, (5) Pax2/5/8, (6) Pax3/7 prd, and (7) Pax1/9 meso. We identified Pax family genes in two analyses—an analysis of PRD Class homeodomains (Fig. 5), which are present in some Pax gene subfamilies [69, 102], and an analysis of PAX domains (Fig. 7). Our analyses suggest that H. exemplaris possesses single orthologs of each Pax gene subfamily except Pax eyg and Pax3/7. We were surprised to not find a pax3/7 ortholog in any nucleotide database for H. exemplaris, since an ortholog of this gene is strongly predicted to have been present in the last common ancestor of Nephrozoa and is present in non-tardigrade panarthropods [69, 90, 102]. Additionally, we did not detect a pax3/7 ortholog in the genome of the tardigrade R. varieornatus. Our results suggest that pax3/7 was deleted somewhere in the tardigrade lineage after it split from the lineage leading to euarthropods and onychophorans, but before the divergence of the H. exemplaris and R. varieornatus lineages.
The conclusion that a Pax eyg ortholog was present in the nephrozoan last common ancestor is a relatively recent idea. This idea stems from a phylogenetic analysis of Pax family genes that recovered a clade of insect pax eyg genes and pax eyg candidates from sea urchins and hemichordates . However, there does not appear to be a pax eyg ortholog in either tardigrades or onychophorans . To our knowledge, pax eyg orthologs have not been identified in any animals outside of sea urchins, hemichordates, and Holometabola. Therefore, we predict that pax eyg evolved in the lineage leading to holometabolous insects, after this lineage split from other insects. In this view, previously identified pax eyg orthologs in sea urchins and hemichordates evolved independently from the pax eyg gene of Holometabola. The nomenclature of Pax genes should be modified to reflect the fact that the previously identified sequences referred to as pax eyg are not directly orthologous across Nephrozoa.
The tripartite brain hypothesis suggests that the forebrain, midbrain, and hindbrain of vertebrates are directly homologous to the proto-, deuto-, and tritocerebral brain segments of flies. Therefore, this hypothesis suggests that the ancestor of vertebrates and flies—the last common ancestor of Nephrozoa—also exhibited a tripartite brain [2, 10, 12,13,14,15]. Based on this hypothesis, other nephrozoans should also exhibit tripartite brains. While all euarthropods are generally interpreted as possessing tripartite brains , developmental  and morphological  data suggest that Onychophora—the sister lineage of Euarthropoda [34, 36]—is characterized by a bipartite brain. Additionally, several stem representatives of both onychophorans and euarthropods exhibited unipartite brains [39, 104]. A comprehensive model of panarthropod brain evolution, based on investigations of extant and extinct representatives of Panarthropoda, presents a more parsimonious view of panarthropod brain evolution compared to the tripartite hypothesis [39, 104]. Under this model, the last common ancestor of Euarthropoda and Onychophora is predicted to have possessed a unipartite brain homologous to the protocerebrum of modern euarthropods. More complex brains evolved independently along the onycophoran and euarthropod lineages, according to models that include data from fossils [39, 104]. In the onychophoran lineage, the central nervous system of the second segment fused to the protocerebral brain, forming a bipartite brain. In the euarthropod lineage, the central nervous system of the second segment fused to the ancestral protocerebral brain, giving rise to the deutocerebral brain segment. The central nervous system of the third segment fused to the deutocerebrum, giving rise to the tritocerebral brain segment. Conceivably, tardigrades could have also independently evolved a multipartite brain—a possibility that draws support from the fact that the evolution of small body size is often associated with the fusion of segmental ganglia . However, our analyses of brain pattering genes in H. exemplaris suggest that tardigrades possess a unipartite brain. Even if the tardigrade brain is currently unipartite, we need to consider the possibility that the tardigrade brain was secondarily simplified from a multipartite brain during the evolution of the tardigrade lineage. Although we cannot completely rule out this possibility, it is a less parsimonious view of brain evolution. Our results place the origin of the unipartite brain in the stem lineage leading to Panarthropoda and suggest that tardigrades retain the ancestral protocerebral-grade brain [39, 104]. This model of brain evolution is inconsistent with an ancestral tripartite nephrozoan brain.
Although the ancestor of Panarthropoda most likely exhibited a unipartite brain, it is clear based on comparisons between H. exemplaris and other panarthropods that its nervous system patterning genes exhibited regionalized expression patterns. Our study supports the model in which regionalized expression patterns of vertebrate and fly brain patterning genes evolved in a nephrozoan ancestor that did not possess a tripartite brain and that these regionalized patterns were co-opted independently for the evolution of tripartite brains in the vertebrate and euarthropod lineages [1, 24, 25, 29, 30].
Hejnol A, Lowe CJ. Embracing the comparative approach: how robust phylogenies and broader developmental sampling impacts the understanding of nervous system evolution. Philos Trans R Soc Lond B Biol Sci. 2015. https://doi.org/10.1098/rstb.2015.0045.
Hirth F. On the origin and evolution of the tripartite brain. Brain Behav Evol. 2010;76(1):3–10.
Arendt D, Denes AS, Jekely G, Tessmar-Raible K. The evolution of nervous system centralization. Philos Trans R Soc Lond B Biol Sci. 2008;363(1496):1523–8.
Moroz LL. On the independent origins of complex brains and neurons. Brain Behav Evol. 2009;74(3):177–90.
Moroz L. Phylogenomics meets neuroscience: how many times might complex brains have evolved? Acta Biol Hung. 2012;63(Supplement 2):3–19.
Strausfeld NJ, Ma X, Edgecombe GD. Fossils and the evolution of the arthropod brain. Curr Biol. 2016;26(20):R989–1000.
Strausfeld NJ, Hirth F. Homology versus convergence in resolving transphyletic correspondences of brain organization. Brain Behav Evol. 2013;82(4):215–9.
Northcutt RG. Evolution of centralized nervous systems: two schools of evolutionary thought. Proc Natl Acad Sci USA. 2012;109(Suppl 1):10626–33.
Holland LZ, Carvalho JE, Escriva H, Laudet V, Schubert M, Shimeld SM, Yu J. Evolution of bilaterian central nervous systems: a single origin? EvoDevo. 2013;4(1):27.
Hirth F, Kammermeier L, Frei E, Walldorf U, Noll M, Reichert H. An urbilaterian origin of the tripartite brain: developmental genetic insights from Drosophila. Development. 2003;130(11):2365–73.
Giribet G. Genomics and the animal tree of life: conflicts and future prospects. Zool Scr. 2016;45(S1):14–21.
Sprecher SG, Reichert H. The urbilaterian brain: developmental insights into the evolutionary origin of the brain in insects and vertebrates. Arthropod Struct Dev. 2003;32(1):141–56.
Lichtneckert R, Reichert H. Insights into the urbilaterian brain: conserved genetic patterning mechanisms in insect and vertebrate brain development. Heredity. 2005;94(5):465–77.
Reichert H. A tripartite organization of the urbilaterian brain: developmental genetic evidence from Drosophila. Brain Res Bull. 2005;66(4):491–4.
Lichtneckert R, Reichert H. Anteroposterior regionalization of the brain: genetic and comparative aspects. In: Technau GM, editor. Brain development in Drosophila melanogaster. New York: Springer; 2008. p. 32–41.
Steinmetz PR, Urbach R, Posnien N, Eriksson J, Kostyuchenko RP, Brena C, Guy K, Akam M, Bucher G, Arendt D. Six3 demarcates the anterior-most developing brain region in bilaterian animals. EvoDevo. 2010;1(1):1–9.
Yang X, Weber M, ZarinKamar N, Posnien N, Friedrich F, Wigand B, Beutel R, Damen WG, Bucher G, Klingler M. Probing the Drosophila retinal determination gene network in Tribolium (II): the Pax6 genes eyeless and twin of eyeless. Dev Biol. 2009;333(1):215–27.
Blackburn DC, Conley KW, Plachetzki DC, Kempler K, Battelle B, Brown NL. Isolation and expression of Pax6 and atonal homologues in the American horseshoe crab, Limulus polyphemus. Dev Dyn. 2008;237(8):2209–19.
Kammermeier L, Leemans R, Hirth F, Flister S, Wenger U, Walldorf U, Gehring WJ, Reichert H. Differential expression and function of the Drosophila Pax6 genes eyeless and twin of eyeless in embryonic central nervous system development. Mech Dev. 2001;103(1–2):71–8.
Tomer R, Denes AS, Tessmar-Raible K, Arendt D. Profiling by image registration reveals common origin of annelid mushroom bodies and vertebrate pallium. Cell. 2010;142(5):800–9.
Strausfeld NJ, Hirth F. Deep homology of arthropod central complex and vertebrate basal ganglia. Science. 2013;340(6129):157–61.
Farries MA. How ‘basal’ are the basal ganglia? Brain Behav Evol. 2013;82(4):211–4.
Brusca RC, Brusca GJ. Invertebrates. 2nd ed. Sunderland: Sinauer Associates; 2003.
Lowe CJ, Wu M, Salic A, Evans L, Lander E, Stange-Thomann N, Gruber CE, Gerhart J, Kirschner M. Anteroposterior patterning in hemichordates and the origins of the chordate nervous system. Cell. 2003;113(7):853–65.
Pani AM, Mullarkey EE, Aronowicz J, Assimacopoulos S, Grove EA, Lowe CJ. Ancient deuterostome origins of vertebrate brain signalling centres. Nature. 2012;483(7389):289–94.
Lowe CJ, Terasaki M, Wu M, Freeman RM Jr, Runft L, Kwan K, Haigo S, Aronowicz J, Lander E, Gruber C. Dorsoventral patterning in hemichordates: insights into early chordate evolution. PLoS Biol. 2006;4(9):e291.
Fritsch M, Wollesen T, Wanninger A. Hox and ParaHox gene expression in early body plan patterning of polyplacophoran mollusks. J Exp Zool B. 2016;326(2):89–104.
Fritsch M, Wollesen T, de Oliveira A, Wanninger A. Unexpected co-linearity of Hox gene expression in an aculiferan mollusk. BMC Evol Biol. 2015;15(1):151.
Redl E, Scherholz M, Wollesen T, Todt C, Wanninger A. Expression of six3 and otx in Solenogastres (Mollusca) supports an ancestral role in bilaterian anterior-posterior axis patterning. Evol Dev. 2018;20(1):17–28.
Wollesen T, Monje SVR, 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(1):231.
Wollesen T, Scherholz M, Monje SVR, Redl E, Todt C, Wanninger A. Brain regionalization genes are co-opted into shell field patterning in Mollusca. Sci Rep. 2017;7(1):5486.
Martín-Durán JM, Pang K, Børve A, Lê HS, Furu A, Cannon JT, Jondelius U, Hejnol A. Convergent evolution of bilaterian nerve cords. Nature. 2018;553(7686):45.
Arendt D. Animal evolution: convergent nerve cords? Curr Biol. 2018;28(5):R225–7.
Campbell LI, Rota-Stabelli O, Edgecombe GD, Marchioro T, Longhorn SJ, Telford MJ, Philippe H, Rebecchi L, Peterson KJ, Pisani D. MicroRNAs and phylogenomics resolve the relationships of Tardigrada and suggest that velvet worms are the sister group of Arthropoda. Proc Natl Acad Sci USA. 2011;108(38):15920–4.
Smith MR, Ortega-Hernández J. Hallucigenia’s onychophoran-like claws and the case for Tactopoda. Nature. 2014;514(7522):363–6.
Rota-Stabelli O, Daley AC, Pisani D. Molecular timetrees reveal a Cambrian colonization of land and a new scenario for ecdysozoan evolution. Curr Biol. 2013;23(5):392–8.
Mayer G, Martin C, Rudiger J, Kauschke S, Stevenson PA, Poprawa I, Hohberg K, Schill RO, Pfluger HJ, Schlegel M. Selective neuronal staining in tardigrades and onychophorans provides insights into the evolution of segmental ganglia in panarthropods. BMC Evol Biol. 2013;13:230-2148-13-230.
Smith FW, Goldstein B. Segmentation in Tardigrada and diversification of segmental patterns in Panarthropoda. Arthropod Struct Dev. 2017;46(3):328–40.
Ortega-Hernández J, Janssen R, Budd GE. Origin and evolution of the panarthropod head—a palaeobiological and developmental perspective. Arthropod Struct Dev. 2016;46(3):354–79.
Kristensen RM. The first record of cyclomorphosis in Tardigrada based on a new genus and species from Arctic meiobenthos. J Zool Syst Evol Res. 1983;20(4):249–70.
Persson DK, Halberg KA, Jørgensen A, Møbjerg N, Kristensen RM. Brain anatomy of the marine tardigrade Actinarctus doryphorus (Arthrotardigrada). J Morphol. 2014;275(2):173–90.
Persson DK, Halberg KA, Jorgensen A, Mobjerg N, Kristensen RM. Neuroanatomy of Halobiotus crispae (Eutardigrada: Hypsibiidae): Tardigrade brain structure supports the clade Panarthropoda. J Morphol. 2012;273(11):1227–45.
Nielsen C. Animal evolution. Interrelationships of the living phyla. 2nd ed. London: Oxford University Press; 2001. p. 226–31.
Gąsiorek P, Stec D, Morek W, Michalczyk Ł. An integrative redescription of Hypsibius dujardini (Doyère, 1840), the nominal taxon for Hypsibioidea (Tardigrada: Eutardigrada). Zootaxa. 2018;4415(1):45–75.
Smith FW, Boothby TC, Giovannini I, Rebecchi L, Jockusch EL, Goldstein B. The compact body plan of tardigrades evolved by the loss of a large body region. Curr Biol. 2016;26(2):224–9.
Gabriel WN, McNuff R, Patel SK, Gregory TR, Jeck WR, Jones CD, Goldstein B. The tardigrade Hypsibius dujardini, a new model for studying the evolution of development. Dev Biol. 2007;312(2):545–59.
Smith FW, Jockusch EL. The metameric pattern of Hypsibius dujardini (Eutardigrada) and its relationship to that of other panarthropods. Front Zool. 2014;11(1):66.
Gabriel WN, Goldstein B. Segmental expression of Pax3/7 and Engrailed homologs in tardigrade development. Dev Genes Evol. 2007;217(6):421–33.
Boothby TC, Tenlen JR, Smith FW, Wang JR, Patanella KA, Nishimura EO, Tintori SC, Li Q, Jones CD, Yandell M, Messina DN, Glasscock J, Goldstein B. Evidence for extensive horizontal gene transfer from the draft genome of a tardigrade. Proc Natl Acad Sci USA. 2015;112(52):15976–81.
Delmont TO, Eren AM. Identifying contamination with advanced visualization and analysis practices: metagenomic approaches for eukaryotic genome assemblies. PeerJ. 2016;4:e1839.
Koutsovoulos G, Kumar S, Laetsch DR, Stevens L, Daub J, Conlon C, Maroon H, Thomas F, Aboobaker AA, Blaxter M. No evidence for extensive horizontal gene transfer in the genome of the tardigrade Hypsibius dujardini. Proc Natl Acad Sci USA. 2016;113(18):5053–8.
Yoshida Y, Koutsovoulos G, Laetsch DR, Stevens L, Kumar S, Horikawa DD, Ishino K, Komine S, Kunieda T, Tomita M. Comparative genomics of the tardigrades Hypsibius dujardini and Ramazzottius varieornatus. PLoS Biol. 2017;15(7):e2002266.
Levin M, Anavy L, Cole AG, Winter E, Mostov N, Khair S, Senderovich N, Kovalev E, Silver DH, Feder M. The mid-developmental transition and the evolution of animal body plans. Nature. 2016;531(7596):637–41.
Boothby TC, Tapia H, Brozena AH, Piszkiewicz S, Smith AE, Giovannini I, Rebecchi L, Pielak GJ, Koshland D, Goldstein B. Tardigrades use intrinsically disordered proteins to survive desiccation. Mol Cell. 2017;65(6):975–84.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.
Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000;17(4):540–52.
Le SQ, Gascuel O. An improved general amino acid replacement matrix. Mol Biol Evol. 2008;25(7):1307–20.
Guindon S, Dufayard J, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59(3):307–21.
Huelsenbeck JP, Ronquist F. Bayesian analysis of molecular evolution using MrBayes. In: Nielsen R, editor. Statistical methods in molecular evolution. New York: Springer; 2005. p. 183–226.
Marchler-Bauer A, Derbyshire MK, Gonzales NR, Lu S, Chitsaz F, Geer LY, Geer RC, He J, Gwadz M, Hurwitz DI. CDD: NCBI’s conserved domain database. Nucleic Acids Res. 2014;43(D1):D222–6.
Kapustin Y, Souvorov A, Tatusova T, Lipman D. Splign: algorithms for computing spliced alignments with identification of paralogs. Biol Direct. 2008;3(1):20.
Hashimoto T, Horikawa DD, Saito Y, Kuwahara H, Kozuka-Hata H, Shin T, Minakuchi Y, Ohishi K, Motoyama A, Aizu T. Extremotolerant tardigrade genome and improved radiotolerance of human cultured cells by tardigrade-unique protein. Nat Commun. 2016;7:12808.
Colombrita C, Silani V, Ratti A. ELAV proteins along evolution: back to the nucleus? Mol Cell Neurosci. 2013;56:447–55.
Smith FW, Bartels PJ, Goldstein B. A hypothesis for the composition of the tardigrade brain and its implications for panarthropod brain evolution. Integr Comp Biol. 2017;57(3):546–59.
Gross V, Mayer G. Neural development in the tardigrade Hypsibius dujardini based on anti-acetylated α-tubulin immunolabeling. EvoDevo. 2015;6(1):12.
Bebenek IG, Gates RD, Morris J, Hartenstein V, Jacobs DK. Sine oculis in basal Metazoa. Dev Genes Evol. 2004;214(7):342–51.
Ryan JF, Burton PM, Mazza ME, Kwong GK, Mullikin JC, Finnerty JR. The cnidarian-bilaterian ancestor possessed at least 56 homeoboxes: evidence from the starlet sea anemone, Nematostella vectensis. Genome Biol. 2006;7(7):R64.
Noll M. Evolution and role of Pax genes. Curr Opin Genet Dev. 1993;3(4):595–605.
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(5):255–7.
Janssen R. A molecular view of onychophoran segmentation. Arthropod Struct Dev. 2017;46(3):341–53.
Mayer G, Kauschke S, Rudiger J, Stevenson PA. Neural markers reveal a one-segmented head in tardigrades (water bears). PLoS ONE. 2013;8(3):e59090.
Schulze C, Schmidt-Rhaesa A. The architecture of the nervous system of Echiniscus testudo (Echiniscoidea, Heterotardigrada). J Limnol. 2013;72(1s):6.
Schulze C, Neves RC, Schmidt-Rhaesa A. Comparative immunohistochemical investigation on the nervous system of two species of Arthrotardigrada (Heterotardigrada, Tardigrada). Zool Anz. 2014;253(3):225–35.
Zantke J, Wolff C, Scholtz G. Three-dimensional reconstruction of the central nervous system of Macrobiotus hufelandi (Eutardigrada, Parachela): implications for the phylogenetic position of Tardigrada. Zoomorphology. 2008;127(1):21–36.
Martin C, Gross V, Pflüger H, Stevenson PA, Mayer G. Assessing segmental versus non-segmental features in the ventral nervous system of onychophorans (velvet worms). BMC Evol Biol. 2017;17(1):3.
Mayer G, Whitington PM. Neural development in Onychophora (velvet worms) suggests a step-wise evolution of segmentation in the nervous system of Panarthropoda. Dev Biol. 2009;335(1):263–75.
Mayer G, Whitington PM, Sunnucks P, Pfluger HJ. A revision of brain composition in Onychophora (velvet worms) suggests that the tritocerebrum evolved in arthropods. BMC Evol Biol. 2010;10:255-2148-10-255.
Whitington PM, Mayer G. The origins of the arthropod nervous system: insights from the Onychophora. Arthropod Struct Dev. 2011;40(3):193–209.
Martin C, Gross V, Hering L, Tepper B, Jahn H, de Sena Oliveira I, Stevenson PA, Mayer G. The nervous and visual systems of onychophorans and tardigrades: learning about arthropod evolution from their closest relatives. J Comp Physiol A. 2017;203(8):565–90.
Damen WG, Hausdorf M, Seyfarth EA, Tautz D. A conserved mode of head segmentation in arthropods revealed by the expression pattern of Hox genes in a spider. Proc Natl Acad Sci USA. 1998;95(18):10665–70.
Jager M, Murienne J, Clabaut C, Deutsch J, Le Guyader H, Manuel M. Homology of arthropod anterior appendages revealed by Hox gene expression in a sea spider. Nature. 2006;441(7092):506–8.
Telford MJ, Thomas RH. Expression of homeobox genes shows chelicerate arthropods retain their deutocerebral segment. Proc Natl Acad Sci USA. 1998;95(18):10671–5.
Eriksson BJ, Tait NN, Budd GE, Janssen R, Akam M. Head patterning and Hox gene expression in an onychophoran and its implications for the arthropod head problem. Dev Genes Evol. 2010;220(3–4):117–22.
Janssen R, Eriksson BJ, Tait NN, Budd GE. Onychophoran Hox genes and the evolution of arthropod Hox gene expression. Front Zool. 2014;11:22.
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.
Li Y, Brown SJ, Hausdorf B, Tautz D, Denell RE, Finkelstein R. Two orthodenticle-related genes in the short-germ beetle Tribolium castaneum. Dev Genes Evol. 1996;206(1):35–45.
Browne WE, Schmid BG, Wimmer EA, Martindale MQ. Expression of otd orthologs in the amphipod crustacean, Parhyale hawaiensis. Dev Genes Evol. 2006;216(10):581–95.
Janssen R, Budd GE, Damen WG. Gene expression suggests conserved mechanisms patterning the heads of insects and myriapods. Dev Biol. 2011;357(1):64–72.
Garwood RJ, Sharma PP, Dunlop JA, Giribet G. A Paleozoic stem group to mite harvestmen revealed through integration of phylogenetics and development. Curr Biol. 2014;24(9):1017–23.
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(1):3–20.
Prpic N. Duplicated Pax6 genes in Glomeris marginata (Myriapoda: Diplopoda), an arthropod with simple lateral eyes. Zoology. 2005;108(1):47–53.
Urbach R. A procephalic territory in Drosophila exhibiting similarities and dissimilarities compared to the vertebrate midbrain/hindbrain boundary region. Neural Dev. 2007;2(1):23.
Sharma PP, Schwager EE, Extavour CG, Giribet G. Hox gene expression in the harvestman Phalangium opilio reveals divergent patterning of the chelicerate opisthosoma. Evol Dev. 2012;14(5):450–63.
Janssen R, Damen WG. The ten Hox genes of the millipede Glomeris marginata. Dev Genes Evol. 2006;216(7–8):451–65.
Averof M, Patel NH. Crustacean appendage evolution associated with changes in Hox gene expression. Nature. 1997;388(6643):682.
Di-Poi N, Montoya-Burgos JI, Miller H, Pourquié O, Milinkovitch MC, Duboule D. Changes in Hox genes’ structure and function during the evolution of the squamate body plan. Nature. 2010;464(7285):99.
Rogers BT, Peterson MD, Kaufman TC. The development and evolution of insect mouthparts as revealed by the expression patterns of gnathocephalic genes. Evol Dev. 2002;4(2):96–110.
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(1):40.
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(3):109–16.
Steinmetz PR, Kostyuchenko RP, Fischer A, Arendt D. The segmental pattern of otx, gbx, and Hox genes in the annelid Platynereis dumerilii. Evol Dev. 2011;13(1):72–9.
Vopalensky P, Kozmik Z. Eye evolution: common use and independent recruitment of genetic components. Philos Trans R Soc Lond B Biol Sci. 2009;364(1531):2819–32.
Keller RG, Desplan C, Rosenberg MI. Identification and characterization of Nasonia Pax genes. Insect Mol Biol. 2010;19(s1):109–20.
Friedrich M, Caravas J. New insights from hemichordate genomes: prebilaterian origin and parallel modifications in the paired domain of the Pax gene eyegone. J Exp Zool B Mol Dev Evol. 2011;316(6):387–92.
Park TS, Kihm J, Woo J, Park C, Lee WY, Smith MP, Harper DA, Young F, Nielsen AT, Vinther J. Brain and eyes of Kerygmachela reveal protocerebral ancestry of the panarthropod head. Nat Commun. 2018;9(1):1019.
Polilov AA. At the size limit-effects of miniaturization in insects. New York: Springer; 2016. p. 240–1.
FWS and BG designed the study. FWS performed all analyses. FWS and MC performed wet lab work. FWS wrote the first draft of the manuscript. All authors read and approved the final manuscript.
We would like to thank three anonymous reviewers for helpful comments that improved the quality of our manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
For non-H. exemplaris sequences, GenBank accession numbers are provided (Additional file 1: Table S1; Additional file 2: Table S2). For H. exemplaris sequences, the source of the publically available genome or transcriptome assembly and the sequence identification numbers to access the sequences are provided (Additional file 1: Table S1; Additional file 2: Table S2).
Consent for publication
Ethics approval and consent to participate
Studies of tardigrades do not require ethics approval or consent to participate.
This study was funded by NSF Grant 1557432 to BG and startup funds provided by the University of North Florida to FWS.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.