FoxK1splice variants show developmental stage-specific plasticity of expression with temperature in the tiger pufferfish

FoxK1 is a member of the highly conserved forkhead/winged helix (Fox)family of transcription factors and it is known to play a key role in mammalian muscle development and myogenic stem cell function. The tiger pufferfish (Takifugu rubripes) orthologue of mammalian FoxK1(TFoxK1) has seven exons and is located in a region of conserved synteny between pufferfish and mouse. TFoxK1 is expressed as three alternative transcripts: TFoxK1-α, TFoxK1-γ and TFoxK1-δ. TFoxK1-α is the orthologue of mouse FoxK1-α, coding for a putative protein of 558 residues that contains the forkhead and forkhead-associated domains typical of Fox proteins and shares 53% global identity with its mammalian homologue. TFoxK1-γ and TFoxK1-δ arise from intron retention events and these transcripts translate into the same 344-amino acid protein with a truncated forkhead domain. Neither are orthologues of mouse FoxK1-β. In adult fish, the TFoxK1 splice variants were differentially expressed between fast and slow myotomal muscle, as well as other tissues, and the FoxK1-α protein was expressed in myogenic progenitor cells of fast myotomal muscle. During embryonic development, TFoxK1 was transiently expressed in the developing somites, heart,brain and eye. The relative expression of TFoxK1-α and the other two alternative transcripts varied with the incubation temperature regime for equivalent embryonic stages and the differences were particularly marked at later developmental stages. The developmental expression pattern of TFoxK1 and its localisation to mononuclear myogenic progenitor cells in adult fast muscle indicate that it may play an essential role in myogenesis in T. rubripes.

Ligand-dependent transcriptional activities of four torafugu pufferfish Takifugu rubripes peroxisome proliferator-activated receptors

Structural and functional properties were investigated for four peroxisome proliferator-activated receptors (PPARs), PPARalpha1, PPARalpha2, PPARbeta, and PPARgamma, from torafugu pufferfish Takifugu rubripes and determined for their transcriptional activity by the reporter assay using reporter plasmids containing three copies of the acyl-CoA oxidase PPAR response element. Although torafugu PPARs showed a high similarity in the primary structure to other vertebrate counterparts, torafugu PPARalpha2 and gamma contained additional sequences of 21 and 28 amino acids, respectively, as in the case of other teleost fish species when compared with African clawed frog counterparts. The transcriptional activity of torafugu PPARalpha1 was enhanced 4.5- and 11.5-fold by Wy-14643 and 5,8,11,14-eicosatetraynoic acid (ETYA) each at 10 microM, respectively, whereas that of PPARalpha2, 4.5- and 7.3-fold at the same concentration of the respective ligands, respectively. The activities of torafugu PPARalpha1 and alpha2 were also enhanced 5.6- and 6.3-fold by ETYA at 1 microM, respectively, but not by Wy-14643 at this concentration. Furthermore, the activities of the two torafugu PPARalphas were enhanced 4.3- and 7.6-fold by arachidonic acid, 4.4- and 5.2-fold by docosahexaenoic acid, and 6.7- and 8.0-fold by eicosapentaenoic acid each at 50 microM, respectively. On the other hand, the activities of torafugu PPARbeta and gamma were not changed by Wy-14643, ETYA, rosiglitazone, nor PUFAs. These results suggest that the activities of torafugu PPARbeta and gamma require undefined ligands. Alternatively, the molecular mechanisms involved in their activation are different from those of other vertebrates.

Fish calcitonin receptor has novel features

Calcitonin (CT), a 32-amino acid peptide, was initially isolated from fish. Fish CT has higher affinity to mammalian CT receptor (CTR), and has activity on calcium homeostasis. Therefore, fish CT has been used as a drug for the treatment of human bone diseases. However, the physiological roles of CT in fish as well as the characteristics of the fish CTR have not been clarified. Here, we cloned and characterized CTR from mefugu (Takifugu obscurus). Full-length cDNA sequencing revealed that mfCTR (mf, mefugu) consists of N-terminal four tandem putative hormone-binding domains (HBDs). Database mining showed that the multiple HBD-containing CTR is a common feature for some other fishes. Detailed pharmacological studies revealed that mfCTR generated cAMP in response to (1) fish CT, (2) calcitonin gene-related peptide (CGRP) in combinations with receptor activity-modifying proteins (mfRAMPs) 1 and 4, and (3) amylin in combinations with mfRAMPs 1-5. Unlike mammalian CTR, mfCTR showed dual affinity sites. Corresponding EC(50) values of those are in close proximity of the in vivo concentration of CT in fish. Analyses of the deletion mutants of mfCTR demonstrated that only the nearmost HBD to the first transmembrane region is functional to the ligands. Although, fish CT has higher affinity to the human CTR, human CT did not bind to the mfCTR. This is the first report that demonstrates the structure and property of fish receptor for CT, CGRP, and amylin. Fish CTR is the first example that has multiple HBD-like sequences.

