Amphioxus functional genomics and the origins of vertebrate gene regulation.

Vertebrates have greatly elaborated the basic chordate body plan and evolved highly distinctive genomes that have been sculpted by two whole-genome duplications. Here we sequence the genome of the Mediterranean amphioxus (Branchiostoma lanceolatum) and characterize DNA methylation, chromatin accessibility, histone modifications and transcriptomes across multiple developmental stages and adult tissues to investigate the evolution of the regulation of the chordate genome. Comparisons with vertebrates identify an intermediate stage in the evolution of differentially methylated enhancers, and a high conservation of gene expression and its cis-regulatory logic between amphioxus and vertebrates that occurs maximally at an earlier mid-embryonic phylotypic period. We analyse regulatory evolution after whole-genome duplications, and find that-in vertebrates-over 80% of broadly expressed gene families with multiple paralogues derived from whole-genome duplications have members that restricted their ancestral expression, and underwent specialization rather than subfunctionalization. Counter-intuitively, paralogues that restricted their expression increased the complexity of their regulatory landscapes. These data pave the way for a better understanding of the regulatory principles that underlie key vertebrate innovations.

The genome of the deep-sea anemone Actinernus sp. contains a mega-array of ANTP-class homeobox genes.

Members of the phylum Cnidaria include sea anemones, corals and jellyfish, and have successfully colonized both marine and freshwater habitats throughout the world. The understanding of how cnidarians adapt to extreme environments such as the dark, high-pressure deep-sea habitat has been hindered by the lack of genomic information. Here, we report the first chromosome-level deep-sea cnidarian genome, of the anemone Actinernus sp., which was 1.39 Gbp in length and contained 44 970 gene models including 14 806 tRNA genes and 30 164 protein-coding genes. Analyses of homeobox genes revealed the longest chromosome hosts a mega-array of Hox cluster, HoxL, NK cluster and NKL homeobox genes; until now, such an array has only been hypothesized to have existed in ancient ancestral genomes. In addition to this striking arrangement of homeobox genes, analyses of microRNAs revealed cnidarian-specific complements that are distinctive for nested clades of these animals, presumably reflecting the progressive evolution of the gene regulatory networks in which they are embedded. Also, compared with other sea anemones, circadian rhythm genes were lost in Actinernus sp., which likely reflects adaptation to living in the dark. This high-quality genome of a deep-sea cnidarian thus reveals some of the likely molecular adaptations of this ecologically important group of metazoans to the extreme deep-sea environment. It also deepens our understanding of the evolution of genome content and organization of animals in general and cnidarians in particular, specifically from the viewpoint of key developmental control genes like the homeobox-encoding genes, where we find an array of genes that until now has only been hypothesized to have existed in the ancient ancestor that pre-dated both the cnidarians and bilaterians.

Amphioxus muscle transcriptomes reveal vertebrate-like myoblast fusion genes and a highly conserved role of insulin signalling in the metabolism of muscle.

BACKGROUND: The formation and functioning of muscles are fundamental aspects of animal biology, and the evolution of 'muscle genes' is central to our understanding of this tissue. Feeding-fasting-refeeding experiments have been widely used to assess muscle cellular and metabolic responses to nutrition. Though these studies have focused on vertebrate models and only a few invertebrate systems, they have found similar processes are involved in muscle degradation and maintenance. Motivation for these studies stems from interest in diseases whose pathologies involve muscle atrophy, a symptom also triggered by fasting, as well as commercial interest in the muscle mass of animals kept for consumption. Experimentally modelling atrophy by manipulating nutritional state causes muscle mass to be depleted during starvation and replenished with refeeding so that the genetic mechanisms controlling muscle growth and degradation can be understood. RESULTS: Using amphioxus, the earliest branching chordate lineage, we address the gap in previous work stemming from comparisons between distantly related vertebrate and invertebrate models. Our amphioxus feeding-fasting-refeeding muscle transcriptomes reveal a highly conserved myogenic program and that the pro-orthologues of many vertebrate myoblast fusion genes were present in the ancestral chordate, despite these invertebrate chordates having unfused mononucleate myocytes. We found that genes differentially expressed between fed and fasted amphioxus were orthologous to the genes that respond to nutritional state in vertebrates. This response is driven in a large part by the highly conserved IGF/Akt/FOXO pathway, where depleted nutrient levels result in activation of FOXO, a transcription factor with many autophagy-related gene targets. CONCLUSION: Reconstruction of these gene networks and pathways in amphioxus muscle provides a key point of comparison between the distantly related groups assessed thus far, significantly refining the reconstruction of the ancestral state for chordate myoblast fusion genes and identifying the extensive role of duplicated genes in the IGF/Akt/FOXO pathway across animals. Our study elucidates the evolutionary trajectory of muscle genes as they relate to the increased complexity of vertebrate muscles and muscle development.

