Optimizing CRISPR/Cas9 approaches in the polymorphic tunicate Ciona intestinalis

CRISPR/Cas9 became a powerful tool for genetic engineering and in vivo knockout also in the invertebrate chordate Ciona intestinalis. Ciona (ascidians, tunicates) is an important model organism because it shares developmental features with the vertebrates, considered the sister group of tunicates, and offers outstanding experimental advantages: a compact genome and an invariant developmental cell lineage that, combined with electroporation mediated transgenesis allows for precise and cell type specific targeting in vivo. A high polymorphism and the mosaic expression of electroporated constructs, however, often hamper the efficient CRISPR knockout, and an optimization in Ciona is desirable. Furthermore, seasonality and artificial maintenance settings can profit from in vitro approaches that would save on animals. Here we present improvements for the CRISPR/Cas9 protocol in silico, in vitro and in vivo. Firstly, in designing sgRNAs, prior sequencing of target genomic regions from experimental animals and alignment with reference genomes of C. robusta and C. intestinalis render a correction possible of subspecies polymorphisms. Ideally, the screening for efficient and non-polymorphic sgRNAs will generate a database compatible for worldwide Ciona populations. Secondly, we challenged in vitro assays for sgRNA validation towards reduced in vivo experimentation and report their suitability but also overefficiency concerning mismatch tolerance. Thirdly, when comparing Cas9 with Cas9:Geminin, thought to synchronize editing and homology-direct repair, we could indeed increase the in vivo efficiency and notably the access to an early expressed gene. Finally, for in vivo CRISPR, genotyping by next generation sequencing (NGS) ex vivo streamlined the definition of efficient single guides. Double CRISPR then generates large deletions and reliable phenotypic excision effects. Overall, while these improvements render CRISPR more efficient in Ciona, they are useful when newly establishing the technique and very transferable to CRISPR in other organisms.

Cellular remodeling and JAK inhibition promote zygotic gene expression in the Ciona germline

Transcription control is a major determinant of cell fate decisions in somatic tissues. By contrast, early germline fate specification in numerous vertebrate and invertebrate species relies extensively on RNA-level regulation, exerted on asymmetrically inherited maternal supplies, with little-to-no zygotic transcription. However delayed, a maternal-to-zygotic transition is nevertheless poised to complete the deployment of pre-gametic programs in the germline. Here, we focus on early germline specification in the tunicate Ciona to study zygotic genome activation. We first demonstrate that a peculiar cellular remodeling event excludes localized postplasmic Pem-1 mRNA, which encodes the general inhibitor of transcription. Subsequently, zygotic transcription begins in Pem-1-negative primordial germ cells (PGCs), as revealed by histochemical detection of elongating RNA Polymerase II, and nascent Mef2 transcripts. In addition, we uncover a provisional antagonism between JAK and MEK/BMPRI/GSK3 signaling, which controls the onset of zygotic gene expression, following cellular remodeling of PGCs. We propose a 2-step model for the onset of zygotic transcription in the Ciona germline and discuss the significance of germ plasm dislocation and remodeling in the context of developmental fate specification.

Cis-regulatory interfaces reveal the molecular mechanisms underlying the notochord gene regulatory network of Ciona

Tissue-specific gene expression is fundamental in development and evolution, and is mediated by transcription factors and by the cis-regulatory regions (enhancers) that they control. Transcription factors and their respective tissue-specific enhancers are essential components of gene regulatory networks responsible for the development of tissues and organs. Although numerous transcription factors have been characterized from different organisms, the knowledge of the enhancers responsible for their tissue-specific expression remains fragmentary. Here we use Ciona to study the enhancers associated with ten transcription factors expressed in the notochord, an evolutionary hallmark of the chordate phylum. Our results illustrate how two evolutionarily conserved transcription factors, Brachyury and Foxa2, coordinate the deployment of other notochord transcription factors. The results of these detailed cis-regulatory analyses delineate a high-resolution view of the essential notochord gene regulatory network of Ciona, and provide a reference for studies of transcription factors, enhancers, and their roles in development, disease, and evolution.

Characterization of a putative orexin receptor in Ciona intestinalis sheds light on the evolution of the orexin/hypocretin system in chordates

Tunicates are evolutionary model organisms bridging the gap between vertebrates and invertebrates. A genomic sequence in Ciona intestinalis (CiOX) shows high similarity to vertebrate orexin receptors and protostome allatotropin receptors (ATR). Here, molecular phylogeny suggested that CiOX is divergent from ATRs and human orexin receptors (hOX1/2). However, CiOX appears closer to hOX1/2 than to ATR both in terms of sequence percent identity and in its modelled binding cavity, as suggested by molecular modelling. CiOX was heterologously expressed in a recombinant HEK293 cell system. Human orexins weakly but concentration-dependently activated its Gq signalling (Ca2+ elevation), and the responses were inhibited by the non-selective orexin receptor antagonists TCS 1102 and almorexant, but only weakly by the OX1-selective antagonist SB-334867. Furthermore, the 5-/6-carboxytetramethylrhodamine (TAMRA)-labelled human orexin-A was able to bind to CiOX. Database mining was used to predict a potential endogenous C. intestinalis orexin peptide (Ci-orexin-A). Ci-orexin-A was able to displace TAMRA-orexin-A, but not to induce any calcium response at the CiOX. Consequently, we suggested that the orexin signalling system is conserved in Ciona intestinalis, although the relevant peptide-receptor interaction was not fully elucidated.

Distribution of cionin, a cholecystokinin/gastrin family peptide, and its receptor in the central nervous system of Ciona intestinalis type A

The cholecystokinin (CCK)/gastrin family peptides are involved in regulation of feeding and digestion in vertebrates. In the ascidian Ciona intestinalis type A (Ciona robusta), cionin, a CCK/gastrin family peptide, has been identified. Cionin is expressed exclusively in the central nervous system (CNS). In contrast, cionin receptor expression has been detected in the CNS, digestive tract, and ovary. Although cionin has been reported to be involved in ovulation, its physiological function in the CNS remains to be investigated. To elucidate its neural function, in the present study, we analyzed the expression of cionin and cionin receptors in the CNS. Cionin was expressed mainly in neurons residing in the anterior region of the cerebral ganglion. In contrast, the gene expressin of the cionin receptor gene CioR1, was detected in the middle part of the cerebral ganglion and showed a similar expression pattern to that of VACHT, a cholinergic neuron marker gene. Moreover, CioR1 was found to be expressed in cholinergic neurons. Consequently, these results suggest that cionin interacts with cholinergic neurons as a neurotransmitter or neuromodulator via CioR1. This study provides insights into a biological role of a CCK/gastrin family peptide in the CNS of ascidians.

Extrinsic apoptosis participates to tail regression during the metamorphosis of the chordate Ciona

Apoptosis is a regulated cell death ubiquitous in animals defined by morphological features depending on caspases. Two regulation pathways are described, currently named the intrinsic and the extrinsic apoptosis. While intrinsic apoptosis is well studied and considered ancestral among metazoans, extrinsic apoptosis is poorly studied outside mammals. Here, we address extrinsic apoptosis in the urochordates Ciona, belonging to the sister group of vertebrates. During metamorphosis, Ciona larvae undergo a tail regression depending on tissue contraction, migration and apoptosis. Apoptosis begin at the tail tip and propagates towards the trunk as a polarized wave. We identified Ci-caspase 8/10 by phylogenetic analysis as homolog to vertebrate caspases 8 and 10 that are the specific initiator of extrinsic apoptosis. We detected Ci-caspase 8/10 expression in Ciona larvae, especially at the tail tip. We showed that chemical inhibition of Ci-caspase 8/10 leads to a delay of tail regression, and Ci-caspase 8/10 loss of function induced an incomplete tail regression. The specificity between apoptotic pathways and initiator caspase suggests that extrinsic apoptosis regulates cell death during the tail regression. Our study presents rare in vivo work on extrinsic apoptosis outside mammals, and contribute to the discussion on its evolutionary history in animals.

