A Nearly Complete Genome of Ciona intestinalis Type A (C. robusta) Reveals the Contribution of Inversion to Chromosomal Evolution in the Genus Ciona

Early transcriptional similarities between two distinct neural lineages during ascidian embryogenesis

In chordates, the central nervous system arises from precursors that have distinct developmental and transcriptional trajectories. Anterior nervous systems are ontogenically associated with ectodermal lineages while posterior nervous systems are associated with mesoderm. Taking advantage of the well-documented cell lineage of ascidian embryos, we asked to what extent the transcriptional states of the different neural lineages become similar during the course of progressive lineage restriction. We performed single-cell RNA sequencing (scRNA-seq) analyses on hand-dissected neural precursor cells of the two distinct lineages, together with those of their sister cell lineages, with a high temporal resolution covering five successive cell cycles from the 16-cell to neural plate stages. A transcription factor binding site enrichment analysis of neural specific genes at the neural plate stage revealed limited evidence for shared transcriptional control between the two neural lineages, consistent with their different ontogenies. Nevertheless, PCA analysis and hierarchical clustering showed that, by neural plate stages, the two neural lineages cluster together. Consistent with this, we identified a set of genes enriched in both neural lineages at the neural plate stage, including miR-124, Celf3.a, Zic.r-b, and Ets1/2. Altogether, the current study has revealed genome-wide transcriptional dynamics of neural progenitor cells of two distinct developmental origins. Our scRNA-seq dataset is unique and provides a valuable resource for future analyses, enabling a precise temporal resolution of cell types not previously described from dissociated embryos.

The alternative enzymes-bearing tunicates lack multiple widely distributed genes coding for peripheral OXPHOS subunits

The respiratory chain alternative enzymes (AEs) NDX and AOX from the tunicate Ciona intestinalis (Ascidiacea) have been xenotopically expressed and characterized in human cells in culture and in the model organisms Drosophila melanogaster and mouse, with the purpose of developing bypass therapies to combat mitochondrial diseases in human patients with defective complexes I and III/IV, respectively. The fact that the genes coding for NDX and AOX have been lost from genomes of evolutionarily successful animal groups, such as vertebrates and insects, led us to investigate if the composition of the respiratory chain of Ciona and other tunicates differs significantly from that of humans and Drosophila, to accommodate the natural presence of AEs. We have failed to identify in tunicate genomes fifteen orthologous genes that code for subunits of the respiratory chain complexes; all of these putatively missing subunits are peripheral to complexes I, III and IV in mammals, and many are important for complex-complex interaction in supercomplexes (SCs), such as NDUFA11, UQCR11 and COX7A. Modeling of all respiratory chain subunit polypeptides of Ciona indicates significant structural divergence that is consistent with the lack of these fifteen clear orthologous subunits. We also provide evidence using Ciona AOX expressed in Drosophila that this AE cannot access the coenzyme Q pool reduced by complex I, but it is readily available to oxidize coenzyme Q molecules reduced by glycerophosphate oxidase, a mitochondrial inner membrane-bound dehydrogenase that is not involved in SCs. Altogether, our results suggest that Ciona AEs might have evolved in a mitochondrial inner membrane environment much different from that of mammals and insects, possibly without SCs; this correlates with the preferential functional interaction between these AEs and non-SC dehydrogenases in heterologous mammalian and insect systems. We discuss the implications of these findings for the applicability of Ciona AEs in human bypass therapies and for our understanding of the evolution of animal respiratory chain.

Conservation of cis-Regulatory Syntax Underlying Deuterostome Gastrulation

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

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.

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.

Ascidian embryonic cells with properties of neural-crest cells and neuromesodermal progenitors of vertebrates

Neural-crest cells and neuromesodermal progenitors (NMPs) are multipotent cells that are important for development of vertebrate embryos. In embryos of ascidians, which are the closest invertebrate relatives of vertebrates, several cells located at the border between the neural plate and the epidermal region have neural-crest-like properties; hence, the last common ancestor of ascidians and vertebrates may have had ancestral cells similar to neural-crest cells. However, these ascidian neural-crest-like cells do not produce cells that are commonly of mesodermal origin. Here we showed that a cell population located in the lateral region of the neural plate has properties resembling those of vertebrate neural-crest cells and NMPs. Among them, cells with Tbx6-related expression contribute to muscle near the tip of the tail region and cells with Sox1/2/3 expression give rise to the nerve cord. These observations and cross-species transcriptome comparisons indicate that these cells have properties similar to those of NMPs. Meanwhile, transcription factor genes Dlx.b, Zic-r.b and Snai, which are reminiscent of a gene circuit in vertebrate neural-crest cells, are involved in activation of Tbx6-related.b. Thus, the last common ancestor of ascidians and vertebrates may have had cells with properties of neural-crest cells and NMPs and such ancestral cells may have produced cells commonly of ectodermal and mesodermal origins.

Extreme genome scrambling in marine planktonic Oikopleura dioica cryptic species

Genome structural variations within species are rare. How selective constraints preserve gene order and chromosome structure is a central question in evolutionary biology that remains unsolved. Our sequencing of several genomes of the appendicularian tunicate Oikopleura dioica around the globe reveals extreme genome scrambling caused by thousands of chromosomal rearrangements, although showing no obvious morphological differences between these animals. The breakpoint accumulation rate is an order of magnitude higher than in ascidian tunicates, nematodes, Drosophila, or mammals. Chromosome arms and sex-specific regions appear to be the primary unit of macrosynteny conservation. At the microsyntenic level, scrambling did not preserve operon structures, suggesting an absence of selective pressure to maintain them. The uncoupling of the genome scrambling with morphological conservation in O. dioica suggests the presence of previously unnoticed cryptic species and provides a new biological system that challenges our previous vision of speciation in which similar animals always share similar genome structures.

