The characterisation of six ADAMTS proteases in the basal chordate Ciona intestinalis provides new insights into the vertebrate ADAMTS family

ADAMTS, constituting a recently discovered family of secreted zinc-dependent metalloproteases, have been shown to have critical physiological roles through identification of a number of natural animal and human gene mutations. The identification of six ADAMTS genes in the basal chordate Ciona intestinalis provides new insight into how, when and in what order the vertebrate orthologues have evolved. The phylogenetic assignments, based on sequences conserved across all genes, are supported by conserved domain structures within defined sub-families. The phylogeny and the frequent localisation of ADAMTS genes in paralogous regions of the genome are consistent with the vertebrate lineages having arisen by large scale or genome duplication. The high level of conservation in the protease active site of vertebrate orthologues within some sub-families suggests subfunctionalisation, whereas the greater divergence in others would favour the evolution of novel substrate specificities and these observations are borne-out where substrate-specificity is known. The expansion and sub-specialization of the ADAMTS family is a component of the increased complexity of extracellular matrix that is associated with the evolution of vertebrates.     

A cDNA resource from the basal chordate Ciona intestinalis

The genome of the basal choradate Ciona intestinalis contains a basic set of genes with less redundancy compared to the vertebrate genome. Extensive EST analyses, cDNA sequencing, and clustering yielded "Ciona intestinalis Gene Collection Release 1," which contains cDNA clones for 13,464 genes, covering nearly 85% of the Ciona mRNA species. This release is ready for use in cDNA cloning, micro/macroarray analysis, and other comprehensive genome-wide analyses for further molecular studies of basal chordates.     

The evolution of the vertebrate metzincins; insights from Ciona intestinalis and Danio rerio

Background

The metzincins are a large gene superfamily of proteases characterized by the presence of a zinc protease domain, and include the ADAM, ADAMTS, BMP1/TLL, meprin and MMP genes. Metzincins are involved in the proteolysis of a wide variety of proteins, including those of the extracellular matrix. The metzincin gene superfamily comprises eighty proteins in the human genome and ninety-three in the mouse. When and how the level of complexity apparent in the vertebrate metzincin gene superfamily arose has not been determined in detail. Here we present a comprehensive analysis of vertebrate metzincins using genes from both Ciona intestinalis and Danio rerio to provide new insights into the complex evolution of this gene superfamily.

Results

We have identified 19 metzincin genes in the ciona genome and 83 in the zebrafish genome. Phylogenetic analyses reveal that the expansion of the metzincin gene superfamily in vertebrates has occurred predominantly by the simple duplication of pre-existing genes rather than by the appearance and subsequent expansion of new metzincin subtypes (the only example of which is the meprin gene family). Despite the number of zebrafish metzincin genes being relatively similar to that of tetrapods (e.g. man and mouse), the pattern of gene retention and loss within these lineages is markedly different. In addition, we have studied the evolution of the related TIMP gene family and identify a single ciona and four zebrafish TIMP genes.

Conclusion

The complexity seen in the vertebrate metzincin gene families was mainly acquired during vertebrate evolution. The metzincin gene repertoire in protostomes and invertebrate deuterostomes has remained relatively stable. The expanded metzincin gene repertoire of extant tetrapods, such as man, has resulted largely from duplication events associated with early vertebrate evolution, prior to the sarcopterygian-actinopterygian split. The teleost repertoire of metzincin genes in part parallels that of tetrapods but has been significantly modified, perhaps as a consequence of a teleost-specific duplication event.     

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.     

Bacterial diversity associated with the tunic of the model chordate Ciona intestinalis

The sea squirt Ciona intestinalis is a well-studied model organism in developmental biology, yet little is known about its associated bacterial community. In this study, a combination of 454 pyrosequencing of 16S ribosomal RNA genes, catalyzed reporter deposition-fluorescence in situ hybridization and bacterial culture were used to characterize the bacteria living inside and on the exterior coating, or tunic, of C. intestinalis adults. The 454 sequencing data set demonstrated that the tunic bacterial community structure is different from that of the surrounding seawater. The observed tunic bacterial consortium contained a shared community of <10 abundant bacterial phylotypes across three individuals. Culture experiments yielded four bacterial strains that were also dominant groups in the 454 sequencing data set, including novel representatives of the classes Alphaproteobacteria and Flavobacteria. The relatively simple bacterial community and availability of dominant community members in culture make C. intestinalis a promising system in which to investigate functional interactions between host-associated microbiota and the development of host innate immunity.     

