The complete mitochondrial genome of the Hwanghae Rockfish Sebastes koreanus (Scorpaenidae, Scorpaeniformes).

The complete mitochondrial genome sequence of Sebastes koreanus was determined by next-generation sequencing. The total genome was 16,499 bp in length, comprising two rRNA, 22 tRNA and 13 protein-coding genes (PCGs), as well as the mtDNA control region. The genome constitutes 27.9% A, 26.6% T, 17.1% G and 28.4% C, showing a slight AT bias (54.5%). Of the 13 mitochondrial PCGs, 12 have an ATG start codon while COX1 had a GTG start codon; the incomplete stop codon (T) was present in the ND3, ND4 and Cytb genes. The anti-G bias was considerably lower at the third codon position of PCGs (8.8%). The complete mitogenome sequence of S. koreanus will contribute to molecular phylogenetic and evolutionary analyses of Sebastinae and to the identification of its genetic populations in the Yellow Sea.

Expressed sequence tags in venomous tissue of Scorpaena plumieri (Scorpaeniformes: Scorpaenidae)

Species of the family Scorpaenidae are responsible for accidents and sporadic casualties by the shore they inhabit. The species Scorpaena plumieri from this family populate the Northeastern and Eastern coast of Brazil causing human envenomation characterized by local and systemic symptoms. In experimental animals the venom induces cardiotoxic, hypotensive, and airway respiratory effects. As first step to identify the venom components we isolated gland mRNA to produce a cDNA library from the fish gland. This report describes the partial sequencing of 356 gland transcripts from S. plumieri. BLAST analysis of transcripts showed that 30% were unknown sequences, 17% hypothetical proteins, 17% related to metabolic enzymes, 14% belonged to signal transducing functions and the remaining groups (7-8%) composed by gene related with expressing proteins, regulatory proteins and structural proteins. A considerable number of these EST were not found in available databases suggesting the existence of new proteins and/or functions yet to be discovered. By screening the library with antibodies against a lectin fraction from S. plumieri venom we identified several clones whose DNA sequence showed similarities with lectins found in fish. In silico analysis of these clones confirm the identity of these molecules in the venom gland of S. plumieri. .

Espécies da família Scorpaenidae são responsáveis por acidentes e mortes esporádicas ao longo da costa que habitam. A espécie Scorpaena plumieri desta família povoam a costa Leste e Nordeste do Brasil, causando envenenamento humano caracterizado por sintomas locais e sistêmicos. Em modelos experimentais animais a peçonha induz cardiotoxicidade, efeitos hipotensivos e alterações nas vias aéreas respiratórias. Como primeiro passo para identificar os componentes da peçonha foram isolados os mRNA das glândulas do peixe para produzir uma biblioteca de cDNAs. Esse artigo descreve o sequenciamento parcial de 356 transcritos das glândulas de S. plumieri. Análises em bancos de dados (BLAST) dos transcritos demonstraram que 30% eram sequências desconhecidas, 17% proteínas hipotéticas, 17% relacionadas às enzimas do metabolismo, 14% pertenciam a funções de transdução de sinais e os demais grupos (7-8%) formados por genes relacionados com a expressão de proteínas, proteínas regulatórias e estruturais. Um número considerável destes EST não foi encontrado em bases de dados disponíveis, sugerindo a existência de novas proteínas e/ou funções ainda a serem descobertas. Ao fazer um barrido da biblioteca com anticorpos produzidos contra uma fração das lectinas do veneno de S. plumieri, identificamos vários clones, cuja sequência de DNA mostram semelhanças com lectinas encontradas em peixes. A análise in silico destes clones confirmam a identidade destas moléculas na glândula de peçonha de S. plumieri.

Toxin gene determination and evolution in scorpaenoid fish.

