The draft genome of Ciona intestinalis: insights into chordate and vertebrate origins

Proteome-wide classification and identification of mammalian-type GPCRs by binary topology pattern

G protein-coupled receptors (GPCRs), a large eukaryotic protein family, have proved difficult to comprehensively detect and functionally identify by homology searches and domain detection, because they are highly divergent and their sequences share strikingly little similarity. Transmembrane (TM) topology pattern analysis has been used to classify TM proteins, and such patterns are conserved within GPCRs of similar function. Here, we developed a stepwise binary topology pattern (BTP) method for GPCR classification and identification and used it to identify and classify mammalian-type GPCRs in the genomes of 10 different eukaryotic species. A binary topology pattern was obtained for each functional class or group by assigning binary loop threshold lengths of "0" (short loop) or "1" (long loop). The GPCR-classification ability of the BTP method had quite high accuracies for classifying GPCR functions at the class level (Classes A, B, C, Frizzled/Smoothened, Non-GPCR, based on the GPCRDB classification scheme), with many classes being classified with 100% accuracy. Sufficiently high accuracies were also maintained at the functional group level, 0.945 over 15 functional groups. Proteome-wide mammalian-type GPCR searches in 10 eukaryotic genomes (H. sapiens, M. musculus, F. rubripes, C. intestinalis, A. thaliana, D. melanogaster, A. gambiae, C. elegans, P. falciparum, S. cerevisiae) using the BTP method showed much higher classification/identification in non-mammalian genomes than typical BLAST searches, in which a higher number of sequences were classified as Non-GPCR. This stepwise BTP method should prove useful for the identification and functional classification of GPCRs from the genomes of a wide range of species.

The revolution of the biology of the genome

Sequence data of entire eukaryotic genomes and their detailed comparison have provided new evidence on genome evolution. The major mechanisms involved in the increase of genome sizes are polyploidization and gene duplication. Subsequent gene silencing or mutations, preferentially in regulatory sequences of genes, modify the genome and permit the development of genes with new properties. Mechanisms such as lateral gene transfer, exon shuffling or the creation of new genes by transposition contribute to the evolution of a genome, but remain of relatively restricted relevance. Mechanisms to decrease genome sizes and, in particular, to remove specific DNA sequences, such as blocks of satellite DNAs, appear to involve the action of RNA interference (RNAi). RNAi mechanisms have been proven to be involved in chromatin packaging related with gene inactivation as well as in DNA excision during the macronucleus development in ciliates.

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.

A functional cellulose synthase from ascidian epidermis

Among animals, urochordates (e.g., ascidians) are unique in their ability to biosynthesize cellulose. In ascidians cellulose is synthesized in the epidermis and incorporated into a protective coat know as the tunic. A putative cellulose synthase-like gene was first identified in the genome sequences of the ascidian Ciona intestinalis. We describe here a cellulose synthase gene from the ascidian Ciona savignyi that is expressed in the epidermis. The predicted C. savignyi cellulose synthase amino acid sequence showed conserved features found in all cellulose synthases, including plants, but was most similar to cellulose synthases from bacteria, fungi, and Dictyostelium discoidium. However, unlike other known cellulose synthases, the predicted C. savignyi polypeptide has a degenerate cellulase-like region near the carboxyl-terminal end. An expression construct carrying the C. savignyi cDNA was found to restore cellulose biosynthesis to a cellulose synthase (CelA) minus mutant of Agrobacterium tumefaciens, showing that the predicted protein has cellulose synthase activity. The lack of cellulose biosynthesis in all other groups of metazoans and the similarity of the C. savignyi cellulose synthase to enzymes from cellulose-producing organisms support the hypothesis that the urochordates acquired the cellulose biosynthetic pathway by horizontal transfer.

A genome-wide survey of the genes for planar polarity signaling or convergent extension-related genes in Ciona intestinalis and phylogenetic comparisons of evolutionary conserved signaling components

Non-canonical Wnt signals similar to planar cell polarity (PCP) signaling in the fly control convergent extension (CE) of the dorsal mesoderm during gastrulation in vertebrates. Using the Ciona complete genome sequence and EST sequence data, we present here an initial and exhaustive search in non-vertebrate chordates, Ciona intestinalis for the family members as well as homologs or orthologs that are involved in PCP/CE signaling cascades. We clarified 7 cardinal gene families, including the MAPK, STE20 group kinase, Rho small GTPase, STAT, Glypican, Fz and Wnt gene families, as well as gene homologs or orthologs for known PCP/CE signaling components with their phylogenetic nature. As a result, we characterized 62 Ciona component genes. Among them, 59 genes were novel and functional genes which were supported by EST expressions and 15 genes belonged to PCP/CE component orthologs of other organisms or common ancestor genes. Moreover, from the phylogenetic point of view, we compared these components genome-widely with the PCP signaling components of fly and the CE signaling components of vertebrates. We then discovered not only that ascidians contain the basic ancestral signaling pathway components in chordates but also that several signaling components have not found in ascidian, indicating that ascidian CE pathway might have several gaps from vertebrate CE pathway. The present study provides an initial step for the subsequent analysis of CE in the non-vertebrate chordates, ascidians. In addition, this phylogenetic approach will help to facilitate understanding of the relationship between fly PCP signaling and the vertebrate CE pathway.

Evolution of the multifunctional protein tyrosine phosphatase family

The protein tyrosine phosphatase (PTP) family plays a central role in signal transduction pathways by controlling the phosphorylation state of serine, threonine, and tyrosine residues. PTPs can be divided into dual specificity phosphatases and the classical PTPs, which can comprise of one or two phosphatase domains. We studied amino acid substitutions at functional sites in the phosphatase domain and identified putative noncatalytic phosphatase domains in all subclasses of the PTP family. The presence of inactive phosphatase domains in all subclasses indicates that they were invented multiple times in evolution. Depending on the domain composition, loss of catalytic activity can result in different consequences for the function of the protein. Inactive single-domain phosphatases can still specifically bind substrate and protect it from dephosphorylation by other phosphatases. The inactive domains of tandem phosphatases can be further subdivided. The first class is more conserved, still able to bind phosphorylated tyrosine residues and might recruit multiphosphorylated substrates for the adjacent active domain. The second has accumulated several variable amino acid substitutions in the catalytic center, indicating a complete loss of tyrosine-binding capabilities. To study the impact of substitutions in the catalytic center to the evolution of the whole domain, we examined the evolutionary rates for each individual site and compared them between the classes. This analysis revealed a release of evolutionary constraint for multiple sites surrounding the catalytic center only in the second class, emphasizing its difference in function compared with the first class. Furthermore, we found a region of higher conservation common to both domain classes, suggesting a new regulatory center. We discuss the influence of evolutionary forces on the development of the phosphatase domain, which has led to additional functions, such as the specific protection of phosphorylated tyrosine residues, substrate recruitment, and regulation of the catalytic activity of adjacent domains.

