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

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.

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

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

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.

Genome-Wide Identification, Comparison, and Expression Analysis of Transcription Factors in Ascidian Styela clava

Tunicates include diverse species, as they are model animals for evolutionary developmental biology study. The embryonic development of tunicates is known to be extensively regulated by transcription factors (TFs). Styela clava, the globally distributed invasive tunicate, exhibits a strong capacity for environmental adaptation. However, the TFs were not systematically identified and analyzed. In this study, we reported 553 TFs categorized into 60 families from S. clava, based on the whole genome data. Comparison of TFs analysis among the tunicate species revealed that the gene number in the zinc finger superfamily displayed the most significant discrepancy, indicating this family was under the highly evolutionary selection and might be related to species differentiation and environmental adaptation. The greatest number of TFs was discovered in the Cys2His2-type zinc finger protein (zf-C2H2) family in S. clava. From the point of temporal view, more than half the TFs were expressed at the early embryonic stage. The expression correlation analysis revealed the existence of a transition for TFs expression from early embryogenesis to the later larval development in S. clava. Eight Hox genes were identified to be located on one chromosome, exhibiting different arrangement and expression patterns, compared to Ciona robusta (C. intestinalis type A). In addition, a total of 23 forkhead box (fox) genes were identified in S. clava, and their expression profiles referred to their potential roles in neurodevelopment and sensory organ development. Our data, thus, provides crucial clues to the potential functions of TFs in development and environmental adaptation in the leathery sea squirt.

Transcription Factors of the bHLH Family Delineate Vertebrate Landmarks in the Nervous System of a Simple Chordate

Tunicates are marine invertebrates whose tadpole-like larvae feature a highly simplified version of the chordate body plan. Similar to their distant vertebrate relatives, tunicate larvae develop a regionalized central nervous system and form distinct neural structures, which include a rostral sensory vesicle, a motor ganglion, and a caudal nerve cord. The sensory vesicle contains a photoreceptive complex and a statocyst, and based on the comparable expression patterns of evolutionarily conserved marker genes, it is believed to include proto-hypothalamic and proto-retinal territories. The evolutionarily conserved molecular fingerprints of these landmarks of the vertebrate brain consist of genes encoding for different transcription factors, and of the gene batteries that they control, and include several members of the bHLH family. Here we review the complement of bHLH genes present in the streamlined genome of the tunicate Ciona robusta and their current classification, and summarize recent studies on proneural bHLH transcription factors and their expression territories. We discuss the possible roles of bHLH genes in establishing the molecular compartmentalization of the enticing nervous system of this unassuming chordate.

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.

The genetic program to specify ectodermal cells in ascidian embryos

The ascidian belongs to the sister group of vertebrates and shares many features with them. The gene regulatory network (GRN) controlling gene expression in ascidian embryonic development leading to the tadpole larva has revealed evolutionarily conserved gene circuits between ascidians and vertebrates. These conserved mechanisms are indeed useful to infer the original developmental programs of the ancestral chordates. Simultaneously, these studies have revealed which gene circuits are missing in the ascidian GRN; these gene circuits may have been acquired in the vertebrate lineage. In particular, the GRN responsible for gene expression in ectodermal cells of ascidian embryos has revealed the genetic programs that regulate the regionalization of the brain, formation of palps derived from placode-like cells, and differentiation of sensory neurons derived from neural crest-like cells. We here discuss how these studies have given insights into the evolution of these traits.

Finding cell-specific expression patterns in the early Ciona embryo with single-cell RNA-seq

Single-cell RNA-seq has been established as a reliable and accessible technique enabling new types of analyses, such as identifying cell types and studying spatial and temporal gene expression variation and change at single-cell resolution. Recently, single-cell RNA-seq has been applied to developing embryos, which offers great potential for finding and characterising genes controlling the course of development along with their expression patterns. In this study, we applied single-cell RNA-seq to the 16-cell stage of the Ciona embryo, a marine chordate and performed a computational search for cell-specific gene expression patterns. We recovered many known expression patterns from our single-cell RNA-seq data and despite extensive previous screens, we succeeded in finding new cell-specific patterns, which we validated by in situ and single-cell qPCR.

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.

Znf703 is a novel RA target in the neural plate border

Znf703 is an RAR- and Wnt-inducible transcription factor that exhibits a complex expression pattern in the developing embryo: Znf703 mRNA is found in the early circumblastoporal ring, then later throughout the neural plate and its border, and subsequently in the mid/hindbrain and somites. We show that Znf703 has a different and separable function in early mesoderm versus neural crest and placode development. Independent of its early knockdown phenotype on Gdf3 and Wnt8, Znf703 disrupts patterning of distinct neural crest migratory streams normally delineated by Sox10, Twist, and Foxd3 and inhibits otocyst formation and otic expression of Sox10 and Eya1. Furthermore, Znf703 promotes massive overgrowth of SOX2+ cells, disrupting the SoxB1 balance at the neural plate border. Despite prominent expression in other neural plate border-derived cranial and sensory domains, Znf703 is selectively absent from the otocyst, suggesting that Znf703 must be specifically cleared or down-regulated for proper otic development. We show that mutation of the putative Groucho-repression domain does not ameliorate Znf703 effects on mesoderm, neural crest, and placodes. We instead provide evidence that Znf703 requires the Buttonhead domain for transcriptional repression.

