Zebrafish Lineage Research Articles

Differential allelic representation (DAR) identifies candidate eQTLs and improves transcriptome analysis.

In comparisons between mutant and wild-type genotypes, transcriptome analysis can reveal the direct impacts of a mutation, together with the homeostatic responses of the biological system. Recent studies have highlighted that, when the effects of homozygosity for recessive mutations are studied in non-isogenic backgrounds, genes located proximal to the mutation on the same chromosome often appear over-represented among those genes identified as differentially expressed (DE). One hypothesis suggests that DE genes chromosomally linked to a mutation may not reflect functional responses to the mutation but, instead, result from an unequal distribution of expression quantitative trait loci (eQTLs) between sample groups of mutant or wild-type genotypes. This is problematic because eQTL expression differences are difficult to distinguish from genes that are DE due to functional responses to a mutation. Here we show that chromosomally co-located differentially expressed genes (CC-DEGs) are also observed in analyses of dominant mutations in heterozygotes. We define a method and a metric to quantify, in RNA-sequencing data, localised differential allelic representation (DAR) between those sample groups subjected to differential expression analysis. We show how the DAR metric can predict regions prone to eQTL-driven differential expression, and how it can improve functional enrichment analyses through gene exclusion or weighting-based approaches. Advantageously, this improved ability to identify probable eQTLs also reveals examples of CC-DEGs that are likely to be functionally related to a mutant phenotype. This supports a long-standing prediction that selection for advantageous linkage disequilibrium influences chromosome evolution. By comparing the genomes of zebrafish (Danio rerio) and medaka (Oryzias latipes), a teleost with a conserved ancestral karyotype, we find possible examples of chromosomal aggregation of CC-DEGs during evolution of the zebrafish lineage. Our method for DAR analysis requires only RNA-sequencing data, facilitating its application across new and existing datasets.

Generation and characterization of a novel gne Knockout Model in Zebrafish.

GNE Myopathy is a rare, recessively inherited neuromuscular worldwide disorder, caused by a spectrum of bi-allelic mutations in the human GNE gene. GNE encodes a bi-functional enzyme responsible for the rate-limiting step of sialic acid biosynthesis pathway. However, the process in which GNE mutations lead to the development of a muscle pathology is not clear yet. Cellular and mouse models for GNE Myopathy established to date have not been informative. Further, additional GNE functions in muscle have been hypothesized. In these studies, we aimed to investigate gne functions using zebrafish genetic and transgenic models, and characterized them using macroscopic, microscopic, and molecular approaches. We first established transgenic zebrafish lineages expressing the human GNE cDNA carrying the M743T mutation, driven by the zebrafish gne promoter. These fish developed entirely normally. Then, we generated a gne knocked-out (KO) fish using the CRISPR/Cas9 methodology. These fish died 8-10days post-fertilization (dpf), but a phenotype appeared less than 24h before death and included progressive body axis curving, deflation of the swim bladder and decreasing movement and heart rate. However, muscle histology uncovered severe defects, already at 5 dpf, with compromised fiber organization. Sialic acid supplementation did not rescue the larvae from this phenotype nor prolonged their lifespan. To have deeper insights into the potential functions of gne in zebrafish, RNA sequencing was performed at 3 time points (3, 5, and 7 dpf). Genotype clustering was progressive, with only 5 genes differentially expressed in gne KO compared to gne WT siblings at 3 dpf. Enrichment analyses of the primary processes affected by the lack of gne also at 5 and 7 dpf point to the involvement of cell cycle and DNA damage/repair processes in the gne KO zebrafish. Thus, we have established a gne KO zebrafish lineage and obtained new insights into gne functions. This is the only model where GNE can be related to clear muscle defects, thus the only animal model relevant to GNE Myopathy to date. Further elucidation of gne precise mechanism-of-action in these processes could be relevant to GNE Myopathy and allow the identification of novel therapeutic targets.

A Highly Divergent Picornavirus Infecting the Gut Epithelia of Zebrafish (Danio rerio) in Research Institutions Worldwide.

