Competition may determine the diversity of transposable elements

Transposable elements are genomic parasites that replicate independently from their hosts. They harm their hosts by causing mutations or genomic rearrangements, and most organisms have evolved various mechanisms to suppress their activity. The evolutionary dynamics of transposons in insects, fish, birds and mammals are dramatically different. Mammalian genomes contain few, very abundant but relatively inactive transposon strains, while Drosophila and fish species harbour diverse strains, which typically have low abundance but are much more virulent. We hypothesise that the variation in the diversity and activity of transposable elements between various animal genomes is caused by the differences in the host defence mechanisms against transposon activity. In recent years RNAi, a mechanism capable of gene, virus and transposon silencing has been discovered. We model RNAi as a density dependant mechanism of defence, which can cause competition among transposons depending on its specificity, and test its predictions using the complete Caenorhabditis elegans, Drosophila melanogaster, Fugu rubripes, chicken, mouse, rat and human genome sequences.

Synuclein proteins of the pufferfish Fugu rubripes: sequences and functional characterization

In humans, three genes encode the related alpha-, beta-, and gamma-synucleins, which function as lipid-binding proteins in vitro. They are being widely studied, mainly because of the central involvement of alpha-synuclein in a number of neurodegenerative diseases, including Parkinson's disease, dementia with Lewy bodies, and multiple system atrophy. In these diseases, the normally soluble alpha-synuclein assembles into abnormal filaments. Here, we have identified and characterized the synuclein gene family from the pufferfish Fugu rubripes. It consists of four genes, which encode alpha-, beta-, gamma1-, and gamma2-synucleins. They range from 113 to 127 amino acids in length and share many of the characteristics of human synucleins, including the presence of imperfect amino-terminal repeats of 11 amino acids, a hydrophobic middle region, and a negatively charged carboxy-terminus. All four synucleins are expressed in the Fugu brain. Recombinant Fugu synucleins exhibited differential liposome binding, which was strongest for alpha-synuclein, followed by beta-, gamma2-, and gamma1-synucleins. In assembly experiments, Fugu alpha-, gamma1-, and gamma2-synucleins formed filaments more readily than human alpha-synuclein. Fugu beta-synuclein, by contrast, failed to assemble in bulk. Filament assembly of synucleins was directly proportional to their degree of hydrophobicity and their tendency to form beta-sheet structure, and correlated inversely with their net charge.

Identification, cloning and characterization of a plasma membrane zinc efflux transporter, TrZnT-1, from fugu pufferfish (Takifugu rubripes)

An orthologue of the mammalian ZnT-1 (zinc transporter-1) gene was cloned from the intestine of the torafugu pufferfish (Takifugu rubripes), demonstrating that this gene predates the evolution of land-living vertebrates. TrZnT-1 (T. rubripes ZnT-1) shares overall topology with other members of the ZnT-1 family of zinc transporters, with six TMs (transmembrane domains) including a large histidine-rich intracellular loop between TM IV and V and intracellular C- and N-termini. Expression of TrZnT-1 in a metallothionein acquiescent cell line suggested that this protein reduces intracellular Zn2+ levels. Manipulation of the transporting media showed that several externally applied hydrominerals had no effect on TrZnT-1 activity. However, addition of N-ethylmaleimide increased TrZnT-1-mediated transport, possibly by increasing intracellular free Zn2+ levels by Zn2+ release from carrier proteins. Generation of a specific antibody and subsequent immunocytochemistry on fixed cells overexpressing TrZnT-1 indicated that the protein is localized to the plasma membrane in these cells. The genomic organization of TrZnT-1 is the same as that in mammals with two exons. The upstream regulatory region of the TrZnT-1 gene contains several putative cis-acting elements, including metal-response elements and an Sp1 site. Analysis of the DNA contigs surrounding the TrZnT-1 gene reveal limited synteny between corresponding regions in the rat, mouse and human; however, this was very low, with only two syntenic genes, ZnT-1 and NEK2 (never in mitosis gene A-related kinase).

ACE2 orthologues in non-mammalian vertebrates (Danio, Gallus, Fugu, Tetraodon and Xenopus)