Antero-posterior ectoderm patterning by canonical Wnt signaling during ascidian development

Wnt/β-catenin signaling is an ancient pathway in metazoans and controls various developmental processes, in particular the establishment and patterning of the embryonic primary axis. In vertebrates, a graded Wnt activity from posterior to anterior endows cells with positional information in the central nervous system. Recent studies in hemichordates support a conserved role for Wnt/β-catenin in ectoderm antero-posterior patterning at the base of the deuterostomes. Ascidians are marine invertebrates and the closest relatives of vertebrates. By combining gain- and loss-of-function approaches, we have determined the role of Wnt/β-catenin in patterning the three ectoderm derivatives of the ascidian Ciona intestinalis, central nervous system, peripheral nervous system and epidermis. Activating Wnt/β-catenin signaling from gastrulation led to a dramatic transformation of the ectoderm with a loss of anterior identities and a reciprocal anterior extension of posterior identities, consistent with studies in other metazoans. Surprisingly, inhibiting Wnt signaling did not produce a reciprocal anteriorization of the embryo with a loss of more posterior identities like in vertebrates and hemichordate. Epidermis patterning was overall unchanged. Only the identity of two discrete regions of the central nervous system, the anteriormost and the posteriormost regions, were under the control of Wnt. Finally, the caudal peripheral nervous system, while being initially Wnt dependent, formed normally. Our results show that the Ciona embryonic ectoderm responds to Wnt activation in a manner that is compatible with the proposed function for this pathway at the base of the deuterostomes. However, possibly because of its fast and divergent mode of development that includes extensive use of maternal determinants, the overall antero-posterior patterning of the Ciona ectoderm is Wnt independent, and Wnt/β-catenin signaling controls the formation of some sub-domains. Our results thus indicate that there has likely been a drift in the developmental systems controlling ectoderm patterning in the lineage leading to ascidians.

The Wnt/β-catenin pathway is a system of cell-cell communication. It has an ancient origin in animals and plays multiple roles during embryogenesis and adult life. In particular, it is involved in determining, in the vertebrate embryo, the identity of the different parts of the body and their relative positions along the antero-posterior axis. We have investigated in an ascidian (or sea squirt) species, a marine invertebrate that is closely related to vertebrates, whether this pathway had a similar role. Like in vertebrates, activating Wnt/β-catenin led to a posteriorization of the embryo with a loss of anterior structures. By contrast, unlike vertebrates, ascidian embryos formed rather normally following Wnt/β-catenin inactivation. Since hemichordates (or acorn worms), earlier divergent invertebrates, use Wnt/β-catenin in a manner comparable to vertebrates, it is in the ascidian lineage that changes have occurred. Consequently, ascidians build an antero-posterior axis, very similarly organized to that of vertebrates, but in a different way.

Positioning a multifunctional basic helix-loop-helix transcription factor within the Ciona notochord gene regulatory network

In a multitude of organisms, transcription factors of the basic helix-loop-helix (bHLH) family control the expression of genes required for organ development and tissue differentiation. The functions of different bHLH transcription factors in the specification of nervous system and paraxial mesoderm have been widely investigated in various model systems. Conversely, the knowledge of the role of these regulators in the development of the axial mesoderm, the embryonic territory that gives rise to the notochord, and the identities of their target genes, remain still fragmentary. Here we investigated the transcriptional regulation and target genes of Bhlh-tun1, a bHLH transcription factor expressed in the developing Ciona notochord as well as in additional embryonic territories that contribute to the formation of both larval and adult structures. We describe its possible role in notochord formation, its relationship with the key notochord transcription factor Brachyury, and suggest molecular mechanisms through which Bhlh-tun1 controls the spatial and temporal expression of its effectors.

Distinct regulation of Snail in two muscle lineages of the ascidian embryo achieves temporal coordination of muscle development

The transcriptional repressor Snail is required for proper differentiation of the tail muscle of ascidian tadpole larvae. Two muscle lineages (B5.1 and B6.4) contribute to the anterior tail muscle cells, and are consecutively separated from a transcriptionally quiescent germ cell lineage at the 16- and 32-cell stages. Concomitantly, cells of these lineages begin to express Tbx6.b (Tbx6-r.b) at the 16- and 32-cell stages, respectively. Meanwhile, Snail expression begins in these two lineages simultaneously at the 32-cell stage. Here, we show that Snail expression is regulated differently between these two lineages. In the B5.1 lineage, Snail was activated through Tbx6.b, which is activated by maternal factors, including Zic-r.a. In the B6.4 lineage, the MAPK pathway was cell-autonomously activated by a constitutively active form of Raf, enabling Zic-r.a to activate Snail independently of Tbx6.b As a result, Snail begins to be expressed at the 32-cell stage simultaneously in these two lineages. Such shortcuts might be required for coordinating developmental programs in embryos in which cells become separated progressively from stem cells, including germline cells.

Organ specific gene expression in the regenerating tail of Macrostomum lignano

