Archetypal organization of the amphioxus Hox gene cluster

40 years of the homeobox: mechanisms of Hox spatial-temporal collinearity in vertebrates

Animal body plans are established during embryonic development by the Hox genes. This patterning process relies on the differential expression of Hox genes along the head-to-tail axis. Hox spatial collinearity refers to the relationship between the organization of Hox genes in clusters and the differential Hox expression, whereby the relative order of the Hox genes within a cluster mirrors the spatial sequence of expression in the developing embryo. In vertebrates, the cluster organization is also associated with the timing of Hox activation, which harmonizes Hox expression with the progressive emergence of axial tissues. Thereby, in vertebrates, Hox temporal collinearity is intimately linked to Hox spatial collinearity. Understanding the mechanisms contributing to Hox temporal and spatial collinearity is thus key to the comprehension of vertebrate patterning. Here, we provide an overview of the main discoveries pertaining to the mechanisms of Hox spatial-temporal collinearity.

3D reconstruction of a gastrulating human embryo

Progress in understanding early human development has been impeded by the scarcity of reference datasets from natural embryos, particularly those with spatial information during crucial stages like gastrulation. We conducted high-resolution spatial transcriptomics profiling on 38,562 spots from 62 transverse sections of an intact Carnegie stage (CS) 8 human embryo. From this spatial transcriptomic dataset, we constructed a 3D model of the CS8 embryo, in which a range of cell subtypes are identified, based on gene expression patterns and positional register, along the anterior-posterior, medial-lateral, and dorsal-ventral axis in the embryo. We further characterized the lineage trajectories of embryonic and extra-embryonic tissues and associated regulons and the regionalization of signaling centers and signaling activities that underpin lineage progression and tissue patterning during gastrulation. Collectively, the findings of this study provide insights into gastrulation and post-gastrulation development of the human embryo.

Hox, homology, and parsimony: An organismal perspective

Hox genes are important regulators in animal development. They often show a mosaic of conserved (e.g., longitudinal axis patterning) and lineage-specific novel functions (e.g., development of skeletal, sensory, or locomotory systems). Despite extensive research over the past decades, it remains controversial at which node in the animal tree of life the Hox cluster evolved. Its presence already in the last common metazoan ancestor has been proposed, although the genomes of both putative earliest extant metazoan offshoots, the ctenophores and the poriferans, are devoid of Hox sequences. The lack of Hox genes in the supposedly "simple"-built poriferans and their low number in cnidarians and the basally branching bilaterians, the xenacoelomorphs, seems to support the classical notion that the number of Hox genes is correlated with the degree of animal complexity. However, the 4-fold increase of the Hox cluster in xiphosurans, a basally branching chelicerate clade, as well as the situation in some teleost fishes that show a multitude of Hox genes compared to, e.g., human, demonstrates, that there is no per se direct correlation between organismal complexity and Hox number. Traditional approaches have tried to base homology on the morphological level on shared expression profiles of individual genes, but recent data have shown that, in particular with respect to Hox and other regulatory genes, complex gene-gene interactions rather than expression signatures of individual genes alone are responsible for shaping morphological traits during ontogeny. Accordingly, for sound homology assessments and reconstructions of character evolution on organ system level, additional independent datasets (e.g., morphological, developmental) need to be included in any such analyses. If supported by solid data, proposed structural homology should be regarded as valid and not be rejected solely on the grounds of non-parsimonious distribution of the character over a given phylogenetic topology.

Spiralian genomics and the evolution of animal genome architecture

Recent developments in sequencing technologies have greatly improved our knowledge of phylogenetic relationships and genomic architectures throughout the tree of life. Spiralia, a diverse clade within Protostomia, is essential for understanding the evolutionary history of parasitism, gene conversion, nervous systems and animal body plans. In this review, we focus on the current hypotheses of spiralian phylogeny and investigate the impact of long-read sequencing on the quality of genome assemblies. We examine chromosome-level assemblies to highlight key genomic features that have driven spiralian evolution, including karyotype, synteny and the Hox gene organization. In addition, we show how chromosome rearrangement has influenced spiralian genomic structures. Although spiralian genomes have undergone substantial changes, they exhibit both conserved and lineage-specific features. We recommend increasing sequencing efforts and expanding functional genomics research to deepen insights into spiralian biology.

Sequential and directional insulation by conserved CTCF sites underlies the Hox timer in stembryos

During development, Hox genes are temporally activated according to their relative positions on their clusters, contributing to the proper identities of structures along the rostrocaudal axis. To understand the mechanism underlying this Hox timer, we used mouse embryonic stem cell-derived stembryos. Following Wnt signaling, the process involves transcriptional initiation at the anterior part of the cluster and a concomitant loading of cohesin complexes enriched on the transcribed DNA segments, that is, with an asymmetric distribution favoring the anterior part of the cluster. Chromatin extrusion then occurs with successively more posterior CTCF sites acting as transient insulators, thus generating a progressive time delay in the activation of more posterior-located genes due to long-range contacts with a flanking topologically associating domain. Mutant stembryos support this model and reveal that the presence of evolutionary conserved and regularly spaced intergenic CTCF sites controls the precision and the pace of this temporal mechanism.

Chromosome-level genome and the identification of sex chromosomes in Uloborus diversus

The orb web is a remarkable example of animal architecture that is observed in families of spiders that diverged over 200 million years ago. While several genomes exist for araneid orb-weavers, none exist for other orb-weaving families, hampering efforts to investigate the genetic basis of this complex behavior. Here we present a chromosome-level genome assembly for the cribellate orb-weaving spider Uloborus diversus. The assembly reinforces evidence of an ancient arachnid genome duplication and identifies complete open reading frames for every class of spidroin gene, which encode the proteins that are the key structural components of spider silks. We identified the 2 X chromosomes for U. diversus and identify candidate sex-determining loci. This chromosome-level assembly will be a valuable resource for evolutionary research into the origins of orb-weaving, spidroin evolution, chromosomal rearrangement, and chromosomal sex determination in spiders.

A new deep learning technique reveals the exclusive functional contributions of individual cancer mutations

Cancers are caused by genomic alterations that may be inherited, induced by environmental carcinogens, or caused due to random replication errors. Postinduction of carcinogenicity, mutations further propagate and drastically alter the cancer genomes. Although a subset of driver mutations has been identified and characterized to date, most cancer-related somatic mutations are indistinguishable from germline variants or other noncancerous somatic mutations. Thus, such overlap impedes appreciation of many deleterious but previously uncharacterized somatic mutations. The major bottleneck arises due to patient-to-patient variability in mutational profiles, making it difficult to associate specific mutations with a given disease outcome. Here, we describe a newly developed technique Continuous Representation of Codon Switches (CRCS), a deep learning-based method that allows us to generate numerical vector representations of mutations, thereby enabling numerous machine learning-based tasks. We demonstrate three major applications of CRCS; first, we show how CRCS can help detect cancer-related somatic mutations in the absence of matched normal samples, which has applications in cell-free DNA–based assessment of tumor mutation burden. Second, the proposed approach also enables identification and exploration of driver genes; our analyses implicate DMD, RSK4, OFD1, WDR44, and AFF2 as potential cancer drivers. Finally, we used CRCS to score individual mutations in a tumor sample, which was found to be predictive of patient survival in bladder urothelial carcinoma, hepatocellular carcinoma, and lung adenocarcinoma. Taken together, we propose CRCS as a valuable computational tool for analysis of the functional significance of individual cancer mutations.