Human SLITRK family genes: genomic organization and expression profiling in normal brain and brain tumor tissue

Slitrk family proteins are characterized as integral membrane proteins that have two leucine-rich repeat (LRR) domains and a carboxy-terminal domain that is partially similar to trk neurotrophin receptor proteins. The LRR domains are similar to those of slit proteins. In a previous study, we showed that mouse Slitrk genes are expressed predominantly in neural tissue and have neurite-modulating activity in cultured neuronal cells. Their expression profiles as well as their functions vary among the family members. In this paper, we characterized the human SLITRK1, SLITRK2, SLITRK3, SLITRK4, SLITRK5, and SLITRK6 genes. The six genes are located in three clusters, on 3q, 13q, and Xq, respectively. Their expression was detected mainly in the brain, but the expression profile of each SLITRK was unique. SLITRK expression was also investigated in various types of brain tumor tissue. The results showed that all SLITRK genes are differentially expressed in brain tumors, including astrocytoma, oligodendroglioma, glioblastoma, medulloblastoma, and supratentorial primitive neuroectodermal tumor (PNET). Particularly interesting findings were that SLITRK3 expression was enhanced in tissue from several different types of tumors and SLITRK6 expression was highly selective. These results suggest that the human SLITRK genes are useful molecular indicators of brain tumor properties.

Generation, annotation, evolutionary analysis, and database integration of 20,000 unique sea urchin EST clusters

Together with the hemichordates, sea urchins represent basal groups of nonchordate invertebrate deuterostomes that occupy a key position in bilaterian evolution. Because sea urchin embryos are also amenable to functional studies, the sea urchin system has emerged as one of the leading models for the analysis of the function of genomic regulatory networks that control development. We have analyzed a total of 107,283 cDNA clones of libraries that span the development of the sea urchin Strongylocentrotus purpuratus. Normalization by oligonucleotide fingerprinting, EST sequencing and sequence clustering resulted in an EST catalog comprised of 20,000 unique genes or gene fragments. Around 7000 of the unique EST consensus sequences were associated with molecular and developmental functions. Phylogenetic comparison of the identified genes to the genome of the urochordate Ciona intestinalis indicate that at least one quarter of the genes thought to be chordate specific were already present at the base of deuterostome evolution. Comparison of the number of gene copies in sea urchins to those in chordates and vertebrates indicates that the sea urchin genome has not undergone extensive gene or complete genome duplications. The established unique gene set represents an essential tool for the annotation and assembly of the forthcoming sea urchin genome sequence. All cDNA clones and filters of all analyzed libraries are available from the resource center of the German genome project at http://www.rzpd.de.

A comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes

We examined functional and evolutionary patterns in the recently constructed set of 5,873 clusters of predicted orthologs from seven eukaryotic genomes. The analysis reveals a conserved core of largely essential eukaryotic genes as well as major diversification and innovation associated with evolution of eukaryotic genomes.

Background

Sequencing the genomes of multiple, taxonomically diverse eukaryotes enables in-depth comparative-genomic analysis which is expected to help in reconstructing ancestral eukaryotic genomes and major events in eukaryotic evolution and in making functional predictions for currently uncharacterized conserved genes.

Results

We examined functional and evolutionary patterns in the recently constructed set of 5,873 clusters of predicted orthologs (eukaryotic orthologous groups or KOGs) from seven eukaryotic genomes: Caenorhabditis elegans, Drosophila melanogaster, Homo sapiens, Arabidopsis thaliana, Saccharomyces cerevisiae, Schizosaccharomyces pombe and Encephalitozoon cuniculi. Conservation of KOGs through the phyletic range of eukaryotes strongly correlates with their functions and with the effect of gene knockout on the organism's viability. The approximately 40% of KOGs that are represented in six or seven species are enriched in proteins responsible for housekeeping functions, particularly translation and RNA processing. These conserved KOGs are often essential for survival and might approximate the minimal set of essential eukaryotic genes. The 131 single-member, pan-eukaryotic KOGs we identified were examined in detail. For around 20 that remained uncharacterized, functions were predicted by in-depth sequence analysis and examination of genomic context. Nearly all these proteins are subunits of known or predicted multiprotein complexes, in agreement with the balance hypothesis of evolution of gene copy number. Other KOGs show a variety of phyletic patterns, which points to major contributions of lineage-specific gene loss and the 'invention' of genes new to eukaryotic evolution. Examination of the sets of KOGs lost in individual lineages reveals co-elimination of functionally connected genes. Parsimonious scenarios of eukaryotic genome evolution and gene sets for ancestral eukaryotic forms were reconstructed. The gene set of the last common ancestor of the crown group consists of 3,413 KOGs and largely includes proteins involved in genome replication and expression, and central metabolism. Only 44% of the KOGs, mostly from the reconstructed gene set of the last common ancestor of the crown group, have detectable homologs in prokaryotes; the remainder apparently evolved via duplication with divergence and invention of new genes.

Conclusions

The KOG analysis reveals a conserved core of largely essential eukaryotic genes as well as major diversification and innovation associated with evolution of eukaryotic genomes. The results provide quantitative support for major trends of eukaryotic evolution noticed previously at the qualitative level and a basis for detailed reconstruction of evolution of eukaryotic genomes and biology of ancestral forms.

Kowalevsky, comparative evolutionary embryology, and the intellectual lineage of evo-devo

Alexander Kowalevsky was one of the most significant 19th century biologists working at the intersection of evolution and embryology. The reinstatement of the Alexander Kowalevsky Medal by the St. Petersburg Society of Naturalists for outstanding contributions to understanding evolutionary relationships in the animal kingdom, evolutionary developmental biology, and comparative zoology is timely now that Evo-devo has emerged as a major research discipline in contemporary biology. Consideration of the intellectual lineage of comparative evolutionary embryology explicitly forces a reconsideration of some current conceptions of the modern emergence of Evo-devo, which has tended to exist in the shadow of experimental embryology throughout the 20th century, especially with respect to the recent success of developmental biology and developmental genetics. In particular we advocate a sharper distinction between the heritage of problems and the heritage of tools for contemporary Evo-devo. We provide brief overviews of the work of N. J. Berrill and D. T. Anderson to illustrate comparative evolutionary embryology in the 20th century, which provides an appropriate contextualization for a conceptual review of our research on the sea urchin genus Heliocidaris over the past two decades. We conclude that keeping research questions rather than experimental capabilities at the forefront of Evo-devo may be an antidote to any repeat of the stagnation experienced by the first group of evolutionary developmental biologists over one hundred years ago and acknowledges Kowalevsky's legacy in evolutionary embryology.

