EST analysis of gene expression in the tentacle of Cyanea capillata

Maintenance of a Protein Structure in the Dynamic Evolution of TIMPs over 600 Million Years

Deciphering the events leading to protein evolution represents a challenge, especially for protein families showing complex evolutionary history. Among them, TIMPs represent an ancient eukaryotic protein family widely distributed in the animal kingdom. They are known to control the turnover of the extracellular matrix and are considered to arise early during metazoan evolution, arguably tuning essential features of tissue and epithelial organization. To probe the structure and molecular evolution of TIMPs within metazoans, we report the mining and structural characterization of a large data set of TIMPs over approximately 600 Myr. The TIMPs repertoire was explored starting from the Cnidaria phylum, coeval with the origins of connective tissue, to great apes and humans. Despite dramatic sequence differences compared with highest metazoans, the ancestral proteins displayed the canonical TIMP fold. Only small structural changes, represented by an α-helix located in the N-domain, have occurred over the evolution. Both the occurrence of such secondary structure elements and the relative solvent accessibility of the corresponding residues in the three-dimensional structures raises the possibility that these sites represent unconserved element prone to accept variations.

Tentacle Transcriptome and Venom Proteome of the Pacific Sea Nettle, Chrysaora fuscescens (Cnidaria: Scyphozoa)

Jellyfish venoms are rich sources of toxins designed to capture prey or deter predators, but they can also elicit harmful effects in humans. In this study, an integrated transcriptomic and proteomic approach was used to identify putative toxins and their potential role in the venom of the scyphozoan jellyfish Chrysaora fuscescens. A de novo tentacle transcriptome, containing more than 23,000 contigs, was constructed and used in proteomic analysis of C. fuscescens venom to identify potential toxins. From a total of 163 proteins identified in the venom proteome, 27 were classified as putative toxins and grouped into six protein families: proteinases, venom allergens, C-type lectins, pore-forming toxins, glycoside hydrolases and enzyme inhibitors. Other putative toxins identified in the transcriptome, but not the proteome, included additional proteinases as well as lipases and deoxyribonucleases. Sequence analysis also revealed the presence of ShKT domains in two putative venom proteins from the proteome and an additional 15 from the transcriptome, suggesting potential ion channel blockade or modulatory activities. Comparison of these potential toxins to those from other cnidarians provided insight into their possible roles in C. fuscescens venom and an overview of the diversity of potential toxin families in cnidarian venoms.

Transcriptomics and venomics: implications for medicinal chemistry

Over the last three decades, transcriptomic studies of venom gland cells have continuously evolved, opening up new possibilities for exploring the molecular diversity of animal venoms, a prerequisite for the discovery of new drug candidates and molecular phylogenetics. The molecular complexity of animal venoms is much greater than initially thought. In this review, we describe the different technologies available for transcriptomic studies of venom, from the original individual cloning approaches to the more recent global Next Generation Sequencing strategies. Our understanding of animal venoms is evolving, with the discovery of complex and diverse bio-optimized cocktails of compounds, including mostly peptides and proteins, which are now beginning to be studied by academic and industrial researchers.

Characterization of small HSPs from Anemonia viridis reveals insights into molecular evolution of alpha crystallin genes among cnidarians

Gene family encoding small Heat-Shock Proteins (sHSPs containing α-crystallin domain) are found both in prokaryotic and eukaryotic organisms; however, there is limited knowledge of their evolution. In this study, two small HSP genes termed AvHSP28.6 and AvHSP27, both organized in one intron and two exons, were characterised in the Mediterranean snakelocks anemone Anemonia viridis. The release of the genome sequence of Hydra magnipapillata and Nematostella vectensis enabled a comprehensive study of the molecular evolution of α-crystallin gene family among cnidarians. Most of the H. magnipapillata sHSP genes share the same gene organization described for AvHSP28.6 and AvHSP27, differing from the sHSP genes of N. vectensis which mainly show an intronless architecture. The different genomic organization of sHSPs, the phylogenetic analyses based on protein sequences, and the relationships among Cnidarians, suggest that the A.viridis sHSPs represent the common ancestor from which H. magnipapillata genes directly evolved through segmental genome duplication. Additionally retroposition events may be considered responsible for the divergence of sHSP genes of N. vectensis from A. viridis. Analyses of transcriptional expression profile showed that AvHSP28.6 was constitutively expressed among different tissues from both ectodermal and endodermal layers of the adult sea anemones, under normal physiological conditions and also under different stress condition. Specifically, we profiled the transcriptional activation of AvHSP28.6 after challenges with different abiotic/biotic stresses showing induction by extreme temperatures, heavy metals exposure and immune stimulation. Conversely, no AvHSP27 transcript was detected in such dissected tissues, in adult whole body cDNA library or under stress conditions. Hence, the involvement of AvHSP28.6 gene in the sea anemone defensome is strongly suggested.