Evolutionary diversification of the canonical Wnt signaling effector TCF/LEF in chordates.

Wnt signaling is essential during animal development and regeneration, but also plays an important role in diseases such as cancer and diabetes. The canonical Wnt signaling pathway is one of the most conserved signaling cascades in the animal kingdom, with the T-cell factor/lymphoid enhancer factor (TCF/LEF) proteins being the major mediators of Wnt/ß-catenin-regulated gene expression. In comparison with invertebrates, vertebrates possess a high diversity of TCF/LEF family genes, implicating this as a possible key change to Wnt signaling at the evolutionary origin of vertebrates. However, the precise nature of this diversification is only poorly understood. The aim of this study is to clarify orthology, paralogy, and isoform relationships within the TCF/LEF gene family within chordates via in silico comparative study of TCF/LEF gene structure, molecular phylogeny, and gene synteny. Our results support the notion that the four TCF/LEF paralog subfamilies in jawed vertebrates (gnathostomes) evolved via the two rounds of whole-genome duplications that occurred during early vertebrate evolution. Importantly, gene structure comparisons and synteny analysis of jawless vertebrate (cyclostome) TCFs suggest that a TCF7L2-like form of gene structure is a close proxy for the ancestral vertebrate structure. In conclusion, we propose a detailed evolutionary path based on a new pre-whole-genome duplication vertebrate TCF gene model. This ancestor gene model highlights the chordate and vertebrate innovations of TCF/LEF gene structure, providing the foundation for understanding the role of Wnt/ß-catenin signaling in vertebrate evolution.

More Than One-to-Four via 2R: Evidence of an Independent Amphioxus Expansion and Two-Gene Ancestral Vertebrate State for MyoD-Related Myogenic Regulatory Factors (MRFs).

The evolutionary transition from invertebrates to vertebrates involved extensive gene duplication, but understanding precisely how such duplications contributed to this transition requires more detailed knowledge of specific cases of genes and gene families. Myogenic differentiation (MyoD) has long been recognized as a master developmental control gene and member of the MyoD family of bHLH transcription factors (myogenic regulatory factors [MRFs]) that drive myogenesis across the bilaterians. Phylogenetic reconstructions within this gene family are complicated by multiple instances of gene duplication and loss in several lineages. Following two rounds of whole-genome duplication (2R WGD) at the origin of the vertebrates, the ancestral function of MRFs is thought to have become partitioned among the daughter genes, so that MyoD and Myf5 act early in myogenic determination, whereas Myog and Myf6 are expressed later, in differentiating myoblasts. Comparing chordate MRFs, we find an independent expansion of MRFs in the invertebrate chordate amphioxus, with evidence for a parallel instance of subfunctionalization relative to that of vertebrates. Conserved synteny between chordate MRF loci supports the 2R WGD events as a major force in shaping the evolution of vertebrate MRFs. We also resolve vertebrate MRF complements and organization, finding a new type of vertebrate MRF gene in the process, which allowed us to infer an ancestral two-gene state in the vertebrates corresponding to the early- and late-acting types of MRFs. This necessitates a revision of previous conclusions about the simple one-to-four origin of vertebrate MRFs.

Micro-RNA Clusters Integrate Evolutionary Constraints on Expression and Target Affinities: The miR-6/5/4/286/3/309 Cluster in Drosophila.

