- Research article
- Open Access
Hox gene cluster of the ascidian, Halocynthia roretzi, reveals multiple ancient steps of cluster disintegration during ascidian evolution
Zoological Lettersvolume 3, Article number: 17 (2017)
The Update to this article has been published in Zoological Letters 2019 5:8
Hox gene clusters with at least 13 paralog group (PG) members are common in vertebrate genomes and in that of amphioxus. Ascidians, which belong to the subphylum Tunicata (Urochordata), are phylogenetically positioned between vertebrates and amphioxus, and traditionally divided into two groups: the Pleurogona and the Enterogona. An enterogonan ascidian, Ciona intestinalis (Ci), possesses nine Hox genes localized on two chromosomes; thus, the Hox gene cluster is disintegrated. We investigated the Hox gene cluster of a pleurogonan ascidian, Halocynthia roretzi (Hr) to investigate whether Hox gene cluster disintegration is common among ascidians, and if so, how such disintegration occurred during ascidian or tunicate evolution.
Our phylogenetic analysis reveals that the Hr Hox gene complement comprises nine members, including one with a relatively divergent Hox homeodomain sequence. Eight of nine Hr Hox genes were orthologous to Ci-Hox1, 2, 3, 4, 5, 10, 12 and 13. Following the phylogenetic classification into 13 PGs, we designated Hr Hox genes as Hox1, 2, 3, 4, 5, 10, 11/12/13.a, 11/12/13.b and HoxX. To address the chromosomal arrangement of the nine Hox genes, we performed two-color chromosomal fluorescent in situ hybridization, which revealed that the nine Hox genes are localized on a single chromosome in Hr, distinct from their arrangement in Ci. We further examined the order of the nine Hox genes on the chromosome by chromosome/scaffold walking. This analysis suggested a gene order of Hox1, 11/12/13.b, 11/12/13.a, 10, 5, X, followed by either Hox4, 3, 2 or Hox2, 3, 4 on the chromosome. Based on the present results and those previously reported in Ci, we discuss the establishment of the Hox gene complement and disintegration of Hox gene clusters during the course of ascidian or tunicate evolution.
The Hox gene cluster and the genome must have experienced extensive reorganization during the course of evolution from the ancestral tunicate to Hr and Ci. Nevertheless, some features are shared in Hox gene components and gene arrangement on the chromosomes, suggesting that Hox gene cluster disintegration in ascidians involved early events common to tunicates as well as later ascidian lineage-specific events.
Hox genes comprise a subset of Antp class homeobox genes  conserved throughout animal phylogeny, and are closely involved in morphogenetic patterning along the anterior-posterior axis. In chordates, Hox genes are classified into 13 paralog groups (PGs) according to homeodomain similarity . Hox genes are often found in a relatively narrow region on one chromosome, forming a Hox gene cluster. It is generally accepted that the Hox gene cluster consists of a subset of the 13 PG Hox genes, which are aligned according to PG number, and that have the same transcription direction, spanning about 100–120 kb on a chromosome. These characteristics of the Hox gene cluster are almost exclusively observed in vertebrate genomes . In contrast, invertebrate Hox gene clusters span the chromosome much more broadly than do their vertebrate counterparts , though the number of Hox genes that constitute a single cluster in invertebrates is at most 13, with the exception of amphioxus and lepidopteran insects, which possess 15 [4, 5] and 14 or more Hox genes , respectively. It is now accepted that Hox gene clusters in vertebrates are exceptionally tightly organized and that the structure of the Hox gene cluster, or the placement of Hox genes on chromosomes, is more variable among invertebrates . The biological significance of such Hox gene clustering in a relatively small portion of the genome has only been explained to a certain extent in vertebrates, while in other taxa it remains poorly understood .
Ascidians belong to the class Ascidiacea, subphylum Tunicata (Urochordata), and phylum Chordata . Ascidians occupy a phylogenetic position between vertebrates and amphioxus [8, 9]. Amphioxus, a basal chordate, has a single Hox gene cluster, consisting of 15 Hox genes with the same transcription direction, spanning about 470 kb . In the ascidian, Ciona intestinalis (Ci), nine Hox genes were identified during draft genome analysis . It was subsequently revealed that the nine Ci Hox genes are on two chromosomes, seven on one and two on the other, exhibiting an unusual gene order as shown using chromosomal FISH analysis by our group . Based on these observations, it was suggested that the Hox cluster has disintegrated in Ci, with the loss of some genes and changes in gene placement on the chromosome . It is thus anticipated that the Hox gene cluster may have also disintegrated in other ascidians. However, substantial evidence to support this speculation has yet to be reported. If such disintegration did in fact occur, the process responsible for its occurrence in ascidian evolution remains enigmatic. In the present study, we address these points.
Ascidians are traditionally divided into two groups (subclasses) , Enterogona (Aplousobranchia and Phlebobranchia) and Pleurogona (Stolidobranchia). Ci is a member of the Phlebobranchia, and Hr belongs to the Stolidobranchia. Both species are widely used in scientific research, especially in developmental studies, and their embryos exhibit very similar development, including most of the same cell lineages . Nevertheless, the non-protein coding regions of their genomes are difficult to align , reflecting a remote phylogenetic relationship between these two ascidians.
In the present study, we analyzed the Hox gene complement and its organization in the Hr genome by chromosomal FISH and chromosome walking to clarify the Hox gene cluster structure. We show that the Hox gene complement consists of nine genes in Hr, as in Ci, but unexpectedly, the nine Hox genes of Hr reside on a single chromosome. We further inferred the Hox gene order on the chromosome. By comparing the information of Hr with that of Ci as well as of other lower chordates, we propose a scenario to explain how the disintegration occurred during ascidian or tunicate evolution.
Isolation of Halocynthia roretzi genomic DNA and Hox gene candidates
Halocynthia roretzi genomic DNAs were prepared individually from single adult animals. Gonads were excised and frozen in liquid nitrogen. Genomic DNA was prepared according to the protocol of Blin and Stafford  and used for PCR as templates. We found that genomic PCR was occasionally not successful with some DNA preparation, probably due to genomic sequence heterogeneity among individuals.
Genomic PCR for isolation of Hox gene candidates was performed as described previously . Additionally, the following degenerate primer sets were used: 5′GARYTNGARAARGARTTY3′ (corresponding to ELEKEF), 5′AARAARMGNCARCCNTAY3′ (KKRQPY) and 5′NCKNCKRTTYTGRAACCA3′ (WFQNRR).
Phylogenetic analysis of Halocynthia roretzi Hox gene candidates
Phylogenetic analysis of Hox gene candidates was done using CLUSTAL W for alignment and the MEGA5 software package  to construct ML trees with 1000 trials. For the reference data set, homeodomains with the flanking 20 N-terminal and seven C-terminal amino acid residues (87 amino acid residues) of 39 mouse (Mus musculus), 35 coelacanth (Latimeria menadoensis), 21 horn shark (Heterodontus francisci) and nine ascidian (Ciona intestinalis; Ci) Hox proteins were used (see Additional file 1: Figure S1).
Construction of a BAC library of Halocynthia roretzi
A Halocynthia roretzi BAC library was constructed from sperm DNA of a single adult H. roretzi, which was obtained at Otsuchi Marine Research Center of the University Tokyo, in Iwate, Japan.
A BAC library was constructed essentially as described previously . The sperm was washed twice with phosphate-buffered saline and then with lysis buffer (10 mM Tris-HCl pH 8.0, 50 mM NaCl, 1% lithium dodecyl sulfate, 100 mM EDTA) for 2 h at 37 °C in a 1.0% agarose gel plug. The plug was stored in 20% NDS solution (0.2% N-lauryl sarcosine, 2 mM Tris-HCl pH 9.0, 0.14 M EDTA). After exchanging the NDS buffer with TE, genomic DNA was digested partially with BamH1, and 150–250 kb DNA fragments were isolated using pulsed field gel electrophoresis. DNA fragments eluted from the gel were ligated into pKS145 and an aliquot of the ligation reaction mixture was used to transform E. coli DH10B. Using a Flexys robot (Genomic Solutions, USA) and a 3D:Biomek FX robot (Beckman Coulter, USA), a total of 20,736 BAC clones were picked up and arrayed in 54 × 384-well microtiter plates in LB medium containing 10% glycerol and 25 μg/mL ampicillin. Plates were incubated overnight at 37 °C and then stored at –80 °C. The BAC library was constructed so as to be amenable to the dimension pooling system for PCR screening.
To estimate the size of the average genomic DNA insert in the BAC library, 20 randomly selected clones were digested with NotI and analyzed by pulsed field gel electrophoresis. This analysis revealed that the average insert size was 110 kbp, and the coverage was estimated as 14.2 ×, assuming that the genome size of H. roretzi is 160 Mbp .
Embryos for chromosomal FISH (fluorescence in situ hybridization)
Fertilized eggs were raised in FSW until the 2-cell stage. When embryos began to divide to the 4-cell stage, they were transferred into Ca++ free seawater. When they reached the 64-cell stage, colchicine (Sigma) was added at a final concentration of 0.025% (w/v). Embryos were cultured for 30 min and fixed with acetic acid:methanol (3:1) overnight, and then transferred to 70% ethanol and kept at –20 °C until use.
Fluorescence in situ hybridization (FISH)
Two-color chromosomal FISH was performed according to a procedure previously described for Ci with some modifications . For preparation of metaphase spreads, 15–30 fixed embryos were de-chorionated manually under a stereomicroscope and embryonic cells were transferred to a microtube. After removal of excess liquid, 100 μL of 60% acetic acid was added to the tube, and the cell suspension was mixed by gentle rolling for 90 s. Then, the mixture was agitated for 30 s by gentle pipetting about 20 times. Immediately after agitation, the mixture was spread gently on a warmed (48 °C) clean slide glass using a pipet. The glass slide was allowed to stand at 48 °C for 2.5 h before being subjected to FISH. Probes for chromosomal FISH were prepared using BAC clones labeled with biotin or digoxigenin using a nick translation kit (Roche).
Chromosome/scaffold walking with BAC library and ANISEED database
Chromosome walking using PCR and BAC library screening was performed using standard methods [18, 19]. Scaffold walking is here referred to as a modification of chromosome walking, using nucleotide sequence information of scaffolds out of database open to public for designing PCR for BAC library screening to identify adjacent scaffolds.
A pair of primers was designed around the starting BAC clone end region. The primers were used for PCR screening of the BAC library. When positive clones were found, their DNAs were prepared and sequenced from both ends of the insert. The resulting nucleotide sequences were used for designing PCR primers. By using the primers and DNAs of isolated clones and of the starting BAC clone, reciprocal PCR was performed to determine the relative positional relationship between isolated clones and the end region of the starting BAC clone, and a desired clone was identified. By using the identified clone, another screening cycle was performed to walk further along the chromosome.
In the case of scaffold walking, the resulting nucleotide sequences were used for BLAST surveys of the Halocynthia roretzi genomic sequences, Halocynthia roretzi MTP2014, of the ANISEED genomic database  (https://www.aniseed.cnrs.fr/aniseed/default/blast_search) to locate the ends of the BAC clone on scaffolds. A pair of primers was designed in the end region of the scaffold, tested for compatibility with PCR using several genomic DNA preparations, and used for the screening of the BAC library. When the positive clones were isolated, the nucleotide sequences were determined for both end regions, and the resulting nucleotide sequence information was used for the BLAST surveys to identify an adjacent scaffold so as to walk further along the chromosome.
Nucleotide sequence determination of the BAC clone end regions
BAC clone DNA was prepared from 5 mL overnight culture using a QIAprep Spin Miniprep kit (QIAGEN). The BAC end-region sequence was determined by a standard method using BigDye ver. 3, upper and lower primers for pKS145 vectors and ABI 3000 sequencers.
Nucleotide sequence determination of a BAC clone insert
Circular BAC clone DNA was prepared using a QIAGEN Large-Construct kit according to the supplier’s procedure, in which an exonuclease digestion procedure is included (QIAGEN). For nucleotide sequence determination of BAC clone DNAs, an Illumina Miseq was used. Libraries were prepared according to a protocol provided by the manufacturer, with slight modifications. Fragmented BAC clone DNA was further purified using Blue Pippin (Sage Science). A paired-end library consisting of clones containing ∼720 bp insert DNA fragment was prepared for the Miseq using a TruSeq DNA PCR-Free LT Sample Prep Kit (Illumina). Adapter sequences were removed from all sequence reads using Trimmomatic-0.30 . Paired-end reads of high quality (quality-value ≥20) were assembled de novo using Newbler 2.9 (GS Assembler) to create a scaffold. From the scaffold, vector sequence was removed, and genomic sequence was extracted.
Hox gene complement in Hr genome
In a previous study, we reported the isolation of Hox1, as well as Hox gene fragments from the Hr genome. Hox1 was identified by alignment to other homeodomain sequences available at that time, and its structure and developmental expression were reported . Since then, we have repeated genomic PCR with various sets of degenerate primers and RT-PCR, using RNA from embryos at various stages, and eventually isolated nine Hox gene candidate sequences, including the previously reported Hox1 from the Hr genome. These candidates were subjected to phylogenetic analysis in the present study.
A phylogenetic tree (Fig. 1) identified nine Hr Hox gene candidates (see also Additional file 1: Figure S1). Eight of the nine Hr Hox genes always clustered with Ci Hox genes, Ci-Hox1, 2, 3, 4, 5, 10, 12 and 13. This suggests that eight Hr Hox genes and their respective counterparts in Ci are orthologous. The remaining Hr Hox gene candidate did not show significant similarity to Ci-Hox6 or any other Ci Hox gene. Accordingly, we tentatively designated the nine Hr Hox genes as Hr Hox1, 2, 3, 4, 5, 10, 12, 13 and X.
Next, we asked to which of the 13 paralog groups (PGs) the nine Hr Hox genes belong. In the phylogenetic tree, Hr candidates for Hox1, Hox2, Hox3, and Hox4 were clearly classified into PGs 1, 2, 3 and 4 with bootstrap values of 95%, 96%, 83% and 91%, respectively (Fig. 1). Although Hr Hox5, 10, 12, 13 and X could not be classified into single PGs in this tree, a clade consisting of Hox genes of PGs 1–8 was supported by bootstrap values of 75% (Fig. 1); both Hox5 and X genes may thus be classified into PGs 1–8. Since Hox genes of PGs 1–4 were clearly identified, another tree was constructed to determine to which PGs the remaining two Hox genes could be classified. In a tree using Hox genes of PGs 5–8, Hr Hox5 and Ci-Hox5 were likely grouped into PG5 with a bootstrap value of 73% (Fig. 2a). By contrast, HoxX was hardly classified into any PG (Fig. 2a). It is also noted here that Ci-Hox6 could hardly be classified into PG6 (see Discussion).
Regarding the three posterior genes, Hox10 was classified into PG10, and Hox12 and 13 into PGs 11–13, albeit with poor bootstrap support (Fig. 1). Assignment of Hox10 as PG10 gene was supported by the conservation of four diagnostic residues (Gly 1, Glu 29, Leu 32, and Asp 42) in the homeodomain and three Lys residues in the flanking C-terminal side region (Additional file 1: Figure S1). These criteria were previously used to assign Ci-Hox10 as such . The remaining two posterior genes were in a clade consisting of Hox genes of PGs11–13 with relatively low bootstrap values (61%, Fig. 1). In another tree using Hox genes of PGs 9–13, the clustering of the two genes in the clade consisting of PGs 11, 12, and 13 was supported by a bootstrap value of 82% (Fig. 2b). Therefore, we propose Hr Hox11/12/13.a and Hox11/12/13.b as counterparts of Ci-Hox12 and Ci-Hox13, respectively.
Thus, the Hr Hox gene complement consists of nine members. These Hox genes are designated Harore Hox1, Hox2, Hox3, Hox4, Hox5, Hox10, Hox11/12/13.a and Hox11/12/13.b according to the newly proposed nomenclature for ascidian genes , and the remaining one gene is tentatively designated Harore HoxX.
Hox gene cluster structural analysis using chromosomal FISH
In order to address the genomic organization of the nine Hox genes using chromosomal FISH, we screened the BAC genomic library developed from sperm of a single individual and obtained clones for all nine Hr Hox genes. Among the clones isolated, clones containing as many as three Hox genes, except for Hox1, were found (data not shown, see next section). Using the isolated BAC clones for probes, we carried out FISH on chromosome spreads prepared from cleavage stage Hr embryos.
In Fig. 3a, red and green spots corresponding to two BAC clones (one containing Hox1 and the other containing Hox2, Hox3 and Hox4, respectively) are shown located on the same chromosome. It is also noted here that Hox1 is localized closer to the chromosome end than Hox2, 3 and 4. In Fig. 3b, a green spot corresponding to the BAC clone containing Hox10, 11/12/13.a and Hox11/12/13.b was localized on the chromosome with a red spot representing Hox1. The red spot was closer to the chromosome end than the green spot (Fig. 3b). Figure 3c shows green and red spots corresponding to two BAC clones, one containing Hox5 and HoxX and the other containing the three posterior genes mentioned above, overlapping on a single chromosome. In Fig. 3d, Hox1 is localized closer to the chromosome end than Hox5, HoxX or the three posterior Hox genes. These results indicate that all of the nine Hox genes are present on a single chromosome in the Hr genome, unlike the Ci genome. In addition, these results suggest that eight of the nine Hox genes are localized closely together on the chromosome, although the positional relationship among the eight genes could not be determined in this analysis. In contrast, Hox1 was localized relatively close to the chromosome end, away from the other Hox genes. This arrangement is somewhat similar to that in the Ci genome (see Discussion), except that Hox12 and 13 genes are on a different chromosome in Ci.
Hox gene order on the chromosome as inferred by chromosome/scaffold walking
Since chromosomal FISH analysis suggested that eight of nine Hox genes, excepting Hox1, may be localized in close proximity on a single chromosome, we examined overlapping of the BAC clones containing at least one of the eight Hox genes by reciprocal PCR using isolated BAC clone DNAs as templates. We found that many clones overlapped one another (data not shown), in a manner suggesting that Hox2, Hox3 and Hox4 form a subcluster and are aligned on the chromosome in this order. Similarly, five Hox genes, Hox5, HoxX, Hox10, Hox11/12/13.a and Hox11/12/13.b, form another subcluster and are aligned in the order, HoxX, Hox5, Hox10, Hox11/12/13.a, Hox11/12/13.b.
In order to resolve the positional relationships between Hox1 and the two subclusters, we carried out chromosome/scaffold walking by utilizing the BAC library and genomic sequence information, Halocynthia roretzi MTP2014, in the genome browser of ANISEED database . Scaffolds including Hox1 and the two subclusters were identified and the chromosome/scaffold walking was started from the ends of these scaffolds. We were successful in connecting the two scaffolds, S11 and S54, which contained Hox1 and Hox11/12/13.b, respectively (Fig. 4). At least five scaffolds were located between Hox1 and Hox11/12/13.b, and the distance between the two genes was about 1.53 Mbp according to calculations based on each scaffold length in the ANISEED database (Fig. 4). Walking distal to HoxX, no adjoining clone was isolated after isolation of a BAC clone (32G9 in Fig. 4), and the chromosome/scaffold walking was aborted. Similarly, walking distal to Hox4 was aborted, because no clone was isolated to connect scaffold S201 with its adjacent scaffold (Fig. 4). On the other hand, chromosome/scaffold walking distal to Hox2 yielded many clones at the end of scaffold S36 (Fig. 4). The clones contained similar, but not identical, nucleotide sequences at one end, and when used for BLAST queries to search the database, every clone hit many short scaffolds containing similar sequences. As a result, the scaffold neighboring S36 was not determined. These observations suggest that the Hox gene order on the chromosome is Hox1, Hox11/12/13.b, Hox11/12/13.a, Hox10, Hox5, HoxX, followed by either Hox4, Hox3, Hox2 or Hox 2, Hox3, Hox4, from the chromosome end to center (Fig. 4). In either case, the nine Hr Hox genes are estimated to span at least ~2.3 Mbp on the chromosome. A mir10 sequence that has been reported to reside in upstream of Hox4 in hemichordates and amphioxus [24,25,26] was not found in the Hr genome using BLAST survey over ANISEED genome browsers (data not shown).
Recent phylogenetic studies suggest that ascidians, comprising a major group of tunicates, are not monophyletic . Tunicates are divided into two branches; one includes Stolidobranchia (Hr) and Appendicularia, and the other includes Phlebobranchia (Ci), Aplousobranchia (another ascidian group) and Thaliacea . In the present study, we analyzed the Hox gene complement and the Hox gene cluster structure of the stolidobranchian ascidian, Halocynthia roretzi, which is phylogenetically remote from the phlebobranchian ascidian, Ciona intestinalis. The Hr Hox gene complement consists of nine members, the same as that of Ci. The nine Hr Hox genes are located on a single chromosome, unlike Ci, in which they reside on two chromosomes .
The Hox gene complement in the last common ancestor of Hr and Ci
The present phylogenetic analysis suggested that eight of the nine Hr Hox gene complement are Hr orthologs of eight Ci Hox genes. When the remaining gene, Hr HoxX, was used to query various ascidian genomes in the ANISEED database, a Hox gene that exhibits high similarity was found only in the genome of a closely related species, Halocynthia aurantium (data not shown). On the other hand, in the ascidian species, Ci and Ciona savignyi, the best-hit Hox genes were Ci-Hox6 and a probable Cs ortholog of Ci-Hox6, respectively (data not shown). In other ascidians, the best-hit gene was difficult to assign to a single PG during phylogenetic analysis (data not shown).
In the previous study, Ci-Hox6 was tentatively designated as such, but without reliable evidence other than that it is localized proximal to Ci-Hox5 (~3 kb apart, according to the ANISEED genome browser) . In the present study, the phylogenetic position of the gene was close to, but not within the clade of PG6 genes, and it exhibited some affinity for PGs 7 and 8 (Fig. 1). The designation, Ci-Hox6, should thus be revisited in future studies.
Based on the above observations, we suggest that the last common ancestor to Hr and Ci may have possessed three central Hox genes, one of PG4, one of PG5 Hox genes, and one out of PGs 6–8 genes, and that the last one may have evolved in the lineages to Hr and Ci, and resulted in extant Harore HoxX and Ci-Hox6, respectively.
The Hox gene complement of ascidians in comparison to that of amphioxus
The number of Hox genes in the Hr and Ci genomes is smaller than in amphioxus, which has 15 Hox genes [4, 5]. It seems reasonable that the ancestral tunicate may have lost several Hox genes after diverging from the evolutionary lineage leading to the vertebrates. However, the Hox genes of amphioxus are designated according to their order on the chromosome, not necessarily based on PGs, which were originally invented for classification of vertebrate Hox genes . In recent years, studies have examined the relationship between amphioxus Hox genes and PGs, employing methods independent of phylogenetic tree construction [28, 29]. The results of these studies, although not necessarily concordant, suggest the common presence of each of PGs 1–5 genes and amphioxus-specific posterior gene paralogs in the genome [28, 29]. As regards the latter, it has been proposed that posterior Hox genes expanded in the amphioxus lineage, and that the last common ancestor of amphioxus and vertebrates may have possessed three ancestral posterior genes, PG9/10, PG11/12 and PG13/14 genes . On the other hand, with respect to the anterior and central Hox genes, it is generally accepted that the last common ancestor of amphioxus and vertebrates possessed three anterior (PGs 1 through 3) and five central (PGs 4 through 8) Hox genes .
When amphioxus Hox protein sequences were analyzed in our phylogenetic tree (Additional file 1: Figure S1 and Additional file 2: Figure S2), amphioxus Hox1, 2, 3 and 4 were clearly classified into PGs1, 2, 3 and 4, respectively. Hox5 may be classified into PG5, although the clustering was not supported by a high bootstrap value. Amphioxus Hox6, 7, and 8 were grouped into a clade consisting of PGs 4–8 genes, but could be excluded from PG4 and PG5; thus, the three genes may be classified into PGs 6–8. Amphioxus Hox10, 11 and 12 apparently seemed to be paralogs (see Additional file 2: Figure S2) and were classified, together with Hox9, into a clade consisting of PGs 9 and 10, which was barely supported by a low bootstrap value. Similarly, amphioxus Hox13 and 14 seemed to be paralogs and were classified into PGs 11–13. By contrast, amphioxus Hox15 was apparently classified into PG13 with a relatively high bootstrap value (84%, Additional file 2: Figure S2).
Considering our observations and those of others, we speculate that the Hox gene cluster of the ancestral amphioxus (and the last common ancestor for amphioxus and vertebrates, too) comprised 11 Hox genes, including three anterior (each of PGs 1-3), five central (each of PGs 4–8), and three posterior (out of PGs 9–13) Hox genes. If this is the case, the last common ancestor of Hr and Ci must have lost two central Hox genes out of PGs 6–8 after divergence from the lineage continuing from the ancestral chordate to the ancestral vertebrate.
Disintegration of the Hox gene cluster during evolution of Hr and Ci
In comparison with Ci, disintegration of the Hr Hox gene cluster seems less extensive, in that all nine Hox genes are on a single chromosome (Fig. 4). Thus, the Hox gene cluster of the last common ancestor of Hr and Ci, appears to have disintegrated differently in these two ascidians’ evolutionary lineages.
Nevertheless, there are some structural features shared by Hr and Ci. First, Hox1 is located away from other Hox genes in both Hr and Ci (Fig. 4, ). Second, Hox11/12/13.a and Hox11/12/13.b are adjacent to each other with reversed orientation. This situation is the same in Ci counterparts (Ci-Hox12 and Ci-Hox13). Third, in both ascidians, Hox2, Hox3, and Hox4 are aligned in the same direction without intervening genes. It should be noted that in both ascidian genomes, the gene immediately adjacent to Hox2 is STAC (SH3 and cysteine-rich domain-containing protein) and two neighboring genes to Hox4 are CHST (carbohydrate sulfotransferase) and NEBL (Nebullet) (Additional file 3: Figure S3). This is the only conserved gene arrangement surrounding the Hox genes that we observed in these two ascidians. This suggests two possibilities. First, after divergence of the two lineages to Hr and Ci, genomic shuffling occurred in each lineage to such an extent that conservation of the gene arrangement surrounding the Hox genes was limited only to one small region, about 170 kb and 140 kb in Hr and Ci, respectively (see Additional file 3: Figure S3). Second, but more importantly, the gene arrangement observed in common between Hr and Ci must have been established prior to the divergence of Hr and Ci.
From these shared structural features, it appears that the disintegration of the Hox gene cluster must have included certain early events, such as translocation of Hox1 or of the Hox2, 3, 4 group, and tail-to-tail location of the Hox11/12/13.a and 11/12/13.b pair. These changes in the Hox gene cluster must have occurred in the last common ancestor of Hr and Ci.
Conclusion: a theoretical scenario for the disintegration of the Hox gene cluster in the ascidian or tunicate evolution
In an appendicularian tunicate, Oikopleura dioica, all central Hox genes and the PG3 Hox gene are missing, and the Hox gene complement in this species is quite different from that of Hr or Ci . Considering this and information about Hox gene cluster of amphioxus, a simple scenario for the disintegration of the Hox gene cluster during the course of ascidian or tunicate evolution is as shown in Fig. 5.
In this scheme, 1) when the ancestral chordate emerged, it had a single Hox gene cluster consisting of three anterior (PGs 1–3), five central (PGs 4–8) and three ancestral posterior (PG9/10, PG11/12 and PG13/14) genes . 2) The ancestral chordate evolved, and the last common ancestor of tunicates and vertebrates diverged from the lineage to cephalochordate. 3) When the ancestral tunicate diverged from the lineage to vertebrates, it must have experienced extensive genomic rearrangement and lost at least one (or two) central Hox genes. At the same time, the ancestral tunicate likely came to possess tunicate characteristics, and a Hox gene complement consisting of nine genes (three each of anterior, central, and posterior Hox genes) was established. Meanwhile, early disintegration events in the Hox gene cluster occurred. Loss of the central Hox genes and disintegration of the Hox gene cluster may be correlated with peculiar way of development of tunicates  and/or limited function of Hox genes as observed in the early development of Ci . 4) The ancestral tunicate evolved and diverged into two distinct lineages, and in turn, ancestral ascidians of the Pleurogona (Stolidobranchia) and Enterogona (Phlebobranchia and Aplousobranchia) diverged from Appendicularia and Thaliacea, respectively. The Hox gene cluster as well as the genome must have experienced further genomic rearrangement. The relatively small conserved gene arrangement between Hr and Ci in the regions surrounding Hox genes may support this part of the scenario.
In the above simple scenario for the disintegration of the ascidian Hox gene cluster, it remains unresolved why one putative central Hox gene has diverged considerably more than other Hox genes in Hr. The evolutionary constraints governing the disintegration of the Hox gene cluster in Hr or Ci, which apparently occurred to a much smaller extent than in Appendicularia, also remain unknown. Answering these questions will further clarify the characteristic features of the Hox gene cluster in ascidians and/or tunicates.
Bürglin TR, Affolter M. Homeodomain proteins: an update. Chromosoma. 2016;125:497–521.
Scott MP. A rational nomenclature for vertebrate homeobox (HOX) genes. Nucleic Acids Res. 1993;21:1687–8.
Duboule D. The rise and fall of Hox gene clusters. Development. 2007;134:2549–60.
Holland L, Albalat R, Azumi K, Benito-Gutiérrez E, Blow M, Bronner-Fraser M, et al. The amphioxus genome illuminates vertebrate origins and cephalochordate biology. Genome Res. 2008;18:1100–11.
Takatori N, Butts T, Candiani S, Pestarino M, Ferrier DEK, Saiga H, et al. Comprehensive survey and classification of homeobox genes in the genome of amphioxus, Branchiostoma Floridae. Dev Genes Evol. 2008;218:579–90.
Ferguson L, Marlétaz F, Carter JM, Taylor WR, Gibbs M, Breuker CJ, et al. Ancient Expansion of the Hox Cluster in Lepidoptera Generated Four Homeobox Genes Implicated in Extra-Embryonic Tissue Formation. PLoS Genet. 2014;10(10):e1004698.
Satoh N. Developmental biology of ascidians. New York: Cambridge University Press; 1994.
Delsuc F, Brinkmann H, Chourrout DPH. Tunicates and not cephalochordates are the closest living relatives of vertebrates. Nature. 2006;439:965–8.
Putnam NH, Butts T, Ferrier DEK, Furlong RF, Hellsten U, Kawashima T, et al. The amphioxus genome and the evolution of the chordate karyotype. Nature. 2008;453:1064–71.
Dehal P, Satou Y, Campbell RK, Chapman J, Degnan B, De Tomaso A, et al. The draft genome of Ciona intestinalis: insights into chordate and vertebrate origins. Science. 2002;298:2157–67.
Ikuta T, Yoshida N, Satoh N, Saiga H. Ciona intestinalis Hox gene cluster: its dispersed structure and residual colinear expression in development. Proc Natl Acad Sci U S A. 2004;101:15118–23.
Ikuta T, Saiga H. Organization of Hox genes in ascidians: present, past, and future. Dev Dyn. 2005;233:382–9.
Lemaire P. Evolutionary crossroads in developmental biology: the tunicates. Development. 2011;138:2143–52.
Blin N, Stafford DW. A general method for isolation of high molecular weight DNA from eukaryotes. Nucleic Acids Res. 1976;3:2303–8.
Katsuyama Y, Wada S, Yasugi S, Saiga H. Expression of the labial group Hox gene HrHox-1 and its alteration induced by retinoic acid in development of the ascidian Halocynthia roretzi. Development. 1995;121:3197–205.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.
Fujiyama A, Watanabe H, Toyoda A, Taylor TD, Itoh T, Tsai SF, Park HS, Yaspo ML, Lehrach H, Chen Z, Fu G, Saitou N, Osoegawa K, de Jong PJ, Suto Y, Hattori MSY. Construction and analysis of a human-chimpanzee comparative clone map. Science. 2002;295:131–4.
Cai L, Taylor JF, Wing RA, Gallagher DS, Woo S-S, Davis SK. Construction and characterization of a bovine bacterial artificial chromosome library. Genomics. 1995;29:413–25.
Kim U-J, Birren BW, Slepak T, Mancino V, Boysen C, Kang H-L, et al. Construction and characterization of a human bacterial artificial chromosome library. Genomics. 1996;34:213–8.
Brozovic M, Martin C, Dantec C, Dauga D, Mendez M, Simion P, et al. ANISEED 2015: a digital framework for the comparative developmental biology of ascidians. Nucleic Acids Res. 2016;44:D808–18.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
Wada S, Tokuoka M, Shoguchi E, Kobayashi K, Di Gregorio A, Spagnuolo A, et al. A genomewide survey of developmentally relevant genes in Ciona intestinalis. II. Genes for homeobox transcription factors. Dev. Genes Evol. 2003;213:222–34.
Stolfi A, Sasakura Y, Chalopin D, Satou Y, Christiaen L, Dantec C, et al. Guidelines for the nomenclature of genetic elements in tunicate genomes. Genesis. 2015;53:1–14.
Tanzer A, Amemiya CT, Kim CB, Stadler PF. Evolution of microRNAs located within Hox gene clusters. J. Exp. Zool. Part B Mol. Dev. Evol. 2005;304B:75–85.
Freeman R, Ikuta T, Wu M, Koyanagi R, Kawashima T, Tagawa K, et al. Identical genomic organization of two hemichordate Hox clusters. Curr Biol. 2012;22:2053–8.
Amemiya CT, Prohaska SJ, Hill-Force A, Cook A, Wasserscheid J, Ferrier DEK, et al. The amphioxus Hox cluster: characterization, comparative genomics, and evolution. J Exp Zool Part B Mol Dev Evol. 2008;310:465–77.
Tsagkogeorga G, Turon X, Hopcroft RR, Tilak M-K, Feldstein T, Shenkar N, et al. An updated 18S rRNA phylogeny of tunicates based on mixture and secondary structure models. BMC Evol Biol. 2009;9:187.
Hueber SD, Weiller GF, Djordjevic MA, Frickey T. Improving Hox protein classification across the major model organisms. PLoS One. 2010;5:5.
Thomas-Chollier M, Ledent V, Leyns L, Vervoort M. A non-tree-based comprehensive study of metazoan Hox and ParaHox genes prompts new insights into their origin and evolution. BMC Evol Biol. 2010;10:73.
Pascual-Anaya J, D’Aniello S, Kuratani S, Garcia-Fernàndez J. Evolution of Hox gene clusters in deuterostomes. BMC Dev Biol. 2013;13:26.
Seo HC, Edvardsen RB, Maeland AD, Bjordal M, Jensen MF, Hansen A, et al. Hox cluster disintegration with persistent anteroposterior order of expression in Oikopleura dioica. Nature. 2004;431:67–71.
Ikuta T, Satoh N, Saiga H. Limited functions of Hox genes in the larval development of the ascidian Ciona intestinalis. Development. 2010;137:1505–13.
We thank Drs. You Katsuyama (Shiga University of Medical Science) and Shuichi Wada (Nagahama Institute of Bio-Science and Technology) for initial contributions to this work. We also thank Drs. Steven Aird (OIST) for technical editing of the manuscript, Masanori Taira (University of Tokyo) for working facilities, Masaru Nonaka (University of Tokyo), and Peter Holland (Oxford University) for valuable discussions.
This work was partly supported by JSPS KAKENHI Grant Numbers 25440113, 22570207, 17018018 and 18370088 to HS and OIST Internal Funds to NS.
Availability of data and materials
All data generated and analyzed in the present study, except for the nucleotide sequences of genes and BAC clones, are included in this article. Nucleotide sequence data of the Hr Hox gene fragments encoding 87 amino acid residues including the homeodomain and BAC clones used for chromosome/scaffold walking have been deposited under accession numbers LC272074-LC272082 (for nine genes), DE999384-DE999451 (for 68 end regions) and AP018119-AP018128 (for 10 whole inserts) in DDBJ.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Amino acid sequences used for the analysis of Hox genes of Halocynthia roretzi (Hr) by construction of ML phylogenetic trees. Amino acid sequences include the homeodomain (60 residues in yellow) and the adjacent 20 N-terminal and seven C-terminal residues. Original accession numbers for these sequences are indicated in brackets. Taxonomic abbreviations are Mm for Mus musculus, Lm for Latimeria menadoensis, Hf for Heterodontus francisci, Ci for Ciona intestinalis, Hr for Halocynthia roretzi, Bl for Branchiostoma lanceolatum and Bf for Branchiostoma floridae. Hr Hox genes, designated according to orthology with Ci-Hox counterparts and according to their classification into paralog groups (PGs) are indicated prior to and in parentheses, respectively. In ascidian sequences, letters in red indicate diagnostic residues for Hox10 homeodomain proteins (see text). (PDF 41 kb)
Phylogenetic analysis of amphioxus Hox genes by constructing an ML tree. The ML tree was constructed using homeodomain sequences and the adjacent 20 N-terminal and seven C-terminal amino acids (Additional file 1: Figure S1) and MEGA5 software. The percentage of 1000 replicated trees, in which gene clustering was supported, is indicated at nodes. Within a clade consisting of only vertebrate Hox genes, the percentage was not indicated at the node. Amphioxus Hox genes are marked by colored circles. Color code and taxonomic abbreviations are the same as in Fig. 1, except that the color code for posterior Hox genes is the same as that shown in Fig. 5. Bl and Bf denote Branchiostoma lanceolatum and Branchiostoma floridae, respectively. (PDF 519 kb)
Conservation of gene arrangements surrounding Hox2, 3, and 4 between Hr and Ci. Genomic regions surrounding Hox2, Hox3, and Hox4 are depicted schematically, based on genomic browser information from the ANISEED database (Halocynthia roretzi MTP 2014, Ciona intestinalis type A (KH2012)). Genes are indicated by thick arrows. The color code for Hox genes is the same as in Fig. 4. Pink arrows downstream of Hox2 indicate the gene encoding SH3 and cysteine-rich domain-containing protein (STAC). Pale pink arrows upstream of Hox4 indicate carbohydrate sulfotransferase (CHST1)/chondroitin 6-O-sulfotransferase (C6ST). Dark pink arrows indicate the NEBL gene encoding the Nebullete protein. Blank arrows indicate genes without positional conservation. Grey arrays of short vertical bars indicate 10 kbp. (PDF 246 kb)