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. Received: 13 October 1998 / Accepted: 25 July 1999     

Generation and Characterization of a genetic zebrafish model of SMA carrying the human SMN2 gene

Background

Animal models of human diseases are essential as they allow analysis of the disease process at the cellular level and can advance therapeutics by serving as a tool for drug screening and target validation. Here we report the development of a complete genetic model of spinal muscular atrophy (SMA) in the vertebrate zebrafish to complement existing zebrafish, mouse, and invertebrate models and show its utility for testing compounds that alter SMN2 splicing.

Results

The human motoneuron disease SMA is caused by low levels, as opposed to a complete absence, of the survival motor neuron protein (SMN). To generate a true model of SMA in zebrafish, we have generated a transgenic zebrafish expressing the human SMN2 gene (hSMN2), which produces only a low amount of full-length SMN, and crossed this onto the smn ^-/- background. We show that human SMN2 is spliced in zebrafish as it is in humans and makes low levels of SMN protein. Moreover, we show that an antisense oligonucleotide that enhances correct hSMN2 splicing increases full-length hSMN RNA in this model. When we placed this transgene on the smn mutant background it rescued the neuromuscular presynaptic SV2 defect that occurs in smn mutants and increased their survival.

Conclusions

We have generated a transgenic fish carrying the human hSMN2 gene. This gene is spliced in fish as it is in humans and mice suggesting a conserved splicing mechanism in these vertebrates. Moreover, antisense targeting of an intronic splicing silencer site increased the amount of full length SMN generated from this transgene. Having this transgene on the smn mutant fish rescued the presynaptic defect and increased survival. This model of zebrafish SMA has all of the components of human SMA and can thus be used to understand motoneuron dysfunction in SMA, can be used as an vivo test for drugs or antisense approaches that increase full-length SMN, and can be developed for drug screening.     

The neural cell adhesion molecule L1: genomic organisation and differential splicing is conserved between man and the pufferfish Fugu

The human gene for the neural cell adhesion molecule L1 is located on Xq28 between the ALD and MeCP2 loci. Mutations in the L1 gene are associated with four related neurological disorders, X-linked hydrocephalus, spastic paraplegia (SPG1), MASA syndrome, and X-linked corpus callosum agenesis. The clinical relevance of L1 has led us to sequence the L1 gene in human and to investigate its conservation in the vertebrate model genome of the pufferfish, Fugu rubripes (Fugu), a species with a compact genome of around 40 Mb. For this purpose we have sequenced a human and a Fugu cosmid clone containing the corresponding L1 genes. For comparison, we have also amplified and sequenced the complete Fugu L1 cDNA. We find that the genomic structure of L1 is conserved. The human and Fugu L1 gene both have 28 exons of nearly identical size. Differential splicing of exons 2 and 27 is conserved over 430 million years, the evolutionary time span between the teleost Fugu and the human L1 gene. In contrast to previously published Fugu genes, many introns are larger in the Fugu L1 gene, making it slightly larger in size despite the compact nature of the Fugu genome. Homology at the amino acid and the nucleotide level with 40% and 51%, respectively, is lower than that of any previously reported Fugu gene. At the level of protein structure, both human and Fugu L1 molecules are composed of six immunoglobulin (Ig)-like domains and five fibronectin (Fn) type III domains, followed by a transmembrane domain and a short cytoplasmic domain. Only the transmembrane and the cytoplasmic domains are significantly conserved in Fugu, supporting their proposed function in intracellular signalling and interaction with cytoskeletal elements in the process of neurite outgrowth and fascicle formation. Our results show that the cytoplasmic domain can be further subdivided into a conserved and a variable region, which may correspond to different functions. Most pathological missense mutations in human L1 affect conserved residues. Fifteen out of 22 reported missense mutations alter amino acids that are identical in both species.     

Characterization of the pufferfish Otx2 cis-regulators reveals evolutionarily conserved genetic mechanisms for vertebrate head specification