Epithelial cell cultures from Botryllus schlosseri palleal buds: accomplishments and challenges

This study focuses on recent improvement in epithelial monolayer cultures originating from whole extirpated Botryllus schlosseri (Urochordata) buds. Buds (n = 2,000) were taken at different ('A' to 'D') blastogenic stages. We tested the suitability of 35 combinations of various substrates and media on attachment, cell spread, epithelial growth frequencies and on monolayer lifespans. Under favorable conditions, cultured buds at blastogenic stages 'B' to 'D' (but not stage 'A') started to attach to the substrates following a 3-day transient period that leads to formation of spheres and attached monolayers. Substrate type is important for the attachment and the development of monolayers. Under various culture conditions, some of stages 'B' and 'C' buds develop (3-20 days) one or more large (1 mm diameter) spheres. Stage 'D' buds develop monolayers (up to 20% of buds) without going through a sphere phase. Neither spheres nor attached monolayers of epithelium were observed in stage 'A' bud cultures. Spheres grew at a rate of 60 microm in diameter per day using specific medium types and did not attach unless the appropriate substrate was present. When attached, epithelial monolayers expanded at a rate of 200 microm in diameter per day, for 3-15 days, and subsequently detached and died. Sixteen types of media were tested. Medium and substrate combinations were found to determine epithelial lifespan. These results revealed significant improvements in the culture of epithelial monolayers from Botryllus palleal buds. However, an early senescence of the developed epithelial sheets (up to two weeks from onset of appearance) may indicate an internal ageing clock that should be taken into consideration in future approaches.

Exclusion of repetitive DNA elements from gnathostomeHox clusters

Despite their homology and analogous function, the Hox gene clusters of vertebrates and invertebrates are subject to different constraints on their structural organization. This is demonstrated by a drastically different distribution of repetitive DNA elements in the Hox cluster regions. While gnathostomes have a strong tendency to exclude repetitive DNA elements from the inside of their Hox clusters, no such trend can be detected in the Hox gene clusters of protostomes. Repeats "invade" the gnathostome Hox clusters from the 5' and 3' ends while the core of the clusters remains virtually free of repetitive DNA. This invasion appears to be correlated with relaxed constraints associated with gene loss after cluster duplications.

Morpholino-based gene knockdown screen of novel genes with developmental function inCiona intestinalis

In the present study, we conducted an extensive analysis to identify novel genes with developmental function among Ciona intestinalis genes discovered by cDNA projects. Translation of a total of 200 genes expressed during embryogenesis was suppressed by using specific morpholino antisense oligonucleotides. Suppression of the translation of any of 40 genes (one-fifth of the genes tested) was thereby shown to cause specific embryonic defects. Most of these genes have counterpart(s) in mouse and human, suggesting that the present approach will be useful for identifying candidate genes essential for the development of vertebrates. Suppression of translation of 14 of these 40 genes resulted in the `disorganized body plan' phenotype characterized by gross morphological abnormalities caused by early defects in embryogenesis. These genes encode zinc-finger, transmembrane or Pbx homeodomain proteins. The morphological features of larvae of this phenotypic class varied according to the gene suppressed, suggesting that a distinct developmental event such as tissue specification or cell cycle progression was affected in each type of larva. Suppression of the remaining 26 genes resulted in the `abnormal tail'phenotype. Some of these genes encode proteins with known functional structures such as Zn-finger and HLH motifs. Twelve genes among them are especially interesting, because their suppression produced defects in the nervous system, as demonstrated by the loss of the sensory pigment cells or palps of the adhesive organ in the knockdown larvae. These results suggest that screening for developmental genes by the reverse genetic approach in Ciona intestinalis embryos is effective for identifying novel genes with developmental functions required for the development of chordates.

Mechanisms of germ cell specification across the metazoans: epigenesis and preformation

Germ cells play a unique role in gamete production, heredity and evolution. Therefore, to understand the mechanisms that specify germ cells is a central challenge in developmental and evolutionary biology. Data from model organisms show that germ cells can be specified either by maternally inherited determinants (preformation) or by inductive signals (epigenesis). Here we review existing data on 28 metazoan phyla, which indicate that although preformation is seen in most model organisms, it is actually the less prevalent mode of germ cell specification, and that epigenetic germ cell specification may be ancestral to the Metazoa.

Hox cluster duplications and the opportunity for evolutionary novelties

Hox genes play a key role in animal body plan development. These genes tend to occur in tightly linked clusters in the genome. Vertebrates and invertebrates differ in their Hox cluster number, with vertebrates having multiple clusters and invertebrates usually having only one. Recent evidence shows that vertebrate Hox clusters are structurally more constrained than invertebrate Hox clusters; they exclude transposable elements, do not undergo tandem duplications, and conserve their intergenic distances and gene order. These constraints are only relaxed after a cluster duplication. In contrast, invertebrate Hox clusters are structurally more plastic; tandem duplications are common, the linkage of Hox genes can change quickly, or they can lose their structural integrity completely. We propose that the constraints on vertebrate Hox cluster structure lead to an association between the retention of duplicated Hox clusters and adaptive radiations. After a duplication the constraints on Hox cluster structure are temporarily lifted, which opens a window of evolvability for the Hox clusters. If this window of evolvability coincides with an adaptive radiation, chances are that a modified Hox cluster becomes recruited in an evolutionary novelty and then both copies of duplicated Hox clusters are retained.

Fugu ESTs: new resources for transcription analysis and genome annotation

The draft Fugu rubripes genome was released in 2002, at which time relatively few cDNAs were available to aid in the annotation of genes. The data presented here describe the sequencing and analysis of 24,398 expressed sequence tags (ESTs) generated from 15 different adult and juvenile Fugu tissues, 74% of which matched protein database entries. Analysis of the EST data compared with the Fugu genome data predicts that approximately 10,116 gene tags have been generated, covering almost one-third of Fugu predicted genes. This represents a remarkable economy of effort. Comparison with the Washington University zebrafish EST assemblies indicates strong conservation within fish species, but significant differences remain. This potentially represents divergence of sequence in the 5' terminal exons and UTRs between these two fish species, although clearly, complete EST data sets are not available for either species. This project provides new Fugu resources, and the analysis adds significant weight to the argument that EST programs remain an essential resource for genome exploitation and annotation. This is particularly timely with the increasing availability of draft genome sequence from different organisms and the mounting emphasis on gene function and regulation.

