Large scale EST analyses in Ciona intestinalis: its application as Northern blot analyses

Developmental Roles and Evolutionary Significance of AMPA-Type Glutamate Receptors

Organogenesis and metamorphosis require the intricate orchestration of multiple types of cellular interactions and signaling pathways. Glutamate (Glu) is an excitatory extracellular signaling molecule in the nervous system, while Ca²⁺ is a major intracellular signaling molecule. The first Glu receptors to be cloned are Ca²⁺ -permeable receptors in mammalian brains. Although recent studies have focused on Glu signaling in synaptic mechanisms of the mammalian central nervous system, it is unclear how this signaling functions in development. Our recent article demonstrated that Ca²⁺ -permeable AMPA-type Glu receptors (GluAs) are essential for formation of a photosensitive organ, development of some neurons, and metamorphosis, including tail absorption and body axis rotation, in ascidian embryos. Based on findings in these embryos and mammalian brains, we formed several hypotheses regarding the evolution of GluAs, the non-synaptic function of Glu, the origin of GluA-positive neurons, and the neuronal network that controls metamorphosis in ascidians.

Evolutionary History of Voltage-Gated Sodium Channels

Every cell within living organisms actively maintains an intracellular Na⁺ concentration that is 10-12 times lower than the extracellular concentration. The cells then utilize this transmembrane Na⁺ concentration gradient as a driving force to produce electrical signals, sometimes in the form of action potentials. The protein family comprising voltage-gated sodium channels (Na_Vs) is essential for such signaling and enables cells to change their status in a regenerative manner and to rapidly communicate with one another. Na_Vs were first predicted in squid and were later identified through molecular biology in the electric eel. Since then, these proteins have been discovered in organisms ranging from bacteria to humans. Recent research has succeeded in decoding the amino acid sequences of a wide variety of Na_V family members, as well as the three-dimensional structures of some. These studies and others have uncovered several of the major steps in the functional and structural transition of Na_V proteins that has occurred along the course of the evolutionary history of organisms. Here we present an overview of the molecular evolutionary innovations that established present-day Na_V α subunits and discuss their contribution to the evolutionary changes in animal bodies.

A cDNA resource for gene expression studies of a hemichordate, Ptychodera flava

Recent investigations into the evolution of deuterostomes and the origin of chordates have paid considerable attention to hemichordates (acorn worms), as hemichordates and echinoderms are the closest chordate relatives. The present study prepared cDNA libraries from Ptychodera flava, to study expression and function of genes involved in development of the hemichordate body plan. Expressed sequence tag (EST) analyses of nine cDNA libraries yielded 18,832 cloned genes expressed in eggs, 18,739 in blastulae, 18,539 in gastrulae, 18,811 in larvae, 18,978 in juveniles, 11,802 in adult proboscis, 17,259 in stomochord, 11,886 in gills, and 11,580 in liver, respectively. A set of 34,159 uni-gene clones of P. flava was obtained. This cDNA resource will be valuable for studying temporal and spatial expression of acorn worm genes during development.

Expression of germline markers in three species of amphioxus supports a preformation mechanism of germ cell development in cephalochordates

Background

In a previous study, we showed that the cephalochordate amphioxus Branchiostoma floridae has localized maternal transcripts of conserved germ cell markers Vasa and Nanos in its early embryos. These results provided strong evidence to support a preformation mechanism for primordial germ cell (PGC) development in B. floridae.

Results

In this study, we further characterize the expression of B. floridae homologs of Piwi and Tudor, which play important roles in germline development in diverse metazoan animals. We show that maternal mRNA of one of the identified Piwi-like homologs, Bf-Piwil1, also colocalizes with Vasa in the vegetal germ plasm and has zygotic expression in both the putative PGCs and the tail bud, suggesting it may function in both germline and somatic stem cells. More interestingly, one Tudor family gene, Bf-Tdrd7, is only expressed maternally and colocalizes with Vasa in germ plasm, suggesting that it may function exclusively in germ cell specification. To evaluate the conservation of the preformation mechanism among amphioxus species, we further analyze Vasa, Nanos, Piwil1, and Tdrd7 expression in two Asian amphioxus species, B. belcheri and B. japonicum. Their maternal transcripts all localize in similar patterns to those seen in B. floridae. In addition, we labeled putative PGCs with Vasa antibody to trace their dynamic distribution in developing larvae.

Conclusions

We identify additional germ plasm components in amphioxus and demonstrate the molecular distinction between the putative germline stem cells and somatic stem cells. Moreover, our results suggest that preformation may be a conserved mechanism for PGC specification among Branchiostoma species. Our Vasa antibody staining results suggest that after the late neurula stage, amphioxus PGCs probably proliferate with the tail bud cells during posterior elongation and are deposited near the forming myomere boundaries. Subsequently, these PGCs would concentrate at the ventral tip of the myoseptal walls to form the gonad anlagen.

A methodical microarray design enables surveying of expression of a broader range of genes in Ciona intestinalis

We provide a new oligo-microarray for Ciona intestinalis, based on the NimbleGen 12-plex×135k format. The array represents 106,285 probes, which is more than double the probe number of the currently available 44k microarray. These probes cover 99.2% of the transcripts in the KyotoHoya (KH) models, published in 2008, and they contain 81.1% of the entries in the UniGene database that are not included in the KH models. In this paper, we show that gene expression levels measured by this new 135k microarray are highly correlated with those obtained by the existing 44k microarray for genes common to both arrays. We also investigated gene expression using samples obtained from the ovary and the neural complex of adult C. intestinalis, showing that the expression of tissue-specific genes is consistent with previous reports. Approximately half of the highly expressed genes identified in the 135k microarray are not included in the previous microarray. The high coverage of gene models by this microarray made it possible to identify splicing variants for a given transcript. The 135k microarray is useful in investigating the functions of genes that are not yet well characterized. Detailed information about this 135k microarray is accessible at no charge from supplemental materials, NCBI Gene Expression Omnibus (GEO), and http://marinegenomics.oist.jp.