Ultraconserved non-coding sequence element controls a subset of spatiotemporal GLI3 expression

The zinc-finger transcription factor GLI3 acts during vertebrate development in a combinatorial, context-dependent fashion as a primary transducer of sonic hedgehog (SHH) signaling. In humans, mutations affecting this key regulator of development are associated with GLI3-morphopathies, a group of congenital malformations in which forebrain and limb development are preferentially affected. We show that a non-coding element from intron two of GLI3, ultraconserved in mammals and highly conserved in the pufferfish Fugu, is a transcriptional enhancer. In transient transfection assays, it activates reporter gene transcription in human cell cultures expressing endogenous GLI3 but not in GLI3 negative cells. The identified enhancer element is predicted to contain conserved binding sites for transcription factors crucial for developmental steps in which GLI3 is involved. The regulatory potential of this element is conserved and was used to direct tissue-specific expression of a green fluorescent protein reporter gene in zebrafish embryos and of a beta-galactosidase reporter in transgenic mouse embryos. Time, location, and quantity of reporter gene expression are congruent with part of the pattern previously reported for endogenous GLI3 transcription.

cDNA cloning of Runx family genes from the pufferfish (Fugu rubripes)

The Runx family genes are involved in hematopoiesis, osteogenesis and neuropoiesis, and mutations in these genes have been frequently associated with human hereditary diseases and cancers. Here we report the cDNA cloning of the full Runx gene family of the pufferfish (Fugu rubripes), which comprises frRunx1, frRunx2, frRunx3, frRunt and frCbfb. Fugu is evolutionarily distant from mammals, thus the annotation of the frRunx family genes greatly facilitates comparative genomics approaches. Protein sequence comparison revealed that the fugu genes show high conservation in the Runt domain and PY and VWRPY motifs. frRunx1 had an extra stretch of eight histidine residues, while frRunx2 lacked the poly-glutamine/-alanine stretch that is a hallmark of the mammalian Runx2 genes. Analysis of the promoter regions revealed high conservation of the binding sites for transcription factors, including Runx sites in the P1 promoters. Abundant CpG dinucleotides in the P2 promoter regions were also detected. The expression patterns of the frRunx family genes in various tissues showed high similarity to those of the mammalian Runx genes. The genomic structures of the fugu and mammalian Runx genes are largely conserved except for a split exon 2 in frRunx1 and an extra exon in the C-terminal region of frRunx3 that is missing in mammalian Runx3 genes. The similarities and differences between the Runx family genes of fugu and mammals will improve our understanding of the functions of these proteins.

Identification and characterization of suppressor of cytokine signaling 1 (SOCS-1) homologues in teleost fish

The suppressor of cytokine signaling 1 (SOCS-1), as the name implies, is a member of proteins that regulate cytokine signaling pathways by inhibiting the events of key tyrosine phosphorylation on cytokine receptors and signaling molecules such as Janus kinase family members. Although mammalian SOCS-1 homologues were characterized in several species, no similar research work has been reported in fish yet. In this paper, we initially cloned the SOCS-1 genes from Tetraodon nigroviridis and Danio rerio, and identified other two SOCS-1 genes from Fugu rubripes and Gasterosteus aculeatus. The results showed that the fish SOCS-1-encoding genes consist of two exons and a single intron, a typical characteristic of SOCS family in gene organization. Moreover, two alternatively spliced transcripts that encoded 220 and 196 amino acids were obtained in T. nigroviridis, proving the distinct existence of alternative splicing in fish SOCS-1 different from higher vertebrates. By reverse transcriptase polymerase chain reaction (PCR) and real-time quantitative PCR, gene expression studies indicated that both two alternatively spliced transcripts of Tetraodon SOCS-1 were expressed extensively in major tissues as we examined and their corresponding expression levels could be strikingly raised at 3 h postinjection with lipopolysaccharide, which strongly suggested that SOCS-1 proteins in fish might be involved in inflammatory responses. This is the first report of cloning and characterization of SOCS-1 complementary deoxyribonucleic acids and genes in fish.

