Amphioxus functional genomics and the origins of vertebrate gene regulation

Assaying Chromatin Accessibility Using ATAC-Seq in Invertebrate Chordate Embryos

Cis-regulatory elements (CREs) are non-coding DNA regions involved in the spatio-temporal regulation of gene expression. Gene regulatory changes drive animal development and play major roles during evolution of animal body plans. Therefore, we believe that determining CREs at different developmental stages and across animal lineages is critical to understand how evolution operates through development. The Assay for Transposase-Accessible Chromatin followed by high-throughput sequencing (ATAC-seq) is a powerful technique for the study of CREs that takes advantage of Tn5 transposase activity. Starting from fewer than 10⁵ cells, in a 1-day procedure, it is possible to detect, at a genome-wide level, CREs located in open chromatin regions with high resolution. Here, we describe a detailed step-by-step ATAC-seq protocol for invertebrate chordate marine embryos. We have successfully applied this technique to amphioxus and two species of tunicate embryos. We also show an easy workflow to analyze data generated with this technique. Moreover, we point out that this method and our bioinformatic pipeline are efficient to detect CREs associated with Wnt signaling pathway by simply using embryos treated with a drug that perturbs this pathway. This approach can be extended to other signaling pathways and also to embryo mutants for critical genes. Our results therefore demonstrate the power of ATAC-seq for the identification of CREs that play essential functions during animal development in a wide range of invertebrate or vertebrate animals.

The origin and evolution of vertebrate neural crest cells

The neural crest is a vertebrate-specific migratory stem cell population that generates a remarkably diverse set of cell types and structures. Because many of the morphological, physiological and behavioural novelties of vertebrates are derived from neural crest cells, it is thought that the origin of this cell population was an important milestone in early vertebrate history. An outstanding question in the field of vertebrate evolutionary-developmental biology (evo-devo) is how this cell type evolved in ancestral vertebrates. In this review, we briefly summarize neural crest developmental genetics in vertebrates, focusing in particular on the gene regulatory interactions instructing their early formation within and migration from the dorsal neural tube. We then discuss how studies searching for homologues of neural crest cells in invertebrate chordates led to the discovery of neural crest-like cells in tunicates and the potential implications this has for tracing the pre-vertebrate origins of the neural crest population. Finally, we synthesize this information to propose a model to explain the origin of neural crest cells. We suggest that at least some of the regulatory components of early stages of neural crest development long pre-date vertebrate origins, perhaps dating back to the last common bilaterian ancestor. These components, originally directing neuroectodermal patterning and cell migration, served as a gene regulatory 'scaffold' upon which neural crest-like cells with limited migration and potency evolved in the last common ancestor of tunicates and vertebrates. Finally, the acquisition of regulatory programmes controlling multipotency and long-range, directed migration led to the transition from neural crest-like cells in invertebrate chordates to multipotent migratory neural crest in the first vertebrates.

Lampreys, the jawless vertebrates, contain three Pax6 genes with distinct expression in eye, brain and pancreas

The transcription factor Pax6 is crucial for the development of the central nervous system, eye, olfactory system and pancreas, and is implicated in human disease. While a single Pax6 gene exists in human and chicken, Pax6 occurs as a gene family in other vertebrates, with two members in elephant shark, Xenopus tropicalis and Anolis lizard and three members in teleost fish such as stickleback and medaka. However, the complement of Pax6 genes in jawless vertebrates (cyclostomes), the sister group of jawed vertebrates (gnathostomes), is unknown. Using a combination of BAC sequencing and genome analysis, we discovered three Pax6 genes in lampreys. Unlike the paired-less Pax6 present in some gnathostomes, all three lamprey Pax6 have a highly conserved full-length paired domain. All three Pax6 genes are expressed in the eye and brain, with variable expression in other tissues. Notably, lamprey Pax6α transcripts are found in the pancreas, a vertebrate-specific organ, indicating the involvement of Pax6 in development of the pancreas in the vertebrate ancestor. Multi-species sequence comparisons revealed only a single conserved non-coding element, in the lamprey Pax6β locus, with similarity to the PAX6 neuroretina enhancer. Using a transgenic zebrafish enhancer assay we demonstrate functional conservation of this element over 500 million years of vertebrate evolution.