Characterizing embryonic gene expression patterns in the mouse using nonredundant sequence-based selection

This article investigates the expression patterns of 160 genes that are expressed during early mouse development. The cDNAs were isolated from 7.5 d postcoitum (dpc) endoderm, a region that comprises visceral endoderm (VE), definitive endoderm, and the node-tissues that are required for the initial steps of axial specification and tissue patterning in the mouse. To avoid examining the same gene more than once, and to exclude potentially ubiquitously expressed housekeeping genes, cDNA sequence was derived from 1978 clones of the Endoderm library. These yielded 1440 distinct cDNAs, of which 123 proved to be novel in the mouse. In situ hybridization analysis was carried out on 160 of the cDNAs, and of these, 29 (18%) proved to have restricted expression patterns.

Structure, expression, and cluster organization of genes encoding gonadotropin-releasing hormone receptors found in the neural complex of the ascidian Ciona intestinalis

Gonadotropin-releasing hormone (GnRH) controls reproductive functions in vertebrates. In the present study, two distinct homologues (Ci-GnRHR1 and Ci-GnRHR2) of the vertebrate GnRH receptor (GnRHR) were identified by cDNA cloning from the neural complex of the ascidian Ciona intestinalis. The Ci-GnRHR1 and Ci-GnRHR2 genes are closely linked to each other and form a cluster with another GnRHR-like gene in the genome. Ci-GnRHR1 and Ci-GnRHR2 are more closely related to vertebrate GnRHRs than to the GnRHR-related protein of Drosophila. The phylogenetic analysis indicates that vertebrates and urochordates independently acquired multiple GnRHRs by gene duplications that occurred during the evolution of each lineage. A voltage clamp of Xenopus oocytes injected with synthetic Ci-GnRHR1 mRNA revealed inward currents in response to an ascidian form of GnRH, suggesting that Ci-GnRHR1 is a bona-fide GnRHR. Expression patterns of Ci-GnRHR1 and Ci-GnRHR2 suggest that a GnRH signaling system is involved in regulation of neuronal and reproductive processes as well as in other physiological functions in ascidians.

Turning the clock back on ancient genome duplication

Complete genome sequence data led rapidly to the conclusion that ancient genome duplications had shaped the genomes of the model organisms Saccharomyces cerevisiae and Arabidopsis thaliana. Recent contributions have gone on to refine date estimates for these duplications and, in the case of Arabidopsis, to infer additional, more ancient, rounds of duplication by reconstructing gene order before the most recent duplication event. It is becoming widely accepted that an ancient duplication occurred before the radiation of the ray-finned fish. However, despite methodological advances and the availability of complete genome sequence data the debate over whether very ancient genome duplications have occurred early in the vertebrate lineage has not yet been fully resolved.

A morphological and genetic characterization of metamorphosis in the ascidian Boltenia villosa

Ascidian metamorphosis is a critical life history stage for exploring chordate evolution and conserved chordate developmental signaling pathways. The vast majority of research on ascidian development has been focused on embryogenesis. Thus there is still little known about the development of ascidian post-larval structures, including differentiation of the chordate pharyngeal gill slits and endostyle along with the heart, blood cells and gut. In this paper, we present our research on metamorphosis in the solitary ascidian Boltenia villosa. Through careful analysis of phalloidin staining in young juveniles, we have discerned a highly coordinated series of developmental events underlying the differentiation of the gut and body wall musculature. Additionally, we have employed subtractive hybridizations to isolate genes that are differentially transcribed during Boltenia metamorphosis. Some of these genes are expressed throughout ascidian development and some appear to be uniquely expressed during metamorphosis. Here we characterize several transcripts with potential developmental functions and discuss their possible roles in the differentiation of adult structures during solitary ascidian metamorphosis.

Eukaryotic domain evolution inferred from genome comparisons

Comparative analyses of eukaryotic genomes are providing insights into the mode and tempo of domain family evolution. Gene duplication, the source of family expansion, far exceeds the rate of emergence of domains from non-coding sequence, and the rate of recruitment of domains into novel architectures. Domain families that appear to be restricted to certain lineages are likely to be the result of gene duplication, coupled with rapid sequence diversification. If such families are evidence of past adaptation, then their functions must relate to the underlying mechanism of selection: competition among organisms.

EST Analysis of the Cnidarian Acropora millepora Reveals Extensive Gene Loss and Rapid Sequence Divergence in the Model Invertebrates

A significant proportion of mammalian genes are not represented in the genomes of Drosophila, Caenorhabditis or Saccharomyces, and many of these are assumed to have been vertebrate innovations. To test this assumption, we conducted a preliminary EST project on the anthozoan cnidarian, Acropora millepora, a basal metazoan. More than 10% of the Acropora ESTs with strong metazoan matches to the databases had clear human homologs but were not represented in the Drosophila or Caenorhabditis genomes; this category includes a surprising diversity of transcription factors and metabolic proteins that were previously assumed to be restricted to vertebrates. Consistent with higher rates of divergence in the model invertebrates, three-way comparisons show that most Acropora ESTs match human sequences much more strongly than they do any Drosophila or Caenorhabditis sequence. Gene loss has thus been much more extensive in the model invertebrate lineages than previously assumed and, as a consequence, some genes formerly thought to be vertebrate inventions must have been present in the common metazoan ancestor. The complexity of the Acropora genome is paradoxical, given that this organism contains apparently few tissue types and the simplest extant nervous system consisting of a morphologically homogeneous nerve net.

Molecular evolutionary analysis of the widespread piggyBac transposon family and related "domesticated" sequences

piggyBac is a short inverted-repeat-type DNA transposable element originally isolated from the genome of the moth Trichoplusia ni. It is currently the gene vector of choice for the transformation of various insect species. A few sequences with similarity to piggyBac have previously been identified from organisms such as humans ( Looper), the pufferfish Takifugu rubripes ( Pigibaku), Xenopus ( Tx), Daphnia ( Pokey), and the Oriental fruit fly Bactrocera dorsalis. We have now identified 50 piggyBac-like sequences from publicly available genome sequences and expressed sequence tags (ESTs). This survey allows the first comparative examination of the distinctive piggyBac transposase, suggesting that it might contain a highly divergent DDD domain, comparable to the widespread DDE domain found in many DNA transposases and retroviral integrases which consists of two absolutely conserved aspartic acids separated by about 70 amino acids with a highly conserved glutamic acid about 35 amino acids further away. Many piggyBac-like sequences were found in the genomes of a phylogenetically diverse range of organisms including fungi, plants, insects, crustaceans, urochordates, amphibians, fishes and mammals. Also, several instances of "domestication" of the piggyBac transposase sequence by the host genome for cellular functions were identified. Novel members of the piggyBac family may be useful in genetic engineering of many organisms.