The neuroendocrine system of Ciona intestinalis Type A, a deuterostome invertebrate and the closest relative of vertebrates

Deuterostome invertebrates, including echinoderms, hemichordates, cephalochordates, and urochordates, exhibit common and species-specific morphological, developmental, physiological, and behavioral characteristics that are regulated by neuroendocrine and nervous systems. Over the past 15 years, omics, genetic, and/or physiological studies on deuterostome invertebrates have identified low-molecular-weight transmitters, neuropeptides and their cognate receptors, and have clarified their various biological functions. In particular, there has been increasing interest on the neuroendocrine and nervous systems of Ciona intestinalis Type A, which belongs to the subphylum Urochordata and occupies the critical phylogenetic position as the closest relative of vertebrates. During the developmental stage, gamma-aminobutylic acid, D-serine, and gonadotropin-releasing hormones regulate metamorphosis of Ciona. In adults, the neuropeptidergic mechanisms underlying ovarian follicle growth, oocyte maturation, and ovulation have been elucidated. This review article provides the most recent and fundamental knowledge of the neuroendocrine and nervous systems of Ciona, and their evolutionary aspects.

A Novel Hemocyte-Derived Peptide and Its Possible Roles in Immune Response of Ciona intestinalis Type A

A wide variety of bioactive peptides have been identified in the central nervous system and several peripheral tissues in the ascidian Ciona intestinalis type A (Ciona robusta). However, hemocyte endocrine peptides have yet to be explored. Here, we report a novel 14-amino-acid peptide, CiEMa, that is predominant in the granular hemocytes and unilocular refractile granulocytes of Ciona. RNA-seq and qRT-PCR revealed the high CiEma expression in the adult pharynx and stomach. Immunohistochemistry further revealed the highly concentrated CiEMa in the hemolymph of the pharynx and epithelial cells of the stomach, suggesting biological roles in the immune response. Notably, bacterial lipopolysaccharide stimulation of isolated hemocytes for 1-4 h resulted in 1.9- to 2.4-fold increased CiEMa secretion. Furthermore, CiEMa-stimulated pharynx exhibited mRNA upregulation of the growth factor (Fgf3/7/10/22), vanadium binding proteins (CiVanabin1 and CiVanabin3), and forkhead and homeobox transcription factors (Foxl2, Hox3, and Dbx) but not antimicrobial peptides (CrPap-a and CrMam-a) or immune-related genes (Tgfbtun3, Tnfa, and Il17-2). Collectively, these results suggest that CiEMa plays roles in signal transduction involving tissue development or repair in the immune response, rather than in the direct regulation of immune response genes. The present study identified a novel Ciona hemocyte peptide, CiEMa, which paves the way for research on the biological roles of hemocyte peptides in chordates.

Anionic omega currents from single countercharge mutants in the voltage-sensing domain of Ci-VSP

The S4 segment of voltage-sensing domains (VSDs) directly responds to voltage changes by reorienting within the electric field as a permion. A narrow hydrophobic "gasket" or charge transfer center at the core of most VSDs focuses the electric field into a narrow region and catalyzes the sequential and reversible translocation of S4 positive gating charge residues across the electric field while preventing the permeation of physiological ions. Mutating specific S4 gating charges can cause ionic leak currents through the VSDs. These gating pores or omega currents play important pathophysiological roles in many diseases of excitability. Here, we show that mutating D129, a key countercharge residue in the Ciona intestinalis voltage-sensing phosphatase (Ci-VSP), leads to the generation of unique anionic omega currents. Neutralizing D129 causes a dramatic positive shift of activation, facilitates the formation of a continuous water path through the VSD, and creates a positive electrostatic potential landscape inside the VSD that contributes to its unique anionic selectivity. Increasing the population or dwell time of the conducting state by a high external pH or an engineered Cd2+ bridge markedly increases the current magnitude. Our findings uncover a new role of countercharge residues in the impermeable VSD of Ci-VSP and offer insights into mechanisms of the conduction of anionic omega currents linked to countercharge residue mutations.

Lhx3/4 initiates a cardiopharyngeal-specific transcriptional program in response to widespread FGF signaling

Individual signaling pathways, such as fibroblast growth factors (FGFs), can regulate a plethora of inductive events. According to current paradigms, signal-dependent transcription factors (TFs), such as FGF/MapK-activated Ets family factors, partner with lineage-determining factors to achieve regulatory specificity. However, many aspects of this model have not been rigorously investigated. One key question relates to whether lineage-determining factors dictate lineage-specific responses to inductive signals or facilitate these responses in collaboration with other inputs. We utilize the chordate model Ciona robusta to investigate mechanisms generating lineage-specific induction. Previous studies in C. robusta have shown that cardiopharyngeal progenitor cells are specified through the combined activity of FGF-activated Ets1/2.b and an inferred ATTA-binding transcriptional cofactor. Here, we show that the homeobox TF Lhx3/4 serves as the lineage-determining TF that dictates cardiopharyngeal-specific transcription in response to pleiotropic FGF signaling. Targeted knockdown of Lhx3/4 leads to loss of cardiopharyngeal gene expression. Strikingly, ectopic expression of Lhx3/4 in a neuroectodermal lineage subject to FGF-dependent specification leads to ectopic cardiopharyngeal gene expression in this lineage. Furthermore, ectopic Lhx3/4 expression disrupts neural plate morphogenesis, generating aberrant cell behaviors associated with execution of incompatible morphogenetic programs. Based on these findings, we propose that combinatorial regulation by signal-dependent and lineage-determinant factors represents a generalizable, previously uncategorized regulatory subcircuit we term "cofactor-dependent induction." Integration of this subcircuit into theoretical models will facilitate accurate predictions regarding the impact of gene regulatory network rewiring on evolutionary diversification and disease ontogeny.

Multi-regional expression of pancreas-related digestive enzyme genes in the intestinal chamber of the ascidian Ciona intestinalis type A

Bilateria share sequential steps in their digestive systems, and digestion occurs in a pre-absorption step within a chamber-like structure. Previous studies on the ascidian Ciona intestinalis type A, an evolutionary research model of vertebrate organs, revealed that Ciona homologs of pancreas-related exocrine digestive enzymes (XDEs) are exclusively expressed in the chamber-like bulging stomach. In the development of the gastrointestinal tract, genes for the pancreas-related transcription factors, namely Ptf1a, Nr5a2, and Pdx, are expressed near the stomach. Recent organ/tissue RNA-seq studies on two Ciona species reported that transcripts of the XDE homologs exist in the intestinal regions, as well as in the stomach. In the present study, we investigated the spatial gene expression of XDE homologs in the gastrointestinal region of the C. intestinalis type A. Whole-mount in situ hybridization using adult and juvenile specimens revealed apparent expression signals of XDE homologs in a small number of gastrointestinal epithelial cells. Furthermore, two pancreas-related transcription factor genes, Nr5a2 and Pdx, exhibited multi-regional expression along the Ciona juvenile intestines. These results imply that ascidians may form multiple digestive regions corresponding to the vertebrate pancreas.

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.

Repetitive and zonal expression profiles of absorption-related genes in the gastrointestinal tract of ascidian Ciona intestinalis type A