AmphiEST: Enabling comparative analysis of ESTs from five developmental stages of amphioxus

Amphioxus has been an important model for understanding the evolution of chordates and origin of vertebrates. Comparative transcriptome analysis can facilitate delineation of gene expression patterns of amphioxus at different developmental stages. So far, however, few such analyses have been performed. Here we have systematically compared amphioxus ESTs from five developmental stages. For the egg, gastrula, neurula, larva and adult stages, amphioxus ESTs were assembled, respectively, into 3364, 3230, 10,299, 4052 and 3866 contigs, and 193, 183, 933, 178 and 151 singlets. 25,796 ORFs were identified, of which 6,529 predicted ORFs (egg, 922; gastrula, 997; neurula, 2,159; larva, 1,014; adult, 1,437) have found matches from the UniProt database. For all five stages, more GO terms were assigned to stage-specific contigs than to repertoire ones. Respectively, 2, 5, 11, 4 and 2 miRNAs were matched with ESTs from the five stages, and 12,232, 10,147, 43,170, 8,049 and 12,811 miRNA target sites were predicted from the 3' UTR sequences. A web interface was developed for visualization of the contigs, miRNAs and annotations (http://www.rich.yunda.org/test/amphioxusest/). The constructed AmphiEST database may serve as an essential resource for future studies of the development and evolution of amphioxus.

Identification of a novel vanadium-binding protein by EST analysis on the most vanadium-rich ascidian, Ascidia gemmata

Ascidians are known to accumulate extremely high levels of vanadium in their blood cells (up to 350 mM). The branchial sac and the intestine are thought to be the first tissues to contact the outer environment and absorb vanadium ions. The concentration of vanadium in the branchial sac and the intestine of the most vanadium-rich ascidian Ascidia gemmata were determined to be 32.4 and 11.9 mM, respectively. Using an expressed sequence tag (EST) analysis of a cDNA library from the intestine of A. gemmata, we determined 960 ESTs and found 55 clones of metal-related gene orthologs, 6 redox-related orthologs, and 18 membrane transporter orthologs. Among them, two genes, which exhibited significant similarity to the vanadium-binding proteins of other vanadium-rich ascidian species, were designated AgVanabin1 and AgVanabin2. Immobilized metal ion affinity chromatography revealed that recombinant AgVanabin1 bound to metal ions with an increasing affinity for Cu(II) > Zn(II) > Co(II) and AgVanabin2 bound to metal ions with an increasing affinity for Cu(II) > Fe(III) > V(IV). To examine the use of AgVanabins for a metal absorption system, we constructed Escherichia coli strains that expressed AgVanabin1 or AgVanabin2 fused to maltose-binding protein and secreted into the periplasmic space. We found that the strain expressing AgVanabin2 accumulated about 13.5 times more Cu(II) ions than the control TB1 strain. Significant accumulation of vanadium was also observed in the AgVanabin2-expressing strain as seen by a 1.5-fold increase.

No chromosomal clustering of housekeeping genes in the marine chordate Ciona intestinalis

Housekeeping genes, widely expressed genes that are required for the basal function of most cell types, are clustered in the human and worm genomes. This arrangement suggests coordinate control of housekeeping gene expression at the chromosomal level. Here we examined whether this notion is applicable to a marine chordate, Ciona intestinalis. Using microarrays, we analyzed genes that were expressed in 11 organs of the adult, including the neural complex, branchial sac, esophagus, stomach, endostyle, intestine, body-wall muscle, heart, blood cells, ovary and testis. This analysis identified 158 genes that are expressed ubiquitously in these organs. These housekeeping genes could be classified into a range of Gene Ontology categories, in particular, ribosomal protein components. Of these 158 genes, we were able to map 141 genes onto the 14 pairs of the C. intestinalis chromosomes. They were distributed rather evenly over all the chromosomes, except for small clusters containing two or three genes. Therefore, the notion of chromosomal clustering of housekeeping genes is not applicable in this chordate.

CIPRO 2.5: Ciona intestinalis protein database, a unique integrated repository of large-scale omics data, bioinformatic analyses and curated annotation, with user rating and reviewing functionality

The Ciona intestinalis protein database (CIPRO) is an integrated protein database for the tunicate species C. intestinalis. The database is unique in two respects: first, because of its phylogenetic position, Ciona is suitable model for understanding vertebrate evolution; and second, the database includes original large-scale transcriptomic and proteomic data. Ciona intestinalis has also been a favorite of developmental biologists. Therefore, large amounts of data exist on its development and morphology, along with a recent genome sequence and gene expression data. The CIPRO database is aimed at collecting those published data as well as providing unique information from unpublished experimental data, such as 3D expression profiling, 2D-PAGE and mass spectrometry-based large-scale analyses at various developmental stages, curated annotation data and various bioinformatic data, to facilitate research in diverse areas, including developmental, comparative and evolutionary biology. For medical and evolutionary research, homologs in humans and major model organisms are intentionally included. The current database is based on a recently developed KH model containing 36 034 unique sequences, but for higher usability it covers 89 683 all known and predicted proteins from all gene models for this species. Of these sequences, more than 10 000 proteins have been manually annotated. Furthermore, to establish a community-supported protein database, these annotations are open to evaluation by users through the CIPRO website. CIPRO 2.5 is freely accessible at http://cipro.ibio.jp/2.5.

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.