Angiotensin-converting enzyme 2 (ACE2), a newly identified member in the renin–angiotensin system (RAS), acts as a negative regulator of ACE. It is mainly expressed in cardiac blood vessels and the tubular epithelia of kidneys and abnormal expression has been implicated in diabetes, hypertension and heart failure. The mechanism and physiological function of this zinc metallopeptidase in mammals are not yet fully understood. Non-mammalian vertebrate models offer attractive and simple alternatives that could facilitate the exploration of ACE2 function. In this paper we report the in silico analysis of Ace2 genes from the Gallus (chicken), Xenopus (frog), Fugu and Tetraodon (pufferfish) genome assembly databases, and from the Danio (zebrafish) cDNA library. Exon ambiguities of Danio and Xenopus Ace2s were resolved by RT-PCR and 3′RACE. Analyses of the exon–intron structures, alignment, phylogeny and hydrophilicity plots, together with the conserved synteny among these vertebrates, support the orthologous relationship between mammalian and non-mammalian ACE2s. The putative promoters of Ace2 from human, Tetraodon and Xenopus tropicalis drove the expression of enhanced green fluorescent protein (EGFP) specifically in the heart tissue of transgenic Xenopus thus making it a suitable model for future functional genomic studies. Additionally, the search for conserved cis-elements resulted in the discovery of WGATAR motifs in all the putative Ace2 promoters from 7 different animals, suggesting a possible role of GATA family transcriptional factors in regulating the expression of Ace2.

Carbohydrate-binding site of a novel mannose-specific lectin from fugu (Takifugu rubripes) skin mucus

Pufflectin-s, identified in the skin mucus of the fugu Takifugu rubripes, is a novel mannose-specific lectin with similar structure to monocotyledonous plant lectins. In the present study, mutational analysis was used to reveal the mannose-binding sites of pufflectin-s. Putative binding sites were mutated as follows: binding site 1; rPL-D32E (Asp32-->Glu32), rPL-N34S (Asn34-->Ser34) and rPL-V36A (Val36-->Ala36) whereas binding site 2; rPL-D61E (Asp61-->Glu61), rPL-N63S (Asn63-->Ser63) and rPL-V65A (Val65-->Ala65). All recombinant proteins were expressed in Escherichia coli, purified with two chromatographic steps, and then subjected to mannose-binding assay by affinity chromatography. Recombinant wild-type pufflectin-s (rPL-wt) as well as three mutants with changes in binding site 2 could bind to mannose, in contrast to the three mutants with changes in binding site 1 in which mannose-binding activity was completely lost. These results clearly demonstrate that, at the least, binding site 1 is critical to mannose-binding activity in pufflectin-s.

Evolution of cis elements in the differential expression of two Hoxa2 coparalogous genes in pufferfish (Takifugu rubripes)

Sequence divergence in cis-regulatory elements is an important mechanism contributing to functional diversity of genes during evolution. Gene duplication and divergence provide an opportunity for selectively preserving initial functions and evolving new activities. Many vertebrates have 39 Hox genes organized into four clusters (Hoxa-Hoxd); however, some ray-finned fishes have extra Hox clusters. There is a single Hoxa2 gene in most vertebrates, whereas fugu (Takifugu rubripes) and medaka (Oryzias latipes) have two coparalogous genes [Hoxa2(a) and Hoxa2(b)]. In the hindbrain, both genes are expressed in rhombomere (r) 2, but only Hoxa2(b) is expressed in r3, r4, and r5. Multiple regulatory modules directing segmental expression of chicken and mouse Hoxa2 genes have been identified, and each module is composed of a series of discrete elements. We used these modules to investigate the basis of differential expression of duplicated Hoxa2 genes, as a model for understanding the divergence of cis-regulatory elements. Therefore, we cloned putative regulatory regions of the fugu and medaka Hoxa2(a) and -(b) genes and assayed their activity. We found that these modules direct reporter expression in a chicken assay, in a manner corresponding to their endogenous expression pattern in fugu. Although sequence comparisons reveal many differences between the two coparalogous genes, specific subtle changes in seven cis elements of the Hoxa2(a) gene restore segmental regulatory activity. Therefore, drift in subsets of the elements in the regulatory modules is responsible for the differential expression of the two coparalogous genes, thus providing insight into the evolution of cis elements.

