Zebrafish comparative genomics and the origins of vertebrate chromosomes

The LN54 radiation hybrid map of zebrafish expressed sequences

To increase the density of a gene map of the zebrafish, Danio rerio, we have placed 3119 expressed sequence tags (ESTs) and cDNA sequences on the LN54 radiation hybrid (RH) panel. The ESTs and genes mapped here join 748 SSLp markers and 459 previously mapped genes and ESTs, bringing the total number of markers on the LN54 RH panel to 4226. Addition of these new markers brings the total LN54 map size to 14,372 cR, with 118 kb/cR. The distribution of ESTs according to linkage groups shows relatively little variation (minimum, 73; maximum, 201). This observation, combined with a relatively uniform size for zebrafish chromosomes, as previously indicated by karyotyping, indicates that there are no especially gene-rich or gene-poor chromosomes in this species. We developed an algorithm to provide a semiautomatic method for the selection of additional framework markers for the LN54 map. This algorithm increased the total number of framework markers to 1150 and permitted the mapping of a high percentage of sequences that could not be placed on a previous version of the LN54 map. The increased concentration of expressed sequences on the LN54 map of the zebrafish genome will facilitate the molecular characterization of mutations in this species.

In situ hybridization screen in zebrafish for the selection of genes encoding secreted proteins

An in situ hybridization expression screen using a signal sequence trap system has been conducted in zebrafish to isolate cDNAs that encode secreted proteins. Random clones (secreted expressed sequence tags; sESTs) were sequenced from zebrafish embryonic (18-24 hr postfertilization) and adult kidney libraries. From the two RNA sources, 627 random sEST cDNAs were identified as being homologous or identical to known genes and 166 clones encode currently unidentified genes. The sESTs represent a broad range of enzymes and other regulatory molecules. Whole-mount in situ hybridization analysis was carried out by using antisense probes generated from 244 selected sESTs, and a range of expression patterns was obtained. Genetic mapping undertaken with sEST sequences demonstrated that assignment of map position was attainable by using 5' primers. The signal sequence trap system used in this work has yielded a range of cDNAs that encode secreted proteins and, together with analysis of patterns of expression and genetic mapping, has the potential to facilitate analysis of signaling pathways central to development and physiology.

Identification and characterization of two distinct GnRH receptor subtypes in a teleost, the medaka Oryzias latipes

We report the identification and characterization of two distinct GnRH receptor (GnRH-R) subtypes, designated GnRH-R1 and GnRH-R2, in a model teleost, the medaka Oryzias latipes. These seven-transmembrane receptors of the medaka contain a cytoplasmic C-terminal tail, which has been found in all other nonmammalian GnRH-Rs cloned to date. The GnRH-R1 gene is composed of three exons separated by two introns, whereas the GnRH-R2 gene has an additional intron and therefore consists of four exons and three introns. The GnRH-R1 and GnRH-R2 genes, both of which exist as single-copy genes in the medaka genome, were mapped to linkage groups 3 and 16, respectively. Inositol phosphate assays using COS-7 cells transfected with GnRH-R1 and GnRH-R2 demonstrated that they had remarkably different ligand sensitivities, although both receptors showed highest preference for chicken-II-type GnRH. Phylogenetic analysis showed the presence of three paralogous lineages for vertebrate GnRH-Rs and indicated that neither GnRH-R1 nor GnRH-R2 is the medaka ortholog to mammalian GnRH-Rs that lack a cytoplasmic tail. This, together with an observation that medaka-type GnRH had low affinity for GnRH-R1 and GnRH-R2, suggests that a third GnRH-R may exist in the medaka.

Comparative genomics of the SOX9 region in human and Fugu rubripes: conservation of short regulatory sequence elements within large intergenic regions

