Genomic structure and sequence of the Fugu rubripes glucose-6-phosphate dehydrogenase gene (G6PD)

A Journey from Blood Cells to Genes and Back

I was attracted to hematology because by combining clinical findings with the use of a microscope and simple laboratory tests, one could often make a diagnosis. I was attracted to genetics when I learned about inherited blood disorders, at a time when we had only hints that somatic mutations were also important. It seemed clear that if we understood not only what genetic changes caused what diseases but also the mechanisms through which those genetic changes contribute to cause disease, we could improve management. Thus, I investigated many aspects of the glucose-6-phosphate dehydrogenase system, including cloning of the gene, and in the study of paroxysmal nocturnal hemoglobinuria (PNH), I found that it is a clonal disorder; subsequently, we were able to explain how a nonmalignant clone can expand, and I was involved in the first trial of PNH treatment by complement inhibition. I was fortunate to do clinical and research hematology in five countries; in all of them, I learned from mentors, from colleagues, and from patients.

Associations between G6PD, OATP1B1 and BLVRA variants and susceptibility to neonatal hyperbilirubinaemia in a Chinese Han population

AIM

Hyperbilirubinaemia is a common disorder in newborns. The aim of this study was to investigate the associations between G6PD 1388 G>A, SLCO1B1 rs4149056 and BLVRA rs699512 variants and the risk of neonatal hyperbilirubinaemia in a Chinese neonate population.

METHODS

A total of 447 Chinese neonates with hyperbilirubinaemia were selected as the study group and 544 healthy subjects were recruited as the control group matched by baseline sex, age, feeding pattern and delivery mode. About 2 mL of peripheral venous blood was taken from all subjects. The single nucleotide polymorphisms (SNPs) of G6PD 1388 G>A, SLCO1B1 rs4149056 and BLVRA rs699512 loci were examined by the polymerase chain reaction and Sanger sequencing technique in the peripheral blood of all subjects.

RESULTS

For the G6PD 1388 G>A SNP, individuals carrying the A-allele were associated with a significantly increased risk of neonatal hyperbilirubinaemia (adjusted odds ratio (OR) = 1.49, P < 0.001, 95% confidence interval (CI): 1.31-1.67). This risk increased significantly in the CC genotype carriers at the rs4149056 locus of the SLCO1B1 gene (OR = 2.17, 95% CI: 1.87-2.33), whereas it decreased significantly in individuals carrying the G-allele at the rs699512 locus of the BLVRA gene (adjusted OR = 0.84, P = 0.01, 95% CI: 0.75-0.95). The G6PD 1388 G>A, SLCO1B1 rs4149056 and BLVRA rs699512 SNPs had a significant impact on serum total bilirubin levels. Moreover, individuals carrying the A-allele of G6PD 1388 G>A and BLVRA rs699512 had a significantly increased risk of developing neonatal hyperbilirubinaemia (OR = 5.01, P < 0.001, 95% CI: 3.42-7.85).

CONCLUSION

Genetic variants of bilirubin metabolism genes, including G6PD 1388 G>A, SLCO1B1 rs4149056 and BLVRA rs699512, are associated with the risk of neonatal hyperbilirubinaemia, and are potential markers for predicting the disorder.

Glucose-6-phosphate dehydrogenase in blunt snout bream Megalobrama amblycephala: molecular characterization, tissue distribution, and the responsiveness to dietary carbohydrate levels

This study aimed to characterize the full-length cDNA of glucose-6-phosphate dehydrogenase (G6PD) from Megalobrama amblycephala with its responses to dietary carbohydrate levels characterized. The cDNA obtained covered 2768 bp with an open reading frame of 1572 bp. Sequence alignment and phylogenetic analysis revealed a high degree of conservation (77-97%) among most fish and other higher vertebrates. The highest transcription of G6PD was observed in kidney followed by liver, whereas relatively low abundance was detected in eye. Then, the transcriptions and activities of G6PD as well as lipid contents were determined in the liver, muscle, and the adipose tissue of fish fed two dietary carbohydrate levels (30 and 42%) for 12 weeks. Hepatic transcriptions of fatty acid synthetase (FAS), acetyl-CoA carboxylase α (ACCα), sterol regulatory element-binding protein-1 (SREBP1), and peroxisome proliferator-activated receptor γ (PPARγ) were also measured to corroborate the lipogenesis derived from carbohydrates. The G6PD expressions and activities in both liver and the adipose tissue as well as the lipid contents in whole-body, liver, and the adipose tissue all increased significantly after high-carbohydrate feeding. Hepatic transcriptions of FAS, ACCα, SREBP1, and PPARγ were also up-regulated remarkably by the intake of a high-carbohydrate diet. These results indicated that the G6PD of M. amblycephala shared a high similarity with that of other vertebrates. Its expressions and activities in tissues were both highly inducible by high-carbohydrate feeding, as also held true for the transcriptions of other enzymes and/or transcription factors involved in lipogenesis, evidencing an enhanced lipogenesis by high dietary carbohydrate levels.