Temporal and spatial characterization of gene expression is a prerequisite for the understanding of cell-, tissue-, and organ-differentiation. In a multifaceted approach to investigate gene expression in the tail plate of the free-living marine flatworm Macrostomum lignano, we performed a posterior-region-specific in situ hybridization screen, RNA sequencing (RNA-seq) of regenerating animals, and functional analyses of selected tail-specific genes. The in situ screen revealed transcripts expressed in the antrum, cement glands, adhesive organs, prostate glands, rhabdite glands, and other tissues. Next we used RNA-seq to characterize temporal expression in the regenerating tail plate revealing a time restricted onset of both adhesive organs and copulatory apparatus regeneration. In addition, we identified three novel previously unannotated genes solely expressed in the regenerating stylet. RNA interference showed that these genes are required for the formation of not only the stylet but the whole male copulatory apparatus. RNAi treated animals lacked the stylet, vesicula granulorum, seminal vesicle, false seminal vesicle, and prostate glands, while the other tissues of the tail plate, such as adhesive organs regenerated normally. In summary, our findings provide a large resource of expression data during homeostasis and regeneration of the morphologically complex tail regeneration and pave the way for a better understanding of organogenesis in M. lignano.

How complexity increases in development: An analysis of the spatial-temporal dynamics of Gene expression in Ciona intestinalis

The increase in complexity in an embryo over developmental time is perhaps one of the most intuitive processes of animal development. It is also intuitive that the embryo becomes progressively compartmentalized over time and space. In spite of this intuitiveness, there are no systematic attempts to quantify how this occurs. Here, we present a quantitative analysis of the compartmentalization and spatial complexity of Ciona intestinalis over developmental time by analyzing thousands of gene expression spatial patterns from the ANISEED database. We measure compartmentalization in two ways: as the relative volume of expression of genes and as the disparity in gene expression between body parts. We also use a measure of the curvature of each gene expression pattern in 3D space. These measures show a similar increase over time, with the most dramatic change occurring from the 112-cell stage to the early tailbud stage. Combined, these measures point to a global pattern of increase in complexity in the Ciona embryo. Finally, we cluster the different regions of the embryo depending on their gene expression similarity, within and between stages. Results from this clustering analysis, which partially correspond to known fate maps, provide a global quantitative overview about differentiation and compartmentalization between body parts at each developmental stage.

Brachyury, Foxa2 and the cis-Regulatory Origins of the Notochord

A main challenge of modern biology is to understand how specific constellations of genes are activated to differentiate cells and give rise to distinct tissues. This study focuses on elucidating how gene expression is initiated in the notochord, an axial structure that provides support and patterning signals to embryos of humans and all other chordates. Although numerous notochord genes have been identified, the regulatory DNAs that orchestrate development and propel evolution of this structure by eliciting notochord gene expression remain mostly uncharted, and the information on their configuration and recurrence is still quite fragmentary. Here we used the simple chordate Ciona for a systematic analysis of notochord cis-regulatory modules (CRMs), and investigated their composition, architectural constraints, predictive ability and evolutionary conservation. We found that most Ciona notochord CRMs relied upon variable combinations of binding sites for the transcription factors Brachyury and/or Foxa2, which can act either synergistically or independently from one another. Notably, one of these CRMs contains a Brachyury binding site juxtaposed to an (AC) microsatellite, an unusual arrangement also found in Brachyury-bound regulatory regions in mouse. In contrast, different subsets of CRMs relied upon binding sites for transcription factors of widely diverse families. Surprisingly, we found that neither intra-genomic nor interspecific conservation of binding sites were reliably predictive hallmarks of notochord CRMs. We propose that rather than obeying a rigid sequence-based cis-regulatory code, most notochord CRMs are rather unique. Yet, this study uncovered essential elements recurrently used by divergent chordates as basic building blocks for notochord CRMs.

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.

Molecular biology approaches in bioadhesion research

The use of molecular biology tools in the field of bioadhesion is still in its infancy. For new research groups who are considering taking a molecular approach, the techniques presented here are essential to unravelling the sequence of a gene, its expression and its biological function. Here we provide an outline for addressing adhesion-related genes in diverse organisms. We show how to gradually narrow down the number of candidate transcripts that are involved in adhesion by (1) generating a transcriptome and a differentially expressed cDNA list enriched for adhesion-related transcripts, (2) setting up a BLAST search facility, (3) perform an in situ hybridization screen, and (4) functional analyses of selected genes by using RNA interference knock-down. Furthermore, latest developments in genome-editing are presented as new tools to study gene function. By using this iterative multi-technologies approach, the identification, isolation, expression and function of adhesion-related genes can be studied in most organisms. These tools will improve our understanding of the diversity of molecules used for adhesion in different organisms and these findings will help to develop innovative bio-inspired adhesives.

Anterior-posterior regionalized gene expression in the Ciona notochord

BACKGROUND

In the simple ascidian chordate Ciona, the signaling pathways and gene regulatory networks giving rise to initial notochord induction are largely understood and the mechanisms of notochord morphogenesis are being systematically elucidated. The notochord has generally been thought of as a non-compartmentalized or regionalized organ that is not finely patterned at the level of gene expression. Quantitative imaging methods have recently shown, however, that notochord cell size, shape, and behavior vary consistently along the anterior-posterior (AP) axis.