A striking feature of micro-RNAs is that they are often clustered in the genomes of animals. The functional and evolutionary consequences of this clustering remain obscure. Here, we investigated a micro-RNA cluster miR-6/5/4/286/3/309 that is conserved across drosophilid lineages. Small RNA sequencing revealed expression of this micro-RNA cluster in Drosophila melanogaster leg discs, and conditional overexpression of the whole cluster resulted in leg appendage shortening. Transgenic overexpression lines expressing different combinations of micro-RNA cluster members were also constructed. Expression of individual micro-RNAs from the cluster resulted in a normal wild-type phenotype, but either the expression of several ancient micro-RNAs together (miR-5/4/286/3/309) or more recently evolved clustered micro-RNAs (miR-6-1/2/3) can recapitulate the phenotypes generated by the whole-cluster overexpression. Screening of transgenic fly lines revealed downregulation of leg-patterning gene cassettes in generation of the leg-shortening phenotype. Furthermore, cell transfection with different combinations of micro-RNA cluster members revealed a suite of downstream genes targeted by all cluster members, as well as complements of targets that are unique for distinct micro-RNAs. Considered together, the micro-RNA targets and the evolutionary ages of each micro-RNA in the cluster demonstrate the importance of micro-RNA clustering, where new members can reinforce and modify the selection forces on both the cluster regulation and the gene regulatory network of existing micro-RNAs. Key words: micro-RNA, cluster, evolution.

Sampling the fish gill microbiome: a comparison of tissue biopsies and swabs.

BACKGROUND: Understanding the influence of methodology on results is an essential consideration in experimental design. In the expanding field of fish microbiology, many best practices and targeted techniques remain to be refined. This study aimed to compare microbial assemblages obtained from Atlantic salmon (Salmo salar) gills by swabbing versus biopsy excision. Results demonstrate the variation introduced by altered sampling strategies and enhance the available knowledge of the fish gill microbiome. RESULTS: The microbiome was sampled using swabs and biopsies from fish gills, with identical treatment of samples for 16S next generation Illumina sequencing. Results show a clear divergence in microbial communities obtained through the different sampling strategies, with swabbing consistently isolating a more diverse microbial consortia, and suffering less from the technical issue of host DNA contamination associated with biopsy use. Sequencing results from biopsy-derived extractions, however, hint at the potential for more cryptic localisation of some community members. CONCLUSIONS: Overall, results demonstrate a divergence in the obtained microbial community when different sampling methodology is used. Swabbing appears a superior method for sampling the microbiota of mucosal surfaces for broad ecological research in fish, whilst biopsies might be best applied in exploration of communities beyond the reach of swabs, such as sub-surface and intracellular microbes, as well as in pathogen diagnosis. Most studies on the external microbial communities of aquatic organisms utilise swabbing for sample collection, likely due to convenience. Much of the ultrastructure of gill tissue in live fish is, however, potentially inaccessible to swabbing, meaning swabbing might fail to capture the full diversity of gill microbiota. This work therefore also provides valuable insight into partitioning of the gill microbiota, informing varied applications of different sampling methods in experimental design for future research.

Reconstruction of ancient homeobox gene linkages inferred from a new high-quality assembly of the Hong Kong oyster (Magallana hongkongensis) genome.

BACKGROUND: Homeobox-containing genes encode crucial transcription factors involved in animal, plant and fungal development, and changes to homeobox genes have been linked to the evolution of novel body plans and morphologies. In animals, some homeobox genes are clustered together in the genome, either as remnants from ancestral genomic arrangements, or due to coordinated gene regulation. Consequently, analyses of homeobox gene organization across animal phylogeny provide important insights into the evolution of genome organization and developmental gene control, and their interaction. However, homeobox gene organization remains to be fully elucidated in several key animal ancestors, including those of molluscs, lophotrochozoans and bilaterians. RESULTS: Here, we present a high-quality chromosome-level genome assembly of the Hong Kong oyster, Magallana hongkongensis (2n = 20), for which 93.2% of the genomic sequences are contained on 10 pseudomolecules (~ 758 Mb, scaffold N50 = 72.3 Mb). Our genome assembly was scaffolded using Hi-C reads, facilitating a larger scaffold size compared to the recently published M. hongkongensis genome of Peng et al. (Mol Ecol Resources, 2020), which was scaffolded using the Crassostrea gigas assembly. A total of 46,963 predicted gene models (45,308 protein coding genes) were incorporated in our genome, and genome completeness estimated by BUSCO was 94.6%. Homeobox gene linkages were analysed in detail relative to available data for other mollusc lineages. CONCLUSIONS: The analyses performed in this study and the accompanying genome sequence provide important genetic resources for this economically and culturally valuable oyster species, and offer a platform to improve understanding of animal biology and evolution more generally. Transposable element content is comparable to that found in other mollusc species, contrary to the conclusion of another recent analysis. Also, our chromosome-level assembly allows the inference of ancient gene linkages (synteny) for the homeobox-containing genes, even though a number of the homeobox gene clusters, like the Hox/ParaHox clusters, are undergoing dispersal in molluscs such as this oyster.