Hagfish intestinal antimicrobial peptides are ancient cathelicidins

Three potent broad-spectrum antimicrobial peptides (HFIAP-1, -2, and -3) isolated from intestinal tissues of Myxine glutinosa (Atlantic hagfish) are identified as ancient members of the cathelicidin family of antimicrobial peptides, hitherto known only from mammals. In situ hybridization reveals that HFIAPs are produced in nests of myeloid cells within the loose connective tissue of the gut wall, a tissue reminiscent of both gut-associated lymphoid tissue (GALT) and vertebrate spleen. We suggest that this tissue organization provides local defense of the hagfish gastrointestinal tract via innate immunity and possibly served as the architectural plan upon which the adaptive immune system evolved.

Neural Tissue in Ascidian Embryos Is Induced by FGF9/16/20, Acting via a Combination of Maternal GATA and Ets Transcription Factors

In chordates, formation of neural tissue from ectodermal cells requires an induction. The molecular nature of the inducer remains controversial in vertebrates. Here, using the early neural marker Otx as an entry point, we dissected the neural induction pathway in the simple embryos of Ciona intestinalis. We first isolated the regulatory element driving Otx expression in the prospective neural tissue, showed that this element directly responds to FGF signaling and that FGF9/16/20 acts as an endogenous neural inducer. Binding site analysis and gene loss of function established that FGF9/16/20 induces neural tissue in the ectoderm via a synergy between two maternal response factors. Ets1/2 mediates general FGF responsiveness, while the restricted activity of GATAa targets the neural program to the ectoderm. Thus, our study identifies an endogenous FGF neural inducer and its early downstream gene cascade. It also reveals a role for GATA factors in FGF signaling.

Novel vanadium-binding proteins (vanabins) identified in cDNA libraries and the genome of the ascidian Ciona intestinalis

Ascidians, especially those belonging to the suborder Phlebobranchia, can accumulate high levels of vanadium. Vanadium-binding proteins (vanabins) were first isolated from a vanadium-accumulating ascidian, Ascidia sydneiensis samea, and then the vanabins were cloned, their expression was studied, and metal-binding assays were conducted. In order to unravel the mechanism of vanadium accumulation, we searched for vanabin-like genes in other animals, including other ascidians. A database search revealed five groups of cDNAs that encoded vanabin-like proteins in another ascidian, Ciona intestinalis. The genes encoding C. intestinalis vanabins, CiVanabin1 to CiVanabin5, were clustered in an 8.4-kb genomic region. The direction of transcription of each gene was identical and each gene had a single intron. All the C. intestinalis vanabins were cysteine-rich, and the repetitive pattern of cysteines closely resembled that of A. sydneiensis samea vanabins. Using immobilized metal ion affinity chromatography (IMAC), we found that a recombinant protein of at least one of the C. intestinalis vanabins (CiVanabin5) bound to vanadium(IV) ions.

Reevaluation of human cytomegalovirus coding potential

The Bio-Dictionary-based Gene Finder was used to reassess the coding potential of the AD169 laboratory strain of human cytomegalovirus and sequences in the Toledo strain that are missing in the laboratory strain of the virus. The gene-finder algorithm assesses the potential of an ORF to encode a protein based on matches to a database of amino acid patterns derived from a large collection of proteins. The algorithm was used to score all human cytomegalovirus ORFs with the potential to encode polypeptides >/=50 aa in length. As a further test for functionality, the genomes of the chimpanzee, rhesus, and murine cytomegaloviruses were searched for orthologues of the predicted human cytomegalovirus ORFs. The analysis indicates that 37 previously annotated ORFs ought to be discarded, and at least nine previously unrecognized ORFs with relatively strong coding potential should be added. Thus, the human cytomegalovirus genome appears to contain approximately 192 unique ORFs with the potential to encode a protein. Support for several of the predictions of our in silico analysis was obtained by sequencing several domains within a clinical isolate of human cytomegalovirus.

Terminal deoxynucleotidyl transferases from elasmobranchs reveal structural conservation within vertebrates

The DNA polymerase (pol) X family is an ancient group of enzymes that function in DNA replication and repair (pol beta), translesion synthesis (pol lambda and pol micro) and terminal addition of non-templated nucleotides. This latter terminal deoxynucleotidyl transferase (TdT) activity performs the unique function of providing diversity at coding joins of immunoglobulin and T-cell receptor genes. The first isolated full-length TdT genes from shark and skate are reported here. Comparisons with the three-dimensional structure of mouse TdT indicate structural similarity with elasmobranch orthologues that supports both a template-independent mode of replication and a lack of strong nucleotide bias. The vertebrate TdTs appear more closely related to pol micro and fungal polymerases than to pol lambda and pol beta. Thus, unlike other molecules of adaptive immunity, TdT is a member of an ancient gene family with a clear gene phylogeny and a high degree of similarity, which implies the existence of TdT ancestors in jawless fishes and invertebrates.

The non-LTR retrotransposons in Ciona intestinalis: new insights into the evolution of chordate genomes

In silico and experimental approaches have been used to identify the non-long terminal repeat retrotransposons of the urochordate Ciona intestinalis providing valuable data for understanding the evolution of early chordate genomes.

Background

Non-long terminal repeat (non-LTR) retrotransposons have contributed to shaping the structure and function of genomes. In silico and experimental approaches have been used to identify the non-LTR elements of the urochordate Ciona intestinalis. Knowledge of the types and abundance of non-LTR elements in urochordates is a key step in understanding their contribution to the structure and function of vertebrate genomes.

Results

Consensus elements phylogenetically related to the I, LINE1, LINE2, LOA and R2 elements of the 14 eukaryotic non-LTR clades are described from C. intestinalis. The ascidian elements showed conservation of both the reverse transcriptase coding sequence and the overall structural organization seen in each clade. The apurinic/apyrimidinic endonuclease and nucleic-acid-binding domains encoded upstream of the reverse transcriptase, and the RNase H and the restriction enzyme-like endonuclease motifs encoded downstream of the reverse transcriptase were identified in the corresponding Ciona families.

Conclusions

The genome of C. intestinalis harbors representatives of at least five clades of non-LTR retrotransposons. The copy number per haploid genome of each element is low, less than 100, far below the values reported for vertebrate counterparts but within the range for protostomes. Genomic and sequence analysis shows that the ascidian non-LTR elements are unmethylated and flanked by genomic segments with a gene density lower than average for the genome. The analysis provides valuable data for understanding the evolution of early chordate genomes and enlarges the view on the distribution of the non-LTR retrotransposons in eukaryotes.

