A pre-vertebrate endodermal origin of calcitonin-producing neuroendocrine cells.

Vertebrate calcitonin-producing cells (C-cells) are neuroendocrine cells that secrete the small peptide hormone calcitonin in response to elevated blood calcium levels. Whereas mouse C-cells reside within the thyroid gland and derive from pharyngeal endoderm, avian C-cells are located within ultimobranchial glands and have been reported to derive from the neural crest. We use a comparative cell lineage tracing approach in a range of vertebrate model systems to resolve the ancestral embryonic origin of vertebrate C-cells. We find, contrary to previous studies, that chick C-cells derive from pharyngeal endoderm, with neural crest-derived cells instead contributing to connective tissue intimately associated with C-cells in the ultimobranchial gland. This endodermal origin of C-cells is conserved in a ray-finned bony fish (zebrafish) and a cartilaginous fish (the little skate, Leucoraja erinacea). Furthermore, we discover putative C-cell homologs within the endodermally-derived pharyngeal epithelium of the ascidian Ciona intestinalis and the amphioxus Branchiostoma lanceolatum, two invertebrate chordates that lack neural crest cells. Our findings point to a conserved endodermal origin of C-cells across vertebrates and to a pre-vertebrate origin of this cell type along the chordate stem.

Conservation of cis-Regulatory Syntax Underlying Deuterostome Gastrulation.

Throughout embryonic development, the shaping of the functional and morphological characteristics of embryos is orchestrated by an intricate interaction between transcription factors and cis-regulatory elements. In this study, we conducted a comprehensive analysis of deuterostome cis-regulatory landscapes during gastrulation, focusing on four paradigmatic species: the echinoderm Strongylocentrotus purpuratus, the cephalochordate Branchiostoma lanceolatum, the urochordate Ciona intestinalis, and the vertebrate Danio rerio. Our approach involved comparative computational analysis of ATAC-seq datasets to explore the genome-wide blueprint of conserved transcription factor binding motifs underlying gastrulation. We identified a core set of conserved DNA binding motifs associated with 62 known transcription factors, indicating the remarkable conservation of the gastrulation regulatory landscape across deuterostomes. Our findings offer valuable insights into the evolutionary molecular dynamics of embryonic development, shedding light on conserved regulatory subprograms and providing a comprehensive perspective on the conservation and divergence of gene regulation underlying the gastrulation process.

Ribosomal protein L17 functions as an antimicrobial protein in amphioxus.

Antimicrobial peptides (AMPs), characterized by their cationic nature and amphiphilic properties, play a pivotal role in inhibiting the biological activity of microbes. Currently, only a fraction of the antimicrobial potential within the ribosomal protein family has been explored, despite its extensive membership and resemblance to AMPs. Herein we demonstrated that amphioxus RPL17 (BjRPL17) exhibited not only upregulated expression upon bacterial stimulation but also possessed bactericidal capabilities against both Gram-negative and -positive bacteria through combined action mechanisms including interaction with cell surface molecules LPS, LTA, and PGN, disruption of cell membrane integrity, promotion of membrane depolarization, and induction of intracellular ROS production. Furthermore, a peptide derived from residues 127-141 of BjRPL17 (termed BjRPL17-1) showed antibacterial activity against Staphylococcus aureus and its methicillin-resistant strain via the same mechanism observed for the full-length protein. Additionally, the rpl17 gene was highly conserved in Metazoa, hinting it may play a universal role in the antibacterial defense system in different animals. Importantly, neither BjRPL17 nor peptide BjRPL17-1 exhibited toxicity towards mammalian cells thereby offering prospects for designing novel AMP agents based on these findings. Collectively, our results establish RPL17 as a novel member of AMPs with remarkable evolutionary conservation.

Functional characterization of growth hormone releasing hormone and its receptor in amphioxus with implication for origin of hypothalamic-pituitary axis.

