Fugu and human sequence comparison identifies novel human genes and conserved non-coding sequences

Exploring the Role of SMPD3 in the lncRNA-miRNA-mRNA Regulatory Network in TBI Progression by Influencing Energy Metabolism

Background

Traumatic brain injury (TBI) is associated with disturbances in energy metabolism. This study aimed to construct a lncRNA-miRNA-mRNA network through bioinformatics methods to explore energy metabolism-related genes in the pathogenesis of TBI.

Methods

Data from datasets GSE171718, GSE131695, and GSE223245 obtained from the Gene Expression Omnibus, were analyzed to identify differentially expressed (DE) genes. Regulatory relationships were investigated through miRDB, miRTarBase, and TargetScan, thereby forming a lncRNA-miRNA-mRNA network. The Molecular Signatures Database (MSigDB) was utilized to identify energy metabolism-related genes, and a protein-protein interaction (PPI) network was established through the STRING database. Functional annotation and enrichment analysis were conducted using GO and KEGG. The TBI mouse model was established to detect the expression levels of GOLGA8B, ZNF367, and SMPD3 in brain tissues.

Results

SMPD3 emerged as the key DE gene linked to energy metabolism in TBI, demonstrating a negative correlation with miR-218-5p and being associated with moderate unconsciousness and female patients. The PPI network revealed SMPD3 interactions with proteins associated with cell death, sphingolipid metabolism, and neurodegenerative diseases such as Alzheimer's disease. In vivo, GOLGA8B, ZNF367, and SMPD3 mRNA levels were significantly lower in TBI mice.

Conclusion

In summary, SMPD3 represents a crucial metabolic gene in the progression of TBI. It potentially provides a new therapeutic target for metabolic disorders caused by traumatic brain injury (TBI) and holds significant theoretical value for further research.

Inhibition of zinc finger protein 367 exerts a tumor suppressive role in colorectal cancer by affecting the activation of oncogenic YAP1 signaling

Zinc finger protein 367 (ZNF367) has been documented as a new cancer-related protein that exerts a pivotal role in the carcinogenesis of multiple cancers. However, whether ZNF367 plays a role in colorectal cancer has not been fully understood. Our data showed that ZNF367 expression was higher in colorectal cancer. Depletion of ZNF367 weakened the proliferative and invasive capabilities of colorectal cancer cells. Up-regulation of ZNF367 enhanced the in vitro malignant features of colorectal cancer cells. Knockdown of ZNF367 impeded the activation of Yes-associated protein (YAP1). Reactivation of YAP1 reversed the ZNF367-knockdown-mediated anticancer effects. Suppression of YAP1 significantly abolished ZNF367-overexpression-induced tumor-promotion effects. Depletion of ZNF367 repressed the tumorigenicity of colorectal cancer cells in vivo. These findings demonstrate that ZNF367 is overexpressed in colorectal cancer and acts as a potential tumor-promoter that contributes to the proliferation and invasion of colorectal cancer by enhancing the activation of YAP1 signaling.

Embryonic development of the grass pufferfish (Takifugu niphobles): From egg to larvae

Tetraodontidae (pufferfish) family members carry the smallest genomes among vertebrates, and these pocket-sized genomes have directly contributed to our understanding of the structure and evolution of higher animals. The grass pufferfish (Takifugu niphobles) could be considered a potential new model organism for comparative genomics and development due to the potential access to embryos, and availability of sequence data for two similar genomes: that of spotted green pufferfish (Tetraodon nigroviridis) and Fugu (Takifugu rubripes). In this study, we provide the first description of the normal embryonic development of T. niphobles, by drawing comparisons with the closely related species cited above. Embryos were obtained by in vitro fertilization of eggs, and subsequent development was monitored at a constant temperature consistent with natural conditions. T. niphobles development was divided into seven periods of embryogenesis: the zygote, cleavage, blastula, gastrula, segmentation, pharyngula, and hatching periods; and stages subdividing these periods are defined based on morphological characteristics. The developmental stage series described in this study aims to provide the utilization of T. niphobles as an experimental model organism for comparative developmental studies.

Chemokine receptors in Atlantic salmon

Teleost sequence data have revealed that many immune genes have evolved differently when compared to other vertebrates. Thus, each gene family needs functional studies to define the biological role of individual members within major species groups. Chemokine receptors, being excellent markers for various leukocyte subpopulations, are one such example where studies are needed to decipher individual gene function. The unique salmonid whole genome duplication that occurred approximately 95 million years ago has provided salmonids with many additional duplicates further adding to the complexity and diversity. Here we have performed a systematic study of these receptors in Atlantic salmon with particular focus on potential inflammatory receptors. Using the preliminary salmon genome data we identified 48 chemokine or chemokine-like receptors including orthologues to the ten receptors previously published in trout. We found expressed support for 40 of the bona fide salmon receptors. Eighteen of the chemokine receptors are duplicated, and when tested against a diploid sister group the majority were shown to be remnants of the 4R whole genome duplication with subsequent high sequence identity. The salmon chemokine receptor repertoire of 40 expressed bona fide genes is comparably larger than that found in humans with 23 receptors. Diversification has been a major driving force for these duplicate genes with the main variability residing in ligand binding and signalling domains.

ZNF367 inhibits cancer progression and is targeted by miR-195

Background

Several members of the zinc finger protein family have been recently shown to have a role in cancer initiation and progression. Zinc finger protein 367 (ZNF367) is a member of the zinc finger protein family and is expressed in embryonic or fetal erythroid tissue but is absent in normal adult tissue.

Methodology/Principal Findings

We show that ZNF367 is overexpressed in adrenocortical carcinoma, malignant pheochromocytoma/paraganglioma and thyroid cancer as compared to normal tissue and benign tumors. Using both functional knockdown and ectopic overexpression in multiple cell lines, we show that ZNF367 inhibits cellular proliferation, invasion, migration, and adhesion to extracellular proteins in vitro and in vivo. Integrated gene and microRNA expression analyses showed an inverse correlation between ZNF367 and miR-195 expression. Luciferase assays demonstrated that miR-195 directly regulates ZNF367 expression and that miR-195 regulates cellular invasion. Moreover, integrin alpha 3 (ITGA3) expression was regulated by ZNF367.

Conclusions/Significance

Our findings taken together suggest that ZNF367 regulates cancer progression.

Expression of Wnt signaling skeletal development genes in the cartilaginous fish, elephant shark (Callorhinchus milii)

