HSP70 gene expression responses to the temperature stress in pufferfish (Takifugu rubripes).

In this study, we isolated and characterized HSP70 cDNA from pufferfish (Takifugu rubripes). The 3053 bp full-length TrHSP70 sequence consisted of a 167 bp 5'-UTR (untranslated region), a 2535 bp open reading frame, and a 351 bp 3'-UTR. BLAST analysis revealed that the TrHSP70 shared high similarity with HSP70 sequences in other species. In our study, we set 3 experimental groups as H1 group (20 °C), H2 group (24 °C), and H3 group (28 °C) for checking the expression level of TrHSP70 in T. rubripes. Tissue-specific gene expression results showed that TrHSP70 had higher expression in the intestines than other tissues of the T. rubripes by RT-qPCR. In the experimental group, we found that the expression of TrHSP70 was upregulated in different tissues in the H3 group. The results show that TrHSP70 is a constitutively expressed gene, which plays an important role in maintaining normal physiological function and coping with stress.

Characterization of the torafugu (Takifugu rubripes) immunoglobulin heavy chain gene locus.

In this study, we investigated the immunoglobulin heavy (IGH) gene locus of torafugu (Takifugu rubripes) from publicly available assembly sequences and presented an annotated locus map, including the IGHV genes, pseudogenes, and IGHC genes. Three new IGHV gene families (IGHV3-IGHV5) were discovered. We observed the interspersion of IGHV1 and IGHV2 family members and that they often intermingled with each other, while other family members were further interspersed. Conservation of the promoter and recombination signal sequences (RSS) was observed in a family-specific manner. In addition to known variable region genes present on chromosome 5 (current torafugu genome assembly), we found 34 additional IGHV genes on scaffold 287 and three novel potentially functional IGHD genes on scaffold 483. In total, the variable region of the torafugu IGH locus consists of at least 48 IGHV genes, seven IGHD genes, and six IGHJ genes. IGHC genes have also been mapped in this study, with three genes encoding immunoglobulin classes: IgT, IgM, and IgD. We confirmed the expression of newly identified IGHV3 family sequences in the spleen and kidney of adult torafugu and found a favorable IGHV segment usage by IgM and IgT. Possible structural variation in the IGHδ locus was observed based on the current torafugu assembly. The complete characterization of the torafugu IGH locus will facilitate detailed studies of large-scale mechanisms associated with the recombination of the variable region genes and will offer insights into the genetic basis of the potential diversity in the antibody response observed in torafugu.

Characterization and functional analysis of two PKR genes in fugu (Takifugu rubripes).

PKR (protein kinase R) is a serine-threonine kinase that inhibits protein synthesis by the phosphorylation of the eukaryotic translation initiation factor 2-alpha (eIF2α), and activates NFκB by inducing NFκB-inducing kinase and IκB (inhibitor of NFκB) kinase. This can lead to antiviral and anti-proliferative effects. In this study, the complete sequence and organization of two fugu PKR genes (fPKRs) were determined by in silico analysis and conventional PCR. The full-length fPKR1 and fPKR2 genes were 3832 bp and 4325 bp, which encoded 523 and 492 amino acids, respectively. Both encoded two dsRNA binding domains and a Serine/Threonine protein kinase domain, and showed very high similarity to green spotted puffer PKRs. Gene expression of the two fPKRs was measured by quantitative real-time PCR on tissue samples from healthy fish and peripheral blood leukocytes stimulated with polyinosinic:polycytidylic acid (PolyI:C) or lipopolysaccharides (LPS). The fPKRs were highly expressed in the skin and fPKR2 was significantly induced in PBLs by PolyI:C but not by LPS. The fPKRs inhibited translation of a luciferase reporter gene in a dose-dependent manner and induced transcriptional activity of a mammalian NFκB luciferase reporter. These results demonstrate that two PKRs in a single species can both be independently, but not equally, functional and support the hypothesis that fish PKRs have roles in the innate immune response similar to those of mammalian PKRs.

