Genome evolution and biodiversity in teleost fish

Genome-wide analysis and expression profiling of the small heat shock proteins in zebrafish

Small Heat Shock Proteins (sHSPs) have important roles in preventing disease and promoting resistance to environmental stressors. Mutations in any one of a number of sHSPs, including HSP27 (HSPB1), HSP22 (HSPB8), alphaA-crystallin (HSPB4), or alphaB-crystallin (HSPB5) can result in neuronal degeneration, myopathy, and/or cataract in humans. Ten sHSPs are known in humans, and thirteen have been identified in teleost fish. Here we report the identification of thirteen zebrafish sHSPs. Using a combination of phylogenetic analysis and analysis of synteny, we have determined that ten are likely orthologs of human sHSPs. We have used quantitative RT-PCR to determine the relative expression levels of all thirteen sHSPs during development and in response to heat shock. Our findings indicate that most of the zebrafish sHSPs are expressed during development, and five of these genes are transcriptionally upregulated by heat shock at one or more stages of development.

Early diversification of the TNF superfamily in teleosts: genomic characterization and expression analysis

The TNF superfamily (TNFSF) of proteins are cytokines involved in diverse immunological and developmental pathways. Little is known about their evolution or expression in lower vertebrate species. Bioinformatic searches of Zebrafish, Tetraodon, and Fugu genome and other teleost expressed sequence tag databases identified 44 novel gene sequences containing a TNF homology domain. This work reveals the following: 1) teleosts possess orthologs of BAFF, APRIL, EDA, TWEAK, 4-1BBL, Fas ligand, LIGHT, CD40L, RANKL, and possibly TL1A; 2) the BAFF-APRIL subfamily is enriched by a third member, BALM, unique to fish; 3) orthologs of lymphotoxins alpha and beta were not clearly identified in teleosts and are substituted by a related ligand, TNF-New; 4) as many as four TRAIL-like genes are present in teleosts, as compared with only one in mammals; and 5) T cell activation ligands OX40L, CD27L, CD30L, and GITRL were not identified in any fish species. Finally, we characterize mRNA expression of TNFSF members CD40L, LIGHT, BALM, APRIL, Fas ligand, RANKL, TRAIL-like, and TNF-New in rainbow trout, Oncorhynchus mykiss, immune and nonimmune tissues. In conclusion, we identified a total of 14 distinct TNFSF members in fishes, indicating expansion of this superfamily before the divergence of bony fish and tetrapods, approximately 360-450 million years ago. Based on these findings, we extend a model of TNFSF evolution and the co-emergence of the vertebrate adaptive immune system.

GnRH and GnRH receptors in metazoa: a historical, comparative, and evolutive perspective

About 50years after Harris's first demonstration of its existence, GnRH has strongly stimulated the interest and imagination of scientists, resulting in a high number of studies in an increasing number of species. For the endocrinologist, GnRH, via its actions on the synthesis and release of pituitary gonadotrophins, is first an essential hormone for the initiation and maintenance of the reproductive axis, but recent data suggest that GnRH emerged in animals lacking a pituitary. In this context, this review intends to explore the current status of knowledge on GnRH and GnRH receptors in metazoa in order to see if it is possible to draw an evolutive scenario according to which GnRH actions progressively evolved from the control of simple basic functions in early metazoa to an indirect mean of controlling gonadal activity in vertebrates through a sophisticated network of finely tuned neurons developing in a rather fascinating way. This review also intends to provide an evolutive scenario based on the recent advances of whole genome sequencing possibly explaining the number of GnRH and GnRH receptor variants according to the 2R and 3R theories accompanied by gene losses.

Emerging roles for zic genes in early development

Members of the Zic family of zinc finger transcription factors play critical roles in a variety of developmental processes. They are involved in development of neural tissues and the neural crest, in left-right axis patterning, in somite development, and in formation of the cerebellum. In addition to their roles in cell-fate specification, zic genes also promote cell proliferation. Further, they are expressed in postmitotic cells of the cerebellum and in retinal ganglion cells. Efforts to determine the role of individual zic genes within an array of developmental and cellular processes are complicated by overlapping patterns of zic gene expression and strong sequence conservation within this gene family. Nevertheless, substantial progress has been made. This review summarizes our knowledge of the molecular events that govern the activities of zic family members, including emerging relationships between upstream signaling pathways and zic genes. In addition, advancements in our understanding of the molecular events downstream of Zic transcription factors are reviewed. Despite significant progress, however, much remains to be learned regarding the mechanisms through which zic genes exert their function in a variety of different contexts.

Comparative sequence analysis of the complete set of 40S ribosomal proteins in the Senegalese sole (Solea senegalensis Kaup) and Atlantic halibut (Hippoglossus hippoglossus L.) (Teleostei: Pleuronectiformes): phylogeny and tissue- and development-specific expression

Background

Ribosomal proteins (RPs) are key components of ribosomes, the cellular organelle responsible for protein biosynthesis in cells. Their levels can vary as a function of organism growth and development; however, some RPs have been associated with other cellular processes or extraribosomal functions. Their high representation in cDNA libraries has resulted in the increase of RP sequences available from different organisms and their proposal as appropriate molecular markers for phylogenetic analysis.

Results