In this study, we determine the toxin genes from both cDNA and genomic DNA of four scorpaenoid fish and reconstruct their evolutionary relationship. The deduced protein sequences of the two toxin subunits in Sebastapistes strongia, Scorpaenopsis oxycephala, and Sebastiscus marmoratus are about 700 amino acid, similar to the sizes of the stonefish (Synanceia horrida, and Synanceia verrucosa) and lionfish (Pterois antennata and Pterois volitans) toxins previously published. The intron positions are highly conserved among these species, which indicate the applicability of gene finding by using genomic DNA template. The phylogenetic analysis shows that the two toxin subunits were duplicated prior to the speciation of Scorpaenoidei. The precedence of the gene duplication over speciation indicates that the toxin genes may be common to the whole family of Scorpaeniform. Furthermore, one additional toxin gene has been determined in the genomic DNA of Dendrochirus zebra. The phylogenetic analysis suggests that an additional gene duplication occurred before the speciation of the lionfish (Pteroinae) and a pseudogene may be generally present in the lineage of lionfish.

DNA microarray analysis on gene candidates possibly related to tetrodotoxin accumulation in pufferfish.

Pufferfish accumulate tetrodotoxin (TTX) at high levels in liver and ovary through the food chain. However, the mechanisms underlying TTX toxification in pufferfish have been poorly understood. In order to search gene candidates involved in TTX accumulation in the torafugu pufferfish Takifugu rubripes, a custom 4x44k oligonucleotide microarray slide was designed by the Agilent eArray program using oligonucleotide probes of 60 bp in length referring to 42,724 predicted transcripts in the publicly available Fugu genome database. DNA microarray analysis was performed with total RNA samples from the livers of two toxic wild specimens in comparison with those from a nontoxic wild specimen and two nontoxic cultured specimens. The mRNA levels of 1108 transcripts were more than 2-fold higher in the toxic specimens than in the nontoxic specimens. The levels of 613 transcripts were remarkably high, and 16 transcripts encoded by 9 genes were up-regulated more than 10-fold. These genes included those encoding forming structural filaments (keratins) and those related to vitamin D metabolism and immunity. It was also noted that the levels of the transcripts encoding serpin peptidase inhibitor clade C member 1, coagulation factor X precursor, complement C2, C3, C5, C8 precursors, and interleukin-6 receptor were high in the toxic liver samples.

Venom evolution through gene duplications.

Venoms contain highly complex mixtures that typically include hundreds of different components and have evolved independently in a diverse range of animals including platypuses, shrews, snakes, lizards, fishes, echinoderms, spiders, wasps, centipedes, sea snails, cephalopods, jellyfish and sea anemones. Many venom genes evolved through gene duplication. Gene duplication occurs in all domains of life and provides the raw substrate from which novel function arise. In this review, we focus on the role that gene duplication has played in the origin and diversification of venom genes. We outline the selective advantages of venom gene duplicates and the role that selection has played in the retention of these duplicates. We use toxin gene intermediates to help trace the evolution of toxin innovation. We also focus on other genomic processes, such as exon and domain duplications, in venom evolution. Finally, we conclude by focusing on the use of high throughput sequencing technology in understanding venom evolution.

Properties and cDNA cloning of a hyaluronidase from the stonefish Synanceia verrucosa venom.

The venoms of two classes of fish, freshwater stingray (members of the genus Potamotrygon) and stonefish (members of the genus Synanceia), contain not only proteinaceous toxins but also hyaluronidases, which are considered as spreading factors that facilitate the tissue diffusion of toxins by degrading hyaluronan. So far, the freshwater stingray Potamotrygon motoro hyaluronidase and the stonefish Synanceia horrida hyaluronidase (SFHYA1) have been purified and characterized, although their spreading activity is still unclear. In this study, a 59 kDa hyaluronidase was partially purified from the stonefish Synanceia verrucosa and shown to be optimally active at pH 6.6, 37 °C and 0.15 M NaCl. Importantly, the S. verrucosa hyaluronidase enhanced the capillary permeability-increasing activity of the S. verrucosa toxin (neoverrucotoxin), providing evidence for its spreading activity. Furthermore, the primary structure of the S. verrucosa hyaluronidase was elucidated by cDNA cloning. The S. verrucosa hyaluronidase (463 amino acid residues) shares as high as 92% sequence identity with SFHYA1 but less than 50% with other hyaluronidases. Nevertheless, one catalytic residue and four substrate positioning residues, which constitute the active site of human hyaluronidases, are conserved in the S. verrucosa hyaluronidase.

