Characterization of trans-spliced chimeric RNAs: insights into the mechanism of trans-splicing.

Trans-splicing is a post-transcriptional processing event that joins exons from separate RNAs to produce a chimeric RNA. However, the detailed mechanism of trans-splicing remains poorly understood. Here, we characterize trans-spliced genes and provide insights into the mechanism of trans-splicing in the tunicate Ciona. Tunicates are the closest invertebrates to humans, and their genes frequently undergo trans-splicing. Our analysis revealed that, in genes that give rise to both trans-spliced and non-trans-spliced messenger RNAs, trans-splice acceptor sites were preferentially located at the first functional acceptor site, and their paired donor sites were weak in both Ciona and humans. Additionally, we found that Ciona trans-spliced genes had GU- and AU-rich 5' transcribed regions. Our data and findings not only are useful for Ciona research community, but may also aid in a better understanding of the trans-splicing mechanism, potentially advancing the development of gene therapy based on trans-splicing.

Retinal Cone Mosaic in sws1-Mutant Medaka (Oryzias latipes), A Teleost.

Purpose: Ablation of short single cones (SSCs) expressing short-wavelength-sensitive opsin (SWS1) is well analyzed in the field of regenerative retinal cells. In contrast with ablation studies, the phenomena caused by the complete deletion of SWS1 are less well-understood. To assess the effects of SWS1 deficiency on retinal structure, we established and analyzed sws1-mutant medaka. Methods: To visualize SWS1, a monoclonal anti-SWS1 antibody and transgenic reporter fish (Tg(sws1:mem-egfp)) were generated. We also developed a CRISPR/Cas-driven sws1-mutant line. Retinal structure of sws1 mutant was visualized using anti-SWS1, 1D4, and ZPR1 antibodies and coumarin derivatives and compared with wild type, Tg(sws1:mem-egfp), and another opsin (lws) mutant. Results: Our rat monoclonal antibody specifically recognized medaka SWS1. Sws1 mutant retained regularly arranged cone mosaic as lws mutant and its SSCs had neither SWS1 nor long wavelength sensitive opsin. Depletion of sws1 did not affect the expression of long wavelength sensitive opsin, and vice versa. ZPR1 antibody recognized arrestin spread throughout double cones and long single cones in wild-type, transgenic, and sws1-mutant lines. Conclusions: Comparative observation of sws1-mutant and wild-type retinas revealed that ZPR1 negativity is not a marker for SSCs with SWS1, but SSCs themselves. Loss of functional sws1 did not cause retinal degeneration, indicating that sws1 is not essential for cone mosaic development in medaka. Our two fish lines, one with visualized SWS1 and the other lacking functional SWS1, offer an opportunity to study neural network synapsing with SSCs and to clarify the role of SWS1 in vision.

Neuronal identities derived by misexpression of the POU IV sensory determinant in a protovertebrate.

The protovertebrate Ciona intestinalis type A (sometimes called Ciona robusta) contains a series of sensory cell types distributed across the head-tail axis of swimming tadpoles. They arise from lateral regions of the neural plate that exhibit properties of vertebrate placodes and neural crest. The sensory determinant POU IV/Brn3 is known to work in concert with regional determinants, such as Foxg and Neurogenin, to produce palp sensory cells (PSCs) and bipolar tail neurons (BTNs), in head and tail regions, respectively. A combination of single-cell RNA-sequencing (scRNA-seq) assays, computational analysis, and experimental manipulations suggests that misexpression of POU IV results in variable transformations of epidermal cells into hybrid sensory cell types, including those exhibiting properties of both PSCs and BTNs. Hybrid properties are due to coexpression of Foxg and Neurogenin that is triggered by an unexpected POU IV feedback loop. Hybrid cells were also found to express a synthetic gene battery that is not coexpressed in any known cell type. We discuss these results with respect to the opportunities and challenges of reprogramming cell types through the targeted misexpression of cellular determinants.

A single motor neuron determines the rhythm of early motor behavior in Ciona.