Comparative venom gland transcriptome analysis of the scorpion Lychas mucronatus reveals intraspecific toxic gene diversity and new venomous components

Background

Lychas mucronatus is one scorpion species widely distributed in Southeast Asia and southern China. Anything is hardly known about its venom components, despite the fact that it can often cause human accidents. In this work, we performed a venomous gland transcriptome analysis by constructing and screening the venom gland cDNA library of the scorpion Lychas mucronatus from Yunnan province and compared it with the previous results of Hainan-sourced Lychas mucronatus.

Results

A total of sixteen known types of venom peptides and proteins are obtained from the venom gland cDNA library of Yunnan-sourced Lychas mucronatus, which greatly increase the number of currently reported scorpion venom peptides. Interestingly, we also identified nineteen atypical types of venom molecules seldom reported in scorpion species. Surprisingly, the comparative transcriptome analysis of Yunnan-sourced Lychas mucronatus and Hainan-sourced Lychas mucronatus indicated that enormous diversity and vastly abundant difference could be found in venom peptides and proteins between populations of the scorpion Lychas mucronatus from different geographical regions.

Conclusions

This work characterizes a large number of venom molecules never identified in scorpion species. This result provides a comparative analysis of venom transcriptomes of the scorpion Lychas mucronatus from different geographical regions, which thoroughly reveals the fact that the venom peptides and proteins of the same scorpion species from different geographical regions are highly diversified and scorpion evolves to adapt a new environment by altering the primary structure and abundance of venom peptides and proteins.

Phospholipases A(2) in the genome of the sea anemone Nematostella vectensis

The genome of the sea anemone Nematostella vectensis (Nv) (Cnidaria, Anthozoa) was sequenced recently (Putnam et al., Science 317: 86, 2007). In the current study, 22 proteins of Nv were identified as putative phospholipases A(2) (PLA(2)) that showed up to 40-50% sequence identity with secreted or intracellular PLA(2)s including those of humans. Nv1-Nv6 PLA(2)s have identity with secreted human group (G)IB and GIIA PLA(2)s and PLA(2)s of the sea anemones Adamsia carciniopados and Urticina crassicornis. Nv7 and Nv8 PLA(2)s have identity with human and bee venom GIII PLA(2)s and Nv9 PLA(2) with GXIIA PLA(2). Nv10-Nv13 PLA(2)s show identity with GIX PLA(2) of Conus magus and bacterial PLA(2)s but no significant identity with any human PLA(2). Nv14 has identity with intracellular GIV PLA(2), Nv15 with GVII PLA(2), Nv16 and Nv17 with GVIII PLA(2), Nv18-Nv20 with GVI PLA(2), and Nv21 and Nv22 with patatin, respectively. The observations indicate that the cnidarian phospholipasome contains a rich array of orthologs of most types of animal PLA(2)s, and that many of the PLA(2)-driven vital functions prevail in these ancient metazoans. Cnidarian PLA(2)s may be considered as evolutionary precursors of PLA(2)s of higher animals.

Analysis of expressed sequence tags of the cyclically parthenogenetic rotifer Brachionus plicatilis

BACKGROUND

