A reporter assay in lamprey embryos reveals both functional conservation and elaboration of vertebrate enhancers

Sea lamprey enlightens the origin of the coupling of retinoic acid signaling to vertebrate hindbrain segmentation

Retinoic acid (RA) is involved in antero-posterior patterning of the chordate body axis and, in jawed vertebrates, has been shown to play a major role at multiple levels of the gene regulatory network (GRN) regulating hindbrain segmentation. Knowing when and how RA became coupled to the core hindbrain GRN is important for understanding how ancient signaling pathways and patterning genes can evolve and generate diversity. Hence, we investigated the link between RA signaling and hindbrain segmentation in the sea lamprey Petromyzon marinus, an important jawless vertebrate model providing clues to decipher ancestral vertebrate features. Combining genomics, gene expression, and functional analyses of major components involved in RA synthesis (Aldh1as) and degradation (Cyp26s), we demonstrate that RA signaling is coupled to hindbrain segmentation in lamprey. Thus, the link between RA signaling and hindbrain segmentation is a pan vertebrate feature of the hindbrain and likely evolved at the base of vertebrates.

Molecular evolution of the heat shock protein family and the role of HSP30 in immune response and wound healing in lampreys (Lethenteron reissneri)

Heat shock proteins (HSPs) are molecular chaperones that ubiquitously exist in various organisms and play essential roles in protein folding, transport, and expression. While most HSPs are highly conserved across species, a few HSPs are evolutionarily distinct in some species and may have unique functions. To explore the evolutionary history of the vertebrate HSP family, we identify members of the HSP family at the genome-wide level in lampreys (Lethenteron reissneri), a living representative of jawless vertebrates diverged from jawed vertebrates over 500 million years ago. The phylogenetic analysis reveals that the lamprey HSP family contains HSP90a1, HSP90a2, HSC70, HSP60, HSP30, HSP27, HSP17, and HSP10, which have a primitive status in the molecular evolution of vertebrate HSPs. Transcriptome analysis reveals the expression distribution of members of the HSP family in various tissues of lampreys. It is shown that HSP30, normally found in birds, amphibians, and fish, is also present in lampreys, with remarkable expansion of HSP30 gene copies in the lamprey genome. The transcription of HSP30 is significantly induced in leukocytes and heart of lampreys during various pathogens or poly(I:C) stimulation, indicating that HSP30 may be involved in the immune defense of lampreys in response to bacterial or viral infection. Immunohistochemistry demonstrates significantly increased HSP30 expression in subcutaneous muscle tissue after skin injury in lamprey models of wound repair. Furthermore, transcriptome analysis shows that ectopic expression of HSP30 in 3T3-L1 fibroblasts affect the expression of genes related to the PI3K-AKT signaling pathway, suggesting that HSP30 could serves as a negative regulator of fibrosis. These results indicate that HSP30 may play a critical role in facilitating the process of lamprey skin repair following injury. This study provides new insights into the origin and evolution of the HSP gene family in vertebrates and offers valuable clues to reveal the important role of HSP30 in immune defense and wound healing of lampreys.

Evolution of the expression and regulation of the nuclear hormone receptor ERR gene family in the chordate lineage

The Estrogen Related Receptor (ERR) nuclear hormone receptor genes have a wide diversity of roles in vertebrate development. In embryos, ERR genes are expressed in several tissues, including the central and peripheral nervous systems. Here we seek to establish the evolutionary history of chordate ERR genes, their expression and their regulation. We examine ERR expression in mollusc, amphioxus and sea squirt embryos, finding the single ERR orthologue is expressed in the nervous system in all three, with muscle expression also found in the two chordates. We show that most jawed vertebrates and lampreys have four ERR paralogues, and that vertebrate ERR genes were ancestrally linked to Estrogen Receptor genes. One of the lamprey paralogues shares conserved expression domains with jawed vertebrate ERRγ in the embryonic vestibuloacoustic ganglion, eye, brain and spinal cord. Hypothesising that conserved expression derives from conserved regulation, we identify a suite of pan-vertebrate conserved non-coding sequences in ERR introns. We use transgenesis in lamprey and chicken embryos to show that these sequences are regulatory and drive reporter gene expression in the nervous system. Our data suggest an ancient association between ERR and the nervous system, including expression in cells associated with photosensation and mechanosensation. This includes the origin in the vertebrate common ancestor of a suite of regulatory elements in the 3' introns that drove nervous system expression and have been conserved from this point onwards.