Diversification and Functional Evolution of HOX Proteins

Gene duplication and divergence is a major contributor to the generation of morphological diversity and the emergence of novel features in vertebrates during evolution. The availability of sequenced genomes has facilitated our understanding of the evolution of genes and regulatory elements. However, progress in understanding conservation and divergence in the function of proteins has been slow and mainly assessed by comparing protein sequences in combination with in vitro analyses. These approaches help to classify proteins into different families and sub-families, such as distinct types of transcription factors, but how protein function varies within a gene family is less well understood. Some studies have explored the functional evolution of closely related proteins and important insights have begun to emerge. In this review, we will provide a general overview of gene duplication and functional divergence and then focus on the functional evolution of HOX proteins to illustrate evolutionary changes underlying diversification and their role in animal evolution.

Regulation of posterior Hox genes by sex steroids explains vertebral variation in inbred mouse strains

A series of elegant embryo transfer experiments in the 1950s demonstrated that the uterine environment could alter vertebral patterning in inbred mouse strains. In the intervening decades, attention has tended to focus on the technical achievements involved and neglected the underlying biological question: how can genetically homogenous individuals have a heterogenous number of vertebrae? Here I revisit these experiments and, with the benefit of knowledge of the molecular-level processes of vertebral patterning gained over the intervening decades, suggest a novel hypothesis for homeotic transformation of the last lumbar vertebra to the adjacent sacral type through regulation of Hox genes by sex steroids. Hox genes are involved in both axial patterning and development of male and female reproductive systems and have been shown to be sensitive to sex steroids in vitro and in vivo. Regulation of these genes by sex steroids and resulting alterations to vertebral patterning may hint at a deep evolutionary link between the ribless lumbar region of mammals and the switch from egg-laying to embryo implantation. An appreciation of the impact of sex steroids on Hox genes may explain some puzzling aspects of human disease, and highlights the spine as a neglected target for in utero exposure to endocrine disruptors.

Essay the (unusual) heuristic value of Hox gene clusters; a matter of time?

Ever since their first report in 1984, Antennapedia-type homeobox (Hox) genes have been involved in such a series of interesting observations, in particular due to their conserved clustered organization between vertebrates and arthropods, that one may legitimately wonder about the origin of this heuristic value. In this essay, I first consider different examples where Hox gene clusters have been instrumental in providing conceptual advances, taken from various fields of research and mostly involving vertebrate embryos. These examples touch upon our understanding of genomic evolution, the revisiting of 19th century views on the relationships between development and evolution and the building of a new framework to understand long-range and pleiotropic gene regulation during development. I then discuss whether the high value of the Hox gene family, when considered as an epistemic object, is related to its clustered structure (and the absence thereof in some animal species) and, if so, what is it in such particular genetic oddities that made them so generous in providing the scientific community with interesting information.

Evolution and expression analysis of STAT family members in blunt snout bream (Megalobrama amblycephala)

The Janus kinase/signal transducer and activator of transcription (JAK/STAT) pathway is involved in regulating the body's immunity, cell proliferation, differentiation, and apoptosis. Members of the STAT family have been extensively studied in different mammalian species. However, there are few studies on the STAT family genes in farmed economic fish. In this study, eight STAT genes including STAT1a, STAT1b, STAT2, STAT3, STAT4, STAT5a, STAT5b and STAT6, in blunt snout bream (Megalobrama amblycephala), an economically important fish in China, were identified and characterized. Analyses of gene location, phylogeny and conserved synteny were conducted to infer the evolutionary origin of these STAT family genes. Furthermore, the evolutionary origin model of STATs was constructed based on the 2R hypothesis and teleost genome duplication (TGD) hypothesis, which clarified the evolutionary origin of the eight STATs in blunt snout bream. Besides, expression of the eight STATs was detected in 10 tissues of healthy blunt snout bream, which showed different expression patterns, and all had the highest level in the blood. In addition, expression of the STATs was significantly induced in the spleen, liver, and kidney after infection of Aeromonas hydrophila, suggesting that they play an important role in protecting the host from pathogens. In general, the evolution of cytokine-related genes parallels that of the immune system, which has likely been a main evolutionary driver. Therefore, the evolutionary model of STAT genes, constructed in the non-model organism pioneeringly, may provide some enlightenment for the evolution of the fish STAT family genes and their involvement in the immune function.

Gene cluster from plant to microbes: Their role in genome architecture, organism's development, specialized metabolism and drug discovery

Plants and microbes fulfil our daily requirements through different high-value chemicals, e.g., nutraceuticals, pharmaceuticals, cosmetics, and through varieties of fruits, crops, vegetables, and many more. Utmost care would therefore be taken for growth, development and sustainability of these important crops and medicinal plants and microbes. Homeobox genes and HOX clusters and their recently characterized expanded family members, including newly discovered homeobox, WOX gene from medicinal herb, Panax ginseng, significantly contributes in the growth and development of these organisms. On the other hand, secondary metabolites produced through secondary metabolism of plants and microbes are used as organisms defense as well as drugs/drug-like molecules for humans. Both the developmental HOX cluster and the biosynthetic gene-cluster (BGC) for secondary metabolites are organised in organisms genome. Genome mining and genomewide analysis of these clusters will definitely identify and characterize many more important molecules from unexplored plants and microbes and underexplored human microbiota and the evolution studies of these clusters will indicate their source of origin. Although genomics revolution now continues at a pace, till date only few hundred plant genome sequences are available. However, next-generation sequencing (NGS) technology now in market and may be applied even for plants with recalcitrant genomes, eventually may discover genomic potential towards production of secondary metabolites of diverse plants and micro-organisms present in the environment and microbiota. Additionally, the development of tools for genome mining e.g., antiSMASH, plantiSMASH, and more and more computational approaches that predicts hundreds of secondary metabolite BGCs will be discussed.

Evolution and genomic organization of the insect sHSP gene cluster and coordinate regulation in phenotypic plasticity

Background

Conserved syntenic gene complexes are rare in Arthropods and likely only retained due to functional constraint. Numerous sHSPs have been identified in the genomes of insects, some of which are located clustered in close proximity. Previous phylogenetic analyses of these clustered sHSP have been limited to a small number of holometabolous insect species and have not determined the pattern of evolution of the clustered sHSP genes (sHSP-C) in insect or Arthropod lineages.

Results

Using eight genomes from representative insect orders and three non-insect arthropod genomes we have identified that a syntenic cluster of sHSPs (sHSP-C) is a hallmark of most Arthropod genomes. Using 11 genomes from Hymenopteran species our phylogenetic analyses have refined the evolution of the sHSP-C in Hymenoptera and found that the sHSP-C is order-specific with evidence of birth-and-death evolution in the hymenopteran lineage. Finally we have shown that the honeybee sHSP-C is co-ordinately expressed and is marked by genomic features, including H3K27me3 histone marks consistent with coordinate regulation, during honeybee ovary activation.