The action of fish peptide Orpotrin analogs on microcirculation

In order to investigate the relationshipbetween theprimary structure of Orpotrin, a vasoactive peptidepreviously isolated from the freshwater stingray Potamotrygon gr. orbignyi, and its microcirculatory effects, three Orpotrin analogs were synthesized. The analogs have a truncated N-terminal with a His residue deletion and two substituted amino acid residues, where one Nle is substituted for one internal Lys residue and the third analog has a substitution of a Pro for an Ala (Orp-desH1, Orp- Nle and Orp-Pro/Ala, respectively). Only Orp-desH1 could induce a lower vasoconstriction effect compared with the natural Orpotrin, indicating that besides the N-terminal, the positive charge of Lys and the Pro residues located at the center of the amino acid chain is crucial for this vasoconstriction effect. Importantly, the suggestions made with bioactive peptides were based on the molecular modeling and dynamics of peptides, the presence of key amino acids and shared activity in microcirculation, characterized by intravital microscopy. Moreover, this study has demonstrated that even subtle changes in the primary structure of Orpotrin alter the biological effects of this native peptide significantly, which could be of interest for biotechnological applications.

Purification, properties and cDNA cloning of neoverrucotoxin (neoVTX), a hemolytic lethal factor from the stonefish Synanceia verrucosa venom.

A proteinaceous toxin with hemolytic and lethal activities, named neoverrucotoxin (neoVTX), was purified from the venom fluid of stonefish Synanceia verrucosa and its primary structure was elucidated by a cDNA cloning technique. NeoVTX is a dimeric 166 kDa protein composed of alpha-subunit (702 amino acid residues) and beta-subunit (699 amino acid residues) and lacks carbohydrate moieties. Its hemolytic activity is inhibited by anionic lipids, especially potently by cardiolipin. These properties are comparable to those of stonustoxin (SNTX) previously purified from S. horrida. Alignment of the amino acid sequences also reveals that the neoVTX alpha- and beta-subunits share as high as 87 and 95% sequence identity with the SNTX alpha- and beta-subunits, respectively. The distinct differences between neoVTX and SNTX are recognized only in the numbers of Cys residues (18 for neoVTX and 15 for SNTX) and free thiol groups (10 for neoVTX and 5 for SNTX). In contrast, neoVTX considerably differs from verrucotoxin (VTX), a tetrameric 322 kDa glycoprotein, previously purified from S. verrucosa. In addition, the sequence identity of the neoVTX beta-subunit with the reported VTX beta-subunit is 90%, being lower than that with the SNTX beta-subunit.

Transcriptome analysis of expressed sequence tags from the venom glands of the fish Thalassophryne nattereri.

Thalassophryne nattereri (niquim) is a venomous fish found on the northern and northeastern coasts of Brazil. Every year, hundreds of humans are affected by the poison, which causes excruciating local pain, edema, and necrosis, and can lead to permanent disabilities. In experimental models, T. nattereri venom induces edema and nociception, which are correlated to human symptoms and dependent on venom kininogenase activity; myotoxicity; impairment of blood flow; platelet lysis and cytotoxicity on endothelial cells. These effects were observed with minute amounts of venom. To characterize the primary structure of T. nattereri venom toxins, a list of transcripts within the venom gland was made using the expressed sequence tag (EST) strategy. Here we report the analysis of 775 ESTs that were obtained from a directional cDNA library of T. nattereri venom gland. Of these ESTs, 527 (68%) were related to sequences previously described. These were categorized into 10 groups according to their biological functions. Sequences involved in gene and protein expression accounted for 14.3% of the ESTs, reflecting the important role of protein synthesis in this gland. Other groups included proteins engaged in the assembly of disulfide bonds (0.5%), chaperones involved in the folding of nascent proteins (1.4%), and sequences related to clusterin (1.5%), as well as transcripts related to calcium binding proteins (1.0%). We detected a large cluster (1.3%) related to cocaine- and amphetamine-regulated transcript (CART), a peptide involved in the regulation of food intake. Surprisingly, several retrotransposon-like sequences (1.0%) were found in the library. It may be that their presence accounts for some of the variation in venom toxins. The toxin category (18.8%) included natterins (18%), which are a new group of kininogenases recently described by our group, and a group of C-type lectins (0.8%). In addition, a considerable number of sequences (32%) was not related to sequences in the databases, which indicates that a great number of new toxins and proteins are still to be discovered from this fish venom gland.