Growth hormone-releasing hormone (GHRH) has been widely shown to stimulate growth hormone (GH) production via binding to GHRH receptor GHRHR in various species of vertebrates, but information regarding the functional roles of GHRH and GHRHR in the protochordate amphioxus remains rather scarce. We showed here that two mature peptides, BjGHRH-1 and BjGHRH-2, encoded by BjGHRH precursor, and a single BjGHRHR protein were identified in the amphioxus Branchiostoma. japonicum. Like the distribution profiles of vertebrate GHRHs and GHRHRs, both the genes Bjghrh and Bjghrhr were widely expressed in the different tissues of amphioxus, including in the cerebral vesicle, Hatschek's pit, neural tube, gill, hepatic caecum, notochord, testis and ovary. Moreover, both BjGHRH-1 and BjGHRH-2 interacted with BjGHRHR, and triggered the cAMP/PKA signal pathway in a dose-dependent manner. Importantly, BjGHRH-1 and BjGHRH-2 were both able to activate the expression of GH-like gene in the cells of Hatschek's pit. These indicate that a functional vertebrate-like GHRH-GHRHR axis had already emerged in amphioxus, which is a seminal innovation making physiological divergence including reproduction, growth, metabolism, stress and osmoregulation possible during the early evolution of vertebrates.

An amphioxus neurula stage cell atlas supports a complex scenario for the emergence of vertebrate head mesoderm.

The emergence of new structures can often be linked to the evolution of novel cell types that follows the rewiring of developmental gene regulatory subnetworks. Vertebrates are characterized by a complex body plan compared to the other chordate clades and the question remains of whether and how the emergence of vertebrate morphological innovations can be related to the appearance of new embryonic cell populations. We previously proposed, by studying mesoderm development in the cephalochordate amphioxus, a scenario for the evolution of the vertebrate head mesoderm. To further test this scenario at the cell population level, we used scRNA-seq to construct a cell atlas of the amphioxus neurula, stage at which the main mesodermal compartments are specified. Our data allowed us to validate the presence of a prechordal-plate like territory in amphioxus. Additionally, the transcriptomic profile of somite cell populations supports the homology between specific territories of amphioxus somites and vertebrate cranial/pharyngeal and lateral plate mesoderm. Finally, our work provides evidence that the appearance of the specific mesodermal structures of the vertebrate head was associated to both segregation of pre-existing cell populations, and co-option of new genes for the control of myogenesis.

A novel C-type lectin protein (BjCTL5) interacts with apoptosis stimulating proteins of p53 (ASPP) to activate NF-κB signaling pathway in primitive chordate.

C-type lectin proteins (CTLs), a group of pattern recognition receptors (PRRs), play pivotal roles in immune responses. However, the signal transduction and regulation of CTLs in cephalochordates have yet to be explored. In this study, we examined the composition of CTLs in Branchiostoma japonicum, identifying a total of 272 CTLs. These CTLs underwent further analysis concerning domain arrangement, tandem and segmental duplication events. A multidomain C-type lectin gene, designated as BjCTL5, encompassing CLECT, KR, CUB, MAM, and SR domains, was the focal point of our investigation. BjCTL5 exhibits ubiquitous expression across all detected tissues and is responsive to stimulation by LPS, mannose, and poly (I:C). The recombinant protein of BjCTL5 can bind to Escherichia coli and Staphylococcus aureus, inducing their agglutination and inhibiting the proliferation of S. aureus. Yeast two-hybrid, CoIP, and confocal immunofluorescence experiments revealed the interaction between BjCTL5 and apoptosis-stimulating proteins of p53, BjASPP. Intriguingly, BjCTL5 was observed to induce the luciferase activity of the NF-κB promoter in HEK293T cells. These results suggested a potential interaction between BjCTL5 and BjASPP, implicating that they involve in the activation of the NF-κB signaling pathway, which provides an evolutionary viewpoint on NF-κB signaling pathway in primitive chordate.

Evolution of the gene regulatory network of body axis by enhancer hijacking in amphioxus.

A central goal of evolutionary developmental biology is to decipher the evolutionary pattern of gene regulatory networks (GRNs) that control embryonic development, and the mechanism underlying GRNs evolution. The Nodal signaling that governs the body axes of deuterostomes exhibits a conserved GRN orchestrated principally by Nodal, Gdf1/3, and Lefty. Here we show that this GRN has been rewired in cephalochordate amphioxus. We found that while the amphioxus Gdf1/3 ortholog exhibited nearly no embryonic expression, its duplicate Gdf1/3-like, linked to Lefty, was zygotically expressed in a similar pattern as Lefty. Consistent with this, while Gdf1/3-like mutants showed defects in axial development, Gdf1/3 mutants did not. Further transgenic analyses showed that the intergenic region between Gdf1/3-like and Lefty could drive reporter gene expression as that of the two genes. These results indicated that Gdf1/3-like has taken over the axial development role of Gdf1/3 in amphioxus, possibly through hijacking Lefty enhancers. We finally demonstrated that, to compensate for the loss of maternal Gdf1/3 expression, Nodal has become an indispensable maternal factor in amphioxus and its maternal mutants caused axial defects as Gdf1/3-like mutants. We therefore demonstrated a case that the evolution of GRNs could be triggered by enhancer hijacking events. This pivotal event has allowed the emergence of a new GRN in extant amphioxus, presumably through a stepwise process. In addition, the co-expression of Gdf1/3-like and Lefty achieved by a shared regulatory region may have provided robustness during body axis formation, which provides a selection-based hypothesis for the phenomena called developmental system drift.