Zebrafish have been extensively used as a model system for research in vertebrate development and pathogen-host interactions. We describe the complete genome of a novel picornavirus identified during a viral metagenomics analysis of zebrafish gut tissue. The closest relatives of this virus showed identity of <20% in their P1 capsids and <36% in their RdRp qualifying zebrafish picornavirus-1 (ZfPV-1) as member of a novel genus with a proposed name of Cyprivirus. Reverse transcription (RT)-PCR testing of zebrafish from North America, Europe, and Asia showed ZfPV-1 to be globally distributed, being detected in 23 of 41 (56%) institutions tested. In situ hybridization of whole zebrafish showed viral RNA was restricted to a subset of enterocytes and cells in the subjacent lamina propria of the intestine and the intestinal mucosa. This naturally occurring and apparently asymptomatic infection (in wild-type zebrafish lineage AB) provides a natural infection system to study picornavirus-host interactions in an advanced vertebrate model organism. Whether ZfPV-1 infection affects any immunological, developmental, or other biological processes in wild-type or mutant zebrafish lineages remains to be determined.

Evolution of gene expression after whole-genome duplication: New insights from the spotted gar genome.

Whole-genome duplications (WGDs) are important evolutionary events. Our understanding of underlying mechanisms, including the evolution of duplicated genes after WGD, however, remains incomplete. Teleost fish experienced a common WGD (teleost-specific genome duplication, or TGD) followed by a dramatic adaptive radiation leading to more than half of all vertebrate species. The analysis of gene expression patterns following TGD at the genome level has been limited by the lack of suitable genomic resources. The recent concomitant release of the genome sequence of spotted gar (a representative of holosteans, the closest-related lineage of teleosts that lacks the TGD) and the tissue-specific gene expression repertoires of over 20 holostean and teleostean fish species, including spotted gar, zebrafish, and medaka (the PhyloFish project), offers a unique opportunity to study the evolution of gene expression following TGD in teleosts. We show that most TGD duplicates gained their current status (loss of one duplicate gene or retention of both duplicates) relatively rapidly after TGD (i.e., prior to the divergence of medaka and zebrafish lineages). The loss of one duplicate is the most common fate after TGD with a probability of approximately 80%. In addition, the fate of duplicate genes after TGD, including subfunctionalization, neofunctionalization, or retention of two "similar" copies occurred not only before but also after the divergence of species tested, in consistency with a role of the TGD in speciation and/or evolution of gene function. Finally, we report novel cases of TGD ohnolog subfunctionalization and neofunctionalization that further illustrate the importance of these processes.

Fluorescent transgenic zebrafish as a biosensor for growth-related effects of methyl parathion

Transgenic fish models are potential alternative subjects in toxicological studies, since they can provide in vivo information on the deleterious effects of different substances. Here, we used a transgenic zebrafish (Danio rerio) lineage, which expresses a destabilized fluorescent protein (DsRED) driven by the myosin light chain promoter (Mylz2), in order to propose a new research tool for environmental biomonitoring. For validating the MYO-RED lineage, we exposed fish to the organophosphorated pesticide methyl parathion (MP). The effect of MP on fish growth was assessed by evaluating weight, length, condition factor and muscle fiber diameter. All factors suffered reduction at both tested concentrations (0.13μM and 13μM of MP). Similarly, fluorescence intensity decreased in a concentration-dependent manner, suggesting muscle protein catabolism. However, DsRED gene expression lowered only at the higher MP concentration. Results indicate that the MYO-RED transgenic zebrafish is an interesting model for detecting the growth-related effects of pollutants. Destabilized proteins such as reporter genes are apparently sensitive biomarkers, since effects were observed even at the lower, environmentally acceptable concentration. Therefore, this transgenic fish is a promising candidate model for sensitive, fast, and easy environmental monitoring.

Comparative Analysis of Teleost Genome Sequences Reveals an Ancient Intron Size Expansion in the Zebrafish Lineage