Transcriptome analysis of expressed sequence tags from the venom glands of the fish Thalassophryne nattereri

Thalassophryne nattereri (niquim) is a venomous fish found on the northern and northeastern coasts of Brazil. Every year, hundreds of humans are affected by the poison, which causes excruciating local pain, edema, and necrosis, and can lead to permanent disabilities. In experimental models, T. nattereri venom induces edema and nociception, which are correlated to human symptoms and dependent on venom kininogenase activity; myotoxicity; impairment of blood flow; platelet lysis and cytotoxicity on endothelial cells. These effects were observed with minute amounts of venom. To characterize the primary structure of T. nattereri venom toxins, a list of transcripts within the venom gland was made using the expressed sequence tag (EST) strategy. Here we report the analysis of 775 ESTs that were obtained from a directional cDNA library of T. nattereri venom gland.

Fish retroposons related to the Penelope element of Drosophila virilis define a new group of retrotransposable elements.

Poseidon and Neptune are two ancient lineages of retroposons related to the Penelope element from Drosophila virilis. They have been identified in various teleost fish species, including the medakafish (Oryzias latipes), and the pufferfishes Fugu rubripes and Tetraodon nigroviridis, whose genomes are currently being sequenced. Some of these elements are highly reiterated in fish genomes. Penelope-related elements were also identified in blood fluke, shrimp, sea urchin, cichlid fish and frog, showing that they are widespread in animals. Penelope-related retroposons were not detected among sequences from the Drosophila melanogaster and human genome projects, suggesting that they have been lost from certain animal lineages. A sequence encoding a putative Uri (also called GIY-YIG) endonuclease domain was detected downstream from the gene for reverse transcriptase. To the best of our knowledge, this type of endonuclease sequence has previously been identified in group I introns and in genes for prokaryotic excinucleases but not in retrotransposable elements. Penelope-related elements are frequently truncated at their 5' ends and can also be flanked by long terminal repeat-like structures. Phylogenetic analysis of the reverse transcriptase domain failed to assign Penelope-related retroposons to one of the major groups of retroelements. Overall, therefore, the evidence strongly suggests that these sequences represent a new group of retrotransposable elements.

Comparative analysis of the ETV6 gene in vertebrate genomes from pufferfish to human.

The ETV6 gene encodes an Ets-like transcription factor that is frequently rearranged in leukemias. While some of the functions of ETV6 have been uncovered recently, little is known about the key structural elements involved. Comparative genome analysis may provide novel insights into gene evolution and functions. In this study, we cloned and sequenced the homologue of ETV6 from the compact genome of the pufferfish Fugu rubripes (fETV6). The genomic structure of the fETV6 gene was investigated by sequence analysis of a contig of genomic clones. The fETV6 gene, composed of eight exons, spans about 15 kb and is 16 times smaller than its human counterpart mainly because of the reduced intron size. Three of the seven introns of fETV are unusually large (more than 2 kb), including the 8.2 kb intron 2. The gene codes for a protein of 465 amino acids that is highly related to its human homologue, exhibiting an overall identity of 58% (72% similarity). To investigate the functional and evolutionary aspects of ETV6, we undertook a comparative analysis of this gene from various vertebrates (human, mouse, chicken, zebrafish and Fugu). As expected, the PNT and ETS domains were highly conserved, with on average 81 and 95% peptide sequence identity, respectively. In addition, we found several new highly conserved regions within the central section of the protein that are likely to represent further functional or structural domains, which may be associated with the transcription repression capacity of this protein. We also found conserved putative regulatory elements in the promoter as well as in the large intron 2 of fETV6. The information derived from this comparative analysis will serve as the basis for more precise functional studies of ETV6 gene regulation and function.