Comparison of Evolutionary Relationships between Branchiostoma floridae, Ciona intestinalis, and Homo sapiens Globins Provide Evidence of Gene Co-Option and Convergent Evolution.

Globins have been studied as model proteins to elucidate the principles of protein evolution. This was achieved by understanding the relationship between amino acid sequence, three-dimensional structure, physicochemical properties, and physiological function. Previous molecular phylogenies of chordate globin genes revealed the monophyletic evolution of urochordate globins and suggested convergent evolution. However, to provide evidence of convergent evolution, it is necessary to determine the physicochemical and functional similarities between vertebrates and urochordate globins. In this study, we determined the expression patterns of Ciona globin genes using real-time RT-PCR. Two genes (Gb-1 and Gb-2) were predominantly expressed in the branchial sac, heart, and hemocytes and were induced under hypoxia. Combined with the sequence analysis, our findings suggest that Gb-1/-2 correspond to vertebrate hemoglobin-α/-ß. However, we did not find a robust similarity between Gb-3, Gb-4, and vertebrate globins. These results suggested that, even though Ciona globins obtained their unique functions differently from vertebrate globins, the two of them shared some physicochemical features and physiological functions. Our findings offer a good example for understanding the molecular mechanisms underlying gene co-option and convergence, which could lead to evolutionary innovations.

Amphioxus as a model to study the evolution of development in chordates.

Cephalochordates and tunicates represent the only two groups of invertebrate chordates, and extant cephalochordates - commonly known as amphioxus or lancelets - are considered the best proxy for the chordate ancestor, from which they split around 520 million years ago. Amphioxus has been an important organism in the fields of zoology and embryology since the 18th century, and the morphological and genomic simplicity of cephalochordates (compared to vertebrates) makes amphioxus an attractive model for studying chordate biology at the cellular and molecular levels. Here we describe the life cycle of amphioxus, and discuss the natural histories and habitats of the different species of amphioxus. We also describe their use as laboratory animal models, and discuss the techniques that have been developed to study different aspects of amphioxus.

Amphioxus adenosine-to-inosine tRNA-editing enzyme that can perform C-to-U and A-to-I deamination of DNA.

Adenosine-to-inosine tRNA-editing enzyme has been identified for more than two decades, but the study on its DNA editing activity is rather scarce. We show that amphioxus (Branchiostoma japonicum) ADAT2 (BjADAT2) contains the active site 'HxE-PCxxC' and the key residues for target-base-binding, and amphioxus ADAT3 (BjADAT3) harbors both the N-terminal positively charged region and the C-terminal pseudo-catalytic domain important for recognition of substrates. The sequencing of BjADAT2-transformed Escherichia coli genome suggests that BjADAT2 has the potential to target E. coli DNA and can deaminate at TCG and GAA sites in the E. coli genome. Biochemical analyses further demonstrate that BjADAT2, in complex with BjADAT3, can perform A-to-I editing of tRNA and convert C-to-U and A-to-I deamination of DNA. We also show that BjADAT2 preferentially deaminates adenosines and cytidines in the loop of DNA hairpin structures of substrates, and BjADAT3 also affects the type of DNA substrate targeted by BjADAT2. Finally, we find that C89, N113, C148 and Y156 play critical roles in the DNA editing activity of BjADAT2. Collectively, our study indicates that BjADAT2/3 is the sole naturally occurring deaminase with both tRNA and DNA editing capacity identified so far in Metazoa.

Characterization of GSDME in amphioxus provides insights into the functional evolution of GSDM-mediated pyroptosis.