Jawed vertebrates (Gnasthostomes) are broadly separated into cartilaginous fishes (Chondricthyes) and bony vertebrates (Osteichthyes). Cartilaginous fishes are divided into chimaeras (e.g. ratfish, rabbit fish and elephant shark) and elasmobranchs (e.g. sharks, rays and skates). Both cartilaginous fish and bony vertebrates are believed to have a common armoured bony ancestor (Class Placodermi), however cartilaginous fish are believed to have lost bone. This study has identified and investigated genes involved in skeletal development in vertebrates, in the cartilaginous fish, elephant shark (Callorhinchus milii). Ctnnb1 (β-catenin), Sfrp (secreted frizzled protein) and a single Sost or Sostdc1 gene (sclerostin or sclerostin domain-containing protein 1) were identified in the elephant shark genome and found to be expressed in a number of tissues, including cartilage. β-catenin was also localized in several elephant shark tissues. The expression of these genes, which belong to the Wnt/β-catenin pathway, is required for normal bone formation in mammals. These findings in the cartilaginous skeleton of elephant shark support the hypothesis that the common ancestor of cartilaginous fishes and bony vertebrates had the potential for making bone.

Zar1 represses translation in Xenopus oocytes and binds to the TCS in maternal mRNAs with different characteristics than Zar2

Maternal mRNAs are translationally regulated during early development. Zar1 and its closely related homolog, Zar2, are both crucial in early development. Xenopus laevis Zygote arrest 2 (Zar2) binds to the Translational Control Sequence (TCS) in maternal mRNAs and regulates translation. The molecular mechanism of Zar1 has not been described. Here we report similarities and differences between Xenopus Zar1 and Zar2. Analysis of Zar sequences in vertebrates revealed two Zar family members with conserved, characteristic amino acid differences in the C-terminal domain. The presence of only two vertebrate Zar proteins was supported by analyzing Zar1 synteny. We propose that the criteria for naming Zar sequences are based on the characteristic amino acids and the chromosomal context. We also propose reclassification of some Zar sequences. We found that Zar1 is expressed throughout oogenesis and is stable during oocyte maturation. The N-terminal domain of Zar1 repressed translation of a reporter construct in immature oocytes. Both Zar1 and Zar2 bound to the TCS in the Wee1 and Mos 3' UTRs using a zinc finger in the C-terminal domain. However, Zar1 had much higher affinity for RNA than Zar2. To show the functional significance of the conserved amino acid substitutions, these residues in Zar2 were mutated to those found in Zar1. We show that these residues contributed to the different RNA binding characteristics of Zar1 compared to Zar2. Our study shows that Zar proteins have generally similar molecular functions in the translational regulation of maternal mRNAs, but they may have different roles in early development.

Repertoire of Protein Kinases Encoded in the Genome of Takifugu rubripes

Takifugu rubripes is teleost fish widely used in comparative genomics to understand the human system better due to its similarities both in number of genes and structure of genes. In this work we survey the fugu genome, and, using sensitive computational approaches, we identify the repertoire of putative protein kinases and classify them into groups and subfamilies. The fugu genome encodes 519 protein kinase-like sequences and this number of putative protein kinases is comparable closely to that of human. However, in spite of its similarities to human kinases at the group level, there are differences at the subfamily level as noted in the case of KIS and DYRK subfamilies which contribute to differences which are specific to the adaptation of the organism. Also, certain unique domain combination of galectin domain and YkA domain suggests alternate mechanisms for immune response and binding to lipoproteins. Lastly, an overall similarity with the MAPK pathway of humans suggests its importance to understand signaling mechanisms in humans. Overall the fugu serves as a good model organism to understand roles of human kinases as far as kinases such as LRRK and IRAK and their associated pathways are concerned.

Dendritic cell-specific transcript: dendritic cell marker and regulator of TLR-induced cytokine production

Dendritic cells (DCs) are the professional APCs of the immune system that dictate the type and course of an immune response. Molecular understanding of DC biology is important for the design of DC-based immunotherapies and optimal clinical applications in vaccination settings. Previously, we isolated and characterized the cDNA-encoding dendritic cell-specific transcript (DC-SCRIPT; also known as ZNF366). DC-SCRIPT mRNA expression in the immune system was confined to DCs and was reported to be an early hallmark of DC differentiation. In this study, we demonstrate IL-4 to be the dominant factor for DC-SCRIPT expression in human monocyte-derived DCs. In addition, to our knowledge, we show for the first time endogenous DC-SCRIPT protein expression in human DCs both in vitro and in situ. DC-SCRIPT protein is detected early upon differentiation of monocytes into DCs and is also present in multiple freshly isolated DC subsets. Maturation of DCs with TLR ligands further increased DC-SCRIPT mRNA expression, suggesting a role in DC maturation. Indeed, small interfering RNA-mediated knockdown of DC-SCRIPT affected the cytokine response upon TLR stimulation. These DCs displayed enhanced IL-10 and decreased IL-12 production, compared with wild-type DCs. Silencing of IL-10 in DC-SCRIPT knockdown DCs rescued IL-12 expression, suggesting a primary role for DC-SCRIPT in the regulation of IL-10 production.

Genome-wide identification of hypoxia-inducible factor binding sites and target genes by a probabilistic model integrating transcription-profiling data and in silico binding site prediction

The transcriptional response driven by Hypoxia-inducible factor (HIF) is central to the adaptation to oxygen restriction. Hence, the complete identification of HIF targets is essential for understanding the cellular responses to hypoxia. Herein we describe a computational strategy based on the combination of phylogenetic footprinting and transcription profiling meta-analysis for the identification of HIF-target genes. Comparison of the resulting candidates with published HIF1a genome-wide chromatin immunoprecipitation indicates a high sensitivity (78%) and specificity (97.8%). To validate our strategy, we performed HIF1a chromatin immunoprecipitation on a set of putative targets. Our results confirm the robustness of the computational strategy in predicting HIF-binding sites and reveal several novel HIF targets, including RE1-silencing transcription factor co-repressor (RCOR2). In addition, mapping of described polymorphisms to the predicted HIF-binding sites identified several single-nucleotide polymorphisms (SNPs) that could alter HIF binding. As a proof of principle, we demonstrate that SNP rs17004038, mapping to a functional hypoxia response element in the macrophage migration inhibitory factor (MIF) locus, prevents induction of this gene by hypoxia. Altogether, our results show that the proposed strategy is a powerful tool for the identification of HIF direct targets that expands our knowledge of the cellular adaptation to hypoxia and provides cues on the inter-individual variation in this response.