The development of large-scale genomics of Senegalese sole (Solea senegalensis) and Atlantic halibut (Hippoglossus hippoglossus), two commercially important flatfish species, has made possible the identification and systematic analysis of the complete set of RP sequences for the small (40S) ribosome subunit. Amino acid sequence comparisons showed a high similarity both between these two flatfish species and with respect to other fish and human. EST analysis revealed the existence of two and four RPS27 genes in Senegalese sole and Atlantic halibut, respectively. Phylogenetic analysis clustered RPS27 in two separate clades with their fish and mammalian counterparts. Steady-state transcript levels for eight RPs (RPS2, RPS3a, RPS15, RPS27-1, RPS27-2, RPS27a, RPS28, and RPS29) in sole were quantitated during larval development and in tissues, using a real-time PCR approach. All eight RPs exhibited different expression patterns in tissues with the lowest levels in brain. On the contrary, RP transcripts increased co-ordinately after first larval feeding reducing progressively during the metamorphic process.

Conclusion

The genomic resources and knowledge developed in this survey will provide new insights into the evolution of Pleuronectiformes. Expression data will contribute to a better understanding of RP functions in fish, especially the mechanisms that govern growth and development in larvae, with implications in aquaculture.

Evolution of pigment synthesis pathways by gene and genome duplication in fish

BACKGROUND

Coloration and color patterning belong to the most diverse phenotypic traits in animals. Particularly, teleost fishes possess more pigment cell types than any other group of vertebrates. As the result of an ancient fish-specific genome duplication (FSGD), teleost genomes might contain more copies of genes involved in pigment cell development than tetrapods. No systematic genomic inventory allowing to test this hypothesis has been drawn up so far for pigmentation genes in fish, and almost nothing is known about the evolution of these genes in different fish lineages.

RESULTS

Using a comparative genomic approach including phylogenetic reconstructions and synteny analyses, we have studied two major pigment synthesis pathways in teleost fish, the melanin and the pteridine pathways, with respect to different types of gene duplication. Genes encoding three of the four enzymes involved in the synthesis of melanin from tyrosine have been retained as duplicates after the FSGD. In the pteridine pathway, two cases of duplicated genes originating from the FSGD as well as several lineage-specific gene duplications were observed. In both pathways, genes encoding the rate-limiting enzymes, tyrosinase and GTP-cyclohydrolase I (GchI), have additional paralogs in teleosts compared to tetrapods, which have been generated by different modes of duplication. We have also observed a previously unrecognized diversity of gchI genes in vertebrates. In addition, we have found evidence for divergent resolution of duplicated pigmentation genes, i.e., differential gene loss in divergent teleost lineages, particularly in the tyrosinase gene family.

CONCLUSION

Mainly due to the FSGD, teleost fishes apparently have a greater repertoire of pigment synthesis genes than any other vertebrate group. Our results support an important role of the FSGD and other types of duplication in the evolution of pigmentation in fish.

Evolution of pigment synthesis pathways by gene and genome duplication in fish

BACKGROUND

Coloration and color patterning belong to the most diverse phenotypic traits in animals. Particularly, teleost fishes possess more pigment cell types than any other group of vertebrates. As the result of an ancient fish-specific genome duplication (FSGD), teleost genomes might contain more copies of genes involved in pigment cell development than tetrapods. No systematic genomic inventory allowing to test this hypothesis has been drawn up so far for pigmentation genes in fish, and almost nothing is known about the evolution of these genes in different fish lineages.

RESULTS

Using a comparative genomic approach including phylogenetic reconstructions and synteny analyses, we have studied two major pigment synthesis pathways in teleost fish, the melanin and the pteridine pathways, with respect to different types of gene duplication. Genes encoding three of the four enzymes involved in the synthesis of melanin from tyrosine have been retained as duplicates after the FSGD. In the pteridine pathway, two cases of duplicated genes originating from the FSGD as well as several lineage-specific gene duplications were observed. In both pathways, genes encoding the rate-limiting enzymes, tyrosinase and GTP-cyclohydrolase I (GchI), have additional paralogs in teleosts compared to tetrapods, which have been generated by different modes of duplication. We have also observed a previously unrecognized diversity of gchI genes in vertebrates. In addition, we have found evidence for divergent resolution of duplicated pigmentation genes, i.e., differential gene loss in divergent teleost lineages, particularly in the tyrosinase gene family.

CONCLUSION

Mainly due to the FSGD, teleost fishes apparently have a greater repertoire of pigment synthesis genes than any other vertebrate group. Our results support an important role of the FSGD and other types of duplication in the evolution of pigmentation in fish.

The evolution of the vertebrate metzincins; insights from Ciona intestinalis and Danio rerio

BACKGROUND

The metzincins are a large gene superfamily of proteases characterized by the presence of a zinc protease domain, and include the ADAM, ADAMTS, BMP1/TLL, meprin and MMP genes. Metzincins are involved in the proteolysis of a wide variety of proteins, including those of the extracellular matrix. The metzincin gene superfamily comprises eighty proteins in the human genome and ninety-three in the mouse. When and how the level of complexity apparent in the vertebrate metzincin gene superfamily arose has not been determined in detail. Here we present a comprehensive analysis of vertebrate metzincins using genes from both Ciona intestinalis and Danio rerio to provide new insights into the complex evolution of this gene superfamily.