RESULTS

Here we screen candidate genes by whole mount in situ hybridization for potential AP asymmetry. We identify 4 genes that show non-uniform expression in the notochord. Ezrin/radixin/moesin (ERM) is expressed more strongly in the secondary notochord lineage than the primary. CTGF is expressed stochastically in a subset of notochord cells. A novel calmodulin-like gene (BCamL) is expressed more strongly at both the anterior and posterior tips of the notochord. A TGF-β ortholog is expressed in a gradient from posterior to anterior. The asymmetries in ERM, BCamL, and TGF-β expression are evident even before the notochord cells have intercalated into a single-file column.

CONCLUSIONS

We conclude that the Ciona notochord is not a homogeneous tissue but instead shows distinct patterns of regionalized gene expression.

Tbx2/3 is an essential mediator within the Brachyury gene network during Ciona notochord development

T-box genes are potent regulators of mesoderm development in many metazoans. In chordate embryos, the T-box transcription factor Brachyury (Bra) is required for specification and differentiation of the notochord. In some chordates, including the ascidian Ciona, members of the Tbx2 subfamily of T-box genes are also expressed in this tissue; however, their regulatory relationships with Bra and their contributions to the development of the notochord remain uncharacterized. We determined that the notochord expression of Ciona Tbx2/3 (Ci-Tbx2/3) requires Ci-Bra, and identified a Ci-Tbx2/3 notochord CRM that necessitates multiple Ci-Bra binding sites for its activity. Expression of mutant forms of Ci-Tbx2/3 in the developing notochord revealed a role for this transcription factor primarily in convergent extension. Through microarray screens, we uncovered numerous Ci-Tbx2/3 targets, some of which overlap with known Ci-Bra-downstream notochord genes. Among the Ci-Tbx2/3 notochord targets are evolutionarily conserved genes, including caspases, lineage-specific genes, such as Noto4, and newly identified genes, such as MLKL. This work sheds light on a large section of the notochord regulatory circuitry controlled by T-box factors, and reveals new components of the complement of genes required for the proper formation of this structure.

Differential gene expression in notochord and nerve cord fate segregation in the Ciona intestinalis embryo

During early embryogenesis, embryonic cells gradually restrict their developmental potential and are eventually destined to give rise to one type of cells. Molecular mechanisms underlying developmental fate restriction are one of the major research subjects within developmental biology. In this article, this subject was addressed by combining blastomere isolation with microarray analysis. During the 6th cleavage of the Ciona intestinalis embryo, from the 32-cell to the 64-cell stage, four mother cells divide into daughter cells with two distinct fates, one giving rise to notochord precursor cells and the other to nerve cord precursors. Approximately 2,200 each of notochord and nerve cord precursor cells were isolated, and their mRNA expression profiles were compared by microarray. This analysis identified 106 and 68 genes, respectively, that are differentially expressed in notochord and nerve cord precursor cells. These included not only genes for transcription factors and signaling molecules but also those with generalized functions observed in many types of cells. In addition, whole-mount in situ hybridization showed dynamic spatial expression profiles of these genes during segregation of the two fates: partitioning of transcripts present in the mother cells into either type of daughter cells, and initiation of preferential gene expression in either type of cells.

Highly conserved elements discovered in vertebrates are present in non-syntenic loci of tunicates, act as enhancers and can be transcribed during development

Co-option of cis-regulatory modules has been suggested as a mechanism for the evolution of expression sites during development. However, the extent and mechanisms involved in mobilization of cis-regulatory modules remains elusive. To trace the history of non-coding elements, which may represent candidate ancestral cis-regulatory modules affirmed during chordate evolution, we have searched for conserved elements in tunicate and vertebrate (Olfactores) genomes. We identified, for the first time, 183 non-coding sequences that are highly conserved between the two groups. Our results show that all but one element are conserved in non-syntenic regions between vertebrate and tunicate genomes, while being syntenic among vertebrates. Nevertheless, in all the groups, they are significantly associated with transcription factors showing specific functions fundamental to animal development, such as multicellular organism development and sequence-specific DNA binding. The majority of these regions map onto ultraconserved elements and we demonstrate that they can act as functional enhancers within the organism of origin, as well as in cross-transgenesis experiments, and that they are transcribed in extant species of Olfactores. We refer to the elements as 'Olfactores conserved non-coding elements'.

Correlating structure and function during the evolution of fibrinogen-related domains

Fibrinogen-related domains (FReDs) are found in a variety of animal proteins with widely different functions, ranging from non-self recognition to clot formation. All appear to have a common surface where binding of one sort or other occurs. An examination of 19 completed animal genomes--including a sponge and sea anemone, six protostomes, and 11 deuterostomes--has allowed phylogenies to be constructed that show where various types of FReP (proteins containing FReDs) first made their appearance. Comparisons of sequences and structures also reveal particular features that correlate with function, including the influence of neighbor-domains. A particular set of insertions in the carboxyl-terminal subdomain was involved in the transition from structures known to bind sugars to those known to bind amino-terminal peptides. Perhaps not unexpectedly, FReDs with different functions have changed at different rates, with ficolins by far the fastest changing group. Significantly, the greatest amount of change in ficolin FReDs occurs in the third subdomain ("P domain"), the very opposite of the situation in most other vertebrate FReDs. The unbalanced style of change was also observed in FReDs from non-chordates, many of which have been implicated in innate immunity.