The Otx2 gene, containing a highly conserved paired-type homeobox, plays a pivotal role in the development of the rostral head throughout vertebrates. Precise regulation of the temporal and spatial expression of Otx2 is likely to be crucial for proper head specification. However, regulatory mechanisms of Otx2 expression remain largely unknown. In this study, the Otx2 genome of the puffer fish Fugu rubripes, which has been proposed as a model vertebrate owing to its highly compact genome, was cloned. Consistently, Fugu Otx2 possesses introns threefold smaller in size than those of the mouse Otx2 gene. Otx2 mRNA was transcribed after MBT, and expressed in the rostral head region throughout the segmentation and pharyngula periods of wild-type Fugu embryos. To elucidate regulatory mechanisms of Otx2 expression, the expression of Otx2-lacZ reporter genes nearly covering the Fugu Otx2 locus, from -30.5 to +38.5 kb, was analyzed, by generating transgenic mice. Subsequently, seven independent cis-regulators were identified over an expanse of 60 kb; these regulators are involved in the mediation of spatiotemporally distinct subdomains of Otx2 expression. Additionally, these expression domains appear to coincide with local signaling centers and developing sense organs. Interestingly, most domains do not overlap with one another, which implies that cis-regulators for redundant expression may be abolished exclusively in the pufferfish so as to reduce its genome size. Moreover, these cis-regions were also able to direct expression in zebrafish embryos equivalent to that observed in transgenic mice. Further comparative sequence analysis of mouse and pufferfish intergenic regions revealed eight highly conserved elements within these cis-regulators. Therefore, we propose that, in vertebrate evolution, the Otx2 promoter acquires multiple, spatiotemporally specific cis-regulators in order to precisely control highly coordinated processes in head development.     

Cloning and sequencing of complement component C9 and its linkage to DOC-2 in the pufferfish Fugu rubripes

The Japanese pufferfish Fugu rubripes has a 400 Mb genome with high gene density and minimal non-coding complexity, and is therefore an ideal vertebrate model for sequence comparison. The identification of regions of conserved synteny between Fugu and humans would greatly accelerate the mapping and ordering of genes. Fugu C9 was cloned and sequenced as a first step in an attempt to characterize the region in Fugu homologous to human chromosome 5p13. The 11 exons of the Fugu C9 gene share 33% identity with human C9 and span 2.9 kb of genomic DNA. By comparison, human C9 spans 90 kb, representing a 30-fold difference in size. We have also determined by cosmid sequence scanning that DOC-2, a tumour suppresser gene which also maps to human 5p13, lies 6–7 kb from C9 in a head-to-head or 5′ to 5′ orientation. These results demonstrate that the Fugu C9/DOC-2 locus is a region of conserved synteny. Sequence scanning of overlapping cosmids has identified two other genes, GAS-1 and FBP, both of which map to human chromosome 9q22, and lie adjacent to the Fugu C9/DOC-2 locus, indicating the boundary between two syntenic regions.     

Genomic structure and sequence of the pufferfish (Fugu rubripes) growth hormone-encoding gene: a comparative analysis of teleost growth hormone genes

A nested polymerase chain reaction (PCR) technique for amplifying a fragment of the gene (GH) encoding teleost growth hormone has been developed. Using this technique, a fragment of the pufferfish, Fugu rubripes and Arothron maculatus; dwarf gourami, Colisa lalia; guppy, Poecilia reticulata; and goldfish, Carassius auratus GH genes were cloned. The Fugu rubripes (Fugu) gene fragment was used to isolate the GH gene from a Fugu genomic library. The complete nucleotide sequence of a 8.5-kb SacI genomic fragment containing the Fugu GH gene has been determined. The GH gene spans 2.5 kb from the first codon to polyadenylation signal, and contains six exons and five introns similar to the GH genes of salmonids, tilapia, barramundi, flounder and yellowtail. The GH introns contain microsatellite and satellite sequences. The microsatellites found in the fifth intron of the GH gene are also present in the corresponding introns of tilapia, barramundi and flounder GH genes. Southern analysis revealed that the GH gene is a single-copy gene in the Fugu. The promoter region of the Fugu GH gene contains conserved sequences that are likely to be involved in the pituitary-specific expression of the gene. A phylogenetic tree of nucleotide (nt) sequences of all known teleost GH genes has been inferred using the distance matrix method. The topology of this tree reflects the major phylogenetic groupings of teleosts. The intron patterns and repetitive sequences of GH genes can serve as useful natural markers for the classification and phylogenetic studies of teleosts.     

In vivo characterization of a vertebrate ultraconserved enhancer