RESULTS

We have identified 19 metzincin genes in the ciona genome and 83 in the zebrafish genome. Phylogenetic analyses reveal that the expansion of the metzincin gene superfamily in vertebrates has occurred predominantly by the simple duplication of pre-existing genes rather than by the appearance and subsequent expansion of new metzincin subtypes (the only example of which is the meprin gene family). Despite the number of zebrafish metzincin genes being relatively similar to that of tetrapods (e.g. man and mouse), the pattern of gene retention and loss within these lineages is markedly different. In addition, we have studied the evolution of the related TIMP gene family and identify a single ciona and four zebrafish TIMP genes.

CONCLUSION

The complexity seen in the vertebrate metzincin gene families was mainly acquired during vertebrate evolution. The metzincin gene repertoire in protostomes and invertebrate deuterostomes has remained relatively stable. The expanded metzincin gene repertoire of extant tetrapods, such as man, has resulted largely from duplication events associated with early vertebrate evolution, prior to the sarcopterygian-actinopterygian split. The teleost repertoire of metzincin genes in part parallels that of tetrapods but has been significantly modified, perhaps as a consequence of a teleost-specific duplication event.

Expression analysis of two Eomesodermin homologues in zebrafish lymphoid tissues and cells

Eomesodermin (Eomes) is a T-box transcription factor that is involved in mesoderm formation in most vertebrates. Eomes is also expressed in CD8+ T cells and NK cells. No information is available on the role of Eomes in the immune system of lower vertebrates to date, although developmental studies on Eomes (Eomes1) have been performed in zebrafish. Here we report the identification of a second Eomes (Eomes2) in zebrafish and compare expression of the two Eomes genes in the immune system. Zebrafish Eomes1 and Eomes2, composed of 661 and 534 amino acids, respectively, share 49.3% amino acid identity in their coding regions and 88.7% amino acid identity in their T-box regions. Conserved synteny between regions of the human and zebrafish genomes, gene organization and phylogenetic analysis all indicate that the zebrafish Eomes2 gene is a homologue of mammalian Eomes, as previously found for zebrafish Eomes1. Eomes1 mRNA was found to be expressed in the gonad, body kidney, spleen and gill, while Eomes2 mRNA was not detected in any of these tissues. However, strong expression of both Eomes mRNAs was detected in the leukocytes from the spleen, followed by those from body kidney and peripheral blood, with expression of Eomes1 always stronger than that of Eomes2. RT-PCR analysis of body kidney cells sorted by FACS revealed that Eomes1 was expressed strongly in lymphocytes, weakly in blast cells, and was not expressed in granulocytes, while Eomes2 was expressed weakly in lymphocytes. These results suggest that both Eomes genes are involved in the zebrafish immune response, particularly in lymphocyte function as has been found in mammals.

Characterization of twin toll-like receptors from rainbow trout (Oncorhynchus mykiss): evolutionary relationship and induced expression by Aeromonas salmonicida salmonicida

Structure and function of factors contributing to the innate immune system of lower vertebrates, including fish are only sparsely characterized. We retrieved with RT-PCR cDNA copies of two closely related Toll-like receptors (TLR) from liver RNA of the rainbow trout (Oncorhynchus mykiss). The cDNA sequences are homologous to 95.6%. The phylogenetic analysis of their deduced amino acid sequences places these twin factors closely to other known TLRs from fish. The twin factors are equally expressed in all tissues analysed, most abundantly in spleen and head kidney and lowest in adipose tissue. Formalin-inactivated Aeromonas salmonicida pathogens induce their expression up to eight-fold in vitro in peripheral blood lymphocytes and in tissues from spleen and head kidney. Our sequence information will be useful to establish expression constructs for these factors necessary to analyse the pathogen specific signal transduction activating the innate immune defence in fish.

Transcriptional regulation of arylalkylamine-N-acetyltransferase-2 gene in the pineal gland of the gilthead seabream

Pineal serotonin-N-acetyltransferase (arylalkylamine-N-acetyltransferase; AANAT) is considered the key enzyme in the generation of circulating melatonin rhythms; the rate of melatonin production is determined by AANAT activity. In all the examined species, AANAT activity is regulated at the post-translational level and, to a variable degree, also at the transcriptional level. Here, the transcriptional regulation of pineal aanat (aanat2) of the gilthead seabream (Sparus aurata) was investigated. Real-time polymerase chain reaction quantification of aanat2 mRNA levels in the pineal gland collected throughout the 24-h cycle revealed a rhythmic expression pattern. In cultured pineal glands, the amplitude was reduced, but the daily rhythmic expression pattern was maintained under constant illumination, indicating a circadian clock-controlled regulation of seabream aanat2. DNA constructs were prepared in which green fluorescent protein was driven by the aanat2 promoters of seabream and Northern pike. In vivo transient expression analyses in zebrafish embryos indicated that these promoters contain the necessary elements to drive enhanced expression in the pineal gland. In the light-entrainable clock-containing PAC-2 zebrafish cell line, a stably transfected seabream aanat2 promoter-luciferase DNA construct exhibited a clock-controlled circadian rhythm of luciferase activity, characteristic for an E-box-driven expression. In NIH-3T3 cells, the seabream aanat2 promoter was activated by a synergistic action of BMAL/CLOCK and orthodenticle homeobox 5 (OTX5). Promoter sequence analyses revealed the presence of the photoreceptor conserved element and an extended E-box (i.e. the binding sites for BMAL/CLOCK and OTX5 that have been previously associated with pineal-specific and rhythmic gene expression). These results suggest that seabream aanat2 is a clock-controlled gene that is regulated by conserved mechanisms.