Calcisponges have a ParaHox gene and dynamic expression of dispersed NK homeobox genes.

Sponges are simple animals with few cell types, but their genomes paradoxically contain a wide variety of developmental transcription factors, including homeobox genes belonging to the Antennapedia (ANTP) class, which in bilaterians encompass Hox, ParaHox and NK genes. In the genome of the demosponge Amphimedon queenslandica, no Hox or ParaHox genes are present, but NK genes are linked in a tight cluster similar to the NK clusters of bilaterians. It has been proposed that Hox and ParaHox genes originated from NK cluster genes after divergence of sponges from the lineage leading to cnidarians and bilaterians. On the other hand, synteny analysis lends support to the notion that the absence of Hox and ParaHox genes in Amphimedon is a result of secondary loss (the ghost locus hypothesis). Here we analysed complete suites of ANTP-class homeoboxes in two calcareous sponges, Sycon ciliatum and Leucosolenia complicata. Our phylogenetic analyses demonstrate that these calcisponges possess orthologues of bilaterian NK genes (Hex, Hmx and Msx), a varying number of additional NK genes and one ParaHox gene, Cdx. Despite the generation of scaffolds spanning multiple genes, we find no evidence of clustering of Sycon NK genes. All Sycon ANTP-class genes are developmentally expressed, with patterns suggesting their involvement in cell type specification in embryos and adults, metamorphosis and body plan patterning. These results demonstrate that ParaHox genes predate the origin of sponges, thus confirming the ghost locus hypothesis, and highlight the need to analyse the genomes of multiple sponge lineages to obtain a complete picture of the ancestral composition of the first animal genome.

Horizons in evolutionary genomics: an interview with David Ferrier.

David Ferrier is a Reader at the University of St Andrews and Deputy Director of the Scottish Oceans Institute, where his lab studies how the diversity of form in the animal kingdom evolved, with an emphasis on using comparative genomics. In this interview, David shares his thoughts on how to escape the 'straitjacket' of traditional model systems, transparency in peer review, and the past and future of genome sequencing.

Two more Posterior Hox genes and Hox cluster dispersal in echinoderms.

BACKGROUND: Hox genes are key elements in patterning animal development. They are renowned for their, often, clustered organisation in the genome, with supposed mechanistic links between the organisation of the genes and their expression. The widespread distribution and comparable functions of Hox genes across the animals has led to them being a major study system for comparing the molecular bases for construction and divergence of animal morphologies. Echinoderms (including sea urchins, sea stars, sea cucumbers, feather stars and brittle stars) possess one of the most unusual body plans in the animal kingdom with pronounced pentameral symmetry in the adults. Consequently, much interest has focused on their development, evolution and the role of the Hox genes in these processes. In this context, the organisation of echinoderm Hox gene clusters is distinctive. Within the classificatory system of Duboule, echinoderms constitute one of the clearest examples of Disorganized (D) clusters (i.e. intact clusters but with a gene order or orientation rearranged relative to the ancestral state). RESULTS: Here we describe two Hox genes (Hox11/13d and e) that have been overlooked in most previous work and have not been considered in reconstructions of echinoderm Hox complements and cluster organisation. The two genes are related to Posterior Hox genes and are present in all classes of echinoderm. Importantly, they do not reside in the Hox cluster of any species for which genomic linkage data is available. CONCLUSION: Incorporating the two neglected Posterior Hox genes into assessments of echinoderm Hox gene complements and organisation shows that these animals in fact have Split (S) Hox clusters rather than simply Disorganized (D) clusters within the Duboule classification scheme. This then has implications for how these genes are likely regulated, with them no longer covered by any potential long-range Hox cluster-wide, or multigenic sub-cluster, regulatory mechanisms.