Rotifers are among the most common non-arthropod animals and are the most experimentally tractable members of the basal assemblage of metazoan phyla known as Gnathifera. The monogonont rotifer Brachionus plicatilis is a developing model system for ecotoxicology, aquatic ecology, cryptic speciation, and the evolution of sex, and is an important food source for finfish aquaculture. However, basic knowledge of the genome and transcriptome of any rotifer species has been lacking.

METHODOLOGY/PRINCIPAL FINDINGS

We generated and partially sequenced a cDNA library from B. plicatilis and constructed a database of over 2300 expressed sequence tags corresponding to more than 450 transcripts. About 20% of the transcripts had no significant similarity to database sequences by BLAST; most of these contained open reading frames of significant length but few had recognized Pfam motifs. Sixteen transcripts accounted for 25% of the ESTs; four of these had no significant similarity to BLAST or Pfam databases. Putative up- and downstream untranslated regions are relatively short and AT rich. In contrast to bdelloid rotifers, there was no evidence of a conserved trans-spliced leader sequence among the transcripts and most genes were single-copy.

CONCLUSIONS/SIGNIFICANCE

Despite the small size of this EST project it revealed several important features of the rotifer transcriptome and of individual monogonont genes. Because there is little genomic data for Gnathifera, the transcripts we found with no known function may represent genes that are species-, class-, phylum- or even superphylum-specific; the fact that some are among the most highly expressed indicates their importance. The absence of trans-spliced leader exons in this monogonont species contrasts with their abundance in bdelloid rotifers and indicates that the presence of this phenomenon can vary at the subphylum level. Our EST database provides a relatively large quantity of transcript-level data for B. plicatilis, and more generally of rotifers and other gnathiferan phyla, and can be browsed and searched at gmod.mbl.edu.

The cnidarian-bilaterian ancestor possessed at least 56 homeoboxes: evidence from the starlet sea anemone, Nematostella vectensis

BACKGROUND

Homeodomain transcription factors are key components in the developmental toolkits of animals. While this gene superclass predates the evolutionary split between animals, plants, and fungi, many homeobox genes appear unique to animals. The origin of particular homeobox genes may, therefore, be associated with the evolution of particular animal traits. Here we report the first near-complete set of homeodomains from a basal (diploblastic) animal.

RESULTS

Phylogenetic analyses were performed on 130 homeodomains from the sequenced genome of the sea anemone Nematostella vectensis along with 228 homeodomains from human and 97 homeodomains from Drosophila. The Nematostella homeodomains appear to be distributed among established homeodomain classes in the following fashion: 72 ANTP class; one HNF class; four LIM class; five POU class; 33 PRD class; five SINE class; and six TALE class. For four of the Nematostella homeodomains, there is disagreement between neighbor-joining and Bayesian trees regarding their class membership. A putative Nematostella CUT class gene is also identified.

CONCLUSION

The homeodomain superclass underwent extensive radiations prior to the evolutionary split between Cnidaria and Bilateria. Fifty-six homeodomain families found in human and/or fruit fly are also found in Nematostella, though seventeen families shared by human and fly appear absent in Nematostella. Homeodomain loss is also apparent in the bilaterian taxa: eight homeodomain families shared by Drosophila and Nematostella appear absent from human (CG13424, EMXLX, HOMEOBRAIN, MSXLX, NK7, REPO, ROUGH, and UNC4), and six homeodomain families shared by human and Nematostella appear absent from fruit fly (ALX, DMBX, DUX, HNF, POU1, and VAX).

Transcriptome analysis of the venom gland of the Mexican scorpion Hadrurus gertschi (Arachnida: Scorpiones)

Background

Scorpions like other venomous animals posses a highly specialized organ that produces, secretes and disposes the venom components. In these animals, the last postabdominal segment, named telson, contains a pair of venomous glands connected to the stinger. The isolation of numerous scorpion toxins, along with cDNA-based gene cloning and, more recently, proteomic analyses have provided us with a large collection of venom components sequences. However, all of them are secreted, or at least are predicted to be secretable gene products. Therefore very little is known about the cellular processes that normally take place inside the glands for production of the venom mixture. To gain insights into the scorpion venom gland biology, we have decided to perform a transcriptomic analysis by constructing a cDNA library and conducting a random sequencing screening of the transcripts.