Campomelic dysplasia (CD), a human skeletal malformation syndrome with XY sex reversal, is caused by heterozygous mutations in and around the gene SOX9. SOX9 has an extended 5' control region, as indicated by CD translocation breakpoints scattered over 1 Mb proximal to SOX9 and by expression data from mice transgenic for human SOX9-spanning yeast artificial chromosomes. To identify long-range regulatory elements within the SOX9 5' control region, we compared approximately 3.7 Mb and 195 kb of sequence around human and Fugu rubripes SOX9, respectively. We identified only seven and five protein-coding genes in the human and F. rubripes sequences, respectively. Four of the F. rubripes genes have been mapped in humans; all reside on chromosome 17 but show extensive intrachromosomal gene shuffling compared with the gene order in F. rubripes. In both species, very large intergenic distances separate SOX9 from its directly flanking genes: 2 Mb and 500 kb on either side of SOX9 in humans, and 68 and 97 kb on either side of SOX9 in F. rubripes. Comparative sequence analysis of the intergenic regions revealed five conserved elements, E1-E5, up to 290 kb 5' to human SOX9 and up to 18 kb 5' to F. rubripes SOX9, and three such elements, E6-E8, 3' to SOX9. Where available, mouse sequences confirm conservation of the elements. From the yeast artificial chromosome transgenic data, elements E3-E5 are candidate enhancers for SOX9 expression in limb and vertebral column, and 8 of 10 CD translocation breakpoints separate these elements from SOX9.

Partitioning of tissue expression accompanies multiple duplications of the Na+/K+ ATPase alpha subunit gene

Vertebrate genomes contain multiple copies of related genes that arose through gene duplication. In the past it has been proposed that these duplicated genes were retained because of acquisition of novel beneficial functions. A more recent model, the duplication-degeneration-complementation hypothesis (DDC), posits that the functions of a single gene may become separately allocated among the duplicated genes, rendering both duplicates essential. Thus far, empirical evidence for this model has been limited to the engrailed and sox family of developmental regulators, and it has been unclear whether it may also apply to ubiquitously expressed genes with essential functions for cell survival. Here we describe the cloning of three zebrafish alpha subunits of the Na(+),K(+)-ATPase and a comprehensive evolutionary analysis of this gene family. The predicted amino acid sequences are extremely well conserved among vertebrates. The evolutionary relationships and the map positions of these genes and of other alpha-like sequences indicate that both tandem and ploidy duplications contributed to the expansion of this gene family in the teleost lineage. The duplications are accompanied by acquisition of clear functional specialization, consistent with the DDC model of genome evolution.

The medaka rs-3 locus required for scale development encodes ectodysplasin-A receptor

The bodies of most teleost fish species are covered with specialized subepithelial structures known as scales. The scale is an epithelial appendage that differentiates from the dermal mesenchyme. Mammals, on the other hand, have no scales, but instead their bodies are covered with hair. Although their appearances are quite different, scales and hair can be considered structurally similar in that both of them are epithelial appendages distributed over the body surface in an orderly pattern. This analogy suggests that they may have the same evolutionary origin. But, to date, no molecular evidence has been presented that links scales and hair. A mutation at the rs-3 locus of medaka (Oryzias latipes) leads to almost complete loss of scales. We demonstrated that the rs-3 locus encodes ectodysplasin-A receptor (EDAR), which is required for the initiation of hair development in mammals. We identified a novel transposon inserted in the first intron of EDAR, which causes aberrant splicing. This work shows that EDAR is required for scale development in fish and suggests that it is an evolutionarily conserved molecule that is required for the development of epithelial appendages in vertebrates.

Genome duplication, divergent resolution and speciation

What are the evolutionary consequences of gene duplication? One answer is speciation, according to a model initially called Reciprocal Silencing and recently expanded and renamed Divergent Resolution. This model shows how the loss of different copies of a duplicated gene in allopatric populations (divergent resolution) can promote speciation by genetically isolating these populations should they become reunited. Genome duplication events produce thousands of duplicated genes. Therefore, lineages with a history of genome duplication might have been especially prone to speciation via divergent resolution.