Multiple transcripts encode glucose 6-phosphate dehydrogenase in the southern cattle tick, Rhipicephalus (Boophilus) microplus

Glucose 6-phosphate dehydrogenase (G6PDH) is an enzyme that plays a critical role in the production of NADPH. Here we describe the identification of four transcripts (G6PDH-A, -B, -C, and -D) that putatively encode the enzyme in the southern cattle tick, Rhipicephalus (Boophilus) microplus. The genomic DNA that is spliced to produce G6PDH-A and -B is 8,600-9,000 bases in length and comprises 12 exons. Comparison of the R. microplus G6PDH gene structure with those available from insects and mammals revealed that the tick gene is most like that of humans. Detection of the four transcripts was evaluated by quantitative RT-PCR using template from larvae, unfed adult females and males, salivary gland tissues from 2- to 3-day-fed adult females and males, and salivary gland tissue of 4- to 5-day-fed adult females. The G6PDH-A and -C transcripts were present in all templates, and both displayed induced expression in salivary gland tissue of fed, adult females but not matched males. The G6PDH-D transcript was detected only in unfed adults and in larvae, a stage in which it was most abundant relative to the other three transcripts. The G6PDH-B transcript, while detectable in all templates, was of low copy number suggesting it is a rare transcript. Induced expression of G6PDH-A and G6PDH-C in fed females may play a role in the tolerance of oxidative stress that is induced upon feeding, and the transcript abundance in fed females may be a function of bloodmeal volume and the time adult females spend on the host relative to adult males.

Use of two transcription starts in the G6PD gene of the bark beetle Ips typographus

The enzyme glucose-6-phosphate dehydrogenase (G6PD) of the bark beetle Ips typographus is derived from a gene that includes eight exons and spans over 7100 nucleotides (nt). By means of two transcription starts, the gene generates two mRNA isoforms that are present in similar amounts in the larvae, pupae and adults. The A isoform includes exon IA of 115 nt, which is followed by intron 1a extending to position 3457 of the gene. The B mRNA isoform begins with exon IB (100 nt) that occupies positions 3291-3390 within the 1a intron. Exons II to VII are included in both mRNA isoforms. The gene contains 31.6% (36.5% in the translated region) of the GC nucleotides. Two transcription starts and the exon/intron organization distinguish bark beetle G6PD from the homologous genes known in other insects. Two enzyme variants were detected in the protein extracts of individual bark beetles but their relationship to the A and B mRNA isoforms is uncertain.

Sublethal ammonia and urea concentrations inhibit rainbow trout (Oncorhynchus mykiss) erythrocyte glucose-6-phosphate dehydrogenase

In vitro and in vivo effects of sublethal ammonia and urea concentrations were assayed on glucose-6-phosphate dehydrogenase (G6PD) of rainbow trout (Oncorhynchus mykiss) erythrocyte. G6PD was purified from erythrocytes with a specific activity of 16.7 EU (mmol NADP+/min)/mg protein and approximately 1600-fold in a yield of approximately 60% by ammonium sulphate precipitation and 2',5'-ADP Sepharose 4B affinity chromatography. The purity of the enzyme was confirmed using SDS polyacrylamide gel electrophoresis. Experiments with ammonia (2.2-5.5 microM) and urea (20-50 microM) showed the inhibitory effects on the enzyme, in vitro. Inhibition effects were determined in vitro by Lineweaver-Burk and regression graphs. The dissociation constant of the enzyme inhibitor complex (Ki) and 50% inhibitory values were 2.26+/-1.21 and 2.86+/-3.51 microM for ammonia and 18.69+/-6.75 and 23.77+/-4.58 microM for urea, respectively. In vivo studies in rainbow trout erythrocytes showed significant (p < 0.01) inhibition of G6PD by ammonia and urea. However, ammonia inhibited more than urea since there were significant differences between the final values of erythrocyte G6PD activities.

Identification and comparative analysis of the RpL14 gene from Takifugu rubripes