Comparative Genomics and Transcriptomics To Analyze Fruiting Body Development in Filamentous Ascomycetes

Many filamentous ascomycetes develop three-dimensional fruiting bodies for production and dispersal of sexual spores. Fruiting bodies are among the most complex structures differentiated by ascomycetes; however, the molecular mechanisms underlying this process are insufficiently understood. Previous comparative transcriptomics analyses of fruiting body development in different ascomycetes suggested that there might be a core set of genes that are transcriptionally regulated in a similar manner across species. Conserved patterns of gene expression can be indicative of functional relevance, and therefore such a set of genes might constitute promising candidates for functional analyses. In this study, we have sequenced the genome of the Pezizomycete Ascodesmis nigricans, and performed comparative transcriptomics of developing fruiting bodies of this fungus, the Pezizomycete Pyronema confluens, and the Sordariomycete Sordaria macrospora With only 27 Mb, the A. nigricans genome is the smallest Pezizomycete genome sequenced to date. Comparative transcriptomics indicated that gene expression patterns in developing fruiting bodies of the three species are more similar to each other than to nonsexual hyphae of the same species. An analysis of 83 genes that are upregulated only during fruiting body development in all three species revealed 23 genes encoding proteins with predicted roles in vesicle transport, the endomembrane system, or transport across membranes, and 13 genes encoding proteins with predicted roles in chromatin organization or the regulation of gene expression. Among four genes chosen for functional analysis by deletion in S. macrospora, three were shown to be involved in fruiting body formation, including two predicted chromatin modifier genes.

Evolution of DNA Methylome Diversity in Eukaryotes

Cytosine DNA methylation (5mC) is a widespread base modification in eukaryotic genomes with critical roles in transcriptional regulation. In recent years, our understanding of 5mC has changed because of advances in 5mC detection techniques that allow mapping of this mark on the whole genome scale. Profiling DNA methylomes from organisms across the eukaryotic tree of life has reshaped our views on the evolution of 5mC. In this review, we explore the macroevolution of 5mC in major eukaryotic groups, and then focus on recent advances made in animals. Genomic 5mC patterns as well as the mechanisms of 5mC deposition tend to be evolutionary labile across large phylogenetic distances; however, some common patterns are starting to emerge. Within the animal kingdom, 5mC diversity has proven to be much greater than anticipated. For example, a previously held common view that genome hypermethylation is a trait exclusive to vertebrates has recently been challenged. Also, data from genome-wide studies are starting to yield insights into the potential roles of 5mC in invertebrate cis regulation. Here we provide an evolutionary perspective of both the well-known and enigmatic roles of 5mC across the eukaryotic tree of life.

The evolutionary landscape of the Rab family in chordates

Intracellular traffic amongst organelles represents a key feature for eukaryotes and is orchestrated principally by members of Rab family, the largest within Ras superfamily. Given that variations in Rab repertoire have been fundamental in animal diversification, we provided the most exhaustive survey regarding the Rab toolkit of chordates. Our findings reveal the existence of 42 metazoan conserved subfamilies exhibiting a univocal intron/exon structure preserved from cnidarians to vertebrates. Since the current view does not capture the Rab complexity, we propose a new Rab family classification in three distinct monophyletic clades. The Rab complement of chordates shows a dramatic diversification due to genome duplications and independent gene duplications and losses with sharp differences amongst cephalochordates, tunicates and gnathostome vertebrates. Strikingly, the analysis of the domain architecture of this family highlighted the existence of chimeric calcium-binding Rabs, which are animal novelties characterized by a complex evolutionary history in gnathostomes and whose role in cellular metabolism is obscure. This work provides novel insights in the knowledge of Rab family: our hypothesis is that chordates represent a hotspot of Rab variability, with many events of gene gains and losses impacting intracellular traffic capabilities. Our results help to elucidate the role of Rab members in the transport amongst endomembranes and shed light on intracellular traffic routes in vertebrates. Then, since the predominant role of Rabs in the molecular communication between different cellular districts, this study paves to way to comprehend inherited or acquired human disorders provoked by dysfunctions in Rab genes.