Genomic sequence comparisons among human, mouse, and pufferfish (Takifugu rubripes (Fugu)) have revealed a set of extremely conserved noncoding sequences. While this high degree of sequence conservation suggests severe evolutionary constraint and predicts a lack of tolerance to change to retain in vivo functionality, such elements have been minimally explored experimentally. In this study, we describe the in-depth characterization of an ancient conserved enhancer, Dc2, located near the dachshund gene, which displays a human-Fugu identity of 84% over 424 basepairs (bp). In addition to this large overall conservation, we find that Dc2 is characterized by the presence of a large block of sequence (144 bp) that is completely identical among human, mouse, chicken, zebrafish, and Fugu. Through the testing of reporter vector constructs in transgenic mice, we observed that the 424-bp Dc2-conserved element is necessary and sufficient for brain tissue enhancer activity. In vivo analyses also revealed that the 144-bp 100% conserved sequence is necessary, but not sufficient, to replicate Dc2 enhancer function. However, the introduction of two separate 16-bp insertions into the highly conserved enhancer core did not cause any detectable modification of its in vivo activity. Our observations indicate that the 144-bp 100% conserved element is tolerant of change at least at the resolution of this transgenic mouse assay and suggest that purifying selection on the Dc2 sequence might not be as strong as we predicted or that some unknown property also constrains this highly conserved enhancer sequence.     

Somatic mutagenesis with a Sleeping Beauty transposon system leads to solid tumor formation in zebrafish

Large-scale sequencing of human cancer genomes and mouse transposon-induced tumors has identified a vast number of genes mutated in different cancers. One of the outstanding challenges in this field is to determine which genes, when mutated, contribute to cellular transformation and tumor progression. To identify new and conserved genes that drive tumorigenesis we have developed a novel cancer model in a distantly related vertebrate species, the zebrafish, Danio rerio. The Sleeping Beauty (SB) T2/Onc transposon system was adapted for somatic mutagenesis in zebrafish. The carp ß-actin promoter was cloned into T2/Onc to create T2/OncZ. Two transgenic zebrafish lines that contain large concatemers of T2/OncZ were isolated by injection of linear DNA into the zebrafish embryo. The T2/OncZ transposons were mobilized throughout the zebrafish genome from the transgene array by injecting SB11 transposase RNA at the 1-cell stage. Alternatively, the T2/OncZ zebrafish were crossed to a transgenic line that constitutively expresses SB11 transposase. T2/OncZ transposon integration sites were cloned by ligation-mediated PCR and sequenced on a Genome Analyzer II. Between 700–6800 unique integration events in individual fish were mapped to the zebrafish genome. The data show that introduction of transposase by transgene expression or RNA injection results in an even distribution of transposon re-integration events across the zebrafish genome. SB11 mRNA injection resulted in neoplasms in 10% of adult fish at ∼10 months of age. T2/OncZ-induced zebrafish tumors contain many mutated genes in common with human and mouse cancer genes. These analyses validate our mutagenesis approach and provide additional support for the involvement of these genes in human cancers. The zebrafish T2/OncZ cancer model will be useful for identifying novel and conserved genetic drivers of human cancers.     

Comparative genomics using fugu: a tool for the identification of conserved vertebrate cis-regulatory elements

With the imminent completion of the whole genome sequence of humans, increasing attention is being focused on the annotation of cis-regulatory elements in the human genome. Comparative genomics approaches based on evolutionary conservation have proved useful in the detection of conserved cis-regulatory elements. The pufferfish, Fugu rubripes, is an attractive vertebrate model for comparative genomics, by virtue of its compact genome and maximal phylogenetic distance from mammals. Fugu has lost a large proportion of nonessential DNA, and retained single orthologs for many duplicate genes that arose in the fish lineage. Non-coding sequences conserved between fugu and mammals have been shown to be functional cis-regulatory elements. Thus, fugu is a model fish genome of choice for discovering evolutionarily conserved regulatory elements in the human genome. Such evolutionarily conserved elements are likely to be shared by all vertebrates, and related to regulatory interactions fundamental to all vertebrates. The functions of these conserved vertebrate elements can be rapidly assayed in mammalian cell lines or in transgenic systems such as zebrafish/medaka and Xenopus, followed by validation of crucial elements in transgenic rodents.     

The first completed genome sequence from a teleost fish (Fugu rubripes) adds significant diversity to the nuclear receptor superfamily