Conclusions

The syntenic sHSP-C is present in most insect genomes, and its conserved coordinate expression and regulation implies that it is an integral genomic component of environmental response in arthropods.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01885-8.

A Preliminary Single-Cell RNA-Seq Analysis of Embryonic Cells That Express Brachyury in the Amphioxus, Branchiostoma japonicum

Among chordate taxa, the cephalochordates diverged earlier than urochordates and vertebrates; thus, they retain unique, primitive developmental features. In particular, the amphioxus notochord has muscle-like properties, a feature not seen in urochordates or vertebrates. Amphioxus contains two Brachyury genes, Bra1 and Bra2. Bra2 is reportedly expressed in the blastopore, notochord, somites, and tail bud, in contrast to a low level of Bra1 expression only in notochord. To distinguish the expression profiles of the two Brachyury genes at the single-cell level, we carried out single-cell RNA-seq (scRNA-seq) analysis using the amphioxus, Branchiostoma japonicum. This scRNA-seq analysis classified B. japonicum embryonic cells into 15 clusters at developmental stages from midgastrula to early swimming larva. Brachyury was expressed in cells of clusters 4, 5, 8, and 9. We first confirmed that cluster 8 comprises cells that form somites since this cluster specifically expresses four myogenic factor genes. Cluster 9 contains a larger number of cells with high levels of Bra2 expression and a smaller number of cells with Bra1 expression. Simultaneous expression in cluster 9 of tool-kit genes, including FoxA, Goosecoid, and hedgehog, showed that this cluster comprises cells that form the notochord. Expression of Bra2, but not Bra1, in cells of clusters 4 and 5 at the gastrula stage together with expression of Wnt1 and Caudal indicates that clusters 4 and 5 comprise cells of the blastopore, which contiguously form the tail bud. In addition, Hox1, Hox3, and Hox4 were highly expressed in Bra2-expressing clusters 4, 5, 8, and 9 in a temporally coordinated manner, suggesting roles of anterior Hox genes in specification of mesodermal organs, including somites, notochord, and tail bud. This scRNA-seq analysis therefore highlights differences between the two Brachyury genes in relation to embryonic regions in which they are expressed and their levels of expression. Bra2 is the ancestral Brachyury in amphioxus, since expression in the blastopore is shared with other deuterostomes. On the other hand, Bra1 is a duplicate copy and likely evolved a supplementary function in notochord and somite formation in the Branchiostoma lineage.

Improved Understanding of the Role of Gene and Genome Duplications in Chordate Evolution With New Genome and Transcriptome Sequences

Comparative approaches to understanding chordate genomes have uncovered a significant role for gene duplications, including whole genome duplications (WGDs), giving rise to and expanding gene families. In developmental biology, gene families created and expanded by both tandem and WGDs are paramount. These genes, often involved in transcription and signalling, are candidates for underpinning major evolutionary transitions because they are particularly prone to retention and subfunctionalisation, neofunctionalisation, or specialisation following duplication. Under the subfunctionalisation model, duplication lays the foundation for the diversification of paralogues, especially in the context of gene regulation. Tandemly duplicated paralogues reside in the same regulatory environment, which may constrain them and result in a gene cluster with closely linked but subtly different expression patterns and functions. Ohnologues (WGD paralogues) often diversify by partitioning their expression domains between retained paralogues, amidst the many changes in the genome during rediploidisation, including chromosomal rearrangements and extensive gene losses. The patterns of these retentions and losses are still not fully understood, nor is the full extent of the impact of gene duplication on chordate evolution. The growing number of sequencing projects, genomic resources, transcriptomics, and improvements to genome assemblies for diverse chordates from non-model and under-sampled lineages like the coelacanth, as well as key lineages, such as amphioxus and lamprey, has allowed more informative comparisons within developmental gene families as well as revealing the extent of conserved synteny across whole genomes. This influx of data provides the tools necessary for phylogenetically informed comparative genomics, which will bring us closer to understanding the evolution of chordate body plan diversity and the changes underpinning the origin and diversification of vertebrates.

Analysis of the lamprey genotype provides insights into caspase evolution and functional divergence

The cysteine-containing aspartate specific proteinase (caspase) family plays important roles in apoptosis and the maintenance of homeostasis in lampreys. We conducted genomic and functional comparisons of six distinct lamprey caspase groups with human counterparts to determine how these expanded molecules evolved to adapt to the changing caspase-mediated signaling pathways. Our results showed that lineage-specific duplication and rearrangement were responsible for expanding lamprey caspases 3 and 7, whereas caspases 1, 6, 8, and 9 maintained a relatively stable genome and protein structure. Lamprey caspase family molecules displayed various expression patterns and were involved in the innate immune response. Caspase 1 and 7 functioned as a pattern recognition receptor with a broad-spectrum of microbial recognition and bactericidal effect. Additionally, caspases 1 and 7 may induce cell apoptosis in a time- and dose-dependent manner; however, apoptosis was inhibited by caspase inhibitors. Thus, these molecules may reflect the original state of the vertebrates caspase family. Our phylogenetic and functional data provide insights into the evolutionary history of caspases and illustrate their functional characteristics in primitive vertebrates.

Deeply conserved synteny resolves early events in vertebrate evolution

Although it is widely believed that early vertebrate evolution was shaped by ancient whole-genome duplications, the number, timing and mechanism of these events remain elusive. Here, we infer the history of vertebrates through genomic comparisons with a new chromosome-scale sequence of the invertebrate chordate amphioxus. We show how the karyotypes of amphioxus and diverse vertebrates are derived from 17 ancestral chordate linkage groups (and 19 ancestral bilaterian groups) by fusion, rearrangement and duplication. We resolve two distinct ancient duplications based on patterns of chromosomal conserved synteny. All extant vertebrates share the first duplication, which occurred in the mid/late Cambrian by autotetraploidization (that is, direct genome doubling). In contrast, the second duplication is found only in jawed vertebrates and occurred in the mid-late Ordovician by allotetraploidization (that is, genome duplication following interspecific hybridization) from two now-extinct progenitors. This complex genomic history parallels the diversification of vertebrate lineages in the fossil record.

In vivo Hox binding specificity revealed by systematic changes to a single cis regulatory module

Hox proteins belong to a family of transcription factors with similar DNA binding specificities that control animal differentiation along the antero-posterior body axis. Hox proteins are expressed in partially overlapping regions where each one is responsible for the formation of particular organs and structures through the regulation of specific direct downstream targets. Thus, explaining how each Hox protein can selectively control its direct targets from those of another Hox protein is fundamental to understand animal development. Here we analyse a cis regulatory module directly regulated by seven different Drosophila Hox proteins and uncover how different Hox class proteins differentially control its expression. We find that regulation by one or another Hox protein depends on the combination of three modes: Hox-cofactor dependent DNA-binding specificity; Hox-monomer binding sites; and interaction with positive and negative Hox-collaborator proteins. Additionally, we find that similar regulation can be achieved by Amphioxus orthologs, suggesting these three mechanisms are conserved from insects to chordates.