Transcriptional enhancers in ascidian development

The study of cis-regulatory DNAs that control developmental gene expression is integral to the modeling of comprehensive genomic regulatory networks for embryogenesis. Ascidian embryos provide a unique opportunity for the analysis of cis-regulatory DNAs with cellular resolution in the context of a simple but typical chordate body plan. Here, we review landmark studies that have laid the foundations for the study of transcriptional enhancers, among other cis-regulatory DNAs, and their roles in ascidian development. The studies using ascidians of the Ciona genus have capitalized on a unique electroporation technique that permits the simultaneous transfection of hundreds of fertilized eggs, which develop rapidly and express transgenes with little mosaicism. Current studies using the ascidian embryo benefit from extensively annotated genomic resources to characterize transcript models in silico. The search for functional noncoding sequences can be guided by bioinformatic analyses combining evolutionary conservation, gene coexpression, and combinations of overrepresented short-sequence motifs. The power of the transient transfection assays has allowed thorough dissection of numerous cis-regulatory modules, which provided insights into the functional constraints that shape enhancer architecture and diversification. Future studies will benefit from pioneering stable transgenic lines and the analysis of chromatin states. Whole genome expression, functional and DNA binding data are being integrated into comprehensive genomic regulatory network models of early ascidian cell specification with a single-cell resolution that is unique among chordate model systems.

Sequential activation of apical and basolateral contractility drives ascidian endoderm invagination

BACKGROUND

Epithelial invagination is a fundamental morphogenetic behavior that transforms a flat cell sheet into a pit or groove. Previous studies of invagination have focused on the role of actomyosin-dependent apical contraction; other mechanisms remain largely unexplored.

RESULTS

We combined experimental and computational approaches to identify a two-step mechanism for endoderm invagination during ascidian gastrulation. During Step 1, which immediately precedes invagination, endoderm cells constrict their apices because of Rho/Rho-kinase-dependent apical enrichment of 1P-myosin. Our data suggest that endoderm invagination itself occurs during Step 2, without further apical shrinkage, via a novel mechanism we call collared rounding: Rho/Rho-kinase-independent basolateral enrichment of 1P-myosin drives apico-basal shortening, whereas Rho/Rho-kinase-dependent enrichment of 1P and 2P myosin in circumapical collars is required to prevent apical expansion and for deep invagination. Simulations show that boundary-specific tension values consistent with these distributions of active myosin can explain the cell shape changes observed during invagination both in normal embryos and in embryos treated with pharmacological inhibitors of either Rho-kinase or Myosin II ATPase. Indeed, we find that the balance of strong circumapical and basolateral tension is the only mechanism based on differential cortical tension that can explain ascidian endoderm invagination. Finally, simulations suggest that mesectoderm cells resist endoderm shape changes during both steps, and we confirm this prediction experimentally.

CONCLUSIONS

Our findings suggest that early ascidian gastrulation is driven by the coordinated apposition of circumapical and lateral endoderm contraction, working against a resisting mesectoderm. We propose that similar mechanisms may operate during other invaginations.

The ANISEED database: digital representation, formalization, and elucidation of a chordate developmental program

Developmental biology aims to understand how the dynamics of embryonic shapes and organ functions are encoded in linear DNA molecules. Thanks to recent progress in genomics and imaging technologies, systemic approaches are now used in parallel with small-scale studies to establish links between genomic information and phenotypes, often described at the subcellular level. Current model organism databases, however, do not integrate heterogeneous data sets at different scales into a global view of the developmental program. Here, we present a novel, generic digital system, NISEED, and its implementation, ANISEED, to ascidians, which are invertebrate chordates suitable for developmental systems biology approaches. ANISEED hosts an unprecedented combination of anatomical and molecular data on ascidian development. This includes the first detailed anatomical ontologies for these embryos, and quantitative geometrical descriptions of developing cells obtained from reconstructed three-dimensional (3D) embryos up to the gastrula stages. Fully annotated gene model sets are linked to 30,000 high-resolution spatial gene expression patterns in wild-type and experimentally manipulated conditions and to 528 experimentally validated cis-regulatory regions imported from specialized databases or extracted from 160 literature articles. This highly structured data set can be explored via a Developmental Browser, a Genome Browser, and a 3D Virtual Embryo module. We show how integration of heterogeneous data in ANISEED can provide a system-level understanding of the developmental program through the automatic inference of gene regulatory interactions, the identification of inducing signals, and the discovery and explanation of novel asymmetric divisions.

Genomic cis-regulatory networks in the early Ciona intestinalis embryo