Defining complete sets of gene family members from diverse species provides the foundation for comparative studies. Using a bioinformatic approach, we have defined the entire nuclear receptor complement within the first available complete sequence of a non-human vertebrate (the teleost fish Fugu rubripes). In contrast to the human set (48 total nuclear receptors), we found 68 nuclear receptors in the Fugu genome. All 68 Fugu receptors had a clear human homolog, thus defining no new nuclear receptor subgroups. A reciprocal analysis showed that each human receptor had one or more Fugu orthologs, excepting CAR (NR1I3) and LXRβ (NR1H2). These 68 receptors add striking diversity to the known nuclear receptor superfamily and provide important comparators to human nuclear receptors. We have compared several pharmacologically relevant human nuclear receptors (FXR, LXRα/β, CAR, PXR, VDR and PPARα/γ/δ) to their Fugu orthologs. This comparison included expression analysis across five Fugu tissue types. All of the Fugu receptors that were analyzed by PCR in this study were expressed, indicating that the majority of the additional Fugu receptors are likely to be functional.     

Analysis of the conservation of synteny between Fugu and human chromosome 12

Background

The pufferfish Fugu rubripes (Fugu) with its compact genome is increasingly recognized as an important vertebrate model for comparative genomic studies. In particular, large regions of conserved synteny between human and Fugu genomes indicate its utility to identify disease-causing genes. The human chromosome 12p12 is frequently deleted in various hematological malignancies and solid tumors, but the actual tumor suppressor gene remains unidentified.

Results

We investigated approximately 200 kb of the genomic region surrounding the ETV6 locus in Fugu (fETV6) in order to find conserved functional features, such as genes or regulatory regions, that could give insight into the nature of the genes targeted by deletions in human cancer cells. Seven genes were identified near the fETV6 locus. We found that the synteny with human chromosome 12 was conserved, but extensive genomic rearrangements occurred between the Fugu and human ETV6 loci.

Conclusion

This comparative analysis led to the identification of previously uncharacterized genes in the human genome and some potentially important regulatory sequences as well. This is a good indication that the analysis of the compact Fugu genome will be valuable to identify functional features that have been conserved throughout the evolution of vertebrates.

     

Tetraodon genome analysis provides further evidence for whole-genome duplication in the ray-finned fish lineage

A whole-genome duplication in the ray-finned fish lineage has been supported by the analyses of the genome sequence of the Japanese pufferfish, Fugu rubripes. Recently, genome sequence of a second teleost fish, the freshwater pufferfish, Tetraodon nigroviridis, was completed. Comparisons of long-range synteny between the Tetraodon and human genomes provided additional evidence for the whole-genome duplication in the ray-finned fish lineage. In the present study, we conducted phylogenetic analysis of the Tetraodon and human proteins to identify ray-finned fish lineage-specific (‘fish-specific’) duplicate genes in the Tetraodon genome. Our analyses provide evidence for 1087 well defined fish-specific duplicate genes in Tetraodon. We also analyzed the Fugu proteome that was predicted in the recent Fugu genome assembly, and identified 346 duplicate genes in addition to the 425 duplicates previously identified. We estimated the ages of duplicate genes using the molecular clock. The ages of duplicate genes in the two pufferfishes independently support a large-scale gene duplication around 380–400 Myr ago. In addition, a burst of recent gene duplications was evident in the Tetraodon lineage. These findings provide further evidence for a whole-genome duplication early in the evolution of ray-finned fishes, and suggest that independent gene duplications have occurred recently in the Tetraodon lineage.     

Genomic Structure and Nucleotide Sequence of the p55 Gene of the Puffer Fish Fugu rubripes

The p55 gene, which codes for a 55-kDa erythrocyte membrane protein, has been cloned and sequenced from the genome of the Japanese puffer fish Fugu rubripes (Fugu). This organism has the smallest recorded vertebrate genome and therefore provides an efficient way to sequence genes at the genomic level. The gene encoding p55 covers 5.5 kb from the beginning to the end of the coding sequence, four to six times smaller than the estimated size of the human gene, and is encoded by 12 exons. The structure of this gene has not been previously elucidated, but from this and other data we would predict a similar or identical structure in mammals. The predicted amino acid sequence of this gene in Fugu, coding for a polypeptide of 467 amino acids, is very similar to that of the human gene with the exception of the first two exons, which differ considerably. The predicted Fugu protein has a molecular weight (52.6 kDa compared with 52.3 kDa) and an isoelectric point very similar to those of human p55. In human, the p55 gene lies in the gene-dense Xq28 region, just 30 kb 3′ to the Factor VIII gene, and is estimated to cover 20-30 kb. Its 5′ end is associated with a CpG island, although there is no evidence that this is the case in Fugu. The small size of genes in Fugu and the high coding homology that they share with their mammalian equivalents, both in structure and sequence, make this compact vertebrate genome an ideal model for genomic studies.     

Something fishy in the rat brain: molecular genetics of the hypothalamo-neurohypophysial system