MGEScan-non-LTR: computational identification and classification of autonomous non-LTR retrotransposons in eukaryotic genomes

Computational methods for genome-wide identification of mobile genetic elements (MGEs) have become increasingly necessary for both genome annotation and evolutionary studies. Non-long terminal repeat (non-LTR) retrotransposons are a class of MGEs that have been found in most eukaryotic genomes, sometimes in extremely high numbers. In this article, we present a computational tool, MGEScan-non-LTR, for the identification of non-LTR retrotransposons in genomic sequences, following a computational approach inspired by a generalized hidden Markov model (GHMM). Three different states represent two different protein domains and inter-domain linker regions encoded in the non-LTR retrotransposons, and their scores are evaluated by using profile hidden Markov models (for protein domains) and Gaussian Bayes classifiers (for linker regions), respectively. In order to classify the non-LTR retrotransposons into one of the 12 previously characterized clades using the same model, we defined separate states for different clades. MGEScan-non-LTR was tested on the genome sequences of four eukaryotic organisms, Drosophila melanogaster, Daphnia pulex, Ciona intestinalis and Strongylocentrotus purpuratus. For the D. melanogaster genome, MGEScan-non-LTR found all known ‘full-length’ elements and simultaneously classified them into the clades CR1, I, Jockey, LOA and R1. Notably, for the D. pulex genome, in which no non-LTR retrotransposon has been annotated, MGEScan-non-LTR found a significantly larger number of elements than did RepeatMasker, using the current version of the RepBase Update library. We also identified novel elements in the other two genomes, which have only been partially studied for non-LTR retrotransposons.     

Direct examination of chromosomal clustering of organ-specific genes in the chordate Ciona intestinalis

One of challenges in the field of developmental biology is to understand how spatially and/or temporally coordinated expression of genes is controlled at the chromosomal level. It remains controversial whether genes expressed in a given tissue are randomly distributed throughout a given animal genome, or instead resolve into clusters. Here we used microarray analysis to identify more than 1,700 genes that are expressed preferentially in each of 11 organs of the chordate Ciona intestinalis adult, and determined the location of these genes on the 14 pairs of Ciona chromosomes. In spite of extensive mapped gene analysis, we only confirmed small clusters containing two or three genes. Our result indicates that organ-specific genes are distributed rather evenly all over chromosomes, suggesting that the notion of clustering of organ-specific genes in animal genomes is not generally applicable to this chordate.     

Efficient targeted mutagenesis of the chordate Ciona intestinalis genome with zinc-finger nucleases

Zinc‐finger nucleases (ZFNs) are engineered nucleases that induce DNA double‐strand breaks (DSBs) at target sequences. They have been used as tools for generating targeted mutations in the genomes of multiple organisms in both animals and plants. The DSB induced by ZFNs is repaired by non‐homologous end joining (NHEJ) or by homologous recombination (HR) mechanisms. Non‐homologous end joining induces some errors because it is independent of a reference DNA sequence. Through the NHEJ mechanism, ZFNs generate insertional or deletional mutations at the target sequence. We examined the usability, specificity and toxicity of ZFNs in the basal chordate Ciona intestinalis. As the target of ZFNs, we chose an enhanced green fluorescent protein (EGFP) gene artificially inserted in the C. intestinalis genome because this locus is neutral for the development and growth of C. intestinalis, and the efficiency of mutagenesis with ZFNs can thus be determined without any bias. We introduced EGFP ‐ZFN mRNAs into the embryos of an EGFP ‐transgenic line and observed the mutation frequency in the target site of EGFP . We also examined the effects of the EGFP ‐ZFNs at off‐target sites resembling the EGFP target sequence in the C. intestinalis genome in order to examine the specificity of ZFNs. We further investigated the influence of ZFNs on embryogenesis, and showed that adequate amounts of ZFNs, which do not disrupt embryogenesis, can efficiently induce mutations on the on‐target site with less effect on the off‐target sites. This suggests that target mutagenesis with ZFNs will be a powerful technique in C. intestinalis.     

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