The ribosomal protein RpL14 gene has been characterized in several species, including, human, rat and fruit fly. Haploinsufficiency for the gene causes the Minute phenotype in Drosophila, and it has been proposed as a regulator in the tumorigenic pathway in human. Several features concerning the gene structure have been studied, and some of these differ between human/rat and Drosophila. To address functional and evolutionary questions about these differences we have isolated and sequenced a cDNA and a genomic clone covering the RpL14 gene from the pufferfish Takifugu rubripes (Fugu). The Fugu RpL14 gene is approximately 2 Kb, with 5 introns, and encodes a protein of 137 amino acids. The protein contains a KOW-motif and a nuclear localization signal, which are conserved among a wide range of RPL14 proteins. On the other hand, a variable amino acid (alanine) repeat observed in human is missing in Takifugu rubripes, and the protein is shorter than its mammalian counterparts. Compared with human, the RpL14 gene in Fugu contains introns localized at identical positions in the gene, and most of them are shorter. A comparison of the RpL14 gene structure from a broad range of organisms indicates that both loss and gain of introns have occurred during the evolution of the gene.

Comparative analysis of vertebrate dystrophin loci indicate intron gigantism as a common feature

The human DMD gene is the largest known to date, spanning > 2000 kb on the X chromosome. The gene size is mainly accounted for by huge intronic regions. We sequenced 190 kb of Fugu rubripes (pufferfish) genomic DNA corresponding to the complete dystrophin gene (FrDMD) and provide the first report of gene structure and sequence comparison among dystrophin genomic sequences from different vertebrate organisms. Almost all intron positions and phases are conserved between FrDMD and its mammalian counterparts, and the predicted protein product of the Fugu gene displays 55% identity and 71% similarity to human dystrophin. In analogy to the human gene, FrDMD presents several-fold longer than average intronic regions. Analysis of intron sequences of the human and murine genes revealed that they are extremely conserved in size and that a similar fraction of total intron length is represented by repetitive elements; moreover, our data indicate that intron expansion through repeat accumulation in the two orthologs is the result of independent insertional events. The hypothesis that intron length might be functionally relevant to the DMD gene regulation is proposed and substantiated by the finding that dystrophin intron gigantism is common to the three vertebrate genes.

The mitochondrial genome of the pufferfish, Fugu rubripes, and ordinal teleostean relationships

The small nuclear genome of the pufferfish, Fugu rubripes (order Tetraodontiformes), makes this species highly interesting for genome research. In order to establish the phylogenetic position of the Tetraodontiformes relative to other teleostean orders that might also have a reduced nuclear genome size, we have sequenced the mitochondrial (mt) genome of the pufferfish. The gene order, nucleotide composition and evolutionary rate of the mt genome of the fugu correspond to those of other teleosts. This suggests that the evolution of this genome has not been affected by the processes that led to the dramatic reduction of the size of the nuclear genome of the fugu. The phylogenetic analyses, which were based on the concatenated amino acid sequences of twelve protein-coding mt genes, placed the fugu among the percomorphs. The affinities between the Tetraodontiformes and either the Perciformes or the Zeiformes were limited, however. The common notion of a separate euteleostean clade remained unsupported. The analyses did not support the traditional systematic understanding that the Clupeiformes constitute a basal teleostean lineage. In addition the findings strongly suggest that three teleostean orders, the Perciformes, Zeiformes and Scorpaeniformes, are paraphyletic.

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

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

Molecular cloning and characterization of the Fugu rubripes MEST/COPG2 imprinting cluster and chromosomal localization in Fugu and Tetraodon nigroviridis

We isolated Fugu genomic clones using the human MEST (Mesoderm-Specific Transcript) cDNA as probe. Sequence analysis revealed the presence of MEST and three additional genes which show homology to plant DNBP (DNA-Binding Protein), vertebrate COPG2 (Coat Protein Gamma 2), as well as to human and mouse UCN (Urocortin). Structures of Fugu and human MEST, COPG2 and UCN genes are very similar. Since MEST and COPG2 are neighboring genes on human chromosome 7q32, we can conclude that we identified their orthologs and that linkage of these genes is evolutionarily conserved in vertebrates. Unlike human MEST which underlies isoform-specific imprinting and is methylated in a parent-of-origin-specific fashion, the CpG island of the Fugu ortholog is completely methylated. The translation start of Fugu MEST is identical to the non-imprinted human isoform which is in good agreement with the assumption that genomic imprinting is restricted to mammals. Comparative mapping of these genes by fluorescence in-situ hybridization to metaphase chromosomes of Fugu rubripes and Tetraodon nigroviridis showed clear signals on one of the smallest acrocentric chromosomal pairs, which in Fugu, can be easily classified by its unique triangular shape.

Two classes of genes in plants