Computational analysis of Ciona intestinalis operons

Operons are clusters of genes that are co-regulated from a common promoter. Operons are typically associated with prokaryotes, although a small number of eukaryotes have been shown to possess them. Among metazoans, operons have been extensively characterized in the nematode Caenorhabditis elegans in which ∼15% of the total genes are organized into operons. The most recent genome assembly for the ascidian Ciona intestinalis placed ∼20% of the genes (2909 total) into 1310 operons. The majority of these operons are composed of two genes, while the largest are composed of six. Here is reported a computational analysis of the genes that comprise the Ciona operons. Gene ontology (GO) terms were identified for about two-thirds of the operon-encoded genes. Using the extensive collection of public EST libraries, estimates of temporal patterns of gene expression were generated for the operon-encoded genes. Lastly, conservation of operons was analyzed by determining how many operon-encoded genes were present in the ascidian Ciona savignyi and whether these genes were organized in orthologous operons. Over 68% of the operon-encoded genes could be assigned one or more GO terms and 697 of the 1310 operons contained genes in which all genes had at least one GO term. Of these 697 operons, GO terms were shared by all of the genes within 146 individual operons, suggesting that most operons encode genes with unrelated functions. An analysis of operon gene expression from nine different EST libraries indicated that for 587 operons, all of the genes that comprise an individual operon were expressed together in at least one EST library, suggesting that these genes may be co-regulated. About 50% (74/146) of the operons with shared GO terms also showed evidence of gene co-regulation. Comparisons with the C. savignyi genome identified orthologs for 1907 of 2909 operon genes. About 38% (504/1310) of the operons are conserved between the two Ciona species. These results suggest that like C. elegans, operons in Ciona are comprised of a variety of genes that are not necessarily related in function. The genes in only 50% of the operons appear to be co-regulated, suggesting that more complex gene regulatory mechanisms are likely operating.

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.

Stress response in the ascidian Ciona intestinalis: transcriptional profiling of genes for the heat shock protein 70 chaperone system under heat stress and endoplasmic reticulum stress

The genome of Ciona intestinalis contains eight genes for HSP70 superfamily proteins, 36 genes for J-proteins, a gene for a J-like protein, and three genes for BAG family proteins. To understand the stress responses of genes in the HSP70 chaperone system comprehensively, the transcriptional profiles of these 48 genes under heat stress and endoplasmic reticulum (ER) stress were studied using real-time reverse transcriptase-polymerase chain reaction (RT-PCR) analysis. Heat stress treatment increased the messenger RNA (mRNA) levels of six HSP70 superfamily genes, eight J-protein family genes, and two BAG family genes. In the cytoplasmic group of the DnaK subfamily of the HSP70 family, Ci-HSPA1/6/7-like was the only heat-inducible gene and Ci-HSPA2/8 was the only constitutively active gene which showed striking simplicity in comparison with other animals that have been examined genome-wide so far. Analyses of the time course and temperature dependency of the heat stress responses showed that the induction of Ci-HSPA1/6/7-like expression rises to a peak after heat stress treatment at 28 degrees C (10 degrees C upshift from control temperature) for 1 h. ER stress treatment with Brefeldin A, a drug that is known to act as ER stress inducer, increased the mRNA levels of four HSP70 superfamily genes and four J-protein family genes. Most stress-inducible genes are conserved between Ciona and vertebrates, as expected from a close evolutionary relationship between them. The present study characterized the stress responses of HSP70 chaperone system genes in Ciona for the first time and provides essential data for comprehensive understanding of the functions of the HSP70 chaperone system.

Ciona intestinalis and Oxycomanthus japonicus, representatives of marine invertebrates

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

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

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

Domain shuffling and the evolution of vertebrates

The evolution of vertebrates has included a number of important events: the development of cartilage, the immune system, and complicated craniofacial structures. Here, we examine domain shuffling as one of the mechanisms that contributes novel genetic material required for vertebrate evolution. We mapped domain-shuffling events during the evolution of deuterostomes with a focus on how domain shuffling contributed to the evolution of vertebrate- and chordate-specific characteristics. We identified approximately 1000 new domain pairs in the vertebrate lineage, including approximately 100 that were shared by all seven of the vertebrate species examined. Some of these pairs occur in the protein components of vertebrate-specific structures, such as cartilage and the inner ear, suggesting that domain shuffling made a marked contribution to the evolution of vertebrate-specific characteristics. The evolutionary history of the domain pairs is traceable; for example, the Xlink domain of aggrecan, one of the major components of cartilage, was originally utilized as a functional domain of a surface molecule of blood cells in protochordate ancestors, and it was recruited by the protein of the matrix component of cartilage in the vertebrate ancestor. We also identified genes that were created as a result of domain shuffling in ancestral chordates. Some of these are involved in the functions of chordate structures, such as the endostyle, Reissner's fiber of the neural tube, and the notochord. Our analyses shed new light on the role of domain shuffling, especially in the evolution of vertebrates and chordates.

The Fox genes of Branchiostoma floridae