Intestinal absorption is essential for heterotrophic bilaterians with a tubular gut. Although the fundamental features of the digestive system were shared among chordates with evolution, the gut morphologies of vertebrates diverged and adapted to different food habitats. The ascidian Ciona intestinalis type A, a genome-wide research model of basal chordates, is used to examine the functional morphology of the intestines because of its transparent juvenile body. In the present study, the characteristic gene expression patterns (GEP) of Ciona absorptive proteins, e.g., brush border membrane enzymes for terminal digestion (lactase, maltase, APA, and APN) and transporters (SGLT1, GLUT5, PEPT1, and B0AT1), were investigated in juveniles and young adults, with a special reference to the absorption of other nutrients by pinocytosis- and phagocytosis-related proteins (megalin, cubilin, amnionless, Dab2, Rab7, LAMP, cathepsins, and MRC1). Whole-mount in situ hybridization revealed that these GEP showed multi-regional and repetitive features along the Ciona gastrointestinal tract, mainly in the stomach and several regions of the intestines. In young adults, many absorption-related genes, including pinocytosis-/phagocytosis-related genes, were also expressed between the stomach and mid-intestine. In the gastrointestinal epithelium, absorption-related genes showed zonal GEP along the epithelial structure. Comparisons of GEP, including other intestinal functions, such as nutrient digestion and intestinal protection, indicated the repetitive assignment of a well-coordinated set of intestinal GEP in the Ciona gastrointestinal tract.

Cereblon influences the timing of muscle differentiation in Ciona tadpoles

Thalidomide has a dark history as a teratogen, but in recent years, its derivates have been shown to function as potent chemotherapeutic agents. These drugs bind cereblon (CRBN), the substrate receptor of an E3 ubiquitin ligase complex, and modify its degradation targets. Despite these insights, remarkably little is known about the normal function of cereblon in development. Here, we employ Ciona, a simple invertebrate chordate, to identify endogenous Crbn targets. In Ciona, Crbn is specifically expressed in developing muscles during tail elongation before they acquire contractile activity. Crbn expression is activated by Mrf, the ortholog of MYOD1, a transcription factor important for muscle differentiation. CRISPR/Cas9-mediated mutations of Crbn lead to precocious onset of muscle contractions. By contrast, overexpression of Crbn delays contractions and is associated with decreased expression of contractile protein genes such as troponin. This reduction is possibly due to reduced Mrf protein levels without altering Mrf mRNA levels. Our findings suggest that Mrf and Crbn form a negative feedback loop to control the precision of muscle differentiation during tail elongation.

Interplays between cis- and trans-Acting Factors for Alternative Splicing in Response to Environmental Changes during Biological Invasions of Ascidians

Alternative splicing (AS), a pivotal biological process contributing to phenotypic plasticity, creates a bridge linking genotypes with phenotypes. Despite its importance, the AS mechanisms underlying environmental response and adaptation have not been well studied, and more importantly, the cis- and trans-acting factors influencing AS variation remain unclear. Using the model invasive congeneric ascidians, Ciona robusta, and Ciona savignyi, we compared their AS responses to environmental changes and explored the potential determinants. Our findings unveiled swift and dynamic AS changes in response to environmental challenges, and differentially alternative spliced genes (DASGs) were functionally enriched in transmembrane transport processes. Interestingly, both the prevalence and level of AS in C. robusta were lower than those observed in C. savignyi. Furthermore, these two indices were higher under temperature stresses compared to salinity stresses in C. savignyi. All the observed patterns underscore the species-specific and environmental context-dependent AS responses to environmental challenges. The dissimilarities in genomic structure and exon/intron size distributions between these two species likely contributed to the observed AS variation. Moreover, we identified a total of 11 and 9 serine/arginine-rich splicing factors (SRSFs) with conserved domains and gene structures in the genomes of C. robusta and C. savignyi, respectively. Intriguingly, our analysis revealed that all detected SRSFs did not exhibit prevalent AS regulations. Instead, we observed AS control over a set of genes related to splicing factors and spliceosome components. Altogether, our results elucidate species-specific and environmental challenge-dependent AS response patterns in closely related invasive ascidians. The identified splicing factors and spliceosome components under AS control offer promising candidates for further investigations into AS-mediated rapid responses to environmental challenges complementary to SRSFs.

Characterization of a novel species-specific 51-amino acid peptide, PEP51, as a caspase-3/7 activator in ovarian follicles of the ascidian, Ciona intestinalis Type A

Invertebrates lack hypothalamic-pituitary-gonadal axis, and have acquired species-specific regulatory systems for ovarian follicle development. Ascidians are marine invertebrates that are the phylogenetically closest living relatives to vertebrates, and we have thus far substantiated the molecular mechanisms underlying neuropeptidergic follicle development of the cosmopolitan species, Ciona intestinalis Type A. However, no ovarian factor has so far been identified in Ciona. In the present study, we identified a novel Ciona-specific peptide, termed PEP51, in the ovary. Immunohistochemical analysis demonstrated the specific expression of PEP51 in oocyte-associated accessory cells, test cells, of post-vitellogenic (stage III) follicles. Immunoelectron microscopy revealed that PEP51 was localized in the cytosol of test cells in early stage III follicles, which lack secretory granules. These results indicate that PEP51 acts as an intracellular factor within test cells rather than as a secretory peptide. Confocal laser microscopy verified that activation of caspase-3/7, the canonical apoptosis marker, was detected in most PEP51-positive test cells of early stage III. This colocalization of PEP51 and the apoptosis marker was consistent with immunoelectron microscopy observations demonstrating that a few normal (PEP51-negative) test cells reside in the aggregates of PEP51-positive apoptotic test cells of early stage III follicles. Furthermore, transfection of the PEP51 gene into COS-7 cells and HEK293MSR cells resulted in activation of caspase-3/7, providing evidence that PEP51 induces apoptotic signaling. Collectively, these results showed the existence of species-specific ovarian peptide-driven cell metabolism in Ciona follicle development. Consistent with the phylogenetic position of Ciona as the closest sister group of vertebrates, the present study sheds new light on the molecular and functional diversity of the regulatory systems of follicle development in the Chordata.

Ciona robusta macrophage migration inhibitory factor (Mif1 and Mif2) genes are differentially regulated in the lipopolysaccharide-challenged pharynx

The effects of lipopolysaccharide (LPS) on Mif (macrophage migration inhibitory factor) gene expression in the pharynx (haemapoetic tissue) of Ciona robusta were investigated using quantitative reverse-transcription PCR (qRT-PCR) and in situ hybridisation (ISH). To verify the induction of an inflammatory response in the pharynx, a qRT-PCR analysis was performed to evaluate the change in the expression of proinflammatory marker genes such as Mbl, Ptx-like, Tnf-α and Nf-kb, which were shown to be upregulated 1 h post LPS challenge. The change in the expression of the two Mif paralogs in the pharynx was assessed before and after stimulation, and qRT-PCR and ISH results showed that, although Mif2 and Mif2 were expressed in clusters of haemocytes in pharynx vessels, only Mif1 expression increased after LPS stimulation. This indicates that the Mif genes are differently regulated and respond to different ambient inputs that need further analysis.

A digital twin reproducing gene regulatory network dynamics of early Ciona embryos indicates robust buffers in the network

How gene regulatory networks (GRNs) encode gene expression dynamics and how GRNs evolve are not well understood, although these problems have been studied extensively. We created a digital twin that accurately reproduces expression dynamics of 13 genes that initiate expression in 32-cell ascidian embryos. We first showed that gene expression patterns can be manipulated according to predictions by this digital model. Next, to simulate GRN rewiring, we changed regulatory functions that represented their regulatory mechanisms in the digital twin, and found that in 55 of 100 cases, removal of a single regulator from a conjunctive clause of Boolean functions did not theoretically alter qualitative expression patterns of these genes. In other words, we found that more than half the regulators gave theoretically redundant temporal or spatial information to target genes. We experimentally substantiated that the expression pattern of Nodal was maintained without one of these factors, Zfpm, by changing the upstream regulatory sequence of Nodal. Such robust buffers of regulatory mechanisms may provide a basis of enabling developmental system drift, or rewiring of GRNs without changing expression patterns of downstream genes, during evolution.