Two classes of genes were identified in three Gramineae (maize, rice, barley) and six dicots (Arabidopsis, soybean, pea, tobacco, tomato, potato). One class, the GC-rich class, contained genes with no, or few, short introns. In contrast, the GC-poor class contained genes with numerous, long introns. The similarity of the properties of each class, as present in the genomes of maize and Arabidopsis, is particularly remarkable in view of the fact that these plants exhibit large differences in genome size, average intron size, and DNA base composition. The functional relevance of the two classes of genes is stressed by (1) the conservation in homologous genes from maize and Arabidopsis not only of the number of introns and of their positions, but also of the relative size of concatenated introns; and (2) the existence of two similar classes of genes in vertebrates; interestingly, the differences in intron sizes and numbers in genes from the GC-poor and GC-rich classes are much more striking in plants than in vertebrates.

Genomic Structure and Comparative Analysis of Nine Fugu Genes: Conservation of Synteny with Human Chromosome Xp22.2–p22.1

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

[Sequence data for this article were deposited with the EMBL/GenBank data libraries under accession nos. AJ011381 and AF094327.]

Genomic organisation of rainbow trout, Oncorhynchus mykiss TGF-beta

The genomic organisation of Oncorhnchus mykiss TGF-beta has been determined through the generation of contiguous clones by PCR. The O. mykiss TGF-beta gene is approximately 3.4 Kb in length and consists of 7 coding exons with no introns in the 5'-UTR. Whilst this is the same number of exons found in TGF-beta genes of amphibians, birds and mammals, in the O. mykiss gene intron 2 of other vertebrates is absent and an additional intron is present at the 3' end of the molecule, splitting exon 7 of the other known TGF-beta genes into two exons (trout exons 6 and 7). Comparison of exon sizes in the coding region support the suggestion that the Xenopus TGF-beta5 and trout TGF-beta sequences are the forerunners of TGF-beta1. Conservation of exons coding for the mature TGF-beta peptide is relatively high (63-73% identity) but other exons show lower identities (37-58%). Comparison of the TGF-beta intron sequences reveals that in general the O. mykiss introns are considerably shorter than the avian homologs. The impact of the teleost TGF-beta gene organisation on theories of the gene evolution of this cytokine family are discussed.

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.

Structures, sequence characteristics, and synteny relationships of the transcription factor E4TF1, the splicing factor U2AF35 and the cystathionine beta synthetase genes from Fugu rubripes

A cosmid containing the beta-amyloid precursor protein (APP) from Fugu rubripes has been completely sequenced. In addition to APP, the cosmid contains the E4TF1-60 transcription factor, the U2AF35 pre-mRNA splicing factor, and the cystathionine beta synthetase (CBS) gene. The human homologues of all four genes map to human chromosome 21 but are not clustered; APP and E4TF1-60 map within 21q21, whereas U2AF35 and CBS map approximately 20Mb distal in 21q22. 3. The protein sequences of the Fugu genes vary in their overall level of similarity to their mammalian homologues, but several regions of functional importance are almost identical. As expected, the intron/exon structures of the homologous pairs of genes are highly conserved, but there are significant differences in the compaction ratios. The introns of APP and E4TF1-60 are 49- and 24-fold smaller in Fugu than in human, and the intergenic distance is compressed at least 100-fold. For U2AF35 and CBS, the introns are compressed only five- to eightfold. These size differences were compared with those for a number of previously reported Fugu genes; in general, levels of compaction of Fugu genes are consistent with the isochore locations of the human homologues.

Analysis of pufferfish homologues of the AT-rich human APP gene

Mutations in the beta-amyloid precursor protein (APP) gene are associated with some forms of Familial Alzheimer's Disease. The human APP gene is large, the 19 exons span approximately 300 kb, and AT-rich, at 40% GC. We have examined the genomic structure and cDNA sequence of the APP gene in the pufferfish Fugu rubripes and Tetraodon fluviatilis, respectively. In contrast to human, the Fugu APP gene spans less than 10 kb of DNA, with the introns compacted 48-fold on average. Two axons, alternatively processed in humans, are absent in both pufferfish. APP is the largest, most AT-rich gene examined in Fugu and is also the most highly compressed. The genomic sequences spanning the human and the Fugu APP genes were analysed with a set of exon and gene prediction programs. Results show that these are highly reliable for the Fugu gene with lower false positive and false negative rates than are seen in the analysis of the human gene. Comparative analysis of Fugu sequences homologous to very AT-rich regions in the human genome may, therefore, be advantageous in gene-finding efforts, both for their highly reduced sizes and their reliable gene predictions.