Quantitative proteome dynamics across embryogenesis in a model chordate

The evolution of gene expression programs underlying the development of vertebrates remains poorly characterized. Here, we present a comprehensive proteome atlas of the model chordate Ciona, covering eight developmental stages and ∼7,000 translated genes, accompanied by a multi-omics analysis of co-evolution with the vertebrate Xenopus. Quantitative proteome comparisons argue against the widely held hourglass model, based solely on transcriptomic profiles, whereby peak conservation is observed during mid-developmental stages. Our analysis reveals maximal divergence at these stages, particularly gastrulation and neurulation. Together, our work provides a valuable resource for evaluating conservation and divergence of multi-omics profiles underlying the diversification of vertebrates.

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.

Temporospatial hierarchy and allele-specific expression of zygotic genome activation revealed by distant interspecific urochordate hybrids

Zygotic genome activation (ZGA) is a universal process in early embryogenesis of metazoan, when the quiescent zygotic nucleus initiates global transcription. However, the mechanisms related to massive genome activation and allele-specific expression (ASE) remain not well understood. Here, we develop hybrids from two deeply diverged (120 Mya) ascidian species to symmetrically document the dynamics of ZGA. We identify two coordinated ZGA waves represent early developmental and housekeeping gene reactivation, respectively. Single-cell RNA sequencing reveals that the major expression wave exhibits spatial heterogeneity and significantly correlates with cell fate. Moreover, allele-specific expression occurs in a species- rather than parent-related manner, demonstrating the divergence of cis-regulatory elements between the two species. These findings provide insights into ZGA in chordates.

A change in cis-regulatory logic underlying obligate versus facultative muscle multinucleation in chordates

Vertebrates and tunicates are sister groups that share a common fusogenic factor, Myomaker (Mymk), that drives myoblast fusion and muscle multinucleation. Yet they are divergent in when and where they express Mymk. In vertebrates, all developing skeletal muscles express Mymk and are obligately multinucleated. In tunicates, Mymk is only expressed in post-metamorphic multinucleated muscles, but is absent from mononucleated larval muscles. In this study, we demonstrate that cis-regulatory sequence differences in the promoter region of Mymk underlie the different spatiotemporal patterns of its transcriptional activation in tunicates and vertebrates. While in vertebrates Myogenic Regulatory Factors (MRFs) like MyoD1 alone are required and sufficient for Mymk transcription in all skeletal muscles, we show that transcription of Mymk in post-metamorphic muscles of the tunicate Ciona requires the combinatorial activity of MRF/MyoD and Early B-Cell Factor (Ebf). This macroevolutionary difference appears to be encoded in cis, likely due to the presence of a putative Ebf binding site adjacent to predicted MRF binding sites in the Ciona Mymk promoter. We further discuss how Mymk and myoblast fusion might have been regulated in the last common ancestor of tunicates and vertebrates, for which we propose two models.

Influence of calcium concentration on larval adhesion in a highly invasive fouling ascidian: From morphological changes to molecular mechanisms

Calcium ion (Ca2+) is involved in the protein-mediated larval adhesion of fouling ascidians, yet the effects of environmental Ca2+ on larval adhesion remain largely unexplored. Here, the larvae of fouling ascidian C. robusta were exposed to different concentrations of Ca2+. Exposures to low-concentration (0 mM and 5 mM) and high-concentration (20 mM and 40 mM) Ca2+ significantly decreased the adhesion rate of larvae, which was primarily attributed to the decreases in adhesive structure length and curvature. Changes in the expressions of genes encoding adhesion-, microvilli-, muscle contraction-, and collagen-related proteins provided a molecular-level explanation for adhesion rate reduction. Additionally, larvae likely prioritized their energy towards immunomodulation in response to Ca2+ stresses, ultimately leading to adhesion reduction. These findings advance our understanding of the influencing mechanisms of environmental Ca2+ on larval adhesion, which are expected to provide references for the development of precise antifouling strategies against ascidians and other fouling species.

Dynamics of Chromatin Opening across Larval Development in the Urochordate Ascidian Ciona savignyi

Ascidian larvae undergo tail elongation and notochord lumenogenesis, making them an ideal model for investigating tissue morphogenesis in embryogenesis. The cellular and mechanical mechanisms of these processes have been studied; however, the underlying molecular regulatory mechanism remains to be elucidated. In this study, assays for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA sequencing (RNA-seq) were applied to investigate potential regulators of the development of ascidian Ciona savignyi larvae. Our results revealed 351 and 138 differentially accessible region genes through comparisons of ATAC-seq data between stages 21 and 24 and between stages 24 and 25, respectively. A joint analysis of RNA-seq and ATAC-seq data revealed a correlation between chromatin accessibility and gene transcription. We further verified the tissue expression patterns of 12 different genes. Among them, Cs-matrix metalloproteinase 24 (MMP24) and Cs-krüppel-like factor 5 (KLF5) were highly expressed in notochord cells. Functional assay results demonstrated that both genes are necessary for notochord lumen formation and expansion. Finally, we performed motif enrichment analysis of the differentially accessible regions in different tailbud stages and summarized the potential roles of these motif-bearing transcription factors in larval development. Overall, our study found a correlation between gene expression and chromatin accessibility and provided a vital resource for understanding the mechanisms of the development of ascidian embryos.