Species discrimination among three kinds of puffer fish using an electronic nose combined with olfactory sensory evaluation.

Species discrimination among three kinds of puffer fish, Takifugu obscurus, Takifugu flavidus and Takifugu rubripes, was conducted using an electronic nose combined with olfactory sensory evaluation. All data were treated by multivariate data processing based on principal component analysis (PCA) and discriminant factor analysis (DFA). The results showed the discriminant model by PCA method and DFA method. Using PCA and DFA, it was shown that the electronic nose was able to reasonably distinguish between each of the eleven puffer fish groups, with a discrimination index of 85. The olfactory sensory evaluation was undertaken in accordance to Sensory analysis-Methodology-Initiation and training of assessors in the detection and recognition of odors (BS ISO 5496-2006), and the results showed that the evaluation was able to identify puffer fish samples according to their species, geographical origin and age. Results from this analysis demonstrate that the E-nose can be used to complement the discrimination of odors by sensory evaluation from the three species of puffer fish studied here.

Different phylogenomic approaches to resolve the evolutionary relationships among model fish species.

Comparative genomics holds the promise to magnify the information obtained from individual genome sequencing projects, revealing common features conserved across genomes and identifying lineage-specific characteristics. To implement such a comparative approach, a robust phylogenetic framework is required to accurately reconstruct evolution at the genome level. Among vertebrate taxa, teleosts represent the second best characterized group, with high-quality draft genome sequences for five model species (Danio rerio, Gasterosteus aculeatus, Oryzias latipes, Takifugu rubripes, and Tetraodon nigroviridis), and several others are in the finishing lane. However, the relationships among the acanthomorph teleost model fishes remain an unresolved taxonomic issue. Here, a genomic region spanning over 1.2 million base pairs was sequenced in the teleost fish Dicentrarchus labrax. Together with genomic data available for the above fish models, the new sequence was used to identify unique orthologous genomic regions shared across all target taxa. Different strategies were applied to produce robust multiple gene and genomic alignments spanning from 11,802 to 186,474 amino acid/nucleotide positions. Ten data sets were analyzed according to Bayesian inference, maximum likelihood, maximum parsimony, and neighbor joining methods. Extensive analyses were performed to explore the influence of several factors (e.g., alignment methodology, substitution model, data set partitions, and long-branch attraction) on the tree topology. Although a general consensus was observed for a closer relationship between G. aculeatus (Gasterosteidae) and Di. labrax (Moronidae) with the atherinomorph O. latipes (Beloniformes) sister taxon of this clade, with the tetraodontiform group Ta. rubripes and Te. nigroviridis (Tetraodontiformes) representing a more distantly related taxon among acanthomorph model fish species, conflicting results were obtained between data sets and methods, especially with respect to the choice of alignment methodology applied to noncoding parts of the genomic region under study. This may limit the use of intergenic/noncoding sequences in phylogenomics until more robust alignment algorithms are developed.

High-Resolution Melting Analysis of COI Sequences Distinguishes Pufferfish Species (Takifugu spp.) in China.

Pufferfish is a traditional, delicious dish in Asia. However, eating wild or improperly processed pufferfish causes serious poisoning. This study aimed to exploit the high-resolution melting (HRM) method for authenticating four species of Takifugu pufferfish (Takifugu xanthopterus, T. fasciatus, T. flavidus, and T. rubripes). Candidate DNA barcodes, including the cytochrome c oxidase subunit I (COI), cytochrome oxidase b (Cytb), and the control region (D-loop), were analyzed, with COI selected as the optimal DNA barcode. An HRM method was developed to identify 57 commercial fish samples in China, including 33 commercial pufferfish products and 24 unlabeled fish products. The findings revealed that the pufferfish products were T. rubripes or T. fasciatus, and four T. xanthopterus samples were detected in unlabeled fish products. These results showed that DNA barcode coupled with HRM analysis was a rapid and efficient tool to identify pufferfish, which might aid in the prevention of consumer fraud or mislabeling of fish products.