Isolation and genomic characterization of the TUPLE1/HIRA gene of the pufferfish Fugu rubripes

In an effort to obtain a small genomic construct for the generation of a HIRA transgenic mouse, we have isolated and sequenced the Fugu TUPLE1/HIRA gene. We have compared the gene organization and the proteins encoded in pufferfish and human and also searched for conserved DNA sequences that might be important in gene regulation. The pufferfish gene spans approx. 9 kb, which is approx. 11 times smaller than the human gene, owing to the reduced size of the introns. Like its human counterpart, it is organized into 25 exons. The majority of the splice sites are in identical positions to those found in the human gene, however, for three internal exons the positions of the splice sites are not directly comparable. The coding regions are almost identical in size and show a high degree of similarity, especially at the amino and carboxy termini. Comparisons of 5' and 3' sequences failed to detect similarities or sequences involved in regulation.

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 40Mb. 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.

Clonability and gene distribution on human chromosome 21: reflections of junk DNA content?

Data from transcriptional mapping of human chromosome 21 have been compiled from a number of sources. Regardless of the gene identification technique used, a consistent picture has developed: the centromere proximal half of 21q, which contains 50% of the DNA (20 Mb), harbors only 10% of the expressed sequences. Because of the variety of gene isolation techniques used, this result is unlikely to arise simply from methodological artefacts, biases in clonability or tissue specificity of expression. This region is known to be AT-rich and to contain APP, the largest gene (spanning 300 kb) currently analyzed on 21q. Interesting preliminary data from analysis of the Fugu rubripes homolog of APP has shown an unusually high, 50-fold, compaction of intron size, raising the intriguing possibility that >90% of the DNA in the human gene may be functionless. Thus, data from a variety of approaches suggest that a large part of 21q very likely has neither coding capacity nor associated regulatory function. By these criteria, it is a good candidate for a repository of junk DNA.

The comparative genomic structure and sequence of the surfeit gene homologs in the puffer fish Fugu rubripes and their association with CpG-rich islands

The puffer fish Fugu rubripes (Fugu) has a compact genome approximately one-seventh the size of man, mainly owing to small intron size and the presence of few dispersed repetitive DNA elements, which greatly facilitates the study of its genes at the genomic level. It has been shown previously that, whereas the Surfeit genes are tightly clustered at a single locus in mammals and birds, the genes are found at three separate loci in the Fugu genome. Here, Fugu gene homologs of all six Surfeit genes (Surf-1 to Surf-6) have been cloned and sequenced, and their gene structure has been compared with that of their mammalian and avian homologs. The predicted protein products of each gene are well conserved between vertebrate species, and in most cases their gene structures are identical to their mammalian and avian homologs except for the Fugu Surf-6 gene, which was found to lack an intron present in the mouse gene. In addition, we have identified conserved regulatory elements at the 5' and 3' ends of the Surf-3/rpL7a gene by comparison with the mammalian and chicken Surf-3/rpL7a gene homologs, including the presence of a polypyrimidine tract at the extreme 5' end of this ribosomal protein gene. The Fugu Surfeit gene homologs appear to be associated with CpG-rich islands, like the Surfeit genes in higher vertebrates, but these Fugu CpG islands are similar to the nonclassical islands characteristic of other fish species. Our observations support the use of the Fugu genome to study vertebrate gene structure, to predict the structure of mammalian genes, and to identify vertebrate regulatory elements. [The sequence data described in this paper have been submitted to the data library under accession nos. Y15170 (Surf-2, Surf-4), Y15171 (Surf-3, Surf-1, Surf-6), and Y15172 (Surf-5.)]

Aberrant processing of the Fugu HD (FrHD) mRNA in mouse cells and in transgenic mice

The puffer fish ( Fugu rubripes ) has a compact genome of 400 Mbp which is approximately 7.5-fold smaller than the human genome. It contains a similar number of genes but is deficient in intergenic, intronic and dispersed repetitive sequences. Fugu is becoming established as the model vertebrate genome for the identification and characterisation of novel human genes and conserved regulatory sequences. It has also been proposed that Fugu genes may provide natural mini-genes for the production of transgenic mice. We have used the Fugu homologue of the Huntington's disease (HD) gene to test this possibility. The human and Fugu HD genes cover 170 kb and 23 kb respectively and have previously been sequenced in their entirety. In Fugu tissue, the Fugu HD gene was found to be expressed as predicted from the gene sequence but three differentially spliced forms were also detected. Despite the absence of conserved promoter sequences, the Fugu promoter was found to be functional in mouse cells. We have generated mice transgenic for the Fugu HD gene and conducted a detailed expression analysis across the entire 10 kb transcript. This revealed the presence of many aberrant splice forms which would be incompatible with the production of the Fugu huntingtin protein. The Fugu HD gene is incorrectly processed in mouse cells both in vitro and in vivo which sheds doubt on the usefulness of Fugu genes for transgenesis.

