Cataloging transcription factor and major signaling molecule genes for functional genomic studies in Ciona intestinalis

Single-cell Transcriptomic Studies Unveil Potential Nodes of the Notochord Gene Regulatory Network

Transcription factors (TFs) are DNA-binding proteins able to modulate the timing, location, and levels of gene expression by binding to regulatory DNA regions. Therefore, the repertoire of TFs present in the genome of a multicellular organism and the expression of variable constellations of TFs in different cellular cohorts determine the distinctive characteristics of developing tissues and organs. The information on tissue-specific assortments of TFs, their cross-regulatory interactions, and the genes/regulatory regions targeted by each TF is summarized in gene regulatory networks (GRNs), which provide genetic blueprints for the specification, development, and differentiation of multicellular structures. In this study, we review recent transcriptomic studies focused on the complement of TFs expressed in the notochord, a distinctive feature of all chordates. We analyzed notochord-specific datasets available from organisms representative of the three chordate subphyla, and highlighted lineage-specific variations in the suite of TFs expressed in their notochord. We framed the resulting findings within a provisional evolutionary scenario, which allows the formulation of hypotheses on the genetic/genomic changes that sculpted the structure and function of the notochord on an evolutionary scale.

Actin Filament in the First Cell Cycle Contributes to the Determination of the Anteroposterior Axis in Ascidian Development

In many animal species, the body axis is determined by the relocalization of maternal determinants, organelles, or unique cell populations in a cytoskeleton-dependent manner. In the ascidian first cell cycle, the myoplasm, including mitochondria, endoplasmic reticulum (ER), and maternal mRNAs, move to the future posterior side concomitantly (called ooplasmic segregation or cytoplasmic and cortical reorganization). This translocation consists of first and second phases depending on the actin and microtubule, respectively. However, the transition from first to second phase, that is, translocation of myoplasmic components from microfilaments to microtubules, has been poorly investigated. In this study, we analyzed the relationship between these cytoskeletons and myoplasmic components during the first cell cycle and their role in morphogenesis by inhibitor experiments. Owing to our improved visualization techniques, there was unexpected F-actin accumulation at the vegetal pole during this transition period. When this F-actin was depolymerized, the microtubule structure was strongly affected, the myoplasmic components, including maternal mRNA, were mislocalized, and the anteroposterior axis formation was disordered. These results suggested the importance of F-actin during the first cell cycle and the existence of interactions between microfilaments and microtubules, implying the enigmatic mechanism of ooplasmic segregation. Solving this mystery leads us to an improved understanding of ascidian early development.

Dynamic changes in the association between maternal mRNAs and endoplasmic reticulum during ascidian early embryogenesis

Axis formation is one of the most important events occurring at the beginning of animal development. In the ascidian egg, the antero-posterior axis is established at this time owing to a dynamic cytoplasmic movement called cytoplasmic and cortical reorganisation. During this movement, mitochondria, endoplasmic reticulum (ER), and maternal mRNAs (postplasmic/PEM RNAs) are translocated to the future posterior side. Although accumulating evidence indicates the crucial roles played by the asymmetrical localisation of these organelles and the translational regulation of postplasmic/PEM RNAs, the organisation of ER has not been described in sufficient detail to date owing to technical difficulties. In this study, we developed three different multiple staining protocols for visualising the ER in combination with mitochondria, microtubules, or mRNAs in whole-mount specimens. We defined the internally expanded “dense ER” using these protocols and described cisterna-like structures of the dense ER using focused ion beam-scanning electron microscopy. Most importantly, we described the dynamic changes in the colocalisation of postplasmic/PEM mRNAs and dense ER; for example, macho-1 mRNA was detached and excluded from the dense ER during the second phase of ooplasmic movements. These detailed descriptions of the association between maternal mRNA and ER can provide clues for understanding the translational regulation mechanisms underlying axis determination during ascidian early embryogenesis.