Forkhead transcription factor O1 (FoxO1) in torafugu pufferfish Takifugu rubripes: Molecular cloning, in vitro DNA binding, and target gene screening in fish metagenome.

The fish insulin/insulin-like growth factor (IGF) pathway has weak control over carbohydrate metabolism. To understand the molecular basis for the metabolic diversity, we characterized the forkhead box transcription factor O1A (FoxO1A), a downstream target of the insulin/IGF pathway, in torafugu Takifugu rubripes. The cloned torafugu FoxO1A cDNA contained all conserved features critical for its transcriptional activity and a unique unspliced intron encoding a poly-glutamine stretch. Torafugu FoxO1A showed the IGF-dependent nuclear exclusion and in vitro binding to the well-conserved FoxO1 binding site, DAF-16 binding element (DBE), but failed to bind to the insulin-responsive element by which mammalian FoxO1 mediates insulin effects. The subsequent in silico genomic screening provided a list of 587 potential torafugu FoxO1A target genes containing the DBE. Some carbohydrate metabolic genes regulated by FoxO1 in mammals were not included in the list. We further identified about 250 potential fish FoxO1 target genes by integrating results of the DBE screening against fish metagenome that contained 262 species. Neuronal processes appeared to be the common major function of fish FoxO1, although further annotation of the potential target genes is required. These results provide a part of the molecular basis underlying the weak association between the insulin/IGF pathway and carbohydrate metabolism in fish.

Classification of Takifugu rubripes, T. chinensis and T. pseudommus by genotyping-by-sequencing.

Takifugu rubripes is more expensive than other species of the genus because of its high protein content and special flavor. However, it is easily confused with imported T. chinensis and T. pseudommus because they have similar morphological characteristics. We identified single nucleotide polymorphism (SNP) markers of T. rubripes by genotyping-by-sequencing (GBS) and evaluated their ability to distinguish among T. rubripes, T. chinensis, and T. pseudommus. In all, 18 polymorphic SNPs were subjected to phylogenetic analyses of the three Takifugu species. Additionally, we subjected a second set of samples to Sanger sequencing to verify that the polymorphic SNPs could be used to evaluate the genetic variation among the three Takifugu species. A phylogenetic tree that included the analyzed sequence of set A, which is referred to as the reference sequence, and a validation sequence of set B with 18 SNPs were produced. Based on this phylogenetic tree and STRUCTURE analyses, T. rubripes, T. chinensis and T. pseudommus have low genetic variation and should be considered the same gene pool. Our findings suggest that further studies are needed to estimate the genetic association of the three Takifugu species.

Chromosomal-level assembly of Takifugu obscurus (Abe, 1949) genome using third-generation DNA sequencing and Hi-C analysis.

The Tetraodontidae family are known to have relatively small and compact genomes compared to other vertebrates. The obscure puffer fish Takifugu obscurus is an anadromous species that migrates to freshwater from the sea for spawning. Thus the euryhaline characteristics of T. obscurus have been investigated to gain understanding of their survival ability, osmoregulation, and other homeostatic mechanisms in both freshwater and seawater. In this study, a high quality chromosome-level reference genome for T. obscurus was constructed using long-read Pacific Biosciences (PacBio) Sequel sequencing and a Hi-C-based chromatin contact map platform. The final genome assembly of T. obscurus is 381 Mb, with a contig N50 length of 3,296 kb and longest length of 10.7 Mb, from a total of 62 Gb of raw reads generated using single-molecule real-time sequencing technology from a PacBio Sequel platform. The PacBio data were further clustered into chromosome-scale scaffolds using a Hi-C approach, resulting in a 373 Mb genome assembly with a contig N50 length of 15.2 Mb and and longest length of 28 Mb. When we directly compared the 22 longest scaffolds of T. obscurus to the 22 chromosomes of the tiger puffer Takifugu rubripes, a clear one-to-one orthologous relationship was observed between the two species, supporting the chromosome-level assembly of T. obscurus. This genome assembly can serve as a valuable genetic resource for exploring fugu-specific compact genome characteristics, and will provide essential genomic information for understanding molecular adaptations to salinity fluctuations and the evolution of osmoregulatory mechanisms.