Identification of an osteocalcin isoform in fish with a large acidic prodomain

Osteocalcin is a small, secreted bone protein whose gene consists of four exons. In the course of analyzing the structure of fish osteocalcin genes, we recently found that the spotted green pufferfish has two possible exon 2 structures, one of 15 bp and the other of 324 bp. Subsequent analysis of the pufferfish cDNA showed that only the transcript with a large exon 2 exists. Exon 2 codes for the osteocalcin propeptide, and exon 2 of pufferfish osteocalcin is approximately 3.4-fold larger than exon 2 previously found in other vertebrate species. We have termed this new pufferfish osteocalcin isoform OC2. Additional studies showed that the OC2 isoform is restricted to a unique fish taxonomic group, the Osteichthyes; OC2 is the only osteocalcin isoform found so far in six Osteichthyes species, whereas both OC1 and OC2 isoforms coexist in zebrafish and rainbow trout. The larger size of the OC2 propeptide is due to an acidic region that is likely to be highly phosphorylated and has no counterpart in the OC1 propeptide. We propose 1) that OC1 and OC2 are encoded by distinct genes that originated from a duplication event that probably occurred in the teleost fish lineage soon after divergence from tetrapods and 2) that the novel OC2 propeptide could be, if secreted, a phosphoprotein that participates in the regulation of biomineralization through its large acidic and phosphorylated propeptide.

Conservation of RET regulatory function from human to zebrafish without sequence similarity

Evolutionary sequence conservation is an accepted criterion to identify noncoding regulatory sequences. We have used a transposon-based transgenic assay in zebrafish to evaluate noncoding sequences at the zebrafish ret locus, conserved among teleosts, and at the human RET locus, conserved among mammals. Most teleost sequences directed ret-specific reporter gene expression, with many displaying overlapping regulatory control. The majority of human RET noncoding sequences also directed ret-specific expression in zebrafish. Thus, vast amounts of functional sequence information may exist that would not be detected by sequence similarity approaches.

Ancient duplicated conserved noncoding elements in vertebrates: a genomic and functional analysis

Fish-mammal genomic comparisons have proved powerful in identifying conserved noncoding elements likely to be cis-regulatory in nature, and the majority of those tested in vivo have been shown to act as tissue-specific enhancers associated with genes involved in transcriptional regulation of development. Although most of these elements share little sequence identity to each other, a small number are remarkably similar and appear to be the product of duplication events. Here, we searched for duplicated conserved noncoding elements in the human genome, using comparisons with Fugu to select putative cis-regulatory sequences. We identified 124 families of duplicated elements, each containing between two and five members, that are highly conserved within and between vertebrate genomes. In 74% of cases, we were able to assign a specific set of paralogous genes with annotation relating to transcriptional regulation and/or development to each family, thus removing much of the ambiguity in identifying associated genes. We find that duplicate elements have the potential to up-regulate reporter gene expression in a tissue-specific manner and that expression domains often overlap, but are not necessarily identical, between family members. Over two thirds of the families are conserved in duplicate in fish and appear to predate the large-scale duplication events thought to have occurred at the origin of vertebrates. We propose a model whereby gene duplication and the evolution of cis-regulatory elements can be considered in the context of increased morphological diversity and the emergence of the modern vertebrate body plan.

Two interleukin (IL)-15 homologues in fish from two distinct origins

Here, we report two distinct genes in teleosts that are homologous to interleukin (IL)-15. The two genes, isolated from fugu (Takifugu rubripes), resemble to mammalian IL-15 but differ from IL-2 and IL-21 in their amino acid sequences, the possessing of an extraordinary long signal peptide and more widespread tissue localization. In addition, multiple out-of-frame AUG codons, the negative translational regulators of mammalian IL-15 genes were also detected in the 5'-UTR of the two genes. Fugu IL-15 homologues also contain four conserved cysteines allowing the formation of two disulfide bridges along with four predicted alpha-helices. Genomic analysis showed that one of the fugu IL-15 homologues possessed six coding exons and exhibited a similar exon-intron organization and synteny structure to that of mammalian and chicken IL-15 genes. Conversely, the other fugu IL-15 homologue possesses four exons and exhibits a different synteny structure with that of IL-15, suggesting that the two genes were derived from two different origins. Moreover, the two genes also differ from each other in tissue localizations and in their expression in response to mitogens. The existence of these two IL-15 homologues in telesots was further supported by their characterization in zebrafish Danio rerio, and the green-spotted pufferfish Tetraodon nigroviridis. The discovery of two distinct IL-15 homologues in fish will assist investigations into the evolution of these genes and their relative contribution to the fish immune system.