Constructing a feedback loop with circadian clock molecules from the silkmoth, Antheraea pernyi

Circadian clocks are important regulators of behavior and physiology. The circadian clock of Drosophila depends on an autoinhibitory feedback loop involving dCLOCK, CYCLE (also called dBMAL, for Drosophila brain and muscle ARNT-like protein), dPERIOD, and dTIMELESS. Recent studies suggest that the clock mechanism in other insect species may differ strikingly from that of Drosophila. We cloned Clock, Bmal, and Timeless homologs (apClock, apBmal, and apTimeless) from the silkmoth Antheraea pernyi, from which a Period homolog (apPeriod) has already been cloned. In Schneider 2 (S2) cell culture assays, apCLOCK:apBMAL activates transcription through an E-box enhancer element found in the 5' region of the apPeriod gene. Furthermore, apPERIOD can robustly inhibit apCLOCK: apBMAL-mediated transactivation, and apTIMELESS can augment this inhibition. Thus, a complete feedback loop, resembling that found in Drosophila, can be constructed from silkmoth CLOCK, BMAL, PERIOD, and TIMELESS. Our results suggest that the circadian autoinhibitory feedback loop discovered in Drosophila is likely to be widespread among insects. However, whereas the transactivation domain in Drosophila lies in the C terminus of dCLOCK, in A. pernyi, it lies in the C terminus of apBMAL, which is highly conserved with the C termini of BMALs in other insects (except Drosophila) and in vertebrates. Our analysis sheds light on the molecular function and evolution of clock genes in the animal kingdom.

Transcriptome analysis of the acoelomate human parasite Schistosoma mansoni

Schistosoma mansoni is the primary causative agent of schistosomiasis, which affects 200 million individuals in 74 countries. We generated 163,000 expressed-sequence tags (ESTs) from normalized cDNA libraries from six selected developmental stages of the parasite, resulting in 31,000 assembled sequences and 92% sampling of an estimated 14,000 gene complement. By analyzing automated Gene Ontology assignments, we provide a detailed view of important S. mansoni biological systems, including characterization of metazoa-specific and eukarya-conserved genes. Phylogenetic analysis suggests an early divergence from other metazoa. The data set provides insights into the molecular mechanisms of tissue organization, development, signaling, sexual dimorphism, host interactions and immune evasion and identifies novel proteins to be investigated as vaccine candidates and potential drug targets.

Resurrecting the ancestral steroid receptor: ancient origin of estrogen signaling

Receptors for sex and adrenal steroid hormones are absent from fully sequenced invertebrate genomes and have not been recovered from other invertebrates. Here we report the isolation of an estrogen receptor ortholog from the mollusk Aplysia californica and the reconstruction, synthesis, and experimental characterization of functional domains of the ancestral protein from which all extant steroid receptors (SRs) evolved. Our findings indicate that SRs are extremely ancient and widespread, having diversified from a primordial gene before the origin of bilaterally symmetric animals, and that this ancient receptor had estrogen receptor-like functionality. This gene was lost in the lineage leading to arthropods and nematodes and became independent of hormone regulation in the Aplysia lineage.

Independent Hox-cluster duplications in lampreys

The analysis of the publicly available Hox gene sequences from the sea lamprey Petromyzon marinus provides evidence that the Hox clusters in lampreys and other vertebrate species arose from independent duplications. In particular, our analysis supports the hypothesis that the last common ancestor of agnathans and gnathostomes had only a single Hox cluster which was subsequently duplicated independently in the two lineages.

Phylogenetic analyses alone are insufficient to determine whether genome duplication(s) occurred during early vertebrate evolution

The widely accepted notion that two whole-genome duplications occurred during early vertebrate evolution (the 2R hypothesis) stems from the fact that vertebrates often possess several genes corresponding to a single invertebrate homolog. However the number of genes predicted by the Human Genome Project is less than twice as many as in the Drosophila melanogaster or Caenorhabditis elegans genomes. This ratio could be explained by two rounds of genome duplication followed by extensive gene loss, by a single genome duplication, by sequential local duplications, or by a combination of any of the above. The traditional method used to distinguish between these possibilities is to reconstruct the phylogenetic relationships of vertebrate genes to their invertebrate orthologs; ratios of invertebrate-to-vertebrate counterparts are then used to infer the number of gene duplication events. The lancelet, amphioxus, is the closest living invertebrate relative of the vertebrates, and unlike protostomes such as flies or nematodes, is therefore the most appropriate outgroup for understanding the genomic composition of the last common ancestor of all vertebrates. We analyzed the relationships of all available amphioxus genes to their vertebrate homologs. In most cases, one to three vertebrate genes are orthologous to each amphioxus gene (median number=2). Clearly this result, and those of previous studies using this approach, cannot distinguish between alternative scenarios of chordate genome expansion. We conclude that phylogenetic analyses alone will never be sufficient to determine whether genome duplication(s) occurred during early chordate evolution, and argue that a "phylogenomic" approach, which compares paralogous clusters of linked genes from complete amphioxus and human genome sequences, will be required if the pattern and process of early chordate genome evolution is ever to be reconstructed.

The COG database: an updated version includes eukaryotes

BACKGROUND

The availability of multiple, essentially complete genome sequences of prokaryotes and eukaryotes spurred both the demand and the opportunity for the construction of an evolutionary classification of genes from these genomes. Such a classification system based on orthologous relationships between genes appears to be a natural framework for comparative genomics and should facilitate both functional annotation of genomes and large-scale evolutionary studies.

RESULTS

We describe here a major update of the previously developed system for delineation of Clusters of Orthologous Groups of proteins (COGs) from the sequenced genomes of prokaryotes and unicellular eukaryotes and the construction of clusters of predicted orthologs for 7 eukaryotic genomes, which we named KOGs after eukaryotic orthologous groups. The COG collection currently consists of 138,458 proteins, which form 4873 COGs and comprise 75% of the 185,505 (predicted) proteins encoded in 66 genomes of unicellular organisms. The eukaryotic orthologous groups (KOGs) include proteins from 7 eukaryotic genomes: three animals (the nematode Caenorhabditis elegans, the fruit fly Drosophila melanogaster and Homo sapiens), one plant, Arabidopsis thaliana, two fungi (Saccharomyces cerevisiae and Schizosaccharomyces pombe), and the intracellular microsporidian parasite Encephalitozoon cuniculi. The current KOG set consists of 4852 clusters of orthologs, which include 59,838 proteins, or approximately 54% of the analyzed eukaryotic 110,655 gene products. Compared to the coverage of the prokaryotic genomes with COGs, a considerably smaller fraction of eukaryotic genes could be included into the KOGs; addition of new eukaryotic genomes is expected to result in substantial increase in the coverage of eukaryotic genomes with KOGs. Examination of the phyletic patterns of KOGs reveals a conserved core represented in all analyzed species and consisting of approximately 20% of the KOG set. This conserved portion of the KOG set is much greater than the ubiquitous portion of the COG set (approximately 1% of the COGs). In part, this difference is probably due to the small number of included eukaryotic genomes, but it could also reflect the relative compactness of eukaryotes as a clade and the greater evolutionary stability of eukaryotic genomes.