Results

From the cDNA library prepared from a single venom gland of the scorpion Hadrurus gertschi, 160 expressed sequence tags (ESTs) were analyzed. These transcripts were further clustered into 68 unique sequences (20 contigs and 48 singlets), with an average length of 919 bp. Half of the ESTs can be confidentially assigned as homologues of annotated gene products. Annotation of these ESTs, with the aid of Gene Ontology terms and homology to eukaryotic orthologous groups, reveals some cellular processes important for venom gland function; including high protein synthesis, tuned posttranslational processing and trafficking. Nonetheless, the main group of the identified gene products includes ESTs similar to known scorpion toxins or other previously characterized scorpion venom components, which account for nearly 60% of the identified proteins.

Conclusion

To the best of our knowledge this report contains the first transcriptome analysis of genes transcribed by the venomous gland of a scorpion. The data were obtained for the species Hadrurus gertschi, belonging to the family Caraboctonidae. One hundred and sixty ESTs were analyzed, showing enrichment in genes that encode for products similar to known venom components, but also provides the first sketch of cellular components, molecular functions, biological processes and some unique sequences of the scorpion venom gland.

Diversity and evolution of conotoxins based on gene expression profiling of Conus litteratus

Cone snails are attracting increasing scientific attention due to their unprecedented diversity of invaluable channel-targeted peptides. As arguably the largest and most successful evolutionary genus of invertebrates, Conus also may become the model system to study the evolution of multigene families and biodiversity. Here, a set of 897 expressed sequence tags (ESTs) derived from a Conus litteratus venom duct was analyzed to illuminate the diversity and evolution mechanism of conotoxins. Nearly half of these ESTs represent the coding sequences of conotoxins, which were grouped into 42 novel conotoxin cDNA sequences (seven superfamilies), with T-superfamily conotoxins being the dominant component. The gene expression profile of conotoxin revealed that transcripts are expressed with order-of-magnitude differences, sequence divergence within a superfamily increases from the N to the C terminus of the open reading frame, and even multiple scaffold-different mature peptides exist in a conotoxin gene superfamily. Most excitingly, we identified a novel conotoxin superfamily and three novel cysteine scaffolds. These results give an initial insight into the C. litteratus transcriptome that will contribute to a better understanding of conotoxin evolution and the study of the cone snail genome in the near future.

Expansion of the SOX gene family predated the emergence of the Bilateria

Members of the SOX gene family are involved in regulating many developmental processes including neuronal determination and differentiation, and in carcinogenesis. So far they have only been identified in species from the Bilateria (deuterostomes and protostomes). To understand the origins of the SOX family, we used a PCR-based strategy to obtain 28 new sequences of SOX gene HMG domains from four non-bilaterian Metazoa: two sponge species, one ctenophore and one cnidarian. One additional SOX sequence was retrieved from EST sequences of the cnidarian species Clytia hemisphaerica. Unexpected SOX gene diversity was found in these species, especially in the cnidarian and the ctenophore. The topology of gene relationships deduced by Maximum Likelihood analysis, although not supported by bootstrap values, suggested that the SOX family started to diversify in the metazoan stem branch prior to the divergence of demosponges, and that further diversification occurred in the eumetazoan branch, as well as later in calcisponges, ctenophores, cnidarians and vertebrates. In contrast, gene loss appears to have occurred in the nematode and probably in other protostome lineages, explaining their lower number of SOX genes.

Silencing of the hydra serine protease inhibitorKazal1gene mimics the humanSPINK1pancreatic phenotype