B-lymphocyte-induced maturation protein-1 (Blimp-1) gene of torafugu (Takifugu rubripes)

Homologous gene of B-lymphocyte-induced maturation protein-1 (Blimp-1) of torafugu (Takifugu rubripes) was identified by tblast search analysis. RT-PCR and 5'RACE clearly defined the sequence of the UTR and coding region which has been ambiguously determined by tblast analysis. Fugu Blimp-1 was mainly expressed in the lymphoid organs. These finding imply that Blimp-1 would take a major role in the terminal differentiation of B-cells to plasma cells in fish.

Cloning and expression analysis of interferon-gamma-inducible-lysosomal thiol reductase gene in large yellow croaker (Pseudosciaena crocea)

In mammals, interferon-gamma-inducible-lysosomal thiol reductase (GILT) has been demonstrated to play a key role in the processing and presentation of MHC class II-restricted antigen (Ag) by catalyzing disulfide bond reduction, thus unfolding native protein Ag and facilitating subsequent cleavage by proteases. Here, we reported the cloning of a GILT gene homologue from the spleen of large yellow croaker, a marine fish (LycGILT). The full-length cDNA of LycGILT gene is 1033 nucleotides (nt) encoding a protein of 256 amino acids (aa), with a putative molecular weight of 28.9 kDa. The deduced protein is highly homologous to that of mammalian and zebrafish GILTs and shares 54.1% sequence identity to that of zebrafish and 43.2-39.2% sequence identity to that of various mammals. The deduced LycGILT possesses the typical structural feature of mammalian GILT, including an active-site CXXC motif, a GILT signature sequence CQHGX2ECX2NX4C, and other six cysteines responsible for the formation of disulfide bonds in the C-terminus. Genomic analysis revealed that LycGILT gene, spanning a 3159nt fragment, contained seven exons interrupted by six introns and exhibited a similar exon-intron organization to human and mouse GILT genes except for a slightly more compact intron arrangement. The LycGILT expression is obviously up-regulated in spleen and kidney after immunization with inactivated trivalent bacterial vaccine consisting of Vibrio alginolyticus, V. paraphaemolyticus, and Aeromonas hydrophila although it also is constitutively expressed in liver, gills, brain, and heart, suggesting that LycGILT may be involved in the immune response to bacterial challenge in large yellow croaker. A search of NCBI sequence data with LycGILT cDNA identified a pufferfish (fugu rubrides) GILT homologue cDNA and its genomic DNA sequence, where two putative interferon-gamma activation sites (GAS) were found within the promoter region. This provided evidence that a fish GILT homologue like mammalian GILT, may also be regulated by interferon-gamma (IFN-gamma) through the JAK-STAT signal pathway. These results indicate that the bony fish GILT is a functional homologue of mammalian GILT.

Organization and structure of hox gene loci in medaka genome and comparison with those of pufferfish and zebrafish genomes

We isolated BAC clones that cover the entire hox gene loci in the medaka fish Oryzias latipes. The BAC clones were characterized by the Southern hybridization with many hox gene probes isolated in our previous study and by PCR using primers designed for selective amplification of respective hox genes. Then, the BAC clones have been subjected to shotgun sequencing. The results revealed the organization of the entire hox gene loci. Forty-six hox genes in total are encoded in seven clusters as follows: 10 hox genes in Aa cluster; 5 in Ab; 9 in Ba; 4 in Bb; 10 in Ca; 6 in Da; and 2 in Db. Together with the information on the hox gene loci registered in the Fugu genome database and in the Danio genome database, the physical maps of three fish genomes were constructed and compared one another. Not only numbers of hox genes but also the distances between the neighboring hox genes are highly similar between medaka and fugu. As for six clusters, Aa, Ab, Ba, Bb, Ca and Da that are commonly present in the three fishes, only few or no differences were found in each cluster. Thus, the hox gene sets should have been well conserved once they had been established in respective species.