Convergent evolution of a vertebrate-like methylome in a marine sponge

Vertebrates have highly methylated genomes at CpG positions, whereas invertebrates have sparsely methylated genomes. This increase in methylation content is considered a major regulatory innovation of vertebrate genomes. However, here we report that a sponge, proposed as the potential sister group to the rest of animals, has a highly methylated genome. Despite major differences in genome size and architecture, we find similarities between the independent acquisitions of the hypermethylated state. Both lineages show genome-wide CpG depletion, conserved strong transcription factor methyl-sensitivity and developmental methylation dynamics at 5-hydroxymethylcytosine enriched regions. Together, our findings trace back patterns associated with DNA methylation in vertebrates to the early steps of animal evolution. Thus, the sponge methylome challenges previous hypotheses concerning the uniqueness of vertebrate genome hypermethylation and its implications for regulatory complexity.

Genome-Wide Identification and Transcriptomic Analysis of MicroRNAs Across Various Amphioxus Organs Using Deep Sequencing

Amphioxus is the closest living invertebrate proxy of the vertebrate ancestor. Systematic gene identification and expression profile analysis of amphioxus organs are thus important for clarifying the molecular mechanisms of organ function formation and further understanding the evolutionary origin of organs and genes in vertebrates. The precise regulation of microRNAs (miRNAs) is crucial for the functional specification and differentiation of organs. In particular, those miRNAs that are expressed specifically in organs (OSMs) play key roles in organ identity, differentiation, and function. In this study, the genome-wide miRNA transcriptome was analyzed in eight organs of adult amphioxus Branchiostoma belcheri using deep sequencing. A total of 167 known miRNAs and 23 novel miRNAs (named novel_mir), including 139 conserved miRNAs, were discovered, and 79 of these were identified as OSMs. Additionally, analyses of the expression patterns of eight randomly selected known miRNAs demonstrated the accuracy of the miRNA deep sequencing that was used in this study. Furthermore, potentially OSM-regulated genes were predicted for each organ type. Functional enrichment of these predicted targets, as well as further functional analyses of known OSMs, was conducted. We found that the OSMs were potentially to be involved in organ-specific functions, such as epidermis development, gonad development, muscle cell development, proteolysis, lipid metabolism, and generation of neurons. Moreover, OSMs with non-organ-specific functions were detected and primarily include those related to innate immunity and response to stimuli. These findings provide insights into the regulatory roles of OSMs in various amphioxus organs.

Exciting times to study the identity and evolution of cell types

The EMBO/EMBL Symposium on ‘The Identity and Evolution of Cell Types’ took place in Heidelberg, Germany, on 15-19 May 2019. The symposium, which brought together a diverse group of speakers addressing a wide range of questions in multiple model systems, provided a platform to discuss how the concept of a cell type should be considered in the era of single cell omics techniques and how cell type evolution can be studied.

Acquisition of multipotent and migratory neural crest cells in vertebrate evolution

The emergence of multipotent and migratory neural crest (NC) cells defines a key evolutionary transition from invertebrates to vertebrates. Studies in vertebrates have identified a complex gene regulatory network that governs sequential stages of NC ontogeny. Comparative analysis has revealed extensive conservation of the overall architecture of the NC gene regulatory network between jawless and jawed vertebrates. Among invertebrates, urochordates express putative NC gene homologs in the neural plate border region, but these NC-like cells do not have migratory capacity, whereas cephalochordates contain no NC cells but its genome contains most homologs of vertebrate NC genes. Whether the absence of migratory NC cells in invertebrates is due to differences in enhancer elements or an intrinsic limitation in potency remains unclear. We provide a brief overview of mechanisms that might explain how ancestral NC-like cells acquired the multipotency and migratory capacity seen in vertebrates.