Balanced gene regulation by an embryonic brain ncRNA is critical for adult hippocampal GABA circuitry

Genomic studies demonstrate that while the majority of the mammalian genome is transcribed, only about 2% of these transcripts are protein coding. We have been investigating how the long, polyadenylated Evf2 non-coding RNA regulates transcription of homeodomain transcription factors DLX5 and DLX6 in the developing mouse forebrain. Here we show that in developing ventral forebrain, Evf2 recruits DLX and MECP2 transcription factors to key DNA regulatory elements in the Dlx 5/6 intergenic region and controls Dlx5, Dlx6, and GAD67 expression through trans and cis-acting mechanisms. Evf2 mouse mutants have reduced numbers of GABAergic interneurons in early post-natal hippocampus and dentate gyrus. Although the numbers of GABAergic interneurons and GAD67 RNA levels return to normal in Evf2 mutant adult hippocampus, reduced synaptic inhibition occurs. These results suggest that non-coding RNA-dependent balanced gene regulation in embryonic brain is critical for proper formation of GABA-dependent neuronal circuitry in adult brain.

Extensive expansion and diversification of the chemokine gene family in zebrafish: identification of a novel chemokine subfamily CX

Background

The chemokine family plays important roles in cell migration and activation. In humans, at least 44 members are known. Based on the arrangement of the four conserved cysteine residues, chemokines are now classified into four subfamilies, CXC, CC, XC and CX3C. Given that zebrafish is an important experimental model and teleost fishes constitute an evolutionarily diverse group that forms half the vertebrate species, it would be useful to compare the zebrafish chemokine system with those of mammals. Prior to this study, however, only incomplete lists of the zebrafish chemokine genes were reported.

Results

We systematically searched chemokine genes in the zebrafish genome and EST databases, and identified more than 100 chemokine genes. These genes were CXC, CC and XC subfamily members, while no CX3C gene was identified. We also searched chemokine genes in pufferfish fugu and Tetraodon, and found only 18 chemokine genes in each species. The majority of the identified chemokine genes are unique to zebrafish or teleost fishes. However, several groups of chemokines are moderately similar to human chemokines, and some chemokines are orthologous to human homeostatic chemokines CXCL12 and CXCL14. Zebrafish also possesses a novel species-specific subfamily consisting of five members, which we term the CX subfamily. The CX chemokines lack one of the two N-terminus conserved cysteine residues but retain the third and the fourth ones. (Note that the XC subfamily only retains the second and fourth of the signature cysteines residues.) Phylogenetic analysis and genome organization of the chemokine genes showed that successive tandem duplication events generated the CX genes from the CC subfamily. Recombinant CXL-chr24a, one of the CX subfamily members on chromosome 24, showed marked chemotactic activity for carp leukocytes. The mRNA was expressed mainly during a certain period of the embryogenesis, suggesting its role in the zebrafish development.

Conclusion

The phylogenic and genomic organization analyses suggest that a substantial number of chemokine genes in zebrafish were generated by zebrafish-specific tandem duplication events. During such duplications, a novel chemokine subfamily termed CX was generated in zebrafish. Only two human chemokines CXCL12 and CXCL14 have the orthologous chemokines in zebrafish. The diversification observed in the numbers and sequences of chemokines in the fish may reflect the adaptation of the individual species to their respective biological environment.

A putative protein structurally related to zygote arrest 1 (Zar1), Zar1-like, is encoded by a novel gene conserved in the vertebrate lineage

Identification and characterization of a bovine cDNA and the corresponding gene coding for a novel protein structurally related to Zar1, therefore called Zar1-like, are here reported for the first time. Structure of Zar1-like is similar to Zar1 gene, nevertheless they are located on distinct chromosomes. We demonstrated that the new gene as well as its genomic context are conserved along the whole vertebrate lineage. Analysis of the deduced protein primary structure showed a high conservation, among vertebrates, of the C-terminal region, where the putative presence of both zinc finger motifs and classical nuclear localization signals is also shared with Zar1. Bovine Zar1-like and the only two other available mRNA leader sequences (human and chicken) exhibit a number of upstream AUGs, suggesting that they are likely to be regulated at translational level. Expression patterns of the cattle transcripts show that Zar1-like is absent in early stages of embryo development, whereas Zar1 is expressed in matured oocytes and in in vitro produced pre-implantation embryos. In adult tissues Zar1-like transcript expression appears to be less restricted than Zar1, nevertheless, at least in bovine, both mRNAs are co-expressed in gonads, raising the question of a possible functional link.

Functional diversification of sonic hedgehog paralog enhancers identified by phylogenomic reconstruction

BACKGROUND

Cis-regulatory modules of developmental genes are targets of evolutionary changes that underlie the morphologic diversity of animals. Little is known about the 'grammar' of interactions between transcription factors and cis-regulatory modules and therefore about the molecular mechanisms that underlie changes in these modules, particularly after gene and genome duplications. We investigated the ar-C midline enhancer of sonic hedgehog (shh) orthologs and paralogs from distantly related vertebrate lineages, from fish to human, including the basal vertebrate Latimeria menadoensis.

RESULTS

We demonstrate that the sonic hedgehog a (shha) paralogs sonic hedgehog b (tiggy winkle hedgehog; shhb) genes of fishes have a modified ar-C enhancer, which specifies a diverged function at the embryonic midline. We have identified several conserved motifs that are indicative of putative transcription factor binding sites by local alignment of ar-C enhancers of numerous vertebrate sequences. To trace the evolutionary changes among paralog enhancers, phylogenomic reconstruction was carried out and lineage-specific motif changes were identified. The relation between motif composition and observed developmental differences was evaluated through transgenic functional analyses. Altering and exchanging motifs between paralog enhancers resulted in reversal of enhancer specificity in the floor plate and notochord. A model reconstructing enhancer divergence during vertebrate evolution was developed.

CONCLUSION

Our model suggests that the identified motifs of the ar-C enhancer function as binary switches that are responsible for specific activity between midline tissues, and that these motifs are adjusted during functional diversification of paralogs. The unraveled motif changes can also account for the complex interpretation of activator and repressor input signals within a single enhancer.

Genomic identification of regulatory elements by evolutionary sequence comparison and functional analysis