Many genes in fish have species-specific asymmetric rates of molecular evolution

BACKGROUND

Gene and genome duplication events increase the amount of genetic material that might then contribute to an increase in the genomic and phenotypic complexity of organisms during evolution. Thus, it has been argued that there is a relationship between gene copy number and morphological complexity and/or species diversity. This hypothesis implies that duplicated genes have subdivided or evolved novel functions compared to their pre-duplication proto-orthologs. Such a functional divergence might be caused by an increase in evolutionary rates in one ortholog, by changes in expression, regulatory evolution, insertion of repetitive elements, or due to positive Darwinian selection in one copy. We studied a set of 2466 genes that were present in Danio rerio, Takifugu rubripes, Tetraodon nigroviridis and Oryzias latipes to test (i) for forces of positive Darwinian selection; (ii) how frequently duplicated genes are retained, and (iii) whether novel gene functions might have evolved.

RESULTS

25% (610) of all investigated genes show significantly smaller or higher genetic distances in the genomes of particular fish species compared to their human ortholog than their orthologs in other fish according to relative rate tests. We identified 49 new paralogous pairs of duplicated genes in fish, in which one of the paralogs is under positive Darwinian selection and shows a significantly higher rate of molecular evolution in one of the four fish species, whereas the other copy apparently did not undergo adaptive changes since it retained the original rate of evolution. Among the genes under positive Darwinian selection, we found a surprisingly high number of ATP binding proteins and transcription factors.

CONCLUSION

The significant rate difference suggests that the function of these rate-changed genes might be essential for the respective fish species. We demonstrate that the measurement of positive selection is a powerful tool to identify divergence rates of duplicated genes and that this method has the capacity to identify potentially interesting candidates for adaptive gene evolution.

Fugu genome does not contain mitochondrial pseudogenes

Contrary to previous observations that fish genomes are devoid of nuclear mitochondrial pseudogenes, a genome-wide survey identified a large number of "recent" and "ancient" nuclear mitochondrial DNA fragments (Numts) in the whole-genome sequences of the fugu (Takifugu rubripes), Tetraodon nigroviridis, and zebrafish (Danio rerio). We have analyzed the latest assembly (v4.0) of the fugu genome and show that, like the Anopheles genome, the fugu nuclear genome does not contain mitochondrial pseudogenes. Fugu assembly v4.0 contains a single scaffold representing the near complete sequence of the fugu mitochondria. The "recent" Numts identified by the previous study in fugu assembly v2.0 are in fact shotgun sequences of mitochondrial DNA that were misassembled with the nuclear sequences, whereas the "ancient" Numts appear to be the result of spurious matches. It is likely that the Numts identified in the genomes of Tetraodon and zebrafish are also similar artifacts. Shotgun sequences of whole genomes often include some mitochondrial sequences. Therefore, any Numts identified in shotgun-sequence assemblies should be verified by Southern hybridization or PCR amplification.

The co-existence of two growth hormone receptors in teleost fish and their differential signal transduction, tissue distribution and hormonal regulation of expression in seabream