Identification of novel genes related to tetrodotoxin intoxication in pufferfish

To investigate the genes related to the biosynthesis or accumulation of tetrodotoxin (TTX) in pufferfish, mRNA expression patterns in the liver from pufferfish, akamefugu Takifugu chrysops and kusafugu Takifugu niphobles, were compared by mRNA arbitrarily primed reverse transcription-polymerase chain reaction (RAP RT-PCR) with fish bearing different concentrations of TTX and its derivatives. RAP RT-PCR provided a 383 bp cDNA fragment and its transcripts were higher in toxic than non-toxic pufferfish liver. Its deduced amino acid sequence was similar to those of fibrinogen-like proteins reported for other vertebrates. Northern blot analysis and rapid amplification of cDNA ends (RACE) revealed that the cDNA fragment of 383 bp was composed of at least three fibrinogen-like protein (flp) genes, flp-1, flp-2 and flp-3. Relative mRNA levels of flp-1, flp-2 and flp-3 showed a linear correlation with toxicity of the liver for two pufferfish species.

Patterns of gene duplication and intron loss in the ENCODE regions suggest a confounding factor

The exon-intron structure of eukaryotic genes allows for phenomena such as alternative splicing, nonsense-mediated decay, and regulation through untranslated regions. However, the evolution of the exon structure of genes is not well elucidated because of limited and phylogenetically sparse data sets. In this study, we use the phylogenetically diverse sequencing of the ENCODE regions to study gene structure evolution in mammalian genomes. This first phylogenetically diverse study of gene structure changes offers insights into the mode and tempo of mammalian gene structure evolution. The genes undergoing structure changes appear to be moderately to highly expressed in germline cells and show levels of selection similar to those of other ENCODE genes. Patterns of gene duplication of the affected genes are more complex than expected. The number of sampled genomes is sufficiently dense to infer that certain gene duplications happened after intron loss. Thus, although gene duplication is highly correlated with intron loss, we conclude that structural changes in genes are not necessarily due to a loss of constraint following gene duplication as previously suggested.

Duplicated members of the Groucho/Tle gene family in fish

The highly conserved Groucho/Tle gene family has widespread functions during embryonic development and in adults. For mammalians, four full-length Tle paralogues are known, whereas the whole spectrum of this gene family in fish species has not been analysed yet. Most detailed data exist for medaka, where 3 Tle genes have been described, Tle1, Tle3, and Tle4. We now isolated 3 additional Tle genes from the medaka genome. Sequence analysis identifies these genes as Tle2a, Tle2b, and Tle3b. Database searches of genomic sequences revealed an identical set of Tle paralogues being present in distantly related fish species, indicating duplicated Tle2 and Tle3 genes for the complete teleost lineage. Like the previously analysed medaka Tle genes, the three new genes show a broad expression pattern during embryogenesis. Nevertheless, a detailed comparison of all six Tle genes reveals critical differences in certain aspects of their expression pattern. In particular, we concentrated on the activity of Tle genes during ear development and found Tle2a and Tle2b expressed in this sensory organ.

Analysis of the exon–intron structures of fish, amphibian, bird and mammalian hatching enzyme genes, with special reference to the intron loss evolution of hatching enzyme genes in Teleostei

Using gene cloning and in silico cloning, we analyzed the structures of hatching enzyme gene orthologs of vertebrates. Comparison led to a hypothesis that hatching enzyme genes of Japanese eel conserve an ancestral structure of the genes of fishes, amphibians, birds and mammals. However, the exon-intron structure of the genes was different from species to species in Teleostei: Japanese eel hatching enzyme genes were 9-exon-8-intron genes, and zebrafish genes were 5-exon-4-intron genes. In the present study, we further analyzed the gene structures of fishes belonging to Acanthopterygii. In the species of Teleostei we examined, diversification of hatching enzyme gene into two paralogous genes for HCE (high choriolytic enzyme) and LCE (low choriolytic enzyme) was found only in the acanthopterygian fishes such as medaka Oryzias latipes, Fundulus heteroclitus, Takifugu rubripes and Tetraodon nigroviridis. In addition, the HCE gene had no intron, while the LCE gene consisted of 8 exons and 7 introns. Phylogenetic analysis revealed that HCE and LCE genes were paralogous to each other, and diverged during the evolutionary lineage to Acanthopterygii. Analysis of gene synteny and cluster structure showed that the syntenic genes around the HCE and LCE genes were highly conserved between medaka and Teraodon, but such synteny was not found around the zebrafish hatching enzyme genes. We hypothesize that the zebrafish hatching enzyme genes were translocated from chromosome to chromosome, and lost some of their introns during evolution.