In hydra, the endodermal epithelial cells carry out the digestive function together with the gland cells that produce zymogens and express the evolutionarily conserved gene Kazal1. To assess the hydra Kazal1 function, we silenced gene expression through double-stranded RNA feeding. A progressive Kazal1 silencing affected homeostatic conditions as evidenced by the low budding rate and the induced animal death. Concomitantly, a dramatic disorganization followed by a massive death of gland cells was observed, whereas the cytoplasm of digestive cells became highly vacuolated. The presence of mitochondria and late endosomes within those vacuoles assigned them as autophagosomes. The enhanced Kazal1 expression in regenerating tips was strongly diminished in Kazal1(–) hydra, and the amputation stress led to an immediate disorganization of the gland cells, vacuolization of the digestive cells and death after prolonged silencing. This first cellular phenotype resulting from a gene knock-down in cnidarians suggests that the Kazal1 serine-protease-inhibitor activity is required to prevent excessive autophagy in intact hydra and to exert a cytoprotective function to survive the amputation stress. Interestingly, these functions parallel the pancreatic autophagy phenotype observed upon mutation within the Kazal domain of the SPINK1 and SPINK3 genes in human and mice, respectively.

StellaBase: the Nematostella vectensis Genomics Database

StellaBase, the Nematostella vectensis Genomics Database, is a web-based resource that will facilitate desktop and bench-top studies of the starlet sea anemone. Nematostella is an emerging model organism that has already proven useful for addressing fundamental questions in developmental evolution and evolutionary genomics. StellaBase allows users to query the assembled Nematostella genome, a confirmed gene library, and a predicted genome using both keyword and homology based search functions. Data provided by these searches will elucidate gene family evolution in early animals. Unique research tools, including a Nematostella genetic stock library, a primer library, a literature repository and a gene expression library will provide support to the burgeoning Nematostella research community. The development of StellaBase accompanies significant upgrades to CnidBase, the Cnidarian Evolutionary Genomics Database. With the completion of the first sequenced cnidarian genome, genome comparison tools have been added to CnidBase. In addition, StellaBase provides a framework for the integration of additional species-specific databases into CnidBase. StellaBase is available at .

Analysis of expressed sequence tags from the venom ducts of Conus striatus: focusing on the expression profile of conotoxins

Cone snails (genus Conus) are predatory marine gastropods that use venom peptides for interacting with prey, predators and competitors. A majority of these peptides, generally known as conotoxins demonstrate striking selectivity in targeting specific subtypes of ion channels and neurotransmitter receptors. So they are not only useful tools in neuroscience to characterize receptors and receptor subtypes, but offer great potential in new drug research and development as well. Here, a cDNA library from the venom ducts of a fish-hunting cone snail species, Conus striatus is described for the generation of expressed sequence tags (ESTs). A total of 429 ESTs were grouped into 137 clusters or singletons. Among these sequences, 221 were toxin sequences, accounting for 52.1% (corresponding to 19 clusters) of all transcripts. A-superfamily (132 ESTs) and O-superfamily conotoxins (80 ESTs) constitute the predominant toxin components. Some non-disulfide-rich Conus peptides were also found. The expression profile of conotoxins also explained to some extent the pharmacological and physiological reactions elicited by this typical piscivorous species. For the first time, a nonstop transcript of conotoxin was identified, which is suggestive that alternative polyadenylation may be a means of post-transcriptional regulation of conotoxin production. A comparison analysis of these conotoxins reveals the different variation and divergence patterns in these two superfamilies. Our investigations indicate that focal hyper-mutation, block substitution and exon shuffling are three main mechanisms leading to the conotoxin diversity in a species. The comprehensive set of Conus gene sequences allowed the identification of the representative classes of conotoxins and related components, which may lay the foundation for further research and development of conotoxins.

T-box and homeobox genes from the ctenophore Pleurobrachia pileus: comparison of Brachyury, Tbx2/3 and Tlx in basal metazoans and bilaterians

Most animals are classified as Bilateria and only four phyla are still extant as outgroups, namely Porifera, Placozoa, Cnidaria and Ctenophora. These non-bilaterians were not considered to have a mesoderm and hence mesoderm-specific genes. However, the T-box gene Brachyury could be isolated from sponges, placozoans and cnidarians. Here, we describe the first Brachyury and a Tbx2/3 homologue from a ctenophore. In addition, analysing T-box and homeobox genes under comparable conditions in all four basal phyla lead to the discovery of novel T-box genes in sponges and cnidarians and a Tlx homeobox gene in the ctenophore Pleurobrachia pileus. The conservation of the T-box and the homeobox genes suggest that distinct subfamilies with different roles in bilaterians were already split in non-bilaterians.