Precise spatiotemporal gene expression during animal development is achieved through gene regulatory networks, in which sequence-specific transcription factors (TFs) bind to cis-regulatory elements of target genes. Although numerous cis-regulatory elements have been identified in a variety of systems, their global architecture in the gene networks that regulate animal development is not well understood. Here, we determined the structure of the core networks at the cis-regulatory level in early embryos of the chordate Ciona intestinalis by chromatin immunoprecipitation (ChIP) of 11 TFs. The regulatory systems of the 11 TF genes examined were tightly interconnected with one another. By combining analysis of the ChIP data with the results of previous comprehensive analyses of expression profiles and knockdown of regulatory genes, we found that most of the previously determined interactions are direct. We focused on cis-regulatory networks responsible for the Ciona mesodermal tissues by examining how the networks specify these tissues at the level of their cis-regulatory architecture. We also found many interactions that had not been predicted by simple gene knockdown experiments, and we showed that a significant fraction of TF-DNA interactions make major contributions to the regulatory control of target gene expression.

Highly divergent gene expression programs can lead to similar chordate larval body plans

The diversity of animal morphologies is thought to result largely from spatial or temporal variations in gene expression. Conversely, we explored here the extent of divergence in transcriptional expression patterns compatible with a common morphological output, the chordate larva. We compared two organisms that share a prototypical tadpole larval body plan but are separated by over half a billion years of divergent evolution: the zebrafish (Danio rerio) and the ascidian Ciona intestinalis, an invertebrate chordate belonging to the sister group of vertebrates. The large databases of whole-mount in situ hybridization expression patterns available for these two species allowed us to carry out a systematic large-scale comparison of spatiotemporal expression patterns of 1103 groups of orthologous genes. We found an extensive overall divergence in gene expression profiles between the two species that was similar at all developmental stages and did not discriminate developmental regulators from their targets. The level of conservation in individual tissues, however, varied. Conservation of tissue-specific expression patterns was highest in tissues involved in locomotion, including muscle, notochord, and the central nervous system. Thus, a broad divergence in gene expression profiles is compatible with the conservation of similar body plans across large evolutionary distances.

Differential regional expression of genes in the developing brain of Ciona intestinalis embryos

Our previous transcriptome analysis identified 565 genes that are preferentially expressed in the developing brain of Ciona intestinalis larvae. Here, we show by in-situ hybridization that the spatial expression patterns of these brain-specific genes fall into different categories depending on the regions where the gene is expressed. For example, Ci-opsin3 and Ci-Dkk3 are expressed in the entire brain, Ci-tyrosinase and Ci-TYRP1 in the dorsal region, and Ci-synaptotagmin3, Ci-ZF399, and Ci-PTFb in the ventral region. Other genes are specific to the posterior, anterior, central, posterior and ventral, or anterior-ventral region of the brain. This regional expression of genes in the Ciona brain is not always associated with cell lineage, suggesting that complex mechanisms control the regionalized expression of brain-specific genes.

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.

Gene-regulatory networks in the Ciona embryos

Ascidians belong to the subphylum Urochordata or Tunicata, which is the sister group of the vertebrates. The simple architecture of the ascidian larva represents the basic chordate body plan. Recent analyses have shown many instances of developmental mechanisms conserved during evolution, while these studies have also revealed a much larger number of instances of divergence. However, to precisely determine the degree of conservation and divergence, that is, how many ways are used to make tadpole-like larvae, we need a systems-level understanding of development. Because animal development is organized by the genome and the minimal functional unit of development is a cell, comprehensiveness and single-cell resolution are necessary for a systems-biological understanding of the development. In the ascidian Ciona intestinalis, gene-regulatory networks responsible for the embryonic development have been studied on a genome-wide scale and at single-cell resolution. The simplicity and compactness of the genome facilitates genome-wide studies. In the Ciona genome, only approximately 670 transcription factor genes are encoded, and their expression profiles during the embryonic development have been analyzed. Gene-knockdown analyses of the transcription factor genes expressed during the embryonic development have been performed. The simplicity of the embryo permits these analyses to be done at single-cell resolution. Actually, these simple embryos are now being modeled in the computer, which allows us to understand the gene-regulatory networks very precisely in three dimensions.

Evolutionary conservation of vertebrate notochord genes in the ascidian Ciona intestinalis

To reconstruct a minimum complement of notochord genes evolutionarily conserved across chordates, we scanned the Ciona intestinalis genome using the sequences of 182 genes reported to be expressed in the notochord of different vertebrates and identified 139 candidate notochord genes. For 66 of these Ciona genes expression data were already available, hence we analyzed the expression of the remaining 73 genes and found notochord expression for 20. The predicted products of the newly identified notochord genes range from the transcription factors Ci-XBPa and Ci-miER1 to extracellular matrix proteins. We examined the expression of the newly identified notochord genes in embryos ectopically expressing Ciona Brachyury (Ci-Bra) and in embryos expressing a repressor form of this transcription factor in the notochord, and we found that while a subset of the genes examined are clearly responsive to Ci-Bra, other genes are not affected by alterations in its levels. We provide a first description of notochord genes that are not evidently influenced by the ectopic expression of Ci-Bra and we propose alternative regulatory mechanisms that might control their transcription.

C2H2 zinc finger genes of the Gli, Zic, KLF, SP, Wilms' tumour, Huckebein, Snail, Ovo, Spalt, Odd, Blimp-1, Fez and related gene families from Branchiostoma floridae