Human chromosome 19 and related regions in mouse: conservative and lineage-specific evolution

To illuminate the function and evolutionary history of both genomes, we sequenced mouse DNA related to human chromosome 19. Comparative sequence alignments yielded confirmatory evidence for hypothetical genes and identified exons, regulatory elements, and candidate genes that were missed by other predictive methods. Chromosome-wide comparisons revealed a difference between single-copy HSA19 genes, which are overwhelmingly conserved in mouse, and genes residing in tandem familial clusters, which differ extensively in number, coding capacity, and organization between the two species. Finally, we sequenced breakpoints of all 15 evolutionary rearrangements, providing a view of the forces that drive chromosome evolution in mammals.

The repertoire of Na,K-ATPase alpha and beta subunit genes expressed in the zebrafish, Danio rerio

We have identified a cohort of zebrafish expressed sequence tags encoding eight Na,K-ATPase alpha subunits and five beta subunits. Sequence comparisons and phylogenetic analysis indicate that five of the zebrafish alpha subunit genes comprise an alpha1-like gene subfamily and two are orthologs of the mammalian alpha3 subunit gene. The remaining alpha subunit clone is most similar to the mammalian alpha2 subunit. Among the five beta subunit genes, two are orthologs of the mammalian beta1 isoform, one represents a beta2 ortholog, and two are orthologous to the mammalian beta3 subunit. Using zebrafish radiation hybrid and meiotic mapping panels, we determined linkage assignments for each alpha and beta subunit gene. Na,K-ATPase genes are dispersed in the zebrafish genome with the exception of four of the alpha1-like genes, which are tightly clustered on linkage group 1. Comparative mapping studies indicate that most of the zebrafish Na,K-ATPase genes localize to regions of conserved synteny between zebrafish and humans. The expression patterns of Na,K-ATPase alpha and beta subunit genes in zebrafish are quite distinctive. No two alpha or beta subunit genes exhibit the same expression profile. Together, our data imply a very high degree of Na,K-ATPase isoenzyme heterogeneity in zebrafish, with the potential for 40 structurally distinct alpha/beta subunit combinations. Differences in expression patterns of alpha and beta subunits suggest that many of the isoenzymes are also likely to exhibit differences in functional properties within specific cell and tissue types. Our studies form a framework for analyzing structure function relationships for sodium pump isoforms using reverse genetic approaches.

Genome duplication events and functional reduction of ploidy levels in sturgeon (Acipenser, Huso and Scaphirhynchus)

Sturgeon (order Acipenserformes) provide an ideal taxonomic context for examination of genome duplication events. Multiple levels of ploidy exist among these fish. In a novel microsatellite approach, data from 962 fish from 20 sturgeon species were used for analysis of ploidy in sturgeon. Allele numbers in a sample of individuals were assessed at six microsatellite loci. Species with approximately 120 chromosomes are classified as functional diploid species, species with approximately 250 chromosomes as functional tetraploid species, and with approximately 500 chromosomes as functional octaploids. A molecular phylogeny of the sturgeon was determined on the basis of sequences of the entire mitochondrial cytochrome b gene. By mapping the estimated levels of ploidy on this proposed phylogeny we demonstrate that (I) polyploidization events independently occurred in the acipenseriform radiation; (II) the process of functional genome reduction is nearly finished in species with approximately 120 chromosomes and more active in species with approximately 250 chromosomes and approximately 500 chromosomes; and (III) species with approximately 250 and approximately 500 chromosomes arose more recently than those with approximately 120 chromosomes. These results suggest that gene silencing, chromosomal rearrangements, and transposition events played an important role in the acipenseriform genome formation. Furthermore, this phylogeny is broadly consistent with previous hypotheses but reveals a highly supported oceanic (Atlantic-Pacific) subdivision within the Acipenser/Huso complex.