Hox proteins are expressed in partially overlapping regions to inform development along the embryo’s head-tail axis. Here the authors analyse a cis regulatory module directly regulated by seven different Drosophila Hox proteins to uncover how different Hox class proteins differentially control its expression.

Evolution of vertebrate spinal cord patterning

The vertebrate spinal cord is organized across three developmental axes, anterior-posterior (AP), dorsal-ventral (DV), and medial-lateral (ML). Patterning of these axes is regulated by canonical intercellular signaling pathways: the AP axis by Wnt, fibroblast growth factor, and retinoic acid (RA), the DV axis by Hedgehog, Tgfβ, and Wnt, and the ML axis where proliferation is controlled by Notch. Developmental time plays an important role in which signal does what and when. Patterning across the three axes is not independent, but linked by interactions between signaling pathway components and their transcriptional targets. Combined this builds a sophisticated organ with many different types of cell in specific AP, DV, and ML positions. Two living lineages share phylum Chordata with vertebrates, amphioxus, and tunicates, while the jawless fish such as lampreys, survive as the most basally divergent vertebrate lineage. Genes and mechanisms shared between lampreys and other vertebrates tell us what predated vertebrates, while those also shared with other chordates tell us what evolved early in chordate evolution. Between these lie vertebrate innovations: genetic and developmental changes linked to evolution of new morphology. These include gene duplications, differences in how signals are received, and new regulatory connections between signaling pathways and their target genes.

Evolutionary transition between invertebrates and vertebrates via methylation reprogramming in embryogenesis

Major evolutionary transitions are enigmas, and the most notable enigma is between invertebrates and vertebrates, with numerous spectacular innovations. To search for the molecular connections involved, we asked whether global epigenetic changes may offer a clue by surveying the inheritance and reprogramming of parental DNA methylation across metazoans. We focused on gametes and early embryos, where the methylomes are known to evolve divergently between fish and mammals. Here, we find that methylome reprogramming during embryogenesis occurs neither in pre-bilaterians such as cnidarians nor in protostomes such as insects, but clearly presents in deuterostomes such as echinoderms and invertebrate chordates, and then becomes more evident in vertebrates. Functional association analysis suggests that DNA methylation reprogramming is associated with development, reproduction and adaptive immunity for vertebrates, but not for invertebrates. Interestingly, the single HOX cluster of invertebrates maintains unmethylated status in all stages examined. In contrast, the multiple HOX clusters show dramatic dynamics of DNA methylation during vertebrate embryogenesis. Notably, the methylation dynamics of HOX clusters are associated with their spatiotemporal expression in mammals. Our study reveals that DNA methylation reprogramming has evolved dramatically during animal evolution, especially after the evolutionary transitions from invertebrates to vertebrates, and then to mammals.

Hox genes: Downstream "effectors" of retinoic acid signaling in vertebrate embryogenesis

One of the major regulatory challenges of animal development is to precisely coordinate in space and time the formation, specification, and patterning of cells that underlie elaboration of the basic body plan. How does the vertebrate plan for the nervous and hematopoietic systems, heart, limbs, digestive, and reproductive organs derive from seemingly similar population of cells? These systems are initially established and patterned along the anteroposterior axis (AP) by opposing signaling gradients that lead to the activation of gene regulatory networks involved in axial specification, including the Hox genes. The retinoid signaling pathway is one of the key signaling gradients coupled to the establishment of axial patterning. The nested domains of Hox gene expression, which provide a combinatorial code for axial patterning, arise in part through a differential response to retinoic acid (RA) diffusing from anabolic centers established within the embryo during development. Hence, Hox genes are important direct effectors of retinoid signaling in embryogenesis. This review focuses on describing current knowledge on the complex mechanisms and regulatory processes, which govern the response of Hox genes to RA in several tissue contexts including the nervous system during vertebrate development.

A Hox-TALE regulatory circuit for neural crest patterning is conserved across vertebrates

In jawed vertebrates (gnathostomes), Hox genes play an important role in patterning head and jaw formation, but mechanisms coupling Hox genes to neural crest (NC) are unknown. Here we use cross-species regulatory comparisons between gnathostomes and lamprey, a jawless extant vertebrate, to investigate conserved ancestral mechanisms regulating Hox2 genes in NC. Gnathostome Hoxa2 and Hoxb2 NC enhancers mediate equivalent NC expression in lamprey and gnathostomes, revealing ancient conservation of Hox upstream regulatory components in NC. In characterizing a lamprey hoxα2 NC/hindbrain enhancer, we identify essential Meis, Pbx, and Hox binding sites that are functionally conserved within Hoxa2/Hoxb2 NC enhancers. This suggests that the lamprey hoxα2 enhancer retains ancestral activity and that Hoxa2/Hoxb2 NC enhancers are ancient paralogues, which diverged in hindbrain and NC activities. This identifies an ancestral mechanism for Hox2 NC regulation involving a Hox-TALE regulatory circuit, potentiated by inputs from Meis and Pbx proteins and Hox auto-/cross-regulatory interactions.

The phylum Vertebrata: a case for zoological recognition

The group Vertebrata is currently placed as a subphylum in the phylum Chordata, together with two other subphyla, Cephalochordata (lancelets) and Urochordata (ascidians). The past three decades, have seen extraordinary advances in zoological taxonomy and the time is now ripe for reassessing whether the subphylum position is truly appropriate for vertebrates, particularly in light of recent advances in molecular phylogeny, comparative genomics, and evolutionary developmental biology. Four lines of current research are discussed here. First, molecular phylogeny has demonstrated that Deuterostomia comprises Ambulacraria (Echinodermata and Hemichordata) and Chordata (Cephalochordata, Urochordata, and Vertebrata), each clade being recognized as a mutually comparable phylum. Second, comparative genomic studies show that vertebrates alone have experienced two rounds of whole-genome duplication, which makes the composition of their gene family unique. Third, comparative gene-expression profiling of vertebrate embryos favors an hourglass pattern of development, the most conserved stage of which is recognized as a phylotypic period characterized by the establishment of a body plan definitively associated with a phylum. This mid-embryonic conservation is supported robustly in vertebrates, but only weakly in chordates. Fourth, certain complex patterns of body plan formation (especially of the head, pharynx, and somites) are recognized throughout the vertebrates, but not in any other animal groups. For these reasons, we suggest that it is more appropriate to recognize vertebrates as an independent phylum, not as a subphylum of the phylum Chordata.

My Favorite Animal, Amphioxus: Unparalleled for Studying Early Vertebrate Evolution

Amphioxus represents the most basally divergent group in chordates and probably the best extant proxy to the ancestor of all chordates including vertebrates. The amphioxus, or lancelets, are benthic filter feeding marine animals and their interest as a model in research is due to their phylogenetic position and their anatomical and genetic stasis throughout their evolutionary history. From the first works in the 19th century to the present day, enormous progress is made mainly favored by technical development at different levels, from spawning induction and husbandry techniques, through techniques for studies of gene function or of the role of different signalling pathways through embryonic development, to functional genomics techniques. Together, these advances foretell a plethora of interesting developments in the world of research with the amphioxus model. Here, the discovery and development of amphioxus as a superb model organism in evolutionary and evolutionary-developmental biology are reviewed.