Sea lamprey enlightens the origin of the coupling of retinoic acid signaling to vertebrate hindbrain segmentation

Retinoic acid (RA) is involved in antero-posterior patterning of the chordate body axis and, in jawed vertebrates, has been shown to play a major role at multiple levels of the gene regulatory network (GRN) regulating hindbrain segmentation. Knowing when and how RA became coupled to the core hindbrain GRN is important for understanding how ancient signaling pathways and patterning genes can evolve and generate diversity. Hence, we investigated the link between RA signaling and hindbrain segmentation in the sea lamprey Petromyzon marinus, an important jawless vertebrate model providing clues to decipher ancestral vertebrate features. Combining genomics, gene expression, and functional analyses of major components involved in RA synthesis (Aldh1as) and degradation (Cyp26s), we demonstrate that RA signaling is coupled to hindbrain segmentation in lamprey. Thus, the link between RA signaling and hindbrain segmentation is a pan vertebrate feature of the hindbrain and likely evolved at the base of vertebrates.

Fishing for Developmental Regulatory Regions: Zebrafish Tissue-Specific ATAC-seq

Interactions between transcription factors and regulatory DNA can be described by gene regulatory networks. These networks provide a systems-level view of embryonic tissue development. Here, we describe a protocol for the isolation, identification, and experimental manipulation of tissue-specific cis-regulatory elements during zebrafish embryonic development using low-input ATAC-seq. With the methods described, genome-wide assessments of regulatory DNA in small populations of developing tissues can be identified, allowing for the construction of gene regulatory networks.

Hmx gene conservation identifies the origin of vertebrate cranial ganglia

The evolutionary origin of vertebrates included innovations in sensory processing associated with the acquisition of a predatory lifestyle¹. Vertebrates perceive external stimuli through sensory systems serviced by cranial sensory ganglia, whose neurons arise predominantly from cranial placodes; however, the understanding of the evolutionary origin of placodes and cranial sensory ganglia is hampered by the anatomical differences between living lineages and the difficulty in assigning homology between cell types and structures. Here we show that the homeobox transcription factor Hmx is a constitutive component of vertebrate sensory ganglion development and that in the tunicate Ciona intestinalis, Hmx is necessary and sufficient to drive the differentiation programme of bipolar tail neurons, cells previously thought to be homologues of neural crest^2,3. Using Ciona and lamprey transgenesis, we demonstrate that a unique, tandemly duplicated enhancer pair regulated Hmx expression in the stem-vertebrate lineage. We also show notably robust vertebrate Hmx enhancer function in Ciona, demonstrating that deep conservation of the upstream regulatory network spans the evolutionary origin of vertebrates. These experiments demonstrate regulatory and functional conservation between Ciona and vertebrate Hmx, and point to bipolar tail neurons as homologues of cranial sensory ganglia.

Analysis of lamprey meis genes reveals that conserved inputs from Hox, Meis and Pbx proteins control their expression in the hindbrain and neural tube