The Evolution of Sexually Selected Traits and Antagonistic Androgen Expression in Actinopterygiian Fishes

Many sexually selected traits in male fishes are controlled by testosterone. Directional selection for male ornaments could theoretically increase male testosterone levels over evolutionary timescales, and when genetically correlated, female testosterone levels as well. Because of the negative fitness consequences of high testosterone, it is plausible that female choice for sexually selected traits in males results in decreased female reproductive fitness. I used comparative analysis to examine the association between male peak testosterone expression and sexually selected ornaments. I also tested for genetic correlation between male and female androgen levels. The presence of sexually selected traits in males was significantly correlated with increased peak androgen levels in males as well as females, and female testosterone levels were significantly correlated with male peak testosterone titers, although the slope was only marginally <1. This suggests that selection to decouple high male and female testosterone levels is either weak or otherwise ineffective.

Globalization and fruitfly invasion and expansion: the medfly paradigm

The phytophagous insects of the Tephritidae family commonly referred to as "true fruit flies" offer different case histories of successful invasions. Mankind has played an important role in altering the distributions of some of the more polyphagous and oligophagous species. However, the question arises why only a few species have become major invaders. The understanding of traits underlying adaptation in different environments is a major topic in invasion biology. Being generalists or specialists, along the K-r gradient of the growth curve, make a difference in term of food resources exploitation and interspecies competition and displacement. The species of the genus Ceratitis are good examples of r-strategists. The genetic and biological data of the most notorious Ceratitis species, the Mediterranean fruit fly Ceratitis capitata (medfly), are reviewed to investigate the traits and behaviours that make the medfly an important invader. It can be learnt from medfly, that invasions in a modern global trade network tend to be due to multiple introductions. This fact allows a maintenance or enhancement of genetic variability in the adventive populations, which in turn increases their potential invasiveness. Our current knowledge of the medfly genome opens the way for future studies on functional genomics.

Phylogenetic conservation of chromosome numbers in Actinopterygiian fishes

The genomes of ray-finned fishes (Actinopterygii) are well known for their evolutionary dynamism as reflected by drastic alterations in DNA content often via regional and whole-genome duplications, differential patterns of gene silencing or loss, shifts in the insertion-to-deletion ratios of genomic segments, and major re-patternings of chromosomes via non-homologous recombination. In sharp contrast, chromosome numbers in somatic karyotypes have been highly conserved over vast evolutionary timescales - a histogram of available counts is strongly leptokurtic with more than 50% of surveyed species displaying either 48 or 50 chromosomes. Here we employ comparative phylogenetic analyses to examine the evolutionary history of alterations in fish chromosome numbers. The most parsimonious ancestral state for major actinopterygiian clades is 48 chromosomes. When interpreted in a phylogenetic context, chromosome numbers evidence many recent instances of polyploidization in various lineages but there is no clear indication of a singular polyploidization event that has been hypothesized to have immediately preceded the teleost radiation. After factoring out evident polyploidizations, a correlation between chromosome numbers and genome sizes across the Actinopterygii is marginally statistically significant (p = 0.012) but exceedingly weak (R (2) = 0.0096). Overall, our phylogenetic analysis indicates a mosaic evolutionary pattern in which the forces that govern labile features of fish genomes must operate largely independently of those that operate to conserve chromosome numbers.

AtEXO70A1, a member of a family of putative exocyst subunits specifically expanded in land plants, is important for polar growth and plant development

The exocyst is a hetero-oligomeric protein complex involved in exocytosis and has been extensively studied in yeast and animal cells. Evidence is now accumulating that the exocyst is also present in plants. Bioinformatic analysis of genes encoding plant homologs of the exocyst subunit, Exo70, revealed that three Exo70 subgroups are evolutionarily conserved among angiosperms, lycophytes and mosses. Arabidopsis and rice contain 22 and approximately 39 EXO70 genes, respectively, which can be classified into nine clusters considered to be ancient in angiosperms (one has been lost in Arabidopsis). We characterized two independent T-DNA insertional mutants of the AtEXO70A1 gene (exo70A1-1 and exo70A1-2). Heterozygous EXO70A1/exo70A1 plants appear to be normal and segregate in a 1:2:1 ratio, suggesting that neither male nor female gametophytes are affected by the EXO70A1 disruption. However, both exo70A1-1 and exo70A1-2 homozygotes exhibit an array of phenotypic defects. The polar growth of root hairs and stigmatic papillae is disturbed. Organs are generally smaller, plants show a loss of apical dominance and indeterminate growth where instead of floral meristems new lateral inflorescences are initiated in a reiterative manner. Both exo70A1 mutants have dramatically reduced fertility. These results suggest that the putative exocyst subunit EXO70A1 is involved in cell and organ morphogenesis.

Coding limits on the number of transcription factors

Background