Members of the gasdermin (GSDM) family are pore-forming effectors that cause membrane permeabilization and pyroptosis, a lytic proinflammatory type of cell death. To reveal the functional evolution of GSDM-mediated pyroptosis at the transition from invertebrates to vertebrates, we conducted functional characterization of amphioxus GSDME (BbGSDME) and found that it can be cleaved by distinct caspase homologs, yielding the N253 and N304 termini with distinct functions. The N253 fragment binds to cell membrane, triggers pyroptosis, and inhibits bacterial growth, while the N304 performs negative regulation of N253-mediated cell death. Moreover, BbGSDME is associated with bacteria-induced tissue necrosis and transcriptionally regulated by BbIRF1/8 in amphioxus. Interestingly, several amino acids that are evolutionarily conserved were found to be important for the function of both BbGSDME and HsGSDME, shedding new lights on the functional regulation of GSDM-mediated inflammation.

Genomic Analysis of Amphioxus Reveals a Wide Range of Fragments Homologous to Viral Sequences.

Amphioxus species are considered living fossils and are important in the evolutionary study of chordates and vertebrates. To explore viral homologous sequences, a high-quality annotated genome of the Beihai amphioxus (Branchiostoma belcheri beihai) was examined using virus sequence queries. In this study, 347 homologous fragments (HFs) of viruses were identified in the genome of B. belcheri beihai, of which most were observed on 21 genome assembly scaffolds. HFs were preferentially located within protein-coding genes, particularly in their CDS regions and promoters. A range of amphioxus genes with a high frequency of HFs is proposed, including histone-related genes that are homologous to the Histone or Histone H2B domains of viruses. Together, this comprehensive analysis of viral HFs provides insights into the neglected role of viral integration in the evolution of amphioxus.

Identification of a novel amphioxus leucine-rich repeat receptor involved in phagocytosis reveals a role for Slit2-N-type LRR in bacterial elimination.

The basal chordate amphioxus is a model for tracing the origin and evolution of vertebrate immunity. To explore the evolution of immunoreceptor signaling pathways, we searched the associated receptors of the amphioxus Branchiostoma belcheri (Bb) homolog of immunoreceptor signaling adaptor protein Grb2. Mass-spectrum analysis of BbGrb2 immunoprecipitates from B. belcheri intestine lysates revealed a folate receptor (FR) domain- and leucine-rich repeat (LRR)-containing protein (FrLRR). Sequence and structural analysis showed that FrLRR is a membrane protein with a predicted curved solenoid structure. The N-terminal Fr domain contains very few folate-binding sites; the following LRR region is a Slit2-type LRR, and a GPI-anchored site was predicted at the C-terminus. RT-PCR analysis showed FrLRR is a transcription-mediated fusion gene of BbFR-like and BbSlit2-N-like genes. Genomic DNA structure analysis implied the B. belcheri FrLRR gene locus and the corresponding locus in Branchiostoma floridae might be generated by exon shuffling of a Slit2-N-like gene into an FR gene. RT-qPCR, immunostaining, and immunoblot results showed that FrLRR was primarily distributed in B. belcheri intestinal tissue. We further demonstrated that FrLRR localized to the cell membrane and lysosomes. Functionally, FrLRR mediated and promoted bacteria-binding and phagocytosis, and FrLRR antibody blocking or Grb2 knockdown inhibited FrLRR-mediated phagocytosis. Interestingly, we found that human Slit2-N (hSlit2-N) also mediated direct bacteria-binding and phagocytosis which was inhibited by Slit2-N antibody blocking or Grb2 knockdown. Together, these results indicate FrLRR and hSlit2-N may function as phagocytotic-receptors to promote phagocytosis through Grb2, implying the Slit2-N-type-LRR-containing proteins play a role in bacterial binding and elimination.

SUMO and PIAS repress NF-κB activation in a basal chordate.

Small ubiquitin-like modifier (SUMO) regulates various biological processes, including the MyD88/TICAMs-IRAKs-TRAF6-NF-κB pathway, one of the core immune pathways. However, its functions are inconsistent between invertebrates and vertebrates and have rarely been investigated in lower chordates, including amphioxus and fishes. Here, we investigated the SUMOylation gene system in the amphioxus, a living basal chordate. We found that amphioxus has a SUMOylation system that has a complete set of genes and preserves several ancestral traits. We proceeded to study their molecular functions using the mammal cell lines. Both amphioxus SUMO1 and SUMO2 were shown to be able to attach to NF-κB Rel and to inhibit NF-κB activation by 50-75% in a dose-dependent fashion. The inhibition by SUMO2 could be further enhanced by the addition of the SUMO E2 ligase UBC9. In comparison, while human SUMO2 inhibited RelA, human SUMO1 slightly activated RelA. We also showed that, similar to human PIAS1-4, amphioxus PIAS could serve as a SUMO E3 ligase and promote its self-SUMOylation. This suggests that amphioxus PIAS is functionally compatible in human cells. Moreover, we showed that amphioxus PIAS is not only able to inhibit NF-κB activation induced by MyD88, TICAM-like, TRAF6 and IRAK4 but also able to suppress NF-κB Rel completely in the presence of SUMO1/2 in a dose-insensitive manner. This suggests that PIAS could effectively block Rel by promoting Rel SUMOylation. In comparison, in humans, only PIAS3, but not PIAS1/2/4, has been reported to promote NF-κB SUMOylation. Taken together, the findings from amphioxus, together with those from mammals and other species, not only offer insights into the functional volatility of the animal SUMO system, but also shed light on its evolutionary transitions from amphioxus to fish, and ultimately to humans.