Heterogeneity in maternal mRNAs within clutches of eggs in response to thermal stress during the embryonic stage

BACKGROUND

The origin of variation is of central interest in evolutionary biology. Maternal mRNAs govern early embryogenesis in many animal species, and we investigated the possibility that heterogeneity in maternal mRNA provisioning of eggs can be modulated by environmental stimuli.

RESULTS

We employed two sibling species of the ascidian Ciona, called here types A and B, that are adapted to different temperature regimes and can be hybridized. Previous study showed that hybrids using type B eggs had higher susceptibility to thermal stress than hybrids using type A eggs. We conducted transcriptome analyses of multiple single eggs from crosses using eggs of the different species to compare the effects of maternal thermal stress on heterogeneity in egg provisioning, and followed the effects across generations. We found overall decreases of heterogeneity of egg maternal mRNAs associated with maternal thermal stress. When the eggs produced by the F1 AB generation were crossed with type B sperm and the progeny ('ABB' generation) reared unstressed until maturation, the overall heterogeneity of the eggs produced was greater in a clutch from an individual with a heat-stressed mother compared to one from a non-heat-stressed mother. By examining individual genes, we found no consistent overall effect of thermal stress on heterogeneity of expression in genes involved in developmental buffering. In contrast, heterogeneity of expression in signaling molecules was directly affected by thermal stress.

CONCLUSIONS

Due to the absence of batch replicates and variation in the number of reads obtained, our conclusions are very limited. However, contrary to the predictions of bet-hedging, the results suggest that maternal thermal stress at the embryo stage is associated with reduced heterogeneity of maternal mRNA provision in the eggs subsequently produced by the stressed individual, but there is then a large increase in heterogeneity in eggs of the next generation, although itself unstressed. Despite its limitations, our study presents a proof of concept, identifying a model system, experimental approach and analytical techniques capable of providing a significant advance in understanding the impact of maternal environment on developmental heterogeneity.

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.

Ovarian Follicle Development in Ascidians

Ovarian follicle development is an essential process for continuation of sexually reproductive animals, and is controlled by a wide variety of regulatory factors such as neuropeptides and peptide hormones in the endocrine, neuroendocrine, and nervous systems. Moreover, while some molecular mechanisms underlying follicle development are conserved, others vary among species. Consequently, follicle development processes are closely related to the evolution and diversity of species. Ciona intestinalis type A (Ciona rubusta) is a cosmopolitan species of ascidians, which are the closest relative of vertebrates. However, unlike vertebrates, ascidians are not endowed with the hypothalamus-pituitary-gonadal axis involving pituitary gonadotropins and sexual steroids. Combined with the phylogenetic position of ascidians as the closest relative of vertebrates, such morphological and endocrine features suggest that ascidians possess both common and species-specific regulatory mechanisms in follicle development. To date, several neuropeptides have been shown to participate in the growth of vitellogenic follicles, oocyte maturation of postvitellogenic follicles, and ovulation of fully mature follicles in a developmental stage-specific fashion. Furthermore, recent studies have shed light on the evolutionary processes of follicle development throughout chordates. In this review, we provide an overview of the neuropeptidergic molecular mechanism in the premature follicle growth, oocyte maturation, and ovulation in Ciona, and comparative views of the follicle development processes of mammals and teleosts.

(Epi)genomic adaptation driven by fine geographical scale environmental heterogeneity after recent biological invasions

Elucidating processes and mechanisms involved in rapid local adaptation to varied environments is a poorly understood but crucial component in management of invasive species. Recent studies have proposed that genetic and epigenetic variation could both contribute to ecological adaptation, yet it remains unclear on the interplay between these two components underpinning rapid adaptation in wild animal populations. To assess their respective contributions to local adaptation, we explored epigenomic and genomic responses to environmental heterogeneity in eight recently colonized ascidian (Ciona intestinalis) populations at a relatively fine geographical scale. Based on MethylRADseq data, we detected strong patterns of local environment-driven DNA methylation divergence among populations, significant epigenetic isolation by environment (IBE), and a large number of local environment-associated epigenetic loci. Meanwhile, multiple genetic analyses based on single nucleotide polymorphisms (SNPs) showed genomic footprints of divergent selection. In addition, for five genetically similar populations, we detected significant methylation divergence and local environment-driven methylation patterns, indicating the strong effects of local environments on epigenetic variation. From a functional perspective, a majority of functional genes, Gene Ontology (GO) terms, and biological pathways were largely specific to one of these two types of variation, suggesting partial independence between epigenetic and genetic adaptation. The methylation quantitative trait loci (mQTL) analysis showed that the genetic variation explained only 18.67% of methylation variation, further confirming the autonomous relationship between these two types of variation. Altogether, we highlight the complementary interplay of genetic and epigenetic variation involved in local adaptation, which may jointly promote populations' rapid adaptive capacity and successful invasions in different environments. The findings here provide valuable insights into interactions between invaders and local environments to allow invasive species to rapidly spread, thus contributing to better prediction of invasion success and development of management strategies.

Single-nucleotide variants within heart enhancers increase binding affinity and disrupt heart development

Transcriptional enhancers direct precise gene expression patterns during development and harbor the majority of variants associated with phenotypic diversity, evolutionary adaptations, and disease. Pinpointing which enhancer variants contribute to changes in gene expression and phenotypes is a major challenge. Here, we find that suboptimal or low-affinity binding sites are necessary for precise gene expression during heart development. Single-nucleotide variants (SNVs) can optimize the affinity of ETS binding sites, causing gain-of-function (GOF) gene expression, cell migration defects, and phenotypes as severe as extra beating hearts in the marine chordate Ciona robusta. In human induced pluripotent stem cell (iPSC)-derived cardiomyocytes, a SNV within a human GATA4 enhancer increases ETS binding affinity and causes GOF enhancer activity. The prevalence of suboptimal-affinity sites within enhancers creates a vulnerability whereby affinity-optimizing SNVs can lead to GOF gene expression, changes in cellular identity, and organismal-level phenotypes that could contribute to the evolution of novel traits or diseases.

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.