The Fox genes are united by encoding a fork head domain, a deoxyribonucleic acid (DNA)-binding domain of the winged-helix type that marks these genes as encoding transcription factors. Vertebrate Fox genes are classified into 23 subclasses named from FoxA to FoxS. We have surveyed the genome of the amphioxus Branchiostoma floridae, identifying 32 distinct Fox genes representing 21 of these 23 subclasses. The missing subclasses, FoxR and FoxS, are specific to vertebrates, and in addition, B. floridae has one further group, FoxAB, that is not found in vertebrates. Hence, we conclude B. floridae has maintained a high level of Fox gene diversity. Expressed sequence tag and complementary DNA sequence data support the expression of 23 genes. Several linkages between B. floridae Fox genes were noted, including some that have evolved relatively recently via tandem duplication in the amphioxus lineage and others that are more ancient.

A cDNA resource for the cephalochordate amphioxus Branchiostoma floridae

Cephalochordates are the basal invertebrate chordates within the phylum Chordata. They are widely used as a model system for research in evolutionary developmental biology (EvoDevo) to understand the basic patterning mechanisms for the chordate body plan and the origin of vertebrates. Recently, the genome of the cephalochordate Branchiostoma floridae was sequenced, which further brings this organism to the front for comparative genomic studies. In this paper, we report the generation of large-scale 5'- and 3'-expressed sequence tags (ESTs) from B. floridae and the complementary deoxyribonucleic acid (cDNA) resource for this species. Both 5'- and 3'-ESTs were sequenced for approximately 140,000 cDNA clones derived from five developmental stages, and the cDNA clones were subsequently grouped into independent clusters using 3'-EST sequences. We identified 21,229 cDNA clusters, and each corresponds to a unique transcript species from B. floridae. We then chose 24,020 cDNA clones representing all of these 21,229 clusters to generate the "Branchiostoma floridae Gene Collection Release 1." We also constructed a database with a searchable interface for this EST dataset and the related information on "Branchiostoma floridae Gene Collection Release 1." This set of cDNA clones along with our cDNA database will serve as an important resource for future research in this basal chordate. This Gene Collection and the original 140,000 individual cDNA clones are available to the research community upon request.

Novel genes involved in canonical Wnt/beta-catenin signaling pathway in early Ciona intestinalis embryos

We report here characterization of five genes for novel components of the canonical Wnt/beta-catenin signaling pathway. These genes were identified in the ascidian Ciona intestinalis through a loss-of-function screening for genes required for embryogenesis with morpholinos, and four of them have counterparts in vertebrates. The five genes we studied are as follows: Ci-PGAP1, a Ciona orthologue of human PGAP1, which encodes GPI (glycosylphosphatidylinositol) inositol-deacylase, Ci-ZF278, a gene encoding a C2H2 zinc-finger protein, Ci-C10orf11, a Ciona orthologue of human C10orf11 that encodes a protein with leucine-rich repeats, Ci-Spatial/C4orf17, a single counterpart for two human genes Spatial and C4orf17, and Ci-FLJ10634, a Ciona orthologue of human FLJ10634 that encodes a member of the J-protein family. Knockdown of each of the genes mimicked beta-catenin knockdown and resulted in suppression of the expression of beta-catenin downstream genes (Ci-FoxD, Ci-Lhx3, Ci-Otx and Ci-Fgf9/16/20) and subsequent endoderm formation. For every gene, defects in knockdown embryos were rescued by overexpression of a constitutively active form, but not wild-type, of Ci-beta-catenin. Dosage-sensitive interactions were found between Ci-beta-catenin and each of the genes. These results suggest that these five genes act upstream of or parallel to Ci-beta-catenin in the Wnt/beta-catenin signaling pathway in early Ciona embryos.

Effects of 5-aza-2'-deoxycytidine on the gene expression profile during embryogenesis of the Ascidian ciona intestinalis: a microarray analysis

DNA methylation is an important epigenetic factor that participates in silencing genes. Genomic approaches to studying DNA methylation promise to be particularly fruitful, since DNA methylation is involved in global control of gene expression in many organisms. With its draft genome completed and a large quantity of available cDNA data, Ciona intestinalis is newly emerging as an invaluable model organism for investigating genome-wide gene expression and function. Here we examine the effects of 5-aza-2'-deoxycytidine (5-aza-CdR), a chemical that blocks CpG methylation, on the gene expression profile of early C. intestinalis embryos, using oligonucleotide-based microarray analysis. Embryos treated with 5-aza-CdR show delayed gastrulation and are developmentally arrested at the neurula stage. They subsequently lose cellular adhesion and finally die. Apoptosis was not detected in these embryos by TUNEL staining at 12 h, indicating that the defects observed did not result from 5-aza-CdR-induced apoptosis. Gene expression profiles of 12-h-old 5-aza-CdR-treated embryos compared to wild-type revealed 91 upregulated genes and 168 downregulated genes. Although nearly half of these encoded proteins with unknown functions, several encoded cell-signaling molecules and transcription factors. In addition, genes associated with the stress response and cell defense were upregulated, whereas genes involved in cell adhesion were downregulated.

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.

Characterization of an 18,166 EST dataset for cassava (Manihot esculenta Crantz) enriched for drought-responsive genes