High-Throughput Protein Production Combined with High- Throughput SELEX Identifies an Extensive Atlas of Ciona robusta Transcription Factor DNA-Binding Specificities

Transcription factors (TFs) control gene transcription, binding to specific DNA motifs located in cis-regulatory elements across the genome. The identification of TF-binding motifs is thus an important aspect to understand the role of TFs in gene regulation. SELEX, Systematic Evolution of Ligands by EXponential enrichment, is an efficient in vitro method, which can be used to determine the DNA-binding specificity of TFs. Thanks to the development of high-throughput (HT) DNA cloning system and protein production technology, the classical SELEX assay has be extended to high-throughput scale (HT-SELEX).We report here the detailed protocol for the cloning, production, and purification of 420 Ciona robusta DNA BD. 263 Ciona robusta TF DNA-binding domain proteins were purified in milligram quantities and analyzed by HT-SELEX. The identification of 139 recognition sequences generates an atlas of protein-DNA-binding specificities that is crucial for the understanding of the gene regulatory network (GRN) of Ciona robusta. Overall, our analysis suggests that the Ciona robusta repertoire of sequence-specific transcription factors comprises less than 500 genes. The protocols for high-throughput protein production and HT-SELEX described in this article for the study of Ciona robusta TF DNA-binding specificity are generic and have been successfully applied to a wide range of TFs from other species, including human, mouse, and Drosophila.

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

Since its initial publication in 2002, the genome of Ciona intestinalis type A (Ciona robusta), the first genome sequence of an invertebrate chordate, has provided a valuable resource for a wide range of biological studies, including developmental biology, evolutionary biology, and neuroscience. The genome assembly was updated in 2008, and it included 68% of the sequence information in 14 pairs of chromosomes. However, a more contiguous genome is required for analyses of higher order genomic structure and of chromosomal evolution. Here, we provide a new genome assembly for an inbred line of this animal, constructed with short and long sequencing reads and Hi-C data. In this latest assembly, over 95% of the 123 Mb of sequence data was included in the chromosomes. Short sequencing reads predicted a genome size of 114-120 Mb; therefore, it is likely that the current assembly contains almost the entire genome, although this estimate of genome size was smaller than previous estimates. Remapping of the Hi-C data onto the new assembly revealed a large inversion in the genome of the inbred line. Moreover, a comparison of this genome assembly with that of Ciona savignyi, a different species in the same genus, revealed many chromosomal inversions between these two Ciona species, suggesting that such inversions have occurred frequently and have contributed to chromosomal evolution of Ciona species. Thus, the present assembly greatly improves an essential resource for genome-wide studies of ascidians.

The notochord gene regulatory network in chordate evolution: Conservation and divergence from Ciona to vertebrates

The notochord is a structure required for support and patterning of all chordate embryos, from sea squirts to humans. An increasing amount of information on notochord development and on the molecular strategies that ensure its proper morphogenesis has been gleaned through studies in the sea squirt Ciona. This invertebrate chordate offers a fortunate combination of experimental advantages, ranging from translucent, fast-developing embryos to a compact genome and impressive biomolecular resources. These assets have enabled the rapid identification of numerous notochord genes and cis-regulatory regions, and provide a rather unique opportunity to reconstruct the gene regulatory network that controls the formation of this developmental and evolutionary chordate landmark. This chapter summarizes the morphogenetic milestones that punctuate notochord formation in Ciona, their molecular effectors, and the current knowledge of the gene regulatory network that ensures the accurate spatial and temporal orchestration of these processes.

Repertoires of G protein-coupled receptors for Ciona-specific neuropeptides