Molecular cloning and beta-naphthoflavone-induced expression of a cytochrome P450 1A (CYP1A) gene from an anadromous river pufferfish, Takifugu obscurus.

In recent years, there has been a decline in the wild populations of river pufferfish, Takifugu obscurus. Besides overexploitation for commercial purposes, environmental pollution is believed to have contributed to its decline. However, almost no information exists about genes involved in metabolism of xenobiotics by this species. Nevertheless, there is interest in fugu fishes, since they possess the smallest genome among vertebrates. We cloned and characterized the full-length cDNA sequence of a cytochrome P450 1A (CYP1A) gene from T. obscurus. Phylogenic relationship of T. obscurus CYP1A was also compared to other fish species. The tissue distribution and time-dependant induction of CYP1A mRNA were studied by real-time PCR in T. obscurus exposed to an aryl hydrocarbon receptor (Ahr) agonist, beta-naphthoflavone (BNF). The greatest basal expression in livers of control as well as BNF-treated individuals. However, brain, gill, gonad, intestine, and kidney also expressed CYP1A. Muscles expressed the least CYP1A. The results of the time-course study revealed induction in brain and gills after 6h and at 12h in most tissues. Except for gills, all other organs retained induced expression of CYP1A mRNA up to 96h.

Artificially designed hybrids facilitate efficient generation of high-resolution linkage maps.

When sequencing eukaryotic genomes, linkage maps are indispensable for building scaffolds to assemble and/or to validate chromosomes. However, current approaches to constructing linkage maps are limited by marker density and cost-effectiveness, especially for wild organisms. We have now devised a new strategy based on artificially generated hybrid organisms to acquire ultrahigh-density genomic markers at reduced cost and build highly accurate linkage maps. We have also developed the novel analysis pipeline Scaffold Extender with Low Depth Linkage Analysis (SELDLA) for data processing to generate linkage maps and draft genomes. Using SELDLA, linkage maps and improved genomes for two species of pufferfish, Takifugu rubripes and Takifugu stictonotus, were obtained simultaneously. The strategy is applicable to a wide range of sexually reproducing organisms, and could, therefore, accelerate the whole genome analysis of various organisms including fish, mollusks, amphibians, insects, plants, and even mammals.

Divergent evolution of the myosin heavy chain gene family in fish and tetrapods: evidence from comparative genomic analysis.

Myosin heavy chain genes (MYHs) are the most important functional domains of myosins, which are highly conserved throughout evolution. The human genome contains 15 MYHs, whereas the corresponding number in teleost appears to be much higher. Although teleosts comprise more than one-half of all vertebrate species, our knowledge of MYHs in teleosts is rather limited. A comprehensive analysis of the torafugu (Takifugu rubripes) genome database enabled us to detect at least 28 MYHs, almost twice as many as in humans. RT-PCR revealed that at least 16 torafugu MYH representatives (5 fast skeletal, 3 cardiac, 2 slow skeletal, 1 superfast, 2 smooth, and 3 nonmuscle types) are actually transcribed. Among these, MYH(M743-2) and MYH(M5) of fast and slow skeletal types, respectively, are expressed during development of torafugu embryos. Syntenic analysis reveals that torafugu fast skeletal MYHs are distributed across five genomic regions, three of which form clusters. Interestingly, while human fast skeletal MYHs form one cluster, its syntenic region in torafugu is duplicated, although each locus contains just a single MYH in torafugu. The results of the syntenic analysis were further confirmed by corresponding analysis of MYHs based on databases from Tetraodon, zebrafish, and medaka genomes. Phylogenetic analysis suggests that fast skeletal MYHs evolved independently in teleosts and tetrapods after fast skeletal MYHs had diverged from four ancestral MYHs.

[Characterization and phylogenetic analysis of the cytochrome oxidase subunit I gene of mitochondrial genome from Takifugu fasciatus].