Establishment of a developmental atlas and transgenetic tools in the ascidian Styela clava

The ascidian Styela clava is an ecologically important species that is distributed along coastal regions worldwide. It has a long history as a model animal for evolutionary and developmental biology research owing to its phylogenetic position between vertebrates and invertebrates, and its classical mosaic expression patterns. However, the standard developmental atlas and protocols and tools for molecular manipulation of this organism are inadequate. In this study, we established a standard developmental table and provided a web-based digital image resource for S. clava embryogenesis at each developmental stage from fertilized eggs to hatching larvae by utilizing confocal laser microscopy and 3D reconstruction images. It takes around 10 h for fertilized eggs to develop into swimming larvae and 20-30 min to complete the tail regression processes at the metamorphic stage. We observed that the notochord cells in S. clava embryos did not produce an extracellular lumen like Ciona robusta, but showed polarized elongation behaviors, providing us an ideal comparative model to study tissue morphogenesis. In addition, we established a chemical-washing procedure to remove the chorion easily from the fertilized eggs. Based on the dechorionation technique, we further realized transgenic manipulation by electroporation and successfully applied tissue-specific fluorescent labeling in S. clava embryos. Our work provides a standard imaging atlas and powerful genetic tools for investigating embryogenesis and evolution using S. clava as a model organism.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00200-2.

Gene expression in notochord and nuclei pulposi: a study of gene families across the chordate phylum

The transition from notochord to vertebral column is a crucial milestone in chordate evolution and in prenatal development of all vertebrates. As ossification of the vertebral bodies proceeds, involutions of residual notochord cells into the intervertebral discs form the nuclei pulposi, shock-absorbing structures that confer flexibility to the spine. Numerous studies have outlined the developmental and evolutionary relationship between notochord and nuclei pulposi. However, the knowledge of the similarities and differences in the genetic repertoires of these two structures remains limited, also because comparative studies of notochord and nuclei pulposi across chordates are complicated by the gene/genome duplication events that led to extant vertebrates. Here we show the results of a pilot study aimed at bridging the information on these two structures. We have followed in different vertebrates the evolutionary trajectory of notochord genes identified in the invertebrate chordate Ciona, and we have evaluated the extent of conservation of their expression in notochord cells. Our results have uncovered evolutionarily conserved markers of both notochord development and aging/degeneration of the nuclei pulposi.

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.

De novo genome assembly and comparative genomics for the colonial ascidian Botrylloides violaceus

Ascidians have the potential to reveal fundamental biological insights related to coloniality, regeneration, immune function, and the evolution of these traits. This study implements a hybrid assembly technique to produce a genome assembly and annotation for the botryllid ascidian, Botrylloides violaceus. A hybrid genome assembly was produced using Illumina, Inc. short and Oxford Nanopore Technologies long-read sequencing technologies. The resulting assembly is comprised of 831 contigs, has a total length of 121 Mbp, N50 of 1 Mbp, and a BUSCO score of 96.1%. Genome annotation identified 13 K protein-coding genes. Comparative genomic analysis with other tunicates reveals patterns of conservation and divergence within orthologous gene families even among closely related species. Characterization of the Wnt gene family, encoding signaling ligands involved in development and regeneration, reveals conserved patterns of subfamily presence and gene copy number among botryllids. This supports the use of genomic data from nonmodel organisms in the investigation of biological phenomena.

The crucial role of CTCF in mitotic progression during early development of sea urchin

CCCTC-binding factor (CTCF), an insulator protein with 11 zinc fingers, is enriched at the boundaries of topologically associated domains (TADs) in eukaryotic genomes. In this study, we isolated and analyzed the cDNAs encoding HpCTCF, the CTCF homolog in the sea urchin Hemicentrotus pulcherrimus, to investigate its expression patterns and functions during the early development of sea urchin. HpCTCF contains nine zinc fingers corresponding to fingers 2-10 of the vertebrate CTCF. Expression pattern analysis revealed that HpCTCF mRNA was detected at all developmental stages and in the entire embryo. Upon expressing the HpCTCF-GFP fusion protein in early embryos, we observed its uniform distribution within interphase nuclei. However, during mitosis, it disappeared from the chromosomes and subsequently reassembled on the chromosome during telophase. Moreover, the morpholino-mediated knockdown of HpCTCF resulted in mitotic arrest during the morula to blastula stage. Most of the arrested chromosomes were not phospholylated at serine 10 of histone H3, indicating that mitosis was arrested at the telophase by HpCTCF depletion. Furthermore, impaired sister chromatid segregation was observed using time-lapse imaging of HpCTCF-knockdown embryos. Thus, HpCTCF is essential for mitotic progression during the early development of sea urchins, especially during the telophase-to-interphase transition. However, the normal development of pluteus larvae in CRISPR-mediated HpCTCF-knockout embryos suggests that disruption of zygotic HpCTCF expression has little effect on embryonic and larval development.

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.

Expression and possible functions of a horizontally transferred glycosyl hydrolase gene, GH6-1, in Ciona embryogenesis

BACKGROUND

