Statistical information characterization of conserved non-coding elements in vertebrates.

Recently, a set of highly conserved non-coding elements (CNEs) has been derived from a comparison between the genomes of the puffer fish, Takifugu or Fugu rubripes, and man. In order to facilitate the identification of these conserved elements in silico, we characterize them by a number of statistical features. We found a pronounced information pattern around CNE borders; although the CNEs themselves are AT rich and have high entropy (complexity), they are flanked by GC-rich regions of low entropy (complexity). We also identified the most abundant motifs within and around of CNEs, and identified those that group around their borders. Like in human promoter regions, the TBP, NF-Y and some other binding motifs are clustered around CNE boundaries, which may suggest a possible transcription regulatory function of CNEs.

Duplicated receptors for VIP and PACAP (VPAC1R and PAC1R) in a teleost fish, Fugu rubripes.

Two principal groups of receptors orthologous with human PAC1R and VPAC1R and were identified and characterised at the genomic level in the teleost fish Fugu rubripes. An additional group orthologous with VPAC2R was also identified and partially characterised. In Fugu, gene duplication of each of the PAC1Rs, VPAC1Rs and VPAC2Rs appears to have occurred. The topology of the tree surrounding the Fugu duplications and other isolated piscine sequences indicates that the duplication events for these six genes clearly preceded the speciation event leading to the Cypriniformes and Tetraodontiformes and is probably teleost-specific. Overall, the combined pattern of gene expression for each pair of duplicated genes mirrored the expression in other vertebrates. However, within each pair of duplicates further specialisation had occurred, with each demonstrating differential tissue distribution profiles suggesting they that may be responsible for the divergent action of the ligands, vasoactive intestinal peptide (VIP) and pituitary adenylate cyclase-activating peptide (PACAP). The Fugu VPAC1R gene regions showed conserved synteny with human chromosome 3p21.3 and also C. elegans chromosome X, indicating that the putative ancestral human chromosome 3 region may be equivalent to chromosome X in Caenorhabditis elegans.

Calcitonin: characterisation and expression in a teleost fish, Fugu rubripes.

The present report describes the structure and expression of the calcitonin gene in Fugu rubripes. It is composed of 4 exons and 3 introns. Splicing of exons 1, 2 and 3 generates the calcitonin pre-proprotein, while splicing of exons 1, 2 and 4 generates calcitonin gene-related protein (CGRP). Exons 1 and 2 encoding the signal sequence and the N-terminal peptide are common in both the gene products and this gene organisation has been conserved in human, rat, chicken and salmon. The gene environment around calcitonin in Fugu has been poorly conserved when compared with human, apart from a small gene cluster. The calcitonin gene in Fugu has a widespread tissue distribution but it is most highly expressed in the brain. The abundance of gene expression in the ultimobranchial gland and the pituitary indicates that these are important sites of production and that the peptide is probably secreted into the circulation and/or acts as a paracrine or autocrine controlling factor. Whilst the function of calcitonin in fish is still largely unknown, the distribution described here suggests that one of the potential functions may be as a neuropeptide.

Characterisation of two topoisomerase 1 genes in the pufferfish (Fugu rubripes).

Eukaryotic DNA topoisomerase I manipulates the higher order structures of DNA. Only one functional topoisomerase 1 (top1) gene has previously been identified in any individual eukaryotic species. Here we report the identification and characterisation of two top1 genes in the pufferfish, Fugu rubripes. This shows that the copy number of top1, like that of other topoisomerases, may vary between eukaryotes. Both Fugu genes have 21 exons; a gene structure similar to that of human TOP1. Despite this conservation of structure, and some non-coding elements, both genes are less than a tenth of the size of the human gene. Sequence and phylogenetic analyses have shown that this duplication is ancient and also affects other species in the fish lineage.

Characterization of the MHC class I region of the Japanese pufferfish (Fugu rubripes).