Differentially expressed chaperone genes reveal a stress response required for unidirectional regeneration in the basal chordate Ciona

BACKGROUND

Unidirectional regeneration in the basal chordate Ciona intestinalis involves the proliferation of adult stem cells residing in the branchial sac vasculature and the migration of progenitor cells to the site of distal injury. However, after the Ciona body is bisected, regeneration occurs in the proximal but not in the distal fragments, even if the latter include a part of the branchial sac with stem cells. A transcriptome was sequenced and assembled from the isolated branchial sacs of regenerating animals, and the information was used to provide insights into the absence of regeneration in distal body fragments.

RESULTS

We identified 1149 differentially expressed genes, which were separated into two major modules by weighted gene correlation network analysis, one consisting of mostly upregulated genes correlated with regeneration and the other consisting of only downregulated genes associated with metabolism and homeostatic processes. The hsp70, dnaJb4, and bag3 genes were among the highest upregulated genes and were predicted to interact in an HSP70 chaperone system. The upregulation of HSP70 chaperone genes was verified and their expression confirmed in BS vasculature cells previously identified as stem and progenitor cells. siRNA-mediated gene knockdown showed that hsp70 and dnaJb4, but not bag3, are required for progenitor cell targeting and distal regeneration. However, neither hsp70 nor dnaJb4 were strongly expressed in the branchial sac vasculature of distal fragments, implying the absence of a stress response. Heat shock treatment of distal body fragments activated hsp70 and dnaJb4 expression indicative of a stress response, induced cell proliferation in branchial sac vasculature cells, and promoted distal regeneration.

CONCLUSIONS

The chaperone system genes hsp70, dnaJb4, and bag3 are significantly upregulated in the branchial sac vasculature following distal injury, defining a stress response that is essential for regeneration. The stress response is absent from distal fragments, but can be induced by a heat shock, which activates cell division in the branchial sac vasculature and promotes distal regeneration. This study demonstrates the importance of a stress response for stem cell activation and regeneration in a basal chordate, which may have implications for understanding the limited regenerative activities in other animals, including vertebrates.

ELMOD3-Rab1A-Flotillin2 cascade regulates lumen formation via vesicle trafficking in Ciona notochord

Lumen development is a crucial phase in tubulogenesis, although its molecular mechanisms are largely unknown. In this study, we discovered an ELMO domain-containing 3 (ELMOD3), which belongs to ADP-ribosylation factor GTPase-activating protein family, was necessary to form the notochord lumen in Ciona larvae. We demonstrated that ELMOD3 interacted with lipid raft protein Flotillin2 and regulated its subcellular localization. The loss-of-function of Flotillin2 prevented notochord lumen formation. Furthermore, we found that ELMOD3 also interacted with Rab1A, which is the regulatory GTPase for vesicle trafficking and located at the notochord cell surface. Rab1A mutations arrested the lumen formation, phenocopying the loss-of-function of ELMOD3 and Flotillin2. Our findings further suggested that Rab1A interactions influenced Flotillin2 localization. We thus identified a unique pathway in which ELMOD3 interacted with Rab1A, which controlled the Flotillin2-mediated vesicle trafficking from cytoplasm to apical membrane, required for Ciona notochord lumen formation.

Diverse logics and grammar encode notochord enhancers

The notochord is a defining feature of all chordates. The transcription factors Zic and ETS regulate enhancer activity within the notochord. We conduct high-throughput screens of genomic elements within developing Ciona embryos to understand how Zic and ETS sites encode notochord activity. Our screen discovers an enhancer located near Lama, a gene critical for notochord development. Reversing the orientation of an ETS site within this enhancer abolishes expression, indicating that enhancer grammar is critical for notochord activity. Similarly organized clusters of Zic and ETS sites occur within mouse and human Lama1 introns. Within a Brachyury (Bra) enhancer, FoxA and Bra, in combination with Zic and ETS binding sites, are necessary and sufficient for notochord expression. This binding site logic also occurs within other Ciona and vertebrate Bra enhancers. Collectively, this study uncovers the importance of grammar within notochord enhancers and discovers signatures of enhancer logic and grammar conserved across chordates.

A novel role for ATP2B in ascidians: Ascidian-specific mutations in ATP2B contribute to sperm chemotaxis

Sperm chemotaxis, in which sperms are attracted to conspecific eggs via species-specific attractants, plays an important role in fertilization. This phenomenon has been observed in various animals and species-specific sperm attractants have been reported in some species. However, the mechanisms involved in the reception and recognition of the species-specific attractant by the sperms is poorly studied. Previously, we found that the plasma membrane-type Ca²⁺ /ATPase (PMCA) is the receptor for the sperm-activating and -attracting factor (SAAF) in the ascidian Ciona intestinalis. To determine the role of PMCA in species-specific sperm chemotaxis, we identified the amino acid sequences of PMCAs derived from six Phlebobranchia species. The testis-specific splice variant of PMCA was found to be present in all the species investigated and the ascidian-specific sequence was detected near the 3'-terminus. Moreover, dN/dS analysis revealed that the extracellular loops 1, 2, and 4 in ascidian PMCA underwent a positive selection. These findings suggest that PMCA recognizes the species-specific structure of SAAF at the extracellular loops 1, 2, and 4, and its testis-specific C-terminal region is involved in the activation and chemotaxis of ascidian sperms.

Differentiation of endostyle cells by Nkx2-1 and FoxE in the ascidian Ciona intestinalis type A: insights into shared gene regulation in glandular- and thyroid-equivalent elements of the chordate endostyle

Due to similarities in iodine concentrations and peroxidase activities, the thyroid in vertebrates is considered to originate from the endostyle of invertebrate chordates even though it is a glandular (mucus-producing) organ for aquatic suspension feeding. Among chordates with an endostyle, urochordates are useful evolutionary research models for the study of vertebrate traits. The ascidian Ciona intestinalis forms an endostyle with specific components of glandular- and thyroid-related elements, and molecular markers have been identified for these components. Since we previously examined a simple endostyle in the larvacean Oikopleura dioica, the expression of the thyroid-related transcription factor genes, Ciona Nkx2-1 and FoxE, was perturbed by TALEN-mediated gene knockout in the present study to elucidate the shared and/or divergent features of a complex ascidian endostyle. The knockout of Ciona Nkx2-1 and FoxE exerted different effects on the morphology of the developing endostyle. The knockout of Nkx2-1 eliminated the expression of both glandular and thyroidal differentiation marker genes, e.g., vWFL1, vWFL2, CiEnds1, TPO, and Duox, while that of FoxE eliminated the expression of the differentiation marker genes, TPO and CiEnds1. The supporting element-related expression of Pax2/5/8a, Pax2/5/8b, FoxQ1, and β-tubulin persisted in the hypoplastic endostyles of Nkx2-1- and FoxE-knockout juveniles. Although the gene regulation of ascidian-specific CiEnds1 remains unclear, these results provide insights into the evolution of the vertebrate thyroid as well as the urochordate endostyle.

Co-expression network analysis of environmental canalization in the ascidian Ciona

BACKGROUND

Canalization, or buffering, is defined as developmental stability in the face of genetic and/or environmental perturbations. Understanding how canalization works is important in predicting how species survive environmental change, as well as deciphering how development can be altered in the evolutionary process. However, how developmental gene expression is linked to buffering remains unclear. We addressed this by co-expression network analysis, comparing gene expression changes caused by heat stress during development at a whole-embryonic scale in reciprocal hybrid crosses of sibling species of the ascidian Ciona that are adapted to different thermal environments.

RESULTS