Retinoic Acid and POU Genes in Developing Amphioxus: A Focus on Neural Development.

POU genes are a family of evolutionarily conserved transcription factors with key functions in cell type specification and neurogenesis. In vitro experiments have indicated that the expression of some POU genes is controlled by the intercellular signaling molecule retinoic acid (RA). In this work, we aimed to characterize the roles of RA signaling in the regulation of POU genes in vivo. To do so, we studied POU genes during the development of the cephalochordate amphioxus, an animal model crucial for understanding the evolutionary origins of vertebrates. The expression patterns of amphioxus POU genes were assessed at different developmental stages by chromogenic in situ hybridization and hybridization chain reaction. Expression was further assessed in embryos subjected to pharmacological manipulation of endogenous RA signaling activity. In addition to a detailed description of the effects of these treatments on amphioxus POU gene expression, our survey included the first description of Pou2 and Pou6 expression in amphioxus embryos. We found that Pit-1, Pou2, Pou3l, and Pou6 expression are not affected by alterations of endogenous RA signaling levels. In contrast, our experiments indicated that Brn1/2/4 and Pou4 expression are regulated by RA signaling in the endoderm and the nerve cord, respectively. The effects of the treatments on Pou4 expression in the nerve cord revealed that, in developing amphioxus, RA signaling plays a dual role by (1) providing anteroposterior patterning information to neural cells and (2) specifying neural cell types. This finding is coherent with a terminal selector function of Pou4 for GABAergic neurons in amphioxus and represents the first description of RA-induced changes in POU gene expression in vivo.

Identification of Nodal-dependent enhancer of amphioxus Chordin sufficient to drive gene expression into the chordate dorsal organizer.

The core molecular mechanisms of dorsal organizer formation during gastrulation are highly conserved within the chordate lineage. One of the key characteristics is that Nodal signaling is required for the organizer-specific gene expression. This feature appears to be ancestral, as evidenced by the presence in the most basally divergent chordate amphioxus. To provide a better understanding of the evolution of organizer-specific gene regulation in chordates, we analyzed the cis-regulatory sequence of amphioxus Chordin in the context of the vertebrate embryo. First, we generated stable zebrafish transgenic lines, and by using light-sheet fluorescent microscopy, characterized in detail the expression pattern of GFP driven by the cis-regulatory sequences of amphioxus Chordin. Next, we performed a 5'deletion analysis and identified an enhancer sufficient to drive the expression of the reporter gene into a chordate dorsal organizer. Finally, we found that the identified enhancer element strongly depends on Nodal signaling, which is consistent with the well-established role of this pathway in the regulation of the expression of dorsal organizer-specific genes across chordates. The enhancer identified in our study may represent a suitable simple system to study the interplay of the evolutionarily conserved regulatory mechanisms operating during early chordate development.

Parallel evolution of amphioxus and vertebrate small-scale gene duplications.