Despite remarkable recent advances in genomics that have enabled us to identify most of the genes in the human genome, comparable efforts to define transcriptional cis-regulatory elements that control gene expression are lagging behind. The difficulty of this task stems from two equally important problems: our knowledge of how regulatory elements are encoded in genomes remains elementary, and there is a vast genomic search space for regulatory elements, since most of mammalian genomes are noncoding. Comparative genomic approaches are having a remarkable impact on the study of transcriptional regulation in eukaryotes and currently represent the most efficient and reliable methods of predicting noncoding sequences likely to control the patterns of gene expression. By subjecting eukaryotic genomic sequences to computational comparisons and subsequent experimentation, we are inching our way toward a more comprehensive catalog of common regulatory motifs that lie behind fundamental biological processes. We are still far from comprehending how the transcriptional regulatory code is encrypted in the human genome and providing an initial global view of regulatory gene networks, but collectively, the continued development of comparative and experimental approaches will rapidly expand our knowledge of the transcriptional regulome.

Comparative genomics for detecting human disease genes

Originally, comparative genomics was geared toward defining the synteny of genes between species. As the human genome project accelerated, there was an increase in the number of tools and means to make comparisons culminating in having the genomic sequence for a large number of organisms spanning the evolutionary tree. With this level of resolution and a long history of comparative biology and comparative genetics, it is now possible to use comparative genomics to build or select better animal models and to facilitate gene discovery. Comparative genomics takes advantage of the functional genetic information from other organisms, (vertebrates and invertebrates), to apply it to the study of human physiology and disease. It allows for the identification of genes and regulatory regions, and for acquiring knowledge about gene function. In this chapter, the current state of comparative genomics and the available tools are discussed in the context of developing animal model systems that reflect the clinical picture.

Identification and characterization of a selenoprotein family containing a diselenide bond in a redox motif

Selenocysteine (Sec, U) insertion into proteins is directed by translational recoding of specific UGA codons located upstream of a stem-loop structure known as Sec insertion sequence (SECIS) element. Selenoproteins with known functions are oxidoreductases containing a single redox-active Sec in their active sites. In this work, we identified a family of selenoproteins, designated SelL, containing two Sec separated by two other residues to form a UxxU motif. SelL proteins show an unusual occurrence, being present in diverse aquatic organisms, including fish, invertebrates, and marine bacteria. Both eukaryotic and bacterial SelL genes use single SECIS elements for insertion of two Sec. In eukaryotes, the SECIS is located in the 3' UTR, whereas the bacterial SelL SECIS is within a coding region and positioned at a distance that supports the insertion of either of the two Sec or both of these residues. SelL proteins possess a thioredoxin-like fold wherein the UxxU motif corresponds to the catalytic CxxC motif in thioredoxins, suggesting a redox function of SelL proteins. Distantly related SelL-like proteins were also identified in a variety of organisms that had either one or both Sec replaced with Cys. Danio rerio SelL, transiently expressed in mammalian cells, incorporated two Sec and localized to the cytosol. In these cells, it occurred in an oxidized form and was not reducible by DTT. In a bacterial expression system, we directly demonstrated the formation of a diselenide bond between the two Sec, establishing it as the first diselenide bond found in a natural protein.

Comparative genomic and functional analyses reveal a novel cis-acting PTEN regulatory element as a highly conserved functional E-box motif deleted in Cowden syndrome

Germline mutations in PTEN, encoding a phosphatase on 10q23, cause Cowden syndrome (CS) and Bannayan-Riley-Ruvalcaba syndrome (BRRS). Approximately, 10% of CS-related PTEN mutations occur in the PTEN promoter and 11% of BRRS-related mutations include large deletions, often favoring the gene's 5' end (exon 1, promoter). In order to better understand the mechanism(s) underlying the deregulation of PTEN in these syndromes, it is important that functional cis-regulatory elements be identified. We employed a comparative genomic approach combined with molecular genetic techniques to identify a highly conserved sequence upstream of the PTEN promoter, sharing 80% sequence identity among Homo sapiens, Mus musculus and Rattus norvegicus. Within this region, we identified a canonical E-box sequence (CACGTG) located at position -2181 to -2176, approximately 800 bp upstream of the PTEN core promoter and more than 1.1 kb upstream of its minimal promoter region (located at -958 to -821). In vitro assays suggest that this motif is recognized by members of the basic region-helix-loop-helix-leucine-zipper (bHLH-LZ) transcription factor family, USF1 and USF2, and reporter assays indicate that this novel E-box is involved in mediating PTEN transcriptional activation. Four of 30 CS/CS-like patients, without previously identified PTEN mutations, were found with germline deletions of the E-box element. Of the four, three had deletions stretching to exon 1, but not 3' of it; importantly, one classic CS patient harbored a germline deletion localizing to this E-box region, further affirming the role of this element in PTEN's regulation and deregulation, and its contribution to the pathogenesis of CS.

Cloning, characterization and expression analysis of pufferfish interleukin-4 cDNA: The first evidence of Th2-type cytokine in fish

Interleukin-4 (IL-4) is one of the key cytokines in Th2 mediated immune responses, which has been shown to regulate the responses of many immune cytokines, such as interferon-gamma (IFN-gamma), interleukin-1 (IL-1) and TNF-alpha. Much work on IL-4 has been done in human and several mammal species while little in fish. In this study, we have cloned and characterized the full-length cDNA of IL-4 in Tetraodon. The Tetraodon IL-4 cDNA is 834bp in length and contains a short 5'UTR of 39bp, a 3'UTR of 375bp and an open reading frame of 420bp translating into a protein of 139aa with a predicted molecular mass of 16.131kDa. The Tetraodon IL-4-encoding gene with the same organization as the mammalians and birds consists of four exons and three introns. The encoded protein shows 11-16% identities to other homologues. RT-PCR was optimized to estimate the expression level of IL-4 in Tetraodon. The results showed that IL-4 is constitutively expressed in all selected tissues, including head kidney, spleen, liver, brain, gill, muscle and heart, although low levels were observed in head kidney, spleen, and liver. The ubiquitous expression of IL-4 is consistent with a postulated role in immune cytokines regulation. Stimulating the fish with a mixed stimulant that contained 2 microg ConA, 2 microg PHA, and 2 microg PMA, significantly up-regulated the expression of IL-4 in most tissues examined, which potentially indicated that IL-4 was involved in the immune inflammatory responses triggered by mitogens. This is the first report of cloning and characterization of IL-4 cDNA and gene in fish.

ZNF366 is an estrogen receptor corepressor that acts through CtBP and histone deacetylases