Towards a reconciliation of the introns early or late views: triosephosphate isomerase genes from insects

The gene encoding the glycolytic enzyme, triosephosphate isomerase (TPI; EC 5.3.1.1), is a favourite model for molecular evolutionists who either subscribe to the theory that introns co-evolved with the ancestral gene, the introns early view, or alternatively, that introns are more recent immigrants. The discovery of an intron in the TPI gene of Culex mosquitoes at a site which was predicted by proponents of the intron early school supported that theory. More recently, the discovery of additional intron sites in several eukaryotes was presented as evidence supporting the introns late school. We have found the 'Culex intron' in two closely related mosquitoes, but not in two more evolutionary primitive Dipterans, suggesting that, if it is an 'ancient intron', loss may be more frequent than that supposed by the intron late school. In addition, we have found that three introns punctuating the TPI gene from the Lepidopteran, Heliothis, appear to be ancestrally related and may be the result of transposable element insertion, 50-90 million years ago. It is argued that both opposing schools in the intron debate be reconciled -- some introns may have been early and certainly others have arrived subsequent to the appearance of the TPI gene.

Molecular cloning of 5-hydroxytryptamine (5-HT) type 1 receptor genes from the Japanese puffer fish, Fugu rubripes

To characterize the structure of Fugu G-protein coupled receptor family and its evolutionary divergence, we have cloned and sequenced the Fugu 5-HT type 1 receptor genes by Polymerase Chain Reaction (PCR) with degenerate primers followed by phage library screening. The analysis of the deduced amino acid sequences showed that F1A alpha and F1A beta have the highest homology to the human 5-HT1A receptor (71.5% and 63.7%, respectively). Another clone, F1D, showed highest (70.5%) homology to the human type 1D receptor. The amino acid residues that are important for ligand binding have been conserved in these Fugu genes. The phylogenetic tree analysis suggests that the duplication event of the Fugu type 1A receptor may have occurred after the divergence of Fugu and the tetrapod lineage.

Testis-specific expression of a functional retroposon encoding glucose-6-phosphate dehydrogenase in the mouse

The X-chromosomal gene glucose-6-phosphate dehydrogenase (G6pd) is known to be expressed in most cell types of mammalian species. In the mouse, we have detected a novel gene, designated G6pd-2, encoding a G6PD isoenzyme. G6pd-2 does not contain introns and appears to represent a retroposed gene. This gene is uniquely transcribed in postmeiotic spermatogenic cells in which the X-encoded G6pd gene is not transcribed. Expression of the G6pd-2 sequence in a bacterial system showed that the encoded product is an active enzyme. Zymogramic analysis demonstrated that recombinant G6PD-2, but not recombinant G6PD-1 (the X-chromosome-encoded G6PD), formed tetramers under reducing conditions. Under the same conditions, G6PD tetramers were also found in extracts of spermatids and spermatozoa, indicating the presence of G6pd-2-encoded isoenzyme in these cell types. G6pd-2 is one of the very few known expressed retroposons encoding a functional protein, and the presence of this gene is probably related to X chromosome inactivation during spermatogenesis.

G protein alpha subunit multigene family in the Japanese puffer fish Fugu rubripes: PCR from a compact vertebrate genome

We compare the complexity and organization of the G protein alpha subunit multigene family in the vertebrate genomes of mammals and the Japanese puffer fish Fugu rubripes. Fourteen Fugu G alpha genes were identified of the 16 genes characterized previously in mammals, including Fugu genes from the four classes of alpha subunits Gs, Gi, Gq, and G12. Fugu and mammalian G alpha coding sequences are highly homologous, and the intron/exon structure of the fish and mammalian orthologs is identical throughout the coding regions. A novel G alpha gene, G alpha p1, was also identified in Fugu rubripes and two other species of puffer fish. The complete sequence of Gnaz and the tandemly duplicated genes Gnai2 and Gnat1 were obtained from a Fugu genomic cosmid library. Introns in the puffer fish G alpha genes lacked repeat DNA sequences, other than simple sequence length repeats, and most introns were significantly shorter in Fugu than in mammalian orthologs. The compact genome of puffer fish provides a unique vertebrate model for characterizing multigene families and identifying novel genes directly from genomic DNA by PCR amplification with degenerate primers. The fact that Fugu encodes most, if not all, of the G protein alpha subunits identified in mammals strongly supports Fugu as a model organism for vertebrate genome research.

