Regulatory blueprint for a chordate embryo

Robust target gene discovery through transcriptome perturbations and genome-wide enhancer predictions in Drosophila uncovers a regulatory basis for sensory specification

CisTarget X is a novel computational method that accurately predicts Atonal governed regulatory networks in the retina of the fruit fly.

A comprehensive systems-level understanding of developmental programs requires the mapping of the underlying gene regulatory networks. While significant progress has been made in mapping a few such networks, almost all gene regulatory networks underlying cell-fate specification remain unknown and their discovery is significantly hampered by the paucity of generalized, in vivo validated tools of target gene and functional enhancer discovery. We combined genetic transcriptome perturbations and comprehensive computational analyses to identify a large cohort of target genes of the proneural and tumor suppressor factor Atonal, which specifies the switch from undifferentiated pluripotent cells to R8 photoreceptor neurons during larval development. Extensive in vivo validations of the predicted targets for the proneural factor Atonal demonstrate a 50% success rate of bona fide targets. Furthermore we show that these enhancers are functionally conserved by cloning orthologous enhancers from Drosophila ananassae and D. virilis in D. melanogaster. Finally, to investigate cis-regulatory cross-talk between Ato and other retinal differentiation transcription factors (TFs), we performed motif analyses and independent target predictions for Eyeless, Senseless, Suppressor of Hairless, Rough, and Glass. Our analyses show that cisTargetX identifies the correct motif from a set of coexpressed genes and accurately predicts target genes of individual TFs. The validated set of novel Ato targets exhibit functional enrichment of signaling molecules and a subset is predicted to be coregulated by other TFs within the retinal gene regulatory network.

Tens of thousands of regulatory elements determine the spatiotemporal expression pattern of protein-coding genes in the metazoan genome. Each regulatory element, when bound by the appropriate transcription factors, can affect the temporal transcription of a nearby target gene in a particular cell type. Annotating the genome for regulatory elements, as well as determining the input transcription factors for each element, is a key challenge in genome biology. In this study, we introduce a computational method, cisTargetX, that predicts transcription factor binding motifs and their target genes through the integration of gene expression data and comparative genomics. We first validate this method in silico using public gene expression data and, then, apply cisTargetX to the developmental program governing photoreceptor neuron specification in the retina of Drosophila melanogaster. Particularly, we perturbed predicted key transcription factors during the initial steps of neurogenesis; measure gene expression by microarrays; identify motifs and predict target genes; validate the predictions in vivo using transgenic animals; and study several functional and evolutionary aspects of the validated regulatory elements for the proneural factor Atonal. Overall, we show that cisTargetX efficiently predicts genetic regulatory interactions and provides mechanistic insight into gene regulatory networks of postembryonic developmental systems.

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.

doublesex/mab3 related-1 (dmrt1) is essential for development of anterior neural plate derivatives in Ciona

Ascidian larvae have a hollow, dorsal central nervous system that shares many morphological features with vertebrate nervous systems yet is composed of very few cells. We show here that a null mutation in the gene dmrt1 in the ascidian Ciona savignyi results in profound abnormalities in the development of the sensory vesicle (brain), as well as other anterior ectodermal derivatives, including the palps and oral siphon primordium (OSP). Although the phenotype of the mutant embryos is variable, the majority have a complete loss of the most anterior structures (palps and OSP) and extensive disruption of sensory structures, such as the light-sensitive ocellus, in the sensory vesicle. dmrt1 is expressed early in the blastula embryo in a small group of presumptive ectodermal cells as they become restricted to anterior neural, OSP and palp fates. Despite the early and restricted expression of dmrt1, we were unable, using several independent criteria, to observe a defect in the mutant embryos until the early tailbud stage. We speculate that the variability and late onset in the phenotype may be due to partially overlapping activities of other gene products.

Macho-1 regulates unequal cell divisions independently of its function as a muscle determinant

The anterior-posterior (A-P) axis in ascidian embryos is established through the posteriorizing activities of a localized egg region known as the posterior vegetal cortex/cytoplasm (PVC). Here we describe a novel function of macho-1, a maternally-localized muscle determinant, in establishment of the A-P axis in the Halocynthia roretzi embryo. Macho-1, in addition to its known function in the formation of posterior tissue such as muscle and mesenchyme, and suppression of the anterior-derived notochord fate, acts independently of its transcriptional activity as a regulator of posterior-specific unequal cell divisions, in cooperation with beta-catenin. Our results suggest that macho-1 and beta-catenin regulate the formation of a microtubule bundle that shortens and pulls the centrosome toward a sub-cellular cortical structure known as centrosome-attracting body (CAB), which is located at the posterior pole of the embryo during unequal cell divisions, and act upstream of PEM, a recently-identified regulator of unequal cell divisions. We also present data that suggest that PEM localization to the CAB may not be required for unequal cleavage regulation. The present study provides an important and novel insight into the role of the zinc-finger-containing transcription factor and indicates that it constitutes a major part of the PVC activity.

Systems biology of embryogenesis