A Teleost Polymeric Ig Receptor Exhibiting Two Ig-Like Domains Transports Tetrameric IgM into the Skin

The skin mucus IgM is an important molecule in the mucosal immune system of teleost skin. However, the transport mechanism associated with this molecule has yet to be clarified. In this study, we isolated a gene encoding a polymeric Ig receptor (pIgR) from a species of teleost fish, Takifugu rubripes (fugu). This gene is known to be an Ig transporter in the intestine of mammals. Our studies further demonstrated that fugu pIgR was expressed in the skin and that a fragment of pIgR bound to tetrameric IgM in the skin mucus. These results indicate that the skin pIgR transports tetrameric IgM into the skin mucus. The fugu pIgR exhibits a unique structure containing only two Ig-like domains corresponding to domain 1 and domain 4/5 of mammalian pIgR. This structure was sufficient for successful binding to tetrameric IgM. Teleost skin thus adopts the same Ig transport system as mammalian intestine via a unique pIgR.

Differential regulation of multiple alternatively spliced transcripts of MyoD

Splice variants of the basic helix-loop-helix myoblast determination factor (myoD) have not been previously found in vertebrates. Here we report the identification and characterization of three alternative transcripts of a myoD paralogue from the tiger pufferfish (Takifugu rubripes). The T. rubripes myoD1 gene (TmyoD1) has 3 exons and 2 introns and it is present on scaffold 104, in a region of conserved synteny with zebrafish. The isoform TMyoD1-alpha is a putative protein of 281 residues that contains the basic, helix-loop-helix and helix III domains and shares 61%, 56%, 51%, 49% and 56% overall identity with zebrafish, Xenopus, mouse, human and chicken MyoD1, respectively. TMyoD1-beta arises from an alternative 3' splice site and differs from TMyoD1-alpha by a 26-residue insertion adjacent to helix III, which is one of the functional domains required for chromatin remodelling. The third alternative transcript, TmyoD1-gamma, retains intron I and has two premature termination codons far from the 3'-most exon-exon junction. TmyoD1-gamma is therefore likely to be degraded by nonsense-mediated decay, an important widespread post-transcriptional mechanism that regulates transcript levels. Analysis of gene expression by qPCR revealed that TmyoD1-alpha was the most abundant transcript in fast and slow myotomal muscle. TmyoD1-alpha expression was 2-fold higher in fast muscle of juvenile fish that were actively producing new myotubes compared to adult stages that had stopped recruiting fast muscle fibres. A similar expression pattern was observed for TmyoD1-alpha in slow muscle but the differences were not significant. Transcript levels of TmyoD1-gamma only varied significantly in fast muscle and were 5-fold higher in adult compared to juvenile stages. Significant differences in expression of TmyoD1 splice variants were also observed during embryonic development. The differential expression of three alternative transcripts of myoD1 in developing and adult myotomal muscle of T. rubripes supports the hypothesis that diversity generated by alternative splicing may be of functional significance in muscle development in this species.

Sequencing and comparative analysis of fugu protocadherin clusters reveal diversity of protocadherin genes among teleosts

Background

The synaptic cell adhesion molecules, protocadherins, are a vertebrate innovation that accompanied the emergence of the neural tube and the elaborate central nervous system. In mammals, the protocadherins are encoded by three closely-linked clusters (α, β and γ) of tandem genes and are hypothesized to provide a molecular code for specifying the remarkably-diverse neural connections in the central nervous system. Like mammals, the coelacanth, a lobe-finned fish, contains a single protocadherin locus, also arranged into α, β and γ clusters. Zebrafish, however, possesses two protocadherin loci that contain more than twice the number of genes as the coelacanth, but arranged only into α and γ clusters. To gain further insight into the evolutionary history of protocadherin clusters, we have sequenced and analyzed protocadherin clusters from the compact genome of the pufferfish, Fugu rubripes.

Results