CONCLUSION

The updated collection of orthologous protein sets for prokaryotes and eukaryotes is expected to be a useful platform for functional annotation of newly sequenced genomes, including those of complex eukaryotes, and genome-wide evolutionary studies.

Chromosomal mapping of ANTP class homeobox genes in amphioxus: piecing together ancestral genomes

Homeobox genes encode DNA-binding proteins, many of which are implicated in the control of embryonic development. Evolutionarily, most homeobox genes fall into two related clades: the ANTP and the PRD classes. Some genes in ANTP class, notably Hox, ParaHox, and NK genes, have an intriguing arrangement into physical clusters. To investigate the evolutionary history of these gene clusters, we examined homeobox gene chromosomal locations in the cephalochordate amphioxus, Branchiostoma floridae. We deduce that 22 amphioxus ANTP class homeobox genes localize in just three chromosomes. One contains the Hox cluster plus AmphiEn, AmphiMnx, and AmphiDll. The ParaHox cluster resides in another chromosome, whereas a third chromosome contains the NK type homeobox genes, including AmphiMsx and AmphiTlx. By comparative analysis we infer that clustering of ANTP class homeobox genes evolved just once, during a series of extensive cis-duplication events of genes early in animal evolution. A trans-duplication event occurred later to yield the Hox and ParaHox gene clusters on different chromosomes. The results obtained have implications for understanding the origin of homeobox gene clustering, the diversification of the ANTP class of homeobox genes, and the evolution of animal genomes.

A Twist-like bHLH gene is a downstream factor of an endogenous FGF and determines mesenchymal fate in the ascidian embryos

Ascidian larvae develop mesenchyme cells in their trunk. A fibroblast growth factor (FGF9/16/20) is essential and sufficient for induction of the mesenchyme in Ciona savignyi. We have identified two basic helix-loop-helix (bHLH) genes named Twist-like1 and Twist-like2 as downstream factors of this FGF. These two genes are phylogenetically closely related to each other, and were expressed specifically in the mesenchymal cells after the 110-cell stage. Gene-knockdown experiments using a specific morpholino oligonucleotide demonstrated that Twist-like1 plays an essential role in determination of the mesenchyme and that Twist-like2 is a downstream factor of Twist-like1. In addition, both overexpression and misexpression of Twist-like1 converts non-mesenchymal cells to mesenchymal cells. We also demonstrate that the upstream regulatory mechanisms of Twist-like1 are different between B-line mesenchymal cells and the A-line mesenchymal cells called 'trunk lateral cells'. FGF9/16/20 is required for the expression of Twist-like1 in B-line mesenchymal precursor cells, whereas FGF, FoxD and another novel bHLH factor called NoTrlc are required for Twist-like1 to be expressed in the A-line mesenchymal precursor cells. Therefore, two different but partially overlapping mechanisms are required for the expression of Twist-like1 in the mesenchymal precursors, which triggers the differentiation of the mesenchyme in Ciona embryos.

Targeted knockdown of an opsin gene inhibits the swimming behaviour photoresponse of ascidian larvae

The sequencing of the ascidian Ciona intestinalis genome was completed at the end of 2002. Effective targeted gene knockdown in this model chordate would greatly enhance understanding of how the genes affect the function of the neurons that underlie behaviour. We show here that antisense morpholinos (MOs) are effective and specific translational inhibitors in C. intestinalis larvae. The larvae developed from eggs injected with a Ci-opsin1 MO lost their photoresponse, which affected their swimming behaviour. Immunohistochemical analysis showed that the expression of Ci-opsin1, the visual pigments in the photoreceptor of the ocellus, had been lost in the larvae with the Ci-opsin1 MO. The inhibition of the photic behaviour of the knockdown larvae is solely due to the suppression of the synthesis of the visual pigments in the photoreceptors. These are the first results showing the knockdown of an ascidian gene and the elucidation of that gene's role in ascidian larval behaviour.

Evidence for the presence of a cellulase gene in the last common ancestor of bilaterian animals

Until recently, the textbook view of cellulose hydrolysis in animals was that gut-resident symbiotic organisms such as bacteria or unicellular eukaryotes are responsible for the cellulases produced. This view has been challenged by the characterization and sequencing of endogenous cellulase genes from some invertebrate animals, including plant-parasitic nematodes, arthropods and a mollusc. Most of these genes are completely unrelated in terms of sequence, and their evolutionary origins remain unclear. In the case of plant-parasitic nematodes, it has been suggested that their ancestor obtained a cellulase gene via horizontal gene transfer from a prokaryote, and similar suggestions have been made about a cellulase gene recently discovered in a sea squirt. To improve understanding about the evolution of animal cellulases, we searched for all known types of these enzymes in GenBank, and performed phylogenetic comparisons. Low phylogenetic resolution was found among most of the sequences examined, however, positional identity in the introns of cellulase genes from a termite, a sea squirt and an abalone provided compelling evidence that a similar gene was present in the last common ancestor of protostomes and deuterostomes. In a different enzyme family, cellulases from beetles and plant-parasitic nematodes were found to cluster together. This result questions the idea of lateral gene transfer into the ancestors of the latter, although statistical tests did not allow this possibility to be ruled out. Overall, our results suggest that at least one family of endogenous cellulases may be more widespread in animals than previously thought.

Membrane proteins with immunoglobulin-like domains—a master superfamily of interaction molecules

Recent genomic analysis has shown that proteins with immunoglobulin superfamily (IgSF) domains are extremely abundant and their number has increased markedly in evolution correlating with the development of the adaptive immune system. The IgSF domain is particularly good at being recognised and is involved in many different kinds of interactions. Thus, analysis of the properties of these domains can act as a paradigm for thinking about the roles of newly identified gene products. This review summarises the identification, function and properties of IgSF domains including, their size, variety of interactions, their strength of binding, role of glycosylation and organisation with other proteins.