Transcription factor proteins bind specific DNA sequences to control the expression of genes. They contain DNA binding domains which belong to several super-families, each with a specific mechanism of DNA binding. The total number of transcription factors encoded in a genome increases with the number of genes in the genome. Here, we examined the number of transcription factors from each super-family in diverse organisms.

Results

We find that the number of transcription factors from most super-families appears to be bounded. For example, the number of winged helix factors does not generally exceed 300, even in very large genomes. The magnitude of the maximal number of transcription factors from each super-family seems to correlate with the number of DNA bases effectively recognized by the binding mechanism of that super-family. Coding theory predicts that such upper bounds on the number of transcription factors should exist, in order to minimize cross-binding errors between transcription factors. This theory further predicts that factors with similar binding sequences should tend to have similar biological effect, so that errors based on mis-recognition are minimal. We present evidence that transcription factors with similar binding sequences tend to regulate genes with similar biological functions, supporting this prediction.

Conclusion

The present study suggests limits on the transcription factor repertoire of cells, and suggests coding constraints that might apply more generally to the mapping between binding sites and biological function.

Cytogenetic mapping of immunoglobulin heavy chain genes in Antarctic fish

The chromosomal location of the IgH locus has been analyzed in several bony fish of the Antarctic perciform group Notothenioidei. Two IgH probes were prepared from the species Trematomus bernacchii (family Nototheniidae, tribe Trematominae) and mapped onto the chromosomes of ten species belonging to the same genus (Trematomus) and in two outgroups, through one-color and two-color FISH. A single location of the IgH locus was found in the majority of the species examined, including the outgroups, whereas in four of them the IgH genes splited to two chromosomal loci. RT-PCR experiments revealed the presence of three allelic sequences in T. newnesi, a species in which the IgH genes were organized in two chromosomal loci. Possible pathways leading to IgH genes duplication during the diversification of trematomine fishes were inferred from the analysis of the FISH patterns in a phylogenetic context. The present work provides the first comprehensive picture of IgH genes organization at chromosomal level in a bony fish group.

Single nucleotide polymorphism (SNP) discovery in duplicated genomes: intron-primed exon-crossing (IPEC) as a strategy for avoiding amplification of duplicated loci in Atlantic salmon (Salmo salar) and other salmonid fishes

Background

Single nucleotide polymorphisms (SNPs) represent the most abundant type of DNA variation in the vertebrate genome, and their applications as genetic markers in numerous studies of molecular ecology and conservation of natural populations are emerging. Recent large-scale sequencing projects in several fish species have provided a vast amount of data in public databases, which can be utilized in novel SNP discovery in salmonids. However, the suggested duplicated nature of the salmonid genome may hamper SNP characterization if the primers designed in conserved gene regions amplify multiple loci.

Results

Here we introduce a new intron-primed exon-crossing (IPEC) method in an attempt to overcome this duplication problem, and also evaluate different priming methods for SNP discovery in Atlantic salmon (Salmo salar) and other salmonids. A total of 69 loci with differing priming strategies were screened in S. salar, and 27 of these produced ~13 kb of high-quality sequence data consisting of 19 SNPs or indels (one per 680 bp). The SNP frequency and the overall nucleotide diversity (3.99 × 10^-4) in S. salar was lower than reported in a majority of other organisms, which may suggest a relative young population history for Atlantic salmon. A subset of primers used in cross-species analyses revealed considerable variation in the SNP frequencies and nucleotide diversities in other salmonids.

Conclusion

Sequencing success was significantly higher with the new IPEC primers; thus the total number of loci to screen in order to identify one potential polymorphic site was six times less with this new strategy. Given that duplication may hamper SNP discovery in some species, the IPEC method reported here is an alternative way of identifying novel polymorphisms in such cases.

Gene loss and evolutionary rates following whole-genome duplication in teleost fishes

Teleost fishes provide the first unambiguous support for ancient whole-genome duplication in an animal lineage. Studies in yeast or plants have shown that the effects of such duplications can be mediated by a complex pattern of gene retention and changes in evolutionary pressure. To explore such patterns in fishes, we have determined by phylogenetic analysis the evolutionary origin of 675 Tetraodon duplicated genes assigned to chromosomes, using additional data from other species of actinopterygian fishes. The subset of genes, which was retained in double after the genome duplication, is enriched in development, signaling, behavior, and regulation functional categories. The evolutionary rate of duplicate fish genes appears to be determined by 3 forces: 1) fish proteins evolve faster than mammalian orthologs; 2) the genes kept in double after genome duplication represent the subset under strongest purifying selection; and 3) following duplication, there is an asymmetric acceleration of evolutionary rate in one of the paralogs. These results show that similar mechanisms are at work in fishes as in yeast or plants and provide a framework for future investigation of the consequences of duplication in fishes and other animals.

Uneven evolutionary rates of bradykinin B1 and B2 receptors in vertebrate lineages