Characterization of the Japanese pufferfish (Takifugu rubripes) T-cell receptor alpha locus reveals a unique genomic organization.

Polymerase chain reactions with degenerate V gene segment primers were used to isolate the putative T-cell receptor alpha-chain gene (TCRA) from Japanese pufferfish (Takifugu rubripes). The putative TCRA chain cDNA is composed of an N-terminus leader peptide followed by the variable region and the constant region. The variable portion of the TCRA gene is encoded by V and J gene segments separated in the germline. As in mammals, the V-J junction sequences are GC rich and highly diversified. Amino acid residues that are required to maintain the function and structural integrity of the TCRA polypeptide, including the conserved Trp-Tyr-Lys and Tyr-Tyr-Cys motifs in the V gene segments, the Lys-Leu-X-Phe-Gly-X-Gly-Thr-X-Leu motif in the J gene segment, the three cysteine residues in the constant region and the charged residues in the transmembrane region are all preserved in the pufferfish. These conserved features suggest that the TCRA gene families in fish and mammals have evolved from a common ancestor.

Genomic sequence analysis of Fugu rubripes CFTR and flanking genes in a 60 kb region conserving synteny with 800 kb of human chromosome 7.

To define control elements that regulate tissue-specific expression of the cystic fibrosis transmembrane regulator (CFTR), we have sequenced 60 kb of genomic DNA from the puffer fish Fugu rubripes (Fugu) that includes the CFTR gene. This region of the Fugu genome shows conservation of synteny with 800-kb sequence of the human genome encompassing the WNT2, CFTR, Z43555, and CBP90 genes. Additionally, the genomic structure of each gene is conserved. In a multiple sequence alignment of human, mouse, and Fugu, the putative WNT2 promoter sequence is shown to contain highly conserved elements that may be transcription factor or other regulatory binding sites. We have found two putative ankyrin repeat-containing genes that flank the CFTR gene. Overall sequence analysis suggests conservation of intron/exon boundaries between Fugu and human CFTR and revealed extensive homology between functional protein domains. However, the immediate 5' regions of human and Fugu CFTR are highly divergent with few conserved sequences apart from those resembling diminished cAMP response elements (CRE) and CAAT box elements. Interestingly, the polymorphic polyT tract located upstream of exon 9 is present in human and Fugu but absent in mouse. Similarly, an intron 1 and intron 9 element common to human and Fugu is absent in mouse. The euryhaline killifish CFTR coding sequence is highly homologous to the Fugu sequence, suggesting that upregulation of CFTR in that species in response to salinity may be regulated transcriptionally.

Characterization and repeat analysis of the compact genome of the freshwater pufferfish Tetraodon nigroviridis.

Tetraodon nigroviridis is a freshwater pufferfish 20-30 million years distant from Fugu rubripes. The genome of both tetraodontiforms is compact, mostly because intergenic and intronic sequences are reduced in size compared to other vertebrate genomes. The previously uncharacterized Tetraodon genome is described here together with a detailed analysis of its repeat content and organization. We report the sequencing of 46 megabases of bacterial artificial chromosome (BAC) end sequences, which represents a random DNA sample equivalent to 13% of the genome. The sequence and location of rRNA gene clusters, centromeric and subtelocentric satellite sequences have been determined. Minisatellites and microsatellites have been cataloged and notable differences were observed in comparison with microsatellites from Fugu. The genome contains homologies to all known families of transposable elements, including Ty3-gypsy, Ty1-copia, Line retrotransposons, DNA transposons, and retroviruses, although their overall abundance is <1%. This structural analysis is an important prerequisite to sequencing the Tetraodon genome.

Random expression of main and vomeronasal olfactory receptor genes in immature and mature olfactory epithelia of Fugu rubripes.