Presence of the delta-MSH sequence in a proopiomelanocortin cDNA cloned from the pituitary of the galeoid shark, Heterodontus portusjacksoni

Since a fourth MSH sequence, delta-MSH, has been detected in the proopiomelanocortin (POMC) gene of a dogfish and a stingray, members of superorder Squalea (class Chondrichthyes), it is possible that this novel MSH sequence might be a feature common to the POMC genes of all modern sharks and rays. As an initial step towards addressing this question, a full-length POMC cDNA was cloned and sequenced from the pituitary of the Port Jackson shark, Heterodontus portusjacksoni. The Port Jackson shark represents one of the oldest lineages in superorder Galea, and this superorder together with superorder Squalea form infraclass Neoselachii (the extant sharks and rays). The Port Jackson shark POMC cDNA has an open reading frame that is 1032 nucleotides in length and encodes the deduced amino acids sequences for beta-endorphin, ACTH/alpha-MSH, beta-MSH, gamma-MSH, and delta-MSH. Port Jackson shark delta-MSH has 83% primary sequence identity with dogfish and stingray delta-MSH, and it appears that the delta-MSH sequence may have been the result of an internal domain duplication and reinsertion of the beta-MSH sequence. The presence of the delta-MSH sequence in the POMC genes of representatives of both superorders of infraclass Neoselachii would indicate that the delta-MSH sequence must have been present in the ancestral euselachian shark that gave rise to the neoselachian radiation.

Schizosaccharomyces pombe as a model for metal homeostasis in plant cells: the phytochelatin-dependent pathway is the main cadmium detoxification mechanism

Sequestration of metal ions by phytochelatins is an important metal tolerance mechanism in a wide range of organisms including plants and certain fungi. Substantial progress in understanding phytochelatin formation at the molecular level has been made in Schizosaccharomyces pombe. The genome of S. pombe has been completely sequenced and all the necessary tools of functional genomics are available. Since most other proteins implicated in plant metal tolerance and homeostasis are also present in this yeast, it represents a very powerful system to elucidate basic mechanisms of metal buffering, sequestration, and toxicity in cells that form phytochelatins. Here, we summarize the work on phytochelatin formation and metal homeostasis in S. pombe. We describe examples of molecular insights obtained from experiments with S. pombe that will be useful in guiding studies with plants. We also provide evidence for the dominance of the phytochelatin pathway in Cd detoxification in S. pombe.

Isolation and characterization of beta-catenin downstream genes in early embryos of the ascidian Ciona savignyi

Nuclear localization of beta-catenin is most likely the first step of embryonic axis formation or embryonic cell specification in a wide variety of animal groups. Therefore, the elucidation of beta-catenin target genes is a key research subject in understanding the molecular mechanisms of the early embryogenesis of animals. In Ciona savignyi embryos, nuclear accumulation of beta-catenin is the first step of endodermal cell specification. Previous subtractive hybridization screens of mRNAs between beta-catenin-overexpressed embryos and nuclear beta-catenin-depleted embryos have resulted in the identification of beta-catenin downstream genes in Ciona embryos. In the present study, I characterize seven additional beta-catenin downstream genes, Cs-cadherinII, Cs-protocadherin, Cs-Eph, Cs-betaCD1, Cs-netrin, Cs-frizzled3/6, and Cs-lefty/antivin. All of these genes were expressed in vegetal blastomeres between the 16-cell and 110-cell stages, although their spatial and temporal expression patterns were different from one another. In situ hybridizations and real-time PCR revealed that the expression of all of these genes was up-regulated in beta-catenin-overexpressed embryos, and down-regulated in beta-catenin-suppressed embryos. Therefore, the accumulation of beta-catenin in the nuclei of vegetal blastomeres activates various vegetally expressed genes with potentially important functions in the specification of these cells.

Comparative genomic analysis reveals independent expansion of a lineage-specific gene family in vertebrates: the class II cytokine receptors and their ligands in mammals and fish

BACKGROUND

The high degree of sequence conservation between coding regions in fish and mammals can be exploited to identify genes in mammalian genomes by comparison with the sequence of similar genes in fish. Conversely, experimentally characterized mammalian genes may be used to annotate fish genomes. However, gene families that escape this principle include the rapidly diverging cytokines that regulate the immune system, and their receptors. A classic example is the class II helical cytokines (HCII) including type I, type II and lambda interferons, IL10 related cytokines (IL10, IL19, IL20, IL22, IL24 and IL26) and their receptors (HCRII). Despite the report of a near complete pufferfish (Takifugu rubripes) genome sequence, these genes remain undescribed in fish.

RESULTS

We have used an original strategy based both on conserved amino acid sequence and gene structure to identify HCII and HCRII in the genome of another pufferfish, Tetraodon nigroviridis that is amenable to laboratory experiments. The 15 genes that were identified are highly divergent and include a single interferon molecule, three IL10 related cytokines and their potential receptors together with two Tissue Factor (TF). Some of these genes form tandem clusters on the Tetraodon genome. Their expression pattern was determined in different tissues. Most importantly, Tetraodon interferon was identified and we show that the recombinant protein can induce antiviral MX gene expression in Tetraodon primary kidney cells. Similar results were obtained in Zebrafish which has 7 MX genes.

CONCLUSION

We propose a scheme for the evolution of HCII and their receptors during the radiation of bony vertebrates and suggest that the diversification that played an important role in the fine-tuning of the ancestral mechanism for host defense against infections probably followed different pathways in amniotes and fish.

The C. elegans Hand gene controls embryogenesis and early gonadogenesis

The C. elegans genome encodes a single Hand bHLH transcription factor. Either hnd-1(RNAi) or a hnd-1 deletion causes partially penetrant defects in viability and gonadogenesis. Dead embryos and young larvae are often misshapen at the posterior end. Our primary focus has been the role of hnd-1 in gonadogenesis. Wild-type C. elegans has two somatic gonadal precursors and two primordial germ cells in stereotyped positions within its four-celled gonadal primordium. The hnd-1 gene affects the presence and position of both the somatic gonadal precursors and primordial germ cells within the primordium, but does not appear to have any role in later gonadogenesis. hnd-1 probably acts within the somatic gonadal precursors or their mesodermal predecessors; defects in primordial germ cells and germ line appear to be secondary. In hnd-1 mutants, somatic gonadal precursors are generated normally, but are not maintained properly and sometimes die. A similar role in controlling the maintenance of precursor fates has been described for other genes governing early organogenesis, including the zebrafish Hand gene hands off. We also report the discovery of two genes, ehn-1 and ehn-3, that have overlapping functions with hnd-1 in embryogenesis and gonadogenesis.