Microsyntenic Clusters Reveal Conservation of lncRNAs in Chordates Despite Absence of Sequence Conservation

Homologous long non-coding RNAs (lncRNAs) are elusive to identify by sequence similarity due to their fast-evolutionary rate. Here we develop LincOFinder, a pipeline that finds conserved intergenic lncRNAs (lincRNAs) between distant related species by means of microsynteny analyses. Using this tool, we have identified 16 bona fide homologous lincRNAs between the amphioxus and human genomes. We characterized and compared in amphioxus and Xenopus the expression domain of one of them, Hotairm1, located in the anterior part of the Hox cluster. In addition, we analyzed the function of this lincRNA in Xenopus, showing that its disruption produces a severe headless phenotype, most probably by interfering with the regulation of the Hox cluster. Our results strongly suggest that this lincRNA has probably been regulating the Hox cluster since the early origin of chordates. Our work pioneers the use of syntenic searches to identify non-coding genes over long evolutionary distances and helps to further understand lncRNA evolution.

Pioneer and repressive functions of p63 during zebrafish embryonic ectoderm specification

The transcription factor p63 is a master regulator of ectoderm development. Although previous studies show that p63 triggers epidermal differentiation in vitro, the roles of p63 in developing embryos remain poorly understood. Here, we use zebrafish embryos to analyze in vivo how p63 regulates gene expression during development. We generate tp63-knock-out mutants that recapitulate human phenotypes and show down-regulated epidermal gene expression. Following p63-binding dynamics, we find two distinct functions clearly separated in space and time. During early development, p63 binds enhancers associated to neural genes, limiting Sox3 binding and reducing neural gene expression. Indeed, we show that p63 and Sox3 are co-expressed in the neural plate border. On the other hand, p63 acts as a pioneer factor by binding non-accessible chromatin at epidermal enhancers, promoting their opening and epidermal gene expression in later developmental stages. Therefore, our results suggest that p63 regulates cell fate decisions during vertebrate ectoderm specification.

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.

Genetic regulation of amphioxus somitogenesis informs the evolution of the vertebrate head mesoderm

The evolution of vertebrates from an ancestral chordate was accompanied by the acquisition of a predatory lifestyle closely associated to the origin of a novel anterior structure, the highly specialized head. While the vertebrate head mesoderm is unsegmented, the paraxial mesoderm of the earliest divergent chordate clade, the cephalochordates (amphioxus), is fully segmented in somites. We have previously shown that fibroblast growth factor signalling controls the formation of the most anterior somites in amphioxus; therefore, unravelling the fibroblast growth factor signalling downstream effectors is of crucial importance to shed light on the evolutionary origin of vertebrate head muscles. By using a comparative RNA sequencing approach and genetic functional analyses, we show that several transcription factors, such as Six1/2, Pax3/7 and Zic, act in combination to ensure the formation of three different somite populations. Interestingly, these proteins are orthologous to key regulators of trunk, and not head, muscle formation in vertebrates. Contrary to prevailing thinking, our results suggest that the vertebrate head mesoderm is of visceral and not paraxial origin and support a multistep evolutionary scenario for the appearance of the unsegmented mesoderm of the vertebrates new 'head'.

Deciphering the three-domain architecture in schlafens and the structures and roles of human schlafen12 and serpinB12 in transcriptional regulation