Cassava (Manihot esculenta Crantz) is a staple food for over 600 million people in the tropics and subtropics and is increasingly used as an industrial crop for starch production. Cassava has a high growth rate under optimal conditions but also performs well in drought-prone areas and on marginal soils. To increase the tools for understanding and manipulating drought tolerance in cassava, we generated expressed sequence tags (ESTs) from normalized cDNA libraries prepared from dehydration-stressed and control well-watered tissues. Analysis of a total of 18,166 ESTs resulted in the identification of 8,577 unique gene clusters (5,383 singletons and 3,194 clusters). Functional categories could be assigned to 63% of the unigenes, while another approximately 11% were homologous to hypothetical genes with unclear functions. The remaining approximately 26% were not significantly homologous to sequences in public databases suggesting that some may be novel and putatively specific to cassava. The dehydration-stressed library uncovered numerous ESTs with recognized roles in drought-responses, including those that encode late-embryogenesis-abundant proteins thought to confer osmoprotective functions during water stress, transcription factors, heat-shock proteins as well as proteins involved in signal transduction and oxidative stress. The unigene clusters were screened for short tandem repeats for further development as microsatellite markers. A total of 592 clusters contained 646 repeats, representing 3.3% of the ESTs queried. The ESTs presented here are the first dehydration stress transcriptome of cassava and can be utilized for the development of microarrays and gene-derived molecular markers to further dissect the molecular basis of drought tolerance in cassava.

Analysis of large scale expression sequenced tags (ESTs) from the anural ascidian, Molgula tectiformis

Anural ascidians show embryogenesis during which tail formation does not take place. This mode of development is a derived character acquired several times independently in ascidian evolution. We identified approximately 20,000 each ESTs (i. e. 10,000 clones each were sequenced from both 5' and 3' ends) of adult gonads, cleaving-embryos, gastrulae/neurulae, embryos before hatching, and hatched larvae of the anural ascidian Molgula tectiformis, in order to comprehensively investigate the molecular mechanism of tailless evolution. Analyses of these ESTs showed that in this species, (1) the expression of embryonic/larval muscle structural genes which are expressed abundantly during embryogenesis of the urodele ascidian Ciona intestinalis, is suppressed; (2) genes that encode proteins with no similarity to known proteins of other organisms are abundantly expressed; (3) genes that show similarity with those up-regulated at metamorphosis in urodele ascidians are up-regulated within several hours after hatching; and (4) 15 of 35 putative orthologues of the downstream components of Brachyury, a key transcription factor for ascidian notochord formation, were found in the ESTs, even though differentiation of notochord is suppressed in this species. We discuss these remarkable results that allow insight into the molecular mechanism(s) responsible for the anural mode of ascidian development.

Origin of the vertebrate visual cycle: III. Distinct distribution of RPE65 and beta-carotene 15,15'-monooxygenase homologues in Ciona intestinalis

We previously identified three genes that encode putative visual cycle proteins that are homologues of retinal G-protein coupled receptor (Ci-opsin3), cellular retinaldehyde-binding protein (Ci-CRALBP) and beta-carotene 15,15'-monooxygenase (Ci-BCO) in the ascidian Ciona intestinalis. Ci-opsin3 and Ci-CRALBP are localized in both ocellus photoreceptor cells and surrounding non-photoreceptor cells in the brain vesicle of the larva. In the present study, we investigated the possible role and evolutionary origin of the BCO/RPE65 family in the visual cycle by analyzing Ci-BCO localization by immunohistochemistry and by identifying a novel gene that encodes a homologue of retinal pigment epithelium-specific 65 kDa protein (Ci-RPE65) in C. intestinalis. In situ hybridization and expressed sequence tag (EST) profiles consistently suggest that Ci-RPE65 is not significantly expressed in the ocellus and brain vesicle of the larva. Ci-RPE65 is expressed in the neural complex, a photoreceptor organ of the adult ascidian, at a level comparable to that of Ci-opsin3 and Ci-CRALBP. Ci-RPE65 is also expressed in various adult tissues, including the gill, body wall and intestine, suggesting that Ci-RPE65 plays a role in addition to that in the visual cycle. In contrast, Ci-BCO is predominantly localized in ocellus photoreceptor cells of the larva. The larval visual cycle seems to use Ci-opsin3 as a photo-isomerase. Our results also suggest that the RPE65-dependent visual cycle is used in the adult photoreceptors of a primitive chordate.

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

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

Functional proteomics inCiona intestinalis: A breakthrough in the exploration of the molecular and cellular mechanism of ascidian development

Ascidians have been providing a unique experimental system for a variety of fields, including reproductive biology, developmental biology, neurobiology, immunology, and evolutional biology. Recent progress in the genome sequencing of Ciona intestinalis has led to the development of a great tool for investigating the gene functions and expressions involved in several biological events in ascidians. The disclosure of genomic information has ushered in the postgenomic era, spearheaded by extensive protein analysis. The characterization of the function, localization, and molecular interaction of cellular proteins results in a more direct description of the molecular mechanism underlying several biological processes. Proteomics in ascidians, however, has just recently appeared and is not well established yet. In this study, we give an outline of the technical processes used in proteomics and review the recent status of ascidian proteomics.

Isolation and characterization of a novel uronic acid-containing acidic glycosphingolipid from the ascidian Halocynthia roretzi

A novel uronic acid-containing glycosphingolipid (UGL-1) was isolated from the ascidian Halocynthia roretzi. UGL-1 was prepared from chloroform-methanol extracts and purified by the use of successive column chromatography on DEAE-Sephadex, Florisil, and Iatrobeads. Chemical structural analysis was performed using methylation analysis, gas chromatography, gas chromatography-mass spectrometry, matrix-assisted laser-desorption/ionization time-of-flight mass spectrometry, and 1H-NMR spectra. The chemical structure of UGL-1 was determined to be a glucuronic acid-containing glycosphingolipid, Galbeta1-4(Fucalpha1-3)GlcAbeta1-1Cer. The ceramide component was composed of C16:0 and C18:0 acids and C16-, C17-, and C18-phytosphingosines as major components.

A Novel Sulfatide, GlcCer-I6 Sulfate, from the Ascidian Ciona intestinalis