Background

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

Results

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

Conclusion

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

The formation and positioning of cilia in Ciona intestinalis embryos in relation to the generation and evolution of chordate left-right asymmetry

In the early mouse embryo monocilia on the ventral node rotate to generate a leftward flow of fluid. This nodal flow is essential for the left-sided expression of nodal and pitx2, and for subsequent asymmetric organ patterning. Equivalent left fluid flow has been identified in other vertebrates, including Xenopus and zebrafish, indicating it is an ancient vertebrate mechanism. Asymmetric nodal and Pitx expression have also been identified in several invertebrates, including the vertebrates' nearest relatives, the urochordates. However whether cilia regulate this asymmetric gene expression remains unknown, and previous studies in urochordates have not identified any cilia prior to the larval stage, when asymmetry is already long established. Here we use Scanning and Transmission Electron Microscopy and immunofluorescence to investigate cilia in the urochordate Ciona intestinalis. We show that single cilia are transiently present on each ectoderm cell of the late neurula/early tailbud stage embryo, a time point just before onset of asymmetric nodal expression. Mapping the position of each cilium on these cells shows they are posteriorly positioned, something also described for mouse node cilia. The C. intestinalis cilia have a 9+0 ring ultrastructure, however we find no evidence of structures associated with motility such as dynein arms, radial spokes or nexin. Furthermore the 9+0 ring structure becomes disorganised immediately after the cilia have exited the cell, indicative of cilia which are not capable of motility. Our results indicate that although cilia are present prior to molecular asymmetries, they are not motile and hence cannot be operating in the same way as the flow-generating cilia of the vertebrate node. We conclude that the cilia may have a role in the development of C. intestinalis left-right asymmetry but that this would have to be in a sensory capacity, perhaps as mechanosensors as hypothesised in two-cilia physical models of vertebrate cilia-driven asymmetry.     

Stomodeal and neurohypophysial placodes in Ciona intestinalis: insights into the origin of the pituitary gland

The ascidian larva has a central nervous system which shares basic characteristics with craniates, such as tripartite organisation and many developmental genes. One difference, at metamorphosis, is that this chordate-like nervous system regresses and the adult's neural complex, composed of the cerebral ganglion and associated neural gland, forms. It is known that neural complex differentiation involves two ectodermal structures, the neurohypophysial duct, derived from the embryonic neural tube, and the stomodeum, i.e. the rudiment of the oral siphon; nevertheless, their precise role remains to be clarified. We have shown that in Ciona intestinalis, the neural complex primordium is the neurohypophysial duct, which in the early larva is a short tube, blind anteriorly, with its lumen in continuity with that of the central nervous system, i.e. the sensory vesicle. The tube grows forwards and fuses with the posterior wall of the stomodeum, a dorsal ectodermal invagination of the larva. The duct then loses posterior communication with the sensory vesicle and begins to grow on the roof of the vesicle itself. The neurohypophysial duct differentiates into the neural gland rudiment; its dorsal wall begins to proliferate neuroblasts, which migrate and converge to build up the cerebral ganglion. The most anterior part of the neural gland organizes into the ciliated duct and funnel, whereas the most posterior part elongates and gives rise to the dorsal strand. The hypothesis that the neurohypophysial duct/stomodeum complex possesses cell populations homologous to the craniate olfactory and adenohypophysial placodes and hypothalamus is discussed.     

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.     

Loss of ancestral genes in the genomic evolution of Ciona intestinalis

Comparison of the predicted protein sets encoded by the complete genomes of two vertebrates (human and pufferfish), the urochordate Ciona intestinalis, three nonchordate animals, and two fungi were used to reconstruct a set of gene families present in the common ancestor of chordates. These ancestral families were much more likely to be lost in Ciona than in either vertebrate. In addition, of 256 duplicate gene pairs that arose by duplication prior to the most recent common ancestor of vertebrates and insects, one of the duplicate genes was four times as likely to be lost in Ciona as in the vertebrates. These results show that the genome of Ciona is not representative of the ancestral chordate genome with respect to gene content but rather shows derived features that may reflect adaptation of the specific ecological niche of urochordates.