Evolution of mammalian genome organization inferred from comparative gene mapping

Comparative genome analyses, including chromosome painting in over 40 diverse mammalian species, ordered gene maps from several representatives of different mammalian and vertebrate orders, and large-scale sequencing of the human and mouse genomes are beginning to provide insight into the rates and patterns of chromosomal evolution on a whole-genome scale, as well as into the forces that have sculpted the genomes of extant mammalian species.

Immune-type receptor genes in zebrafish share genetic and functional properties with genes encoded by the mammalian leukocyte receptor cluster

An extensive, highly diversified multigene family of novel immune-type receptor (nitr) genes has been defined in Danio rerio (zebrafish). The genes are predicted to encode type I transmembrane glycoproteins consisting of extracellular variable (V) and V-like C2 (V/C2) domains, a transmembrane region and a cytoplasmic tail. All of the genes examined encode immunoreceptor tyrosine-based inhibition motifs in the cytoplasmic tail. Radiation hybrid panel mapping and analysis of a deletion mutant line (b240) indicate that a minimum of approximately 40 nitr genes are contiguous in the genome and span approximately 0.6 Mb near the top of zebrafish linkage group 7. One flanking region of the nitr gene complex shares conserved synteny with a region of mouse chromosome 7, which shares conserved synteny with human 19q13.3-q13.4 that encodes the leukocyte receptor cluster. Antibody-induced crosslinking of Nitrs that have been introduced into a human natural killer cell line inhibits the phosphorylation of mitogen-activated protein kinase that is triggered by natural killer-sensitive tumor target cells. Nitrs likely represent intermediates in the evolution of the leukocyte receptor cluster.

Conserved vertebrate chromosome segments in the large salamander genome

Urodele amphibians (salamanders) are important models for embryological, physiological, and natural history research and are also a biomedically important group because they are the only vertebrates capable of regenerating entire organ systems. To enhance the utility of salamanders for biomedical research and for understanding genome evolution, genetic linkage analysis was used to identify chromosome segments that are homologous between ambystomatid salamanders and distantly related vertebrate model organisms. A total of 347 loci (AFLPs, RAPDs, and protein-coding loci) were mapped using an interspecific meiotic mapping panel (Ambystoma mexicanum and A. tigrinum tigrinum; family Ambystomatidae). Genome size in Ambystoma was estimated to be 7291 cM, the largest linkage map estimate reported for any organism. However, the relatively large size of the salamander genome did not hinder efforts to map and identify conserved syntenies from a small sample of 24 protein-coding loci. Chromosomal segments that are conserved between fishes and mammals are also conserved in these salamanders. Thus, comparative gene mapping appears to be an efficient strategy for identifying orthologous loci between ambystomatid salamanders and genomically well-characterized vertebrate model organisms.

Fishing for lymphoid genes

Thymic organogenesis and T-cell lymphopoiesis are crucial interdependent processes that establish a functional vertebrate immune system. The current understanding of vertebrate thymic development during embryogenesis remains incomplete and would benefit from novel approaches. The zebrafish Danio rerio is a powerful developmental and genetic system for the dissection of early events in the ontogeny of the immune system. Forward genetic screens have uncovered genes involved in hematopoiesis, and specific screens are being designed to examine the genes that regulate T-cell development and the origin of the thymus. Studies of the zebrafish should improve our understanding of lymphoid development in vertebrates.

A microsatellite-based genetic linkage map for channel catfish, Ictalurus punctatus