The regulation of gene expression by estrogen receptor-alpha (ERalpha) requires the coordinated and temporal recruitment of diverse sets of transcriptional co-regulator complexes, which mediate nucleosome remodelling and histone modification. Using ERalpha as bait in a yeast two-hybrid screen, we have identified a novel ERalpha-interacting protein, ZNF366, which is a potent corepressor of ERalpha activity. The interaction between ZNF366 and ERalpha has been confirmed in vitro and in vivo, and is mediated by the zinc finger domains of the two proteins. Further, we show that ZNF366 acts as a corepressor by interacting with other known ERalpha corepressors, namely RIP140 and CtBP, to inhibit expression of estrogen-responsive genes in vivo. Together, our results indicate that ZNF366 may play an important role in regulating the expression of genes in response to estrogen.

The Evf-2 noncoding RNA is transcribed from the Dlx-5/6 ultraconserved region and functions as a Dlx-2 transcriptional coactivator

The identification of ultraconserved noncoding sequences in vertebrates has been associated with developmental regulators and DNA-binding proteins. One of the first of these was identified in the intergenic region between the Dlx-5 and Dlx-6 genes, members of the Dlx/dll homeodomain-containing protein family. In previous experiments, we showed that Sonic hedgehog treatment of forebrain neural explants results in the activation of Dlx-2 and the novel noncoding RNA (ncRNA), Evf-1. In this report, we show that the Dlx-5/6 ultraconserved region is transcribed to generate an alternatively spliced form of Evf-1, the ncRNA Evf-2. Evf-2 specifically cooperates with Dlx-2 to increase the transcriptional activity of the Dlx-5/6 enhancer in a target and homeodomain-specific manner. A stable complex containing the Evf-2 ncRNA and the Dlx-2 protein forms in vivo, suggesting that the Evf-2 ncRNA activates transcriptional activity by directly influencing Dlx-2 activity. These experiments identify a novel mechanism whereby transcription is controlled by the cooperative actions of an ncRNA and a homeodomain protein. The possibility that a subset of vertebrate ultraconserved regions may function at both the DNA and RNA level to control key developmental regulators may explain why ultraconserved sequences exhibit 90% or more conservation even after 450 million years of vertebrate evolution.

Identification and characterization of DC-SCRIPT, a novel dendritic cell-expressed member of the zinc finger family of transcriptional regulators

Dendritic cells (DC) compose a heterogeneous population of cells that hold a leading role in initiating and directing immune responses. Although their function in recognizing, capturing, and presenting Ags is well defined, the molecular mechanisms that control their differentiation and immune functions are still largely unknown. In this study, we report the isolation and characterization of DC-SCRIPT, a novel protein encoded by an 8-kb mRNA that is preferentially expressed in DC. DC-SCRIPT is expressed in multiple DC subsets in vivo, including myeloid DC, plasmacytoid DC, and Langerhans cells. At the protein level, DC-SCRIPT consists of a proline-rich region, 11 C2H2-type zinc fingers, and an acidic region. Localization studies reveal that DC-SCRIPT resides in the nucleus and that nuclear localization is critically dependent on the zinc fingers. The protein displays no transcriptional activation properties according to assorted transactivation assays, but interacts with the corepressor C-terminal binding protein 1. Taken together, our results show that we have isolated a novel DC marker that could be involved in transcriptional repression. In contrast to other DC molecules, DC-SCRIPT identifies all DC subsets tested to date.

Molecular characterization of the murine homologue of the DC-derived protein DC-SCRIPT

Dendritic cell-specific transcript (DC-SCRIPT) is a putative DC zinc (Zn) finger-type transcription factor described recently in humans. Here, we illustrate that DC-SCRIPT is highly conserved in evolution and report the initial characterization of the murine ortholog of DC-SCRIPT, which is also preferentially expressed in DC as shown by real-time quantitative polymerase chain reaction, and its distribution resembles that of its human counterpart. Studies undertaken in human embryonic kidney 293 cells depict its nuclear localization and reveal that the Zn finger domain of the protein is mainly responsible for nuclear import. The human and the mouse genes are located in syntenic chromosomal regions and exhibit a similar genomic organization with numerous common transcription factor-binding sites in their promoter region, including sites for many factors implicated in haematopoiesis and DC biology, such as Gfi, GATA-1, Spi-B, and c-Rel. Taken together, these data show that DC-SCRIPT is well-conserved in evolution and that the mouse homologue is more than 80% homologous to the human protein. Therefore, mouse models can be used to elucidate the function of this novel DC marker.

Identification and characterization of an interleukin-15 homologue from Tetraodon nigroviridis

Interleukin-15 (IL-15) plays an important role in adaptive immune systems in vertebrates with similar bioactivities to interleukin-2 (IL-2). Here we report molecular cloning, sequence analysis and distribution of an IL-15 homologue from a pufferfish (Tetraodon nigroviridis). It is located within a 3,088 bp genomic fragment, transcribed into a 1,056 bp mRNA including 158 bp 5'UTR (untranslated region), 519 bp ORF (open reading frame) and 379 bp 3'UTR. T. nigroviridis IL-15 is constitutively detectable in tissues and organs selected. Levels of transcripts were observed after various stimulations. Gene organization is similar to mammals and birds, and a high degree of conservation of chromosome synteny exists between them. Systematic genomics search against Takifugu rubripes genome supports our conclusions. The T. nigroviridis IL-15 precursor with 172aa (amino acids) contains a putative 53aa signal peptide, while the mature peptide has a calculated molecular mass of 13.36 kDa and a theoretical pI of 4.67. The protein sequence shares 13.3-62.1% identity with reported IL-15s. Phylogenetic analysis grouped Tetraodon with other fish on a separated branch, excluded from mammalian and avian IL-15s. In addition, our analysis on another annotated T. nigroviridis IL-15 demonstrated that it may be a paralogue of IL-15. To differentiate it from the known IL-15s, we described it as IL-15x.

Expression and genomic organization of zonadhesin-like genes in three species of fish give insight into the evolutionary history of a mosaic protein

Background

The mosaic sperm protein zonadhesin (ZAN) has been characterized in mammals and is implicated in species-specific egg-sperm binding interactions. The genomic structure and testes-specific expression of zonadhesin is known for many mammalian species. All zonadhesin genes characterized to date consist of meprin A5 antigen receptor tyrosine phosphatase mu (MAM) domains, mucin tandem repeats, and von Willebrand (VWD) adhesion domains. Here we investigate the genomic structure and expression of zonadhesin-like genes in three species of fish.