Fugu intron oversize reveals the presence of U15 snoRNA coding sequences in some introns of the ribosomal protein S3 gene

We present here the analysis of the genomic organization of the Fugu gene coding for ribosomal protein S3 and its intron encoded U15 RNA, and compare it with the homologous human and Xenopus genes. Only two of the six Fugu S3 gene introns do not contain the U15 sequence and are in fact shorter than 100 nucleotides, as most Fugu introns. The other four introns are somewhat longer and contain sequences homologous to U15 RNA; two of these represent functional copies, as shown by microinjections of Fugu transcripts into Xenopus oocytes, whereas the other two appear to be nonfunctional pseudocopies. Thus Fugu turns out to be ideal for the study of intron encoded snoRNAs, partly because of the reduced cloning and sequencing workload, and partly because the intron length per se can be an indication of the presence of a snoRNA coding sequence.

Conserved linkage between the puffer fish (Fugu rubripes) and human genes for platelet-derived growth factor receptor and macrophage colony-stimulating factor receptor

We have cloned and sequenced the teleost homologs of the human genes encoding platelet-derived growth factor receptor-beta (PDGFR beta) and macrophage colony-stimulating factor 1 receptor (CSFIR) from the puffer fish Fugu rubripes. The Fugu PDGFR beta and CSFIR genes each consist of 21 coding exons similar to the human CSFI gene, but are considerably smaller than their human counterparts because of the smaller introns. Furthermore, the two Fugu genes are linked tandemly in a head-to-tail array similar to their human homologs with 2.2 kb of intergenic sequence. Amino acid sequences of the Fugu and human PDGFR beta and CSFIR genes show an overall homology of 45% and 39%, respectively, with the kinase domains showing a much higher degree of conservation. Dot-matrix analysis revealed several short stretches of conserved sequences in the 3' untranslated regions of the PDGFR beta genes and the adjacent promoter regions of the CSFIR genes. These conserved sequences may have a role in the regulation of expression of either or both of these closely linked genes.

Quality not quantity: the pufferfish genome

The genome of the pufferfish, Fugu rubripes (Fugu) is compact. With a similar gene complement to mammals and a genome size of just 400 Mb, gene density is high averaging one every 6-7 kb. Initial characterization of this genome has shown that although genes are much smaller and more densely spaced, their intron/exon structure is conserved with the resulting introns being small. There is little repetitive DNA in the genome and this greatly facilitates comparative genomic studies. The coding content of genes is highly conserved as are critical regulatory elements of some genes. Other DNA is not, however, and this allows the identification of homologous coding sequence between Fugu and mammalian genes. Although the genome of Fugu is 7.5 times smaller than the human genome, not all genes are reduced proportionately. Some regions of the genome show conserved synteny with mammalian genomes, although at the present time only short physical distances have been examined. The structure of the genome is also being studied. Initial data suggest that this may be different to that found in mammals. It is not clear that the same kind of isochore structure is present in this early vertebrate genome. Patterns of methylation may be different resulting in a different distribution of CpG islands. An attempt is being made to centralize both resources and data from the genome of Fugu so that everything may be integrated into a single, publicly accessible database which in turn, may be integrated with databases from other organisms.

A functional role for some Fugu introns larger than the typical short ones: the example of the gene coding for ribosomal protein S7 and snoRNA U17

The compact genome of Fugu rubripes, with its very small introns, appears to be particularly suitable to study intron-encoded functions. We have analyzed the Fugu gene for ribosomal protein S7 (formerly S8, see Note), whose Xenopus homolog contains in its introns the coding sequences for the small nucleolar RNA U17. Except for intron length, the organization of the Fugu S7 gene is very similar to that of the Xenopus counterpart. The total length of the Fugu S7 gene is 3930 bp, compared with 12691 bp for Xenopus. This length difference is uniquely due to smaller introns. Although short, the six introns are longer than the approximately 100 bp size of most Fugu introns, as they host U17 RNA coding sequences. While four of the six U17 sequences are 'canonical', the remaining two represent diverged U17 pseudocopies. In fact, microinjection in Xenopus oocytes of in vitro synthesized Fugu transcripts containing the 'canonical' U17f sequence results in efficient production of mature U17 RNA, while injection of a transcript containing the U17 psi b sequence does not.

Translocation events in the evolution of aminoacyl-tRNA synthetases