Meis genes are known to play important roles in the hindbrain and neural crest cells of jawed vertebrates. To explore the roles of Meis genes in head development during evolution of vertebrates, we have identified four meis genes in the sea lamprey genome and characterized their patterns of expression and regulation, with a focus on the hindbrain and pharynx. Each of the lamprey meis genes displays temporally and spatially dynamic patterns of expression, some of which are coupled to rhombomeric domains in the developing hindbrain and select pharyngeal arches. Studies of Meis loci in mouse and zebrafish have identified enhancers that are bound by Hox and TALE (Meis and Pbx) proteins, implicating these factors in the direct regulation of Meis expression. We examined the lamprey meis loci and identified a series of cis-elements conserved between lamprey and jawed vertebrate meis genes. In transgenic reporter assays we demonstrated that these elements act as neural enhancers in lamprey embryos, directing reporter expression in appropriate domains when compared to expression of their associated endogenous meis gene. Sequence alignments reveal that these conserved elements are in similar relative positions of the meis loci and contain a series of consensus binding motifs for Hox and TALE proteins. This suggests that ancient Hox and TALE-responsive enhancers regulated expression of ancestral vertebrate meis genes in segmental domains in the hindbrain and have been retained in the meis loci during vertebrate evolution. The presence of conserved Meis, Pbx and Hox binding sites in these lamprey enhancers links Hox and TALE factors to regulation of lamprey meis genes in the developing hindbrain, indicating a deep ancestry for these regulatory interactions prior to the divergence of jawed and jawless vertebrates.

Generation of knock-in lampreys by CRISPR-Cas9-mediated genome engineering

The lamprey represents the oldest group of living vertebrates and has been a key organism in various research fields such as evolutionary developmental biology and neuroscience. However, no knock-in technique for this animal has been established yet, preventing application of advanced genetic techniques. Here, we report efficient generation of F₀ knock-in lampreys by CRISPR-Cas9-mediated genome editing. A donor plasmid containing a heat-shock promoter was co-injected with a short guide RNA (sgRNA) for genome digestion, a sgRNA for donor plasmid digestion, and Cas9 mRNA. Targeting different genetic loci, we succeeded in generating knock-in lampreys expressing photoconvertible protein Dendra2 as well as those expressing EGFP. With its simplicity, design flexibility, and high efficiency, we propose that the present method has great versatility for various experimental uses in lamprey research and that it can also be applied to other “non-model” organisms.

Programmed DNA Elimination in Vertebrates

Over the last few decades, an increasing number of vertebrate taxa have been identified that undergo programmed genome rearrangement, or programmed DNA loss, during development. In these organisms, the genome of germ cells is often reproducibly different from the genome of all other cells within the body. Although we clearly have not identified all vertebrate taxa that undergo programmed genome loss, the list of species known to undergo loss now represents ∼10% of vertebrate species, including several basally diverging lineages. Recent studies have shed new light on the targets and mechanisms of DNA loss and their association with canonical modes of DNA silencing. Ultimately, expansion of these studies into a larger collection of taxa will aid in reconstructing patterns of shared/independent ancestry of programmed DNA loss in the vertebrate lineage, as well as more recent evolutionary events that have shaped the structure and content of eliminated DNA.

Evolutionary and Developmental Associations of Neural Crest and Placodes in the Vertebrate Head: Insights From Jawless Vertebrates

Neural crest and placodes are key innovations of the vertebrate clade. These cells arise within the dorsal ectoderm of all vertebrate embryos and have the developmental potential to form many of the morphological novelties within the vertebrate head. Each cell population has its own distinct developmental features and generates unique cell types. However, it is essential that neural crest and placodes associate together throughout embryonic development to coordinate the emergence of several features in the head, including almost all of the cranial peripheral sensory nervous system and organs of special sense. Despite the significance of this developmental feat, its evolutionary origins have remained unclear, owing largely to the fact that there has been little comparative (evolutionary) work done on this topic between the jawed vertebrates and cyclostomes—the jawless lampreys and hagfishes. In this review, we briefly summarize the developmental mechanisms and genetics of neural crest and placodes in both jawed and jawless vertebrates. We then discuss recent studies on the role of neural crest and placodes—and their developmental association—in the head of lamprey embryos, and how comparisons with jawed vertebrates can provide insights into the causes and consequences of this event in early vertebrate evolution.