A BAC map of the Japanese pufferfish (Fugu) MHC class I region was constructed using a mixture of sequence scanning and sequence-tagged site mapping methodologies. The Fugu MHC class Ia genes are linked to genes which are found within the human classical MHC class II and extended class II regions, a situation which has been found in the MHC of all teleosts mapped so far. The 300-kb contig comprises 24 MHC-related genes and is bounded by six non-MHC genes, which are thought to represent an evolutionary breakpoint within the region. Comparative analysis with both human and zebrafish MHC maps indicates two blocks of genes (KNSL2, ZNF297, DAXX, TAPBP, FLOTILLIN; and PSMB8, PSMB10, PSMB9, ABCB3, FABGL, BRD2, COL11A2, RXRB) which have remained linked over 400 million years and may represent an ancestral arrangement of the vertebrate MHC. Zebrafish and Fugu diverged between 100-200 million years ago and differences exist between these two fish species. The position and number of MHC class Ia genes is not conserved between species, there is an inversion of a block of nine genes centering on the PSMB cluster, and additional genes are present in zebrafish coding for a transport-associated protein and a beta proteasome subunit. The extent of these differences has implications for the extrapolation of fish model organism data to commercial aquaculture species. The data presented here represent the most extensive analysis of a fish MHC class Ia region described so far and clearly delimit the extent of this region in Fugu and, potentially, all teleosts.

SAND, a new protein family: from nucleic acid to protein structure and function prediction.

As a result of genome, EST and cDNA sequencing projects, there are huge numbers of predicted and/or partially characterised protein sequences compared with a relatively small number of proteins with experimentally determined function and structure. Thus, there is a considerable attention focused on the accurate prediction of gene function and structure from sequence by using bioinformatics. In the course of our analysis of genomic sequence from Fugu rubripes, we identified a novel gene, SAND, with significant sequence identity to hypothetical proteins predicted in Saccharomyces cerevisiae, Schizosaccharomyces pombe, Caenorhabditis elegans, a Drosophila melanogaster gene, and mouse and human cDNAs. Here we identify a further SAND homologue in human and Arabidopsis thaliana by use of standard computational tools. We describe the genomic organisation of SAND in these evolutionarily divergent species and identify sequence homologues from EST database searches confirming the expression of SAND in over 20 different eukaryotes. We confirm the expression of two different SAND paralogues in mammals and determine expression of one SAND in other vertebrates and eukaryotes. Furthermore, we predict structural properties of SAND, and characterise conserved sequence motifs in this protein family.

Use of the Japanese pufferfish (Fugu rubripes) in comparative genomics.

With the draft sequence of the human genome available and an increasing number of organisms being sequenced, attention is becoming focused on sequence interpretation and functional analysis. Comparative genomics will play an important role in evaluating these data. At the molecular level, roles for uncharacterized proteins can be hypothesized by identifying conserved protein domains and putative noncoding regulatory elements can be defined from direct sequence comparisons of evolutionarily distant organisms. At a higher level, questions, such as the importance of gene order positioning, conservation of linkage, and genome evolution, can begin to be answered by collecting map data from different organisms. This minireview, centering on Fugu regions sharing synteny with human chromosomes 11p, 20q, and 6p21.3, details some of the ways in which the Japanese pufferfish can contribute to the study of comparative genomics and evaluation of sequence data from the genome programs.

Comparative mapping of the human 9q34 region in Fugu rubripes.

Twenty-seven genes have been cloned and mapped in Fugu which have orthologues within the human chromosome 9q34 region. The genes are arranged into five cosmid and BAC contigs which physically map to two different Fugu chromosomes. Considering the gene content of these contigs, it is more probable that a translocation event took place early in the Fugu lineage to split the ancestral 9q34 region onto two chromosomes rather than the alternative hypothesis of a large-scale duplication of the region into two chromosomes with subsequent rapid and dramatic gene loss. There are considerable differences in gene order between the two species, which would appear to be the result of a series of complex chromosome inversions; thus suggesting that there have been no positional constraints on this particular gene set.