The Takifugu fasciatus mitochondrial cytochrome oxidase I gene (COI) and its associated tRNA genes were sequenced by PCR. The open reading frame of COI gene contains 1,546 bp nucleotides, encoding a putative protein of 515 amino acid residues. The pattern of codon usage of the COI gene is less biased toward A+T. The COI gene of T. fasciatus shows a high degree of homology with that from the other 14 fish species recorded in the GenBank and has 97.6% homology with Takifugu rubripes, 76.5% with Masturus lanceolatus and 75.4% with Mola mola. The phylogenetic trees show that the relationships based on the homology is consistent with the morphological and taxonomic results. Predicted secondary structures of the tRNA genes suggest that they have the classical cloverleaf structures.

Identification and expression of a novel carbonic anhydrase isozyme in the pufferfish Takifugu vermicularis.

Carbonic anhydrase (CA) is a key element for maintaining acid base balance in fish. In our present experiment, novel CA isozymes were identified from the pear puffer (Takifugu vermicularis). Based on the high homology of two predicted CA sequences of the tiger puffer (Takifugu rubripes), a 1715bp novel cDNA was obtained from T. vermicularis. The open reading frame showed a complete coding sequence of 552bp with a deduced peptide sequence of 183 amino acids that exhibited highest (97%) identity with pufferfish putative CA III and CA IV-like sequences. In addition, this translated protein sequence showed 36-37% identity with zebrafish CA IV-like, CA XVa, CA XVb, and CA XVc proteins. Phylogenetic analysis revealed that the pufferfish novel protein (pCAn) was a membrane-bound CA protein. Alignment of multiple CA sequences illustrated that most of the putative active site residues of the pCAn isozyme were situated at highly conserved regions of the CA sequences. Examination of motif distribution suggested that the pCAn isozyme was very similar to the puffer predicted CA IV-like isozyme. Reverse transcription-polymerase chain reaction (PCR) analysis showed highly differential expression in the brain, gills, kidney, and muscle, whereas CA mRNA expression was almost absent in heart, liver, and intestine. Quantitative PCR expression of CA mRNA abundance suggested several-fold higher expression of pCAn isozymes in the gills compared to other tissues tested. Our results suggest that the pCAn isozyme might be related to CA IV-like isozymes. Further functional studies are needed to investigate the function of the pCAn isozyme in T. vermicularis.

Takifugu obscurus is a euryhaline fugu species very close to Takifugu rubripes and suitable for studying osmoregulation.

BACKGROUND: The genome sequence of the pufferfish Takifugu rubripes is an enormously useful tool in the molecular physiology of fish. Euryhaline fish that can survive both in freshwater (FW) and seawater (SW) are also very useful for studying fish physiology, especially osmoregulation. Recently we learned that there is a pufferfish, Takifugu obscurus, common name "mefugu" that migrates into FW to spawn. If T. obscurus is indeed a euryhaline fish and shares a high sequence homology with T. rubripes, it will become a superior animal model for studying the mechanism of osmoregulation. We have therefore determined its euryhalinity and phylogenetic relationship to the members of the Takifugu family. RESULTS: The following six Takifugu species were used for the analyses: T. obscurus, T. rubripes, T. niphobles, T. pardalis, T. poecilonotus, and T. porphyreus. When transferred to FW, only T. obscurus could survive while the others could not survive more than ten days in FW. During this course of FW adaptation, serum Na+ concentration of T. obscurus decreased only slightly, but a rapid and large decrease occurred even in the case of T. niphobles, a peripheral fresh water species that is often seen in brackish river mouths. Phylogenetic analysis using nucleotide sequences of the mitochondrial 16S ribosomal RNA gene of each species indicated that the six Takifugu species are very closely related with each other. CONCLUSION: T. obscurus is capable of adapting to both FW and SW. Its genomic sequence shares a very high homology with those of the other Takifugu species such that the existing Takifugu genomic information resources can be utilized. These properties make "mefugu", which has drawn little attention from animal physiologists until this study, a useful model animal for studying the molecular mechanism of maintaining body fluid homeostasis.