Fugu contains two unlinked protocadherin loci, Pcdh1 and Pcdh2, that collectively consist of at least 77 genes. The fugu Pcdh1 locus has been subject to extensive degeneration, resulting in the complete loss of Pcdh1γ cluster. The fugu Pcdh genes have undergone lineage-specific regional gene conversion processes that have resulted in a remarkable regional sequence homogenization among paralogs in the same subcluster. Phylogenetic analyses show that most protocadherin genes are orthologous between fugu and zebrafish either individually or as paralog groups. Based on the inferred phylogenetic relationships of fugu and zebrafish genes, we have reconstructed the evolutionary history of protocadherin clusters in the teleost fish lineage.

Conclusion

Our results demonstrate the exceptional evolutionary dynamism of protocadherin genes in vertebrates in general, and in teleost fishes in particular. Besides the 'fish-specific' whole genome duplication, the evolution of protocadherin genes in teleost fishes is influenced by lineage-specific gene losses, tandem gene duplications and regional sequence homogenization. The dynamic protocadherin clusters might have led to the diversification of neural circuitry among teleosts, and contributed to the behavioral and physiological diversity of teleosts.

The sex-determining locus in the tiger pufferfish, Takifugu rubripes

The tiger pufferfish (fugu), Takifugu rubripes, is a model fish that has had its genome entirely sequenced. By performing genomewide linkage analyses, we show that the sex of fugu is determined by a single chromosomal region on linkage group 19 in an XX-XY system.

Characterization of two paralogous muscleblind-like genes from the tiger pufferfish (Takifugu rubripes)

Muscleblind-like (Mbnl) proteins are required for terminal muscle differentiation in mammals. In this study we have identified two mbnl paralogues from the tiger pufferfish, tmbnl2a and tmbnl3, which are the first examples of non-mammalian mbnl genes. Tmbnl2a and tmbnl3 were found in regions of conserved synteny and had a high degree of global conservation with their mammalian homologues. Phylogenetic analysis showed that the T. rubripes genome contains one mbnl3 gene and two copies of mbnl1 and mbnl2. Moreover, the mbnl1 and mbnl3 paralogues are derived from duplication of a common ancestral gene. The average rates of synonymous substitutions between T. rubripes, mouse and human mbnl2 and mbnl3 genes were much higher than the corresponding rates of non-synonymous mutations, suggesting that Mbnl2 and Mbnl3 are subjected to strong purifying selection. Quantitation of tmbnl2a and tmbnl3 transcripts by real-time PCR revealed that these two paralogues are differentially expressed in fast and slow myotomal muscle, heart, liver, skin, brain and testes. Tmbnl2a was expressed at similar levels in all tissues examined, as was the mouse orthologue. Tmbnl3 was expressed at higher levels than tmbnl2a, with a ubiquitous tissue distribution. Expression of tmbnl3 remained high in adult pufferfish muscle whereas the mouse orthologue was down-regulated in adults, perhaps reflecting the indeterminate and determinate growth patterns of these taxa, respectively.

Identification, cDNA cloning, and mRNA localization of a zebrafish ortholog of leukemia inhibitory factor

Leukemia inhibitory factor (LIF) maintains embryonic stem cells in an undifferentiated state. To establish stable cultures of zebrafish embryonic stem cells, recombinant zebrafish LIF (zLIF) is needed because the LIF sequence varies greatly between species. In the current study, we identified the zebrafish (Danio rerio) and pufferfish (Tetraodon nigroviridis) orthologs of lif from genomic databases, and we isolated a cDNA encoding zLIF. Synteny analysis and comparison of sequences identified zebrafish and Tetraodon orthologs of human LIF. The cDNA for zLIF encoded a predicted 215-amino acid protein with a putative 32-amino acid signal peptide, two disulfide bonds, and two N-linked glycosylation sites. We found that transcription of zlif starts at the hatching period during embryogenesis and is present in the brain, visceral organs, bone, and skin.

Ancient noncoding elements conserved in the human genome

Cartilaginous fishes represent the living group of jawed vertebrates that diverged from the common ancestor of human and teleost fish lineages about 530 million years ago. We generated approximately 1.4x genome sequence coverage for a cartilaginous fish, the elephant shark (Callorhinchus milii), and compared this genome with the human genome to identify conserved noncoding elements (CNEs). The elephant shark sequence revealed twice as many CNEs as were identified by whole-genome comparisons between teleost fishes and human. The ancient vertebrate-specific CNEs in the elephant shark and human genomes are likely to play key regulatory roles in vertebrate gene expression.