Recent work in tunicate supports the similarity between the motor circuits of vertebrates and basal deuterostome lineages. To understand how the rhythmic activity in motor circuits is acquired during development of protochordate Ciona, we investigated the coordination of the motor response by identifying a single pair of oscillatory motor neurons (MN2/A10.64). The MN2 neurons had Ca2+ oscillation with an ~80-s interval that was cell autonomous even in a dissociated single cell. The Ca2+ oscillation of MN2 coincided with the early tail flick (ETF). The spikes of the membrane potential in MN2 gradually correlated with the rhythm of ipsilateral muscle contractions in ETFs. The optogenetic experiments indicated that MN2 is a necessary and sufficient component of ETFs. These results indicate that MN2 is indispensable for the early spontaneous rhythmic motor behavior of Ciona. Our findings shed light on the understanding of development and evolution of chordate rhythmical locomotion.

Evolution of Developmental Programs for the Midline Structures in Chordates: Insights From Gene Regulation in the Floor Plate and Hypochord Homologues of Ciona Embryos.

In vertebrate embryos, dorsal midline tissues, including the notochord, the prechordal plate, and the floor plate, play important roles in patterning of the central nervous system, somites, and endodermal tissues by producing extracellular signaling molecules, such as Sonic hedgehog (Shh). In Ciona, hedgehog.b, one of the two hedgehog genes, is expressed in the floor plate of the embryonic neural tube, while none of the hedgehog genes are expressed in the notochord. We have identified a cis-regulatory region of hedgehog.b that was sufficient to drive a reporter gene expression in the floor plate. The hedgehog.b cis-regulatory region also drove ectopic expression of the reporter gene in the endodermal strand, suggesting that the floor plate and the endodermal strand share a part of their gene regulatory programs. The endodermal strand occupies the same topographic position of the embryo as does the vertebrate hypochord, which consists of a row of single cells lined up immediately ventral to the notochord. The hypochord shares expression of several genes with the floor plate, including Shh and FoxA, and play a role in dorsal aorta development. Whole-embryo single-cell transcriptome analysis identified a number of genes specifically expressed in both the floor plate and the endodermal strand in Ciona tailbud embryos. A Ciona FoxA ortholog FoxA.a is shown to be a candidate transcriptional activator for the midline gene battery. The present findings suggest an ancient evolutionary origin of a common developmental program for the midline structures in Olfactores.

The complete cell lineage and MAPK- and Otx-dependent specification of the dopaminergic cells in the Ciona brain.

The Ciona larva has served as a unique model for understanding the development of dopaminergic cells at single-cell resolution owing to the exceptionally small number of neurons in its brain and its fixed cell lineage during embryogenesis. A recent study suggested that the transcription factors Fer2 and Meis directly regulate the dopamine synthesis genes in Ciona, but the dopaminergic cell lineage and the gene regulatory networks that control the development of dopaminergic cells have not been fully elucidated. Here, we reveal that the dopaminergic cells in Ciona are derived from a bilateral pair of cells called a9.37 cells at the center of the neural plate. The a9.37 cells divide along the anterior-posterior axis, and all of the descendants of the posterior daughter cells differentiate into the dopaminergic cells. We show that the MAPK pathway and the transcription factor Otx are required for the expression of Fer2 in the dopaminergic cell lineage. Our findings establish the cellular and molecular framework for fully understanding the commitment to dopaminergic cells in the simple chordate brain.

Cellular identity and Ca2+ signaling activity of the non-reproductive GnRH system in the Ciona intestinalis type A (Ciona robusta) larva.

Tunicate larvae have a non-reproductive gonadotropin-releasing hormone (GnRH) system with multiple ligands and receptor heterodimerization enabling complex regulation. In Ciona intestinalis type A larvae, one of the gnrh genes, gnrh2, is conspicuously expressed in the motor ganglion and nerve cord, which are homologous structures to the hindbrain and spinal cord, respectively, of vertebrates. The gnrh2 gene is also expressed in the proto-placodal sensory neurons, which are the proposed homologue of vertebrate olfactory neurons. Tunicate larvae occupy a non-reproductive dispersal stage, yet the role of their GnRH system remains elusive. In this study, we investigated neuronal types of gnrh2-expressing cells in Ciona larvae and visualized the activity of these cells by fluorescence imaging using a calcium sensor protein. Some cholinergic neurons and dopaminergic cells express gnrh2, suggesting that GnRH plays a role in controlling swimming behavior. However, none of the gnrh2-expressing cells overlap with glycinergic or GABAergic neurons. A role in motor control is also suggested by a relationship between the activity of gnrh2-expressing cells and tail movements. Interestingly, gnrh2-positive ependymal cells in the nerve cord, known as a kind of glia cells, actively produced Ca2+ transients, suggesting that active intercellular signaling occurs in the glia cells of the nerve cord.

