Intact cluster and chordate-like expression of ParaHox genes in a sea star
© Annunziata et al.; licensee BioMed Central Ltd. 2013
Received: 1 February 2013
Accepted: 29 April 2013
Published: 27 June 2013
The ParaHox genes are thought to be major players in patterning the gut of several bilaterian taxa. Though this is a fundamental role that these transcription factors play, their activities are not limited to the endoderm and extend to both ectodermal and mesodermal tissues. Three genes compose the ParaHox group: Gsx, Xlox and Cdx. In some taxa (mostly chordates but to some degree also in protostomes) the three genes are arranged into a genomic cluster, in a similar fashion to what has been shown for the better-known Hox genes. Sea urchins possess the full complement of ParaHox genes but they are all dispersed throughout the genome, an arrangement that, perhaps, represented the primitive condition for all echinoderms. In order to understand the evolutionary history of this group of genes we cloned and characterized all ParaHox genes, studied their expression patterns and identified their genomic loci in a member of an earlier branching group of echinoderms, the asteroid Patiria miniata.
We identified the three ParaHox orthologs in the genome of P. miniata. While one of them, PmGsx is provided as maternal message, with no zygotic activation afterwards, the other two, PmLox and PmCdx are expressed during embryogenesis, within restricted domains of both endoderm and ectoderm. Screening of a Patiria bacterial artificial chromosome (BAC) library led to the identification of a clone containing the three genes. The transcriptional directions of PmGsx and PmLox are opposed to that of the PmCdx gene within the cluster.
The identification of P. miniata ParaHox genes has revealed the fact that these genes are clustered in the genome, in contrast to what has been reported for echinoids. Since the presence of an intact cluster, or at least a partial cluster, has been reported in chordates and polychaetes respectively, it becomes clear that within echinoderms, sea urchins have modified the original bilaterian arrangement. Moreover, the sea star ParaHox domains of expression show chordate-like features not found in the sea urchin, confirming that the dynamics of gene expression for the respective genes and their putative regulatory interactions have clearly changed over evolutionary time within the echinoid lineage.
KeywordsPatiria miniata Asteroid ParaHox Cluster Gsx Xlox Cdx Colinearity Gut
Homeobox-containing genes regulate many aspects of development in Bilateria. Based on sequence similarities, the presence of specific protein motifs, genomic arrangement and other characteristics these genes have been classified into several families. Among the best-characterized families are the Hox and the ParaHox. While the Hox family is best characterized in segmented animals (arthropods and vertebrata), where they play patterning roles in the three germ layers, little is known about the role of ParaHox in most taxa, although it is agreed that they serve an integral role in endoderm development.
A key property of Hox genes in most Bilateria is their genomic arrangement in clusters. The position of different genes within the cluster is related to their relative domains of expression along the major body axis. This property, called spatial colinearity is observed in most animals, even in those with a clearly derived morphology (that is, echinoderms ). Some taxa also exhibit a different form of colinear expression of genes, temporal colinearity, in which genes located in different positions within the cluster are activated progressively during development. While temporal colinearity seems to be strictly dependent on the presence of an intact cluster, spatial colinearity seems to be more permissive with cluster breaks. In fact, there are cases of extreme cluster disintegration, for instance in larvaceans  or acoelomorphs , where the relative spatial domains are still well conserved (with respect to the genomic locations that their clustered paralogs exhibit in other groups). This form of colinearity, in the absence of a genomic cluster, has been named trans-colinearity .
The general features exhibited by the ParaHox genes in animals are less clear. While many ParaHox genes have been identified in different animal groups, little is known about their genomic arrangements and developmental roles. A general agreement is that the ParaHox genes, Gsx, Xlox and Cdx, are an ancient sister group to the Hox genes. Both groups probably evolved from a primitive ProtoHox cluster that was duplicated, giving rise to the Hox and ParaHox clusters . Surprisingly, while the presence of a Hox cluster has been demonstrated in many animal groups, ‘intact’ ParaHox clusters are only known from Chordata. It is presently unclear whether there are complete ParaHox clusters in any protostome group (though a partially intact cluster has been described in Platynereis dumerilii: ) or in any non-chordate deuterostome. It is interesting to note that in some groups, such as sea urchins, the absence of cluster organization has not prevented the ParaHox genes from showing signs of correlative expression (both spatial and temporal ). Moreover, in that system, two of the genes (Xlox and Cdx) are clearly interlinked within a gene regulatory network that controls endoderm regionalization and exhibit mutual regulatory interactions . Whether this property is a product of an older (evolutionary) cluster association is not yet known.
The role of bilaterian ParaHox genes (with the exception of Gsx that is almost entirely expressed in the central nervous system (CNS)) seems to relate to both the CNS and endoderm patterning processes. The orthologs of the gene Gsh, the most ‘anterior’ of the ParaHox genes, are expressed in different domains of the bilaterian nervous system. A glance at the expression of the three ParaHox genes within the Bilateria shows their broad commonalities, most probably a reflection of ancestral roles. For instance, the mice Gsh1 and Gsh2 paralogs are mostly active in the developing brain [9, 10]. The related gene in amphioxus, AmphiGsx, appears first in four cells in the neural tube at the level of somite five, and later only in the cerebral vesicle ; likewise, in the ascidian Ciona intestinalis, its ortholog is expressed in the sensory vesicle . In the sea urchin embryo the gene SpGsx is detected in a small ectodermal domain, probably neural . Within the Protostomia, Gsx has been analyzed in a few phyla. In insects (where it receives the name ind) the gene is expressed in the intermediate columns of the CNS . In nereid worms (P. dumerilii is the best studied example) Gsx is expressed in two domains, one in the CNS and another in the stomadeum plus some cell clusters of the posterior foregut and midgut . Interestingly, in other polychaete worms, such as Capitella teleta, the gut domain has been lost . Dual domains of expression, in the CNS and the gut, are similarly observed in the mollusk Gibbula varia. Interestingly, in this last case it has been observed that after torsion, the domains of Gsx expression are localized to the mouth and foregut, somewhat reminiscent of the pattern seen in Platynereis. All these patterns observed in different bilaterian taxa suggest a primitive function of Gsx in patterning the CNS and the most anterior gut.
The central ParaHox gene, Xlox, has been also described in different groups of Bilateria. Strikingly, the gene seems to be lost from all insect genomes (not present in any of the genomes sequenced to date). In vertebrates, Xlox is expressed in both the developing gut [16, 17] and the CNS . In the chordate amphioxus, AmphiXlox is detected in the gut and in two cells of the neural tube . Similar domains of expression (gut plus neural ectoderm) are described for the sea urchin Xlox ortholog [7, 19]. In protostomes the patterns are, again, very similar. In the mollusk Gibbula, Xlox is expressed in the middle part of the digestive tract and in some neuroectodermal cells . Also in the polychaetes Nereis virens and P. dumerilii this gene is expressed in the midgut and in several cells of the neuroectoderm [6, 20]. In some other protostomian taxa the domains are more restricted. In the leech (Helobdella triserialis) it has been reported that Xlox is only expressed in the midgut  while in the platyhelminth Schmidtea polichroa the ortholog gene is only expressed in the nervous system . We see again that Xlox seems to be dedicated to the patterning of the gut (primarily) and some areas of the CNS.
Cdx (or caudal in insects) is the most ‘posterior’ of the ParaHox genes. It has been cloned and studied in many taxa. Three mouse paralogs (Cdx1, Cdx2 and Cdx4) are expressed in the posterior part of the gut and some areas of the CNS , domains that are shared by most vertebrates. The amphioxus ortholog, AmphiCdx, is also expressed in the posterior part of the gut and in the developing neural tube . In the sea urchin, SpCdx is active in the hindgut . Urochordates express Cdx in cells of the neural plate and the posterior gut [24, 25], although, in at least one case (in Herdmania curvata), both domains appear at different phases of their life cycle . In several arthropods, the caudal gene patterns the posterior end of the animal [26–28]. In polychaetes Cdx is expressed both in the posterior region, including the gut, in addition to areas of the CNS. The mollusk Gibbula expresses its Cdx gene in the posterior end of the digestive tract, whereas in the veliger larva, expression is localized specifically in the hindgut and rectum . Here, again, it is demonstrated that the Cdx genes are involved, as observed for Gsx and Xlox, in the patterning of structures located within the endoderm (mostly the posterior gut) and the neuroectoderm.
The commonalities shown in the expression of ParaHox genes have prompted some investigators to propose an ancestral role for this group of genes in the patterning of endoderm [29, 30] (although we know that the similarities extend to other tissues). An important issue remains to be clarified: how do these commonalities together with the taxon specificities, reflect the history of the original cluster in every clade. The problem is accentuated by the paucity of information related to the genomic organization of these genes in most phyla.
Previous analysis of a sea urchin cluster revealed that the echinoid ParaHox genes were dispersed in the genome. This finding correlates with the rearrangement that occurred within the Hox cluster in the sea urchin genome and represents a clearly derived condition for ParaHox genes. In order to understand the origin of this disruption of genomic links we decided to analyze the history of echinoderm ParaHox clusters and the roles of its genes. Here we present a detailed investigation of ParaHox gene expression and genome organization in a representative of an earlier branching group of echinoderms, the sea star P. miniata. Asteroids and echinoids diverged around 480 million years ago (Mya) . We find interesting differences in the genomic organization of the ParaHox genes in both groups. While, as previously shown, echinoids have all genes dispersed in the genome , here, in asteroids we find that all of them form a single, compact, cluster. Moreover, we find partial conservation and some remarkable divergent features in ParaHox gene expression patterns when compared with those of chordates and sea urchins.
Results and discussion
Identification of the Patiria miniata ParaHox genes: PmGsx, PmLox, and PmCdx
The case of PmLox is particularly interesting since the hexapeptide is followed by a long stretch of 23 amino acids that seem to be unique to sea star Xlox [see Additional file 1: Figure S1], since we do not find it in the published orthologs of the sea urchin Strongylocentrotus purpuratus, the mouse or the mollusk G. varia. However, a peptide of similar size is located in the same position in the cephalopod Sepia officinalis ortholog (A. Cole; personal communication), suggesting the possibility that the presence of this feature might be an old one, and perhaps has been overlooked in the Xlox genes of other animals. This sea star short peptide is encoded in an extra exon, which is not found in the equivalent position of the sea urchin genome. Since the distance between the hexapeptide and the homeodomain (the linker ) seems to be critical for the function of the protein, we suggest that the presence of this long peptide within the linker may have an important impact on its biological function (see also ). It is also possible, moreover, that via alternative splicing, the two forms, with and without the extra exon, are both produced in the embryo (a fact that has not been investigated so far).
PmGsx, PmLox and PmCdx are in cluster in the sea star Patiria miniatagenome
Sea star ParaHox genes are located in a continuous region of 87 kb (depicted in Figure 1B), without any intervening ORFs in between, as determined by sequencing of clones from a P. miniata BAC library. The intergenic distances are 31 kb between PmGsx and PmLox and 13 kb between PmLox and PmCdx. As in other deuterostome clusters, the transcriptional orientation of PmGsx and PmLox are the same while that for PmCdx is reversed. As shown in Figure 1B, PmGsx is composed of two exons while PmLox and PmCdx contain three. In the case of the Cdx gene, the homeobox is split by one intron (see Figure 1A and B). In PmLox and PmGsx, the homeobox is located integrally in the last exon (see asterisks in Figure 1B). All exon (protein-coding) and intron sizes are listed in the Additional file 1: Table S1.
Our results demonstrate that Gsx, Xlox and Cdx are clustered in the sea star genome, a situation radically different from what we have seen in the sea urchins . How does this arrangement relate to the situations in other bilaterian taxa? While ParaHox genes are present in all bilaterian phyla, their cluster arrangement has been demonstrated for a limited number of taxa. Intact clusters appear only in some deuterostomes (Hemichordata and Chordata; [5, 36]). However, a non-intact ParaHox cluster is clearly present in the annelid P. dumerilii. In this latter case, the orthologs of the Gsx, Xlox and Cdx are located within the same chromosome but while the Pdu-Gsx and Pdu-Xlox genes are neighbors (46 kb apart), the Pdu-Cdx gene is located at the opposite end of the chromosome. Interestingly, Hui and collaborators  have shown that genes located in the neighborhood of the Platynereis ParaHox genes can also be found in close proximity to the ParaHox genes in humans, indicating a unique origin of the ParaHox clusters within the Bilateria (see also ). In fact, a thorough analysis of syntenic regions around the ParaHox cluster suggests that the ParaHox cluster was already present in the ancestor of the bilaterians plus the cnidarians . Nonetheless, the ParaHox cluster has been broken repeatedly in different lineages, perhaps due to the activity of transposable elements  or other sources of genomic instability (see, for instance ). Within the Protostomes we find that, for instance, while Platynereis has a (partially) intact cluster, other polychaetes such as Capitella teleta have lost any signs of clustering among their ParaHox genes. It is possible that the ancestral condition within the Lophotrochozoa was that of a pair of tightly linked Gsx and Xlox genes and a dispersed Cdx. This is supported by the fact that the genome of the limpet Lottia gigantea also shows this ParaHox arrangement (DOE Joint Genome Institute ). Within the Ecdysozoa the situation is more simplified, with genes lost in many lineages (Xlox in all of them) and the genomic dispersion of Gsx and Cdx orthologs.
The Deuterostomia, a monophyletic group, comprises four phyla (Hemichordata, Echinodermata, Urochordata and Chordata), with the Xenoacoelomorpha ostensibly as its putative fifth phylum . While, as mentioned, clusters of ParaHox genes have been detected in all these phyla (with the exception of Xenoacelomorpha), there are multiple simplifications occurring in the respective lineages. Intact clusters are seen in vertebrates such as humans, mouse and Xenopus (a maximum of one intact cluster remains in each studied species) but in teleost fishes the ParaHox cluster was apparently lost . Interestingly, in the bichir Polypterus senegalus, the most basal extant ray-finned fish, the cluster is intact, as it is in the immediate outgroup of the teleost fishes (for instance, the bowfin Amia calva; ). In all these cases, the gene organization and composition of ParaHox clusters resemble the condition described originally in the amphioxus genome . The transcriptional orientations of the Gsx and Xlox paralogs are the same, and opposite to that of Cdx. A case of disintegration in the cluster is observed in the urochordate Ciona intestinalis where the three genes are dispersed in two chromosomes, Gsx on chromosome 2q and Xlox and Cdx on chromosome 14q. Xlox and Cdx are separated, however, by 240 kb with many intervening genes in between, a sign of progressive dispersion . Hemichordates have kept their cluster intact, at least in the Ptychoderidae (reported in [36, 45]).
Temporal expression pattern of the three ParaHox genes during sea star embryonic development
Clearly, during the time frame of these experiments, which spans embryogenesis in its entirety, PmLox has a single peak of maximum accumulation at day three, while PmCdx seems to accumulate in two different waves, with peaks of accumulation at days one and five, respectively. The levels of PmGsx mRNA are always low, due most probably to solely maternal contribution. Based on the QPCR data, and understanding that PmGsx is only maternal, although detectable until 24 hpf, we would suggest that PmCdx is the first ParaHox gene to be activated during sea star development, followed by PmLox, which is turned on more than one day later. Moreover, a possible role for PmGsx in advanced larval stages or in adult body formation would indicate that Gsx is the last ParaHox gene to be activated. Our data show that the sea star ParaHox transcription factor temporal activation resembles the sequence of activation observed in chordates  while it is inverted for the sea urchin orthologs . Moreover, the sea star ParaHox cluster characteristics support the theory that mechanisms producing temporal colinearity are likely the major constraining forces on gene cluster maintenance [46–48]. In the sea star, both the temporal order of expression (at least for Xlox and Cdx) and the cluster organization, typical of Chordata, seem to be conserved, while in the sea urchin, the temporal colinearity (assuming the cluster arrangement in chordates) is inverted and the cluster is broken.
The spatial expression domains of the ParaHox genes during sea star embryonic development
Considering all our results together, some very interesting features can be highlighted. First of all, conservation in the expression domains of the sea urchin and sea star orthologs is of note. Further, the expression of the sea star Cdx in the blastula stage, well before gastrulation has been initiated, implicates a function for Cdx that is not present in sea urchins. However, an early activation of Cdx in development, prior to its recruitment in posterior gut patterning, appears as a common feature among chordates: in murine embryos the first embryonic territory of expression for Cdx genes is the posterior primitive streak ; moreover, in Xenopus tropicalis the three orthologs Cad1, Cad2 and Cad3 are first expressed in the early gastrula around the blastopore and later in the posterior embryo, including the gut ; finally, the first expression of the amphioxus Cdx orthologue is detected in a ring of cells surrounding the blastopore opening . The conservation of Cdx gene expression around the primitive streak of the early mouse and the blastopore of frog, amphioxus and sea star embryos, suggests the existence of an ancestral role for Cdx in the early stages of deuterostome development. This function has been lost in the sea urchin embryo where only a late role for Cdx in patterning posterior structures, in this case the gut, has been found .
PmLox and PmCdxdynamic relative expression along the sea star developing gut
PmLoxectodermal domain of expression
Xlox expression in neural territories has been reported in many animals, protostomes and deuterostomes alike [6, 11, 18, 20, 57]. Within the echinoderms it has been shown that in the sea urchin Strongylocentrotus purpuratus SpLox is specifically expressed in cells localized below the ciliary band, in a position corresponding to the lateral ganglion .
Evolutionary implications: changing the genomic organization and the activities of ParaHox genes within the Echinodermata
Over the last few years it has been shown that ParaHox genes are arranged in genomic clusters in several bilaterian taxa, a fact that might have important consequences for their regulation. However, this sometimes tight arrangement is broken in many lineages. One clear example of this evolutionary phenomenon is shown in the echinoderms. We have shown in the past that the three ParaHox genes are not clustered in the sea urchin genome . Here, instead, we demonstrate that a member of the asteroids (P. miniata) has kept the ancestral condition, the presence of a tight cluster of genes , a condition retained in the sister group of all echinoderms, the hemichordates . The evolutionary history of these genes and their putative ancestral roles merit some further comments, which are given in the next paragraphs.
When evaluating the dynamics of some genomic arrangements, we think it is interesting to note that the partial break of ParaHox clusters, or their eventual disintegration, seems to occur in parallel, in the various taxa where these events take place, to the break and dispersion of its sister cluster, the HOX. The modification of the cluster seems to be related to the average substitution rate of the different taxa, as seen in some metazoan phylograms , such that faster evolving clades are more prone to cluster modification. An example in this context would be the sea urchins (echinoids) where the Hox cluster has been reorganized through breaks and the loss of genes (Hox4) or the fusion of transcriptional units (Hox5). In parallel, a dispersion of the three ParaHox genes within the genome has occurred. Interestingly, one break point in the echinoid Hox cluster occurs at the Hox4 locus, with the consequence that this gene is lost . In sea stars, Hox 4 is retained  suggesting that, perhaps, this break is not present in asteroids. If this were the case, and taking into account what we now observe in the Patiria ParaHox, we would be facing another case of parallel Hox/ParaHox clusters evolutionary histories. Hox and ParaHox clusters would be retained in asteroids but broken (to different degrees) in echinoids.
Given the presence of an intact cluster in deuterostomes and of a partially intact cluster in protostomes, we can infer that the common ancestor of all Bilateria had a ParaHox cluster composed of three genes, Gsx, Xlox and Cdx. This fact is supported by the shared presence of similar genes in the genomic neighborhood of Platynereis and human ParaHox genes . The cluster itself might have originated from a single precursor gene that was duplicated in cis, or, alternatively, it was derived, as a whole, from the duplication of a ProtoHox ancestral cluster . This tight genome arrangement was kept in several lineages although in others, perhaps as a result of a relaxed selection, it was broken into pieces. Within the echinoderms, we have a clear example of the dynamics of the cluster in evolutionary time. We show that asteroids (sea stars) maintain a tight chromosomal arrangement of the three genes while in the echinoids (sea urchins) the genes are dispersed in the genome (this paper and ). Interestingly, and independently of the genomic arrangements, there are some aspects of ParaHox gene regulation that are kept, even in the absence of an intact cluster. The most important of these aspects is the nested nature of their endodermal expression domains. When it comes to spatial domains within the gut, in almost all cases studied the expression patterns of the genes (but fundamentally Xlox and Cdx) are maintained, with Xlox being expressed anteriorly to Cdx. The presence of Gsx in the bilaterian guts is less conserved. It is very possible, as it has been shown for the Hox genes, that a tight cluster arrangement is a prerequisite for a strict temporal control of the genes expression, while the relative spatial domains can be preserved in the absence of tight clustering. This suggested hypothesis has not been tested for the ParaHox genes, since most of the observations are purely circumstantial at this point. No clear kinetic studies allow us to confirm or reject such an assumption at the moment.
Given the use of ParaHox genes in gut patterning, and their historical birth in the cis-duplication of a precursor gene, it is reasonable to assume a scenario in which the bilaterian ancestor used a single proto-parahox gene (one possible candidate being the Placozoa Trox-2 gene ) to control the specification of the gut tissue (although not excluding its use in other germ layers). This precursor gene was duplicated in cis, giving rise to the three extant ParaHox genes, Gsx, Xlox and Cdx, linked in the genome. These genes evolved their functions by sub-functionalization  and, perhaps in parallel, giving rise to different functional domains within the gut. In some cases the domains were maintained by cross-regulatory interactions (sea urchins, for instance), although this could also be maintained even in the event of cluster disintegration. The cluster in itself proved quite flexible, in terms of size, as can be seen clearly in Figure 6. Strikingly, the changes in gene spacing within the ParaHox cluster occurred without the incorporation of new ORFs inside it. However, the few analyses we have on the occurrence and nature of the intervening sequences have demonstrated the presence of many repetitive elements, or rests of transposable elements in these regions, suggesting that one mechanism controlling the cluster size might be through the expansion of repetitive sequences or the inclusion (transposition) of new mobile elements . At the same time these elements could provide the structural basis for the later cluster disintegration in specific lineages (as seen in the Hox cluster; ).
We should emphasize, once more, that it is becoming clear, from the analysis of different animal genomes, that there is a parallel disintegration of Hox and ParaHox clusters in some lineages. This might not be so surprising being most probably the reflection of the evolutionary rates of change in the DNAs of each lineage (seen, for instance, as branch lengths in phylogenetic analysis). The cases of tunicates, or acoels, are paradigmatic, but it is a fact also demonstrated here for the echinoderms, where a disintegration of the ParaHox cluster has accompanied the parallel rupture of the Hox clusters. This leads us to the inescapable conclusion (not surprising, though) that the dynamics of cluster evolution are a consequence of the global dynamics followed by the genome over evolutionary time.
Recent advances in genome sequencing technologies have produced a flurry of genomic data from many, previously inaccessible, animal systems. This has generated a renewed interest in the dynamics of genome evolution and in the relationship between changes in the genome and different morphological transformations.
Among those features that are now being thoroughly investigated, the relationship between gene clustering and regulation is of special relevance. Hox and ParaHox genes have been classical examples of vectorial patterning systems with complex regulation. However, most of the data supporting current models of Hox and ParaHox function have been obtained in a relatively small group of animals, mainly insects and chordates. It is clear that a full understanding of how these gene families (and their functions) have evolved over time will be possible only when a wider selection of animal models is considered.
Here, we present a thorough investigation of one such group of genes, the ParaHox. We have cloned and thoroughly characterized the ParaHox group of genes in the sea star P. miniata. Moreover, we have studied their genomic organization and the expression domains of each gene, in space and time, during embryogenesis. We show that sea stars organize the ParaHox genes in a single, compact cluster, reflecting their basal position within the Echinodermata and highlighting the more derived nature of the related group of echinoids. We show that both gene cluster disruption and the changing patterns of gene expression are happening within one single phylum, a clear sign of the dynamic nature of genomes during evolution.
Cloning of PmGsx, PmLox and PmCdx
Universal degenerate primers amplifying the ANTP family homeoboxes were used on cDNA synthesized from a mixture of 3- and 4-day old P. miniata embryo RNA, following the protocol published in Martinez et al. . The PCR product mix was cloned into the pGEM-T Easy vector (Promega, Madison, WI, USA) and a large number of colonies were screened leading to the isolation of the clones containing PmLox and PmCdx homeoboxes. Using as probe the PmLox homeobox, a clone containing a 1,454 bp PmLox fragment was successfully isolated from a P. miniata (3-day embryo) cDNA library. A 3′-RACE strategy was followed to obtain a longer portion of PmCdx. In particular, a mix of 3- and 4-day RNA was used as template for the 3′-RACE performed using the kit 3′-RACE System for Rapid Amplification of cDNA Ends (Invitrogen, Carlsberg, CA, USA); two primers, a first forward primer and a second nested primer were designed based on the available sequence (PmCdx3′-RACE: ACATCACCATCAGACGCAAG, PmCdx3′n-RACE: GGGACTATCCGAGAGACAGG) and the PCRs were conducted following the manufacturer’s instructions. The amplified product was cloned into Topo-TA cloning vector (Invitrogen). A fragment of 1,221 bp of PmGsx was cloned from P. miniata genomic DNA designing specific primers on the sequence of the BAC containing the three P. miniata ParaHox genes (PmGsxF1: AAAACACCGAAAATTGCAAAG, PmGsxR1: AGTTTTGCGGCCACTTTCTA).
cDNA and BAC filter screenings
P. miniata cDNA library arrayed filters were screened using as probes PmLox and PmCdx homeoboxes, while the BAC library filters were screened using the PmLox cDNA clone and PmCdx fragment obtained by 3′-RACE. The screenings were performed following the protocol described in Martinez et al. .
Sea star embryos, in situhybridization and immunostaining experiments
Adult P. miniata were obtained from Patrick Leahy (Kerchoff Marine Laboratory, California Institute of Technology, Pasadena, CA, USA) and housed in circulating sea water aquaria in the Stazione Zoologica Anton Dohrn of Naples. Gametes were obtained following the procedure described in Hinman et al.  and embryos were raised at 15°C in filtered sea water diluted 9:1 with de-ionized water. After day 4, larvae were fed daily with a mix of Isochrysis galbana and Rhodomonas lens algae. For in situ hybridization, probes were transcribed from purified PCR amplified template DNA, using digoxygenin-11-UTP (Roche, Indianapolis, IN, USA) or labeled with dinitrophenol (DNP) (Mirus, Madison, WI, USA) following kit instructions. Whole mount in situ hybridization experiments were performed as described in Hinman et al. , with the modifications suggested in Yankura et al.  for the double fluorescent in situ hybridization procedure. For the acetylated tubulin immunostaining, embryos and larvae were fixed overnight at 4°C in PBS with 2% paraformaldehyde (PFA), washed in phosphate buffered saline with Tween (PBST) several times, incubated in the blocking solution (5% sheep serum, 1X PBST) for one hour at room temperature and then incubated overnight in the blocking solution with a dilution 1:250 of the mouse monoclonal anti acetylated tubulin antibody (Sigma-Aldrich, St Louis, MO, USA). The next day the embryos were rinsed several times with PBST, incubated in the secondary antibody, the AlexaFluor 488 goat anti- mouse immunoglobulin G (IgG) (Molecular Probes, Eugene, OR, USA), and after one hour washed to remove the unbound antibody. For the fluorescent in situ hybridization and the immunostaining, embryos were examined and images obtained with a Zeiss confocal laser-scanning microscope LSM 510. For the colorimetric in situ hybridization, embryos were observed and pictures made with the use of a Zeiss digital camera (Axiocam) mounted on a Zeiss Axioimage 2 MOT microscope operating in DIC mode.
Quantitative polymerase chain reaction
Total RNA was collected from 300 embryos for each analyzed stage with the use of the Ambion® RNAqueous-Micro Kit (Life Technologies, Carlsberg, CA, USA) and cDNA synthesized with the SuperScript® VILO™ cDNA Synthesis Kit (Invitrogen) following the manufacturer’s instructions. Temporal accumulation of messages was monitored using real-time quantitative PCR (QPCR). Specific primer sets for each gene were designed on separate exons (for PmLox and PmGsx two sets of primers were designed and used in the QPCR experiments). Reactions were performed using the ViiA 7 REAL TIME PCR with SYBR Green chemistry (Applied Biosystems, Foster City, CA, USA). Ubiquitin was used to normalize gene expression data following Hinman et al. . Sequence of the primers used for QPCR: PmLoxQF1: GCCGCATCATCATCATACAC; PmLoxQR1: ACATATGAGCGTGCGATTTG; PmLoxQF2: CACGCTCATATGTGGAAAGC; PmLoxQR2: CGCTTGTTCTCGTCAAAGTCT; PmCdxQF: ACCGGAGATGGTCCTGAAC; PmCdxQR: GCATGAAGACAGGGCAGTTT; PmGsxQF1: AGACCCGAGGAGACTCCAAT; PmGsxQR1: ACATGTTGGAGGCGAACTCT; PmGsxQF2: CGGCGGGTCAAGTACAAG; PmGsxQR2: CAATCTCTCTCCGGCTGC).
BAC sequencing and assembling
A mixture of the BAC clones, CeCl-purified, positive for both PmLox and PmCdx, was sequenced at the European Molecular Biology Laboratory (EMBL) Genomics Core Facility using a massive parallel sequencing platform (Illumina, San Diego, CA, USA). Sequence gaps were filled using specific primers and targeted capillary electrophoresis (Abi3730DNA analyzer, Life Technologies, Carlsberg, CA, USA) at the Molecular Biology Service of the Stazione Zoologica Anton Dohrn of Naples. The assembly of the individual reads was performed by the EMBL bioinformatics service, and further assembly including the sequences obtained through the walking strategy was performed using the SeqMan Pro Lasergene suite program. A prediction of the exon-intron junctions was performed using Genescan software and validated with the available cDNA sequences. The sequences of the complete coding region of the three genes have been deposited in Genbank [KC551919, KC551920, KC551921] and are reported in Additional file 1.
Bacterial artificial chromosome
Central nervous system
Days post fertilization
Hours post fertilization
Million years ago
Open reading frame
Phosphate buffered saline with Tween-20
Polymerase chain reaction
Post oral ciliary band
- 3' RACE:
Rapid amplification of cDNA ends.
The regular research visits of PM to the Stazione Zoologica have been funded by the European Union program ‘ASSEMBLE’. We thank Veronica Hinman and Brenna McCauley for clones and assistance with Patiria miniata manipulation; Andrew Cameron for BAC clones; Salvatore D’Aniello for support with phylogenetic analysis; Alessandro Amoroso and Davide Caramiello for technical assistance; the SZN Molecular Biology Service and the EMBL GeneCore facility for BAC sequencing; Patrick Leahy and the Southern California Sea Urchin Company for animal supply; Chris Amemiya and Alison Cole for kindly revising our manuscript. This work was supported in part by the ‘Molmorph’ (contract number MEST-CT-2005–020542, fellowship to RA) and the ‘Evonet’ (project 215781) EU-Marie Curie Early Training Network. The reviewers provided thorough and constructive criticisms helping us to improve the manuscript.
- Arenas-Mena C, Cameron AR, Davidson EH: Spatial expression of Hox cluster genes in the ontogeny of a sea urchin. Development. 2000, 127: 4631-4643.PubMed
- Seo HC, Edvardsen RB, Maeland AD, Bjordal M, Jensen MF, Hansen A, Flaat M, Weissenbach J, Lehrach H, Wincker P, Reinhardt R, Chourrout D: Hox cluster disintegration with persistent anteroposterior order of expression in Oikopleura dioica. Nature. 2004, 431: 67-71. 10.1038/nature02709.View ArticlePubMed
- Moreno E, Nadal M, Baguñà J, Martínez P: Tracking the origins of the bilaterian Hox patterning system: insights from the acoel flatworm Symsagittifera roscoffensis. Evol Dev. 2009, 11: 574-581. 10.1111/j.1525-142X.2009.00363.x.View ArticlePubMed
- Duboule D: The rise and fall of Hox gene clusters. Development. 2007, 134: 2549-2560. 10.1242/dev.001065.View ArticlePubMed
- Brooke NM, Garcia-Fernandez J, Holland PW: The ParaHox gene cluster is an evolutionary sister of the Hox gene cluster. Nature. 1998, 392: 920-922. 10.1038/31933.View ArticlePubMed
- Hui JH, Raible F, Korchagina N, Dray N, Samain S, Magdelenat G, Jubin C, Segurens B, Balavoine G, Arendt D, Ferrier DE: Features of the ancestral bilaterian inferred from Platynereis dumerilii ParaHox genes. BMC Biol. 2009, 7: 43-10.1186/1741-7007-7-43.PubMed CentralView ArticlePubMed
- Arnone MI, Rizzo F, Annunciata R, Cameron RA, Peterson KJ, Martínez P: Genetic organization and embryonic expression of the ParaHox genes in the sea urchin S. purpuratus: insights into the relationship between clustering and colinearity. Dev Biol. 2006, 300: 63-73. 10.1016/j.ydbio.2006.07.037.View ArticlePubMed
- Cole AG, Rizzo F, Martinez P, Fernandez-Serra M, Arnone MI: Two ParaHox genes, SpLox and SpCdx, interact to partition the posterior endoderm in the formation of a functional gut. Development. 2009, 136: 541-549. 10.1242/dev.029959.View ArticlePubMed
- Valerius MT, Li H, Stock JL, Weinstein M, Kaur S, Singh G, Potter SS: Gsh-1: a novel murine homeobox gene expressed in the central nervous system. Dev Dyn. 1995, 203: 337-351. 10.1002/aja.1002030306.View ArticlePubMed
- Hsieh-Li HM, Witte DP, Szucsik JC, Weinstein M, Li H, Potter SS: Gsh-2, a murine homeobox gene expressed in the developing brain. Mech Dev. 1995, 50: 177-186. 10.1016/0925-4773(94)00334-J.View ArticlePubMed
- Osborne PW, Benoit G, Laudet V, Schubert M, Ferrier DE: Differential regulation of ParaHox genes by retinoic acid in the invertebrate chordate amphioxus (Branchiostoma floridae). Dev Biol. 2009, 327: 252-262. 10.1016/j.ydbio.2008.11.027.View ArticlePubMed
- Hudson C, Lemaire P: Induction of anterior neural fates in the ascidian Ciona intestinalis. Mech Dev. 2001, 100: 189-203. 10.1016/S0925-4773(00)00528-1.View ArticlePubMed
- Weiss JB, Von Ohlen T, Mellerick DM, Dressler G, Doe CQ, Scott MP: Dorsoventral patterning in the Drosophila central nervous system: the intermediate neuroblasts defective homeobox gene specifies intermediate column identity. Genes Dev. 1998, 12: 3591-3602. 10.1101/gad.12.22.3591.PubMed CentralView ArticlePubMed
- Frobius AC, Seaver EC: ParaHox gene expression in the polychaete annelid Capitella sp I. Dev Genes Evol. 2006, 216: 81-88. 10.1007/s00427-005-0049-0.View ArticlePubMed
- Samadi L, Steiner G: Conservation of ParaHox genes’ function in patterning of the digestive tract of the marine gastropod Gibbula varia. BMC Dev Biol. 2010, 10: 74-10.1186/1471-213X-10-74.PubMed CentralView ArticlePubMed
- Wright CV, Cho KW, Oliver G, De Robertis EM: Vertebrate homeodomain proteins: families of region-specific transcription factors. Trends Biochem Sci. 1989, 14: 52-56. 10.1016/0968-0004(89)90043-1.View ArticlePubMed
- Jonsson J, Carlsson L, Edlund T, Edlund H: Insulin-promoter-factor 1 is required for pancreas development in mice. Nature. 1994, 371: 606-609. 10.1038/371606a0.View ArticlePubMed
- Perez-Villamil B, Schwartz PT, Vallejo M: The pancreatic homeodomain transcription factor IDX1/IPF1 is expressed in neural cells during brain development. Endocrinology. 1999, 140: 3857-3860. 10.1210/en.140.8.3857.View ArticlePubMed
- Cole AG, Arnone MI: Fluorescent in situ hybridization reveals multiple expression domains for SpBrn1/2/4 and identifies a unique ectodermal cell type that co-expresses the ParaHox gene SpLox. Gene Expr Patterns. 2009, 9: 324-328. 10.1016/j.gep.2009.02.005.View ArticlePubMed
- Kulakova MA, Cook CE, Andreeva TF: ParaHox gene expression in larval and postlarval development of the polychaete Nereis virens (Annelida, Lophotrochozoa). BMC Dev Biol. 2008, 8: 61-10.1186/1471-213X-8-61.PubMed CentralView ArticlePubMed
- Wedeen CJ, Shankland M: Mesoderm is required for the formation of a segmented endodermal cell layer in the leech Helobdella. Dev Biol. 1997, 191: 202-214. 10.1006/dbio.1997.8735.View ArticlePubMed
- Martin-Duran JM, Romero R: Evolutionary implications of morphogenesis and molecular patterning of the blind gut in the planarian Schmidtea polychroa. Dev Biol. 2011, 352: 164-176. 10.1016/j.ydbio.2011.01.032.View ArticlePubMed
- Young T, Deschamps J: Hox, Cdx, and anteroposterior patterning in the mouse embryo. Curr Top Dev Biol. 2009, 88: 235-255.View ArticlePubMed
- Hinman VF, Becker E, Degnan BM: Neuroectodermal and endodermal expression of the ascidian Cdx gene is separated by metamorphosis. Dev Genes Evol. 2000, 210: 212-216. 10.1007/s004270050306.View ArticlePubMed
- Mita K, Fujiwara S: Nodal regulates neural tube formation in the Ciona intestinalis embryo. Dev Genes Evol. 2007, 217: 593-601. 10.1007/s00427-007-0168-x.View ArticlePubMed
- Mlodzik M, Fjose A, Gehring WJ: Isolation of caudal, a Drosophila homeo box-containing gene with maternal expression, whose transcripts form a concentration gradient at the pre-blastoderm stage. EMBO J. 1985, 4: 2961-2969.PubMed CentralPubMed
- Copf T, Rabet N, Celniker SE, Averof M: Posterior patterning genes and the identification of a unique body region in the brine shrimp Artemia franciscana. Development. 2003, 130: 5915-5927. 10.1242/dev.00835.View ArticlePubMed
- Schulz C, Schröder R, Hausdorf B, Wolff C, Tautz D: A caudal homologue in the short germ band beetle Tribolium shows similarities to both, the Drosophila and the vertebrate caudal expression patterns. Dev Genes Evol. 1998, 208: 283-289. 10.1007/s004270050183.View ArticlePubMed
- Holland PW: Beyond the Hox: how widespread is homeobox gene clustering?. J Anat. 2001, 199: 13-23. 10.1046/j.1469-7580.199.parts1-2.3.x.PubMed CentralView ArticlePubMed
- Holland PW: Evolution of homeobox genes. WIREs Dev Biol. 2013, 2: 31-45. 10.1002/wdev.78.View Article
- Pisani D, Feuda R, Peterson KJ, Smith AB: Resolving phylogenetic signal from noise when divergence is rapid: a new look at the old problem of echinoderm class relationships. Mol Phylogenet Evol. 2012, 62: 27-34. 10.1016/j.ympev.2011.08.028.View ArticlePubMed
- Martinez P, Lee JC, Davidson EH: Complete sequence of SpHox8 and its linkage in the single Hox gene cluster of Strongylocentrotus purpuratus. J Mol Evol. 1997, 44: 371-377. 10.1007/PL00006156.View ArticlePubMed
- Ferrier DE, Dewar K, Cook A, Chang JL, Hill-Force A, Amemiya C: The chordate ParaHox cluster. Curr Biol. 2005, 15: R820-R822. 10.1016/j.cub.2005.10.014.View ArticlePubMed
- In der Rieden PM, Mainguy G, Woltering JM, Durston AJ: Homeodomain to hexapeptide or PBC-interaction-domain distance: size apparently matters. Trends Genet. 2004, 20: 76-79. 10.1016/j.tig.2003.12.001.View ArticlePubMed
- Gebelein B, Culi J, Ryoo HD, Zhang W, Mann RS: Specificity of Distalless repression and limb primordia development by abdominal Hox proteins. Dev Cell. 2002, 3: 487-498. 10.1016/S1534-5807(02)00257-5.View ArticlePubMed
- Ikuta T: Evolution of invertebrate deuterostomes and Hox/ParaHox genes. Genomics Proteomics Bioinformatics. 2011, 9: 77-96. 10.1016/S1672-0229(11)60011-9.View ArticlePubMed
- Mendivil Ramos O, Barker D, Ferrier DE: Ghost loci imply Hox and ParaHox existence in the last common ancestor of animals. Curr Biol. 2012, 22: 1951-1956. 10.1016/j.cub.2012.08.023.View ArticlePubMed
- Hui JH, Holland PW, Ferrier DE: Do cnidarians have a ParaHox cluster? Analysis of synteny around a Nematostella homeobox gene cluster. Evol Dev. 2008, 10: 725-730. 10.1111/j.1525-142X.2008.00286.x.View ArticlePubMed
- Osborne PW, Ferrier DE: Chordate Hox and ParaHox gene clusters differ dramatically in their repetitive element content. Mol Biol Evol. 2010, 27: 217-220. 10.1093/molbev/msp235.View ArticlePubMed
- Lupski JR, Stankiewicz P: Genomic disorders: molecular mechanisms for rearrangements and conveyed phenotypes. PLoS Genet. 2005, 1: e49-10.1371/journal.pgen.0010049.PubMed CentralView ArticlePubMed
- JGI Genome Portal: Lottia gigantea v1.0. [http://genome.jgi-psf.org/Lotgi1/Lotgi1.home.html]
- Philippe H, Brinkmann H, Copley RR, Moroz LL, Nakano H, Poustka AJ, Wallberg A, Peterson KJ, Telford MJ: Acoelomorph flatworms are deuterostomes related to Xenoturbella. Nature. 2011, 470: 255-258. 10.1038/nature09676.PubMed CentralView ArticlePubMed
- Mulley JF, Chiu CH, Holland PW: Breakup of a homeobox cluster after genome duplication in teleosts. Proc Natl Acad Sci U S A. 2006, 103: 10369-10372. 10.1073/pnas.0600341103.PubMed CentralView ArticlePubMed
- Shoguchi E, Kawashima T, Satou Y, Hamaguchi M, Sin-I T, Kohara Y, Putnam N, Rokhsar DS, Satoh N: Chromosomal mapping of 170 BAC clones in the ascidian Ciona intestinalis. Genome Res. 2006, 16: 297-303.PubMed CentralView ArticlePubMed
- Ikuta T, Chen Y-C, Annunziata R, Ting H-C, Tung C-H, Koyanagi R, Tagawa K, Humphreys T, Fujiyama A, Saiga H, Satoh N, Yu J-K, Arnone MI, Su Y-H: Identification of an intact ParaHox cluster with temporal colinearity but residual spatial colinearity in the hemichordate Ptychodera flava. BMC Evol Biol. 2013, 13: 129-10.1186/1471-2148-13-129.PubMed CentralView ArticlePubMed
- Ferrier DE, Holland PW: Ciona intestinalis ParaHox genes: evolution of Hox/ParaHox cluster integrity, developmental mode, and temporal colinearity. Mol Phylogenet Evol. 2002, 24: 412-417. 10.1016/S1055-7903(02)00204-X.View ArticlePubMed
- Ferrier DE, Minguillon C: Evolution of the Hox/ParaHox gene clusters. Int J Dev Biol. 2003, 47: 605-611.PubMed
- Tschopp P, Tarchini B, Spitz F, Zakany J, Duboule D: Uncoupling time and space in the collinear regulation of Hox genes. PLoS Genet. 2009, 5: e1000398-10.1371/journal.pgen.1000398.PubMed CentralView ArticlePubMed
- Kuraishi R, Osanai K: Cell movements during gastrulation of starfish larvae. Biol Bull. 1992, 183: 258-268. 10.2307/1542213.View Article
- Chawengsaksophak K, de Graaff W, Rossant J, Deschamps J, Beck F: Cdx2 is essential for axial elongation in mouse development. Proc Natl Acad Sci USA. 2004, 101: 7641-7645. 10.1073/pnas.0401654101.PubMed CentralView ArticlePubMed
- Reece-Hoyes JS, Keenan ID, Isaacs HV: Cloning and expression of the Cdx family from the frog Xenopus tropicalis. Dev Dyn. 2002, 223: 134-140. 10.1002/dvdy.1234.View ArticlePubMed
- Gamer LW, Wright CV: Murine Cdx-4 bears striking similarities to the Drosophila caudal gene in its homeodomain sequence and early expression pattern. Mech Dev. 1993, 43: 71-81. 10.1016/0925-4773(93)90024-R.View ArticlePubMed
- Marom K, Shapira E, Fainsod A: The chicken caudal genes establish an anterior-posterior gradient by partially overlapping temporal and spatial patterns of expression. Mech Dev. 1997, 64: 41-52. 10.1016/S0925-4773(97)00043-9.View ArticlePubMed
- Morales AV, de la Rosa EJ, de Pablo F: Expression of the cCdx-B homeobox gene in chick embryo suggests its participation in rostrocaudal axial patterning. Dev Dyn. 1996, 206: 343-353. 10.1002/(SICI)1097-0177(199608)206:4<343::AID-AJA1>3.0.CO;2-I.View ArticlePubMed
- Meyer BI, Gruss P: Mouse Cdx-1 expression during gastrulation. Development. 1993, 117: 191-203.PubMed
- Grainger S, Savory JG, Lohnes D: Cdx2 regulates patterning of the intestinal epithelium. Dev Biol. 2010, 339: 155-165. 10.1016/j.ydbio.2009.12.025.View ArticlePubMed
- Corrado M, Aniello F, Fucci L, Branno M: Ci-IPF1, the pancreatic homeodomain transcription factor, is expressed in neural cells of Ciona intestinalis larva. Mech Dev. 2001, 102: 271-274. 10.1016/S0925-4773(01)00311-2.View ArticlePubMed
- Nakajima Y, Kaneko H, Murray G, Burke RD: Divergent patterns of neural development in larval echinoids and asteroids. Evol Dev. 2004, 6: 95-104. 10.1111/j.1525-142X.2004.04011.x.View ArticlePubMed
- Yankura KA, Martik ML, Jennings CK, Hinman VF: Uncoupling of complex regulatory patterning during evolution of larval development in echinoderms. BMC Biol. 2010, 8: 143-10.1186/1741-7007-8-143.PubMed CentralView ArticlePubMed
- Hejnol A, Martindale MQ: Coordinated spatial and temporal expression of Hox genes during embryogenesis in the acoel Convolutriloba longifissura. BMC Biol. 2009, 7: 65-10.1186/1741-7007-7-65.PubMed CentralView ArticlePubMed
- Cameron RA, Rowen L, Nesbitt R, Bloom S, Rast JP, Berney K, Arenas-Mena C, Martinez P, Lucas S, Richardson PM, Davidson EH, Peterson KJ, Hood L: Unusual gene order and organization of the sea urchin hox cluster. J Exp Zool B Mol Dev Evol. 2006, 306: 45-58.View ArticlePubMed
- Long S, Martinez P, Chen WC, Thorndyke M, Byrne M: Evolution of echinoderms may not have required modification of the ancestral deuterostome HOX gene cluster: first report of PG4 and PG5 Hox orthologues in echinoderms. Dev Genes Evol. 2003, 213: 573-576. 10.1007/s00427-003-0355-3.View ArticlePubMed
- Kriks S, Lanuza GM, Mizuguchi R, Nakafuku M, Goulding M: Gsh2 is required for the repression of Ngn1 and specification of dorsal interneuron fate in the spinal cord. Development. 2005, 132: 2991-3002. 10.1242/dev.01878.View ArticlePubMed
- Furlong RF, Younger R, Kasahara M, Reinhardt R, Thorndyke M, Holland PW: A degenerate ParaHox gene cluster in a degenerate vertebrate. Mol Biol Evol. 2007, 24: 2681-2686. 10.1093/molbev/msm194.View ArticlePubMed
- Garcia-Fernandez J: The genesis and evolution of homeobox gene clusters. Nat Rev Genet. 2005, 6: 881-892.View ArticlePubMed
- Jakob W, Sagasser S, Dellaporta S, Holland P, Kuhn K, Schierwater B: The Trox-2 Hox/ParaHox gene of Trichoplax (Placozoa) marks an epithelial boundary. Dev Genes Evol. 2004, 214: 170-175. 10.1007/s00427-004-0390-8.View ArticlePubMed
- Cresko WA, Yan YL, Baltrus DA, Amores A, Singer A, Rodríguez-Marí A, Postlethwait JH: Genome duplication, subfunction partitioning, and lineage divergence: Sox9 in stickleback and zebrafish. Dev Dyn. 2003, 228: 480-489. 10.1002/dvdy.10424.View ArticlePubMed
- Martinez P, Rast JP, Arenas-Mena C, Davidson EH: Organization of an echinoderm Hox gene cluster. Proc Natl Acad Sci U S A. 1999, 96: 1469-1474. 10.1073/pnas.96.4.1469.PubMed CentralView ArticlePubMed
- Hinman VF, Nguyen AT, Davidson EH: Expression and function of a starfish Otx ortholog, Am Otx: a conserved role for Otx proteins in endoderm development that predates divergence of the eleutherozoa. Mech Dev. 2003, 120: 1165-1176. 10.1016/j.mod.2003.08.002.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.