We have developed a bioinformatics pipeline for the comparative evolutionary analysis of Ensembl genomes and have used it to analyze the introns of the five available teleost fish genomes. We show our pipeline to be a powerful tool for revealing variation between genomes that may otherwise be overlooked with simple summary statistics. We identify that the zebrafish, Danio rerio, has an unusual distribution of intron sizes, with a greater number of larger introns in general and a notable peak in the frequency of introns of approximately 500 to 2,000 bp compared with the monotonically decreasing frequency distributions of the other fish. We determine that 47% of D. rerio introns are composed of repetitive sequences, although the remainder, over 331 Mb, is not. Because repetitive elements may be the origin of the majority of all noncoding DNA, it is likely that the remaining D. rerio intronic sequence has an ancient repetitive origin and has since accumulated so many mutations that it can no longer be recognized as such. To study such an ancient expansion of repeats in the Danio, lineage will require further comparative analysis of fish genomes incorporating a broader distribution of teleost lineages.

GH overexpression causes muscle hypertrophy independent from local IGF-I in a zebrafish transgenic model

The aim of the present study was to analyse the morphology of white skeletal muscle in males and females from the GH-transgenic zebrafish (Danio rerio) lineage F0104, comparing the expression of genes related to the somatotrophic axis and myogenesis. Histological analysis demonstrated that transgenic fish presented enhanced muscle hypertrophy when compared to non-transgenic fish, with transgenic females being more hypertrophic than transgenic males. The expression of genes related to muscle growth revealed that transgenic hypertrophy is independent from local induction of insulin-like growth factor 1 gene (igf1). In addition, transgenic males exhibited significant induction of myogenin gene (myog) expression, indicating that myog may mediate hypertrophic growth in zebrafish males overexpressing GH. Induction of the α-actin gene (acta1) in males, independently from transgenesis, also was observed. There were no significant differences in total protein content from the muscle. Our results show that muscle hypertrophy is independent from muscle igf1, and is likely to be a direct effect of excess circulating GH and/or IGF1 in this transgenic zebrafish lineage.

SOCS1 and SOCS3 are the main negative modulators of the somatotrophic axis in liver of homozygous GH-transgenic zebrafish ( Danio rerio)

Homozygote individuals (HO) of the GH-transgenic zebrafish lineage (F0104), despite expressing double the amount of growth hormone (GH) in relation to the hemizygote (HE) individuals, presented smaller growth in relation to the last, and similar to the non-transgenic (NT) group. Through the analysis of the expression of genes of the somatotrophic axis in the livers of HO and NT individuals, it was verified that GHR, JAK2 and STAT5.1 did not present significant differences among the analyzed genotypes (NT and HO). However, in the IGF-I gene expression, an accentuated decrease was observed in group HO ( p < 0.01), suggesting a resistance effect to excess GH. This resistance could be related to the insufficient amount of energy for supporting the accelerated metabolic demand caused by excess circulating GH. Analysis of the genes involved in the regulation of GH signalization by dephosphorylation (PTP-H1 and PTP-1B) did not show any significant alteration when comparing groups HO and NT. However, the analysis of the SOCS1 and SOCS3 genes showed an induction in homozygotes of 2.5 times ( p < 0.01) and 4.3 times ( p < 0.05), respectively, in relation to non-transgenics. The results of the present work demonstrate that, in homozygotes, GH signaling is reduced by the action of the SOCS1 and SOCS3 proteins.

Japanese medaka Hox paralog group 2: insights into the evolution of Hox PG2 gene composition and expression in the Osteichthyes