Bradykinin acts through two receptor subtypes in mammals and generates a variety of responses including pain, inflammation and hypotension. The evolutionary history of the bradykinin system has been unclear due to shortage of information outside mammals. We describe here two receptor subtypes and the bradykinin precursor in three species of bony fish (the zebrafish Danio rerio, the Japanese pufferfish Takifugu rubripes, and the green spotted pufferfish Tetraodon nigroviridis) and chicken and analyze the relationships to mammals by a combination of phylogeny, conserved synteny and exon-intron organization. All of these species have two receptor genes located close to each other in a tandem formation, with the B2 gene 5' to the B1 gene, in chromosomal regions displaying conserved synteny between the species (albeit conservation of synteny in zebrafish is still unclear due to poor genome assembly). The evolutionary rate differs between the two genes as well as between lineages leading to differing pharmacological properties for both B1 and B2 across vertebrate classes. Also the bradykinin precursor gene was identified in all of these species in a chromosome region with conserved synteny. The tissue distribution of mRNA in T. rubripes is similar for B1 and B2, suggesting more similar regulation for the two genes than in mammals. In conclusion, the receptor tandem duplication predates the divergence of ray-finned fish and tetrapods and no additional duplicates of the receptors or bradykinin seem to have survived the ray-finned fish tetraploidization.

Contrast between extensive variation of 28S rDNA and stability of 5S rDNA and telomeric repeats in the diploid-polyploid Squalius alburnoides complex and in its maternal ancestor Squalius pyrenaicus (Teleostei, Cyprinidae)

The diploid-polyploid Squalius alburnoides complex resulted from interspecific hybridization. The chromosomal mapping of 28S and 5S ribosomal genes and of (TTAGGG)n telomeric repeats was performed on specimens from the complex and from the sympatric bisexual species S. pyrenaicus (the complex maternal ancestor) as part of an investigation of the evolutionary relationships between genomic constitutions and the consequences of the ongoing polyploidization process in terms of chromosome reshaping. Contrasting results were obtained. While results with 5S rDNA and telomeric probes gave an impression of genomic stability, the variability detected with 28S rDNA probe suggested quite the opposite. The 5S rDNA probe mapped constantly to three chromosomes per haploid genome with apparently conserved locations in morphologically similar chromosomes; conversely, prominent intra- and inter-individual variations of 28S rDNA and of syntenic sites with 5S rDNA were detected with regard to number, size and location. Hypotheses for the causes of such polymorphisms are discussed. The terminal position of most 28S rDNA sites and the absence of detectable interstitial telomeric sequences suggest a mechanism that does not involve major chromosomal rearrangements. These fishes share similar patterns for the studied cytogenetic markers which may be taken as evidence of an apparent stability that may be hiding extensive and subtle genome variations that are possibly related to an ongoing evolutionary process of genome tetraploidization and speciation.

Identification of novel rainbow trout (Onchorynchus mykiss) chemokines, CXCd1 and CXCd2: mRNA expression after Yersinia ruckeri vaccination and challenge

Chemokines play important roles in controlling leukocyte trafficking under normal and inflammatory conditions. Sixteen CXC chemokines have been identified in the human and mouse genomes, while considerably fewer teleost fish CXC chemokines have been reported. Here, we describe a novel clade of trout (Onchorynchus mykiss) CXC chemokines, designated Onmy CXCd, and we identify a novel gene, CXCd1, and a putative duplicate, CXCd2. The trout CXCd proteins contain 112 amino acids and the CXCd1 gene is comprised of four exons and three introns. Constitutive CXCd mRNA expression was detected in skin, gill, visceral fat, and posterior kidney tissues, while low transcript levels were present in the anterior kidney and spleen. Spleen CXCd transcript abundance increased 1 day after bath vaccination (fourfold) and subsided to basal levels by 7 days postvaccination. Challenge with viable Yersinia ruckeri induced expression of trout CXCd RNA up to ninefold in the spleen. The number of viable Y. ruckeri were significantly correlated with CXCd gene transcript abundance (P = 0.0051, Spearman correlation 0.497, n = 30 fish), and fish with the highest bacterial loads had the highest CXCd expression. In contrast, pro-inflammatory cytokine IL-1-beta2 mRNA levels were elevated in fish infected with low numbers of Y. ruckeri, while diminishing in heavily infected fish. CXCd mRNA expression was not increased in rainbow trout infected with infectious hematopoietic necrosis virus, suggesting that up-regulation may be pathogen-specific. Taken together, these results indicate that CXCd transcript elevation follows the pro-inflammatory cytokine response to Y. ruckeri and may be a relevant immunological marker of exposure.

Identification and comparative expression analysis of a second wt1 gene in zebrafish

The Wilms' tumor suppressor gene wt1 encodes a zinc-finger transcription factor that plays an important role in the development of the mammalian genitourinary system. Mutations in WT1 in humans lead to anomalies of kidney and gonad development and cause Wilms' tumor, a pediatric kidney cancer. The inactivation of both wt1 alleles in mice gives rise to multiple organ defects, among them agenesis of kidney, spleen, and gonads. In zebrafish, an ortholog of wt1 has been described that is expressed in the pronephric field and is later restricted to the podocytes. Here, we report the existence of a second wt1 gene in zebrafish, which we have named wt1b (we named the initial gene wt1a). The overall sequence identity of the two Wt1 proteins is 70% and 92% between the zinc-finger regions, respectively. In contrast to wt1a, wt1b is expressed from the earliest stages of development onward, albeit at low levels. Both wt1a and wt1b are expressed in the intermediate mesoderm, with wt1b being restricted to a smaller area lying at the caudal end of the wt1a expression domain. In adult fish, high expression levels for both genes can be found in gonads, kidney, heart, spleen, and muscle.