Xenacoelomorpha Survey Reveals That All 11 Animal Homeobox Gene Classes Were Present in the First Bilaterians

Homeodomain transcription factors are involved in many developmental processes across animals and have been linked to body plan evolution. Detailed classifications of these proteins identified 11 distinct classes of homeodomain proteins in animal genomes, each harboring specific sequence composition and protein domains. Although humans contain the full set of classes, Drosophila melanogaster and Caenorhabditis elegans each lack one specific class. Furthermore, representative previous analyses in sponges, ctenophores, and cnidarians could not identify several classes in those nonbilaterian metazoan taxa. Consequently, it is currently unknown when certain homeodomain protein classes first evolved during animal evolution. Here, we investigate representatives of the sister group to all remaining bilaterians, the Xenacoelomorpha. We analyzed three acoel, one nemertodermatid, and one Xenoturbella transcriptomes and identified their expressed homeodomain protein content. We report the identification of representatives of all 11 classes of animal homeodomain transcription factors in Xenacoelomorpha and we describe and classify their homeobox genes relative to the established animal homeodomain protein families. Our findings suggest that the genome of the last common ancestor of bilateria contained the full set of these gene classes, supporting the subsequent diversification of bilaterians.

Wnt evolution and function shuffling in liberal and conservative chordate genomes

BACKGROUND

What impact gene loss has on the evolution of developmental processes, and how function shuffling has affected retained genes driving essential biological processes, remain open questions in the fields of genome evolution and EvoDevo. To investigate these problems, we have analyzed the evolution of the Wnt ligand repertoire in the chordate phylum as a case study.

RESULTS

We conduct an exhaustive survey of Wnt genes in genomic databases, identifying 156 Wnt genes in 13 non-vertebrate chordates. This represents the most complete Wnt gene catalog of the chordate subphyla and has allowed us to resolve previous ambiguities about the orthology of many Wnt genes, including the identification of WntA for the first time in chordates. Moreover, we create the first complete expression atlas for the Wnt family during amphioxus development, providing a useful resource to investigate the evolution of Wnt expression throughout the radiation of chordates.

CONCLUSIONS

Our data underscore extraordinary genomic stasis in cephalochordates, which contrasts with the liberal and dynamic evolutionary patterns of gene loss and duplication in urochordate genomes. Our analysis has allowed us to infer ancestral Wnt functions shared among all chordates, several cases of function shuffling among Wnt paralogs, as well as unique expression domains for Wnt genes that likely reflect functional innovations in each chordate lineage. Finally, we propose a potential relationship between the evolution of WntA and the evolution of the mouth in chordates.

Evolution of dopamine receptors: phylogenetic evidence suggests a later origin of the DRD2l and DRD4rs dopamine receptor gene lineages

Dopamine receptors are integral membrane proteins whose endogenous ligand is dopamine. They play a fundamental role in the central nervous system and dysfunction of dopaminergic neurotransmission is responsible for the generation of a variety of neuropsychiatric disorders. From an evolutionary standpoint, phylogenetic relationships among the DRD₁ class of dopamine receptors are still a matter of debate as in the literature different tree topologies have been proposed. In contrast, phylogenetic relationships among the DRD₂ group of receptors are well understood. Understanding the time of origin of the different dopamine receptors is also an issue that needs further study, especially for the genes that have restricted phyletic distributions (e.g., DRD_2l and DRD_4rs). Thus, the goal of this study was to investigate the evolution of dopamine receptors, with emphasis on shedding light on the phylogenetic relationships among the D₁ class of dopamine receptors and the time of origin of the DRD_2l and DRD_4rs gene lineages. Our results recovered the monophyly of the two groups of dopamine receptors. Within the DRD₁ group the monophyly of each paralog was recovered with strong support, and phylogenetic relationships among them were well resolved. Within the DRD₁ class of dopamine receptors we recovered the sister group relationship between the DRD_1C and DRD_1E, and this clade was recovered sister to a cyclostome sequence. The DRD₁ clade was recovered sister to the aforementioned clade, and the group containing DRD₅ receptors was sister to all other DRD₁ paralogs. In agreement with the literature, among the DRD₂ class of receptors, DRD₂ was recovered sister to DRD₃, whereas DRD₄ was sister to the DRD₂/DRD₃ clade. According to our phylogenetic tree, the DRD_2l and DRD_4rs gene lineages would have originated in the ancestor of gnathostomes between 615 and 473 mya. Conservation of sequences required for dopaminergic neurotransmission and small changes in regulatory regions suggest a functional refinement of the dopaminergic pathways along evolution.

Coupling the roles of Hox genes to regulatory networks patterning cranial neural crest

The neural crest is a transient population of cells that forms within the developing central nervous system and migrates away to generate a wide range of derivatives throughout the body during vertebrate embryogenesis. These cells are of evolutionary and clinical interest, constituting a key defining trait in the evolution of vertebrates and alterations in their development are implicated in a high proportion of birth defects and craniofacial abnormalities. In the hindbrain and the adjacent cranial neural crest cells (cNCCs), nested domains of Hox gene expression provide a combinatorial'Hox-code' for specifying regional properties in the developing head. Hox genes have been shown to play important roles at multiple stages in cNCC development, including specification, migration, and differentiation. However, relatively little is known about the underlying gene-regulatory mechanisms involved, both upstream and downstream of Hox genes. Furthermore, it is still an open question as to how the genes of the neural crest GRN are linked to Hox-dependent pathways. In this review, we describe Hox gene expression, function and regulation in cNCCs with a view to integrating these genes within the emerging gene regulatory network for cNCC development. We highlight early roles for Hox1 genes in cNCC specification, proposing that this may be achieved, in part, by regulation of the balance between pluripotency and differentiation in precursor cells within the neuro-epithelium. We then describe what is known about the regulation of Hox gene expression in cNCCs and discuss this from the perspective of early vertebrate evolution.

Segmental arithmetic: summing up the Hox gene regulatory network for hindbrain development in chordates

Organization and development of the early vertebrate hindbrain are controlled by a cascade of regulatory interactions that govern the process of segmentation and patterning along the anterior-posterior axis via Hox genes. These interactions can be assembled into a gene regulatory network that provides a framework to interpret experimental data, generate hypotheses, and identify gaps in our understanding of the progressive process of hindbrain segmentation. The network can be broadly separated into a series of interconnected programs that govern early signaling, segmental subdivision, secondary signaling, segmentation, and ultimately specification of segmental identity. Hox genes play crucial roles in multiple programs within this network. Furthermore, the network reveals properties and principles that are likely to be general to other complex developmental systems. Data from vertebrate and invertebrate chordate models are shedding light on the origin and diversification of the network. Comprehensive cis-regulatory analyses of vertebrate Hox gene regulation have enabled powerful cross-species gene regulatory comparisons. Such an approach in the sea lamprey has revealed that the network mediating segmental Hox expression was present in ancestral vertebrates and has been maintained across diverse vertebrate lineages. Invertebrate chordates lack hindbrain segmentation but exhibit conservation of some aspects of the network, such as a role for retinoic acid in establishing nested Hox expression domains. These comparisons lead to a model in which early vertebrates underwent an elaboration of the network between anterior-posterior patterning and Hox gene expression, leading to the gene-regulatory programs for segmental subdivision and rhombomeric segmentation. WIREs Dev Biol 2017, 6:e286. doi: 10.1002/wdev.286 For further resources related to this article, please visit the WIREs website.