The genome size evolution of medaka (Oryzias latipes) and fugu (Takifugu rubripes)

Evolution of the genome size in eukaryotes is often affected by changes in the noncoding sequences, for which insertions and deletions (indels) of small nucleotide sequences and amplification of repetitive elements are considered responsible. In this study, we compared the genomic DNA sequences of two kinds of fish, medaka (Oryzias latipes) and fugu (Takifugu rubripes), which show two-fold difference in the genome size (800 Mb vs. 400 Mb). We selected a contiguous DNA sequence of 790 kb from the medaka chromosome LG22 (linkage group 22), and made a precise comparison with the sequence (387 kb) of the corresponding region of Takifugu. The sequence of 178 kb in total was aligned common between two fishes, and the remaining sequences (612 kb for medaka and 209 kb for fugu) were found abundant in various repetitive elements including many types of unclassified low copy repeats, all of which accounted for more than a half (54%) of the genome size difference. Furthermore, we identified a significant difference in the length ratio of the unaligned sequences that locate between the aligned sequences (USBAS), particularly after eliminating known repetitive elements. These USBAS with no repetitive elements (USBAS-nr) located within the intron and intergenic region. These results strongly indicated that amplification of repetitive elements and compilation of indels are major driving forces to facilitate changes in the genome size.

Cracking the genome's second code: enhancer detection by combined phylogenetic footprinting and transgenic fish and frog embryos

Genes involved in vertebrate development are unusually enriched for highly conserved non-coding sequence elements. These regions are readily detected in silico, by genome-wide sequence comparisons between different vertebrates, from mammals to fish (phylogenetic footprinting). It follows that sequence conservation must be the result of positive selection for an essential physiological role. An obvious possibility is that these conserved sequences possess regulatory or structural functions important for gene expression and, thus, an in vivo assay becomes necessary. We have developed a rapid testing system using zebrafish and Xenopus laevis embryos that allows us to assign transcriptional regulatory functions to conserved non-coding sequence elements. The sequences are cloned into a vector containing a minimal promoter and the GFP reporter, and are assayed for their putative cis-regulatory activity in zebrafish or Xenopus transgenic experiments. Vectors used include plasmid DNA and the Tol2 transposon system in fish and X. laevis. We have followed this logic to detect and analyze conserved elements in an intergenic region present in the Iroquois (Irx) gene clusters of zebrafish, Xenopus tropicalis, Fugu rubripes and mouse. We have assayed approximately 50 of these conserved elements and shown that the majority behave as modular positive regulatory elements (enhancers) that contribute to specific temporal and spatial domains that are part of the endogenous gene expression pattern. Moreover, comparison of the activity of cognate Irx enhancers from different organisms demonstrates that conservation of sequence is accompanied by in vivo functional conservation across species. Finally, for some of the most conserved elements, we have been able to identify a critical core sequence, essential for correct enhancer function.

Characterization and tissue distribution of multiple agouti-family genes in pufferfish, Takifugu rubripes

Four types of agouti-family genes (AGRP1, AGRP2, ASIP1 and ASIP2) were obtained from torafugu, Takifugu rubripes. Their characterization and structure were analyzed to elucidate the relationship among the torafugu agouti-family genes. Both AGRP1 and AGRP2 showed genomic synteny with the human AGRP gene. Phylogenetic tree analysis showed that AGRP1 formed a cluster with human AGRP. We inferred that torafugu AGRP1 and AGRP2 are orthologs of human AGRP and that they are paralogous genes derived from genome duplication occurred in the teleost phylogeny. Torafugu ASIP1 showed genomic synteny with the human ASIP, but ASIP2 did not. The ASIP1 expression level was about five times higher in the white ventral skin than in the black dorsal skin. Therefore, we concluded that torafugu ASIP1 is an ortholog of human ASIP, nevertheless, we are unable to determine if torafugu ASIP2 is a paralog of ASIP1 or not.

In vivo enhancer analysis of human conserved non-coding sequences