Asymmetric evolution in two fish-specifically duplicated receptor tyrosine kinase paralogons involved in teleost coloration

The occurrence of a fish-specific genome duplication (FSGD) in the lineage leading to teleost fishes is widely accepted, but the consequences of this event remain elusive. Teleosts, and the cichlid fishes from the species flocks in the East African Great Lakes in particular, evolved a unique complexity and diversity of body coloration and color patterning. Several genes involved in pigment cell development have been retained in duplicate copies in the teleost genome after the FSGD. Here we investigate the evolutionary fate of one of these genes, the type III receptor tyrosine kinase (RTK) colony-stimulating factor 1 receptor (csf1r). We isolated and shotgun sequenced two paralogous csf1r genes from a bacterial artificial chromosome library of the cichlid fish Astatotilapia burtoni that are both linked to paralogs of the pdgfr beta gene, another type III RTK. Two pdgfr beta-csf1r paralogons were also identified in the genomes of pufferfishes and medaka, and our phylogenetic analyses suggest that the pdgfr beta-csf1r locus was duplicated during the course of the FSGD. Comparisons of teleosts and tetrapods suggest asymmetrical divergence at different levels of genomic organization between the teleost-specific pdgfr beta-csf1r paralogons, which seem to have evolved as coevolutionary units. The high-evolutionary rate in the teleost B-paralogon, consisting of csf1rb and pdgfr betab, further suggests neofunctionalization by functional divergence of the extracellular, ligand-binding region of these cell-surface receptors. Finally, we hypothesize that genome duplications and the associated expansion of the RTK family might be causally linked to the evolution of coloration in vertebrates and teleost fishes in particular.

Comparative analysis of her genes during fish somitogenesis suggests a mouse/chick-like mode of oscillation in medaka

Somitogenesis is the key developmental step, which divides the vertebrate body axis into segmentally repeated structures. It requires an intricate process of pre-patterning, which is driven by an oscillator mechanism consisting of the Delta-Notch pathway and various hairy- and Enhancer of split-related (her) genes. The subset of her genes, which are necessary to set up the segmentation clock, reveal a complex scenario of interactions. To understand which her genes are essential core players in this process, we compared the expression patterns of somitogenesis-relevant her genes in zebrafish and medaka (Oryzias latipes). Most of the respective medaka genes (Ol-her) are duplicated like what has been shown for zebrafish (Dr-her) and pufferfish genes (Fr-her). However, zebrafish genes show some additional copies and significant differences in expression patterns. For the paralogues Dr-her1 and Dr-her11, only one copy exists in the medaka (Ol-her1/11), which combines the expression patterns found for both zebrafish genes. In contrast to Dr-her5, the medaka orthologue appears to play a role in somitogenesis because it is expressed in the presomitic mesoderm (PSM). PSM expression also suggests a role for both Ol-her13 genes, homologues of mouse Hes6 (mHes6), in this process, which would be consistent with a conserved mHes6 homologue gear in the segmentation clock exclusively in lower vertebrates. Members of the mHes5 homologue group seem to be involved in somite formation in all vertebrates (e.g. Dr- and Ol-her12), although different paralogues are additionally recruited in zebrafish (e.g. Dr-her15) and medaka (e.g. Ol-her4). We found that the linkage between duplicates is strongly conserved between pufferfish and medaka and less well conserved in zebrafish. Nevertheless, linkage and orientation of several her duplicates are identical in all three species. Therefore, small-scale duplications must have happened before whole genome duplication occurred in a fish ancestor. Expression of multiple stripes in the intermediate PSM, characteristic for the zebrafish orthologues, is absent in all somitogenesis-related her genes of the medaka. In fact, the expression mode of Ol-her1/11 and Ol-her5 indicates dynamism similar to the hairy clock genes in chicken and mouse. This suggests that Danio rerio shows a rather derived clock mode when compared to other fish species and amniotes or that, alternatively, the clock mode evolved independently in zebrafish, medaka and mouse or chicken.

Sex chromosomes and male ornaments: a comparative evaluation in ray-finned fishes

Theory predicts that the mechanism of genetic sex determination can substantially influence the evolution of sexually selected traits. For example, female heterogamety (ZZ/ZW) can favour the evolution of extreme male traits under Fisher's runaway model of sexual selection. We empirically test whether the genetic system of sex determination has played a role in the evolution of exaggerated male ornaments in actinopterygiian fishes, a clade in which both female-heterogametic and male-heterogametic systems of sex determination have evolved multiple times. Using comparative methods both uncorrected and corrected for phylogenetic non-independence, we detected no significant correlation between sex-chromosome systems and sexually selected traits in males. Results suggest that sex-determination mechanism is at best a relatively minor factor affecting the outcomes of sexual selection in ray-finned fishes.

The first radiation hybrid map of a perch-like fish: the gilthead seabream (Sparus aurata L)