The Tunicata or Urochordata is the only animal group with the ability to synthesize cellulose directly and cellulose is a component of the tunic that covers the entire tunicate body. The genome of Ciona intestinalis type A contains a cellulose synthase gene, CesA, that it acquired via an ancient, horizontal gene transfer. CesA is expressed in embryonic epidermal cells and functions in cellulose production. Ciona CesA is composed of both a glycosyltransferase domain, GT2, and a glycosyl hydrolase domain, GH6, which shows a mutation at a key position and seems functionless. Interestingly, the Ciona genome contains a glycosyl hydrolase gene, GH6-1, in which the GH6 domain seems intact. This suggests expression and possible functions of GH6-1 during Ciona embryogenesis. Is GH6-1 expressed during embryogenesis? If so, in what tissues is the gene expressed? Does GH6-1 serve a function? If so, what is it? Answers to these questions may advance our understanding of evolution of this unique animal group.

RESULTS

Quantitative reverse transcription PCR and in situ hybridization revealed that GH6-1 is expressed in epidermis of tailbud embryos and in early swimming larvae, a pattern similar to that of CesA. Expression is downregulated at later stages and becomes undetectable in metamorphosed juveniles. The GH6-1 expression level is higher in the anterior-trunk region and caudal-tip regions of late embryos. Single-cell RNA sequencing analysis of the late tailbud stage showed that cells of three clusters with epidermal identity express GH6-1, and that some of them co-express CesA. TALEN-mediated genome editing was used to generate GH6-1 knockout Ciona larvae. Around half of TALEN-electroporated larvae showed abnormal development of adhesive papillae and altered distribution of surface cellulose. In addition, three-fourths of TALEN-electroporated animals failed to complete larval metamorphosis.

CONCLUSIONS

This study showed that tunicate GH6-1, a gene that originated by horizontal gene transfer of a prokaryote gene, is recruited into the ascidian genome, and that it is expressed and functions in epidermal cells of ascidian embryos. Although further research is required, this observation demonstrates that both CesA and GH6-1 are involved in tunicate cellulose metabolism, impacting tunicate morphology and ecology.

BMP signaling is required to form the anterior neural plate border in ascidian embryos

Cranial neurogenic placodes have been considered vertebrate innovations. However, anterior neural plate border (ANB) cells of ascidian embryos share many properties with vertebrate neurogenic placodes; therefore, it is now believed that the last common ancestor of vertebrates and ascidians had embryonic structures similar to neurogenic placodes of vertebrate embryos. Because BMP signaling is important for specifying the placode region in vertebrate embryos, we examined whether BMP signaling is also involved in gene expression in the ANB region of ascidian embryos. Our data indicated that Admp, a divergent BMP family member, is mainly responsible for BMP signaling in the ANB region, and that two BMP-antagonists, Noggin and Chordin, restrict the domain, in which BMP signaling is activated, to the ANB region, and prevent it from expanding to the neural plate. BMP signaling is required for expression of Foxg and Six1/2 at the late gastrula stage, and also for expression of Zf220, which encodes a zinc finger transcription factor in late neurula embryos. Because Zf220 negatively regulates Foxg, when we downregulated Zf220 by inhibiting BMP signaling, Foxg was upregulated, resulting in one large palp instead of three palps (adhesive organs derived from ANB cells). Functions of BMP signaling in specification of the ANB region give further support to the hypothesis that ascidian ANB cells share an evolutionary origin with vertebrate cranial placodes.

Three amphioxus reference genomes reveal gene and chromosome evolution of chordates

The slow-evolving invertebrate amphioxus has an irreplaceable role in advancing our understanding of the vertebrate origin and innovations. Here we resolve the nearly complete chromosomal genomes of three amphioxus species, one of which best recapitulates the 17 chordate ancestor linkage groups. We reconstruct the fusions, retention, or rearrangements between descendants of whole-genome duplications, which gave rise to the extant microchromosomes likely existed in the vertebrate ancestor. Similar to vertebrates, the amphioxus genome gradually establishes its three-dimensional chromatin architecture at the onset of zygotic activation and forms two topologically associated domains at the Hox gene cluster. We find that all three amphioxus species have ZW sex chromosomes with little sequence differentiation, and their putative sex-determining regions are nonhomologous to each other. Our results illuminate the unappreciated interspecific diversity and developmental dynamics of amphioxus genomes and provide high-quality references for understanding the mechanisms of chordate functional genome evolution.

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.

Improved Genome Editing in the Ascidian Ciona with CRISPR/Cas9 and TALEN

The ascidian Ciona intestinalis type A (or Ciona robusta) is an important organism for elucidating the mechanisms that make the chordate body plan. CRISPR/Cas9 and TAL effector nuclease (TALEN) are widely used to quickly address genetic functions in Ciona. Our previously reported method of CRISPR/Cas9-mediated mutagenesis in this animal has inferior mutation rates compared to those of TALENs. We here describe an updated way to effectively mutate genes with CRISPR/Cas9 in Ciona. Although the construction of TALENs is much more laborious than that of CRISPR/Cas9, this technique is useful for tissue-specific knockouts that are not easy even by the optimized CRISPR/Cas9 method.

Evolution of 3D chromatin organization at different scales

Most animal genomes fold in 3D chromatin domains called topologically associated domains (TADs) that facilitate interactions between cis-regulatory elements (CREs) and promoters. Owing to their critical role in the control of developmental gene expression, we explore how TADs have shaped animal evolution. In the light of recent studies that profile TADs in disparate animal lineages, we discuss their phylogenetic distribution and the mechanisms that underlie their formation. We present evidence indicating that TADs are plastic entities composed of genomic strata of different ages: ancient cores are combined with newer regions and brought into extant TADs through genomic rearrangements. We highlight that newly incorporated TAD strata enable the establishment of new CRE-promoter interactions and in turn new expression patterns that can drive phenotypical innovation. We further highlight how subtle changes in chromatin folding may fine-tune the expression levels of developmental genes and hold a potential for evolutionary significance.