BACKGROUND: Amphioxus are non-vertebrate chordates characterized by a slow morphological and molecular evolution. They share the basic chordate body-plan and genome organization with vertebrates but lack their 2R whole-genome duplications and their developmental complexity. For these reasons, amphioxus are frequently used as an outgroup to study vertebrate genome evolution and Evo-Devo. Aside from whole-genome duplications, genes continuously duplicate on a smaller scale. Small-scale duplicated genes can be found in both amphioxus and vertebrate genomes, while only the vertebrate genomes have duplicated genes product of their 2R whole-genome duplications. Here, we explore the history of small-scale gene duplications in the amphioxus lineage and compare it to small- and large-scale gene duplication history in vertebrates. RESULTS: We present a study of the European amphioxus (Branchiostoma lanceolatum) gene duplications thanks to a new, high-quality genome reference. We find that, despite its overall slow molecular evolution, the amphioxus lineage has had a history of small-scale duplications similar to the one observed in vertebrates. We find parallel gene duplication profiles between amphioxus and vertebrates and conserved functional constraints in gene duplication. Moreover, amphioxus gene duplicates show levels of expression and patterns of functional specialization similar to the ones observed in vertebrate duplicated genes. We also find strong conservation of gene synteny between two distant amphioxus species, B. lanceolatum and B. floridae, with two major chromosomal rearrangements. CONCLUSIONS: In contrast to their slower molecular and morphological evolution, amphioxus' small-scale gene duplication history resembles that of the vertebrate lineage both in quantitative and in functional terms.

Localization of trypsin, chymotrypsin and elastase in the digestive tract of amphioxus Branchiostoma japonicum with implications to the origin of vertebrate pancreas.

The organ pancreas is characteristic of all vertebrates, but its evolutionary origin remains largely unknown. The serine proteinases trypsin, chymotrypsin and elastase are produced primarily in the pancreas in vertebrates. Here we clearly demonstrate by enzymatic activity assay, histochemical staining and in situ hybridization that trypsin, chymotrypsin and elastase are widely distributed in the digestive tract of the amphioxus Branchiostoma japonicum, especially in the hepatic caecum and mid-gut. This suggests that an extensive region of the amphioxus digestive tract including the hepatic caecum is homologous to the vertebrate exocrine pancreatic cells, providing a new angle for the study of the origin and evolution of vertebrate pancreas.

Joint profiling of gene expression and chromatin accessibility during amphioxus development at single-cell resolution.

Vertebrate evolution was accompanied by two rounds of whole-genome duplication followed by functional divergence in terms of regulatory circuits and gene expression patterns. As a basal and slow-evolving chordate species, amphioxus is an ideal paradigm for exploring the origin and evolution of vertebrates. Single-cell sequencing has been widely used to construct the developmental cell atlas of several representative species of vertebrates (human, mouse, zebrafish, and frog) and tunicates (sea squirts). Here, we perform single-nucleus RNA sequencing (snRNA-seq) and single-cell assay for transposase accessible chromatin sequencing (scATAC-seq) for different stages of amphioxus (covering embryogenesis and adult tissues). With the datasets generated, we constructed a developmental tree for amphioxus cell fate commitment and lineage specification and characterize the underlying key regulators and genetic regulatory networks. The data are publicly available on the online platform AmphioxusAtlas.

Highly distinct genetic programs for peripheral nervous system formation in chordates.

BACKGROUND: Vertebrates develop their peripheral nervous system (PNS) from transient unique embryonic structures, the neural crest, and the ectodermal placodes that are located at the border of the forming central nervous system. By contrast, in the invertebrate chordates, amphioxus and ascidians, a large part of the PNS originates at the opposite of the embryo, in the ventral ectoderm. In both groups, a biphasic mechanism regulates ventral PNS formation: high BMP levels specify a neurogenic territory within which glutamatergic epidermal sensory neuron formation is controlled by the Notch pathway. Given these similarities and the phylogenetic relationships within chordates, it is likely that ventral PNS is an ancestral feature in chordates and that it has been lost in vertebrates. RESULTS: In order to get insights into the molecular control of ventral PNS formation and to test the hypothesis of their homology and potential contribution to the emergence of vertebrate PNS, we undertook a close comparison of ventral PNS formation in the ascidian Phallusia mammillata and the amphioxus Branchiostoma lanceolatum. Using timed RNA-seq series, we identified novel markers of the ventral PNS during different phases of its development in both species. By extensively determining the expression of paralogous and orthologous genes, we observed that only a minority of genes have a shared expression in the ventral PNS. However, a large fraction of ventral PNS orthologous genes are expressed in the dorsally forming PNS of vertebrates. CONCLUSIONS: Our work has significantly increased the molecular characterization of ventral PNS formation in invertebrate chordates. The low observed conservation of gene expression in the ventral PNS suggests that the amphioxus and ascidian ventral PNS are either not homologous, or alternatively extensive drift has occurred in their regulatory mechanisms following a long period (600 My) of separate evolution and accelerated evolution in the ascidian lineage. The homology to genes expressed in the dorsally forming PNS of vertebrates suggests that ancestral sensory neurons gene networks have been redeployed in vertebrates.