Skip to main content
Log in

Comprehensive survey and classification of homeobox genes in the genome of amphioxus, Branchiostoma floridae

  • Original Article
  • Published:
Development Genes and Evolution Aims and scope Submit manuscript

Abstract

The homeobox genes comprise a large and diverse gene superfamily, many of which encode transcription factors with pivotal roles in the embryonic development of animals. We searched the assembled draft genome sequence of an amphioxus, Branchiostoma floridae, for genes possessing homeobox sequences. Phylogenetic analysis was used to divide these into gene families and classes. The 133 amphioxus homeobox genes comprise 60 ANTP class genes, 29 PRD genes (excluding Pon and Pax1/9), nine TALE genes, seven POU genes, seven LIM genes, five ZF genes, four CUT genes, four HNF genes, three SINE genes, one CERS gene, one PROS gene, and three unclassified genes. Ten of the 11 homeobox gene classes are less diverse in amphioxus than humans, as a result of gene duplication on the vertebrate lineage. Amphioxus possesses at least one member for all of the 96 homeobox gene families inferred to be present in the common ancestor of chordates, including representatives of the Msxlx, Bari, Abox, Nk7, Ro, and Repo gene families that have been lost from tunicates and vertebrates. We find duplication of several homeobox genes in the cephalochordate lineage (Mnx, Evx, Emx, Vent, Nk1, Nedx, Uncx, Lhx2/9, Hmbox, Pou3, and Irx) and several divergent genes that probably originated by extensive sequence divergence (Hx, Ankx, Lcx, Acut, Atale, Azfh, Ahbx, Muxa, Muxb, Aprd1–6, and Ahnf). The analysis reveals not only the repertoire of amphioxus homeobox genes but also gives insight into the evolution of chordate homeobox genes.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Fig. 1
Fig. 2

Similar content being viewed by others