The C2H2 zinc finger is one of the most common domains encoded by animal genomes and has been implicated in DNA binding as well as protein-protein interactions and RNA binding. Genes encoding C2H2 zinc finger domains include not only well-studied conserved transcription factors such as Gli and Snail but also include a large diversity of more rapidly evolving genes. Here, I focus on the description of amphioxus members of families and super-families of C2H2 zinc finger genes that have been the subject of functional studies in other species, specifically the Gli, Zic, Glis, Snail, Scratch, Krox, Wilms' tumour, Huckebein, SP, KLF, Ovo, Spalt, Blimp-1, Odd and Fez genes. Surveys of the Branchiostoma floridae genome reveal members of all of these groups of genes. Genes are named according to molecular phylogenetic analyses, such that the nomenclature reflects pre-existing gene names in the context of gene families that have descended from a single common ancestral gene in the common ancestor of chordates and insects. In total, this comprises 28 B. floridae C2H2 zinc finger genes, representing at least 15 gene families. For 17 of these genes, expressed sequence tag clusters and associated clone identification codes relating to the B. floridae gene collection are given.

Postplasmic/PEM RNAs: a class of localized maternal mRNAs with multiple roles in cell polarity and development in ascidian embryos

Ascidian is a good model to understand the cellular and molecular mechanisms responsible for mRNA localization with the discovery of a large family of localized maternal mRNAs, called postplasmic/PEM RNAs, which includes more than 40 members in three different ascidian species (Halocynthia roretzi, Ciona intestinalis, and C. savignyi). Among these mRNAs, two types (Type I and Type II) have been identified and show two different localization patterns from fertilization to the eight-cell stage. At the eight-cell stage, both types concentrate to a macromolecular cortical structure called CAB (for Centrosome Attracting Body) in the posterior-vegetal B4.1 blastomeres. The CAB is responsible for unequal cleavages and the partitioning of postplasmic/PEM RNAs at the posterior pole of embryos during cleavage stages. It has also been suggested that the CAB region could contain putative germ granules. In this review, we discuss recent data obtained on the distribution of Type I postplasmic/PEM RNAs from oogenesis to late development, in relation to their localization and translational control. We have first regrouped localization patterns for Type I and Type II into a comparative diagram and included all important definitions in the field. We also have made an exhaustive classification of their embryonic expression profiles (Type I or Type II), and analyzed their functions after knockdown and/or overexpression experiments and the role of the 3'-untranslated region (3'UTR) controlling both their localization and translation. Finally, we propose a speculative model integrating recent data, and we also discuss the relationship between postplasmic/PEM RNAs, posterior specification, and germ cell formation in ascidians.

Novel genes involved in Ciona intestinalis embryogenesis: characterization of gene knockdown embryos

The sequenced genome of the urochordate ascidian Ciona intestinalis contains nearly 2,500 genes that have vertebrate homologues, but their functions are as yet unknown. To identify novel genes involved in early chordates embryogenesis, we previously screened 200 Ciona genes by knockdown experiments using specific morpholino oligonucleotides and found that suppression of the translation of 40 genes caused embryonic defects (Yamada et al. [2003] Development 130:6485-6495). We have since examined an additional 304 genes, that is, screening 504 genes overall, and a total of 111 genes showed morphological defects when gene function was suppressed. We further examined the role of these genes in the differentiation of six major tissues of the embryo: endoderm, muscle, epidermis, neural tissue, mesenchyme, and notochord. Based on the similarity of phenotypes of gene knockdown embryos, genes were categorized into several groups, with the suggestion that the genes within a given group are involved in similar developmental processes. For example, five were shown to be novel genes that are likely involved in beta-catenin-mediated endoderm formation. The type of large-scale screening used is, therefore, a powerful approach to identify novel genes with significant developmental functions, the details of which will be determined in future studies.

An amphioxus Gli gene reveals conservation of midline patterning and the evolution of hedgehog signalling diversity in chordates

BACKGROUND

Hedgehog signalling, interpreted in receiving cells by Gli transcription factors, plays a central role in the development of vertebrate and Drosophila embryos. Many aspects of the signalling pathway are conserved between these lineages, however vertebrates have diverged in at least one key aspect: they have evolved multiple Gli genes encoding functionally-distinct proteins, increasing the complexity of the hedgehog-dependent transcriptional response. Amphioxus is one of the closest living relatives of the vertebrates, having split from the vertebrate lineage prior to the widespread gene duplication prominent in early vertebrate evolution.

PRINCIPAL FINDINGS

We show that amphioxus has a single Gli gene, which is deployed in tissues adjacent to sources of hedgehog signalling derived from the midline and anterior endoderm. This shows the duplication and divergence of the Gli gene family, and hence the origin of vertebrate Gli functional diversity, was specific to the vertebrate lineage. However we also show that the single amphioxus Gli gene produces two distinct transcripts encoding different proteins. We utilise three tests of Gli function to examine the transcription regulatory capacities of these different proteins, demonstrating one has activating activity similar to Gli2, while the other acts as a weak repressor, similar to Gli3.