Schlafen proteins are important in cell differentiation and defense against viruses, and yet this family of vertebrate proteins is just beginning to be understood at the molecular level. Here, the three-dimensional architecture and molecular interfaces of human schlafen12 (hSLFN12), which promotes intestinal stem cell differentiation, are analyzed by sequence conservation and structural modeling in light of the functions of its homologs and binding partners. Our analysis shows that the schlafen or divergent AAA ATPase domain described in the N-terminal region of schlafens in databases and the literature is a misannotation. This N-terminal region is conclusively an AlbA_2 DNA/RNA binding domain, forming the conserved core of schlafens and their sequence homologs from bacteria through mammals. Group III schlafens additionally contain a AAA NTPase domain in their C-terminal helicase region. In hSLFN12, we have uncovered a domain matching rho GTPases, which directly follows the AlbA_2 domain in all group II-III schlafens. Potential roles for the GTPase-like domain include antiviral activity and cytoskeletal interactions that contribute to nucleocytoplasmic shuttling and cell polarization during differentiation. Based on features conserved with rSlfn13, the AlbA_2 region in hSLFN12 is likely to bind RNA, possibly as a ribonuclease. We hypothesize that RNA binding by hSLFN12 contributes to an RNA-induced transcriptional silencing/E3 ligase complex, given the functions of hSLFN12's partners, SUV39H1, JMJD6, and PDLIM7. hSLFN12's partner hSerpinB12 may contribute to heterochromatin formation, based on its homology to MENT, or directly regulate transcription via its binding to RNA polymerase II. The analysis presented here provides clear architectural and transcriptional regulation hypotheses to guide experimental design for hSLFN12 and the thousands of schlafens that share its motifs.

yap1b, a divergent Yap/Taz family member, cooperates with yap1 in survival and morphogenesis via common transcriptional targets

Yap1/Taz are well-known Hippo effectors triggering complex transcriptional programs controlling growth, survival and cancer progression. Here, we describe yap1b, a new Yap1/Taz family member with a unique transcriptional activation domain that cannot be phosphorylated by Src/Yes kinases. We show that yap1b evolved specifically in euteleosts (i.e. including medaka but not zebrafish) by duplication and adaptation of yap1. Using DamID-seq, we generated maps of chromatin occupancy for Yap1, Taz (Wwtr1) and Yap1b in gastrulating zebrafish and medaka embryos. Our comparative analyses uncover the genetic programs controlled by Yap family proteins during early embryogenesis, and show largely overlapping targets for Yap1 and Yap1b. CRISPR/Cas9-induced mutation of yap1b in medaka does not result in an overt phenotype during embryogenesis or adulthood. However, yap1b mutation strongly enhances the embryonic malformations observed in yap1 mutants. Thus yap1 ^-/-; yap1b ^-/- double mutants display more severe body flattening, eye misshaping and increased apoptosis than yap1 ^-/- single mutants, thus revealing overlapping gene functions. Our results indicate that, despite its divergent transactivation domain, Yap1b cooperates with Yap1 to regulate cell survival and tissue morphogenesis during early development.

TET enzymes, DNA demethylation and pluripotency

Ten-eleven translocation (TET) methylcytosine dioxygenases (TET1, TET2, TET3) actively cause demethylation of 5-methylcytosine (5mC) and produce and safeguard hypomethylation at key regulatory regions across the genome. This 5mC erasure is particularly important in pluripotent embryonic stem cells (ESCs) as they need to maintain self-renewal capabilities while retaining the potential to generate different cell types with diverse 5mC patterns. In this review, we discuss the multiple roles of TET proteins in mouse ESCs, and other vertebrate model systems, with a particular focus on TET functions in pluripotency, differentiation, and developmental DNA methylome reprogramming. Furthermore, we elaborate on the recently described non-catalytic roles of TET proteins in diverse biological contexts. Overall, TET proteins are multifunctional regulators that through both their catalytic and non-catalytic roles carry out myriad functions linked to early developmental processes.

De novo assembly of the goldfish (Carassius auratus) genome and the evolution of genes after whole-genome duplication