Two genomic contigs of putative growth hormone receptors (GHRs) were identified in fugu and zebrafish genomes by in silico analysis, suggesting the presence of two GHR subtypes in a single teleost species. We have tested this hypothesis by cloning the full-length cDNA sequence of a second GHR subtype from the black seabream in which the first GHR subtype had been previously reported by us. In addition, we had also cloned the sequences of both GHR subtypes from two other fish species, namely the Southern catfish and the Nile tilapia. Phylogenetic analysis of known GHR sequences from various vertebrates revealed that fish GHRs cluster into two distinct clades, viz. GHR1 and GHR2. One clade (GHR1), containing 6 to 7 extracellular cysteine residues, is structurally more akin to the non-teleost GHRs. The other clade (GHR2), containing only 4 to 5 extracellular cysteine residues, is unique to teleosts and is structurally more divergent from the non-teleost GHRs. In addition, we had examined the biological activities of both GHR subtypes from seabream using a number of reporter transcription assays in cultured eukaryotic cells and demonstrated that both of them were able to activate the Spi 2.1 and beta-casein promoters upon receptor stimulation in a ligand specific manner. In contrast, only GHR1 but not GHR2 in seabream could trigger the c-fos promoter activity, indicating that the two GHR subtypes possess some differences in their signal transduction mechanisms. Also, the expression of GHR2 is significantly higher than GHR1 in many tissues of the seabream including the gonad, kidney, muscle, pituitary and spleen. In vivo hormone treatment data indicated that cortisol upregulated hepatic GHR1 expression in seabream but not GHR2, whereas testosterone decreased hepatic GHR2 expression but not GHR1. On the other hand, hepatic expression of both GHR1 and GHR2 in seabream was decreased by estradiol treatment.

Oct2 transcription factors in fish--a comparative genomic analysis

The Oct2 transcription factor is important in driving expression of the IgH locus of the channel catfish, Ictalurus punctatus. Two isoforms, catfish Oct2alpha and Oct2beta, have been characterized at the level of expression and function, but little is known of the structure of the Oct2 gene in catfish. To gain insight into the diversity of Oct2 gene structure and expression in the teleost fish, a comparative genomic analysis of Oct2 was undertaken in the pufferfish (Fugu rubripes) and the zebrafish (Danio rerio). The orthologues of zebrafish and Fugu Oct2 were identified, and share with catfish Oct2 the expression of a limited number (two in zebrafish, three in Fugu) of isotypes produced by alternative pathways of RNA processing. The alternatively spliced variants of catfish Oct2 showed a different pattern of exon use from those of Fugu and zebrafish. The analysis also identified a novel homologue of Oct2 in both zebrafish and Fugu. This homologue, termed Oct2x, shares similarities to both Oct1 and Oct2. A phylogenetic analysis of the relationships of Oct2x gave an unexpected result, with Oct2x occupying a position basal to the Oct gene families of both vertebrates and Drosophila.

Ray-Fin Fish Tetraploidization Gave Rise to Pufferfish Duplicates of NPY and PYY, but Zebrafish NPY Duplicate Was Lost

We have used sequence information and gene location to identify NPY family genes in the pufferfish, Takifugu rubripes (fugu), and zebrafish. Fugu has two copies of NPY, presumably resulting from the ray-fin fish tetraploidization. Zebrafish has probably lost one of the copies. Both species have two copies of PYY, the second of which was previously named PY. The two fugu NPY genes are predominantly expressed in brain. The two PYY genes are expressed in a broad range of tissues including brain and gonads. Thus, the NPY system appears to be more complex in teleosts than in tetrapods.

Unusual Genomic Structure: Melanocortin Receptors inFugu

The melanocortin (MC) receptors are found in five subtypes in mammals and chicken, while recent studies have shown that the Fugu (Takifugu rubripes) genome has only four MC receptors and the zebrafish genome has six subtypes. The MC3 receptor seems to be missing from the two closely related pufferfishes, Fugu and Tetraodon (Tetraodon nigroviridis). The MC2 and MC5 receptors in the pufferfish have introns. Moreover, these two receptors are found in a tandem that is remarkably conserved in several vertebrate species. Here, we speculate about the genomic origin of the MC receptors.

Nested genes in the human genome

Here we studied one special type of gene, i.e., the nested gene, in the human genome. We collected 373 reliably annotated nested genes. Two-thirds of them were on the strand opposite that of their host gene. About 58% coding nested gene pairs were conserved in mouse and some were even maintained in chicken and fish, while nested pseudogenes were poorly conserved. Ka/Ks analysis revealed that nested genes were under strong selection, although they did not demonstrate greater conservation than other genes. With microarray data we observed that two partners of one nested pair seemed to be expressed reciprocally. A significant proportion of nested genes were tissue-specifically expressed. Gene ontology analysis demonstrated that quite a number of nested genes participated in cellular signal transduction. Based on these observations, we think that nested genes are a group of genes with important physiological functions.