References

  • Abascal F, Zardoya R, Posada D (2005) ProtTest: selection of best-fit models of protein evolution. Bioinformatics 21:2104–2105

    Article  PubMed  CAS  Google Scholar 

  • Booth HAF, Holland PWH (2004) Eleven daughters of NANOG. Genomics 84:229–238

    Article  PubMed  CAS  Google Scholar 

  • Booth HAF, Holland PWH (2007) Annotation, nomenclature and evolution of four novel homeobox genes expressed in the human germ line. Gene 387:7–14

    Article  PubMed  CAS  Google Scholar 

  • Bourlat SJ, Juliusdottir T, Lowe CJ, Freeman R, Aronowicz J, Kirschner M, Lander ES, Thorndyke M, Nakano H, Kohn AB et al (2006) Deuterostome phylogeny reveals monophyletic chordates and the new phylum Xenoturbellida. Nature 444:85–88

    Article  PubMed  CAS  Google Scholar 

  • Brooke NM, Garcia-Fernàndez J, Holland PWH (1998) The ParaHox gene cluster is an evolutionary sister of the Hox gene cluster. Nature 392:920–922

    Article  PubMed  CAS  Google Scholar 

  • Bürglin TR (1994) A comprehensive classification of homeobox genes. In: Duboue D (ed) Guidebook to the homeobox genes. Oxford University Press, Oxford, pp 25–71

    Google Scholar 

  • Bürglin TR (1995) The evolution of homeobox genes. In: Arai R, Kato M, Doi Y (eds) Biodiversity and evolution. The National Science Museum Foundation, Tokyo, pp 291–336

    Google Scholar 

  • Bürglin TR (1997) Analysis of TALE superclass homeobox genes (MEIS, PBC, KNOX, Iroquois, TGIF) reveals a novel domain conserved between plants and animals. Nucleic Acids Res 25:4173–4180

    Article  PubMed  Google Scholar 

  • Bürglin TR (1998) The PBC domain contains a MEINOX domain: coevolution of Hox and TALE homeobox genes. Dev Genes Evol 208:113–116

    Article  PubMed  Google Scholar 

  • Bürglin TR, Cassata G (2002) Loss and gain of domains during evolution of cut superclass homeobox genes. Int J Dev Biol 46:115–123

    PubMed  Google Scholar 

  • Butts T, Holland PWH, Ferrier DEK (2008) The Urbilaterian Super-Hox cluster. Trends in Genetics 24:259–262

    Article  PubMed  CAS  Google Scholar 

  • Candiani S, Castagnola P, Oliveri D, Pestarino M (2002) Cloning and developmental expression of AmphiBrn1/2/4, a POU III gene in amphioxus. Mech Dev 116:231–234

    Article  PubMed  CAS  Google Scholar 

  • Candiani S, Oliveri D, Parodi M, Bertini E, Pestarino M (2006) Expression of AmphiPOU-IV in the developing neural tube and epidermal sensory neural precursors in amphioxus supports a conserved role of class IV POU genes in the sensory cells development. Dev Genes Evol 216:623–633

    Article  PubMed  CAS  Google Scholar 

  • Candiani S, Holland ND, Oliveri D, Parodi M, Pestarino M (2008) Expression of the amphioxus Pit-1 gene (AmphiPOU1F1/Pit-1) exclusively in the developing preoral organ, a putative homolog of the vertebrate adenohypophysis. Brain Res Bull 75:324–330

    Article  PubMed  CAS  Google Scholar 

  • Chi YI, Frantz JD, Oh BC, Hansen L, Dhe-Paganon S, Shoelson SE (2002) Diabetes mutations delineate an atypical POU domain in HNF-1alpha. Mol Cell 10:1129–1137

    Article  PubMed  CAS  Google Scholar 

  • Chu-Lagraff Q, Wright DM, McNeil LK, Doe CQ (1991) The prospero gene encodes a divergent homeodomain protein that controls neuronal identity in Drosophila. Development Suppl 2:79–85

    Google Scholar 

  • Dehal P, Boore JL (2005) Two rounds of whole genome duplication in the ancestral vertebrate. PLoS Biol 3:e314

    Article  PubMed  Google Scholar 

  • Delsuc F, Brinkmann H, Chourrout D, Philippe H (2006) Tunicates and not cephalochordates are the closest living relatives of vertebrates. Nature 439:965–968

    Article  PubMed  CAS  Google Scholar 

  • Dildrop R, Ruther U (2004) Organization of Iroquois genes in fish. Dev Genes Evol 214:267–276

    Article  PubMed  CAS  Google Scholar 

  • Doe CQ, Chu-LaGraff Q, Wright DM, Scott MP (1991) The prospero gene specifies cell fates in the Drosophila central nervous system. Cell 65:451–464

    Article  PubMed  CAS  Google Scholar 

  • Dozier C, Kagoshima H, Niklaus G, Cassata G, Bürglin TR (2001) The Caenorhabditis elegans Six/sine oculis class homeobox gene ceh-32 is required for head morphogenesis. Dev Biol 236:289–303

    Article  PubMed  CAS  Google Scholar 

  • Felsenstein J (1993) PHYLIP (Phylogeny Inference Package) version 3.5c. Department of Genetics, University of Washington, Seattle

    Google Scholar 

  • Ferrier DEK, Minguillon C, Holland PWH, Garcia-Fernàndez J (2000) The amphioxus Hox cluster: deuterostome posterior flexibility and Hox14. Evolution and Development 2:284–293

    Article  PubMed  CAS  Google Scholar 

  • Galliot B, de Vargas C, Miller D (1999) Evolution of homeobox genes: Q50 Paired-like genes founded the Paired class. Dev Genes Evol 209:186–197

    Article  PubMed  CAS  Google Scholar 

  • Garcia-Fernàndez J, Holland PWH (1994) Archetypal organization of the amphioxus Hox gene cluster. Nature 370:563–566

    Article  PubMed  Google Scholar 

  • Glardon S, Holland LZ, Gehring WJ, Holland ND (1998) Isolation and developmental expression of the amphioxus Pax-6 gene (AmphiPax-6): insights into eye and photoreceptor evolution. Development 125:2701–2710

    PubMed  CAS  Google Scholar 

  • Gomez-Skarmeta JL, Modolell J (2002) Iroquois genes: genomic organization and function in vertebrate neural development. Curr Opin Genet Dev 12:403–408

    Article  PubMed  CAS  Google Scholar 

  • Guindon S, Gascuel O (2003) A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 52:696–704

    Article  PubMed  Google Scholar 

  • Guindon S, Lethiec F, Duroux P, Gascuel O (2005) PHYML Online—a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res 33:W557–W559

    Article  PubMed  CAS  Google Scholar 

  • Hobert O, Westphal H (2000) Functions of LIM-homeobox genes. Trends Genet 16:75–83

    Article  PubMed  CAS  Google Scholar 

  • Holland PWH, Garcia-Fernàndez J, Williams NA, Sidow A (1994) Gene duplications and the origins of vertebrate development. Development 1994 (Suppl):125–133

  • Holland ND, Holland LZ, Kozmik Z (1995) An amphioxus Pax gene, AmphiPax-1, expressed in embryonic endoderm, but not in mesoderm: implications for the evolution of class I paired box genes. Mol Mar Biol Biotechnol 4:206–214

    PubMed  CAS  Google Scholar 

  • Holland LZ, Schubert M, Kozmik Z, Holland ND (1999) AmphiPax-3/7, an amphioxus paired box gene: insights into chordate myogenesis, neurogenesis, and the possible evolutionary precursor of definitive vertebrate neural crest. Evol Dev 1:153–165

    Article  PubMed  CAS  Google Scholar 

  • Holland PWH, Booth HAF, Bruford EA (2007) Classification and nomenclature of all human homeobox genes. BMC Biol 5:47

    Article  PubMed  Google Scholar 

  • Holland LZ, Albalat R, Azumi K, Benito-Gutiérrez È, Blow MJ, Bronner-Fraser M, Brunet F, Butts T, Candiani S, Dishaw LJ, Ferrier DEK, Garcia-Fernàndez J, Gibson-Brown JJ, Gissi C, Godzik A, Hallböök F, Hirose D, Hosomichi K, Ikuta T, Inoko H, Kasahara M, Kasamatsu J, Kawashima T, Kimura A, Kobayashi M, Kozmik Z, Kubokawa K, Laudet V, Litman GW, McHardy AC, Meulemans D, Nonaka M, Olinski RP, Pancer Z, Pennacchio LA, Pestarino M, Rast JP, Rigoutsos I, Robinson-Rechavi M, Roch G, Saiga H, Sasakura Y, Satake M, Satou Y, Schubert M, Sherwood N, Shiina T, Takatori N, Tello J, Vopalensky P, Wada S, Xu A, Ye Y, Yoshida K, Yoshizaki F, Yu J-K, Zhang Q, Zmasek CM, Putnam NH, Rokhsar DS, Satoh N, Holland PWH (2008) The amphioxus genome illuminates vertebrate origins and cephalochordate biology. Genome Res 18:1100–1111

    Article  PubMed  CAS  Google Scholar 

  • Howard-Ashby M, Materna SC, Brown CT, Chen L, Cameron RA, Davidson EH (2006) Identification and characterisation of homeobox transcription factor genes in Strongylocentrotus purpuratus, and their expression in embryonic development. Dev Biol 300:74–89

    Article  PubMed  CAS  Google Scholar 

  • Irimia M, Maeso I, Garcia-Fernandez J (2008) Convergent evolution of clustering of Iroquois homeobox genes across metazoans. Mol Biol Evol 25:1521–1525

    Article  PubMed  CAS  Google Scholar 

  • Jackman WR, Langeland JA, Kimmel CB (2000) Islet reveals segmentation in the amphioxus hindbrain homolog. Dev Biol 220:16–26

    Article  PubMed  CAS  Google Scholar 

  • Kozmik Z, Holland ND, Kalousova A, Paces J, Schubert M, Holland LZ (1999) Characterization of an amphioxus paired box gene, AmphiPax2/5/8: developmental expression patterns in optic support cells, nephridium, thyroid-like structures and pharyngeal gill slits, but not in the midbrain-hindbrain boundary region. Development 126:1295–304

    PubMed  CAS  Google Scholar 

  • Kumar S, Tamura K, Nei M (2004) MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform 5:150–163

    Article  PubMed  CAS  Google Scholar 

  • Letunic I, Copley RR, Pils B, Pinkert S, Schultz J, Bork P (2006) SMART 5: domains in the context of genomes and networks. Nucleic Acids Res 34:D257–D260

    Article  PubMed  CAS  Google Scholar 

  • Liu M, Su M, Lyons GE, Bodmer R (2006) Functional conservation of zinc-finger homeodomain gene zfh1/SIP1 in Drosophila heart development. Dev Genes Evol 216:683–693

    Article  PubMed  CAS  Google Scholar 

  • Luke GN, Castro LFC, McLay K, Bird C, Coulson A, Holland PWH (2003) Dispersal of NK homeobox gene clusters in amphioxus and humans. Proc Natl Acad Sci USA 100:5292–5295

    Article  PubMed  CAS  Google Scholar 

  • Mesika A, Ben-Dor S, Laviad EL, Futerman AH (2007) A new functional motif in Hox domain-containing ceramide synthases: identification of a novel region flanking the Hox and TLC domains essential for activity. J Biol Chem 282:27366–27373

    Article  PubMed  CAS  Google Scholar 

  • Minguillón C, Ferrier DEK, Cebrián C, Garcia-Fernàndez J (2002) Gene duplications in the prototypical cephalochordate amphioxus. Gene 287:121–128

    Article  PubMed  Google Scholar 

  • Mizutani Y, Kihara A, Igarashi Y (2005) Mammalian Lass6 and its related family members regulate synthesis of specific ceramides. Biochem J 390:263–271

    Article  PubMed  CAS  Google Scholar 

  • Mukherjee K, Bürglin TR (2007) Comprehensive analysis of animal TALE homeobox genes: new conserved motifs and cases of accelerated evolution. J Mol Evol 65:137–153

    Article  PubMed  CAS  Google Scholar 

  • Nishijima I, Ohtoshi A (2006) Characterization of a novel prospero-related homeobox gene, Prox2. Mol Genet Genom 275:471–478

    Article  CAS  Google Scholar 

  • Ohshima K, Hattori M, Yada T, Gojobori T, Sakaki Y, Okada N (2003) Whole-genome screening indicates a possible burst of formation of processed pseudogenes and Alu repeats by particular L1 subfamilies in ancestral primates. Genome Biol 4:R74

    Article  PubMed  Google Scholar 

  • Peters T, Dildrop R, Ausmeier K, Ruther U (2000) Organization of mouse Iroquois homeobox genes in two clusters suggests a conserved regulation and function in vertebrate development. Genome Res 10:1453–1462

    Article  PubMed  CAS  Google Scholar 

  • Pewzner-Jung Y, Ben-Dor S, Futerman AH (2006) When do Lasses (longevity assurance genes) become CerS (ceramide synthases)?: insights into the regulation of ceramide synthesis. J Biol Chem 281:25001–25005

    Article  PubMed  CAS  Google Scholar 

  • Putnam NH, Butts T, Ferrier DEK, Furlong RF, Hellsten U, Kawashima T, Robinson-Rechavi M, Shoguchi E, Terry A, Yu J-K, Benito-Gutiérrez E, Dubchak I, Garcia-Fernàndez J, Grigoriev IV, Horton AC, de Jong PJ, Jurka J, Kapitonov V, Kohara Y, Kuroki Y, Lindquist E, Lucas S, Kazutoyo O, Pennacchio LA, Salamov AA, Satou Y, Sauka-Spengler T, Schmutz J, Shin I-T, Toyoda A, Gibson-Brown JJ, Bronner-Fraser M, Fujiyama A, Holland LZ, Holland PWH, Satoh N, Rokhsar DS (2008) The amphioxus genome and the evolution of the chordate karyotype. Nature 453:1064–1071

    Article  PubMed  CAS  Google Scholar 

  • Rambaut A (2006) FigTree. http://tree.bio.ed.ac.uk/software/figtree/

  • Ryan JF, Burton PM, Mazza ME, Kwong GK, Mullikin JC, Finnerty JF (2006) The cnidarian-bilaterian ancestor possessed at least 56 homeoboxes: evidence from the starlet sea anemone, Nematostella vectensis. Genome Biol 7:R64

    Article  PubMed  Google Scholar 

  • Schultz J, Milpetz F, Bork P, Ponting CP (1998) SMART, a simple modular architecture research tool: identification of signaling domains. Proc Natl Acad Sci USA 95:5857–5864

    Article  PubMed  CAS  Google Scholar 

  • Seo HC, Curtiss J, Mlodzik M, Fjose A (1999) Six class homeobox genes in Drosophila belong to three distinct families and are involved in head development. Mech Dev 83:127–139

    Article  PubMed  CAS  Google Scholar 

  • Takatori N, Saiga H (2008) Evolution of CUT class homeobox genes in the lineage to the vertebrates: insights from the analysis of the genomic sequence of amphioxus, Branchiostoma floridae. Int J Dev Biol 52(7) (in press). doi:10.1387/ijdb.072541nt

  • Wada S, Tokuoka M, Shoguchi E, Kobayashi K, Di Gregorio A, Spagnuolo A, Branno M, Kohara Y, Rokhsar D, Levine M et al (2003) A genomewide survey of developmentally relevant genes in Ciona intestinalis. II. Genes for homeobox transcription factors. Dev Genes Evol 213:222–234

    Article  PubMed  CAS  Google Scholar 

  • Yousef MS, Matthews BW (2005) Structural basis of Prospero-DNA interaction: implications for transcription regulation in developing cells. Structure 13:601–607

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgments

We thank Nik Putnam, Dan Rokshar, and other members of the amphioxus team at the Joint Genome Institute, Walnut Creek, California, for their very considerable efforts in determining the B. floridae genome sequence and making it available to the research community, and Linda Holland, Nori Satoh and Jeremy Gibson-Brown for important contributions to promotion and coordination of the project. We also thank members of the Developmental Program Laboratory, Department of Biological Sciences, Graduate School of Science and Engineering, Tokyo Metropolitan University, members of the J.W. Jenkinson Laboratory of Evolution and Development, Department of Zoology, University of Oxford, and Jordi Garcia-Fernàndez for discussions and advice. This work was supported by KAKENHI (Grant-in-Aid for Scientific Research) on Priority Area ‘Comparative genomics’ from the Ministry of Education, Culture, Sports, Science and Technology of Japan, Grant-in-Aid for Scientific Research B from JSPS (H.S.), MIUR Italy, FIRB 2001 BAU01WAFY and PRIN 2006 PRIN2006058952 (S.C., M.P.). and BBSRC (T.B, P.W.H.H., D.E.K.F.).

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Hidetoshi Saiga or Peter W. H. Holland.

Additional information

Communicated by: J. J. Gibson-Brown

N. Takatori and T. Butts contributed equally to this work.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Fig. 1