Hox paralog group 2 (PG2) genes function to specify the development of the hindbrain and pharyngeal arch-derived structures in the Osteichthyes. In this article, we describe the cDNA cloning and embryonic expression analysis of Japanese medaka (Oryzias latipes) Hox PG2 genes. We show that there are only two functional canonical Hox genes, hoxa2a and b2a, and that a previously identified hoxa2b gene is a transcribed pseudogene, psihoxa2b. The functional genes, hoxa2a and b2a, were expressed in developing rhombomeres and pharyngeal arches in a manner that was relatively well conserved compared with zebrafish (Danio rerio) but differed significantly from orthologous striped bass (Morone saxatilis) and Nile tilapia (Oreochromis niloticus) genes, which, we suggest, may be owing to effects of post-genome duplication loss of a Hox PG2 gene in the medaka and zebrafish lineages. psihoxa2b was expressed at readily detectable levels in several noncanonical Hox expression domains, including the ventral aspect of the neural tube, the pectoral fin buds and caudal-most region of the embryonic trunk, indicative that regulatory control elements needed for spatio-temporal expression have diverged from their ancestral counterparts. Comparative expression analyses showed medaka hoxa2a and b2a expression in the 2nd pharyngeal arch (PA2) beyond the onset of chondrogenesis, which, according to previous hypotheses, suggests these genes function redundantly as selector genes of PA2 identity. We conclude that Hox PG2 gene composition and expression have diverged significantly during osteichthyan evolution and that this divergence in teleosts may be related to lineage-dependent differential gene loss following an actinopterygian-specific whole genome duplication.

Fish specific duplication of Dmrt2: Characterization of zebrafish Dmrt2b

The protein families with a conserved zinc finger-like DNA binding DM domain are putative transcription factors related to the sexual regulator Dsx of Drosophila and Mab-3 of Caenorhabditis elegans. Although several members have been cloned, there are still other members need to be identified, and origin and evolution of the gene family also remain unclear. We report here cloning, expression and synteny analysis of a duplicated copy of zebrafish Dmrt2a gene, the Dmrt2b which is fish specific, whereas Dmrt2a exists in all vertebrates. During embryogenesis, the Dmrt2b expression increased gradually from shield stage to hatching stage, and mainly localized in branchial arches from 24 hpf to 40 hpf, indicating that it has a potential role in the development of branchial arches. The duplicated Dmrt2b and Dmrt2a with structural variation and expression diversification during development revealed that a process of functional diploidization of gene function occurred during zebrafish lineage. DNA binding experiment indicated that Dmrt2b recognized similar DNA sequences to those of both DSX and MAB-3, indicating a conserved regulatory function. Synteny analysis among chromosomes containing Dmrt2a and Dmrt2b showed that zebrafish Dmrt2b and at least nine genes on chromosome 6 have respective homologues on chromosome 5 containing Dmrt2a. Further synteny search from genome information showed that Dmrt2b and its neighborhood existed only in the genome of teleosts. Dmrt2a and Dmrt2b were duplicated from the duplication event, which might be part of the third genome duplication, occurred during the evolution of ray-finned fishes, probably before the emergence of osteichthyes around 350 Myr. The duplicated Dmrt2b and Dmrt2a with structural variation and expression diversification suggested their diverse roles: Dmrt2b in specification of branchial arches while Dmrt2a in somitogenesis. These analyses undoubtedly help understanding functional divergence and evolution of the DM genes following gene duplication in fishes.

The effect of GH overexpression on GHR and IGF-I gene regulation in different genotypes of GH-transgenic zebrafish

Most biological actions of growth hormone (GH) are mediated by the insulin-like growth factor I (IGF-I) that is produced after the interaction of the hormone with a specific cell surface receptor, the GH receptor (GHR). Even though the GH excess on fish metabolism is poorly known, several species have been genetically engineered for this hormone in order to improve growth for aquaculture. In some GH-transgenic fish growth has been dramatically increased, while in others high levels of transgene expression have shown inhibition of the growth response. In this study, we used for the first time different genotypes (hemizygous and homozygous) of a GH-transgenic zebrafish ( Danio rerio) lineage as a model for studying the GH resistance induced by different GH transgene expression levels. The results obtained here demonstrated that homozygous fish did not grow as expected and have a lower condition factor, which implies a catabolic state. These findings are explained by a decreased IGF-I and GHR gene expression as a consequence of GH resistance. Together, our results demonstrated that homozygous GH-transgenic fish showed similar characteristics to the starvation-induced fish and could be an interesting model for studying the regulation of the GH/GHR/IGF-I axis in fish.