We have characterized hisS, the gene encoding the histidyl-tRNA synthetase (HisRS) from the tetraodontoid fish Fugu rubripes. The hisS gene is about 3.5 kbp long and contains 13 exons and 12 introns of 172 bp, on average. The Fugu hisS gene encodes a putative protein of 519 amino acids with the three motifs identified as signatures of class 2 aminoacyl-tRNA synthetases. A model for the shifting of intron 8 between Fugu and hamster is proposed based on the successive appearance of a cryptic splicing site followed by an insertion mutation that created a new acceptor site. In addition, sequence comparisons suggest that the hisS gene has undergone a translocation through the first intron. As a result, the Fugu HisRS has an N-terminal sequence markedly different from that in the human and hamster enzymes. We propose that similar events have been responsible for variations at the N-terminal end of other aminoacyl-tRNA synthetases. Our analysis suggests that this involves exchanges through introns of two exons encoding an ancestral 32-amino acid motif.

Human lipoprotein lipase last exon is not translated, in contrast to lower vertebrates

We have sequenced the first fish (zebrafish, Brachydanio rerio) lipoprotein lipase (LPL) cDNA clone. Similarities were found in mammalian LPL cDNA, but the codon spanning the last two exons (which is thus split by the last intron) is AGA (Arg) as opposed to TGA in mammals. Exon 10 is thus partially translated. These results were confirmed with rainbow trout (Oncorhynchus mykiss). We also investigated whether mammal TGA coded for selenocystein (SeCys), the 21st amino acid, but found that this was not the case: TGA does not encode SeCys but is a stop codon. It thus appears that the sense codon AGA (fish) has been transformed into a stop codon TGA (human) during the course of evolution. It remains to be determined if the "loss" of the C-terminal end of mammalian LPL protein has conferred an advantage in terms of LPL activity or, on the contrary, a disadvantage (e.g., susceptibility to diabetes or atherosclerosis).

Small is beautiful: comparative genomics with the pufferfish (Fugu rubripes)

As the Human Genome Project advances, it is clear that the emphasis will switch from accumulation of data to their interpretation. Comparative genomics provides a powerful way in which to analyse sequence data. Indeed, there is already a long list of 'model' organisms, which allow comparative analyses in a variety of ways. The very small vertebrate genome of the pufferfish provides a simple and economical way of comparing sequence data from mammals and fish, representing a large evolutionary divergence and so permitting the identification of essential elements that are still present in both species. These elements include genes and the associated machinery that controls their expression; elements that, in many cases, have survived the test of time.

Conservation of synteny between the genome of the pufferfish (Fugu rubripes) and the region on human chromosome 14 (14q24.3) associated with familial Alzheimer disease (AD3 locus)

The genome of the pufferfish (Fugu rubripes) (400 Mb) is approximately 7.5 times smaller than the human genome, but it has a similar gene repertoire to that of man. If regions of the two genomes exhibited conservation of gene order (i.e., were syntenic), it should be possible to reduce dramatically the effort required for identification of candidate genes in human disease loci by sequencing syntenic regions of the compact Fugu genome. We have demonstrated that three genes (dihydrolipoamide succinyltransferase, S31iii125, and S20i15), which are linked to FOS in the familial Alzheimer disease focus (AD3) on human chromosome 14, have homologues in the Fugu genome adjacent to Fugu cFOS. The relative gene order of cFOS, S31iii125, and S20i15 was the same in both genomes, but in Fugu these three genes lay within a 12.4-kb region, compared to >600 kb in the human AD3 locus. These results demonstrate the conservation of synteny between the genomes of Fugu and man and highlight the utility of this approach for sequence-based identification of genes in human disease loci.

Comparative analysis and genomic structure of the tuberous sclerosis 2 (TSC2) gene in human and pufferfish

Germ-line mutations of the TSC2 tumour suppressor gene have been identified in humans with tuberous sclerosis and in the Eker rat. Tuberin, the human TSC2 gene product, has a small region of homology with rap1GAP and stimulates rap1 GTPase activity in vitro, suggesting that one of its cellular roles is to function as a GTPase activating protein (GAP). We have undertaken a comparative analysis of the TSC2 gene in human and the pufferfish, Fugu rubripes. In addition to the GAP domain, three other regions of the proteins are highly conserved (peptide sequence similarity > 80%). These regions are likely to represent further functional domains. To facilitate analysis of mutations within these domains we have determined the genomic structure of the human TSC2 gene. It comprises 41 exons, including exon 31 which was absent from the originally described spliceoform of the human TSC2 transcript and was identified following exon prediction from Fugu genomic sequence. These findings support the proposal of the Fugu genome as a tool for human gene analysis.