Arbitrarily rooted phylogenetic tree of ANTP class genes generated by maximum likelihood based on homeodomain sequences. Gene family support nodes are shown, except in cases where these are not recovered as monophyletic (Hox families, NK2, NK4, and Vent) due to uneven rates of sequence evolution and complex histories. The two Vent genes are represented by a single branch, since they have identical homeodomains. Amphioxus homeodomains identified in this study are marked with closed circles. The two B. floridae ro alleles differ in predicted amino acid sequence; Protein ID 290436 was used. Scale bar indicates evolutionary distance of 0.5 amino acid substitutions per position. The following notes apply to all supplementary figures: numbers at nodes indicate bootstrap support values given as a percentage of 500 bootstrap pseudoreplications of the data. Names on branches indicate protein names, with species name in abbreviated form and accession ID. Species codes are: Bb Branchiostoma belcheri; Bf Branchiostoma floridae; Ce Caenorhabditis elegans; Cf Canis familiaris; Ci Ciona intestinalis; Dm Drosophila melanogaster; Dr Danio rerio; Gg Gallus gallus; He Heliocidaris erythrogramma; Hr Halocynthia roretzi; Hs Homo sapiens; Mm Mus musculus; Nv Nematostella vectensis; Od Oikopleura dioica; Pd Platynereis dumerilii; Sk Saccoglossus kowalevskii; Sp Strongylocentrotus purpuratus;; Tf Takifugu rubripes; and Xl Xenopus laevis. (PDF 168 KB)