Spatio-temporal regulation of Rx and mitotic patterns shape the eye-cup of the photoreceptor cells in Ciona.

The ascidian larva has a pigmented ocellus comprised of a cup-shaped array of approximately 30 photoreceptor cells, a pigment cell, and three lens cells. Morphological, physiological and molecular evidence has suggested evolutionary kinship between the ascidian larval photoreceptors and vertebrate retinal and/or pineal photoreceptors. Rx, an essential factor for vertebrate photoreceptor development, has also been suggested to be involved in the development of the ascidian photoreceptor cells, but a recent revision of the photoreceptor cell lineage raised a crucial discrepancy between the reported expression patterns of Rx and the cell lineage. Here, we report spatio-temporal expression patterns of Rx at single-cell resolution along with mitotic patterns up to the final division of the photoreceptor-lineage cells in Ciona. The expression of Rx commences in non-photoreceptor a-lineage cells on the right side of the anterior sensory vesicle at the early tailbud stage. At the mid tailbud stage, Rx begins to be expressed in the A-lineage photoreceptor cell progenitors located on the right side of the posterior sensory vesicle. Thus, Rx is specifically but not exclusively expressed in the photoreceptor-lineage cells in the ascidian embryo. Two cis-regulatory modules are shown to be important for the photoreceptor-lineage expression of Rx. The cell division patterns of the photoreceptor-lineage cells rationally explain the generation of the cup-shaped structure of the pigmented ocellus. The present findings demonstrate the complete cell lineage of the ocellus photoreceptor cells and provide a framework elucidating the molecular and cellular mechanisms of photoreceptor development in Ciona.

Regulatory cocktail for dopaminergic neurons in a protovertebrate identified by whole-embryo single-cell transcriptomics.

The CNS of the protovertebrate Ciona intestinalis contains a single cluster of dopaminergic (DA) neurons, the coronet cells, which have been likened to the hypothalamus of vertebrates. Whole-embryo single-cell RNA sequencing (RNA-seq) assays identified Ptf1a as the most strongly expressed cell-specific transcription factor (TF) in DA/coronet cells. Knockdown of Ptf1a activity results in their loss, while misexpression results in the appearance of supernumerary DA/coronet cells. Photoreceptor cells and ependymal cells are the most susceptible to transformation, and both cell types express high levels of Meis Coexpression of both Ptf1a and Meis caused the wholesale transformation of the entire CNS into DA/coronet cells. We therefore suggest that the reiterative use of functional manipulations and single-cell RNA-seq assays is an effective means for the identification of regulatory cocktails underlying the specification of specific cell identities.

The Use of cis-Regulatory DNAs as Molecular Tools.

Ascidians possess relatively small and compact genomes. This feature enables us to easily isolate cis-regulatory DNAs of genes of interest. Particularly, cis-regulatory DNAs of genes showing tissue- or cell-type-specific expression are routinely used for the artificial induction of gene expression. This strategy helps us to label cells, tissues, and organs of interest, and to investigate gene functions through overexpression, ectopic expression, and the disruption of functions by dominant-negative forms. Thus, cis-regulatory DNAs provide a powerful tool for tissue-specific genetic manipulation in studies of ascidian development and physiology. This chapter summarizes the types of cis-regulatory DNAs as a genetic manipulation tool, describes the methods used for isolating cis-regulatory DNAs, and provide reported examples of the use of cis-regulatory DNAs as molecular tools for investigating gene functions.

ANISEED 2017: extending the integrated ascidian database to the exploration and evolutionary comparison of genome-scale datasets.

ANISEED (www.aniseed.cnrs.fr) is the main model organism database for tunicates, the sister-group of vertebrates. This release gives access to annotated genomes, gene expression patterns, and anatomical descriptions for nine ascidian species. It provides increased integration with external molecular and taxonomy databases, better support for epigenomics datasets, in particular RNA-seq, ChIP-seq and SELEX-seq, and features novel interactive interfaces for existing and novel datatypes. In particular, the cross-species navigation and comparison is enhanced through a novel taxonomy section describing each represented species and through the implementation of interactive phylogenetic gene trees for 60% of tunicate genes. The gene expression section displays the results of RNA-seq experiments for the three major model species of solitary ascidians. Gene expression is controlled by the binding of transcription factors to cis-regulatory sequences. A high-resolution description of the DNA-binding specificity for 131 Ciona robusta (formerly C. intestinalis type A) transcription factors by SELEX-seq is provided and used to map candidate binding sites across the Ciona robusta and Phallusia mammillata genomes. Finally, use of a WashU Epigenome browser enhances genome navigation, while a Genomicus server was set up to explore microsynteny relationships within tunicates and with vertebrates, Amphioxus, echinoderms and hemichordates.