Cloning, expression and antioxidant activity of a thioredoxin peroxidase from Branchiostoma belcheri tsingtaunese

Peroxiredoxins (Prxs) are ubiquitous antioxidant enzymes that catalyze the thioredoxin- dependent reduction of hydroperoxides. In this study, a novel thioredoxin peroxidase (Bbt-TPx1), a member of the peroxiredoxin superfamily, was found by EST sequence analysis of a cDNA library of Branchiostoma belcheri tsingtaunese ovary. The sequence of a full-length cDNA clone contained an open reading frame encoding a polypeptide of 198 amino acid residues, with a calculated molecular weight of 22,150 Da. The expression patterns of the protein at different developmental stages and adult amphioxus tissues indicate that this enzyme may play important roles in anti-oxidation and innate immunity. The recombinant Bbt-TPx1 protein was expressed with a polyhistidine-tag in Escherichia coli and purified using Ni chromatography followed by SP cation exchange chromatography. The rBbt-TPx1 protein existed as a dimer under non-reducing conditions, and was dissociated into monomers by dithiothreitol (DTT); it might predominantly exist in oligomeric form. The rBbt-TPx1 protein showed a significant thiol-dependent peroxidase activity, removing hydrogen peroxide in the presence of dithiothreitol (DTT), but not glutathione (GSH). Protection of plasmid DNA and the thiol-protein from damage by metal-catalyzed oxidation (MCO) in vitro was also revealed.

Molecular regionalization of the developing amphioxus neural tube challenges major partitions of the vertebrate brain

All vertebrate brains develop following a common Bauplan defined by anteroposterior (AP) and dorsoventral (DV) subdivisions, characterized by largely conserved differential expression of gene markers. However, it is still unclear how this Bauplan originated during evolution. We studied the relative expression of 48 genes with key roles in vertebrate neural patterning in a representative amphioxus embryonic stage. Unlike nonchordates, amphioxus develops its central nervous system (CNS) from a neural plate that is homologous to that of vertebrates, allowing direct topological comparisons. The resulting genoarchitectonic model revealed that the amphioxus incipient neural tube is unexpectedly complex, consisting of several AP and DV molecular partitions. Strikingly, comparison with vertebrates indicates that the vertebrate thalamus, pretectum, and midbrain domains jointly correspond to a single amphioxus region, which we termed Di-Mesencephalic primordium (DiMes). This suggests that these domains have a common developmental and evolutionary origin, as supported by functional experiments manipulating secondary organizers in zebrafish and mice.

Molecular cytogenetic differentiation of paralogs of Hox paralogs in duplicated and re-diploidized genome of the North American paddlefish (Polyodon spathula)

Background

Acipenseriformes is a basal lineage of ray-finned fishes and comprise 27 extant species of sturgeons and paddlefishes. They are characterized by several specific genomic features as broad ploidy variation, high chromosome numbers, presence of numerous microchromosomes and propensity to interspecific hybridization. The presumed palaeotetraploidy of the American paddlefish was recently validated by molecular phylogeny and Hox genes analyses. A whole genome duplication in the paddlefish lineage was estimated at approximately 42 Mya and was found to be independent from several genome duplications evidenced in its sister lineage, i.e. sturgeons. We tested the ploidy status of available chromosomal markers after the expected rediploidization. Further we tested, whether paralogs of Hox gene clusters originated from this paddlefish specific genome duplication are cytogenetically distinguishable.

Results

We found that both paralogs HoxA alpha and beta were distinguishable without any overlapping of the hybridization signal - each on one pair of large metacentric chromosomes. Of the HoxD, only the beta paralog was unequivocally identified, whereas the alpha paralog did not work and yielded only an inconclusive diffuse signal. Chromosomal markers on three diverse ploidy levels reflecting different stages of rediploidization were identified: quadruplets retaining their ancestral tetraploid condition, semi-quadruplets still reflecting the ancestral tetraploidy with clear signs of advanced rediploidization, doublets were diploidized with ancestral tetraploidy already blurred. Also some of the available microsatellite data exhibited diploid allelic band patterns at their loci whereas another locus showed more than two alleles.

Conclusions

Our exhaustive staining of paddlefish chromosomes combined with cytogenetic mapping of ribosomal genes and Hox paralogs and with microsatellite data, brings a closer look at results of the process of rediploidization in the course of paddlefish genome evolution. We show a partial rediploidization represented by a complex mosaic structure comparable with segmental paleotetraploidy revealed in sturgeons (Acipenseridae). Sturgeons and paddlefishes with their high propensity for whole genome duplication thus offer suitable animal model systems to further explore evolutionary processes that were shaping the early evolution of all vertebrates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0484-8) contains supplementary material, which is available to authorized users.

Identification and conservation of gene loss events of Hox gene clusters in the marine medaka (Oryzias melastigma)

In this study, the identification of the whole Hox gene clusters (46 Hox genes) in the marine medaka Oryzias melastigma was investigated using genome assembly and RNA-seq information. Moreover, the gene loss events of Hox gene clusters, which may occur during fish evolution, were examined for a better understanding of the evolutionary status of the gene lost events of the Hox gene cluster across fish species, particularly in the genus Oryzias.

T-Box Genes and Developmental Gene Regulatory Networks in Ascidians

Ascidians are invertebrate chordates with a biphasic life cycle characterized by a dual body plan that displays simplified versions of chordate structures, such as a premetamorphic 40-cell notochord topped by a dorsal nerve cord and postmetamorphic pharyngeal slits. These relatively simple chordates are characterized by rapid development, compact genomes and ease of transgenesis, and thus provide the opportunity to rapidly characterize the genomic organization, developmental function, and transcriptional regulation of evolutionarily conserved gene families. This review summarizes the current knowledge on members of the T-box family of transcription factors in Ciona and other ascidians. In both chordate and nonchordate animals, these genes control a variety of morphogenetic processes, and their mutations are responsible for malformations and developmental defects in organisms ranging from flies to humans. In ascidians, T-box transcription factors are required for the formation and specialization of essential structures, including notochord, muscle, heart, and differentiated neurons. In recent years, the experimental advantages offered by ascidian embryos have allowed the rapid accumulation of a wealth of information on the molecular mechanisms that regulate the expression of T-box genes. These studies have also elucidated the strategies employed by these transcription factors to orchestrate the appropriate spatial and temporal deployment of the numerous target genes that they control.

Brain evolution by brain pathway duplication