Results

The cDNA and corresponding genomic locus of a zonadhesin-like gene (zlg) in Atlantic salmon (Salmo salar) were sequenced. Zlg is similar in adhesion domain content to mammalian zonadhesin; however, the domain order is altered. Analysis of puffer fish (Takifugu rubripes) and zebrafish (Danio rerio) sequence data identified zonadhesin (zan) genes that share the same domain order, content, and a conserved syntenic relationship with mammalian zonadhesin. A zonadhesin-like gene in D. rerio was also identified. Unlike mammalian zonadhesin, D. rerio zan and S. salar zlg were expressed in the gut and not in the testes.

Conclusion

We characterized likely orthologs of zonadhesin in both T. rubripes and D. rerio and uncovered zonadhesin-like genes in S. salar and D. rerio. Each of these genes contains MAM, mucin, and VWD domains. While these domains are associated with several proteins that show prominent gut expression, their combination is unique to zonadhesin and zonadhesin-like genes in vertebrates. The expression patterns of fish zonadhesin and zonadhesin-like genes suggest that the reproductive role of zonadhesin evolved later in the mammalian lineage.

A genomic approach to reveal novel genes associated with myotube formation in the model teleost,Takifugu rubripes

Little is known about the transcriptional networks that regulate myotube production in vertebrates. In the present study, we have used a genomic approach to discover novel genes associated with myotube formation in fast muscle of the tiger puffer fish, Takifugu rubripes. The number of fast muscle fibers per myotome increased until 1.2 kg body mass, and subsequent growth was by fiber hypertrophy alone. Forward and reverse subtracted cDNA libraries were prepared from a 180-g (myotube +) and a 3.4-kg (myotube −) fish, and 1,452 expressed sequence tags (ESTs) were obtained. After these ESTs were grouped into nonredundant clusters and housekeeping and structural genes were eliminated, 57 genes were selected and quantitative PCR was used to investigate their expression levels in different tissues from independent groups of myotube(−) and myotube(+) fish acclimated to the same environmental conditions and diet. Eleven novel genes were found to be consistently differentially expressed, but only four showed appropriate tissue-specific expression. These four genes were upregulated 5–25 times in fast muscle of myotube(−) relative to myotube(+) growth stages, while their expression remained unchanged in the other tissues studied. The novel genes identified, which are also present in other vertebrate genomes, may play a role in inhibiting myotube formation in vertebrate muscle.

Molecular cloning and characterization of ZFF29: a protein containing a unique Cys2His2 zinc-finger motif

We have cloned a gene, ZFF29 (zinc-finger protein of human fetal liver erythroid cells 29), from human fetal liver erythroid cells. Two types of mature mRNA were identified and designated ZFF29a and ZFF29b. In human genome the ZFF29 gene is on chromosome 9q, and the two forms are splice variants. There is a unique transcription start site, which predicts major mRNAs composed of 2485 bases for ZFF29a and 1801 bases for ZFF29b. The anticipated mRNAs were demonstrated in K562 cells, but not in any adult human tissues examined by Northern blotting. In the mouse, reverse transcription-PCR revealed that the ZFF29 mRNA is present in adult bone marrow and ovary at a higher level than in any other tissues examined. These findings suggest that ZFF29 proteins are expressed in embryonic/fetal erythroid tissues. The deduced polypeptide chains of ZFF29a and ZFF29b are composed of 306 and 350 amino acids respectively. A unique zinc-finger motif composed of two contiguous Cys(2)His(2)-type fingers is common to both forms of ZFF29. They are nuclear proteins and ZFF29b, but not ZFF29a, is an activator of erythroid gene promoters.

Mulan: multiple-sequence local alignment and visualization for studying function and evolution