CONCLUSIONS

These data show that vertebrates and amphioxus have evolved functionally-similar repertoires of Gli proteins using parallel molecular routes; vertebrates via gene duplication and divergence, and amphioxus via alternate splicing of a single gene. Our results demonstrate that similar functional complexity of intercellular signalling can be achieved via different evolutionary pathways.

The dystrotelin, dystrophin and dystrobrevin superfamily: new paralogues and old isoforms

BACKGROUND

Dystrophins and dystrobrevins are distantly related proteins with important but poorly understood roles in the function of metazoan muscular and neuronal tissues. Defects in them and their associated proteins cause a range of neuromuscular disorders. Members of this superfamily have been discovered in a relatively serendipitous way; we set out to compile a comprehensive description of dystrophin- and dystrobrevin-related sequences from available metazoan genome sequences, validated in representative organisms by RT-PCR, or acquired de novo from key species.

RESULTS

Features of the superfamily revealed by our survey include: a) Dystrotelin, an entirely novel branch of the superfamily, present in most vertebrates examined. Dystrotelin is expressed in the central nervous system, and is a possible orthologue of Drosophila DAH. We describe the preliminary characterisation of its function, evolution and expression. b) A novel vertebrate member of the dystrobrevin family, gamma-dystrobrevin, an ancient branch now extant only in fish, but probably present in our own ancestors. Like dystrophin, zebrafish gamma-dystrobrevin mRNA is localised to myosepta. c) The extent of conservation of alternative splicing and alternative promoter use in the dystrophin and dystrobrevin genes; alternative splicing of dystrophin exons 73 and 78 and alpha-dystrobrevin exon 13 are conserved across vertebrates, as are the use of the Dp116, Dp71 and G-utrophin promoters; the Dp260 and Dp140 promoters are tetrapod innovations. d) The evolution of the unique N-terminus of DRP2 and its relationship to Dp116 and G-utrophin. e) A C-terminally truncated common ancestor of dystrophin and utrophin in cyclostomes. f) A severely restricted repertoire of dystrophin complex components in ascidians.

CONCLUSION

We have refined our understanding of the evolutionary history and isoform diversity of the five previously reported vertebrate superfamily members and describe two novel members, dystrotelin and gamma-dystrobrevin. Dystrotelins, dystrophins and dystrobrevins are roughly equally related to each other. Vertebrates therefore have a repertoire of seven superfamily members (three dystrophins, three dystrobevins, and one dystrotelin), with one lost in tetrapods. Most invertebrates studied have one member from each branch. Although the basic shared function which is implied by the common architecture of these distantly related proteins remains unclear, it clearly permeates metazoan biology.

The C2H2 zinc finger genes of Strongylocentrotus purpuratus and their expression in embryonic development

The C2H2 zinc finger is one of the most abundant protein domains and is thought to have been extensively replicated in diverse animal clades. Some well-studied proteins that contain this domain are transcriptional regulators. As part of an attempt to delineate all transcription factors encoded in the Strongylocentrotus purpuratus genome, we identified the C2H2 zinc finger genes indicated in the sequence, and examined their involvement in embryonic development. We found 377 zinc finger genes in the sea urchin genome, about half the number found in mice or humans. Their expression was measured by quantitative PCR. Up to the end of gastrulation less than a third of these genes is expressed, and about 75% of the expressed genes are maternal; both parameters distinguish these from all other classes of regulatory genes as measured in other studies. Spatial expression pattern was determined by whole mount in situ hybridization for 43 genes transcribed at a sufficient level, and localized expression was observed in diverse embryonic tissues. These genes may execute important regulatory functions in development. However, the functional meaning of the majority of this large gene family remains undefined.

Embryonic expression profiles and conserved localization mechanisms of pem/postplasmic mRNAs of two species of ascidian, Ciona intestinalis and Ciona savignyi

In many animals, the first cue for development is transcripts and/or proteins that are provided maternally and are localized at specific regions of fertilized eggs and early embryos. The ascidian is known to exhibit a mosaic mode of development, which is largely dependent on localized maternal factors. In early Ciona intestinalis embryos, the posterior-most localization appears to be the major specialized pattern of maternal transcripts. The present study examined the temporal and spatial expression pattern of 40 genes known as pem/postplasmic genes, for which maternal mRNAs are localized at the posterior-most region during early Ciona embryogenesis. Ten of these genes showed redistribution to B8.12-line cells, which are known to give rise to germ cells in ascidians. In addition 23 orthologues were newly identified in a related ascidian species, Ciona savignyi, and 16 of them showed the mRNA localization pattern at the posterior-most region. Furthermore, the localized pattern of exogenous mRNA, which comprised the 3' UTR of C. intestinalis pem/postplasmic genes conjugated with the LacZ ORF, showed the localization at the posterior-most region in C. savignyi embryos. Likewise, the 3' UTR of C. savignyi pem/postplasmic genes conjugated with the LacZ ORF showed localization at the posterior most region in C. intestinalis embryos, suggesting that localization mechanisms are conserved between the two species. The present study therefore provides basic information for future functional analyses of these pem/postplasmic genes and for exploring the mechanisms of localization of mRNAs.