Since our previous work showed that developmental buffering in this group is maternally inherited, we first identified maternal developmental buffering genes (MDBGs) in which the expression level in embryos is both correlated to the level of environmental canalization and also differentially expressed depending on the species' gender roles in hybrid crosses. We found only 15 MDBGs, all of which showed high correlation coefficient values for expression with a large number of other genes, and 14 of these belonged to a single co-expression module. We then calculated correlation coefficients of expression between MDBGs and transcription factors in the central nervous system (CNS) developmental gene network that had previously been identified experimentally. We found that, compared to the correlation coefficients between MDBGs, which had an average of 0.96, the MDBGs are loosely linked to the CNS developmental genes (average correlation coefficient 0.45). Further, we investigated the correlation of each developmental to MDBGs, showing that only four out of 62 CNS developmental genes showed correlation coefficient > 0.9, comparable to the values between MDBGs, and three of these four genes were signaling molecules: BMP2/4, Wnt7, and Delta-like.

CONCLUSIONS

We show that the developmental pathway is not centrally located within the buffering network. We found that out of 62 genes in the developmental gene network, only four genes showed correlation coefficients as high as between MDBGs. We propose that loose links to MDBGs stabilize spatiotemporally dynamic development.

Ecological Review of the Ciona Species Complex

AbstractThe set of four closely related solitary ascidians Ciona spp. were once considered a single cosmopolitan species, Ciona intestinalis, but are now recognized as genetically and morphologically distinct species. The possibility of ecological differences between the species was not widely considered in studies preceding the schism of Ciona spp. Consequently, there may be an over-generalization of the ecology of Ciona spp., with potential implications for the broad range of studies targeting these species, encompassing the evolution, development, genomics, and invasion biology of Ciona spp. We completed a comprehensive review of the ecology of Ciona spp. to establish the similarities and differences between the widely distributed Ciona robusta and C. intestinalis (and what little is known of the two other species, Ciona sp. C and Ciona sp. D). When necessary, we used study locations and the species' geographic ranges to infer the species in each study in the review. As expected, ecological similarities are the norm between the two species, spanning both abiotic and biotic interactions. However, there are also important differences that have potential implications for other aspects of the biology of Ciona spp. For example, differences in temperature and salinity tolerances likely correspond with the disparities in the geographic distribution of the species. Asymmetries in topics studied in each species diminish our ability to fully compare several aspects of the ecology of Ciona spp. and are priority areas for future research. We anticipate that our clarification of common and unique aspects of each species' ecology will help to provide context for future research in many aspects of the biology of Ciona spp.

Xbp1 and Brachyury establish an evolutionarily conserved subcircuit of the notochord gene regulatory network

Gene regulatory networks coordinate the formation of organs and structures that compose the evolving body plans of different organisms. We are using a simple chordate model, the Ciona embryo, to investigate the essential gene regulatory network that orchestrates morphogenesis of the notochord, a structure necessary for the proper development of all chordate embryos. Although numerous transcription factors expressed in the notochord have been identified in different chordates, several of them remain to be positioned within a regulatory framework. Here, we focus on Xbp1, a transcription factor expressed during notochord formation in Ciona and other chordates. Through the identification of Xbp1-downstream notochord genes in Ciona, we found evidence of the early co-option of genes involved in the unfolded protein response to the notochord developmental program. We report the regulatory interplay between Xbp1 and Brachyury, and by extending these results to Xenopus, we show that Brachyury and Xbp1 form a cross-regulatory subcircuit of the notochord gene regulatory network that has been consolidated during chordate evolution.

Omics Studies for the Identification of Ascidian Peptides, Cognate Receptors, and Their Relevant Roles in Ovarian Follicular Development

Omics studies contribute to the elucidation of genomes and profiles of gene expression. In the ascidian Ciona intestinalis Type A (Ciona robusta), mass spectrometry (MS)-based peptidomic studies have detected numerous Ciona-specific (nonhomologous) neuropeptides as well as Ciona homologs of typical vertebrate neuropeptides and hypothalamic peptide hormones. Candidates for cognate G protein-coupled receptors (GPCRs) for these peptides have been found in the Ciona transcriptome by two ways. First, Ciona homologous GPCRs of vertebrate counterparts have been detected by sequence homology searches of cognate transcriptomes. Second, the transcriptome-derived GPCR candidates have been used for machine learning-based systematic prediction of interactions not only between Ciona homologous peptides and GPCRs but also between novel Ciona peptides and GPCRs. These data have ultimately led to experimental evidence for various Ciona peptide-GPCR interactions. Comparative transcriptomics between the wildtype and Ciona vasopressin (CiVP) gene-edited Ciona provide clues to the biological functions of CiVP in ovarian follicular development and whole body growth. Furthermore, the transcriptomes of follicles treated with peptides, such as Ciona tachykinin and cionin (a Ciona cholecystokinin homolog), have revealed key regulatory genes for Ciona follicle growth, maturation, and ovulation, eventually leading to the verification of essential and novel molecular mechanisms underlying these biological events. These findings indicate that omics studies, combined with artificial intelligence and single-cell technologies, pave the way for investigating in greater details the nervous, neuroendocrine, and endocrine systems of ascidians and the molecular and functional evolution and diversity of peptidergic regulatory networks throughout chordates.

Transcriptomic Analyses Reveal 2 and 4 Family Members of Cytochromes P450 (CYP) Involved in LPS Inflammatory Response in Pharynx of Ciona robusta

Cytochromes P450 (CYP) are enzymes responsible for the biotransformation of most endogenous and exogenous agents. The expression of each CYP is influenced by a unique combination of mechanisms and factors including genetic polymorphisms, induction by xenobiotics, and regulation by cytokines and hormones. In recent years, Ciona robusta, one of the closest living relatives of vertebrates, has become a model in various fields of biology, in particular for studying inflammatory response. Using an in vivo LPS exposure strategy, next-generation sequencing (NGS) and qRT-PCR combined with bioinformatics and in silico analyses, compared whole pharynx transcripts from naïve and LPS-exposed C. robusta, and we provide the first view of cytochrome genes expression and miRNA regulation in the inflammatory response induced by LPS in a hematopoietic organ. In C. robusta, cytochromes belonging to 2B,2C, 2J, 2U, 4B and 4F subfamilies were deregulated and miRNA network interactions suggest that different conserved and species-specific miRNAs are involved in post-transcriptional regulation of cytochrome genes and that there could be an interplay between specific miRNAs regulating both inflammation and cytochrome molecules in the inflammatory response in C. robusta.

Comparative analysis of transcriptomic profiles among ascidians, zebrafish, and mice: Insights from tissue-specific gene expression

Tissue/organ-specific genes (TSGs) are important not only for understanding organ development and function, but also for investigating the evolutionary lineages of organs in animals. Here, we investigate the TSGs of 9 adult tissues of an ascidian, Ciona intestinalis Type A (Ciona robusta), which lies in the important position of being the sister group of vertebrates. RNA-seq and qRT-PCR identified the Ciona TSGs in each tissue, and BLAST searches identified their homologs in zebrafish and mice. Tissue distributions of the vertebrate homologs were analyzed and clustered using public RNA-seq data for 12 zebrafish and 30 mouse tissues. Among the vertebrate homologs of the Ciona TSGs in the neural complex, 48% and 63% showed high expression in the zebrafish and mouse brain, respectively, suggesting that the central nervous system is evolutionarily conserved in chordates. In contrast, vertebrate homologs of Ciona TSGs in the ovary, pharynx, and intestine were not consistently highly expressed in the corresponding tissues of vertebrates, suggesting that these organs have evolved in Ciona-specific lineages. Intriguingly, more TSG homologs of the Ciona stomach were highly expressed in the vertebrate liver (17-29%) and intestine (22-33%) than in the mouse stomach (5%). Expression profiles for these genes suggest that the biological roles of the Ciona stomach are distinct from those of their vertebrate counterparts. Collectively, Ciona tissues were categorized into 3 groups: i) high similarity to the corresponding vertebrate tissues (neural complex and heart), ii) low similarity to the corresponding vertebrate tissues (ovary, pharynx, and intestine), and iii) low similarity to the corresponding vertebrate tissues, but high similarity to other vertebrate tissues (stomach, endostyle, and siphons). The present study provides transcriptomic catalogs of adult ascidian tissues and significant insights into the evolutionary lineages of the brain, heart, and digestive tract of chordates.