Understanding the mechanisms of evolution of brain pathways for complex behaviours is still in its infancy. Making further advances requires a deeper understanding of brain homologies, novelties and analogies. It also requires an understanding of how adaptive genetic modifications lead to restructuring of the brain. Recent advances in genomic and molecular biology techniques applied to brain research have provided exciting insights into how complex behaviours are shaped by selection of novel brain pathways and functions of the nervous system. Here, we review and further develop some insights to a new hypothesis on one mechanism that may contribute to nervous system evolution, in particular by brain pathway duplication. Like gene duplication, we propose that whole brain pathways can duplicate and the duplicated pathway diverge to take on new functions. We suggest that one mechanism of brain pathway duplication could be through gene duplication, although other mechanisms are possible. We focus on brain pathways for vocal learning and spoken language in song-learning birds and humans as example systems. This view presents a new framework for future research in our understanding of brain evolution and novel behavioural traits.

Reorganisation of Hoxd regulatory landscapes during the evolution of a snake-like body plan

Within land vertebrate species, snakes display extreme variations in their body plan, characterized by the absence of limbs and an elongated morphology. Such a particular interpretation of the basic vertebrate body architecture has often been associated with changes in the function or regulation of Hox genes. Here, we use an interspecies comparative approach to investigate different regulatory aspects at the snake HoxD locus. We report that, unlike in other vertebrates, snake mesoderm-specific enhancers are mostly located within the HoxD cluster itself rather than outside. In addition, despite both the absence of limbs and an altered Hoxd gene regulation in external genitalia, the limb-associated bimodal HoxD chromatin structure is maintained at the snake locus. Finally, we show that snake and mouse orthologous enhancer sequences can display distinct expression specificities. These results show that vertebrate morphological evolution likely involved extensive reorganisation at Hox loci, yet within a generally conserved regulatory framework.

Location, location, location: the evolutionary history of CD1 genes and the NKR-P1/ligand systems

CD1 genes encode cell surface molecules that present lipid antigens to various kinds of T lymphocytes of the immune system. The structures of CD1 genes and molecules are like the major histocompatibility complex (MHC) class I system, the loading of antigen and the tissue distribution for CD1 molecules are like those in the class II system, and phylogenetic analyses place CD1 between class I and class II sequences, altogether leading to the notion that CD1 is a third ancient system of antigen presentation molecules. However, thus far, CD1 genes have only been described in mammals, birds and reptiles, leaving major questions as to their origin and evolution. In this review, we recount a little history of the field so far and then consider what has been learned about the structure and functional attributes of CD1 genes and molecules in marsupials, birds and reptiles. We describe the central conundrum of CD1 evolution, the genomic location of CD1 genes in the MHC and/or MHC paralogous regions in different animals, considering the three models of evolutionary history that have been proposed. We describe the natural killer (NK) receptors NKR-P1 and ligands, also found in different genomic locations for different animals. We discuss the consequence of these three models, one of which includes the repudiation of a guiding principle for the last 20 years, that two rounds of genome-wide duplication at the base of the vertebrates provided the extra MHC genes necessary for the emergence of adaptive immune system of jawed vertebrates.

Composition and genomic organization of arthropod Hox clusters

Background

The ancestral arthropod is believed to have had a clustered arrangement of ten Hox genes. Within arthropods, Hox gene mutations result in transformation of segment identities. Despite the fact that variation in segment number/character was common in the diversification of arthropods, few examples of Hox gene gains/losses have been correlated with morphological evolution. Furthermore, a full appreciation of the variation in the genomic arrangement of Hox genes in extant arthropods has not been recognized, as genome sequences from each major arthropod clade have not been reported until recently. Initial genomic analysis of the chelicerate Tetranychusurticae suggested that loss of Hox genes and Hox gene clustering might be more common than previously assumed. To further characterize the genomic evolution of arthropod Hox genes, we compared the genomic arrangement and general characteristics of Hox genes from representative taxa from each arthropod subphylum.

Results

In agreement with others, we find arthropods generally contain ten Hox genes arranged in a common orientation in the genome, with an increasing number of sampled species missing either Hox3 or abdominal-A orthologs. The genomic clustering of Hox genes in species we surveyed varies significantly, ranging from 0.3 to 13.6 Mb. In all species sampled, arthropod Hox genes are dispersed in the genome relative to the vertebrate Mus musculus. Differences in Hox cluster size arise from variation in the number of intervening genes, intergenic spacing, and the size of introns and UTRs. In the arthropods surveyed, Hox gene duplications are rare and four microRNAs are, in general, conserved in similar genomic positions relative to the Hox genes.

Conclusions

The tightly clustered Hox complexes found in the vertebrates are not evident within arthropods, and differential patterns of Hox gene dispersion are found throughout the arthropods. The comparative genomic data continue to support an ancestral arthropod Hox cluster of ten genes with a shared orientation, with four Hox gene-associated miRNAs, although the degree of dispersion between genes in an ancestral cluster remains uncertain. Hox3 and abdominal-A orthologs have been lost in multiple, independent lineages, and current data support a model in which inversions of the Abdominal-B locus that result in the loss of abdominal-A correlate with reduced trunk segmentation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-016-0048-4) contains supplementary material, which is available to authorized users.

Conservation of Hox gene clusters in the self-fertilizing fish Kryptolebias marmoratus (Cyprinodontiformes; Rivulidae)

In this study, whole Hox gene clusters in the self-fertilizing mangrove killifish Kryptolebias marmoratus (Cyprinodontiformes; Rivulidae), a unique hermaphroditic vertebrate in which both sex organs are functional at the same time, were identified from whole genome and transcriptome sequences. The aim was to increase the understanding of the evolutionary status of conservation of this Hox gene cluster across fish species.

Neofunctionalization of Androgen Receptor by Gain-of-Function Mutations in Teleost Fish Lineage

Steroid hormone receptor family provides an example of evolution of diverse transcription factors through whole-genome duplication (WGD). However, little is known about how their functions have been evolved after the duplication. Teleosts present a good model to investigate an accurate evolutionary history of protein function after WGD, because a teleost-specific WGD (TSGD) resulted in a variety of duplicated genes in modern fishes. This study focused on the evolution of androgen receptor (AR) gene, as two distinct paralogs, ARα and ARβ, have evolved in teleost lineage after TSGD. ARα showed a unique intracellular localization with a higher transactivation response than that of ARβ. Using site-directed mutagenesis and computational prediction of protein-ligand interactions, we identified two key substitutions generating a new functionality of euteleost ARα. The substitution in the hinge region contributes to the unique intracellular localization of ARα. The substitution on helices 10/11 in the ligand-binding domain possibly modulates hydrogen bonds that stabilize the receptor-ligand complex leading to the higher transactivation response of ARα. These substitutions were conserved in Acanthomorpha (spiny-rayed fish) ARαs, but not in an earlier branching lineage among teleosts, Japanese eel. Insertion of these substitutions into ARs from Japanese eel recapitulates the evolutionary novelty of euteleost ARα. These findings together indicate that the substitutions generating a new functionality of teleost ARα were fixed in teleost genome after the divergence of the Elopomorpha lineage. Our findings provide a molecular explanation for an adaptation process leading to generation of the hyperactive AR subtype after TSGD.

Cdx ParaHox genes acquired distinct developmental roles after gene duplication in vertebrate evolution