Constrained vertebrate evolution by pleiotropic genes.

Despite morphological diversification of chordates over 550 million years of evolution, their shared basic anatomical pattern (or 'bodyplan') remains conserved by unknown mechanisms. The developmental hourglass model attributes this to phylum-wide conserved, constrained organogenesis stages that pattern the bodyplan (the phylotype hypothesis); however, there has been no quantitative testing of this idea with a phylum-wide comparison of species. Here, based on data from early-to-late embryonic transcriptomes collected from eight chordates, we suggest that the phylotype hypothesis would be better applied to vertebrates than chordates. Furthermore, we found that vertebrates' conserved mid-embryonic developmental programmes are intensively recruited to other developmental processes, and the degree of the recruitment positively correlates with their evolutionary conservation and essentiality for normal development. Thus, we propose that the intensively recruited genetic system during vertebrates' organogenesis period imposed constraints on its diversification through pleiotropic constraints, which ultimately led to the common anatomical pattern observed in vertebrates.

Evolutionary steps involving counterion displacement in a tunicate opsin.

Ci-opsin1 is a visible light-sensitive opsin present in the larval ocellus of an ascidian, Ciona intestinalis This invertebrate opsin belongs to the vertebrate visual and nonvisual opsin groups in the opsin phylogenetic tree. Ci-opsin1 contains candidate counterions (glutamic acid residues) at positions 113 and 181; the former is a newly acquired position in the vertebrate visual opsin lineage, whereas the latter is an ancestral position widely conserved among invertebrate opsins. Here, we show that Glu113 and Glu181 in Ci-opsin1 act synergistically as counterions, which imparts molecular properties to Ci-opsin1 intermediate between those of vertebrate- and invertebrate-type opsins. Synergy between the counterions in Ci-opsin1 was demonstrated by E113Q and E181Q mutants that exhibit a pH-dependent spectral shift, whereas only the E113Q mutation in vertebrate rhodopsin yields this spectral shift. On absorbing light, Ci-opsin1 forms an equilibrium between two intermediates with protonated and deprotonated Schiff bases, namely the MI-like and MII-like intermediates, respectively. Adding G protein caused the equilibrium to shift toward the MI-like intermediate, indicating that Ci-opsin1 has a protonated Schiff base in its active state, like invertebrate-type opsins. Ci-opsin1's G protein activation efficiency is between the efficiencies of vertebrate- and invertebrate-type opsins. Interestingly, the E113Y and E181S mutations change the molecular properties of Ci-opsin1 into those resembling invertebrate-type or bistable opsins and vertebrate ancient/vertebrate ancient-long or monostable opsins, respectively. These results strongly suggest that acquisition of counterion Glu113 changed the molecular properties of visual opsin in a vertebrate/tunicate common ancestor as a crucial step in the evolution of vertebrate visual opsins.

Revised lineage of larval photoreceptor cells in Ciona reveals archetypal collaboration between neural tube and neural crest in sensory organ formation.

The Ciona intestinalis larva has two distinct photoreceptor organs, a conventional pigmented ocellus and a nonpigmented ocellus, that are asymmetrically situated in the brain. The ciliary photoreceptor cells of these ocelli resemble visual cells of the vertebrate retina. Precise elucidation of the lineage of the photoreceptor cells will be key to understanding the developmental mechanisms of these cells as well as the evolutionary relationships between the photoreceptor organs of ascidians and vertebrates. Photoreceptor cells of the pigmented ocellus have been thought to develop from anterior animal (a-lineage) blastomeres, whereas the developmental origin of the nonpigmented ocellus has not been determined. Here, we show that the photoreceptor cells of both ocelli develop from the right anterior vegetal hemisphere: those of the pigmented ocellus from the right A9.14 cell and those of the nonpigmented ocellus from the right A9.16 cell. The pigmented ocellus is formed by a combination of two lineages of cells with distinct embryonic origins: the photoreceptor cells originate from a medial portion of the A-lineage neural plate, while the pigment cell originates from the lateral edge of the a-lineage neural plate. In light of the recently proposed close evolutionary relationship between the ocellus pigment cell of ascidians and the cephalic neural crest of vertebrates, the ascidian ocellus may represent a prototypic contribution of the neural crest to a cranial sensory organ.