Multiple-sequence alignment analysis is a powerful approach for understanding phylogenetic relationships, annotating genes, and detecting functional regulatory elements. With a growing number of partly or fully sequenced vertebrate genomes, effective tools for performing multiple comparisons are required to accurately and efficiently assist biological discoveries. Here we introduce Mulan (http://mulan.dcode.org/), a novel method and a network server for comparing multiple draft and finished-quality sequences to identify functional elements conserved over evolutionary time. Mulan brings together several novel algorithms: the TBA multi-aligner program for rapid identification of local sequence conservation, and the multiTF program for detecting evolutionarily conserved transcription factor binding sites in multiple alignments. In addition, Mulan supports two-way communication with the GALA database; alignments of multiple species dynamically generated in GALA can be viewed in Mulan, and conserved transcription factor binding sites identified with Mulan/multiTF can be integrated and overlaid with extensive genome annotation data using GALA. Local multiple alignments computed by Mulan ensure reliable representation of short- and large-scale genomic rearrangements in distant organisms. Mulan allows for interactive modification of critical conservation parameters to differentially predict conserved regions in comparisons of both closely and distantly related species. We illustrate the uses and applications of the Mulan tool through multispecies comparisons of the GATA3 gene locus and the identification of elements that are conserved in a different way in avians than in other genomes, allowing speculation on the evolution of birds. Source code for the aligners and the aligner-evaluation software can be freely downloaded from http://www.bx.psu.edu/miller_lab/.

Prediction of nuclear hormone receptor response elements

The nuclear receptor (NR) class of transcription factors controls critical regulatory events in key developmental processes, homeostasis maintenance, and medically important diseases and conditions. Identification of the members of a regulon controlled by a NR could provide an accelerated understanding of development and disease. New bioinformatics methods for the analysis of regulatory sequences are required to address the complex properties associated with known regulatory elements targeted by the receptors because the standard methods for binding site prediction fail to reflect the diverse target site configurations. We have constructed a flexible Hidden Markov Model framework capable of predicting NHR binding sites. The model allows for variable spacing and orientation of half-sites. In a genome-scale analysis enabled by the model, we show that NRs in Fugu rubripes have a significant cross-regulatory potential. The model is implemented in a web interface, freely available for academic researchers, available at http://mordor.cgb.ki.se/NHR-scan.

Evolution of vertebrate genes related to prion and Shadoo proteins--clues from comparative genomic analysis

Recent findings of new genes in fish related to the prion protein (PrP) gene PRNP, including our recent report of SPRN coding for Shadoo (Sho) protein found also in mammals, raise issues of their function and evolution. Here we report additional novel fish genes found in public databases, including a duplicated SPRN gene, SPRNB, in Fugu, Tetraodon, carp, and zebrafish encoding the Sho2 protein, and we use comparative genomic analysis to analyze the evolutionary relationships and to infer evolutionary trajectories of the complete data set. Phylogenetic footprinting performed on aligned human, mouse, and Fugu SPRN genes to define candidate regulatory promoter regions, detected 16 conserved motifs, three of which are known transcription factor-binding sites for a receptor and transcription factors specific to or associated with expression in brain. This result and other homology-based (VISTA global genomic alignment; protein sequence alignment and phylogenetics) and context-dependent (genomic context; relative gene order and orientation) criteria indicate fish and mammalian SPRN genes are orthologous and suggest a strongly conserved basic function in brain. Whereas tetrapod PRNPs share context with the analogous stPrP-2-coding gene in fish, their sequences are diverged, suggesting that the tetrapod and fish genes are likely to have significantly different functions. Phylogenetic analysis predicts the SPRN/SPRNB duplication occurred before divergence of fish from tetrapods, whereas that of stPrP-1 and stPrP-2 occurred in fish. Whereas Sho appears to have a conserved function in vertebrate brain, PrP seems to have an adaptive role fine-tuned in a lineage-specific fashion. An evolutionary model consistent with our findings and literature knowledge is proposed that has an ancestral prevertebrate SPRN-like gene leading to all vertebrate PrP-related and Sho-related genes. This provides a new framework for exploring the evolution of this unusual family of proteins and for searching for members in other fish branches and intermediate vertebrate groups.

eShadow: a tool for comparing closely related sequences

Primate sequence comparisons are difficult to interpret due to the high degree of sequence similarity shared between such closely related species. Recently, a novel method, phylogenetic shadowing, has been pioneered for predicting functional elements in the human genome through the analysis of multiple primate sequence alignments. We have expanded this theoretical approach to create a computational tool, eShadow, for the identification of elements under selective pressure in multiple sequence alignments of closely related genomes, such as in comparisons of human-to-primate or mouse-to-rat DNA. This tool integrates two different statistical methods and allows for the dynamic visualization of the resulting conservation profile. eShadow also includes a versatile optimization module capable of training the underlying Hidden Markov Model to differentially predict functional sequences. This module grants the tool high flexibility in the analysis of multiple sequence alignments and in comparing sequences with different divergence rates. Here, we describe the eShadow comparative tool and its potential uses for analyzing both multiple nucleotide and protein alignments to predict putative functional elements.

Comparative genomic analysis as a tool for biological discovery

The recent completion of the human genome sequence has enabled the identification of a large fraction of our gene catalogue and their physical chromosomal position. However, current efforts lag at defining the cis-regulatory sequences that control the spatial and temporal patterns of each gene's expression. This task remains difficult due to our lack of knowledge of the vocabulary controlling gene regulation and the vast genomic search space, with greater than 95% of our genome being noncoding. Recent comparative genomic-based strategies are beginning to aid in the identification of functional sequences based on their high levels of evolutionary conservation. This has proven successful for comparisons between closely related species such as human-primate or human-mouse, but also holds true for distant evolutionary comparisons, such as human-fish or human-bird. In this review we provide support for the utility of cross-species sequence comparisons by illustrating several applications of this strategy, including the identification of new genes and functional non-coding sequences. We also discuss emerging concepts as this field matures, such as how to properly select which species for comparison, which may differ significantly between independent studies.

Analysis of multiple genomic sequence alignments: a web resource, online tools, and lessons learned from analysis of mammalian SCL loci

Comparative analysis of genomic sequences is becoming a standard technique for studying gene regulation. However, only a limited number of tools are currently available for the analysis of multiple genomic sequences. An extensive data set for the testing and training of such tools is provided by the SCL gene locus. Here we have expanded the data set to eight vertebrate species by sequencing the dog SCL locus and by annotating the dog and rat SCL loci. To provide a resource for the bioinformatics community, all SCL sequences and functional annotations, comprising a collation of the extensive experimental evidence pertaining to SCL regulation, have been made available via a Web server. A Web interface to new tools specifically designed for the display and analysis of multiple sequence alignments was also implemented. The unique SCL data set and new sequence comparison tools allowed us to perform a rigorous examination of the true benefits of multiple sequence comparisons. We demonstrate that multiple sequence alignments are, overall, superior to pairwise alignments for identification of mammalian regulatory regions. In the search for individual transcription factor binding sites, multiple alignments markedly increase the signal-to-noise ratio compared to pairwise alignments.

Identification of 5'-upstream region of pufferfish ribosomal protein L29 gene as a strong constitutive promoter to drive GFP expression in zebrafish

The genomic structure of Tetraodon fluviatilis L29 gene was determined and its promoter activity was analyzed in COS-1 cells and zebrafish embryos. The TfL29 gene comprises four exons and three introns, spanning approximately 1.7kb. The 5(')-upstream 2.2-kb of the first exon contains 10 E-boxes and many putative binding motifs for transcription factors GATA-1, AML-1a, c-Myb, Oct-1, CdxA, and NRF-2. Promoter activity assay showed that the distal 2.2-kb fragment not only had high luciferase activity in COS-1 cells, but also strong and ubiquitous GFP expression in a variety of tissues in zebrafish embryos. On the other hand, there are no TATA or CAAT boxes within a 300-bp region upstream from the transcription initiation site. Although this region has high luciferase activity in COS-1 cells, it is not sufficient to drive GFP expression in zebrafish embryos. In this proximal 300-bp region, there are two E-boxes, two CdxA sites, and one NRF-2 site that is immediately downstream of the transcription start site.

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.

Whole-genome discovery of transcription factor binding sites by network-level conservation

Comprehensive identification of DNA cis-regulatory elements is crucial for a predictive understanding of transcriptional network dynamics. Strong evidence suggests that these DNA sequence motifs are highly conserved between related species, reflecting strong selection on the network of regulatory interactions that underlie common cellular behavior. Here, we exploit a systems-level aspect of this conservation-the network-level topology of these interactions-to map transcription factor (TF) binding sites on a genomic scale. Using network-level conservation as a constraint, our algorithm finds 71% of known TF binding sites in the yeast Saccharomyces cerevisiae, using only 12% of the sequence of a phylogenetic neighbor. Most of the novel predicted motifs show strong features of known TF binding sites, such as functional category and/or expression profile coherence of their corresponding genes. Network-level conservation should provide a powerful constraint for the systematic mapping of TF binding sites in the larger genomes of higher eukaryotes.

Compact intergenic regions of the pufferfish genome facilitate isolation of gene promoters: characterization ofFugu3′-phosphoadenosine 5′-phosphosulfate synthase 2 (fPapss2) gene promoter function in transgenicXenopus

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.

Shadoo, a new protein highly conserved from fish to mammals and with similarity to prion protein

We report evidence from cDNA isolation and expression analysis as well as analyses of genome, expressed sequence tag (EST), cDNA and expression databases for a new gene named SPRN (shadow of prion protein). SPRN comprises two exons, with the open reading frame (ORF) contained within exon 2, and codes for a protein of 130-150 amino acids named Shadoo (Japanese shadow), predicted to be extracellular and GPI-anchored. The SPRN gene was found in fish (zebrafish, Fugu) and mammals (mouse, rat, human). Conservation of order and transcription orientation of two proximal genes between fishes and mammals strongly indicates gene orthology. Sequence comparison shows: a highly conserved N-terminal signal sequence; Arg-rich basic region containing up to six tetrarepeats of consensus XXRG (where X is G, A or S); a hydrophobic region of 20 residues with strong homology to PrP; a less conserved C-terminal domain containing a conserved glycosylation motif; and a C-terminal peptide predicted to be a signal sequence for glycophosphotidylinositol (GPI)-anchor attachment. Fish Shadoos (Sho) show well conserved sequences (identity 54%) over 106 amino acids (zebrafish length), and conservation among the mammalian sequences is very high (identity 81-96%). The fish and mammalian sequences are also well conserved, particularly for zebrafish, to beyond the end of the hydrophobic sequence (identity 41-53%, 78 amino acids, zebrafish length). The overall structure appears closely related to prion proteins (PrPs), although the C-terminal domains of Shos are quite different from those of PrPs, for which conformational changes in mammals are implicated in disease. The structural similarity is particularly interesting given recent reports of three new genes with similarities to PrPs found in Fugu (PrP-like, PrP-461/stPrP-1 and stPrP-2) and other fish, but for which direct evolution to higher vertebrate PrPs is unlikely and for which no other mammalian homologues have been found. Database information indicates expression of SPRN in embryo, brain and retina of mouse and rat, hippocampus of human, and in embryo and retina of zebrafish, and we directly confirmed a strikingly specific expression of the mammalian (human, mouse, rat) transcripts in whole brain. This result together with some common structural features led to the suggestive hypothesis of a possible functional link between mammalian PrP and Sho proteins.

Neurone-specific expression and regulation of the pufferfish isotocin and vasotocin genes in transgenic mice

We used comparative genetics to investigate the location, structure and evolution of the oxytocin and vasopressin gene regulatory regions. The pufferfish, Fugu rubripes, is an attractive vertebrate model for comparison because of its maximal evolutionary distance from mammals and short intergenic regions. To determine whether regulatory DNA is conserved between oxytocin and vasopressin, and their Fugu homologs, isotocin and vasotocin, we generated transgenic mice bearing overlapping Fugu cosmids that contained the isotocin and/or vasotocin genes as well as short isotocin (5 kb) and vasotocin (9 kb) constructs. Our study shows that the Fugu isotocin and vasotocin genes express specifically in the mouse oxytocinergic and vasopressinergic neurones, respectively, and that the cis-regulatory elements which mediate neurone-specific expression are located within the short transgene constructs tested. Thus, the neurone-specific expression of the oxytocin and vasopressin gene families, and the mechanisms mediating the cell-specificity, evolved before the divergence of the fish and mammalian lineages. Salt-loading of transgenic mice induced an increase in abundance of isotocin, but not vasotocin mRNA in the cognate neurones. It appears that either the vasotocin gene does not respond to osmotic perturbations or the vasotocin transgene construct tested lacks osmotic response elements. Comparisons of homologous flanking sequences of the Fugu and mouse genes identified several short matching sequences, which are candidate regulatory elements.

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.

Insights from human/mouse genome comparisons

Large-scale public genomic sequencing efforts have provided a wealth of vertebrate sequence data poised to provide insights into mammalian biology. These include deep genomic sequence coverage of human, mouse, rat, zebrafish, and two pufferfish ( Fugu rubripes and Tetraodon nigroviridis) (Aparicio et al. 2002; Lander et al. 2001; Venter et al. 2001; Waterston et al. 2002). In addition, a high-priority has been placed on determining the genomic sequence of chimpanzee, dog, cow, frog, and chicken (Boguski 2002). While only recently available, whole genome sequence data have provided the unique opportunity to globally compare complete genome contents. Furthermore, the shared evolutionary ancestry of vertebrate species has allowed the development of comparative genomic approaches to identify ancient conserved sequences with functionality. Accordingly, this review focuses on the initial comparison of available mammalian genomes and describes various insights derived from such analysis.

Comparative genomics: genome-wide analysis in metazoan eukaryotes

The increasing number of complete and nearly complete metazoan genome sequences provides a significant amount of material for large-scale comparative genomic analysis. Finding new effective methods to analyse such enormous datasets has been the object of intense research. Three main areas in comparative genomics have recently shown important developments: whole-genome alignment, gene prediction and regulatory-region prediction. Each of these areas improves the methods of deciphering long genomic sequences and uncovering what lies hidden in them.

Comparative DNA sequence analysis of mouse and human CC chemokine gene clusters

The CC chemokines are a closely related subfamily of the chemokine superfamily. Most of the CC chemokine genes form a cluster on chromosome 11 in mice and chromosome 17 in humans. To date, 11 and 16 functional genes have been localized within the mouse and human clusters, respectively. Notably, some of the genes within these clusters appear to have no counterparts between the two species, and the orthologous relationships of some of the genes are difficult to establish solely on the basis of amino acid similarity. In this study, we have taken a comparative genomic approach to reveal some of the features that may be involved in the dynamic evolution of these gene clusters. We sequenced a 122-kb region containing five chemokine genes of the mouse CC cluster. This mouse sequence was combined with those determined by the Mouse Genome Sequencing Project, and the entire sequence of the mouse CC cluster was compared with that of the corresponding cluster in the human genome by percent identity plot and dot-plot analyses. Although no additional chemokine genes have been found in these clusters, our analysis has revealed that numerous gene rearrangements have occurred even after the diversification of rodents and primates, resulting in several species-specific chemokine genes and pseudogenes. In addition, phylogenetic analysis and comparison of the genomic sequences unambiguously identified the orthologous relationships of some of the chemokine genes in the mouse and human CC gene clusters.