Flounder and fugu have a single lefty gene that covers the functions of lefty1 and lefty2 of zebrafish during L–R patterning

The lefty gene encodes a member of the TGF-β superfamily that regulates L–R axis formation during embryogenesis via antagonistic activity against Nodal, another TGF-β superfamily member. Both mouse and zebrafish have two lefty genes, lefty1 and lefty2. Interestingly, the expression domains of mouse and zebrafish lefty are different from one another. At present, the orthology and functional diversity of the mouse and zebrafish lefty genes are not clear. Here, we report that flounder and two fugu species, Takifugu and Tetraodon, have a single lefty gene in their genomes. In addition, we provide evidence that the mouse lefty genes were duplicated on a single chromosome but the zebrafish lefty genes arose from a whole-genome duplication that occurred early in the divergence of ray-finned fishes. These independent origins likely explain the difference in the expression domains of the mouse and zebrafish lefty gene pairs. Furthermore, we found that the duplication corresponding to the zebrafish lefty2 gene was lost from the fugu genome, suggesting that loss of lefty2 in the fugu/flounder lineage occurred after its divergence from the zebrafish lineage. During L–R patterning, the single lefty gene of flounder covers two expression domains, the left side of the dorsal diencephalon and the left LPM, which are regulated separately by lefty1 and lefty2 in zebrafish. We infer that the lefty genes of the ray-finned fishes and mammals underwent independent gene duplication events that resulted in independent regulation of lefty expression.

Expression of Dlx genes during the development of the zebrafish pharyngeal dentition: evolutionary implications

In order to investigate similarities and differences in genetic control of development among teeth within and between species, we determined the expression pattern of all eight Dlx genes of the zebrafish during development of the pharyngeal dentition and compared these data with that reported for mouse molar tooth development. We found that (i) dlx1a and dlx6a are not expressed in teeth, in contrast to their murine orthologs, Dlx1 and Dlx6; (ii) the expression of the six other zebrafish Dlx genes overlaps in time and space, particularly during early morphogenesis; (iii) teeth in different locations and generations within the zebrafish dentition differ in the number of genes expressed; (iv) expression similarities and differences between zebrafish Dlx genes do not clearly follow phylogenetic and linkage relationships; and (v) similarities and differences exist in the expression of zebrafish and mouse Dlx orthologs. Taken together, these results indicate that the Dlx gene family, despite having been involved in vertebrate tooth development for over 400 million years, has undergone extensive diversification of expression of individual genes both within and between dentitions. The latter type of difference may reflect the highly specialized dentition of the mouse relative to that of the zebrafish, and/or genome duplication in the zebrafish lineage facilitating a redistribution of Dlx gene function during odontogenesis.

Gene duplication, gene loss and evolution of expression domains in the vertebrate nuclear receptor NR5A (Ftz-F1) family

Fushi tarazu factor 1 (Ftz-F1, NR5A) is a zinc-finger transcription factor that belongs to the nuclear receptor superfamily and regulates genes that are involved in sterol and steroid metabolism in gonads, adrenals, liver and other tissues. To understand the evolutionary origins and developmental genetic relationships of the Ftz-F1 genes, we have cloned four homologous Ftz-f1 genes in zebrafish, called ff1a, ff1b, ff1c and ff1d. These four genes have different temporal and spatial expression patterns during development, indicating that they have distinct mechanisms of genetic regulation. Among them, the ff1a expression pattern is similar to mammalian Nr5a2, while the ff1b pattern is similar to that of mammalian Nr5a1. Genetic mapping experiments show that these four ff1 genes are located on chromosome segments conserved between the zebrafish and human genomes, indicating a common ancestral origin. Phylogenetic and conserved synteny analysis show that ff1a is the orthologue of NR5A2, and that ff1b and ff1d genes are co-orthologues of NR5A1 that arose by a gene-duplication event, probably a whole-genome duplication, in the ray-fin lineage, and each gene is located next to an NR6A1 co-orthologue as in humans, showing that the tandem duplication occurred before the divergence of human and zebrafish lineages. ff1c does not have a mammalian counterpart. Thus we have characterized the phylogenetic relationships, expression patterns and chromosomal locations of these Ftz-F1 genes, and have demonstrated their identities as NR5A genes in relation to the orthologous genes in other species.