Identification and characterization of a beta proteasome subunit cluster in the Japanese pufferfish (Fugu rubripes).

The low molecular mass polypeptide (LMP2, LMP7, and MECL-1) genes code for beta-type subunits of the proteasome, a multimeric complex that degrades proteins into peptides as part of the MHC class I-mediated Ag-presenting pathway. These gene products are up-regulated in response to infection by IFN-gamma and replace the corresponding constitutively expressed subunits (X, Y, and Z) during the immune response. In humans, the LMP2 and LMP7 genes both reside within the class II region of the MHC (6p21.3), while MECL-1 is located at 16q22.1. In the present study, we have identified all three IFN-gamma-regulated beta-type proteasome subunits in Fugu, which are present as a cluster within the Fugu MHC class I region. We show that in this species, LMP7, LMP2, and MECL-1 are linked. Also within this cluster is an LMP2-like subunit (which seems specific to all teleosts tested to date) and a closely linked LMP7 pseudogene, indicating that within Fugu and potentially other teleosts, there has been an additional regional duplication involving these genes.

Characterization of the Fugu rubripes NLK and FN5 genes flanking the NF1 (Neurofibromatosis type 1) gene in the 5' direction and mapping of the human counterparts.

To complete the analysis of the Neurofibromatosis type 1 (NF1) gene region in Fugu rubripes, we characterized the upstream flanking region of the NF1 gene and identified the FN5 (flanking the Fugu NF1 gene in 5' direction) gene and the NLK (Nemo-like kinase) gene as its flanking genes. The FN5 gene spans 3807bp and encompasses four exons, three of which belong to the expanded 5' UTR. Only 11% of the FN5 transcript is protein-coding. The function of the FN5 protein spanning 59 amino acids is unknown. We also characterized the human and the mouse FN5 transcripts and found 85% and 83% similarity of deduced amino acid sequences compared with Fugu. Two copies of the human FN5 gene were identified, one on chromosome 17q21.3-q22 several megabases distal to the NF1 gene at 17q11.2. The second copy of the FN5 gene was mapped to 11q13.3-q23.3. In Fugu, the FN5 gene is flanked by the NLK gene, which spans 4513bp from the translation start to the stop codon and encompasses 11 exons. Comparing the deduced amino acid sequences, 82% overall similarity was observed between Fugu and mouse or human NLK and 67% similarity between the Fugu NLK and the highly related LIT-1 kinase of Caenorhabditis elegans, which has been shown, like the vertebrate counterpart, to be involved in the Wnt signalling pathway. We mapped the human NLK gene to 17q11.2 between markers D17S935 and D17S120, more than 1Mb proximal to the NF1 gene. The characterization of the 5' flanking region presented here, together with that of the 3' region, demonstrates the profound differences between Fugu and human considering the gene content within the region flanking the NF1 gene.

Genomic structure and expression of parathyroid hormone-related protein gene (PTHrP) in a teleost, Fugu rubripes.

In this study we describe the isolation and characterisation of the parathyroid hormone-related protein (PTHrP) gene from the teleost Fugu rubripes. The gene has a relatively simple structure, compared with tetrapod PTHrP genes, composed of three exons and two introns, encompassing 2.25kb of genomic DNA. The gene encodes a protein of 163 amino acids, with a putative signal peptide of 37 amino acids and a mature peptide of 126 amino acids. The overall homology with known tetrapod PTHrP proteins is low (36%), with a novel sequence inserted between positions 38 and 65, the absence of the conserved pentapeptide (TRSAW) and shortened C-terminal domain. The N-terminus shows greater conservation (62%), suggesting that it may have a hypercalcaemic function similar to that of tetrapod PTHrP. In situ localisation and RT-PCR have demonstrated the presence of PTHrP in a wide range of tissues with varying levels of expression. Sequence scanning of overlapping cosmids has identified three additional genes, TMPO, LDHB and KCNA1, which map to human chromosome 12, with the latter two mapping to 12p12-11.2. PTHrP in human also maps to this chromosome 12 sub-region, thus demonstrating conservation of synteny between human and Fugu.

Identification and analysis of two snail genes in the pufferfish (Fugu rubripes) and mapping of human SNA to 20q.

All members of the snail gene family are zinc-finger transcription factors expressed early in embryonic development and are involved in the formation of tissues such as mesoderm and presumptive neural crest. Here, we report the identification and structural organisation of two snail genes in the compact genome of the pufferfish Fugu rubripes, and examine the phylogenetic relationships between these and other members of the snail gene family. Both genes have a three exon, two intron structure similar to that previously reported for human SLUG. While human SLUG has been mapped to 8q (Cohen, M.E., Yin, M., Paznekas, W.A., Schertzer, M., Wood, S., Jabs, E.W., 1998. Human SLUG organisation, expression and chromosome map location on 8q. Genomics 51, 468-471), the human sna gene SNA, was previously unmapped. We have used sequence similarity to the Fugu genes to identify a human SNA EST and mapped this by radiation hybrid and physical mapping to the distal end of human 20q. This is likely to be the mapping location of the human sna gene (SNA).

Comparative analysis of the PCOLCE region in Fugu rubripes using a new automated annotation tool.

The Japanese pufferfish Fugu rubripes with a genome of about 400 Mb is becoming increasingly recognized as a vertebrate model organism for comparative gene analysis (see Elgar 1996 for review). We have isolated and sequenced two Fugu cosmids spanning a genomic region of 66 kb containing the Fugu homolog to the human PCOLCE-I (Glöckner et al. 1998). We then examined if RUMMAGE-DP, a newly developed analysis tool for gene discovery which was designed for human and mouse genomic DNA, can be used for automatic annotation of Fugu genomic sequence. The exon prediction programs contained in RUMMAGE-DP performed better overall for the human sequence than for the Fugu contig. The GENSCAN program was the only exon prediction programme that performed equally well for both organisms. We show that RUMMAGE-DP is very useful in automatic analysis of Fugu sequences. Comparative analysis of the genomic structure of the PCOLCE-I genes in Fugu and human reveals that the exon/intron structure throughout the protein coding region is almost identical. We defined an additional domain based on the high degree of similarity of 26 aa between mammals and Fugu. The PCOLCE-I protein in both organisms contains two highly conserved CUB domains. Exons 6 and 7 are the only coding exons that differ in length between the two species. We assume that these exons do not code for any catalytic domain of the protein. Analysis of the remaining five Fugu genes within the 66 kb interval revealed no conserved synteny with the corresponding human 7q22 region.

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.

Isolation and characterisation of the retinoic acid receptor-alpha gene in the Japanese pufferfish, F. rubripes.

Nuclear hormone receptors (NRs) are ligand-inducible transcription factors that mediate critical functions in many species. The majority of novel NRs have hitherto been cloned from cDNA libraries by virtue of their homology to previously identified receptors. In this study, we validate a genomic DNA-based approach to isolating NRs by cloning the retinoic acid receptor-alpha (RARalpha) gene from the genome of the Japanese pufferfish, Fugu rubripes. The fRARalpha gene is more compact than its human and murine counterparts and demonstrates a highly conserved genomic organisation and amino acid sequence, generating two isoforms (fRARalpha1 and fRARalpha2) with divergent aminoterminal domains. In addition, a conserved regulatory element containing a retinoic acid response element was identified upstream of the fRARalpha2-specific exon, implying that retinoid induction of this isoform is evolutionarily conserved and critical to its function in vivo. We propose two uses for the Fugu genome in the study of NRs: the isolation of novel NRs that exhibit restricted spatio-temporal expression from genomic DNA and the identification of evolutionarily conserved promoter or intragenic regulatory DNA elements.

Characterization of three genes, AKAP84, BAW and WSB1, located 3' to the neurofibromatosis type 1 locus in Fugu rubripes.

Sequence analysis of cosmid clones was instrumental to identify three genes in the region flanking the Fugu rubripes NF1 gene in the 3' direction: the AKAP84 gene (A-kinase anchor protein 84), the WSB1 gene (WD-40-repeat protein with a SOCS box) and the BAW gene of yet unknown function located between the AKAP84 and the WSB1 genes. The human homologues of these genes are not located in the immediate vicinity of the NF1 gene at 17q11.2. Although synteny of the NF1, AKAP84, BAW and WSB1 genes is conserved between Fugu and human, the gene order is not conserved, and more than a simple inversion would have been necessary to explain the difference in gene order. The mammalian homologue of the Fugu BAW gene or protein has not yet been characterized. As deduced from the respective cDNAs, the Fugu AKAP84, WSB1 and BAW proteins vary concerning the overall degree of similarity to their mammalian counterparts. Whereas the overall similarity of AKAP84 between Fugu and mouse is low, three regions of known functional importance show considerable conservation. These are the N-terminal anchoring domain mediating the insertion of AKAP84 in the outer mitochondrial membrane, the binding site of the regulatory subunit (RI or RII) of protein kinase A, and the C-terminal domain present in the alternatively spliced isoform AKAP121 with an hnRNP K homology domain involved in RNA binding. A higher overall similarity of deduced protein sequences between Fugu and mouse was observed comparing the BAW gene products (74.1%) and the WSB1 proteins (77.2%).

Characterization of the transcription factor MTF-1 from the Japanese pufferfish (Fugu rubripes) reveals evolutionary conservation of heavy metal stress response.

The pufferfish Fugu rubripes was recently introduced as a new model organism for genomic studies, since it contains a full set of vertebrate genes but only 13% as much DNA as a mammal. Fugu genes tend to be smaller and densely spaced due to shortening of introns and intergenic spacers. We isolated the Fugu gene for the metal-responsive transcription factor MTF-1 (MTF1), a mediator of heavy metal regulation and oxidative stress response previously characterized in mammals. In addition, most of the cDNA sequence was also determined. The 780 amino acid MTF-1 protein of Fugu is very similar to that of mouse and human, with 90% amino acid identity in the DNA binding zinc finger domain and 57% overall identity. Expression of the pufferfish cDNA in mammalian cells shows that Fugu MTF-1 has the same DNA binding specificity as its mammalian counterpart and also induces transcription in response to zinc and cadmium. The protein-coding part of the Fugu MTF-1 gene spans 6.4 kb and consists of 11 exons. Upstream region and first exon constitute a CpG island. The distance between stop codon and polyadenylation motifs is >2 kb, suggesting a very long 3' untranslated mRNA region, followed by another CpG island which may represent the promoter of the next gene downstream. Part of the MTF-1 genomic structure was also determined in the mouse, and some striking similarities were found: for example, the upstream adjacent gene in both species is INPP5P, encoding a phosphatase. The mouse MTF-1 promoter is also embedded in a CpG island, which however shares no sequence similarity to the one of Fugu. The Fugu CpG island is shorter than the one of the mouse and has no elevated [G+C] content; these and other data indicate that CpG islands of fish may represent a primordial stage of CpG island evolution.

Three receptor genes for plasminogen related growth factors in the genome of the puffer fish Fugu rubripes.

Plasminogen related growth factors (PRGFs) and their receptors play major roles in embryogenesis, tissue regeneration and neoplasia. In order to investigate the complexity and evolution of the PRGF receptor family we have cloned and sequenced three receptors for PRGFs in the teleost fish Fugu rubripes, a model vertebrate with a compact genome. One of the receptor genes isolated encodes the orthologue of mammalian MET, whilst the other two may represent Fugu rubripes orthologues of RON and SEA. This is the first time three PRGF receptors have been identified in a single species.