Neuropeptides play pivotal roles in various biological events in the nervous, neuroendocrine, and endocrine systems, and are correlated with both physiological functions and unique behavioral traits of animals. Elucidation of functional interaction between neuropeptides and receptors is a crucial step for the verification of their biological roles and evolutionary processes. However, most receptors for novel peptides remain to be identified. Here, we show the identification of multiple G protein-coupled receptors (GPCRs) for species-specific neuropeptides of the vertebrate sister group, Ciona intestinalis Type A, by combining machine learning and experimental validation. We developed an original peptide descriptor-incorporated support vector machine and used it to predict 22 neuropeptide-GPCR pairs. Of note, signaling assays of the predicted pairs identified 1 homologous and 11 Ciona-specific neuropeptide-GPCR pairs for a 41% hit rate: the respective GPCRs for Ci-GALP, Ci-NTLP-2, Ci-LF-1, Ci-LF-2, Ci-LF-5, Ci-LF-6, Ci-LF-7, Ci-LF-8, Ci-YFV-1, and Ci-YFV-3. Interestingly, molecular phylogenetic tree analysis revealed that these receptors, excluding the Ci-GALP receptor, were evolutionarily unrelated to any other known peptide GPCRs, confirming that these GPCRs constitute unprecedented neuropeptide receptor clusters. Altogether, these results verified the neuropeptide-GPCR pairs in the protochordate and evolutionary lineages of neuropeptide GPCRs, and pave the way for investigating the endogenous roles of novel neuropeptides in the closest relatives of vertebrates and the evolutionary processes of neuropeptidergic systems throughout chordates. In addition, the present study also indicates the versatility of the machine-learning-assisted strategy for the identification of novel peptide-receptor pairs in various organisms.

Gene regulatory systems that control gene expression in the Ciona embryo

Transcriptional control of gene expression is one of the most important regulatory systems in animal development. Specific gene expression is basically determined by combinatorial regulation mediated by multiple sequence-specific transcription factors. The decoding of animal genomes has provided an opportunity for us to systematically examine gene regulatory networks consisting of successive layers of control of gene expression. It remains to be determined to what extent combinatorial regulation encoded in gene regulatory networks can explain spatial and temporal gene-expression patterns. The ascidian Ciona intestinalis is one of the animals in which the gene regulatory network has been most extensively studied. In this species, most specific gene expression patterns in the embryo can be explained by combinations of upstream regulatory genes encoding transcription factors and signaling molecules. Systematic scrutiny of gene expression patterns and regulatory interactions at the cellular resolution have revealed incomplete parts of the network elucidated so far, and have identified novel regulatory genes and novel regulatory mechanisms.

The identification of transcription factors expressed in the notochord of Ciona intestinalis adds new potential players to the brachyury gene regulatory network

The notochord is the distinctive characteristic of chordates; however, the knowledge of the complement of transcription factors governing the development of this structure is still incomplete. Here we present the expression patterns of seven transcription factor genes detected in the notochord of the ascidian Ciona intestinalis at various stages of embryonic development. Four of these transcription factors, Fos-a, NFAT5, AFF and Klf15, have not been directly associated with the notochord in previous studies, while the others, including Spalt-like-a, Lmx-like, and STAT5/6-b, display evolutionarily conserved expression in this structure as well as in other domains. We examined the hierarchical relationships between these genes and the transcription factor Brachyury, which is necessary for notochord development in all chordates. We found that Ciona Brachyury regulates the expression of most, although not all, of these genes. These results shed light on the genetic regulatory program underlying notochord formation in Ciona and possibly other chordates.

Comparative cytogenetic mapping of Sox2 and Sox14 in cichlid fishes and inferences on the genomic organization of both genes in vertebrates

To better understand the genomic organization and evolution of Sox genes in vertebrates, we cytogenetically mapped Sox2 and Sox14 genes in cichlid fishes and performed comparative analyses of their orthologs in several vertebrate species. The genomic regions neighboring Sox2 and Sox14 have been conserved during vertebrate diversification. Although cichlids seem to have undergone high rates of genomic rearrangements, Sox2 and Sox14 are linked in the same chromosome in the Etroplinae Etroplus maculatus that represents the sister group of all remaining cichlids. However, these genes are located on different chromosomes in several species of the sister group Pseudocrenilabrinae. Similarly, the ancestral synteny of Sox2 and Sox14 has been maintained in several vertebrates, but this synteny has been broken independently in all major groups as a consequence of karyotype rearrangements that took place during the vertebrate evolution.