Among Teleosts, Perciformes are the largest order of fishes and include numerous species of commercial importance. Perciformes also comprise species of primary interest for evolutionary studies and analysis of the sex determination systems and sex chromosome plasticity. Unfortunately, genomics tools and resources for Perciformes remain to be developed. Here, we report the production of a seabream whole-genome radiation hybrid (RH) panel in which quality was ascertained by the construction of a 2-Mb-resolution RH map. The map encompasses 440 markers (288 microsatellites, 82 gene-based markers, and 70 STS) suitable for linkage analysis and comparative mapping studies. Achievement of a RH panel and a whole-genome RH map should contribute to establishing seabream as a fish model among the Perciformes and should be of importance in aquaculture for marker-assisted selection, improvement of growth performance, and disease management. Development of RH maps in a cost-effective manner for other fishes with the described methodology will offer a powerful approach in aquaculture and will provide extended capabilities for comparing vertebrate genome evolution.

The evolution of metazoan axial properties

Renewed interest in the developmental basis of organismal complexity, and the emergence of new molecular tools, is improving our ability to study the evolution of metazoan body plans. The most substantial changes in body-plan organization occurred early in metazoan evolution; new model systems for studying basal metazoans are now being developed, and total-genome-sequencing initiatives are underway for at least three of the four most important taxa. The elucidation of how the gene networks that are involved in axial organization, germ-layer formation and cell differentiation are used differently during development is generating a more detailed understanding of the events that have led to the current diversity of multicellular life.

Evolutionary conservation and diversification of Rh family genes and proteins

Rhesus (Rh) proteins were first identified in human erythroid cells and recently in other tissues. Like ammonia transporter (Amt) proteins, their only homologues, Rh proteins have the 12 transmembrane-spanning segments characteristic of transporters. Many think Rh and Amt proteins transport the same substrate, NH(3)/NH(4)(+), whereas others think that Rh proteins transport CO(2) and Amt proteins NH(3). In the latter view, Rh and Amt are different biological gas channels. To reconstruct the phylogeny of the Rh family and study its coexistence with and relationship to Amt in depth, we analyzed 111 Rh genes and 260 Amt genes. Although Rh and Amt are found together in organisms as diverse as unicellular eukaryotes and sea squirts, Rh genes apparently arose later, because they are rare in prokaryotes. However, Rh genes are prominent in vertebrates, in which Amt genes disappear. In organisms with both types of genes, Rh had apparently diverged away from Amt rapidly and then evolved slowly over a long period. Functionally divergent amino acid sites are clustered in transmembrane segments and around the gas-conducting lumen recently identified in Escherichia coli AmtB, in agreement with Rh proteins having new substrate specificity. Despite gene duplications and mutations, the Rh paralogous groups all have apparently been subject to strong purifying selection indicating functional conservation. Genes encoding the classical Rh proteins in mammalian red cells show higher nucleotide substitution rates at nonsynonymous codon positions than other Rh genes, a finding that suggests a possible role for these proteins in red cell morphogenetic evolution.

Molecular evolution of ankyrin: gain of function in vertebrates by acquisition of an obscurin/titin-binding-related domain

Ankyrins form a family of modular adaptor proteins that link between integral membrane proteins and the cytoskeleton. They evolved within the Metazoa as an adaptation for organizing membrane microstructure and directing membrane traffic. Molecular cloning has identified one Caenorhabditis elegans (unc-44), two Drosophila (Dank1, Dank2), and three mammalian (Ank1, Ank2, Ank3) genes. We have previously identified a 76-amino acid (aa) alternatively spliced sequence that is present in muscle polypeptides encoded by the rat Ank3 gene. A closely related sequence in a muscle Ank1 product binds the cytoskeletal muscle proteins obscurin and titin. This obscurin/titin-binding-related domain (OTBD) contains repeated modules of 18 aa: three are encoded by Ank1 and Ank2, two by Ank3; this pattern is conserved throughout vertebrate ankyrin genes. The C. elegans ankyrin, UNC-44, contains one 18-aa module as does the ankyrin gene in the urochordate Ciona intestinalis, but the insect ankyrins contain none. Our data indicate that an ancestral ankyrin acquired an 18-aa module which was preserved in the Ecdysozoa/deuterostome divide, but it was subsequently lost from arthropods. Successive duplications of the module led to a gain of function in vertebrates as it acquired obscurin/titin-binding activity. We suggest that the OTBD represents an adaptation of the cytoskeleton that confers muscle cells with resilience to the forces associated with vertebrate life.

From 2R to 3R: evidence for a fish-specific genome duplication (FSGD)

An important mechanism for the evolution of phenotypic complexity, diversity and innovation, and the origin of novel gene functions is the duplication of genes and entire genomes. Recent phylogenomic studies suggest that, during the evolution of vertebrates, the entire genome was duplicated in two rounds (2R) of duplication. Later, approximately 350 mya, in the stem lineage of ray-finned (actinopterygian) fishes, but not in that of the land vertebrates, a third genome duplication occurred-the fish-specific genome duplication (FSGD or 3R), leading, at least initially, to up to eight copies of the ancestral deuterostome genome. Therefore, the sarcopterygian (lobe-finned fishes and tetrapods) genome possessed originally only half as many genes compared to the derived fishes, just like the most-basal and species-poor lineages of extant fishes that diverged from the fish stem lineage before the 3R duplication. Most duplicated genes were secondarily lost, yet some evolved new functions. The genomic complexity of the teleosts might be the reason for their evolutionary success and astounding biological diversity.