Genome-wide identification and characterization of transcription start sites and promoters in the tunicate Ciona intestinalis.

The tunicate Ciona intestinalis, an invertebrate chordate, has recently emerged as a powerful model organism for gene regulation analysis. However, few studies have been conducted to identify and characterize its transcription start sites (TSSs) and promoters at the genome-wide level. Here, using TSS-seq, we identified TSSs at the genome-wide scale and characterized promoters in C. intestinalis. Specifically, we identified TSS clusters (TSCs), high-density regions of TSS-seq tags, each of which appears to originate from an identical promoter. TSCs were found not only at known TSSs but also in other regions, suggesting the existence of many unknown transcription units in the genome. We also identified candidate promoters of 79 ribosomal protein (RP) genes, each of which had the major TSS in a polypyrimidine tract and showed a sharp TSS distribution like human RP gene promoters. Ciona RP gene promoters, however, did not appear to have typical TATA boxes, unlike human RP gene promoters. In Ciona non-RP promoters, two pyrimidine-purine dinucleotides, CA and TA, were frequently used as TSSs. Despite the absence of CpG islands, Ciona TATA-less promoters showed low expression specificity like CpG-associated human TATA-less promoters. By using TSS-seq, we also predicted trans-spliced gene TSSs and found that their downstream regions had higher G+T content than those of non-trans-spliced gene TSSs. Furthermore, we identified many putative alternative promoters, some of which were regulated in a tissue-specific manner. Our results provide valuable information about TSSs and promoter characteristics in C. intestinalis and will be helpful in future analysis of transcriptional regulation in chordates.

The pre-vertebrate origins of neurogenic placodes.

The sudden appearance of the neural crest and neurogenic placodes in early branching vertebrates has puzzled biologists for over a century. These embryonic tissues contribute to the development of the cranium and associated sensory organs, which were crucial for the evolution of the vertebrate "new head". A previous study suggests that rudimentary neural crest cells existed in ancestral chordates. However, the evolutionary origins of neurogenic placodes have remained obscure owing to a paucity of embryonic data from tunicates, the closest living relatives to those early vertebrates. Here we show that the tunicate Ciona intestinalis exhibits a proto-placodal ectoderm (PPE) that requires inhibition of bone morphogenetic protein (BMP) and expresses the key regulatory determinant Six1/2 and its co-factor Eya, a developmental process conserved across vertebrates. The Ciona PPE is shown to produce ciliated neurons that express genes for gonadotropin-releasing hormone (GnRH), a G-protein-coupled receptor for relaxin-3 (RXFP3) and a functional cyclic nucleotide-gated channel (CNGA), which suggests dual chemosensory and neurosecretory activities. These observations provide evidence that Ciona has a neurogenic proto-placode, which forms neurons that appear to be related to those derived from the olfactory placode and hypothalamic neurons of vertebrates. We discuss the possibility that the PPE-derived GnRH neurons of Ciona resemble an ancestral cell type, a progenitor to the complex neuronal circuit that integrates sensory information and neuroendocrine functions in vertebrates.

Nonreproductive role of gonadotropin-releasing hormone in the control of ascidian metamorphosis.

BACKGROUND: Gonadotropin-releasing hormones (GnRHs) are neuropeptides that play central roles in the reproduction of vertebrates. In the ascidian Ciona intestinalis, GnRHs and their receptors are expressed in the nervous systems at the larval stage, when animals are not yet capable of reproduction, suggesting that the hormones have non-reproductive roles. RESULTS: We showed that GnRHs in Ciona are involved in the animal's metamorphosis by regulating tail absorption and adult organ growth. Absorption of the larval tail and growth of the adult organs are two major events in the metamorphosis of ascidians. When larvae were treated with GnRHs, they completed tail absorption more frequently than control larvae. cAMP was suggested to be a second messenger for the induction of tail absorption by GnRHs. tGnRH-3 and tGnRH-5 (the "t" indicates "tunicate") inhibited the growth of adult organs by arresting cell cycle progression in parallel with the promotion of tail absorption. CONCLUSIONS: This study provides new insights into the molecular mechanisms of ascidian metamorphosis conducted by non-reproductive GnRHs.