A novel sulfated glycosphingolipid, SGL-1, was isolated from the ascidian Ciona intestinalis, prepared from chloroform/methanol extracts and fractionated successively on DEAE Sephadex-A25, Florisil and Iatrobeads column chromatographies. Chemical structural analysis was performed using methylation analysis, gas-liquid chromatography, combined gas-liquid chromatography-mass spectrometry, matrix-assisted laser desorption/ionization time-of-flight mass spectrometry and (1)H-NMR spectroscopy. This chemical structure is presented as GlcCer I(6)-Sulfate. The ceramide moiety was specified by t16:0, t17:0, br,t17:0, t18:0 and br,t18:0 as sphingoids, and 2-hydroxy, saturated fatty acids as represented by docosanoic and tetracosanoic acids.

Isolation and phylogeny of novel cytochrome P450 genes from tunicates (Ciona spp.): A CYP3 line in early deuterostomes?

Cytochromes P450 (CYPs) form a gene superfamily involved in the biotransformation of numerous endogenous and exogenous natural and synthetic compounds. In humans, CYP3A4 is regarded as one of the most important CYPs due to its abundance in liver and its capacity to metabolize more than 50% of all clinically used drugs. It has been suggested that all CYP3s arose from a common ancestral gene lineage that diverged between 800 and 1100 million years ago, before the deuterostome-protostome split. While CYP3s are well known in mammals and have been described in lower vertebrates, they have not been reported in non-vertebrate deuterostomes. Members of the genus Ciona belong to the tunicates, whose lineage is thought to be the most basal among the chordates, and from which the vertebrate line diverged. Here we describe the cloning, exon-intron structure, phylogeny, and estimated expression of four novel genes from Ciona intestinalis. We also describe the gene structure and phylogeny of homologous genes in Ciona savignyi. Comparing these genes with other members of the CYP clan 3, show that the Ciona sequences bear remarkable similarity to vertebrate CYP3A genes, and may be an early deuterostome CYP3 line.

Genomic overview of mRNA 5'-leader trans-splicing in the ascidian Ciona intestinalis

Although spliced leader (SL) trans-splicing in the chordates was discovered in the tunicate Ciona intestinalis there has been no genomic overview analysis of the extent of trans-splicing or the make-up of the trans-spliced and non-trans-spliced gene populations of this model organism. Here we report such an analysis for Ciona based on the oligo-capping full-length cDNA approach. We randomly sampled 2078 5'-full-length ESTs representing 668 genes, or 4.2% of the entire genome. Our results indicate that Ciona contains a single major SL, which is efficiently trans-spliced to mRNAs transcribed from a specific set of genes representing approximately 50% of the total number of expressed genes, and that individual trans-spliced mRNA species are, on average, 2-3-fold less abundant than non-trans-spliced mRNA species. Our results also identify a relationship between trans-splicing status and gene functional classification; ribosomal protein genes fall predominantly into the non-trans-spliced category. In addition, our data provide the first evidence for the occurrence of polycistronic transcription in Ciona. An interesting feature of the Ciona polycistronic transcription units is that the great majority entirely lack intercistronic sequences.

An integrated database of the ascidian, Ciona intestinalis: towards functional genomics

An integrated genome database is essential for future studies of functional genomics. In this study, we update cDNA and genomic resources of the ascidian, Ciona intestinalis, and provide an integrated database of the genomic and cDNA data by extending a database published previously. The updated resources include over 190,000 ESTs (672,396 in total together with the previous ESTs) and over 1,000 full-insert sequences (6,773 in total). In addition, results of mapping information of the determined scaffolds onto chromosomes, ESTs from a full-length enriched cDNA library for indication of precise 5'-ends of genes, and comparisons of SNPs and indels among different individuals are integrated into this database, all of these results being reported recently. These advances continue to increase the utility of Ciona intestinalis as a model organism whilst the integrated database will be useful for researchers in comparative and evolutionary genomics.

Ci-Rga, a gene encoding an MtN3/saliva family transmembrane protein, is essential for tissue differentiation during embryogenesis of the ascidian Ciona intestinalis

A novel gene (Ci-Rga) essential for tissue differentiation during embryogenesis of the ascidian Ciona intestinalis is reported here. This gene was identified through functional screening of Ciona genes required for development by translational inhibition experiments with morpholino antisense oligonucleotides. The deduced protein of Ci-Rga contains two copies of a domain with unknown function called the MtN3/saliva domain. Phylogenetic analysis showed that Ci-Rga belongs to the MtN3/saliva family of genes conserved among metazoans and plants, and is an ortholog of mouse Rga (Recombination-activating gene 1 gene activation). During Ciona embryogenesis, both maternal and zygotic transcripts of Ci-Rga were expressed. Translational inhibition of Ci-Rga with specific morpholino resulted in abnormal embryos in which the cleavage pattern became atypical and expression of marker genes for each of the six major tissues, namely the endoderm, muscle, mesenchyme, notochord, neural tissue, and epidermis, was lost or suppressed at the tailbud stage. Although differentiation of all the six major tissues was affected by Ci-Rga knock-down, the degree of abnormalities and the timing of appearance of abnormalities were different among tissues. Expression analysis of developmentally important genes involved in the fate specification, such as Ci-Bra, Ci-Twist-like1a, Ci-Otx, Ci-Fgf9/16/20, Ci-Lhx3, Ci-FoxD, and Ci-Tbx6b, showed that an initial step of the fate specification of notochord, mesenchyme, and neural tissue, but not of endoderm or muscle, is impaired in the knock-down embryo. These results showed that Ci-Rga is a multifunctional gene essential for tissue differentiation during embryogenesis, and is primarily required for the fate specification of notochord, mesenchyme, and neural tissue, and provide some insights into the function of this little-known group of genes.