Ciona intestinalis and Oxycomanthus japonicus, representatives of marine invertebrates

The study of marine invertebrates is useful in various biological research fields. However, genetic analyses of these animals are limited, mainly due to difficulties in culturing them, and the genetic resources of marine invertebrates have not been organized. Recently, advances have been made in the study of two deuterostomes, an ascidian Ciona intestinalis and a feather star Oxycomanthus japonicus. The draft genome sequence of Ciona intestinalis has been determined, and its compact genome, which has less redundancy of genes compared with vertebrates, provides us with a useful experimental system for analyzing the functions of genes during development. The life cycle of Ciona intestinalis is approximately 2-3 months, and the genetic techniques including a perfect inland culture system, germline transformation with a transposon Minos, enhancer detection and insertional mutagenesis, have been established. The feather star Oxycomanthus japonicus conserves the characteristics of the basic echinoderm body plan with a segmented mesoderm, which is a fascinating characteristic for understanding the evolution of echinoderms. Oxycomanthus japonicus shows strong regeneration ability and is a suitable subject for analysis of the mechanisms of regeneration. In consideration of these features, the National BioResource Project (NBRP) has started to support the supply of wild-types, transgenic lines and inbred lines of Ciona intestinalis and Oxycomanthus japonicus.

Early zygotic expression of transcription factors and signal molecules in fully dissociated embryonic cells of Ciona intestinalis: A microarray analysis

Specification of early embryonic cells of animals is established by maternally provided factors and interactions of neighboring cells. The present study addressed a question of autonomous versus non-autonomous specification of embryonic cells by using the Ciona intestinalis embryo, in particular the genetic cascade of zygotic expression of transcription factor genes responsible for notochord specification. To examine this issue, we combined the classic experiment of continuous dissociation of embryonic cells with the modern technique of oligonucleotide-based microarrays. We measured early zygotic expression of 389 core transcription factors genes and 118 major signal molecule genes in embryonic cells that were fully dissociated from the first cleavage. Our results indicated that even if cells are free from contact with neighbors, the major transcription factor genes that have primary roles in embryonic cell specification commence their zygotic expression at the same time as in normal embryos. Dissociation of embryonic cells did not affect extracellular signal-regulated kinases (ERK) activity. Although normal embryos treated with U0126 failed to express Bra and Twist-like-1, dissociated embryonic cells treated with U0126 expressed the genes. These results are discussed in relation to the grade of autonomous versus non-autonomous genetic cascades that are responsible for the specification of early Ciona embryonic cells.

Evolution of non-coding regulatory sequences involved in the developmental process: reflection of differential employment of paralogous genes as highlighted by Sox2 and group B1 Sox genes

In higher vertebrates, the expression of Sox2, a group B1 Sox gene, is the hallmark of neural primordial cell state during the developmental processes from embryo to adult. Sox2 is regulated by the combined action of many enhancers with distinct spatio-temporal specificities. DNA sequences for these enhancers are conserved in a wide range of vertebrate species, corresponding to a majority of highly conserved non-coding sequences surrounding the Sox2 gene, corroborating the notion that the conservation of non-coding sequences mirrors their functional importance. Among the Sox2 enhancers, N-1 and N-2 are activated the earliest in embryogenesis and regulate Sox2 in posterior and anterior neural plates, respectively. These enhancers differ in their evolutionary history: the sequence and activity of enhancer N-2 is conserved in all vertebrate species, while enhancer N-1 is fully conserved only in amniotes. In teleost embryos, Sox19a/b play the major pan-neural role among the group B1 Sox paralogues, while strong Sox2 expression is limited to the anterior neural plate, reflecting the absence of posterior CNS-dedicated enhancers, including N-1. In Xenopus, neurally expressed SoxD is the orthologue of Sox19, but Sox3 appears to dominate other B1 paralogues. In amniotes, however, Sox19 has lost its group B1 Sox function and transforms into group G Sox15 (neofunctionalization), and Sox2 assumes the dominant position by gaining enhancer N-1 and other enhancers for posterior CNS. Thus, the gain and loss of specific enhancer elements during the evolutionary process reflects the change in functional assignment of particular paralogous genes, while overall regulatory functions attributed to the gene family are maintained.