For over a thousand years, the common goldfish (Carassius auratus) was raised throughout Asia for food and as an ornamental pet. As a very close relative of the common carp (Cyprinus carpio), goldfish share the recent genome duplication that occurred approximately 14 million years ago in their common ancestor. The combination of centuries of breeding and a wide array of interesting body morphologies provides an exciting opportunity to link genotype to phenotype and to understand the dynamics of genome evolution and speciation. We generated a high-quality draft sequence and gene annotations of a "Wakin" goldfish using 71X PacBio long reads. The two subgenomes in goldfish retained extensive synteny and collinearity between goldfish and zebrafish. However, genes were lost quickly after the carp whole-genome duplication, and the expression of 30% of the retained duplicated gene diverged substantially across seven tissues sampled. Loss of sequence identity and/or exons determined the divergence of the expression levels across all tissues, while loss of conserved noncoding elements determined expression variance between different tissues. This assembly provides an important resource for comparative genomics and understanding the causes of goldfish variants.

Molecular and evolutionary aspects of the protochordate digestive system

The digestive system is a functional unit consisting of an endodermal tubular structure (alimentary canal) and accessory organs that function in nutrition processing in most triploblastic animals. Various morphologies and apparatuses are formed depending on the phylogenetical relationship and food habits of the specific species. Nutrition processing and morphogenesis of the alimentary canal and accessory organs have both been investigated in vertebrates, mainly humans and mammals. When attempting to understand the evolutionary processes that led to the vertebrate digestive system, however, it is useful to examine other chordates, specifically protochordates, which share fundamental functional and morphogenetic molecules with vertebrates, which also possess non-duplicated genomes. In protochordates, basic anatomical and physiological studies have mainly described the characteristic traits of suspension feeders. Recent progress in genome sequencing has allowed researchers to comprehensively detail protochordate genes and has compared the genetic backgrounds among chordate nutrition processing and alimentary canal/accessory organ systems based on genomic information. Gene expression analyses have revealed spatiotemporal gene expression profiles in protochordate alimentary canals. Additionally, to investigate the basis of morphological diversity in the chordate alimentary canal and accessory organs, evolutionary developmental research has examined developmental transcription factors related to morphogenesis and anterior-posterior pattering of the alimentary canal and accessory organs. In this review, we summarize the current knowledge of molecules involved in nutrition processing and the development of the alimentary canal and accessory organs with innate immune and endocrine roles in protochordates and we explore the molecular basis for understanding the evolution of the chordate digestive system.

Functional lability of RNA-dependent RNA polymerases in animals

RNA interference (RNAi) requires RNA-dependent RNA polymerases (RdRPs) in many eukaryotes, and RNAi amplification constitutes the only known function for eukaryotic RdRPs. Yet in animals, classical model organisms can elicit RNAi without possessing RdRPs, and only nematode RNAi was shown to require RdRPs. Here we show that RdRP genes are much more common in animals than previously thought, even in insects, where they had been assumed not to exist. RdRP genes were present in the ancestors of numerous clades, and they were subsequently lost at a high frequency. In order to probe the function of RdRPs in a deuterostome (the cephalochordate Branchiostoma lanceolatum), we performed high-throughput analyses of small RNAs from various Branchiostoma developmental stages. Our results show that Branchiostoma RdRPs do not appear to participate in RNAi: we did not detect any candidate small RNA population exhibiting classical siRNA length or sequence features. Our results show that RdRPs have been independently lost in dozens of animal clades, and even in a clade where they have been conserved (cephalochordates) their function in RNAi amplification is not preserved. Such a dramatic functional variability reveals an unexpected plasticity in RNA silencing pathways.