The brain peptides vasopressin and oxytocin play crucial roles in the regulation of salt and water balance. The genes encoding these neurohormones are regulated by cell-specific and physiological cues, but the molecular mechanisms remain obscure. New strategies, involving the introduction of rat transgenes into rats, are being used to address these issues,⁽¹⁾ but the complexity of the rat genome has hampered progress. By contrast, the pufferfish, Fugu rubripes, has a “junk-free” genome.⁽²⁾ The oxytocin homologue from Fugu, isotocin, has been introduced into rats⁽³⁾ and is expressed in oxytocin neurons, where it is upregulated by physiological perturbations that upregulate the oxytocin gene. The Fugu and rat lineages separated 400 million years ago, yet the mechanisms that regulate the isotocin and oxytocin genes have been conserved. Fugu genome analysis and transgenesis in the physiologically tractable rat host are a powerful combination that will enable the identification of fundamental components of the neural systems that control homeostasis. BioEssays 20 :741–749, 1998. © 1998 John Wiley & Sons, Inc.     

Compact intergenic regions of the pufferfish genome facilitate isolation of gene promoters: characterization of Fugu 3'-phosphoadenosine 5'-phosphosulfate synthase 2 (fPapss2) gene promoter function in transgenic Xenopus

The highly compact nature of the pufferfish (Fugu rubripes) genome renders it a useful tool not only for annotating coding regions within vertebrate genomes, but also for the identification of sequences important to gene regulation. Indeed, owing to this compaction it will be feasible in many instances to initiate analyses using entire intergenic regions when mapping gene promoters; a strategy that is very rarely feasible with the expanded genomes of other species. Stemming from our interest in studying promoters expressed in chondrocytes, we selected for study the intergenic region upstream of Fugu 3'-phosphoadenosine 5'-phosphosulfate synthase 2, fPapss2, a gene required for the normal development of cartilage extracellular matrix. Functional characterization of the entire fPapss2 5' intergenic region was carried out by monitoring expression of the enhanced green fluorescent protein (EGFP) gene reporter in the developing cartilage of transgenic Xenopus laevis. By evaluating a series of 5' intergenic region deletions we defined a minimal fPapss2 sequence of approximately 300 bp that was essential for EGFP expression in tadpole cartilage. This functional analysis of an entire Fugu intergenic region, combined with the efficiency of Xenopus transgenesis, serves as a model for the rapid characterization of evolutionarily-conserved regulatory regions of other pufferfish genes.     

Generation and characterization of a zebrafish muscle specific inducible Cre line

Zebrafish transgenic lines provide valuable insights into gene functions, cell lineages and cell behaviors during development. Spatiotemporal control over transgene expression is a critical need in many experimental approaches, with applications in loss- and gain-of-function expression, ectopic expression and lineage tracing experiments. The Cre/loxP recombination system is a powerful tool to provide this control and the demand for validated Cre and loxP zebrafish transgenics is high. One of the major challenges to widespread application of Cre/loxP technology in zebrafish is comparatively small numbers of established tissue-specific Cre or CreERT2 lines. We used Tol2-mediated transgenesis to generate Tg(CrymCherry;-1.9mylz2:CreERT2) which provides an inducible CreERT2 source driven by muscle-specific mylz2 promoter. The transgenic specifically labels the trunk and tail skeletal muscles. We assessed the temporal responsiveness of the transgenic by screening with a validated loxP reporter transgenic ubi:Switch. Further, we evaluated the recombination efficiency in the transgenic with varying concentrations of 4-OHT, for different induction time periods and at different stages of embryogenesis and observed that higher recombination efficiency is achieved when embryos are induced with 10 μM 4-OHT from 10-somites or 24 hpf till 48 or 72 hpf. The transgenic is an addition to currently available zebrafish transgenesis toolbox and a significant tool to advance muscle biology studies in zebrafish.     

Fugu: a compact vertebrate reference genome

At 400 Mb, the Japanese pufferfish, Fugu rubripes, has the smallest vertebrate genome but has a similar gene repertoire to other vertebrates. Its genes are densely packed with short intergenic and intronic sequences devoid of repetitive elements. It likely has a mutational bias towards DNA elimination and is probably close to a 'minimal' vertebrate genome. As such it is a useful reference genome for gene discovery and gene validation in other vertebrates. Its usefulness in the discovery of conserved regulatory elements has already been demonstrated. The Fugu genome sequence is a good complement to genetic studies in other vertebrates.