Transcriptional co-regulation of evolutionarily conserved microRNA/cone opsin gene pairs: implications for photoreceptor subtype specification.

The vertebrate retina contains two types of photoreceptor cells, rods and cones, which use distinct types of opsins and phototransduction proteins. Cones can be further divided into several subtypes with differing wavelength sensitivity and morphology. Although photoreceptor development has been extensively studied in a variety of vertebrate species, the mechanism by which photoreceptor subtypes are established is still largely unknown. Here we report two microRNAs (miRNAs), miR-726 and miR-729, which are potentially involved in photoreceptor subtype specification. In the medaka Oryzias latipes, the genes encoding miR-726 and miR-729 are located upstream of the red-sensitive opsin gene LWS-A and the UV-sensitive opsin gene SWS1, respectively, and are transcribed in the opposite direction from the respective opsin genes. The miR-726/LWS pair is conserved between teleosts and tetrapods, and the miR-729/SWS1 pair is conserved among teleosts. in situ hybridization analyses and fluorescence reporter assays suggest that these miRNAs are co-expressed with the respective opsins in specific cone subtypes. Potential targets of miR-726 and miR-729 predicted in silico include several transcription factors that regulate photoreceptor development. Functional analyses of cis-regulatory sequences in vivo suggest that transcription of the paired microRNA and opsin genes is co-regulated by common cis-regulatory modules. We propose an evolutionarily conserved mechanism that controls photoreceptor subtype identity through coupling between transcriptional and post-transcriptional regulations.

Characterization of the compact bicistronic microRNA precursor, miR-1/miR-133, expressed specifically in Ciona muscle tissues.

Muscle-specific miR-1/206 and miR-133 families have been suggested to play fundamental roles in skeletal and cardiac myogenesis in vertebrates. To gain insights into the relationships between the divergence of these miRs and muscular tissue types, we investigated the expression patterns of miR-1 and miR-133 in two ascidian Ciona species and compared their genomic structures with those of other chordates. We found that Ciona intestinalis and Ciona savignyi each possess a single copy of the miR-1/miR-133 cluster, which is only 350 nucleotide long. During embryogenesis, Ciona miR-1 and miR-133 are generated as a single continuous primary transcript accumulated in the nuclei of the tail muscle cells, starting at the gastrula stage. In adults, mature miR-133 and miR-1 are differentially expressed in the heart and body wall muscle. Expression of the reporter gene linked to the 850-bp upstream region of the predicted transcription start site confirmed that this region drives the muscle-specific expression of the primary transcript of miR-1/miR-133. In many deuterostome lineages, including that of Ciona, the miR-1/133 cluster is located in the same intron of the mind bomb (mib) gene in reverse orientation. Our results suggest that the origin of genomic organization and muscle-specific regulation of miR-1/133 can be traced back to the ancestor of chordates. Duplication of this miR cluster might have led to the remarkable elaboration in the morphology and function of skeletal muscles in the vertebrate lineage.

A conserved non-reproductive GnRH system in chordates.

Gonadotropin-releasing hormone (GnRH) is a neuroendocrine peptide that plays a central role in the vertebrate hypothalamo-pituitary axis. The roles of GnRH in the control of vertebrate reproductive functions have been established, while its non-reproductive function has been suggested but less well understood. Here we show that the tunicate Ciona intestinalis has in its non-reproductive larval stage a prominent GnRH system spanning the entire length of the nervous system. Tunicate GnRH receptors are phylogenetically closest to vertebrate GnRH receptors, yet functional analysis of the receptors revealed that these simple chordates have evolved a unique GnRH system with multiple ligands and receptor heterodimerization enabling complex regulation. One of the gnrh genes is conspicuously expressed in the motor ganglion and nerve cord, which are homologous structures to the hindbrain and spinal cord of vertebrates. Correspondingly, GnRH receptor genes were found to be expressed in the tail muscle and notochord of embryos, both of which are phylotypic axial structures along the nerve cord. Our findings suggest a novel non-reproductive role of GnRH in tunicates. Furthermore, we present evidence that GnRH-producing cells are present in the hindbrain and spinal cord of the medaka, Oryzias latipes, thereby suggesting the deep evolutionary origin of a non-reproductive GnRH system in chordates.