A Genome-Wide Survey of Genes for Enzymes Involved in Pigment Synthesis in an Ascidian, Ciona intestinalis

The draft genome sequence and a large quantity of EST and cDNA information are now available for the ascidian Ciona intestinalis. In the present study, genes involved in pigment synthesis pathways were identified in the decoded genome of Ciona, and information about these genes was obtained from available EST and cDNA sequences. It was found that the Ciona genome contains orthologous genes for each enzyme of the melanin, pteridine, ommochrome, papiliochrome, and heme synthesis pathways. Several appear as independent duplications in the Ciona genome. Because cDNA clones for all but two of these genes have already been isolated by the cDNA project, C. intestinalis will provide an experimental system to explore molecular mechanisms underlying color patterns, through future genome-wide studies.

Maternal determinants and mRNAs in the cortex of ascidian oocytes, zygotes and embryos

The peripheral region of ascidian oocytes and zygotes contains five determinants for morphogenesis and differentiation of the embryo. The determinant for the 24 primary muscle cells of the tadpole, macho1, is one of several cortical mRNAs localized in a gradient along the animal-vegetal axis in the oocyte. After fertilization these mRNAs, together with cortical endoplasmic reticulum (cER) and a subcortical mitochondria-rich domain (myoplasm), relocate in two major reorganization phases forming the posterior plasm (postplasm) of the zygote. At the 8-cell stage cortical mRNAs concentrate in a macroscopic cortical structure called the centrosome-attracting body (CAB), forming a characteristic posterior end mark (PEM) in the two posterior vegetal blastomeres. We propose to call the numerous mRNAs showing this particular cortical localization in the posterior region of the embryo postplasmic/PEM RNAs and suggest a nomemclature. We do not know how postplasmic/PEM RNAs reach their polarized distribution in the oocyte cortex but at least PEM1 and macho1 (and probably others) bind to the network of cER retained in isolated cortical fragments. We propose that after fertilization, these postplasmic/PEM mRNAs move in the zygote cortex together with the cER network (cER/mRNA domain) via microfilament- and microtubule-driven translocations. The cER/mRNA domain is localized posteriorly at the time of first cleavage and distributed equally between the first two blastomeres. After the third cleavage, the cER/mRNA domain and dense particles compact to form the CAB in posterior vegetal blastomeres of the 8-cell stage. We discuss the identity of postplasmic/PEM RNAs, how they localize, anchor, relocate and may be translated. We also examine their roles in unequal cleavage and as a source of posterior morphogenetic and differentiation factors.

openSputnik--a database to ESTablish comparative plant genomics using unsaturated sequence collections

The public expressed sequence tag collections are continually being enriched with high-quality sequences that represent an ever-expanding range of taxonomically diverse plant species. While these sequence collections provide biased insight into the populations of expressed genes available within individual species and their associated tissues, the information is conceivably of wider relevance in a comparative context. When we consider the available expressed sequence tag (EST) collections of summer 2004, most of the major plant taxonomic clades are at least superficially represented. Investigation of the five million available plant ESTs provides a wealth of information that has applications in modelling the routes of plant genome evolution and the identification of lineage-specific genes and gene families. Over four million ESTs from over 50 distinct plant species have been collated within an EST analysis pipeline called openSputnik. The ESTs were resolved down into approximately one million unigene sequences. These have been annotated using orthology-based annotation transfer from reference plant genomes and using a variety of contemporary bioinformatics methods to assign peptide, structural and functional attributes. The openSputnik database is available at http://sputnik.btk.fi.

Computational discovery of DNA motifs associated with cell type-specific gene expression in Ciona

Temporally and spatially co-expressed genes are expected to be regulated by common transcription factors and therefore to share cis-regulatory elements. In the ascidian Ciona intestinalis, the whole-genome sequences and genome-scale gene expression profiles allow the use of computational techniques to investigate cis-elements that control transcription. We collected 5' flanking sequences of 50 tissue-specific genes from genome databases of C. intestinalis and a closely related species Ciona savignyi. We searched for DNA motifs over-represented in upstream regions of a group of co-expressed genes. Several motifs were distributed predominantly in upstream regions of photoreceptor, pan-neuronal, or muscle-specific gene groups. One muscle-specific motif, M2, was distributed preferentially in regions from -200 to -100 bp relative to the translational start sites. Promoters of muscle-specific genes of C. intestinalis were isolated, connected with a green fluorescent protein gene (GFP), and introduced into C. intestinalis embryos. In muscle cells, these promoters specifically drove GFP expression, which mutations of the M2 sites greatly reduced. When M2 sites were located upstream of a basal promoter, the reporter GFP was specifically expressed in muscle cells. These results suggest the validity of our computational prediction of cis-regulatory elements. Thus, bioinformatics can help identify cis-regulatory elements involved in chordate development.

Three distinct lineages of mesenchymal cells in Ciona intestinalis embryos demonstrated by specific gene expression