The development of a complete organism from a single cell involves extraordinarily complex orchestration of biological processes that vary intricately across space and time. Systems biology seeks to describe how all elements of a biological system interact in order to understand, model and ultimately predict aspects of emergent biological processes. Embryogenesis represents an extraordinary opportunity (and challenge) for the application of systems biology. Systems approaches have already been used successfully to study various aspects of development, from complex intracellular networks to four-dimensional models of organogenesis. Going forward, great advancements and discoveries can be expected from systems approaches applied to embryogenesis and developmental biology.

Regulation of notochord-specific expression of Ci-Bra downstream genes in Ciona intestinalis embryos

Brachyury, a T-box transcription factor, is expressed in ascidian embryos exclusively in primordial notochord cells and plays a pivotal role in differentiation of notochord cells. Previously, we identified approximately 450 genes downstream of Ciona intestinalis Brachyury (Ci-Bra), and characterized the expression profiles of 45 of these in differentiating notochord cells. In this study, we looked for cisregulatory sequences in minimal enhancers of 20 Ci-Bra downstream genes by electroporating region within approximately 3 kb upstream of each gene fused with lacZ. Eight of the 20 reporters were expressed in notochord cells. The minimal enchancer for each of these eight genes was narrowed to a region approximately 0.5-1.0-kb long. We also explored the genome-wide and coordinate regulation of 43 Ci-Bra-downstream genes. When we determined their chromosomal localization, it became evident that they are not clustered in a given region of the genome, but rather distributed evenly over 13 of the 14 pairs of chromosomes, suggesting that gene clustering does not contribute to coordinate control of the Ci-Bra downstream gene expression. Our results might provide Insights Into the molecular mechanisms underlying notochord formation in chordates.

Neuropeptides, hormone peptides, and their receptors in Ciona intestinalis: an update

The critical phylogenetic position of ascidians leads to the presumption that neuropeptides and hormones in vertebrates are highly likely to be evolutionarily conserved in ascidians, and the cosmopolitan species Ciona intestinalis is expected to be an excellent deuterostome Invertebrate model for studies on neuropeptides and hormones. Nevertheless, molecular and functional characterization of Ciona neuropeptides and hormone peptides was restricted to a few peptides such as a cholecystokinin (CCK)/gastrin peptide, cionin, and gonadotropin-releasing hormones (GnRHs). In the past few years, mass spectrometric analyses and database searches have detected Ciona orthologs or prototypes of vertebrate peptides and their receptors, including tachykinin, insulin/relaxin, calcitonin, and vasopressin. Furthermore, studies have shown that several Ciona peptides, including vasopressin and a novel GnRH-related peptide, have acquired ascidian-specific molecular forms and/or biological functions. These findings provided indisputable evidence that ascidians, unlike other invertebrates (including the traditional protostome model animals), possess neuropeptides and hormone peptides structurally and functionally related to vertebrate counterparts, and that several peptides have uniquely diverged in ascidian evolutionary lineages. Moreover, recent functional analyses of Ciona tachykinin in the ovary substantiated the novel tachykininergic protease-assoclated oocyte growth pathway, which could not have been revealed in studies on vertebrates. These findings confirm the outstanding advantages of ascidians in understanding the neuroscience, endocrinology, and evolution of vertebrate neuropeptides and hormone peptides. This article provides an overview of basic findings and reviews new knowledge on ascidian neuropeptides and hormone peptides.

Identification of genes downstream of nodal in the Ciona intestinalis embryo

Nodal, a growth factor belonging to the TGF-beta superfamily, is required for the formation of the neural tube in Ciona intestinalis. Previous studies have revealed many genes whose expression is controlled by Nodal in the Ciona embryo; however, all of them encode transcription factors and signaling molecules. In the present study, we identified five genes upregulated or downregulated by the overexpression of Nodal in embryos of C. intestinalis. The upregulated genes included those encoding type IV collagen 1/3/5, laminin-alpha5, and Prickle. The downregulated genes included those encoding glypican and delta1-protocadherln-like. Many of these genes were expressed in the neural plate at the late gastrula stage. The present study revealed candidate effector genes that directly regulate, in response to Nodal, the morphogenesis of the neural tube in Ciona intestinalis.

Computational identification of Ciona intestinalis microRNAs

MicroRNAs (miRNAs) are conserved non-coding small RNAs with potent post-transcriptional gene regulatory functions. Recent computational approaches and sequencing of small RNAs had indicated the existence of about 80 Ciona intestinalis miRNAs, although it was not clear whether other miRNA genes were present in the genome. We undertook an alternative computational approach to look for Ciona miRNAs. Conserved non-coding sequences from the C. intestinalis genome were extracted and computationally folded to identify putative hairpin-like structures. After applying additional criteria, we obtained 458 miRNA candidates whose sequences were used to design a custom microarray. Over 100 of our predicted hairpins were identified in this array when probed with RNA from various Ciona stages. We also compared our predictions to recently deposited sequences of Ciona small RNAs and report that 170 of our predicted hairpins are represented in this data set. Altogether, about 250 of our 458 predicted miRNAs were represented in either our array data or the small-RNA sequence database. These results suggest that Ciona has a large number of genomically encoded miRNAs that play an important role in modulating gene activity in developing embryos and adults.

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.