An FGF timer for zygotic genome activation

Zygotic genome activation has been extensively studied in a variety of systems including flies, frogs, and mammals. However, there is comparatively little known about the precise timing of gene induction during the earliest phases of embryogenesis. Here we used high-resolution in situ detection methods, along with genetic and experimental manipulations, to study the timing of zygotic activation in the simple model chordate Ciona with minute-scale temporal precision. We found that two Prdm1 homologs in Ciona are the earliest genes that respond to FGF signaling. We present evidence for a FGF timing mechanism that is driven by ERK-mediated derepression of the ERF repressor. Depletion of ERF results in ectopic activation of FGF target genes throughout the embryo. A highlight of this timer is the sharp transition in FGF responsiveness between the eight- and 16-cell stages of development. We propose that this timer is an innovation of chordates that is also used by vertebrates.

Regulators specifying cell fate activate cell cycle regulator genes to determine cell numbers in ascidian larval tissues

In animal development, most cell types stop dividing before terminal differentiation; thus, cell cycle control is tightly linked to cell differentiation programmes. In ascidian embryos, cell lineages do not vary among individuals, and rounds of the cell cycle are determined according to cell lineages. Notochord and muscle cells stop dividing after eight or nine rounds of cell division depending on their lineages. In the present study, we showed that a Cdk inhibitor, Cdkn1.b, is responsible for stopping cell cycle progression in these lineages. Cdkn1.b is also necessary for epidermal cells to stop dividing. In contrast, mesenchymal and endodermal cells continue to divide even after hatching, and Myc is responsible for maintaining cell cycle progression in these tissues. Expression of Cdkn1.b in notochord and muscle is controlled by transcription factors that specify the developmental fate of notochord and muscle. Likewise, expression of Myc in mesenchyme and endoderm is under control of transcription factors that specify the developmental fate of mesenchyme and endoderm. Thus, cell fate specification and cell cycle control are linked by these transcription factors.

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.

Molecular Characterization and Expression Analysis of Putative Class C (Glutamate Family) G Protein-Coupled Receptors in Ascidian Styela clava

In this study, we performed the genome-wide domain analysis and sequence alignment on the genome of Styela clava, and obtained a repertoire of 204 putative GPCRs, which exhibited a highly reduced gene number compared to vertebrates and cephalochordates. In this repertoire, six Class C GPCRs, including four metabotropic glutamate receptors (Sc-GRMs), one calcium-sensing receptor (Sc-CaSR), and one gamma-aminobutyric acid (GABA) type B receptor 2-like (Sc-GABABR2-like) were identified, with the absence of type 1 taste and vomeronasal receptors. All the Sc-GRMs and Sc-CaSR contained the typical "Venus flytrap" and cysteine-rich domains required for ligand binding and subsequent propagation of conformational changes. In swimming larvae, Sc-grm3 and Sc-casr were mainly expressed at the junction of the sensory vesicle and tail nerve cord while the transcripts of Sc-grm4, Sc-grm7a, and Sc-grm7b appeared at the anterior trunk, which suggested their important functions in neurotransmission. The high expression of these Class C receptors at tail-regression and metamorphic juvenile stages hinted at their potential involvement in regulating metamorphosis. In adults, the transcripts were highly expressed in several peripheral tissues, raising the possibility that S. clava Class C GPCRs might function as neurotransmission modulators peripherally after metamorphosis. Our study systematically characterized the ancestral chordate Class C GPCRs to provide insights into the origin and evolution of these receptors in chordates and their roles in regulating physiological and morphogenetic changes relevant to the development and environmental adaption.

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.

Protein kinases and protein phosphatases encoded in the Ciona robusta genome

Protein kinases (PKs) and protein phosphatases (PPs) regulate the phosphorylation of proteins that are involved in a variety of biological processes. To study such biological processes systematically, it is important to know the whole repertoire of PKs and PPs encoded in a genome. In the present study, we surveyed the genome of an ascidian (Ciona robusta or Ciona intestinalis type A) to comprehensively identify the genes that encoded PKs and PPs. Because ascidians belong to the sister group of vertebrates, a comparison of the whole repertoire of PKs and PPs of ascidians with those of vertebrates may help to delineate the complements of these proteins that were present in the last common ancestor of these two groups of animals. Our results show that the repertory of PPs was much more expanded in vertebrates than the repertory of PKs. We also showed that approximately 75% of PKs and PPs were expressed during development from eggs to larvae. Thus, the present study provides catalogs for PKs and PPs encoded in the ascidian genome. These catalogs will be useful for systematic studies of biological processes that involve phosphorylation and for evolutionary studies of the origin of vertebrates.

Deeply conserved synteny and the evolution of metazoan chromosomes

Animal genomes show networks of deeply conserved gene linkages whose phylogenetic scope and chromosomal context remain unclear. Here, we report chromosome-scale conservation of synteny among bilaterians, cnidarians, and sponges and use comparative analysis to reconstruct ancestral chromosomes across major animal groups. Comparisons among diverse metazoans reveal the processes of chromosome evolution that produced contemporary karyotypes from their Precambrian progenitors. On the basis of these findings, we introduce a simple algebraic representation of chromosomal change and use it to establish a unified systematic framework for metazoan chromosome evolution. We find that fusion-with-mixing, a previously unappreciated mode of chromosome change, has played a central role. We find that relicts of several metazoan chromosomal units are preserved in unicellular eukaryotes. These conserved pre-metazoan linkages include the chromosomal unit that encodes the most diverse set of metazoan homeobox genes, suggesting a candidate genomic context for the early diversification of this key gene family.

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.