Main olfactory receptor genes were isolated from a seawater fish, Fugu rubripes (pufferfish), and characterized. Two subfamilies of genes encoding seven transmembrane receptors were identified; one consists of five or more members, termed FOR1-1 to 5 of FOR1 subfamily, and the other appears to be a single copy gene, termed the FOR2 subfamily. FOR1 members show extremely high amino acid sequence similarities of about 95% to one another, and are distantly related to catfish-1 with the highest similarity of 37%. FOR2 shows 43% similarity to goldfish-A28. Phylogenically, both FOR members are categorized among pedigrees of the fish main olfactory receptor family outside the mammalian receptor family, although similarities between Fugu receptors and those of fresh-water fishes are lower than those among fresh-water fishes. In situ hybridization shows that both subfamilies of receptor genes are expressed randomly over the olfactory epithelium throughout all developmental stages, and no segregation of the signals was found. On the other hand, when three members of a vomeronasal olfactory receptor gene family, related to the Ca(2+)-sensing receptor, were used as probes, they were also randomly expressed over the same epithelium as the main olfactory receptors. This is in contrast to the expression profiles observed for zebrafish and goldfish, where the main or vomeronasal olfactory receptors are expressed in segregated patterns. It is thus suggested that the expression pattern of fish olfactory receptors varies depending on the species, although fish olfactory receptors are highly related to one another in their primary structures, and are phylogenically distinct from those of mammals.

Generation and analysis of 25 Mb of genomic DNA from the pufferfish Fugu rubripes by sequence scanning.

We have generated and analyzed >50,000 shotgun clones from 1059 Fugu cosmid clones. All sequences have been minimally edited and searched against protein and DNA databases. These data are all displayed on a searchable, publicly available web site at. With an average of 50 reads per cosmid, this is virtually nonredundant sequence skimming, covering 30%-50% of each clone. This essentially random data set covers nearly 25 Mb (>6%) of the Fugu genome and forms the basis of a series of whole genome analyses which address questions regarding gene density and distribution in the Fugu genome and the similarity between Fugu and mammalian genes. The Fugu genome, with eight times less DNA but a similar gene repertoire, is ideally suited to this type of study because most cosmids contain more than one identifiable gene. General features of the genome are also discussed. We have made some estimation of the syntenic relationship between mammals and Fugu and looked at the efficacy of ORF prediction from short, unedited Fugu genomic sequences. Comparative DNA sequence analyses are an essential tool in the functional interpretation of complex vertebrate genomes. This project highlights the utility of using the Fugu genome in this kind of study.

Molecular cloning of EDG-3 and N-Shc genes from the puffer fish, Fugu rubripes, and conservation of synteny with the human genome.

EDG-3 is a receptor for sphingosine-1-phosphate mapped on human chromosome 9q22.1-q22.2. We used the compact Fugu genome for its linkage analysis. The Fugu EDG-3 was composed of one intron and two exons, encoding a 384 amino acid protein that has 56.9% homology with the human EDG-3. Approximately 3 kb apart, a neuronal Shc (N-Shc) gene was identified. It spans 7 kb containing 12 coding exons, and has an overall 53.4% similarity with the human protein. We mapped the human N-Shc gene to chromosome 9q21.3-q22.2. This is the first report of the genomic structure and the linkage of these two genes conserved between Fugu and human.

Genomic structure and comparative analysis of nine Fugu genes: conservation of synteny with human chromosome Xp22.2-p22.1.

The pufferfish Fugu rubripes has a compact 400-Mb genome that is approximately 7.5 times smaller than the human genome but contains a similar number of genes. Focusing on the distal short arm of the human X chromosome, we have studied the evolutionary conservation of gene orders in Fugu and man. Sequencing of 68 kb of Fugu genomic DNA identified nine genes in the following order: (SCML2)-STK9, XLRS1, PPEF-1, KELCH2, KELCH1, PHKA2, AP19, and U2AF1-RS2. Apart from an evolutionary inversion separating AP19 and U2AF1-RS2 from PHKA2, gene orders are identical in Fugu and man, and all nine human homologs map to the Xp22 band. All Fugu genes were found to be smaller than their human counterparts, but gene structures were mostly identical. These data suggest that genomic sequencing in Fugu is a powerful and economical strategy to predict gene orders in the human genome and to elucidate the structure of human genes.