High-throughput sequence analysis of Ciona intestinalis SL trans-spliced mRNAs: alternative expression modes and gene function correlates

Pre-mRNA 5' spliced-leader (SL) trans-splicing occurs in some metazoan groups but not in others. Genome-wide characterization of the trans-spliced mRNA subpopulation has not yet been reported for any metazoan. We carried out a high-throughput analysis of the SL trans-spliced mRNA population of the ascidian tunicate Ciona intestinalis by 454 Life Sciences (Roche) pyrosequencing of SL-PCR-amplified random-primed reverse transcripts of tailbud embryo RNA. We obtained approximately 250,000 high-quality reads corresponding to 8790 genes, approximately 58% of the Ciona total gene number. The great depth of this data revealed new aspects of trans-splicing, including the existence of a significant class of "infrequently trans-spliced" genes, accounting for approximately 28% of represented genes, that generate largely non-trans-spliced mRNAs, but also produce trans-spliced mRNAs, in part through alternative promoter use. Thus, the conventional qualitative dichotomy of trans-spliced versus non-trans-spliced genes should be supplanted by a more accurate quantitative view recognizing frequently and infrequently trans-spliced gene categories. Our data include reads representing approximately 80% of Ciona frequently trans-spliced genes. Our analysis also revealed significant use of closely spaced alternative trans-splice acceptor sites which further underscores the mechanistic similarity of cis- and trans-splicing and indicates that the prevalence of +/-3-nt alternative splicing events at tandem acceptor sites, NAGNAG, is driven by spliceosomal mechanisms, and not nonsense-mediated decay, or selection at the protein level. The breadth of gene representation data enabled us to find new correlations between trans-splicing status and gene function, namely the overrepresentation in the frequently trans-spliced gene class of genes associated with plasma/endomembrane system, Ca(2+) homeostasis, and actin cytoskeleton.

The molecular machinery of germ line specification

Germ cells occupy a unique position in animal reproduction, development, and evolution. In sexually reproducing animals, only they can produce gametes and contribute genetically to subsequent generations. Nonetheless, germ line specification during embryogenesis is conceptually the same as the specification of any somatic cell type: germ cells must activate a specific gene regulatory network in order to differentiate and go through gametogenesis. While many genes with critical roles in the germ line have been characterized with respect to expression pattern and genetic interactions, it is the molecular interactions of the relevant gene products that are ultimately responsible for germ cell differentiation. This review summarizes the current state of knowledge on the molecular functions and biochemical connections between germ line gene products. We find that homologous genes often interact physically with the same conserved molecular partners across the metazoans. We also point out cases of nonhomologous genes from different species whose gene products play analogous biological roles in the germ line. We suggest a preliminary molecular definition of an ancestral "pluripotency module" that could have been modified during metazoan evolution to become specific to the germ line.

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.

The Motif Tool Assessment Platform (MTAP) for sequence-based transcription factor binding site prediction tools

Predicting transcription factor binding sites (TFBS) from sequence is one of the most challenging problems in computational biology. The development of (semi-)automated computer-assisted prediction methods is needed to find TFBS over an entire genome, which is a first step in reconstructing mechanisms that control gene activity. Bioinformatics journals continue to publish diverse methods for predicting TFBS on a monthly basis. To help practitioners in deciding which method to use to predict for a particular TFBS, we provide a platform to assess the quality and applicability of the available methods. Assessment tools allow researchers to determine how methods can be expected to perform on specific organisms or on specific transcription factor families. This chapter introduces the TFBS detection problem and reviews current strategies for evaluating algorithm effectiveness. In this chapter, a novel and robust assessment tool, the Motif Tool Assessment Platform (MTAP), is introduced and discussed.

Practical computational methods for regulatory genomics: a cisGRN-Lexicon and cisGRN-browser for gene regulatory networks