Improved genome assembly and evidence-based global gene model set for the chordate Ciona intestinalis: new insight into intron and operon populations

An improved assembly of the Ciona intestinalis genome reveals that it contains non-canonical introns and that about 20% of Ciona genes reside in operons.

Background

The draft genome sequence of the ascidian Ciona intestinalis, along with associated gene models, has been a valuable research resource. However, recently accumulated expressed sequence tag (EST)/cDNA data have revealed numerous inconsistencies with the gene models due in part to intrinsic limitations in gene prediction programs and in part to the fragmented nature of the assembly.

Results

We have prepared a less-fragmented assembly on the basis of scaffold-joining guided by paired-end EST and bacterial artificial chromosome (BAC) sequences, and BAC chromosomal in situ hybridization data. The new assembly (115.2 Mb) is similar in length to the initial assembly (116.7 Mb) but contains 1,272 (approximately 50%) fewer scaffolds. The largest scaffold in the new assembly incorporates 95 initial-assembly scaffolds. In conjunction with the new assembly, we have prepared a greatly improved global gene model set strictly correlated with the extensive currently available EST data. The total gene number (15,254) is similar to that of the initial set (15,582), but the new set includes 3,330 models at genomic sites where none were present in the initial set, and 1,779 models that represent fusions of multiple previously incomplete models. In approximately half, 5'-ends were precisely mapped using 5'-full-length ESTs, an important refinement even in otherwise unchanged models.

Conclusion

Using these new resources, we identify a population of non-canonical (non-GT-AG) introns and also find that approximately 20% of Ciona genes reside in operons and that operons contain a high proportion of single-exon genes. Thus, the present dataset provides an opportunity to analyze the Ciona genome much more precisely than ever.

Logic of gene regulatory networks

Regulatory networks of transcription factors and signaling molecules lie at the heart of development. Their architecture implements logic functions whose execution propels cells from one regulatory state to the next, thus driving development forward. As an example of a subcircuit that translates transcriptional input into developmental output, we consider a particularly simple case, the regulatory processes underlying pigment cell formation in sea urchin embryos. The regulatory events in this process can be represented as elementary logic functions.

PlanTAPDB, a phylogeny-based resource of plant transcription-associated proteins

Diversification of transcription-associated protein (TAP) families during land plant evolution is a key process yielding increased complexity of plant life. Understanding the evolutionary relationships between these genes is crucial to gain insight into plant evolution. We have determined a substantial set of TAPs that are focused on, but not limited to, land plants using PSI-BLAST searches and subsequent filtering and clustering steps. Phylogenies were created in an automated way using a combination of distance and maximum likelihood methods. Comparison of the data to previously published work confirmed their accuracy and usefulness for the majority of gene families. Evidence is presented that the flowering plant apical stem cell regulator WUSCHEL evolved from an ancestral homeobox gene that was already present after the water-to-land transition. The presence of distinct expanded gene families, such as COP1 and HIT in moss, is discussed within the evolutionary backdrop. Comparative analyses revealed that almost all angiosperm transcription factor families were already present in the earliest land plants, whereas many are missing among unicellular algae. A global analysis not only of transcription factors but also of transcriptional regulators and novel putative families is presented. A wealth of data about plant TAP families and all data accrued throughout their automated detection and analysis are made available via the PlanTAPDB Web interface. Evolutionary relationships of these genes are readily accessible to the nonexpert at a mouse-click. Initial analyses of selected gene families revealed that PlanTAPDB can easily be exerted for knowledge discovery.