Vasopressin Promoter Transgenic and Vasopressin Gene-Edited Ascidian, Ciona intestinalis Type A (Ciona robusta): Innervation, Gene Expression Profiles, and Phenotypes

Oxytocin (OT) and vasopressin (VP) superfamily neuropeptides are distributed in not only vertebrates but also diverse invertebrates. However, no VPergic innervation of invertebrates has ever been documented. In the ascidian, Ciona intestinalis Type A (Ciona robusta), an OT/VP superfamily peptide was identified, and the Ciona vasopressin (CiVP) induces oocyte maturation and ovulation. In the present study, we characterize the innervation and phenotypes of genetically modified Ciona: CiVP promoter-Venus transgenic and CiVP mutants. CiVP promoter-Venus transgenic Ciona demonstrated that CiVP gene was highly expressed in the cerebral ganglion and several nerves. Fluorescence was also detected in the ovary of young CiVP promoter-Venus transgenic ascidians, suggesting that the CiVP gene is also expressed temporarily in the ovary of young ascidians. Furthermore, a marked decrease of post-vitellogenic (stage III) follicles was observed in the ovary of CiVP mutants, whereas pre-vitellogenic (stage I) and vitellogenic (stage II) follicles were increased in the mutant ovary, compared with that of wildtype Ciona. Gene expression profiles showed that the expression of various genes, including genes related to ovarian follicle growth, was altered in the ovary of CiVP mutants. Altogether, these results indicated that CiVP, mainly as a neuropeptide, plays pivotal roles in diverse biological functions, including growth of early-stage ovarian follicles via regulation of the expression of a wide variety of genes. This is the first report describing a VP gene promoter-transgenic and VP gene-edited invertebrate and also on its gene expression profiles and phenotypes.

Stress granules in Ciona robusta: First evidences of TIA-1-related nucleolysin and tristetraprolin gene expression under metal exposure

Stress granules are non-membranous cytoplasmic foci, composed of non-translating messenger ribonucleoproteins, translational initiation factors and other additional proteins. They represent a primary mechanism to rapidly modulate gene expression when cells are subjected to adverse environmental conditions. Very few works have been devoted to study the presence of the molecular components of stress granules in invertebrates. In this work, we characterized the transcript sequences for two important protein components of stress granules, TIA-1-related nucleolysin (TIAR) and tristetraprolin (TTP), in the solitary ascidian Ciona robusta, an invertebrate chordate, and carried out the first studies on their gene expression under stress conditions induced by metals (Cu, Zn and Cd). Data on mRNA expression levels, provided by qRT-PCR analyses, show a generalized decrease at the second day of metal-exposure for both tiar and ttp, suggesting that metal accumulation induces acute stress and the inhibition of the transcription for the two studied proteins. In-situ hybridization analyses demonstrate that TIAR and TTP antisense riboprobes recognize circulating granular amoebocytes in the hemolymph, in both blood lacunae and tunic. The results obtained in this work increase our knowledge on the evolution of anti-stress proteins in metazoans and emphasize the importance of the transcription of tiar and ttp, which represents an efficient physiological response allowing organisms to survive in the environment under stress conditions.

Evolutionary and transcriptional analyses of a pentraxin-like component family involved in the LPS inflammatory response of Ciona robusta

Pentraxins (PTXs) are a superfamily of conserved proteins which are components of the humoral arm of innate immunity. They are considered to be functional ancestors of antibodies and are classified into short and long types. In this study, we show that a pentraxin-like component (Ptx-like) with a C-terminal PTX domain, highly homologous to the short PTX of H. sapiens CRP, and a long N-terminal domain typical of long PTXs, is involved in the inflammatory response of Ciona robusta under LPS exposure in vivo. Analyses of protein domains as well as 3D modelling and phylogenetic tree supported the close relationship of Ptx-like with mammalian CRP, suggesting that C. robusta Ptx-like shares a common ancestor in the chordate lineages. qRT-PCR analysis showed that Ptx-like was transcriptionally upregulated during the inflammatory process induced by LPS inoculation and that it is involved in the initial phase as well as the secondary phase of the inflammatory response in which matrix remodelling and the achievement of homeostasis occur. In situ hybridisation assays revealed that gene transcription was upregulated in the pharynx post-LPS challenge in vivo, and that Ptx-like was expressed by clusters of haemocytes, mainly granulocytes, inside the pharynx vessels. We also found transcript-expressing granulocytes flowing in the musculature and in the lacunae of the circulatory system. These data supported that Ptx-like is a potential molecule of the acute-phase response in C. robusta immune defence systems against bacterial infection.

Brachyury controls Ciona notochord fate as part of a feed-forward network

The notochord is a defining feature of the chordates. The transcription factor Brachyury (Bra) is a key regulator of notochord fate but here we show that it is not a unitary master regulator in the model chordate Ciona Ectopic Bra expression only partially reprograms other cell types to a notochord-like transcriptional profile and a subset of notochord-enriched genes is unaffected by CRISPR Bra disruption. We identify Foxa.a and Mnx as potential co-regulators, and find that combinatorial cocktails are more effective at reprogramming other cell types than Bra alone. We reassess the network relationships between Bra, Foxa.a and other components of the notochord gene regulatory network, and find that Foxa.a expression in the notochord is regulated by vegetal FGF signaling. It is a direct activator of Bra expression and has a binding motif that is significantly enriched in the regulatory regions of notochord-enriched genes. These and other results indicate that Bra and Foxa.a act together in a regulatory network dominated by positive feed-forward interactions, with neither being a classically defined master regulator.

Ciona Brachyury proximal and distal enhancers have different FGF dose-response relationships

Many genes are regulated by two or more enhancers that drive similar expression patterns. Evolutionary theory suggests that these seemingly redundant enhancers must have functionally important differences. In the simple ascidian chordate Ciona, the transcription factor Brachyury is induced exclusively in the presumptive notochord downstream of lineage specific regulators and FGF-responsive Ets family transcription factors. Here we exploit the ability to finely titrate FGF signaling activity via the MAPK pathway using the MEK inhibitor U0126 to quantify the dependence of transcription driven by different Brachyury reporter constructs on this direct upstream regulator. We find that the more powerful promoter-adjacent proximal enhancer and a weaker distal enhancer have fundamentally different dose-response relationships to MAPK inhibition. The Distal enhancer is more sensitive to MAPK inhibition but shows a less cooperative response, whereas the Proximal enhancer is less sensitive and more cooperative. A longer construct containing both enhancers has a complex dose-response curve that supports the idea that the proximal and distal enhancers are moderately super-additive. We show that the overall expression loss from intermediate doses of U0126 is not only a function of the fraction of cells expressing these reporters, but also involves graded decreases in expression at the single-cell level. Expression of the endogenous gene shows a comparable dose-response relationship to the full length reporter, and we find that different notochord founder cells are differentially sensitive to MAPK inhibition. Together, these results indicate that although the two Brachyury enhancers have qualitatively similar expression patterns, they respond to FGF in quantitatively different ways and act together to drive high levels of Brachyury expression with a characteristic input/output relationship. This indicates that they are fundamentally not equivalent genetic elements.