Amphioxus SYCP1: a case of retrogene replacement and co-option of regulatory elements adjacent to the ParaHox cluster.

Retrogenes are formed when an mRNA is reverse-transcribed and reinserted into the genome in a location unrelated to the original locus. If this retrocopy inserts into a transcriptionally favourable locus and is able to carry out its original function, it can, in rare cases, lead to retrogene replacement. This involves the original, often multi-exonic, parental copy being lost whilst the newer single-exon retrogene copy 'replaces' the role of the ancestral parent gene. One example of this is amphioxus SYCP1, a gene that encodes a protein used in synaptonemal complex formation during meiosis and which offers the opportunity to examine how a retrogene evolves after the retrogene replacement event. SYCP1 genes exist as large multi-exonic genes in most animals. AmphiSYCP1, however, contains a single coding exon of ~ 3200 bp and has inserted next to the ParaHox cluster of amphioxus, whilst the multi-exonic ancestral parental copy has been lost. Here, we show that AmphiSYCP1 has not only replaced its parental copy, but also has evolved additional regulatory function by co-opting a bidirectional promoter from the nearby AmphiCHIC gene. AmphiSYCP1 has also evolved a de novo, multi-exonic 5'untranslated region that displays distinct regulatory states, in the form of two different isoforms, and has evolved novel expression patterns during amphioxus embryogenesis in addition to its ancestral role in meiosis. The absence of ParaHox-like expression of AmphiSYCP1, despite its proximity to the ParaHox cluster, also suggests that this gene is not influenced by any potential pan-cluster regulatory mechanisms, which are seemingly restricted to only the ParaHox genes themselves.

Light-sheet microscopy using an Airy beam.

Light-sheet microscopy facilitates rapid, high-contrast, volumetric imaging with minimal sample exposure. However, the rapid divergence of a traditional Gaussian light sheet restricts the field of view (FOV) that provides innate subcellular resolution. We show that the Airy beam innately yields high contrast and resolution up to a tenfold larger FOV. In contrast to the Bessel beam, which also provides an increased FOV, the Airy beam's characteristic asymmetric excitation pattern results in all fluorescence contributing positively to the contrast, enabling a step change for light-sheet microscopy.

The first myriapod genome sequence reveals conservative arthropod gene content and genome organisation in the centipede Strigamia maritima.

Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history.

TCF/Lef regulates the Gsx ParaHox gene in central nervous system development in chordates.

BACKGROUND: The ParaHox genes play an integral role in the anterior-posterior (A-P) patterning of the nervous system and gut of most animals. The ParaHox cluster is an ideal system in which to study the evolution and regulation of developmental genes and gene clusters, as it displays similar regulatory phenomena to its sister cluster, the Hox cluster, but offers a much simpler system with only three genes. RESULTS: Using Ciona intestinalis transgenics, we isolated a regulatory element upstream of Branchiostoma floridae Gsx that drives expression within the central nervous system of Ciona embryos. The minimal amphioxus enhancer region required to drive CNS expression has been identified, along with surrounding sequence that increases the efficiency of reporter expression throughout the Ciona CNS. TCF/Lef binding sites were identified and mutagenized and found to be required to drive the CNS expression. Also, individual contributions of TCF/Lef sites varied across the regulatory region, revealing a partial division of function across the Bf-Gsx-Up regulatory element. Finally, when all TCF/Lef binding sites are mutated CNS expression is not only abolished, but a latent repressive function is also unmasked. CONCLUSIONS: We have identified a B. floridae Gsx upstream regulatory element that drives CNS expression within transgenic Ciona intestinalis, and have shown that this CNS expression is dependent upon TCF/Lef binding sites. We examine the evolutionary and developmental implications of these results, and discuss the possibility of TCF/Lef not only as a regulator of chordate Gsx, but as a deeply conserved regulatory factor controlling all three ParaHox genes across the Metazoa.