Developmental biology of the larvacean Oikopleura dioica: Genome resources, functional screening, and imaging

The larvacean Oikopleura dioica is a cosmopolitan planktonic chordate and is closely related to vertebrates. It is characterized by a tadpole-shaped morphology with notochord flanked by muscle in the tail and brain on the dorsal side, a short life cycle of five days, a compact genome of approximately 56 Mb, a simple and transparent body with a small number of cells (~4000 in functional juveniles), invariant embryonic cell lineages, and fast development that ensures complete morphogenesis and organ formation 10 h after fertilization. With these features, this marine chordate is a promising and advantageous animal model in which genetic manipulation is feasible. In this review, we introduce relevant resources and modern techniques that have been developed: (1) Genome and transcriptomes. Oikopleura dioica has the smallest genome among non-parasitic metazoans. Its genome databases have been generated using three geographically distant O. dioica populations, and several intra-species sequence differences are becoming evident; (2) Functional genetic knockdown techniques. Comprehensive screening of genes is feasible using ovarian microinjection and double-strand DNA-induced gene knockdown; and (3) Live imaging of embryos and larvae. Application of these techniques has uncovered novel aspects of development, including meiotic cell arrest, left-right patterning, epidermal cell patterning, and mouth formation involving the connection of ectoderm and endoderm sheets. Oikopleura dioca has become very useful for developmental and evolutionary studies in chordates.

Tbx15/18/22 shares a binding site with Tbx6-r.b to maintain expression of a muscle structural gene in ascidian late embryos

The ascidian larval tail contains muscle cells for swimming. Most of these muscle cells differentiate autonomously. The genetic program behind this autonomy has been studied extensively and the genetic cascade from maternal factors to initiation of expression of a muscle structural gene, Myl.c, has been uncovered; Myl.c expression is directed initially by transcription factor Tbx6-r.b at the 64-cell stage and then by the combined actions of Tbx6-r.b and Mrf from the gastrula to early tailbud stages. In the present study, we showed that transcription of Myl.c continued in late tailbud embryos and larvae, although a fusion protein of Tbx6-r.b and GFP was hardly detectable in late tailbud embryos. A knockdown experiment, reporter assay, and in vitro binding assay indicated that an essential cis-regulatory element of Myl.c that bound Tbx6-r.b in early embryos bound Tbx15/18/22 in late embryos to maintain expression of Myl.c. We also found that Tbx15/18/22 was controlled by Mrf, which constitutes a regulatory loop with Tbx6-r.b. Therefore, our data indicated that Tbx15/18/22 was activated initially under control of this regulatory loop as in the case of Myl.c, and then Tbx15/18/22 maintained the expression of Myl.c after Tbx6-r.b had disappeared. RNA-sequencing of Tbx15/18/22 morphant embryos revealed that many muscle structural genes were regulated similarly by Tbx15/18/22. Thus, the present study revealed the mechanisms of maintenance of transcription of muscle structural genes in late embryos in which Tbx15/18/22 takes the place of Tbx6-r.b.

The Genome of the "Sea Vomit" Didemnum vexillum

Tunicates are the sister group of vertebrates and thus occupy a key position for investigations into vertebrate innovations as well as into the consequences of the vertebrate-specific genome duplications. Nevertheless, tunicate genomes have not been studied extensively in the past, and comparative studies of tunicate genomes have remained scarce. The carpet sea squirt Didemnum vexillum, commonly known as “sea vomit”, is a colonial tunicate considered an invasive species with substantial ecological and economical risk. We report the assembly of the D. vexillum genome using a hybrid approach that combines 28.5 Gb Illumina and 12.35 Gb of PacBio data. The new hybrid scaffolded assembly has a total size of 517.55 Mb that increases contig length about eightfold compared to previous, Illumina-only assembly. As a consequence of an unusually high genetic diversity of the colonies and the moderate length of the PacBio reads, presumably caused by the unusually acidic milieu of the tunic, the assembly is highly fragmented (L50 = 25,284, N50 = 6539). It is sufficient, however, for comprehensive annotations of both protein-coding genes and non-coding RNAs. Despite its shortcomings, the draft assembly of the “sea vomit” genome provides a valuable resource for comparative tunicate genomics and for the study of the specific properties of colonial ascidians.

Cell geometry, signal dampening, and a bimodal transcriptional response underlie the spatial precision of an ERK-mediated embryonic induction

Precise control of lineage segregation is critical for the development of multicellular organisms, but our quantitative understanding of how variable signaling inputs are integrated to activate lineage-specific gene programs remains limited. Here, we show how precisely two out of eight ectoderm cells adopt neural fates in response to ephrin and FGF signals during ascidian neural induction. In each ectoderm cell, FGF signals activate ERK to a level that mirrors its cell contact surface with FGF-expressing mesendoderm cells. This gradual interpretation of FGF inputs is followed by a bimodal transcriptional response of the immediate early gene, Otx, resulting in its activation specifically in the neural precursors. At low levels of ERK, Otx is repressed by an ETS family transcriptional repressor, ERF2. Ephrin signals are critical for dampening ERK activation levels across ectoderm cells so that only neural precursors exhibit above-threshold levels, evade ERF repression, and "switch on" Otx transcription.

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.