An evolutionarily ancient mechanism for regulation of hemoglobin expression in vertebrate red cells

The oxygen transport function of hemoglobin (HB) is thought to have arisen ∼500 million years ago, roughly coinciding with the divergence between jawless (Agnatha) and jawed (Gnathostomata) vertebrates. Intriguingly, extant HBs of jawless and jawed vertebrates were shown to have evolved twice, and independently, from different ancestral globin proteins. This raises the question of whether erythroid-specific expression of HB also evolved twice independently. In all jawed vertebrates studied to date, one of the HB gene clusters is linked to the widely expressed NPRL3 gene. Here we show that the nprl3-linked hb locus of a jawless vertebrate, the river lamprey (Lampetra fluviatilis), shares a range of structural and functional properties with the equivalent jawed vertebrate HB locus. Functional analysis demonstrates that an erythroid-specific enhancer is located in intron 7 of lamprey nprl3, which corresponds to the NPRL3 intron 7 MCS-R1 enhancer of jawed vertebrates. Collectively, our findings signify the presence of an nprl3-linked multiglobin gene locus, which contains a remote enhancer that drives globin expression in erythroid cells, before the divergence of jawless and jawed vertebrates. Different globin genes from this ancestral cluster evolved in the current NPRL3-linked HB genes in jawless and jawed vertebrates. This provides an explanation of the enigma of how, in different species, globin genes linked to the same adjacent gene could undergo convergent evolution.

Functional genetic analysis in a jawless vertebrate, the sea lamprey: insights into the developmental evolution of early vertebrates

Lampreys and hagfishes are the only surviving relicts of an ancient but ecologically dominant group of jawless fishes that evolved in the seas of the Cambrian era over half a billion years ago. Because of their phylogenetic position as the sister group to all other vertebrates (jawed vertebrates), comparisons of embryonic development between jawless and jawed vertebrates offers researchers in the field of evolutionary developmental biology the unique opportunity to address fundamental questions related to the nature of our earliest vertebrate ancestors. Here, we describe how genetic analysis of embryogenesis in the sea lamprey (Petromyzon marinus) has provided insight into the origin and evolution of developmental-genetic programs in vertebrates. We focus on recent work involving CRISPR/Cas9-mediated genome editing to study gene regulatory mechanisms involved in the development and evolution of neural crest cells and new cell types in the vertebrate nervous system, and transient transgenic assays that have been instrumental in dissecting the evolution of cis-regulatory control of gene expression in vertebrates. Finally, we discuss the broad potential for these functional genomic tools to address previously unanswerable questions related to the evolution of genomic regulatory mechanisms as well as issues related to invasive sea lamprey population control.

A genome-wide assessment of the ancestral neural crest gene regulatory network

The neural crest (NC) is an embryonic cell population that contributes to key vertebrate-specific features including the craniofacial skeleton and peripheral nervous system. Here we examine the transcriptional and epigenomic profiles of NC cells in the sea lamprey, in order to gain insight into the ancestral state of the NC gene regulatory network (GRN). Transcriptome analyses identify clusters of co-regulated genes during NC specification and migration that show high conservation across vertebrates but also identify transcription factors (TFs) and cell-adhesion molecules not previously implicated in NC migration. ATAC-seq analysis uncovers an ensemble of cis-regulatory elements, including enhancers of Tfap2B, SoxE1 and Hox-α2 validated in the embryo. Cross-species deployment of lamprey elements identifies the deep conservation of lamprey SoxE1 enhancer activity, mediating homologous expression in jawed vertebrates. Our data provide insight into the core GRN elements conserved to the base of the vertebrates and expose others that are unique to lampreys.

An atlas of anterior hox gene expression in the embryonic sea lamprey head: Hox-code evolution in vertebrates