Supplementary Fig. 2

Arbitrarily rooted phylogenetic tree of Pax genes generated by maximum likelihood based on paired domain sequences. Scale bar indicates evolutionary distance of 0.05 amino acid substitutions per position. In this and all subsequent supplementary figures, amphioxus proteins identified in this study are represented without an accession ID and are marked with closed circles; amphioxus sequences previously reported are represented with accession ID and marked with open circles. The B. floridae genome Pax6 sequence used corresponds to Protein ID 291575. (GIF 109 KB)

High resolution image file (TIFF 28 MB)

Supplementary Fig. 3

Arbitrarily rooted phylogenetic tree of PRD class genes generated by maximum likelihood based on homeodomain sequences. Scale bar indicates evolutionary distance of 0.1 amino acid substitutions per position. (GIF 193 KB)

High resolution image file (TIFF 32 MB)

Supplementary Fig. 4

Arbitrarily rooted phylogenetic tree of LIM class genes generated by maximum likelihood based on homeodomain sequences. Scale bar indicates evolutionary distance of 0.1 amino acid substitutions per position. (GIF 122 KB)

High resolution image file (TIFF 31 MB)

Supplementary Fig. 5

Arbitrarily rooted phylogenetic tree of POU class genes generated by maximum likelihood based on homeodomain sequences. The highly divergent Hdx gene is not included. Scale bar indicates evolutionary distance of 0.1 amino acid substitutions per position. (GIF 73 KB)