Another biomineralising protostome with an msp130 gene and conservation of msp130 gene structure across Bilateria.

Msp130 genes are known for their association with biomineralisation, principally in echinoderm skeletogenesis. Recently, msp130 genes were shown to exist more widely across the animal kingdom, including in molluscs, and a hypothesis was formed that the genes had arisen independently in the deuterostome and mollusc lineages via horizontal gene transfer, thus facilitating the evolution of biomineralisation in these distinct lineages (Ettensohn, 2014). Here we show that another biomineralising protostome, the polychaete Spirobranchus (formerly Pomatoceros) lamarcki also possesses an msp130 gene, and expresses it during a biomineralisation process. However, based on analysis of gene structure, we hypothesize that the protostome and deuterostome msp130 genes did not originate from independent horizontal gene transfers, but instead are descended from a gene already present in the bilaterian ancestor, with the gene being secondarily lost from several lineages.

Evolutionary crossroads in developmental biology: annelids.

Annelids (the segmented worms) have a long history in studies of animal developmental biology, particularly with regards to their cleavage patterns during early development and their neurobiology. With the relatively recent reorganisation of the phylogeny of the animal kingdom, and the distinction of the super-phyla Ecdysozoa and Lophotrochozoa, an extra stimulus for studying this phylum has arisen. As one of the major phyla within Lophotrochozoa, Annelida are playing an important role in deducing the developmental biology of the last common ancestor of the protostomes and deuterostomes, an animal from which >98% of all described animal species evolved.

The amphioxus genome and the evolution of the chordate karyotype.

Lancelets ('amphioxus') are the modern survivors of an ancient chordate lineage, with a fossil record dating back to the Cambrian period. Here we describe the structure and gene content of the highly polymorphic approximately 520-megabase genome of the Florida lancelet Branchiostoma floridae, and analyse it in the context of chordate evolution. Whole-genome comparisons illuminate the murky relationships among the three chordate groups (tunicates, lancelets and vertebrates), and allow not only reconstruction of the gene complement of the last common chordate ancestor but also partial reconstruction of its genomic organization, as well as a description of two genome-wide duplications and subsequent reorganizations in the vertebrate lineage. These genome-scale events shaped the vertebrate genome and provided additional genetic variation for exploitation during vertebrate evolution.

Time is of the essence for ParaHox homeobox gene clustering.

ParaHox genes, and their evolutionary sisters the Hox genes, are integral to patterning the anterior-posterior axis of most animals. Like the Hox genes, ParaHox genes can be clustered and exhibit the phenomenon of colinearity - gene order within the cluster matching gene activation. Two new instances of ParaHox clustering provide the first examples of intact clusters outside chordates, with gene expression lending weight to the argument that temporal colinearity is the key to understanding clustering.

Extensive chordate and annelid macrosynteny reveals ancestral homeobox gene organization.

Genes with the homeobox motif are crucial in developmental biology and widely implicated in the evolution of development. The Antennapedia (ANTP)-class is one of the two major classes of animal homeobox genes, and includes the Hox genes, renowned for their role in patterning the anterior-posterior axis of animals. The origin and evolution of the ANTP-class genes are a matter of some debate. A principal guiding hypothesis has been the existence of an ancient gene Mega-cluster deep in animal ancestry. This hypothesis was largely established from linkage data from chordates, and the Mega-cluster hypothesis remains to be seriously tested in protostomes. We have thus mapped ANTP-class homeobox genes to the chromosome level in a lophotrochozoan protostome. Our comparison of gene organization in Platynereis dumerilii and chordates indicates that the Mega-cluster, if it did exist, had already been broken up onto four chromosomes by the time of the protostome-deuterostome ancestor (PDA). These results not only elucidate an aspect of the genome organization of the PDA but also reveal high levels of macrosynteny between P. dumerilii and chordates. This implies a very low rate of interchromosomal genome rearrangement in the lineages leading to P. dumerilii and the chordate ancestor since the time of the PDA.