In the hindbrain and the adjacent cranial neural crest (NC) cells of jawed vertebrates (gnathostomes), nested and segmentally-restricted domains of Hox gene expression provide a combinatorial Hox-code for specifying regional properties during head development. Extant jawless vertebrates, such as the sea lamprey (Petromyzon marinus), can provide insights into the evolution and diversification of this Hox-code in vertebrates. There is evidence for gnathostome-like spatial patterns of Hox expression in lamprey; however, the expression domains of the majority of lamprey hox genes from paralogy groups (PG) 1-4 are yet to be characterized, so it is unknown whether they are coupled to hindbrain segments (rhombomeres) and NC. In this study, we systematically describe the spatiotemporal expression of all 14 sea lamprey hox genes from PG1-PG4 in the developing hindbrain and pharynx to investigate the extent to which their expression conforms to the archetypal gnathostome hindbrain and pharyngeal hox-codes. We find many similarities in Hox expression between lamprey and gnathostome species, particularly in rhombomeric domains during hindbrain segmentation and in the cranial neural crest, enabling inference of aspects of Hox expression in the ancestral vertebrate embryonic head. These data are consistent with the idea that a Hox regulatory network underlying hindbrain segmentation is a pan vertebrate trait. We also reveal differences in hindbrain domains at later stages, as well as expression in the endostyle and in pharyngeal arch (PA) 1 mesoderm. Our analysis suggests that many Hox expression domains that are observed in extant gnathostomes were present in ancestral vertebrates but have been partitioned differently across Hox clusters in gnathostome and cyclostome lineages after duplication.

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.

I-SceI Meganuclease-mediated transgenesis in the acorn worm, Saccoglossus kowalevskii

Hemichordates are a phylum of marine invertebrate deuterostomes that are closely related to chordates, and represent one of the most promising models to provide insights into early deuterostome evolution. The genome of the hemichordate, Saccoglossus kowalevskii, reveals an extensive set of non-coding elements conserved across all three deuterostome phyla. Functional characterization and cross-phyla comparisons of these putative regulatory elements will enable a better understanding of enhancer evolution, and subsequently how changes in gene regulation give rise to morphological innovation. Here, we describe an efficient method of transgenesis for the characterization of non-coding elements in S. kowalevskii. We first test the capacity of an I-SceI transgenesis system to drive ubiquitous or regionalized gene expression, and to label specific cell types. Finally, we identified a minimal promoter that can be used to test the capacity of putative enhancers in S. kowalevskii. This work demonstrates that this I-SceI transgenesis technique, when coupled with an understanding of chromatin accessibility, can be a powerful tool for studying how evolutionary changes in gene regulatory mechanisms contributed to the diversification of body plans in deuterostomes.

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.

Conserved Noncoding Elements in the Most Distant Genera of Cephalochordates: The Goldilocks Principle

Cephalochordates, the sister group of vertebrates + tunicates, are evolving particularly slowly. Therefore, genome comparisons between two congeners of Branchiostoma revealed so many conserved noncoding elements (CNEs), that it was not clear how many are functional regulatory elements. To more effectively identify CNEs with potential regulatory functions, we compared noncoding sequences of genomes of the most phylogenetically distant cephalochordate genera, Asymmetron and Branchiostoma, which diverged approximately 120-160 million years ago. We found 113,070 noncoding elements conserved between the two species, amounting to 3.3% of the genome. The genomic distribution, target gene ontology, and enriched motifs of these CNEs all suggest that many of them are probably cis-regulatory elements. More than 90% of previously verified amphioxus regulatory elements were re-captured in this study. A search of the cephalochordate CNEs around 50 developmental genes in several vertebrate genomes revealed eight CNEs conserved between cephalochordates and vertebrates, indicating sequence conservation over >500 million years of divergence. The function of five CNEs was tested in reporter assays in zebrafish, and one was also tested in amphioxus. All five CNEs proved to be tissue-specific enhancers. Taken together, these findings indicate that even though Branchiostoma and Asymmetron are distantly related, as they are evolving slowly, comparisons between them are likely optimal for identifying most of their tissue-specific cis-regulatory elements laying the foundation for functional characterizations and a better understanding of the evolution of developmental regulation in cephalochordates.