Molecular cloning, gene expression and characterization of the third estrogen receptor of the Nile tilapia, Oreochromis niloticus

Estrogens are essential for many reproductive and non-reproductive functions. In teleosts, it is well-known that several subtypes of estrogen receptors are required for the precise action of estrogens. Present study describes the cloning of the third estrogen receptor, ER- beta2, from the Nile tilapia by EST sequencing coupled microarray. The cloned ER-beta2 showed 77.7% amino acid identity with the reported Atlantic croaker ER-beta. Three ERs, ER-alpha, ER-beta1 and ER-beta2, from the fugu genome were also isolated to analyze their gene structures. Comparison of the intron/exon boundaries and exon numbers of fugu, tilapia, rainbow trout and zebrafish, and phylogenetic analysis of 63 ER sequences revealed that ER-beta probably underwent two successive lineage-specific duplications in teleost. The former took place only in zebrafish lineage, and the latter took place in advanced teleosts without the zebrafish lineage, whereas no duplication of the ER-alpha gene has been detected. Tissue distribution analysis by RT-PCR revealed that tilapia ER-alpha and ER-beta1 were expressed ubiquitously, whereas ER-beta2 is expressed only in the pituitary, liver, intestine, kidney and gonads, with the highest expression in the testis and the lowest level in the ovary. Northern blot analysis detected a single transcript of about 3.4 kb in the testis but not in the ovary mRNAs. In transient transfection assays using human embryonic kidney 293 (HEK293) cells, tilapia ER-beta2 showed estrodiol-17beta dependent transactivation.

The ontogeny and homology of the Weberian apparatus in the zebrafish Danio rerio (Ostariophysi: Cypriniformes)

The ontogeny of the Weberian apparatus was examined in the zebrafish, Danio rerio, using both cleared and stained specimens and histology. Over 300 individuals from four independent zebrafish lineages, ranging in size from 3 to 28 mm TL, were examined for this study. Results provide a basic understanding of the development of the Weberian apparatus in the wild-type zebrafish. Our results, in conjunction with those already published, point to substantial variation in the development of the Weberian apparatus among otophysans and new interpretations of the homology of certain ossicles (e.g. tripus and claustrum). Hypotheses of homology among various Weberian ossicles are considered and represent an important step in understanding the evolution of sound transmission in ostariophysan fishes. © 2004 The Linnean Society of London, Zoological Journal of the Linnean Society, 2004, 140, 241–254.

Sequence relationships and expression patterns of zebrafish zic2 and zic5 genes

The zinc finger motif forms a DNA binding domain that is found in a wide variety of proteins. Among them, the members of the zic gene family are highly conserved throughout metazoans. We report here the isolation of two new members of this gene family in zebrafish, zic2.2 and zic5, isolated during random screening for tissue-specific genes. Zic2.2 is closely related to the previously reported zic2 gene, which we propose to rename zic2.1; these two genes form a subfamily with other vertebrate zic2 genes. We compare here the expression patterns of zic2.1, zic2.2, and zic5. All three genes showed dynamic expression patterns starting after the initiation of zygotic transcription, predominantly in the developing neural tube. Compared to zic2.1, zic2.2 was expressed in a similar but distinct manner during early development, particularly in the retina and the forming somites. A zic2.2 ortholog has not been identified in other vertebrate species, suggesting that the zic2.1/ zic2.2 pair resulted from a genome duplication event during the evolution of the zebrafish lineage.