Invertebrate Muscles: Muscle Specific Genes and Proteins

This is the first of a projected series of canonic reviews covering all invertebrate muscle literature prior to 2005 and covers muscle genes and proteins except those involved in excitation-contraction coupling (e.g., the ryanodine receptor) and those forming ligand- and voltage-dependent channels. Two themes are of primary importance. The first is the evolutionary antiquity of muscle proteins. Actin, myosin, and tropomyosin (at least, the presence of other muscle proteins in these organisms has not been examined) exist in muscle-like cells in Radiata, and almost all muscle proteins are present across Bilateria, implying that the first Bilaterian had a complete, or near-complete, complement of present-day muscle proteins. The second is the extraordinary diversity of protein isoforms and genetic mechanisms for producing them. This rich diversity suggests that studying invertebrate muscle proteins and genes can be usefully applied to resolve phylogenetic relationships and to understand protein assembly coevolution. Fully achieving these goals, however, will require examination of a much broader range of species than has been heretofore performed.

Toxic polypeptides of the hydra—a bioinformatic approach to cnidarian allomones

Cnidarians such as hydrae and sea anemones are sessile, predatory, soft bodied animals which depend on offensive and defensive allomones for prey capture and survival. These allomones are distributed throughout the entire organism both in specialized stinging cells (nematocytes) and in the body tissues. The cnidarian allomonal system is composed of neurotoxins, cytolysins and toxic phospholipapses. The present bioinformatic survey was motivated by the fact that while hydrae are the most studied model cnidarian, little is known about their allomones. A large-scale EST database from Hydra magnipapillata was searched for orthologs of known cnidarian allomones, as well as for allomones found in other venomous organisms. We show that the hydrae express orthologs of cnidarian phospholipase A2 toxins and cytolysins belonging to the actinoporin family, but could not find orthologs of the 'classic' short chain neurotoxins affecting sodium and potassium conductance. Hydrae also express proteins similar to elapid-like phospholipases, CRISP proteins, Prokineticin-like polypeptides and toxic deoxyribonucleases. Our results illustrate a high level of complexity in the hydra allomonal system, suggest that several toxins represent a basal component of all cnidarian allomones, and raise the intriguing possibility that similar proteins may fulfill both endogenous and allomonal roles in cnidaria.

Rising starlet: the starlet sea anemone, Nematostella vectensis

In recent years, a handful of model systems from the basal metazoan phylum Cnidaria have emerged to challenge long-held views on the evolution of animal complexity. The most-recent, and in many ways most-promising addition to this group is the starlet sea anemone, Nematostella vectensis. The remarkable amenability of this species to laboratory manipulation has already made it a productive system for exploring cnidarian development, and a proliferation of molecular and genomic tools, including the currently ongoing Nematostella genome project, further enhances the promise of this species. In addition, the facility with which Nematostella populations can be investigated within their natural ecological context suggests that this model may be profitably expanded to address important questions in molecular and evolutionary ecology. In this review, we explore the traits that make Nematostella exceptionally attractive as a model organism, summarize recent research demonstrating the utility of Nematostella in several different contexts, and highlight a number of developments likely to further increase that utility in the near future.

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.

The evolution of multicellularity and early animal genomes

Several independent molecular datasets, including complete mtDNA sequence, indicate that Choanozoa are most closely related to multicellular animals. There is still confusion concerning basal animal phylogeny, although recent data indicate that Placozoa are not degenerate cnidarians and hence (along with sponges) occupy a pivotal position. The transition in evolution from diploblast to bilaterian animals is becoming better understood, with gene expression data arguing that cnidarians have forerunners of the anteroposterior and dorsoventral body axes, and even a putative homologue of mesoderm. The homeobox and kinase gene families have been further analysed in basal animals, although more data are required to enable detailed comparison with Bilateria.