Expression of the AmphiTcf gene in amphioxus: insights into the evolution of the TCF/LEF gene family during vertebrate evolution

T-cell factor (TCF) and lymphoid enhancer factors (LEF) genes encode proteins that are transcription factors mediating beta-catenin/Wnt signaling. Whereas mammals have four such genes, the Florida amphioxus (Branchiostoma floridae) apparently has only one such gene (AmphiTcf). From cleavage through early gastrula, cytoplasmic maternal transcripts of this gene are localized toward the animal pole. In gastrulae, AmphiTcf expression begins in the mesendoderm. In neurulae, there is expression in the pharynx, hindgut, anterior notochord, somites, and at the anterior end of the neural plate. In early larvae, expression is detectable in the floor of the diencephalon, notochord, tail bud, forming somites, pharynx, and ciliated pit (a presumed homolog of the vertebrate adenohypophysis). Phylogenetic analysis of TCF/LEF proteins placed AmphiTcf as the sister group of a clade comprising vertebrate Tcf1, Lef1, Tcf3, and Tcf4. Comparison of developmental expression for amphioxus AmphiTcf and vertebrate TCF/LEF genes indicates that this gene family has undergone extensive subfunctionalization and neofunctionalization during vertebrate evolution.

Microarray analysis of zygotic expression of transcription factor genes and cell signaling molecule genes in early Ciona intestinalis embryos

In ascidians, specification of embryonic cells takes place very early at the 16-, 32- and 64-cell stages, and this developmental event involves zygotic expression of various genes, some encoding transcription factors and some encoding cell signaling molecules. Previous studies have demonstrated that approximately 50 transcription factor genes and 25 signaling molecule genes commence their zygotic expression by the 64-cell stage of Ciona intestinalis embryos. With the aid of oligonucleotide-based microarray, we examined the zygotic expression profiles of developmental genes in early Ciona embryos. Although the microarray method had a tendency to barely detect zygotic expression of genes that are expressed maternally, the present results confirmed the results of previous studies. In addition, the present analysis demonstrated the zygotic expression of four genes that were not identified in previous studies, and this result was confirmed by whole-mount in situ hybridization. Our results therefore provide further information on the developmental genes that are zygotically expressed in early Ciona embryos, and demonstrate that the microarray is a powerful tool for future studies of the gene regulatory network in Ciona, a basal chordate with a body plan similar to that of vertebrates.

Systematic analysis of embryonic expression profiles of zinc finger genes in Ciona intestinalis

The recent decoding of a number of animal genomes has provided unprecedented information regarding evolution and gene structures, but this information must be supplemented with precise gene annotations and the temporal and spatial expression patterns of individual genes. In the present study, we systematically identified and characterized 566 zinc finger genes in the genome of Ciona intestinalis, an emerging model system for genome-wide studies of development and evolution. Of these genes, 356 genes encoded a potential transcription factor based on putative nucleic acid binding activity or domains of unknown function. We further examined the expression patterns of 225 genes during embryogenesis, and, when considered with a previous study [Imai, K.S., Hino, K., Yagi, K., Satoh, N., Satou, Y., 2004. Gene expression profiles of transcription factors and signaling molecules in the ascidian embryo: towards a comprehensive understanding of gene networks. Development 131, 4047-4058], we have characterized the developmental expression patterns of nearly 85% of the potential zinc finger-containing transcription factors. Overall, zinc finger genes are preferentially maternally expressed with little larval expression during development. The present study provides a valuable reference for genome-wide studies in this species and for future studies wishing to examine zinc finger gene expression patterns in other animals.