When and where genes are expressed is controlled by regulatory DNA regions known as enhancers. Genes often have multiple enhancers that control expression in different cell types or embryonic regions, but there are also genes that have multiple enhancers that control near-identical expression patterns. These ‘shadow’ enhancers are common features of many animal genomes, but it is unclear to what extent they are truly identical in function. Here we studied a pair of shadow enhancers for the notochord-specific gene Brachyury in the simple model chordate Ciona that are both directly activated by the same signaling pathway. We titrated the activity of this pathway with graded doses of a pharmacological inhibitor and measured the effects in quantitative enhancer assays. We found that the two enhancers had significant differences in sensitivity and cooperativity to the same shared regulator and are thus not identical in function when assessed quantitatively. We also identified subtle differences in sensitivity to this upstream signal between different notochord precursor cells.

A single-cell analysis of the molecular lineage of chordate embryogenesis

Progressive unfolding of gene expression cascades underlies diverse embryonic lineage development. Here, we report a single-cell RNA sequencing analysis of the complete and invariant embryonic cell lineage of the tunicate Ciona savignyi from fertilization to the onset of gastrulation. We reconstructed a developmental landscape of 47 cell types over eight cell cycles in the wild-type embryo and identified eight fate transformations upon fibroblast growth factor (FGF) inhibition. For most FGF-dependent asymmetric cell divisions, the bipotent mother cell displays the gene signature of the default daughter fate. In convergent differentiation of the two notochord lineages, we identified additional gene pathways parallel to the master regulator T/Brachyury Last, we showed that the defined Ciona cell types can be matched to E6.5-E8.5 stage mouse cell types and display conserved expression of limited number of transcription factors. This study provides a high-resolution single-cell dataset to understand chordate early embryogenesis and cell lineage differentiation.

Genomic basis of environmental adaptation in the leathery sea squirt (Styela clava)

Tunicates occupy the evolutionary position at the boundary of invertebrates and vertebrates. It exhibits adaptation to broad environmental conditions and is distributed globally. Despite hundreds of years of embryogenesis studies, the genetic basis of the invasive habits of ascidians remains largely unknown. The leathery sea squirt, Styela clava, is an important invasive species. We used the chromosomal-level genome and transcriptome of S. clava to explore its genomic- and molecular-network-based mechanisms of adaptation to environments. Compared with Ciona intestinalis type A (C. robusta), the size of the S. clava genome was expanded by 2-fold, although the gene number was comparable. An increase in transposon number and variation in dominant types were identified as potential expansion mechanisms. In the S. clava genome, the number of genes encoding the heat-shock protein 70 family and members of the complement system was expanded significantly, and cold-shock protein genes were transferred horizontally into the S. clava genome from bacteria. The expanded gene families potentially play roles in the adaptation of S. clava to its environments. The loss of key genes in the galactan synthesis pathway might explain the distinct tunic structure and hardness compared with the ascidian Ciona species. We demonstrated further that the integrated thyroid hormone pathway participated in the regulation of larval metamorphosis that provides S. clava with two opportunities for adapting to their environment. Thus, our report of the chromosomal-level leathery sea squirt genome provides a comprehensive genomic basis for the understanding of environmental adaptation in tunicates.

Rimbp, a New Marker for the Nervous System of the Tunicate Ciona robusta

Establishment of presynaptic mechanisms by proteins that regulate neurotransmitter release in the presynaptic active zone is considered a fundamental step in animal evolution. Rab3 interacting molecule-binding proteins (Rimbps) are crucial components of the presynaptic active zone and key players in calcium homeostasis. Although Rimbp involvement in these dynamics has been described in distantly related models such as fly and human, the role of this family in most invertebrates remains obscure. To fill this gap, we defined the evolutionary history of Rimbp family in animals, from sponges to mammals. We report, for the first time, the expression of the two isoforms of the unique Rimbp family member in Ciona robusta in distinct domains of the larval nervous system. We identify intronic enhancers that are able to drive expression in different nervous system territories partially corresponding to Rimbp endogenous expression. The analysis of gene expression patterns and the identification of regulatory elements of Rimbp will positively impact our understanding of this family of genes in the context of Ciona embryogenesis.

The EJC component Magoh in non-vertebrate chordates

Earliest craniates possess a newly enlarged, elaborated forebrain with new cell types and neuronal networks. A key question in vertebrate evolution is when and how this cerebral expansion took place. The exon-junction complex (EJC) plays an essential role in mRNA processing of all Eukarya. Recently, it has been proposed that the EJC represses recursive RNA splicing in Deuterostomes, with implication in human brain diseases like microcephaly and depression. However, the EJC or EJC subunit contribution to brain development in non-vertebrate Deuterostomes remained unknown. Being interested in the evolution of chordate characters, we focused on the model species, Branchiostoma lanceolatum (Cephalochordata) and Ciona robusta (Tunicata), with the aim to investigate the ancestral and the derived expression state of Magoh orthologous genes. This study identifies that Magoh is part of a conserved syntenic group exclusively in vertebrates and suggests that Magoh has experienced duplication and loss events in mammals. During early development in amphioxus and ascidian, maternal contribution and zygotic expression of Magoh genes in various types of progenitor cells and tissues are consistent with the condition observed in other Bilateria. Later in development, we also show expression of Magoh in the brain of cephalochordate and ascidian larvae. Collectively, these results provide a basis to further define what functional role(s) Magoh exerted during nervous system development and evolution.

Exploring miR-9 Involvement in Ciona intestinalis Neural Development Using Peptide Nucleic Acids

The microRNAs are small RNAs that regulate gene expression at the post-transcriptional level and can be involved in the onset of neurodegenerative diseases and cancer. They are emerging as possible targets for antisense-based therapy, even though the in vivo stability of miRNA analogues is still questioned. We tested the ability of peptide nucleic acids, a novel class of nucleic acid mimics, to downregulate miR-9 in vivo in an invertebrate model organism, the ascidian Ciona intestinalis, by microinjection of antisense molecules in the eggs. It is known that miR-9 is a well-conserved microRNA in bilaterians and we found that it is expressed in epidermal sensory neurons of the tail in the larva of C. intestinalis. Larvae developed from injected eggs showed a reduced differentiation of tail neurons, confirming the possibility to use peptide nucleic acid PNA to downregulate miRNA in a whole organism. By identifying putative targets of miR-9, we discuss the role of this miRNA in the development of the peripheral nervous system of ascidians.

Three multi-allelic gene pairs are responsible for self-sterility in the ascidian Ciona intestinalis

Many hermaphroditic organisms possess a self-incompatibility system to avoid inbreeding. Although the mechanisms of self-incompatibility in flowering plants are well known, little is known about the mechanisms of self-sterility in hermaphroditic marine invertebrates. Ascidians are hermaphroditic sessile marine invertebrates that release sperm and eggs into the surrounding seawater. Several species, including Ciona intestinalis type A (Ciona robusta), exhibit strict self-sterility. In a previous study, we found that the candidate genes responsible for self-sterility in Ciona reside in chromosome 2q (locus A) and chromosome 7q (locus B). Two pairs of multi-allelic genes, named s(sperm)-Themis-A and v(vitelline-coat)-Themis-A in locus A and s-Themis-B and v-Themis-B in locus B, are responsible for self-sterility. In this study, we identified a third multi-allelic gene pair, s-Themis-B2 and v-Themis-B2, within locus B that is also involved in this system. Genetic analysis revealed that the haplotypes of s/v-Themis-A, s/v-Themis-B and s/v-Themis-B2 play essential roles in self-sterility. When three haplotypes were matched between s-Themis and v-Themis, fertilization never occurred even in nonself crossing. Interestingly, gene targeting of either s/v-Themis-B/B2 or s/v-Themis-A by genome editing enabled self-fertilization. These results indicate that s/v-Themis-A, -B and -B2 are S-determinant genes responsible for self-sterility in the ascidian C. intestinalis type A.

Effector gene expression underlying neuron subtype-specific traits in the Motor Ganglion of Ciona