RNA interference (RNAi) is a conserved gene regulation system in eukaryotes. In non-animal eukaryotes, it necessitates RNA-dependent RNA polymerases (“RdRPs”). Among animals, only nematodes appear to require RdRPs for RNAi. Yet additional animal clades have RdRPs and it is assumed that they participate in RNAi. Here, we find that RdRPs are much more common in animals than previously thought, but their genes were independently lost in many lineages. Focusing on a species with RdRP genes (a cephalochordate), we found that it does not use them for RNAi. While RNAi is the only known function for eukaryotic RdRPs, our results suggest additional roles. Eukaryotic RdRPs thus have a complex evolutionary history in animals, with frequent independent losses and apparent functional diversification.

Evolution of cis-regulatory modules for the head organizer gene goosecoid in chordates: comparisons between Branchiostoma and Xenopus

BACKGROUND

In cephalochordates (amphioxus), the notochord runs along the dorsal to the anterior tip of the body. In contrast, the vertebrate head is formed anterior to the notochord, as a result of head organizer formation in anterior mesoderm during early development. A key gene for the vertebrate head organizer, goosecoid (gsc), is broadly expressed in the dorsal mesoderm of amphioxus gastrula. Amphioxus gsc expression subsequently becomes restricted to the posterior notochord from the early neurula. This has prompted the hypothesis that a change in expression patterns of gsc led to development of the vertebrate head during chordate evolution. However, molecular mechanisms of head organizer evolution involving gsc have never been elucidated.

RESULTS

To address this question, we compared cis-regulatory modules of vertebrate organizer genes between amphioxus, Branchiostoma japonicum, and frogs, Xenopus laevis and Xenopus tropicalis. Here we show conservation and diversification of gene regulatory mechanisms through cis-regulatory modules for gsc, lim1/lhx1, and chordin in Branchiostoma and Xenopus. Reporter analysis using Xenopus embryos demonstrates that activation of gsc by Nodal/FoxH1 signal through the 5' upstream region, that of lim1 by Nodal/FoxH1 signal through the first intron, and that of chordin by Lim1 through the second intron, are conserved between amphioxus and Xenopus. However, activation of gsc by Lim1 and Otx through the 5' upstream region in Xenopus are not conserved in amphioxus. Furthermore, the 5' region of amphioxus gsc recapitulated the amphioxus-like posterior mesoderm expression of the reporter gene in transgenic Xenopus embryos.

CONCLUSIONS

On the basis of this study, we propose a model, in which the gsc gene acquired the cis-regulatory module bound with Lim1 and Otx at its 5' upstream region to be activated persistently in anterior mesoderm, in the vertebrate lineage. Because Gsc globally represses trunk (notochord) genes in the vertebrate head organizer, this cooption of gsc in vertebrates appears to have resulted in inhibition of trunk genes and acquisition of the head organizer and its derivative prechordal plate.

Recapitulation-like developmental transitions of chromatin accessibility in vertebrates

The relationship between development and evolution has been a central theme in evolutionary developmental biology. Across the vertebrates, the most highly conserved gene expression profiles are found at mid-embryonic, organogenesis stages, whereas those at earlier and later stages are more diverged. This hourglass-like pattern of divergence does not necessarily rule out the possibility that gene expression profiles that are more evolutionarily derived appear at later stages of development; however, no molecular-level evidence of such a phenomenon has been reported. To address this issue, we compared putative gene regulatory elements among different species within a phylum. We made a genome-wide assessment of accessible chromatin regions throughout embryogenesis in three vertebrate species (mouse, chicken, and medaka) and estimated the evolutionary ages of these regions to define their evolutionary origins on the phylogenetic tree. In all the three species, we found that genomic regions tend to become accessible in an order that parallels their phylogenetic history, with evolutionarily newer gene regulations activated at later developmental stages. This tendency was restricted only after the mid-embryonic, phylotypic periods. Our results imply a phylogenetic hierarchy of putative regulatory regions, in which their activation parallels the phylogenetic order of their appearance. One evolutionary mechanism that may explain this phenomenon is that newly introduced regulatory elements are more likely to survive if activated at later stages of embryogenesis. Possible relationships between this phenomenon and the so-called recapitulation are discussed.