The CYRENE Project focuses on the study of cis-regulatory genomics and gene regulatory networks (GRN) and has three components: a cisGRN-Lexicon, a cisGRN-Browser, and the Virtual Sea Urchin software system. The project has been done in collaboration with Eric Davidson and is deeply inspired by his experimental work in genomic regulatory systems and gene regulatory networks. The current CYRENE cisGRN-Lexicon contains the regulatory architecture of 200 transcription factors encoding genes and 100 other regulatory genes in eight species: human, mouse, fruit fly, sea urchin, nematode, rat, chicken, and zebrafish, with higher priority on the first five species. The only regulatory genes included in the cisGRN-Lexicon (CYRENE genes) are those whose regulatory architecture is validated by what we call the Davidson Criterion: they contain functionally authenticated sites by site-specific mutagenesis, conducted in vivo, and followed by gene transfer and functional test. This is recognized as the most stringent experimental validation criterion to date for such a genomic regulatory architecture. The CYRENE cisGRN-Browser is a full genome browser tailored for cis-regulatory annotation and investigation. It began as a branch of the Celera Genome Browser (available as open source at http://sourceforge.net/projects/celeragb /) and has been transformed to a genome browser fully devoted to regulatory genomics. Its access paradigm for genomic data is zoom-to-the-DNA-base in real time. A more recent component of the CYRENE project is the Virtual Sea Urchin system (VSU), an interactive visualization tool that provides a four-dimensional (spatial and temporal) map of the gene regulatory networks of the sea urchin embryo.

Mechanism of DNA replication-dependent transcriptional activation of the acetylcholinesterase gene in the Ciona intestinalis embryo

The acetylcholinesterase-encoding gene in the ascidian Ciona intestinalis (Ci-AChE) is expressed in tail muscle cells from the gastrula stage. When the embryo was continuously treated with aphidicolin from the 32-cell stage, Ci-AChE was not expressed even when control embryos reached the tailbud stage. This result suggests that Ci-AChE acquires the competence to be transcribed after passing through a certain number of DNA replication cycles. A lacZ reporter gene containing the 5' flanking region of Ci-AChE was expressed in the tail muscle cells. Aphidicolin treatment from the 32-cell stage affected, but did not completely suppress, the expression of lacZ. A bisulfite sequencing analysis was carried out to examine the methylation status of four regions within the 5' flanking sequence and the first exon. However, all of these regions remained unmethylated from the 16-cell to 110-cell stages. The results suggested that the DNA of the Ci-AChE locus is not responsible for counting the rounds of replication. We examined the expression of the C. intestinalis MyoD (Ci-MyoD), a transcription factor that activates Ci-AChE. Aphidicolin treatment from the 32-cell stage suppressed the expression of Ci-MyoD, even when control embryos reached the gastrula stage. These results suggest that a lack of Ci-MyoD is critical to the suppression of Ci-AChE in aphidicolin-treated embryos.

The sea squirt Ciona intestinalis

Sea squirts (Ciona intestinalis) are tunicates (or urochordates), the closest living relatives of the vertebrates. Although the adults are simple, sessile filter feeders, the embryos and larvae possess clear chordate features including a prominent notochord and dorsal, hollow neural tube. Tail-bud-stage embryos and mature swimming tadpoles are composed of approximately 1000 and 2600 cells, respectively, with complete lineage information. This cellular simplicity is coupled with a streamlined genome that has not undergone the duplications seen in vertebrates. A variety of molecular tools have been applied to understanding Ciona embryogenesis. Comparisons of the C. intestinalis genome and the related but divergent Ciona savignyi genome have facilitated the identification of conserved non-coding DNAs, including regulatory DNAs such as tissue-specific enhancers. Systematic in situ hybridization assays and gene-disruption experiments using specific morpholino antisense oligonucleotides have led to the elaboration of provisional gene regulatory networks underlying the specification of key chordate tissues, including the notochord, neural tube, and beating heart. These networks provide a foundation for understanding the mechanistic basis of more complex cell-specification processes in vertebrates, and for understanding the evolutionary origins of distinctive vertebrate characteristics such as the neural crest. Because tunicates and vertebrates are sister groups, there is every indication that the developmental mechanisms revealed in the simple Ciona model will be applicable to comparable processes in vertebrates.

Direct activation of a notochord cis-regulatory module by Brachyury and FoxA in the ascidian Ciona intestinalis

The notochord is a defining feature of the chordate body plan. Experiments in ascidian, frog and mouse embryos have shown that co-expression of Brachyury and FoxA class transcription factors is required for notochord development. However, studies on the cis-regulatory sequences mediating the synergistic effects of these transcription factors are complicated by the limited knowledge of notochord genes and cis-regulatory modules (CRMs) that are directly targeted by both. We have identified an easily testable model for such investigations in a 155-bp notochord-specific CRM from the ascidian Ciona intestinalis. This CRM contains functional binding sites for both Ciona Brachyury (Ci-Bra) and FoxA (Ci-FoxA-a). By combining point mutation analysis and misexpression experiments, we demonstrate that binding of both transcription factors to this CRM is necessary and sufficient to activate transcription. To gain insights into the cis-regulatory criteria controlling its activity, we investigated the organization of the transcription factor binding sites within the 155-bp CRM. The 155-bp sequence contains two Ci-Bra binding sites with identical core sequences but opposite orientations, only one of which is required for enhancer activity. Changes in both orientation and spacing of these sites substantially affect the activity of the CRM, as clusters of identical sites found in the Ciona genome with different arrangements are unable to activate transcription in notochord cells. This work presents the first evidence of a synergistic interaction between Brachyury and FoxA in the activation of an individual notochord CRM, and highlights the importance of transcription factor binding site arrangement for its function.

Transcriptional regulation of the retinoic acid receptor in the dorsal midline epidermis in the Ciona intestinalis embryo

Retinoic acid regulates the spatial pattern of gene expression mainly in the epidermis in the protochordate ascidian, Ciona intestinalis. Our previous study characterized the enhancer element responsible for the activation of the retinoic acid receptor (RAR) in the dorsal midline epidermis (DME). In the present study, detailed analysis identified two small sequence elements required for the transcriptional activation in the DME. Deletion of either one of these sequences resulted in suppression of the DME enhancer, suggesting that both sites are necessary. The nucleotide sequences of these two sites were similar to the consensus recognition sequences for the Msx and Sox transcription factors, respectively. These transcription factors are expressed in the DME lineage blastomeres (b7.9 and b7.10 blastomere pairs of bilaterally symmetrical embryos of this species) at the 64-cell stage. Electrophoretic mobility shift assays revealed that recombinant Msxb and SoxB1 proteins specifically bound to these sequences. These results suggest that Msxb and SoxB1 synergistically activate the enhancer in the DME.

The repertoire of heterotrimeric G proteins and RGS proteins in Ciona intestinalis

Background

Heterotrimeric G proteins and regulators of G protein signaling (RGS) proteins are key downstream interacting partners in the G protein coupled receptor (GPCR) signaling pathway. The highly versatile GPCR transmembrane signaling system is a consequence of the coupling of a diverse set of receptors to downstream partners that include multiple subforms of G proteins and regulatory proteins including RGS proteins, among others. While the GPCR repertoire of Ciona intestinalis, representing the basal chordate is known, the repertoire of the heterotrimeric G proteins and RGS proteins is unknown.

Methodology/Principal Findings

In the present study, we performed an in-silico genome-wide search of C. intestinalis for its complement of G proteins and RGS proteins. The identification of several one-to-one orthologs of human G proteins at the levels of families, subfamilies and types and of homologs of the human RGS proteins suggests an evolutionarily conserved structure function relationship of the GPCR signaling mechanism in the chordates.

Conclusions

The C. intestinalis genome encodes a highly conserved, albeit, limited repertoire of the heterotrimeric G protein complexes with the size of subunit types comparable with that in lower eukaryotes.

Early lineage specification of long-lived germline precursors in the colonial ascidian Botryllus schlosseri

In many taxa, germline precursors segregate from somatic lineages during embryonic development and are irreversibly committed to gametogenesis. However, in animals that can propagate asexually, germline precursors can originate in adults. Botryllus schlosseri is a colonial ascidian that grows by asexual reproduction, and on a weekly basis regenerates all somatic and germline tissues. Embryonic development in solitary ascidians is the classic example of determinative specification, and we are interested in both the origins and the persistence of stem cells responsible for asexual development in colonial ascidians. In this study, we characterized vasa as a putative marker of germline precursors. We found that maternally deposited vasa mRNA segregates early in development to a posterior lineage of cells, suggesting that germline formation is determinative in colonial ascidians. In adults, vasa expression was observed in the gonads, as well as in a population of mobile cells scattered throughout the open circulatory system, consistent with previous transplantation/reconstitution results. vasa expression was dynamic during asexual development in both fertile and infertile adults, and was also enriched in a population of stem cells. Germline precursors in juveniles could contribute to gamete formation immediately upon transplantation into fertile adults, thus vasa expression is correlated with the potential for gamete formation, which suggests that it is a marker for embryonically specified, long-lived germline progenitors. Transient vasa knockdown did not have obvious effects on germline or somatic development in adult colonies, although it did result in a profound heterochrony, suggesting that vasa might play a homeostatic role in asexual development.

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.

Global networks of functional coupling in eukaryotes from comprehensive data integration

No single experimental method can discover all connections in the interactome. A computational approach can help by integrating data from multiple, often unrelated, proteomics and genomics pipelines. Reconstructing global networks of functional coupling (FC) faces the challenges of scale and heterogeneity--how to efficiently integrate huge amounts of diverse data from multiple organisms, yet ensuring high accuracy. We developed FunCoup, an optimized Bayesian framework, to resolve these issues. Because interactomes comprise functional coupling of many types, FunCoup annotates network edges with confidence scores in support of different kinds of interactions: physical interaction, protein complex member, metabolic, or signaling link. This capability boosted overall accuracy. On the whole, the constructed framework was comprehensively tested to optimize the overall confidence and ensure seamless, automated incorporation of new data sets of heterogeneous types. Using over 50 data sets in seven organisms and extensively transferring information between orthologs, FunCoup predicted global networks in eight eukaryotes. For the Ciona intestinalis network, only orthologous information was used, and it recovered a significant number of experimental facts. FunCoup predictions were validated on independent cancer mutation data. We show how FunCoup can be used for discovering candidate members of the Parkinson and Alzheimer pathways. Cross-species pathway conservation analysis provided further support to these observations.

A gene regulatory network directed by zebrafish No tail accounts for its roles in mesoderm formation

Using chromatin immunoprecipitation combined with genomic microarrays we have identified targets of No tail (Ntl), a zebrafish Brachyury ortholog that plays a central role in mesoderm formation. We show that Ntl regulates a downstream network of other transcription factors and identify an in vivo Ntl binding site that resembles the consensus T-box binding site (TBS) previously identified by in vitro studies. We show that the notochord-expressed gene floating head (flh) is a direct transcriptional target of Ntl and that a combination of TBSs in the flh upstream region are required for Ntl-directed expression. Using our genome-scale data we have assembled a preliminary gene regulatory network that begins to describe mesoderm formation and patterning in the early zebrafish embryo.

Similar regulatory logic in Ciona intestinalis for two Wnt pathway modulators, ROR and SFRP-1/5

Anteroposterior patterning of the ectoderm in the invertebrate chordate Ciona intestinalis first relies on key zygotic activators, such as FoxA, and later on the coordinated responses to inducing signals from the underlying mesendoderm. Here, we focus on a mechanism of coordination of these responses by looking at the cis-regulatory logics of Ci-Rora and Ci-Rorb, which are coding for putative non-canonical transmembrane Wnt receptors, and are present in tandem along the C. intestinalis chromosome 08q. We showed that during cleavage stages, both Ci-Rora and Ci-Rorb genes are initially expressed in all blastomeres of the anterior ectoderm (a-line), as sFRP1/5 (Lamy, C., Rothbächer, U., Caillol, D., Lemaire, P., 2006. Ci-FoxA-a is the earliest zygotic determinant of the ascidian anterior ectoderm and directly activates Ci-sFRP1/5. Development 133, 2835-2844.). We then carried out a functional analysis of cis-regulatory regions and showed that both genes have elements enriched in Ci-FoxA-a binding sites. We dissected one of these early enhancers, and showed that it is directly activated by Ci-FoxA-a, as one sFRP1/5 cis-regulatory element. We also showed that although FoxA binding sites are abundant in genomes, dense clusters of these sites are found upstream from very few genes, and have a high predictive value of a-line expression. These data indicate an important role for FoxA in anterior specification, via the transcriptional regulation of target genes belonging to various signalling pathways.

Spatio-temporal intersection of Lhx3 and Tbx6 defines the cardiac field through synergistic activation of Mesp

Mesp encodes a bHLH transcription factor required for specification of the cardiac mesoderm in Ciona embryos. The activities of Macho-1 and beta-catenin, two essential maternal determinants, are required for Mesp expression in the B7.5 blastomeres, which constitute the heart field. The T-box transcription factor Tbx6 functions downstream of Macho-1 as a direct activator of Mesp expression. However, Tbx6 cannot account for the restricted expression of Mesp in the B7.5 lineage since it is expressed throughout the presumptive tail muscles. Here we present evidence that the LIM-homeobox gene Lhx3, a direct target of beta-catenin, is essential for localized Mesp expression. Lhx3 is expressed throughout the presumptive endoderm and B7.5 blastomeres. Thus, the B7.5 blastomeres are the only cells to express sustained levels of the Tbx6 and Lhx3 activators. Like mammalian Lhx3 genes, Ci-Lhx3 encodes two isoforms with distinct N-terminal peptides. The Lhx3a isoform appears to be expressed both maternally and zygotically, while the Lhx3b isoform is exclusively zygotic. Misexpression of Lhx3b is sufficient to induce ectopic Mesp activation in cells expressing Tbx6b. Injection of antisense morpholino oligonucleotides showed that the Lhx3b isoform is required for endogenous Mesp expression. Mutations in the Lhx3 half-site of Tbx6/Lhx3 composite elements strongly reduced the activity of a minimal Mesp enhancer. We discuss the delineation of the heart field by the synergistic action of muscle and gut determinants.

Evolution of gene regulatory network architectures: examples of subcircuit conservation and plasticity between classes of echinoderms

Developmental gene regulatory networks (GRNs) explain how regulatory states are established in particular cells during development and how these states then determine the final form of the embryo. Evolutionary changes to the sequence of the genome will direct reorganization of GRN architectures, which in turn will lead to the alteration of developmental programs. A comparison of GRN architectures must consequently reveal the molecular basis for the evolution of developmental programs among different organisms. This review highlights some of the important findings that have emerged from the most extensive direct comparison of GRN architectures to date. Comparison of the orthologous GRNs for endomesodermal specification in the sea urchin and sea star, provides examples of several discrete, functional GRN subcircuits and shows that they are subject to diverse selective pressures. This demonstrates that different regulatory linkages may be more or less amenable to evolutionary change. One of the more surprising findings from this comparison is that GRN-level functions may be maintained while the factors performing the functions have changed, suggesting that GRNs have a high capacity for compensatory changes involving transcription factor binding to cis regulatory modules.

Evolution of the neural crest viewed from a gene regulatory perspective

Neural crest cells are a vertebrate innovation and form a wide variety of embryonic cell types as diverse as peripheral neurons and facial skeleton. They undergo complex migration and differentiation processes from their site of origin in the developing central nervous system to their final destinations in the periphery. In this review, we summarize recent data on the current formulation of a gene regulatory network underlying neural crest formation and its roots at the base of the vertebrate lineage. Analyzing neural crest formation from a gene regulatory viewpoint provides insights into both the developmental mechanisms and evolutionary origins of this vertebrate-specific cell type.

Interaction of notochord-derived fibrinogen-like protein with Notch regulates the patterning of the central nervous system of Ciona intestinalis embryos

The midline organ the notochord and its overlying dorsal neural tube are the most prominent features of the chordate body plan. Although the molecular mechanisms involved in the formation of the central nervous system (CNS) have been studied extensively in vertebrate embryos, none of the genes that are expressed exclusively in notochord cells has been shown to function in this process. Here, we report a gene in the urochordate Ciona intestinalis encoding a fibrinogen-like protein that plays a pivotal role in the notochord-dependent positioning of neuronal cells. While this gene (Ci-fibrn) is expressed exclusively in notochord cells, its protein product is not confined to these cells but is distributed underneath the CNS as fibril-like protrusions. We demonstrated that Ci-fibrn interacts physically and functionally with Ci-Notch that is expressed in the central nervous system, and that the correct distribution of Ci-fibrn protein is dependent on Notch signaling. Disturbance of the Ci-fibrn distribution caused an abnormal positioning of neuronal cells and an abnormal track of axon extension. Therefore, it is highly likely that the interaction between the notochord-based fibrinogen-like protein and the neural tube-based Notch signaling plays an essential role in the proper patterning of CNS.

Brachyury null mutant-induced defects in juvenile ascidian endodermal organs

We report the isolation of a recessive ENU-induced short-tailed mutant in the ascidian Ciona intestinalis that is the product of a premature stop in the brachyury gene. Notochord differentiation and morphogenesis are severely disrupted in the mutant line. At the larval stage, variable degrees of ectopic endoderm staining were observed in the homozygous mutants, indicating that loss of brachyury results in stochastic fate transformation. In post-metamorphosis mutants, a uniform defect in tail resorption was observed, together with variable defects in digestive tract development. Some cells misdirected from the notochord lineage were found to be incorporated into definitive endodermal structures, such as stomach and intestine.

A novel and critical role for Oct4 as a regulator of the maternal-embryonic transition

Background

Compared to the emerging embryonic stem cell (ESC) gene network, little is known about the dynamic gene network that directs reprogramming in the early embryo. We hypothesized that Oct4, an ESC pluripotency regulator that is also highly expressed at the 1- to 2-cell stages in embryos, may be a critical regulator of the earliest gene network in the embryo.

Methodology/Principal Findings

Using antisense morpholino oligonucleotide (MO)-mediated gene knockdown, we show that Oct4 is required for development prior to the blastocyst stage. Specifically, Oct4 has a novel and critical role in regulating genes that encode transcriptional and post-transcriptional regulators as early as the 2-cell stage. Our data suggest that the key function of Oct4 may be to switch the developmental program from one that is predominantly regulated by post-transcriptional control to one that depends on the transcriptional network. Further, we propose to rank candidate genes quantitatively based on the inter-embryo variation in their differential expression in response to Oct4 knockdown. Of over 30 genes analyzed according to this proposed paradigm, Rest and Mta2, both of which have established pluripotency functions in ESCs, were found to be the most tightly regulated by Oct4 at the 2-cell stage.

Conclusions/Significance

We show that the Oct4-regulated gene set at the 1- to 2-cell stages of early embryo development is large and distinct from its established network in ESCs. Further, our experimental approach can be applied to dissect the gene regulatory network of Oct4 and other pluripotency regulators to deconstruct the dynamic developmental program in the early embryo.

Gene regulatory networks underlying the compartmentalization of the Ciona central nervous system

The tripartite organization of the central nervous system (CNS) may be an ancient character of the bilaterians. However, the elaboration of the more complex vertebrate brain depends on the midbrain-hindbrain boundary (MHB) organizer, which is absent in invertebrates such as Drosophila. The Fgf8 signaling molecule expressed in the MHB organizer plays a key role in delineating separate mesencephalon and metencephalon compartments in the vertebrate CNS. Here, we present evidence that an Fgf8 ortholog establishes sequential patterns of regulatory gene expression in the developing posterior sensory vesicle, and the interleaved ;neck' region located between the sensory vesicle and visceral ganglion of the simple chordate Ciona intestinalis. The detailed characterization of gene networks in the developing CNS led to new insights into the mechanisms by which Fgf8/17/18 patterns the chordate brain. The precise positioning of this Fgf signaling activity depends on an unusual AND/OR network motif that regulates Snail, which encodes a threshold repressor of Fgf8 expression. Nodal is sufficient to activate low levels of the Snail repressor within the neural plate, while the combination of Nodal and Neurogenin produces high levels of Snail in neighboring domains of the CNS. The loss of Fgf8 patterning activity results in the transformation of hindbrain structures into an expanded mesencephalon in both ascidians and vertebrates, suggesting that the primitive MHB-like activity predates the vertebrate CNS.

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.

Microarray analysis of Foxa2 mutant mouse embryos reveals novel gene expression and inductive roles for the gastrula organizer and its derivatives

Background

The Spemann/Mangold organizer is a transient tissue critical for patterning the gastrula stage vertebrate embryo and formation of the three germ layers. Despite its important role during development, there are still relatively few genes with specific expression in the organizer and its derivatives. Foxa2 is a forkhead transcription factor that is absolutely required for formation of the mammalian equivalent of the organizer, the node, the axial mesoderm and the definitive endoderm (DE). However, the targets of Foxa2 during embryogenesis, and the molecular impact of organizer loss on the gastrula embryo, have not been well defined.

Results

To identify genes specific to the Spemann/Mangold organizer, we performed a microarray-based screen that compared wild-type and Foxa2 mutant embryos at late gastrulation stage (E7.5). We could detect genes that were consistently down-regulated in replicate pools of mutant embryos versus wild-type, and these included a number of known node and DE markers. We selected 314 genes without previously published data at E7.5 and screened for expression by whole mount in situ hybridization. We identified 10 novel expression patterns in the node and 5 in the definitive endoderm. We also found significant reduction of markers expressed in secondary tissues that require interaction with the organizer and its derivatives, such as cardiac mesoderm, vasculature, primitive streak, and anterior neuroectoderm.

Conclusion

The genes identified in this screen represent novel Spemann/Mangold organizer genes as well as potential Foxa2 targets. Further investigation will be needed to define these genes as novel developmental regulatory factors involved in organizer formation and function. We have placed these genes in a Foxa2-dependent genetic regulatory network and we hypothesize how Foxa2 may regulate a molecular program of Spemann/Mangold organizer development. We have also shown how early loss of the organizer and its inductive properties in an otherwise normal embryo, impacts on the molecular profile of surrounding tissues.

A protocol describing the principles of cis-regulatory analysis in the sea urchin

cis-Regulatory analysis (CRA) is the precise identification of the cis-acting genomic sequences regulating gene transcription. As such, CRA provides essential mechanistic insight into key biological processes such as development. The first phase of this protocol involves identification of a large (approximately 100 kb) clone of genomic material surrounding the gene of interest and use of this clone to establish a reliable and unambiguous reporter assay. In the second phase, phylogenetic footprinting is used to identify candidate regulatory modules; these genomic sequences are then recursively tested for reporter activity. In the final phase, potential transcription factor binding sites are identified and disrupted in reporter constructs for individual testing. The strengths of this method reflect the use of large clones containing all relevant genomic regulatory sequences to establish a reporter assay with high fidelity. Given these foundational elements, predicted or suspected regulatory inputs can be rigorously tested and novel regulatory inputs identified. Although the expected time line varies greatly with the depth of information required, results may be obtained in as little as 4-6 months, but more detailed analyses will require several years to complete.

Mapping the gene expression universe

Methods to globally survey gene expression provide valuable insights into gene function during development. In particular, comprehensive in situ hybridization studies have demonstrated that gene expression patterns are extraordinarily diverse and new imaging techniques have been introduced to capture these patterns with higher resolution at the tissue, cellular, and subcellular levels. The analysis of massive image databases can be greatly facilitated by computer vision techniques once annotated image sets reach the crucial mass sufficient to train the computer in pattern recognition. Ultimately, genome-wide atlases of gene expression during development will record gene activity in living animals with at least cellular resolution and in the context of morphogenetic events. These emerging datasets will lead to great advances in the field of comparative genomics and revolutionize our ability to decipher and model developmental processes for a variety of organisms.

Trunk lateral cells are neural crest-like cells in the ascidian Ciona intestinalis: insights into the ancestry and evolution of the neural crest

Neural crest-like cells (NCLC) that express the HNK-1 antigen and form body pigment cells were previously identified in diverse ascidian species. Here we investigate the embryonic origin, migratory activity, and neural crest related gene expression patterns of NCLC in the ascidian Ciona intestinalis. HNK-1 expression first appeared at about the time of larval hatching in dorsal cells of the posterior trunk. In swimming tadpoles, HNK-1 positive cells began to migrate, and after metamorphosis they were localized in the oral and atrial siphons, branchial gill slits, endostyle, and gut. Cleavage arrest experiments showed that NCLC are derived from the A7.6 cells, the precursors of trunk lateral cells (TLC), one of the three types of migratory mesenchymal cells in ascidian embryos. In cleavage arrested embryos, HNK-1 positive TLC were present on the lateral margins of the neural plate and later became localized adjacent to the posterior sensory vesicle, a staging zone for their migration after larval hatching. The Ciona orthologues of seven of sixteen genes that function in the vertebrate neural crest gene regulatory network are expressed in the A7.6/TLC lineage. The vertebrate counterparts of these genes function downstream of neural plate border specification in the regulatory network leading to neural crest development. The results suggest that NCLC and neural crest cells may be homologous cell types originating in the common ancestor of tunicates and vertebrates and support the possibility that a putative regulatory network governing NCLC development was co-opted to produce neural crest cells during vertebrate evolution.