Conservative and Atypical Ferritins of Sponges

Ferritins comprise a conservative family of proteins found in all species and play an essential role in resistance to redox stress, immune response, and cell differentiation. Sponges (Porifera) are the oldest Metazoa that show unique plasticity and regenerative potential. Here, we characterize the ferritins of two cold-water sponges using proteomics, spectral microscopy, and bioinformatic analysis. The recently duplicated conservative HdF1a/b and atypical HdF2 genes were found in the Halisarca dujardini genome. Multiple related transcripts of HpF1 were identified in the Halichondria panicea transcriptome. Expression of HdF1a/b was much higher than that of HdF2 in all annual seasons and regulated differently during the sponge dissociation/reaggregation. The presence of the MRE and HRE motifs in the HdF1 and HdF2 promotor regions and the IRE motif in mRNAs of HdF1 and HpF indicates that sponge ferritins expression depends on the cellular iron and oxygen levels. The gel electrophoresis combined with specific staining and mass spectrometry confirmed the presence of ferric ions and ferritins in multi-subunit complexes. The 3D modeling predicts the iron-binding capacity of HdF1 and HpF1 at the ferroxidase center and the absence of iron-binding in atypical HdF2. Interestingly, atypical ferritins lacking iron-binding capacity were found in genomes of many invertebrate species. Their function deserves further research.

Nkx2-1 and FoxE regionalize glandular (mucus-producing) and thyroid-equivalent traits in the endostyle of the chordate Oikopleura dioica

The endostyle is a ventral pharyngeal organ used for internal filter feeding of basal chordates and is considered homologous to the follicular thyroid of vertebrates. It contains mucus-producing (glandular) and thyroid-equivalent regions organized along the dorsoventral (DV) axis. Although thyroid-related genes (Nkx2-1, FoxE, and thyroid peroxidase (TPO)) are known to be expressed in the endostyle, their roles in establishing regionalization within the organ have not been demonstrated. We report that Nkx2-1 and FoxE are essential for establishing DV axial identity in the endostyle of Oikopleura dioica. Genome and expression analyses showed von Willebrand factor-like (vWFL) and TPO/dual oxidase (Duox)/Nkx2-1/FoxE as orthologs of glandular and thyroid-related genes, respectively. Knockdown experiments showed that Nkx2-1 is necessary for the expression of glandular and thyroid-related genes, whereas FoxE is necessary only for thyroid-related genes. Moreover, Nkx2-1 expression is necessary for FoxE expression in larvae during organogenesis. The results demonstrate the essential roles of Nkx2-1 and FoxE in establishing regionalization in the endostyle, including (1) the Nkx2-1-dependent glandular region, and (2) the Nkx2-1/FoxE-dependent thyroid-equivalent region. DV axial regionalization may be responsible for organizing glandular and thyroid-equivalent traits of the pharynx along the DV axis.

Two distinct motifs for Zic-r.a drive specific gene expression in two cell lineages

Zic-r.a, a maternal transcription factor, specifies posterior fate in ascidian embryos. However, its direct target, Tbx6-r.b, does not contain typical Zic-r.a-binding sites in its regulatory region. Using an in vitro selection assay, we found that Zic-r.a binds to sites dissimilar to the canonical motif, by which it activates Tbx6-r.b in a sub-lineage of muscle cells. These sites with non-canonical motifs have weak affinity for Zic-r.a; therefore, it activates Tbx6-r.b only in cells expressing Zic-r.a abundantly. Meanwhile, we found that Zic-r.a expressed zygotically in late embryos activates neural genes through canonical sites. Because different zinc-finger domains of Zic-r.a are important for driving reporters with canonical and non-canonical sites, it is likely that the non-canonical motif is not a divergent version of the canonical motif. In other words, our data indicate that the non-canonical motif represents a motif distinct from the canonical motif. Thus, Zic-r.a recognizes two distinct motifs to activate two sets of genes at two timepoints in development. This article has an associated 'The people behind the papers' interview.

The gene regulatory system for specifying germ layers in early embryos of the simple chordate

In animal embryos, gene regulatory networks control the dynamics of gene expression in cells and coordinate such dynamics among cells. In ascidian embryos, gene expression dynamics have been dissected at the single-cell resolution. Here, we revealed mathematical functions that represent the regulatory logics of all regulatory genes expressed at the 32-cell stage when the germ layers are largely specified. These functions collectively explain the entire mechanism by which gene expression dynamics are controlled coordinately in early embryos. We found that regulatory functions for genes expressed in each of the specific lineages contain a common core regulatory mechanism. Last, we showed that the expression of the regulatory genes became reproducible by calculation and controllable by experimental manipulations. Thus, these regulatory functions represent an architectural design for the germ layer specification of this chordate and provide a platform for simulations and experiments to understand the operating principles of gene regulatory networks.

3D genomics across the tree of life reveals condensin II as a determinant of architecture type

We investigated genome folding across the eukaryotic tree of life. We find two types of three-dimensional (3D) genome architectures at the chromosome scale. Each type appears and disappears repeatedly during eukaryotic evolution. The type of genome architecture that an organism exhibits correlates with the absence of condensin II subunits. Moreover, condensin II depletion converts the architecture of the human genome to a state resembling that seen in organisms such as fungi or mosquitoes. In this state, centromeres cluster together at nucleoli, and heterochromatin domains merge. We propose a physical model in which lengthwise compaction of chromosomes by condensin II during mitosis determines chromosome-scale genome architecture, with effects that are retained during the subsequent interphase. This mechanism likely has been conserved since the last common ancestor of all eukaryotes.