Microsatellite loci were identified in channel catfish gene sequences or random clones from a small insert genomic DNA library. Outbred populations of channel catfish contained an average of eight alleles per locus and an average heterozygosity of 0.70. A genetic linkage map of the channel catfish genome (N = 29) was constructed from two reference families. A total of 293 microsatellite loci were polymorphic in one or both families, with an average of 171 informative meioses per locus. Nineteen type I loci, 243 type II loci, and one EST were placed in 32 multipoint linkage groups covering 1958 cM. Nine more type II loci were contained in three two-point linkage groups covering 24.5 cM. Twenty-two type II loci remained unlinked. Multipoint linkage groups ranged in size from 11.9 to 110.5 cM with an average intermarker distance of 8.7 cM. Seven microsatellite loci were closely linked with the sex-determining locus. The microsatellite loci and genetic linkage map will increase the efficiency of selective breeding programs for channel catfish.

Yesterday's polyploids and the mystery of diploidization

Thirty years after Susumu Ohno proposed that vertebrate genomes are degenerate polyploids, the extent to which genome duplication contributed to the evolution of the vertebrate genome, if at all, is still uncertain. Sequence-level studies on model organisms whose genomes show clearer evidence of ancient polyploidy are invaluable because they indicate what the evolutionary products of genome duplication can look like. The greatest mystery is the molecular basis of diploidization, the evolutionary process by which a polyploid genome turns into a diploid one.

A comparative map of the zebrafish genome

Zebrafish mutations define the functions of hundreds of essential genes in the vertebrate genome. To accelerate the molecular analysis of zebrafish mutations and to facilitate comparisons among the genomes of zebrafish and other vertebrates, we used a homozygous diploid meiotic mapping panel to localize polymorphisms in 691 previously unmapped genes and expressed sequence tags (ESTs). Together with earlier efforts, this work raises the total number of markers scored in the mapping panel to 2119, including 1503 genes and ESTs and 616 previously characterized simple-sequence length polymorphisms. Sequence analysis of zebrafish genes mapped in this study and in prior work identified putative human orthologs for 804 zebrafish genes and ESTs. Map comparisons revealed 139 new conserved syntenies, in which two or more genes are on the same chromosome in zebrafish and human. Although some conserved syntenies are quite large, there were changes in gene order within conserved groups, apparently reflecting the relatively frequent occurrence of inversions and other intrachromosomal rearrangements since the divergence of teleost and tetrapod ancestors. Comparative mapping also shows that there is not a one-to-one correspondence between zebrafish and human chromosomes. Mapping of duplicate gene pairs identified segments of 20 linkage groups that may have arisen during a genome duplication that occurred early in the evolution of teleosts after the divergence of teleost and mammalian ancestors. This comparative map will accelerate the molecular analysis of zebrafish mutations and enhance the understanding of the evolution of the vertebrate genome.

A Comparative Map of the Zebrafish Genome

Zebrafish mutations define the functions of hundreds of essential genes in the vertebrate genome. To accelerate the molecular analysis of zebrafish mutations and to facilitate comparisons among the genomes of zebrafish and other vertebrates, we used a homozygous diploid meiotic mapping panel to localize polymorphisms in 691 previously unmapped genes and expressed sequence tags (ESTs). Together with earlier efforts, this work raises the total number of markers scored in the mapping panel to 2119, including 1503 genes and ESTs and 616 previously characterized simple-sequence length polymorphisms. Sequence analysis of zebrafish genes mapped in this study and in prior work identified putative human orthologs for 804 zebrafish genes and ESTs. Map comparisons revealed 139 new conserved syntenies, in which two or more genes are on the same chromosome in zebrafish and human. Although some conserved syntenies are quite large, there were changes in gene order within conserved groups, apparently reflecting the relatively frequent occurrence of inversions and other intrachromosomal rearrangements since the divergence of teleost and tetrapod ancestors. Comparative mapping also shows that there is not a one-to-one correspondence between zebrafish and human chromosomes. Mapping of duplicate gene pairs identified segments of 20 linkage groups that may have arisen during a genome duplication that occurred early in the evolution of teleosts after the divergence of teleost and mammalian ancestors. This comparative map will accelerate the molecular analysis of zebrafish mutations and enhance the understanding of the evolution of the vertebrate genome.