The ascidian embryonic mesenchyme, comprising about 900 cells, forms mesodermal tissues or organs of the adult body after metamorphosis. The mesenchyme originates from the A7.6 [trunk lateral cells (TLCs)], B7.7, and B8.5 blastomeres of the 110-cell stage embryo. Previous studies showed that FGF9/16/20 is required for specification of the mesenchyme in Ciona embryos and that two different (A7.6 and B8.5/B7.7) but partially overlapping molecular mechanisms are associated with the expression of a basic helix-loop-helix (bHLH) transcription factor gene, Twist-like1, in the mesenchymal precursors, which triggers the differentiation process of mesenchyme cells. In the present study, we examined whether the three embryonic lineages express the same mesenchyme-specific structural genes under the control of a common mechanism or whether the three lineages are characterized by the expression of genes specific to each of the lineages. We characterized nine mesenchyme-specific genes in Ciona embryos and found that five were expressed in A7.6/B8.5/B7.7, two in B8.5/B7.7, and two in B7.7 only. FGF9/16/20 and Twist-like1 were required for the expression of all the mesenchyme-specific genes, except for three A7.6/B8.5/B7.7-specific genes in A7.6 progenitors. Overexpression of FGF9/16/20 or Twist-like1 upregulated the expression of A7.6/B8.5/B7.7- and B8.5/B7.7-specific genes, while it downregulated the expression of B7.7-specific genes. These results provide evidence that the differentiation of each of the three mesenchyme lineages of Ciona embryos is characterized by the expression of a specific set of genes, whose expression is controlled differentially.

An enhancer trap in the ascidian Ciona intestinalis identifies enhancers of its Musashi orthologous gene

The enhancer trap technique, established in Drosophila melanogaster, is a very sophisticated tool. Despite its usefulness, however, there have been very few reports on enhancer traps in other animals. The ascidian Ciona intestinalis, a splendid experimental system for developmental biology, provides good material for developmental genetics. Recently, germline transgenesis of C. intestinalis has been achieved using the Tc1/mariner superfamily transposon Minos. During the course of that study, one Minos insertion line that showed a different GFP expression pattern from other lines was isolated. One fascinating possibility is that an enhancer trap event occurred in this line. Here we show that a Minos insertion in the Ci-Musashi gene was responsible for the altered GFP expression. Ci-Musashi showed a similar expression pattern to GFP. In addition, introns of Ci-Musashi have enhancer activity that can alter the expression pattern of nearby genes to resemble that of GFP in this line. These results clearly demonstrate that an enhancer trap event that entrapped enhancers of Ci-Musashi occurred in C. intestinalis.

The neurobiology of the ascidian tadpole larva: recent developments in an ancient chordate

With little more than 330 cells, two thirds within the sensory vesicle, the CNS of the tadpole larva of the ascidian Ciona intestinalis provides us with a chordate nervous system in miniature. Neurulation, neurogenesis and its genetic bases, as well as the gene expression territories of this tiny constituency of cells all follow a chordate plan, giving rise in some cases to frank structural homologies with the vertebrate brain. Recent advances are fueled by the release of the genome and EST expression databases and by the development of methods to transfect embryos by electroporation. Immediate prospects to test the function of neural genes are based on the isolation of mutants by classical genetics and insertional mutagenesis, as well as by the disruption of gene function by morpholino antisense oligo-nucleotides. Coupled with high-speed video analysis of larval swimming, optophysiological methods offer the prospect to analyze at single-cell level the function of a CNS built on a vertebrate plan.

T-box genes in the ascidianCiona intestinalis: Characterization of cDNAs and spatial expression

Members of the T-box family of transcription factors share an evolutionarily conserved DNA-binding domain and play significant roles in various processes of embryonic development. Vertebrate T-box genes are categorized into the following five major subfamilies (eight groups), depending on sequence similarities: Brachyury, Tbx1 (Tbx1/10, Tbx15/18/22, Tbx20), Tbx2/3/4/5 (Tbx2/3 and Tbx4/5), Tbx6, and Tbr/Eomes/TBX21. Ascidians are primitive chordates, and their tadpole larva are considered to represent the simplified and basic body plan of vertebrates. In addition, it has been revealed that the ascidian genome contains the basic ancestral complement of genes involved in development. The present characterization of cDNAs and survey of the Ciona intestinalis draft genome demonstrated that the Ciona genome contains a single copy gene for each of the Brachyury, Tbx1/10, Tbx15/18/22, Tbx20, Tbx2/3, and Tbr/Eomes/TBX21 groups, and at least three copies of the Tbx6 subfamily. Each of the Ciona T-box genes shows a characteristic expression pattern, although that of Tbx20 was not determined in the present study. These results provide basic information that will be useful for future studies of the function of each gene, genetic cascades of different T-box genes, and genome-wide surveys of evolutionary changes in the T-box gene structure and organization in this primitive chordate.

The evolutionary origin of animal cellulose synthase

Urochordates are the only animals that produce cellulose, a polysaccharide existing primarily in the extracellular matrices of plant, algal, and bacterial cells. Here we report a Ciona intestinalis homolog of cellulose synthase, which is the core catalytic subunit of multi-enzyme complexes where cellulose biosynthesis occurs. The Ciona cellulose synthase gene, Ci-CesA, is a fusion of a cellulose synthase domain and a cellulase (cellulose-hydrolyzing enzyme) domain. Both the domains have no animal homologs in public databases. Exploiting this fusion of atypical genes, we provided evidence of a likely lateral transfer of a bacterial cellulose synthase gene into the urochordate lineage. According to fossil records, this likely lateral acquisition of the cellulose synthase gene may have occurred in the last common ancestor of extant urochordates more than 530 million years ago. Whole-mount in situ hybridization analysis revealed the expression of Ci-CesA in C. intestinalis embryos, and the expression pattern of Ci-CesA was spatiotemporally consistent with observed cellulose synthesis in vivo. We propose here that urochordates may use a laterally acquired "homologous" gene for an analogous process of cellulose synthesis.