The vertebrate Hox gene regulatory network for hindbrain segmentation: Evolution and diversification: Coupling of a Hox gene regulatory network to hindbrain segmentation is an ancient trait originating at the base of vertebrates

Hindbrain development is orchestrated by a vertebrate gene regulatory network that generates segmental patterning along the anterior-posterior axis via Hox genes. Here, we review analyses of vertebrate and invertebrate chordate models that inform upon the evolutionary origin and diversification of this network. Evidence from the sea lamprey reveals that the hindbrain regulatory network generates rhombomeric compartments with segmental Hox expression and an underlying Hox code. We infer that this basal feature was present in ancestral vertebrates and, as an evolutionarily constrained developmental state, is fundamentally important for patterning of the vertebrate hindbrain across diverse lineages. Despite the common ground plan, vertebrates exhibit neuroanatomical diversity in lineage-specific patterns, with different vertebrates revealing variations of Hox expression in the hindbrain that could underlie this diversification. Invertebrate chordates lack hindbrain segmentation but exhibit some conserved aspects of this network, with retinoic acid signaling playing a role in establishing nested domains of Hox expression.

Lampreys as Diverse Model Organisms in the Genomics Era

Lampreys, one of the two surviving groups of ancient vertebrates, have become important models for study in diverse fields of biology. Lampreys (of which there are approximately 40 species) are being studied, for example, (a) to control pest sea lamprey in the North American Great Lakes and to restore declining populations of native species elsewhere; (b) in biomedical research, focusing particularly on the regenerative capability of lampreys; and (c) by developmental biologists studying the evolution of key vertebrate characters. Although a lack of genetic resources has hindered research on the mechanisms regulating many aspects of lamprey life history and development, formerly intractable questions are now amenable to investigation following the recent publication of the sea lamprey genome. Here, we provide an overview of the ways in which genomic tools are currently being deployed to tackle diverse research questions and suggest several areas that may benefit from the availability of the sea lamprey genome.

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.

A Hox regulatory network of hindbrain segmentation is conserved to the base of vertebrates

A defining feature governing head patterning of jawed vertebrates is a highly conserved gene regulatory network that integrates hindbrain segmentation with segmentally restricted domains of Hox gene expression. Although non-vertebrate chordates display nested domains of axial Hox expression, they lack hindbrain segmentation. The sea lamprey, a jawless fish, can provide unique insights into vertebrate origins owing to its phylogenetic position at the base of the vertebrate tree. It has been suggested that lamprey may represent an intermediate state where nested Hox expression has not been coupled to the process of hindbrain segmentation. However, little is known about the regulatory network underlying Hox expression in lamprey or its relationship to hindbrain segmentation. Here, using a novel tool that allows cross-species comparisons of regulatory elements between jawed and jawless vertebrates, we report deep conservation of both upstream regulators and segmental activity of enhancer elements across these distant species. Regulatory regions from diverse gnathostomes drive segmental reporter expression in the lamprey hindbrain and require the same transcriptional inputs (for example, Kreisler (also known as Mafba), Krox20 (also known as Egr2a)) in both lamprey and zebrafish. We find that lamprey hox genes display dynamic segmentally restricted domains of expression; we also isolated a conserved exonic hox2 enhancer from lamprey that drives segmental expression in rhombomeres 2 and 4. Our results show that coupling of Hox gene expression to segmentation of the hindbrain is an ancient trait with origin at the base of vertebrates that probably led to the formation of rhombomeric compartments with an underlying Hox code.