Novel bioactive parathyroid hormone and related peptides in teleost fish

We report the identification, gene expression and biological activity of two parathyroid hormones (PTH; PTHA and PTHB), two PTH-related peptides (PTHrP; PTHrPA and PTHrPB) and a PTH-like ligand (PTH-L) with hybrid characteristics in puffer fishes (Takifugu rubripes and Tetraodon fluviatilis). Experimental data are consistent with PTH-L and PTHrPA having calciotropic activities equivalent, respectively, to tetrapod PTH and PTHrP. We hypothesise on the basis of phylogenetic and functional analysis that PTH-L could be a fish relic of an ancestral PTH/PTHrP gene.

Comparative genomics of the human and Fugu voltage-gated calcium channel alpha1-subunit gene family reveals greater diversity in Fugu

Extensive search for the orthologs of 10 human voltage-gated calcium channel (VGCC) alpha(1)-subunit genes in the Fugu genome sequence revealed 21 alpha(1)-subunit genes in the compact genome of Fugu. Subtype classification of the identified Fugu alpha(1) orthologs based on phylogenetic analysis, genomic organization and sequence comparison of the most divergent II/III loop and the C-terminal regions of the alpha(1)-subunits indicated extra copies of alpha(1S)-, alpha(1D)-, alpha(1F)-, alpha(1A)-, alpha(1E)-, alpha(1H)- and alpha(1G)-subunit genes. Phylogenetic analysis reveals that this is likely due to fish lineage specific alpha(1)-subunit subtype duplication. Sequence comparison shows that many of the structural features characteristic of VGCC and specific channel subtypes are also present in the Fugu alpha(1)-subunits. All the Fugu alpha(1)-subunits showed similar expression profile to that of the mammalian alpha(1)-subunits except for Fugu alpha(1S), alpha(1A), alpha(1B) and alpha(1H) which have a more widespread tissue distribution. These results indicate that Fugu, a lower vertebrate, has more extensive channel heterogeneity compared to human.

The discovery, positioning and verification of a set of transcription-associated motifs in vertebrates

We have developed several new methods to investigate transcriptional motifs in vertebrates. We developed a specific alignment tool appropriate for regions involved in transcription control, and exhaustively enumerated all possible 12-mers for involvement in transcription by virtue of their mammalian conservation. We then used deeper comparative analysis across vertebrates to identify the active instances of these motifs. We have shown experimentally in Medaka fish that a subset of these predictions is involved in transcription.

Defining a genomic radius for long-range enhancer action: duplicated conserved non-coding elements hold the key

Many conserved non-coding elements (CNEs) in vertebrate genomes have been shown to function as tissue-specific enhancers. However, the target genes of most CNEs are unknown. Here we show that the target genes of duplicated CNEs can be predicted by considering their neighbouring paralogous genes. This enables us to provide the first systematic estimate of the genomic range for distal cis-regulatory interactions in the human genome: half of CNEs are >250 kb away from their associated gene.

Discovery of a new class of immunoglobulin heavy chain from fugu

In teleosts, the genomic organization of the immunoglobulin (Ig) heavy (H)-chain locus was thought to follow a typical translocon-type multigene structure; however, recent studies have indicated a variation in the structure and this might be teleost specific. Isotypes of the Ig H-chain, namely IgM, IgD, IgZ and IgT, have been identified. In this study, we report the discovery of a new class of IgH from fugu. This isotype was first identified from the genomic sequence of the fugu IgH locus. This novel IgH gene is composed of two constant (C) domains, a hinge region, and two exons encoding membrane regions. Surprisingly, the new IgH gene is present between the variable (V)H and Cmu regions of the locus. The C domains of the new isotype do not show any significant similarity to mammalian or fish IgH genes. The cloned cDNA from the new isotype has typical Ig H-chain characteristics and is expressed as both secretory and membrane form. Transcript analyses suggest that the new IgH from fugu might only use the joining (J)H segments present in front of the new CH domains and that the usage of DH and JH segments is specific to the isotype expressed. The expression pattern of the gene has been confirmed by in situ hybridization and PCR studies.