Background

The functional consequences of whole genome duplications in vertebrate evolution are not fully understood. It remains unclear, for instance, why paralogues were retained in some gene families but extensively lost in others. Cdx homeobox genes encode conserved transcription factors controlling posterior development across diverse bilaterians. These genes are part of the ParaHox gene cluster. Multiple Cdx copies were retained after genome duplication, raising questions about how functional divergence, overlap, and redundancy respectively contributed to their retention and evolutionary fate.

Results

We examined the degree of regulatory and functional overlap between the three vertebrate Cdx genes using single and triple morpholino knock-down in Xenopus tropicalis followed by RNA-seq. We found that one paralogue, Cdx4, has a much stronger effect on gene expression than the others, including a strong regulatory effect on FGF and Wnt genes. Functional annotation revealed distinct and overlapping roles and subtly different temporal windows of action for each gene. The data also reveal a colinear-like effect of Cdx genes on Hox genes, with repression of Hox paralogy groups 1 and 2, and activation increasing from Hox group 5 to 11. We also highlight cases in which duplicated genes regulate distinct paralogous targets revealing pathway elaboration after whole genome duplication.

Conclusions

Despite shared core pathways, Cdx paralogues have acquired distinct regulatory roles during development. This implies that the degree of functional overlap between paralogues is relatively low and that gene expression pattern alone should be used with caution when investigating the functional evolution of duplicated genes. We therefore suggest that developmental programmes were extensively rewired after whole genome duplication in the early evolution of vertebrates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0165-x) contains supplementary material, which is available to authorized users.

A Simple Predictive Enhancer Syntax for Hindbrain Patterning Is Conserved in Vertebrate Genomes

Background

Determining the function of regulatory elements is fundamental for our understanding of development, disease and evolution. However, the sequence features that mediate these functions are often unclear and the prediction of tissue-specific expression patterns from sequence alone is non-trivial. Previous functional studies have demonstrated a link between PBX-HOX and MEIS/PREP binding interactions and hindbrain enhancer activity, but the defining grammar of these sites, if any exists, has remained elusive.

Results

Here, we identify a shared sequence signature (syntax) within a heterogeneous set of conserved vertebrate hindbrain enhancers composed of spatially co-occurring PBX-HOX and MEIS/PREP transcription factor binding motifs. We use this syntax to accurately predict hindbrain enhancers in 89% of cases (67/75 predicted elements) from a set of conserved non-coding elements (CNEs). Furthermore, mutagenesis of the sites abolishes activity or generates ectopic expression, demonstrating their requirement for segmentally restricted enhancer activity in the hindbrain. We refine and use our syntax to predict over 3,000 hindbrain enhancers across the human genome. These sequences tend to be located near developmental transcription factors and are enriched in known hindbrain activating elements, demonstrating the predictive power of this simple model.

Conclusion

Our findings support the theory that hundreds of CNEs, and perhaps thousands of regions across the human genome, function to coordinate gene expression in the developing hindbrain. We speculate that deeply conserved sequences of this kind contributed to the co-option of new genes into the hindbrain gene regulatory network during early vertebrate evolution by linking patterns of hox expression to downstream genes involved in segmentation and patterning, and evolutionarily newer instances may have continued to contribute to lineage-specific elaboration of the hindbrain.

Hox genes pattern the anterior-posterior axis of the juvenile but not the larva in a maximally indirect developing invertebrate, Micrura alaskensis (Nemertea)

BACKGROUND

The pilidium larva is a novel body plan that arose within a single clade in the phylum Nemertea - the Pilidiophora. While the sister clade of the Pilidiophora and the basal nemerteans develop directly, pilidiophorans have a long-lived planktotrophic larva with a body plan distinctly different from that of the juvenile. Uniquely, the pilidiophoran juvenile develops inside the larva from several discrete rudiments. The orientation of the juvenile with respect to the larval body varies within the Pilidiophora, which suggests that the larval and juvenile anteroposterior (AP) axes are patterned differently. In order to gain insight into the evolutionary origins of the pilidium larva and the mechanisms underlying this implied axial uncoupling, we examined the expression of the Hox genes during development of the pilidiophoran Micrura alaskensis.

RESULTS

We identified sequences of nine Hox genes and the ParaHox gene caudal through a combination of transcriptome analysis and molecular cloning, and determined their expression pattern during development using in situ hybridization in whole-mounted larvae. We found that Hox genes are first expressed long after the pilidium is fully formed and functional. The Hox genes are expressed in apparently overlapping domains along the AP axis of the developing juvenile in a subset of the rudiments that give rise to the juvenile trunk. Hox genes are not expressed in the larval body at any stage of development.

CONCLUSIONS

While the Hox genes pattern the juvenile pilidiophoran, the pilidial body, which appears to be an evolutionary novelty, must be patterned by some mechanism other than the Hox genes. Although the pilidiophoran juvenile develops from separate rudiments with no obvious relationship to the embryonic formation of the larva, the Hox genes appear to exhibit canonical expression along the juvenile AP axis. This suggests that the Hox patterning system can maintain conserved function even when widely decoupled from early polarity established in the egg.

Enigmatic orthology relationships between Hox clusters of the African butterfly fish and other teleosts following ancient whole-genome duplication

Numerous ancient whole-genome duplications (WGD) have occurred during eukaryote evolution. In vertebrates, duplicated developmental genes and their functional divergence have had important consequences for morphological evolution. Although two vertebrate WGD events (1R/2R) occurred over 525 Ma, we have focused on the more recent 3R or TGD (teleost genome duplication) event which occurred approximately 350 Ma in a common ancestor of over 26,000 species of teleost fishes. Through a combination of whole genome and bacterial artificial chromosome clone sequencing we characterized all Hox gene clusters of Pantodon buchholzi, a member of the early branching teleost subdivision Osteoglossomorpha. We find 45 Hox genes organized in only five clusters indicating that Pantodon has suffered more Hox cluster loss than other known species. Despite strong evidence for homology of the five Pantodon clusters to the four canonical pre-TGD vertebrate clusters (one HoxA, two HoxB, one HoxC, and one HoxD), we were unable to confidently resolve 1:1 orthology relationships between four of the Pantodon clusters and the eight post-TGD clusters of other teleosts. Phylogenetic analysis revealed that many Pantodon genes segregate outside the conventional "a" and "b" post-TGD orthology groups, that extensive topological incongruence exists between genes physically linked on a single cluster, and that signal divergence causes ambivalence in assigning 1:1 orthology in concatenated Hox cluster analyses. Out of several possible explanations for this phenomenon we favor a model which keeps with the prevailing view of a single TGD prior to teleost radiation, but which also considers the timing of diploidization after duplication, relative to speciation events. We suggest that although the duplicated hoxa clusters diploidized prior to divergence of osteoglossomorphs, the duplicated hoxb, hoxc, and hoxd clusters concluded diploidization independently in osteoglossomorphs and other teleosts. We use the term "tetralogy" to describe the homology relationship which exists between duplicated sequences which originate through a shared WGD, but which diploidize into distinct paralogs from a common allelic pool independently in two lineages following speciation.