The central nervous system of the Ciona larva contains only 177 neurons. The precise regulation of neuron subtype-specific morphogenesis and differentiation observed during the formation of this minimal connectome offers a unique opportunity to dissect gene regulatory networks underlying chordate neurodevelopment. Here we compare the transcriptomes of two very distinct neuron types in the hindbrain/spinal cord homolog of Ciona, the Motor Ganglion (MG): the Descending decussating neuron (ddN, proposed homolog of Mauthner Cells in vertebrates) and the MG Interneuron 2 (MGIN2). Both types are invariantly represented by a single bilaterally symmetric left/right pair of cells in every larva. Supernumerary ddNs and MGIN2s were generated in synchronized embryos and isolated by fluorescence-activated cell sorting for transcriptome profiling. Differential gene expression analysis revealed ddN- and MGIN2-specific enrichment of a wide range of genes, including many encoding potential "effectors" of subtype-specific morphological and functional traits. More specifically, we identified the upregulation of centrosome-associated, microtubule-stabilizing/bundling proteins and extracellular guidance cues part of a single intrinsic regulatory program that might underlie the unique polarization of the ddNs, the only descending MG neurons that cross the midline. Consistent with our predictions, CRISPR/Cas9-mediated, tissue-specific elimination of two such candidate effectors, Efcab6-related and Netrin1, impaired ddN polarized axon outgrowth across the midline.

Alternative NADH dehydrogenase extends lifespan and increases resistance to xenobiotics in Drosophila

Mitochondrial alternative NADH dehydrogenase (aNDH) was found to extend lifespan when expressed in the fruit fly. We have found that fruit flies expressing aNDH from Ciona intestinalis (NDX) had 17–71% lifespan prolongation on media with different protein-tocarbohydrate ratios except NDX-expressing males that had 19% shorter lifespan than controls on a high protein diet. NDX-expressing flies were more resistant to organic xenobiotics, 2,4-dichlorophenoxyacetic acid and alloxan, and inorganic toxicant potassium iodate, and partially to sodium molybdate treatments. On the other hand, NDX-expressing flies were more sensitive to catechol and sodium chromate. Enzymatic analysis showed that NDX-expressing males had higher glucose 6-phosphate dehydrogenase activity, whilst both sexes showed increased glutathione S-transferase activity.

miR-7 Knockdown by Peptide Nucleic Acids in the Ascidian Ciona intestinalis

Peptide Nucleic Acids (PNAs) are synthetic mimics of natural oligonucleotides, which bind complementary DNA/RNA strands with high sequence specificity. They display numerous advantages, but in vivo applications are still rare. One of the main drawbacks of PNAs application is the poor cellular uptake that could be overcome by using experimental models, in which microinjection techniques allow direct delivery of molecules into eggs. Thus, in this communication, we investigated PNAs efficiency in miR-7 downregulation and compared its effects with those obtained with the commercially available antisense molecule, Antagomir (Dharmacon) in the ascidian Ciona intestinalis. Ascidians are marine invertebrates closely related to vertebrates, in which PNA techniques have not been applied yet. Our results suggested that anti-miR-7 PNAs were able to reach their specific targets in the developing ascidian embryos with high efficiency, as the same effects were obtained with both PNA and Antagomir. To the best of our knowledge, this is the first evidence that unmodified PNAs can be applied in in vivo knockdown strategies when directly injected into eggs.

Initial characterization of Wnt-Tcf functions during Ciona heart development

In vertebrate embryos, the cardiopharyngeal mesoderm gives rise to both cardiac and branchiomeric head muscles. The canonical Wnt signaling pathway regulates many aspects of cardiomyocyte specification, and modulates a balance between skeletal and cardiac myogenesis during vertebrate head muscle development. However, the role of Wnt signaling during ascidian cardiopharyngeal development remains elusive. Here, we documented the expression of Wnt pathway components during cardiopharyngeal development in Ciona, and generated tools to investigate potential roles for Wnt signaling, and its transcriptional effector Tcf, on heart vs. pharyngeal muscle fate specification. Neither focused functional analyses nor lineage-specific transcriptome profiling uncovered a significant role for Tcf during early cardiac vs. pharyngeal muscle fate choice. By contrast, Wnt gene expression patterns of Frizzled4 and Lrp4/8 and CRISPR/Cas9-mediated Tcf knock-down suggested a later requirement for Wnt signaling during heart morphogenesis and/or cardiomyocyte differentiation. This study provides a provisional set of reagents to study Wnt signaling function in Ciona, and promising insights for future analyses of Wnt functions during heart organogenesis.

Single-cell transcriptome profiling of the Ciona larval brain

The tadpole-type larva of Ciona has emerged as an intriguing model system for the study of neurodevelopment. The Ciona intestinalis connectome has been recently mapped, revealing the smallest central nervous system (CNS) known in any chordate, with only 177 neurons. This minimal CNS is highly reminiscent of larger CNS of vertebrates, sharing many conserved developmental processes, anatomical compartments, neuron subtypes, and even specific neural circuits. Thus, the Ciona tadpole offers a unique opportunity to understand the development and wiring of a chordate CNS at single-cell resolution. Here we report the use of single-cell RNAseq to profile the transcriptomes of single cells isolated by fluorescence-activated cell sorting (FACS) from the whole brain of Ciona robusta (formerly intestinalis Type A) larvae. We have also compared these profiles to bulk RNAseq data from specific subsets of brain cells isolated by FACS using cell type-specific reporter plasmid expression. Taken together, these datasets have begun to reveal the compartment- and cell-specific gene expression patterns that define the organization of the Ciona larval brain.

The Ciona myogenic regulatory factor functions as a typical MRF but possesses a novel N-terminus that is essential for activity

Electroporation-based assays were used to test whether the myogenic regulatory factor (MRF) of Ciona intestinalis (CiMRF) interferes with endogenous developmental programs, and to evaluate the importance of its unusual N-terminus for muscle development. We found that CiMRF suppresses both notochord and endoderm development when it is expressed in these tissues by a mechanism that may involve activation of muscle-specific microRNAs. Because these results add to a large body of evidence demonstrating the exceptionally high degree of functional conservation among MRFs, we were surprised to discover that non-ascidian MRFs were not myogenic in Ciona unless they formed part of a chimeric protein containing the CiMRF N-terminus. Equally surprising, we found that despite their widely differing primary sequences, the N-termini of MRFs of other ascidian species could form chimeric MRFs that were also myogenic in Ciona. This domain did not rescue the activity of a Brachyury protein whose transcriptional activation domain had been deleted, and so does not appear to constitute such a domain. Our results indicate that ascidians have previously unrecognized and potentially novel requirements for MRF-directed myogenesis. Moreover, they provide the first example of a domain that is essential to the core function of an important family of gene regulatory proteins, one that, to date, has been found in only a single branch of the family.

The gelatinase MMP-9like is involved in regulation of LPS inflammatory response in Ciona robusta

Matrix metalloproteinases (MMPs) are a family of endopeptidases collectively able to degrade the components of the extracellular matrix (ECM), with important roles in many biological processes, such as embryogenesis, normal tissue remodelling, angiogenesis and wound healing. New views on the function of MMPs reveal that they regulate inflammatory response and therefore might represent an early step in the evolution of the immune system. MMPs can affect the activity of cytokines involved in inflammation including TGF-β and TNF-α. MMPs are widely distributed in all kingdoms of life and have likely evolved from a single-domain protein which underwent successive rounds of duplications. In this study, we focused on the Ciona robusta (formerly known as Ciona intestinalis) MMP gelatinase homologue. Gene organization, phylogenetic analysis and 3D modeling supported the closest correlation of C. robusta gelatinase with the human MMP-9. Real-time PCR analysis and zymographic assay showed a prompt expression induced by LPS inoculation and an upregulation of enzymatic activity. Furthermore, we showed that before of the well-known increase of TGF-β and TNF-α levels, a MMP-9like boost occurred, suggesting a possible involvement of MMP-9like in regulating inflammatory response in C. robusta.