High resolution image file (TIFF 30 MB)

Supplementary Fig. 6

Arbitrarily rooted phylogenetic tree of HNF class genes generated by maximum likelihood based on homeodomain sequences. The highly divergent Ahnf gene is not included. Scale bar indicates evolutionary distance of 0.1 amino acid substitutions per position. (GIF 58 KB)

High resolution image file (TIFF 26 MB)

Supplementary Fig. 7

Arbitrarily rooted phylogenetic tree of SINE class genes generated by maximum likelihood based on homeodomain sequences. Scale bar indicates evolutionary distance of 0.05 amino acid substitutions per position. (GIF 59 KB)

High resolution image file (TIFF 28 MB)

Supplementary Fig. 8

Arbitrarily rooted phylogenetic tree of TALE class genes generated by maximum likelihood based on homeodomain sequences. Scale bar indicates evolutionary distance of 0.1 amino acid substitutions per position. (GIF 103 KB)

High resolution image file (TIFF 31 MB)

Supplementary Fig. 9

Arbitrarily rooted phylogenetic tree of PROS class genes generated by maximum likelihood based on homeodomain sequences. Scale bar indicates evolutionary distance of 0.1 amino acid substitutions per position. (GIF 28 KB)

High resolution image file (TIFF 20 MB)

Supplementary Fig. 10

Arbitrarily rooted phylogenetic tree of homeodomains from ZF class genes generated by the maximum likelihood. Multiple homeodomains from a single gene are numbered according to their position from the N terminus. Scale bar indicates evolutionary distance of 0.2 amino acid substitutions per position. (GIF 156 KB)

High resolution image file (TIFF 33 MB)

Supplementary Fig. 11

Arbitrarily rooted phylogenetic tree of CERS class genes generated by maximum likelihood based on homeodomain sequences. Scale bar indicates evolutionary distance of 0.1 amino acid substitutions per position. (GIF 63 KB)

High resolution image file (TIFF 23 MB)

Supplementary Table 1

Chromosomal location and putative synteny of putative human homologues of amphioxus genes neighboring Pou3L. (DOC 38 KB)

Rights and permissions

Reprints and permissions

About this article

Cite this article

Takatori, N., Butts, T., Candiani, S. et al. Comprehensive survey and classification of homeobox genes in the genome of amphioxus, Branchiostoma floridae . Dev Genes Evol 218, 579–590 (2008). https://doi.org/10.1007/s00427-008-0245-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s00427-008-0245-9

Keywords

Navigation