Identifying the sequences that direct the spatial and temporal expression of genes and defining their function in vivo remains a significant challenge in the annotation of vertebrate genomes. One major obstacle is the lack of experimentally validated training sets. In this study, we made use of extreme evolutionary sequence conservation as a filter to identify putative gene regulatory elements, and characterized the in vivo enhancer activity of a large group of non-coding elements in the human genome that are conserved in human-pufferfish, Takifugu (Fugu) rubripes, or ultraconserved in human-mouse-rat. We tested 167 of these extremely conserved sequences in a transgenic mouse enhancer assay. Here we report that 45% of these sequences functioned reproducibly as tissue-specific enhancers of gene expression at embryonic day 11.5. While directing expression in a broad range of anatomical structures in the embryo, the majority of the 75 enhancers directed expression to various regions of the developing nervous system. We identified sequence signatures enriched in a subset of these elements that targeted forebrain expression, and used these features to rank all approximately 3,100 non-coding elements in the human genome that are conserved between human and Fugu. The testing of the top predictions in transgenic mice resulted in a threefold enrichment for sequences with forebrain enhancer activity. These data dramatically expand the catalogue of human gene enhancers that have been characterized in vivo, and illustrate the utility of such training sets for a variety of biological applications, including decoding the regulatory vocabulary of the human genome.

Flounder and fugu have a single lefty gene that covers the functions of lefty1 and lefty2 of zebrafish during L-R patterning

The lefty gene encodes a member of the TGF-beta superfamily that regulates L-R axis formation during embryogenesis via antagonistic activity against Nodal, another TGF-beta superfamily member. Both mouse and zebrafish have two lefty genes, lefty1 and lefty2. Interestingly, the expression domains of mouse and zebrafish lefty are different from one another. At present, the orthology and functional diversity of the mouse and zebrafish lefty genes are not clear. Here, we report that flounder and two fugu species, Takifugu and Tetraodon, have a single lefty gene in their genomes. In addition, we provide evidence that the mouse lefty genes were duplicated on a single chromosome but the zebrafish lefty genes arose from a whole-genome duplication that occurred early in the divergence of ray-finned fishes. These independent origins likely explain the difference in the expression domains of the mouse and zebrafish lefty gene pairs. Furthermore, we found that the duplication corresponding to the zebrafish lefty2 gene was lost from the fugu genome, suggesting that loss of lefty2 in the fugu/flounder lineage occurred after its divergence from the zebrafish lineage. During L-R patterning, the single lefty gene of flounder covers two expression domains, the left side of the dorsal diencephalon and the left LPM, which are regulated separately by lefty1 and lefty2 in zebrafish. We infer that the lefty genes of the ray-finned fishes and mammals underwent independent gene duplication events that resulted in independent regulation of lefty expression.

Comparative study of family 2 GPCRs in Fugu rubripes

In this study, members of family 2 GPCRs, one of the largest families of receptors in vertebrates, were isolated and characterized in the genome of the Japanese pufferfish, Fugu rubripes, and compared with the orthologous genes in other vertebrates. Phylogenetic analysis carried out with all vertebrate family 2 GPCR members indicated that CALR/CGRPR and CRF are the most divergent receptor group within this family and that the remaining members appear to originate from a common ancestral gene precursor.

[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.

Genetic basis of tetrodotoxin resistance in pufferfishes

Tetrodotoxin (TTX) is a highly potent neurotoxin that selectively binds to the outer vestibule of voltage-gated sodium channels. Pufferfishes accumulate extremely high concentrations of TTX without any adverse effect. A nonaromatic amino acid (Asn) residue present in domain I of the pufferfish, Takifugu pardalis, Na v1.4 channel has been implicated in the TTX resistance of pufferfishes . However, the effect of this residue on TTX sensitivity has not been investigated, and it is not known if this residue is conserved in all pufferfishes. We have investigated the genetic basis of TTX resistance in pufferfishes by comparing the sodium channels from two pufferfishes (Takifugu rubripes [fugu] and Tetraodon nigroviridis) and the TTX-sensitive zebrafish. Although all three fishes contain duplicate copies of Na v1.4 channels (Na v1.4a and Na v1.4b), several substitutions were found in the TTX binding outer vestibule of the two pufferfish channels. Electrophysiological studies showed that the nonaromatic residue (Asn in fugu and Cys in Tetraodon) in domain I of Na v1.4a channels confers TTX resistance. The Glu-to-Asp mutation in domain II of Tetraodon channel Na v1.4b is similar to that in the saxitoxin- and TTX-resistant Na+ channels of softshell clams . Besides helping to deter predators, TTX resistance enables pufferfishes to selectively feed on TTX-bearing organisms.