Gene Duplication of Zebrafish JAK2 Homologs Is Accompanied by Divergent Embryonic Expression Patterns: Only jak2a Is Expressed During Erythropoiesis

Members of the JAK family of protein tyrosine kinase (PTK) proteins are required for the transmission of signals from a variety of cell surface receptors, particularly those of the cytokine receptor family. JAK function has been implicated in hematopoiesis and regulation of the immune system, and recent data suggest that the vertebrate JAK2 gene may play a role in leukemia. We have isolated and characterized jak cDNAs from the zebrafish Danio rerio. The zebrafish genome possesses 2 jak2 genes that occupy paralogous chromosome segments in the zebrafish genome, and these segments conserve syntenic relationships with orthologous genes in mammalian genomes, suggesting an ancient duplication in the zebrafish lineage. The jak2a gene is expressed at high levels in erythroid precursors of primitive and definitive waves and at a lower level in early central nervous system and developing fin buds. jak2b is expressed in the developing lens and nephritic ducts, but not in hematopoietic tissue. The expression of jak2a was examined in hematopoietic mutants and found to be disrupted in cloche and spadetail, suggesting an early role in hematopoiesis. Taken together with recent gene knockout data in the mouse, we suggest that jak2a may be functionally equivalent to mammalian Jak2, with a role in early erythropoiesis.

Zebrafish contains two Pax6 genes involved in eye development

The Pax6 genes of both vertebrates and invertebrates are expressed in the developing eye and in the central nervous system. These genes encode transcription factors with two DNA-binding domains, an N-terminal paired domain and a homeodomain separated by a flexible linker region. Ectopic eye structures are obtained upon targeted expression of Drosophila, squid, ascidian or mouse Pax6 genes in various imaginal disc primordia of Drosophila. We have previously cloned a Pax6 cDNA from zebrafish. Here we report the cloning of a novel Pax6 homolog from zebrafish denoted Pax6.2. The coding sequences of the two genes show 82% identity whereas the deduced amino acid sequences are 95% identical with complete conservation of the paired- and homeodomains. The embryonic expression patterns of Pax6.1 and Pax6.2 reveal both overlapping and discrete expression domains suggesting a division of labor between these two very similar gene products during development of brain and eye structures. Both Pax6.1 and Pax6.2 can act as transcriptional activators with Pax6.2 being more efficient than Pax6.1. Both Pax6.1 and Pax6.2 are able to induce ectopic eyes in Drosophila, while Pax2 is not, suggesting that eye induction is not a general feature of Pax family genes but a distinct characteristic of Pax6 and its direct homologs. Attempts to detect Pax6. 2 homologs in chick, mice or humans proved unsuccessful suggesting that this gene either was lost during evolution of higher vertebrates or, more likely, arose as part of a larger scale duplication of chromosome segments occurring in the zebrafish lineage.

Cloning of a neuropeptide Y/peptide YY receptor from the atlantic cod: the Yb receptor

Neuropeptide Y (NPY) belongs to a family of structurally related neuroendocrine peptides that bind to G protein-coupled receptors. Five different receptor subtypes have recently been cloned in mammals and we have found another three receptor genes in the zebrafish, called zYa, zYb, and zYc, that appear to be distinct subtypes as deduced from their widely different sequences. To elucidate the evolutionary relationships between the mammalian and zebrafish receptors, we have used the zebrafish probes to isolate genomic clones from another teleost fish, the Atlantic cod, Gadus morhua. We present here the sequence of the cod Yb gene, whose deduced protein sequence is equally identical to the zebrafish Yb (69%) and Yc proteins (66%). The two zebrafish receptors are 76% identical to each other, suggesting that they arose by gene duplication in the zebrafish lineage after divergence from the cod lineage. The five cloned mammalian NPY-family receptors and the three cloned zebrafish NPY receptors indicate that this is the largest receptor family among all peptide receptors that belong to the superfamily of G protein-coupled receptors.