DANA elements: a family of composite, tRNA-derived short interspersed DNA elements associated with mutational activities in zebrafish (Danio rerio)

Application of Repetitive Sequences in Fish Cell Depletion as a Target for the CRISPR/Cas9 System

Specific cell depletion is a common means to study the physiological function of cell lineages and tissue regeneration. However, 100% depletion is difficult to achieve with existing cell depletion strategies. With the increasing maturity of CRISPR/Cas9 technology, it is increasingly used for the depletion of various cells. However, even with this technology, it is difficult to complete the depletion of specific gene knockout cells. For this reason, cell depletion with the use of repetitive sequences as the target of CRISPR/Cas9 was explored using zebrafish. All cells were used as the target cells for the first set of experiments. The results showed that injection of a mixture of DANA-gRNA and Cas9 mRNA into zygotes resulted in substantial cell apoptosis. Cells are almost invisible in the embryonic animal pole during the dome stage. The activities of the caspase-3 and caspase-9 proteins and the mRNA level of the P53 gene were significantly increased. Then, primordial germ cells (PGCs) in embryos were used as the target cells in subsequent experiments. To specifically knock out PGCs, we injected the mix of DANA-gRNA, pkop: Cas9 plasmid (the kop promotor allows Cas9 expression only in PGCs), and eGFP-nos3'UTR mRNA into zebrafish fertilized eggs. The results revealed that the activity of the caspase-3 protein was significantly increased, and the mRNA levels of P53, ku70, and ku80 were significantly upregulated, while the number of PGCs decreased gradually. Few PGCs labeled with GFP could be seen 20 h post-fertilization (hpf), and no PGCs could be seen at the germinal ridge 24 hpf. Therefore, the combination of CRISPR/Cas9 technology and repetitive sequences can achieve efficient cell depletion regardless of whether there is generalized expression or expression in specific cells. These results indicate that it is feasible to eliminate cells by using repeat sequences as CRISPR/Cas9 system target sites.

PM2.5-induced extensive DNA methylation changes in the heart of zebrafish embryos and the protective effect of folic acid

We previously found that folic acid (FA) attenuated cardiac defects in zebrafish embryos exposed to extractable organic matter (EOM) from PM2.5, but the underlining mechanisms remain to be elucidated. Since DNA methylation is crucial to cardiac development, we hypothesized that EOM-induced aberrant DNA methylation changes could be diminished by FA supplementation. In this study, zebrafish embryos were exposed to EOM in the absence or presence of FA. Genomic-wide DNA methylation analysis identified both DNA hypo- and hyper-methylation changes in CCGG sites in zebrafish embryos exposed to EOM, which were attenuated by FA supplementation. We identified a total of 316 genes with extensive DNA methylation changes in EOM samples but little or no DNA methylation changes in EOM plus FA samples. The genes were involved in critical cellular processes and signaling pathways important for embryo development. In addition, the EOM-decreased SAM/SAH ratio was counteracted by FA supplementation. Furthermore, FA attenuated the EOM-induced changes in the expression of genes involved in the regulation of DNA methylation and in folate biosynthesis. In conclusion, our data suggest that FA supplementation protected zebrafish embryos from the cardiac developmental toxicity of PM2.5 by alleviating EOM-induced DNA methylation changes.

Lsh Is Essential for Maintaining Global DNA Methylation Levels in Amphibia and Fish and Interacts Directly with Dnmt1

Eukaryotic genomes are methylated at cytosine bases in the context of CpG dinucleotides, a pattern which is maintained through cell division by the DNA methyltransferase Dnmt1. Dramatic methylation losses are observed in plant and mouse cells lacking Lsh (lymphoid specific helicase), predominantly at repetitive sequences and gene promoters. However, the mechanism by which Lsh contributes to the maintenance of DNA methylation is unknown. Here we show that DNA methylation is lost in Lsh depleted frog and fish embryos, both of which exhibit developmental delay. Additionally, we show that both Lsh and Dnmt1 are associated with chromatin and that Lsh knockdown leads to a decreased Dnmt1-chromatin association. Coimmunoprecipitation experiments reveal that Lsh and Dnmt1 are found in the same protein complex, and pulldowns show this interaction is direct. Our data indicate that Lsh is usually diffuse in the nucleus but can be recruited to heterochromatin in a HP1α-dependent manner. These data together (a) show that the role of Lsh in DNA methylation is conserved in plants, amphibian, fish, and mice and (b) support a model in which Lsh contributes to Dnmt1 binding to chromatin, explaining how its loss can potentially lead to perturbations in DNA methylation maintenance.

Decrease in cytosine methylation at CpG island shores and increase in DNA fragmentation during zebrafish aging

Age-related changes in DNA methylation have been demonstrated in mammals, but it remains unclear as to the generality of this phenomenon in vertebrates, which is a criterion for the fundamental cause of senescence. Here we showed that the zebrafish genome gradually and clearly lost methylcytosine in somatic cells, but not in male germ cells during aging, and that age-dependent hypomethylation preferentially occurred at a particular domain called the CpG island shore, which is associated with vertebrates' genes and has been shown to be hypomethylated in humans with age. We also found that two CpG island shores hypomethylated in zebrafish oocytes were de novo methylated in fertilized eggs, which suggests that the zebrafish epigenome is reset upon fertilization, enabling new generations to restart with a heavily methylated genome. Furthermore, we observed an increase in cleavage of the zebrafish genome to an oligonucleosome length in somatic cells from the age of 12 months, which is suggestive of an elevated rate of apoptosis in the senescent stage.

Inverse PCR-based method for isolating novel SINEs from genome

Short interspersed elements (SINEs) are moderately repetitive DNA sequences in eukaryotic genomes. Although eukaryotic genomes contain numerous SINEs copy, it is very difficult and laborious to isolate and identify them by the reported methods. In this study, the inverse PCR was successfully applied to isolate SINEs from Opsariichthys bidens genome in Eastern Asian Cyprinid. A group of SINEs derived from tRNA(Ala) molecular had been identified, which were named Opsar according to Opsariichthys. SINEs characteristics were exhibited in Opsar, which contained a tRNA(Ala)-derived region at the 5' end, a tRNA-unrelated region, and AT-rich region at the 3' end. The tRNA-derived region of Opsar shared 76 % sequence similarity with tRNA(Ala) gene. This result indicated that Opsar could derive from the inactive or pseudogene of tRNA(Ala). The reliability of method was tested by obtaining C-SINE, Ct-SINE, and M-SINEs from Ctenopharyngodon idellus, Megalobrama amblycephala, and Cyprinus carpio genomes. This method is simpler than the previously reported, which successfully omitted many steps, such as preparation of probes, construction of genomic libraries, and hybridization.

RNA-Mediated Gene Duplication and Retroposons: Retrogenes, LINEs, SINEs, and Sequence Specificity

A substantial number of “retrogenes” that are derived from the mRNA of various intron-containing genes have been reported. A class of mammalian retroposons, long interspersed element-1 (LINE1, L1), has been shown to be involved in the reverse transcription of retrogenes (or processed pseudogenes) and non-autonomous short interspersed elements (SINEs). The 3′-end sequences of various SINEs originated from a corresponding LINE. As the 3′-untranslated regions of several LINEs are essential for retroposition, these LINEs presumably require “stringent” recognition of the 3′-end sequence of the RNA template. However, the 3′-ends of mammalian L1s do not exhibit any similarity to SINEs, except for the presence of 3′-poly(A) repeats. Since the 3′-poly(A) repeats of L1 and Alu SINE are critical for their retroposition, L1 probably recognizes the poly(A) repeats, thereby mobilizing not only Alu SINE but also cytosolic mRNA. Many flowering plants only harbor L1-clade LINEs and a significant number of SINEs with poly(A) repeats, but no homology to the LINEs. Moreover, processed pseudogenes have also been found in flowering plants. I propose that the ancestral L1-clade LINE in the common ancestor of green plants may have recognized a specific RNA template, with stringent recognition then becoming relaxed during the course of plant evolution.

Analysis of the use of codon pairs in the HE gene of the ISA virus shows a correlation between bias in HPR codon-pair use and mortality rates caused by the virus

BACKGROUND

Segment 6 of the ISA virus codes for hemoagglutinin-esterase (HE). This segment is highly variable, with more than 26 variants identified. The major variation is observed in what is called the high polymorphism region (HPR). The role of the different HPR zones in the viral cycle or evolution remains unknown. However viruses that present the HPR0 are avirulent, while viruses with important deletions in this region have been responsible for outbreaks with high mortality rates. In this work, using bioinformatic tools, we examined the influence of different HPRs on the adaptation of HE genes to the host translational machinery and the relationship to observed virulence.

METHODS

Translational efficiency of HE genes and their HPR were estimated analyzing codon-pair bias (CPB), adaptation to host codon use (codon adaptation index-CAI) and the adaptation to available tRNAs (tAI). These values were correlated with reported mortality for the respective ISA virus and the ΔG of RNA folding. tRNA abundance was inferred from tRNA gene numbers identified in the Salmo salar genome using tRNAScan-SE. Statistical correlation between data was performed using a non-parametric test.

RESULTS

We found that HPR0 contains zones with codon pairs of low frequency and low availability of tRNA with respect to salmon codon-pair usage, suggesting that HPR modifies HE translational efficiency. Although calculating tAI was impossible because one third of tRNAs (~60.000) were tRNA-ala, translational efficiency measured by CPB shows that as HPR size increases, the CPB value of the HE gene decreases (P = 2x10⁻⁷, ρ = -0.675, n = 63) and that these values correlate positively with the mortality rates caused by the virus (ρ = 0.829, P = 2x10⁻⁷, n = 11). The mortality associated with different virus isolates or their corresponding HPR sizes were not related with the ΔG of HPR RNA folding, suggesting that the secondary structure of HPR RNA does not modify virulence.

CONCLUSIONS

Our results suggest that HPR size affects the efficiency of gene translation, which modulates the virulence of the virus by a mechanism similar to that observed in production of live attenuated vaccines through deoptimization of codon-pair usage.

SINEBase: a database and tool for SINE analysis

SINEBase (http://sines.eimb.ru) integrates the revisited body of knowledge about short interspersed elements (SINEs). A set of formal definitions concerning SINEs was introduced. All available sequence data were screened through these definitions and the genetic elements misidentified as SINEs were discarded. As a result, 175 SINE families have been recognized in animals, flowering plants and green algae. These families were classified by the modular structure of their nucleotide sequences and the frequencies of different patterns were evaluated. These data formed the basis for the database of SINEs. The SINEBase website can be used in two ways: first, to explore the database of SINE families, and second, to analyse candidate SINE sequences using specifically developed tools. This article presents an overview of the database and the process of SINE identification and analysis.

Polypteridae (Actinopterygii: Cladistia) and DANA-SINEs insertions

SINE sequences are interspersed throughout virtually all eukaryotic genomes and greatly outnumber the other repetitive elements. These sequences are of increasing interest for phylogenetic studies because of their diagnostic power for establishing common ancestry among taxa, once properly characterized. We identified and characterized a peculiar family of composite tRNA-derived short interspersed SINEs, DANA-SINEs, associated with mutational activities in Danio rerio, in a group of species belonging to one of the most basal bony fish families, the Polypteridae, in order to investigate their own inner specific phylogenetic relationships. DANA sequences were identified, sequenced and then localized, by means of fluorescent in situ hybridization (FISH), in six Polypteridae species (Polypterus delhezi, P. ornatipinnis, P. palmas, P. buettikoferi P. senegalus and Erpetoichthys calabaricus) After cloning, the sequences obtained were aligned for phylogenetic analysis, comparing them with three Dipnoan lungfish species (Protopterus annectens, P. aethiopicus, Lepidosiren paradoxa), and Lethenteron reissneri (Petromyzontidae)was used as outgroup. The obtained overlapping MP, ML and NJ tree clustered together the species belonging to the two taxonomically different Osteichthyans groups: the Polypteridae, by one side, and the Protopteridae by the other, with the monotypic genus Erpetoichthys more distantly related to the Polypterus genus comprising three distinct groups: P. palmas and P. buettikoferi, P. delhezi and P. ornatipinnis and P. senegalus. In situ hybridization with DANA probes marked along the whole chromosome arms in the metaphases of all the Polypteridae species examined.

The role of transposable elements in the evolution of non-mammalian vertebrates and invertebrates

Background

Transposable elements (TEs) have played an important role in the diversification and enrichment of mammalian transcriptomes through various mechanisms such as exonization and intronization (the birth of new exons/introns from previously intronic/exonic sequences, respectively), and insertion into first and last exons. However, no extensive analysis has compared the effects of TEs on the transcriptomes of mammals, non-mammalian vertebrates and invertebrates.

Results

We analyzed the influence of TEs on the transcriptomes of five species, three invertebrates and two non-mammalian vertebrates. Compared to previously analyzed mammals, there were lower levels of TE introduction into introns, significantly lower numbers of exonizations originating from TEs and a lower percentage of TE insertion within the first and last exons. Although the transcriptomes of vertebrates exhibit significant levels of exonization of TEs, only anecdotal cases were found in invertebrates. In vertebrates, as in mammals, the exonized TEs are mostly alternatively spliced, indicating that selective pressure maintains the original mRNA product generated from such genes.

Conclusions

Exonization of TEs is widespread in mammals, less so in non-mammalian vertebrates, and very low in invertebrates. We assume that the exonization process depends on the length of introns. Vertebrates, unlike invertebrates, are characterized by long introns and short internal exons. Our results suggest that there is a direct link between the length of introns and exonization of TEs and that this process became more prevalent following the appearance of mammals.

Loss of Dnmt1 catalytic activity reveals multiple roles for DNA methylation during pancreas development and regeneration

Developmental mechanisms regulating gene expression and the stable acquisition of cell fate direct cytodifferentiation during organogenesis. Moreover, it is likely that such mechanisms could be exploited to repair or regenerate damaged organs. DNA methyltransferases (Dnmts) are enzymes critical for epigenetic regulation, and are used in concert with histone methylation and acetylation to regulate gene expression and maintain genomic integrity and chromosome structure. We carried out two forward genetic screens for regulators of endodermal organ development. In the first, we screened for altered morphology of developing digestive organs, while in the second we screed for the lack of terminally differentiated cell types in the pancreas and liver. From these screens, we identified two mutant alleles of zebrafish dnmt1. Both lesions are predicted to eliminate dnmt1 function; one is a missense mutation in the catalytic domain and the other is a nonsense mutation that eliminates the catalytic domain. In zebrafish dnmt1 mutants, the pancreas and liver form normally, but begin to degenerate after 84 h post fertilization (hpf). Acinar cells are nearly abolished through apoptosis by 100 hpf, though neither DNA replication, nor entry into mitosis is halted in the absence of detectable Dnmt1. However, endocrine cells and ducts are largely spared. Surprisingly, dnmt1 mutants and dnmt1 morpholino-injected larvae show increased capacity for pancreatic beta cell regeneration in an inducible model of pancreatic beta cell ablation. Thus, our data suggest that Dnmt1 is dispensable for pancreatic duct or endocrine cell formation, but not for acinar cell survival. In addition, Dnmt1 may influence the differentiation of pancreatic beta cell progenitors or the reprogramming of cells toward the pancreatic beta cell fate.

Transposon tools hopping in vertebrates

In the past decade, tools derived from DNA transposons have made major contributions to vertebrate genetic studies from gene delivery to gene discovery. Multiple, highly complementary systems have been developed, and many more are in the pipeline. Judging which DNA transposon element will work the best in diverse uses from zebrafish genetic manipulation to human gene therapy is currently a complex task. We have summarized the major transposon vector systems active in vertebrates, comparing and contrasting known critical biochemical and in vivo properties, for future tool design and new genetic applications.

Inactivation of zebrafish mrf4 leads to myofibril misalignment and motor axon growth disorganization

Mrf4 is a basic helix-loop-helix (bHLH) transcription factor associated with myogenesis. Two mrf4 transcripts, mrf4_tv1 and mrf4_tv2, were identified in zebrafish generated by alternative splicing. To study their biological functions, we separately injected the Mrf4-morpholinos, including MO1 (mrf4_tv1:mrf4_tv2 knockdown), MO2+MO3 (mrf4_tv1:mrf4_tv2 knockdown), MO3 (mrf4_tv1 knockdown), and MO4 (mrf4_tv2 knockdown), into zebrafish embryos to observe mrf4 gene knockdown phenotypes. No phenotypic abnormalities were observed following injection with 0.5 ng of MO1 but those injected with 4.5, 9, or 13.5 ng displayed curved-body phenotypes, such as indistinct somite boundaries, and a lack of uniformly sized cell blocks. Similar results were also observed in the (MO2+MO3)-, MO3-, and MO4-injected groups. To further investigate the molecular mechanisms that lead to curved-body phenotypes, we stained embryos with alpha-bungrotoxin and specific monoclonal antibodies F59, Znp1, and Zn5 to detect morphological changes in acetyl-choline receptor (AChR) clusters, muscle fibers, common path of the primary neurons, and secondary neurons axonal projections, respectively. Our results show that the muscle fibers of mrf4_(tv1:tv2)-morphant aligned disorderly and lost their integrity and attachment, while the defects became milder in either mrf4_tv1-morphant or mrf4_tv2-morphant. On the other hand, reduced axonal projections and AChR clusters were found in both mrf4_tv2-morphant and mrf4_(tv1:tv2)-morphant but distributed normally in the mrf4_tv1-morphant. We conclude that Mrf4_tv2 is involved in alignment of muscle fibers, and Mrf4_tv1 might have cooperative function with Mrf4_tv2 in muscle fiber alignment, without affecting the muscle-nerve connection.

Retroposons of salmonoid fishes (Actinopterygii: Salmonoidei) and their evolution

Short and long retroposons, or non-LTR retrotransposons (SINEs and LINEs, respectively) are two groups of interspersed repetitive elements amplifying in the genome via RNA and cDNA-mediated reverse transcription. In this process, SINEs entirely depend on the enzymatic machinery of autonomous LINEs. The impact of retroposons on the host genome is difficult to overestimate: their sequences account for significant portion of the eukaryotic genome, while propagation of their active copies gradually reshapes it. In this way, the retropositional activity plays a role of important evolutionary factor. More than 100 LINE and nearly 100 SINE families have been described to date from the genomes of various eukaryotes, and it is salmonoid fishes (Actinopterygii: Salmonoidei) that are particularly noticeable for the diversity of transposons they host in their genomes, including two LINE and seven SINE families. Moreover, this group of ray-finned fish represents an excellent opportunity to study such a rare evolutionary phenomenon as lateral gene transfer, due to a great variety of transposons and other sequences salmons share with a blood fluke, Schistosoma japonicum (Trematoda: Strigeiformes)--a parasitic helminth infecting various vertebrates. The aim of the present review is to structure all knowledge accumulated about salmonoid retroposons by now, as well as to complement it with the new data pertaining to the distribution of some SINE families.

Mrf4 (myf6) is dynamically expressed in differentiated zebrafish skeletal muscle

Mrf4 (Myf6) is a member of the basic helix-loop-helix (bHLH) myogenic regulatory transcription factor (MRF) family, which also contains Myod, Myf5 and myogenin. Mrf4 is implicated in commitment of amniote cells to skeletal myogenesis and is also abundantly expressed in many adult muscle fibres. The specific role of Mrf4 is unclear both because mrf4 null mice are viable, suggesting redundancy with other MRFs, and because of genetic interactions at the complex mrf4/myf5 locus. We report the cloning and expression of an mrf4 gene from zebrafish, Danio rerio, which shows conservation of linkage to myf5. Mrf4 mRNA accumulates in a subset of terminally differentiated muscle fibres in parallel with myosin protein in the trunk and fin. Although most, possibly all, trunk muscle expresses mrf4, the level of mRNA is dynamically regulated. No expression is detected in muscle precursor cell populations prior to myosin accumulation. Moreover, mrf4 expression is not detected in head muscles, at least at early stages. As fish mature, mrf4 expression is pronounced in the region of slow muscle fibres.

Transposon tools and methods in zebrafish

Zebrafish is an excellent model animal to study vertebrate development by genetic approaches. Hundreds of mutations affecting various processes of development have been isolated by chemical mutagenesis and insertional mutagenesis using a pseudotyped retrovirus. However, useful transposon tools and methods had not been available in zebrafish. This is mainly because no active transposable element has been found from the zebrafish genome. Recently, efficient transgenesis, gene trap, and enhancer trap methods have been developed in zebrafish by using the Tol2 and the Sleeping Beauty transposon systems. These methods should increase the usefulness of zebrafish as a model vertebrate and facilitate the study of developmental biology, genetics, and genomics.

Identification of a gene required for de novo DNA methylation of the zebrafish no tail gene

The zebrafish no tail gene (ntl) is indispensable for tail and notochord development. We have shown previously that ntl is de novo methylated during early embryogenesis. To find the gene that de novo methylates ntl and understand the meaning of this methylation, we cloned seven genes that encode the conserved catalytic domain of methyltransferases. We found that injection of antisense morpholino oligonucleotides against one of them, termed dnmt7, into eggs significantly reduced the level of ntl methylation, although no apparent phenotype was induced by the injection. Inhibition of Dnmt7 activity did not change the level of genome-wide methylation nor did it affect de novo methylation of injected plasmid DNA, indicating that Dnmt7 specifically methylates ntl in the genome.

Comparative evolution history of SINEs in Arabidopsis thaliana and Brassica oleracea: evidence for a high rate of SINE loss

Brassica oleracea and Arabidopsis thaliana belong to the Brassicaceae(Cruciferae) family and diverged 16 to 19 million years ago. Although the genome size of B. oleracea (approximately 600 million base pairs) is more than four times that of A. thaliana (approximately 130 million base pairs), their gene content is believed to be very similar with more than 85% sequence identity in the coding region. Therefore, this important difference in genome size is likely to reflect a different rate of non-coding DNA accumulation. Transposable elements (TEs) constitute a major fraction of non-coding DNA in plant species. A different rate in TE accumulation between two closely related species can result in significant genome size variations in a short evolutionary period. Short interspersed elements (SINEs) are non-autonomous retroposons that have invaded the genome of most eukaryote species. Several SINE families are present in B. oleracea and A. thaliana and we found that two of them (called RathE1 and RathE2) are present in both species. In this study, the tempo of evolution of RathE1 and RathE2 SINE families in both species was compared. We observed that most B. oleracea RathE2 SINEs are "young" (close to the consensus sequence) and abundant while elements from this family are more degenerated and much less abundant in A. thaliana. However, the situation is different for the RathE1 SINE family for which the youngest elements are found in A. thaliana. Surprisingly, no SINE was found to occupy the same (orthologous) genomic locus in both species suggesting that either these SINE families were not amplified at a significant rate in the common ancestor of the two species or that older elements were lost and only the recent (lineage-specific) insertions remain. To test this latter hypothesis, loci containing a recently inserted SINE in the A. thaliana col-0 ecotype were selected and characterized in several other A. thaliana ecotypes. In addition to the expected SINE containing allele and the pre-integrative allele (i.e. the "empty" allele), we observed in the different ecotypes, alleles with truncated portions of the SINE (up to the complete loss of the element) and of the immediate genomic flanking sequences. The absence of SINEs in orthologous positions between B. oleracea and A. thaliana and the presence in recently diverged A. thaliana ecotypes of alleles containing severely truncated SINEs suggest a very high rate of SINE loss in these species.

SINEs and LINEs: symbionts of eukaryotic genomes with a common tail

Many SINEs and LINEs have been characterized to date, and examples of the SINE and LINE pair that have the same 3' end sequence have also increased. We report the phylogenetic relationships of nearly all known LINEs from which SINEs are derived, including a new example of a SINE/LINE pair identified in the salmon genome. We also use several biological examples to discuss the impact and significance of SINEs and LINEs in the evolution of vertebrate genomes.

Genome evolution and biodiversity in teleost fish

Teleost fish, which roughly make up half of the extant vertebrate species, exhibit an amazing level of biodiversity affecting their morphology, ecology and behaviour as well as many other aspects of their biology. This huge variability makes fish extremely attractive for the study of many biological questions, particularly of those related to evolution. New insights gained from different teleost species and sequencing projects have recently revealed several peculiar features of fish genomes that might have played a role in fish evolution and speciation. There is now substantial evidence that a round of tetraploidization/rediploidization has taken place during the early evolution of the ray-finned fish lineage, and that hundreds of duplicate pairs generated by this event have been maintained over hundreds of millions of years of evolution. Differential loss or subfunction partitioning of such gene duplicates might have been involved in the generation of fish variability. In contrast to mammalian genomes, teleost genomes also contain multiple families of active transposable elements, which might have played a role in speciation by affecting hybrid sterility and viability. Finally, the amazing diversity of sex determination systems and the plasticity of sex chromosomes observed in teleost might have been involved in both pre- and postmating reproductive isolation. Comparison of data generated by current and future genome projects as well as complementary studies in other species will allow one to approach the molecular and evolutionary mechanisms underlying genome diversity in fish, and will certainly significantly contribute to our understanding of gene evolution and function in humans and other vertebrates.

Insertional mutagenesis in zebrafish

Insertional mutagenesis is a method for identifying genes essential for a given biological process by using the integration of DNA as the mutagen, thereby facilitating the cloning of the mutated gene. The use of retrovirus-mediated insertional mutagenesis in zebrafish has led to the mutation and rapid identification of hundreds of genes required for embryonic development and cell viability and growth, revealing the diversity of gene products required for the development of this vertebrate. Here, I will review the methodology of this approach and the results to date, as well as other potential ways to use insertional mutagenesis for genetic screens.

Efficient gene delivery and gene expression in zebrafish using the Sleeping Beauty transposon

We used the Tc1/mariner family transposable element Sleeping Beauty (SB) for transgenesis and long-term expression studies in the zebrafish (Danio rerio), a popular organism for clinical disease, vertebrate patterning, and cell biology applications. SB transposase enhanced the transgenesis and expression rate sixfold (from 5 to 31%) and more than doubled the total number of tagged chromosomes over standard, plasmid injection-based transgenesis methods. Molecular analysis of these loci demonstrated a precise integration of these elements into recipient chromosomes with genetic footprints diagnostic of transposition. GFP expression from transposase-mediated integrants was Mendelian through the eighth generation. A blue-shifted GFP variant (BFP) and a red fluorescent protein (DsRed) were also useful transgenesis markers, indicating that multiple reporters are practical for use with SB in zebrafish. We showed that SB is suitable for tissue-specific transgene applications using an abbreviated gamma-crystallin GFP cassette. Finally, we describe a general utility transposon vector for chromosomal engineering and molecular genetics experiments in zebrafish. Together, these data indicate that SB is an efficient tool for transgenesis and expression in zebrafish, and that the transposon will be useful for gene expression in cell biology applications as well as an insertional mutagen for gene discovery during development.

The vasa locus in zebrafish: multiple RGG boxes from duplications

Vasa is the most extensively characterized germ cell-specific marker in metazoans. We determined the sequence of the zebrafish vasa locus, which - together with the flanking regions - is 25 kb long. It contains 26 + 1 exons, with an average length of 83 bp. The 5' end of zebrafish vasa is rich in repeats; it includes a peculiar tandem repeat with three units spanning a total of eight exons, each unit containing several tandemly arranged RNA-binding motifs (RGG boxes). The presence of various, nonmicrosatellite type repeats in vasa seems to be universal among the species studied, and could have functional importance during the evolution of the gene, due to the increase in the number of RGG boxes located in these repeats. Analysis of vasa transcripts in zebrafish identified numerous isoforms resulting from alternative splicing (seven) and polyadenylation (four). Mapping to two radiation hybrid panels strengthened the position of the gene on LG10 of zebrafish.

In silico and in situ characterization of the zebrafish (Danio rerio) gnrh3 (sGnRH) gene

BACKGROUND

Gonadotropin releasing hormone (GnRH) is responsible for stimulation of gonadotropic hormone (GtH) in the hypothalamus-pituitary-gonadal axis (HPG). The regulatory mechanisms responsible for brain specificity make the promoter attractive for in silico analysis and reporter gene studies in zebrafish (Danio rerio).

RESULTS

We have characterized a zebrafish [Trp7, Leu8] or salmon (s) GnRH variant, gnrh3. The gene includes a 1.6 Kb upstream regulatory region and displays the conserved structure of 4 exons and 3 introns, as seen in other species. An in silico defined enhancer at -976 in the zebrafish promoter, containing adjacent binding sites for Oct-1, CREB and Sp1, was predicted in 2 mammalian and 5 teleost GnRH promoters. Reporter gene studies confirmed the importance of this enhancer for cell specific expression in zebrafish. Interestingly the promoter of human GnRH-I, known as mammalian GnRH (mGnRH), was shown capable of driving cell specific reporter gene expression in transgenic zebrafish.

CONCLUSIONS

The characterized zebrafish Gnrh3 decapeptide exhibits complete homology to the Atlantic salmon (Salmo salar) GnRH-III variant. In silico analysis of mammalian and teleost GnRH promoters revealed a conserved enhancer possessing binding sites for Oct-1, CREB and Sp1. Transgenic and transient reporter gene expression in zebrafish larvae, confirmed the importance of the in silico defined zebrafish enhancer at -976. The capability of the human GnRH-I promoter of directing cell specific reporter gene expression in zebrafish supports orthology between GnRH-I and GnRH-III.

Classical and molecular cytogenetics of the zebrafish, Danio rerio (Cyprinidae, Cypriniformes): an overview

The zebrafish, Danio rerio, has recently become the model system for the genetic analysis of vertebrate development. This paper reviews the advances in zebrafish cytogenetics, obtained through classical and molecular techniques, which will lead to the assignment of specific linkage groups to specific chromosome pairs in the zebrafish genome project. Several chromosome pairs of the 50-chromosome karyotype of D. rerio were differentially stained by classical staining techniques and additional information has been obtained by molecular cytogenetics. Indeed, the analysis of constitutive heterochromatin by C-banding and base-specific fluorochrome staining had suggested a differential composition of peri- and paracentromeric constitutive heterochromatin. The chromosome mapping of distinct AT- and GC-rich zebrafish satellite DNAs by means of PRINS (Primed in situ) and multicolor FISH (Fluorescence in situ Hybridization) has confirmed this hypothesis, which therefore provided the chromosome localization of 10% of the zebrafish genome. The analysis of nucleolus organizer regions (NORs) by silver staining and by FISH with 18S rDNA has also revealed the existence of variable and inactive NORs, in addition to those on the terminal regions of the long arms of the three NOR-bearing chromosome pairs. Other multicopy genes, such as minor ribosomal genes, or multicopy repeats, such as telomere specific sequences, have now been mapped on zebrafish chromosomes. The latest advancement in zebrafish molecular cytogenetics is the chromosome mapping of single locus genes. Single-copy genes from each of the 25 genetic linkage groups are now being mapped on zebrafish chromosomes by using PAC clones.

V-SINEs: a new superfamily of vertebrate SINEs that are widespread in vertebrate genomes and retain a strongly conserved segment within each repetitive unit

We have identified a new superfamily of vertebrate short interspersed repetitive elements (SINEs), designated V-SINEs, that are widespread in fishes and frogs. Each V-SINE includes a central conserved domain preceded by a 5'-end tRNA-related region and followed by a potentially recombinogenic (TG)(n) tract, with a 3' tail derived from the 3' untranslated region (UTR) of the corresponding partner long interspersed repetitive element (LINE) that encodes a functional reverse transcriptase. The central domain is strongly conserved and is even found in SINEs in the lamprey genome, suggesting that V-SINEs might be approximately 550 Myr old or older in view of the timing of divergence of the lamprey lineage from the bony fish lineage. The central conserved domain might have been subject to some form of positive selection. Although the contemporary 3' tails of V-SINEs differ from one another, it is possible that the original 3' tail might have been replaced, via recombination, by the 3' tails of more active partner LINEs, thereby retaining retropositional activity and the ability to survive for long periods on the evolutionary time scale. It seems plausible that V-SINEs may have some function(s) that have been maintained by the coevolution of SINEs and LINEs during the evolution of vertebrates.

Use of the zebrafish (Danio rerio) to define hematopoiesis

Hematopoiesis in the vertebrate is characterized by the induction of ventral mesoderm to form hematopoietic stem cells and the eventual differentiation of these progenitors to form the peripheral blood lineages. Several genes have been implicated in the differentiation and development of hematopoietic and vascular progenitor cells, yet our understanding of the discrete steps involved in the induction of these cells from the ventral mesoderm is still incomplete. One method of delineating these processes is based on the use of lower vertebrates. The zebrafish (Danio rerio) is an especially robust vertebrate system for both isolating and characterizing genes involved in these processes. Hematopoietic mutants have been generated with defects in many of the steps of both the primitive and definitive hematopoietic programs. Cloning of the genes that underlie these mutations should yield valuable details of hematopoiesis and may have therapeutic implications for bone marrow transplantation and stem cell gene therapy.

Zebrafish acetylcholinesterase is encoded by a single gene localized on linkage group 7. Gene structure and polymorphism; molecular forms and expression pattern during development

We cloned and sequenced the acetylcholinesterase gene and cDNA of zebrafish, Danio rerio. We found a single gene (ache) located on linkage group LG7. The relative organization of ache, eng2, and shh genes is conserved between zebrafish and mammals and defines a synteny. Restriction fragment length polymorphism analysis was allowed to identify several allelic variations. We also identified two transposable elements in non-coding regions of the gene. Compared with other vertebrate acetylcholinesterase genes, ache gene contains no alternative splicing at 5' or 3' ends where only a T exon is present. The translated sequence is 60-80% identical to acetylcholinesterases of the vertebrates and exhibits an extra loop specific to teleosts. Analysis of molecular forms showed a transition, at the time of hatching, from the globular G4 form to asymmetric A12 form that becomes prominent in adults. In situ hybridization and enzymatic activity detection on whole embryos confirmed early expression of the acetylcholinesterase gene in nervous and muscular tissues. We found no butyrylcholinesterase gene or activity in Danio. These findings make zebrafish a promising model to study function of acetylcholinesterase during development and regulation of molecular forms assembly in vivo.

Birth of a retroposon: the Twin SINE family from the vector mosquito Culex pipiens may have originated from a dimeric tRNA precursor

SINEs are short interspersed repetitive elements found in many eukaryotic genomes and are believed to propagate by retroposition. Almost all SINEs reported to date have a composite structure made of a 5' tRNA-related region followed by a tRNA-unrelated region. Here, we describe a new type of tRNA-derived SINEs from the genome of the mosquito Culex pipiens. These elements, called TWINs, are approximately 220 bp long and reiterated at approximately 500 copies per haploid genome. TWINs have a unique structure compared with other tRNA-SINEs described so far. They consist of two tRNA(Arg)-related regions separated by a 39-bp spacer. Other tRNA-unrelated sequences include a 5-bp leader preceding the left tRNA-like unit and a short trailer located downstream of the right tRNA-like region. This 3' trailer is a 10-bp sequence that is ended by a TTTT motif and followed by a polyA tract of variable length. The right tRNA-like unit also contains a 16-bp sequence which is absent in the left one and appears to be located in the ancestral anticodon stem precisely at a position expected for a nuclear tRNA intron. According to this singular structure, we hypothesize that the TWIN: SINE family originated from an unprocessed polymerase III transcript containing two tRNA sequences. We suggest that some peculiar properties acquired by this dicistronic transcript, such as a polyA tail and a 3' stem-loop secondary structure, promote its retroposition by increasing its chances of being recognized by a reverse transcriptase encoded elsewhere in the C. pipiens genome.

Highly repetitive elements from Chinese bitterlings (genus Rhodeus, Cyprinidae)

We have isolated and characterized several highly repetitive DNA elements from two species of Chinese bitterlings, Rhodeus atremius suigensis and R. ocellatus ocellatus. They comprise a partly interspersed and partly tandem repetitive family of about 1.0 to 1.3 kb in length. Individual elements showed considerable length variation, but genomic Southern blotting revealed two major length groups. Their restricted presence of these elements among related species and relative copy number differences indicated rapid change of genome structure in this group of fish. The isolated elements may be useful landmarks for further chromosomal studies.

Sequence analysis of transposable elements in the sea squirt, Ciona intestinalis

A systematic search of 1 Mb of genomic sequences from the sea squirt, Ciona intestinalis, revealed the presence of six families of transposable elements. The Cigr-1 retrotransposon contains identical 245-bp long terminal repeats (LTRs) and a 3,630-bp open reading frame (ORF) encoding translation products in the same order as the domains characteristic of gypsy/Ty3-type LTR retrotransposons. The closest homologs of the reverse transcriptase domain were in gypsy elements from Drosophila and the sushi element from the pufferfish. However, the capsid-nucleocapsid region shows the clearest homology to an echinoderm element, Tgr1. Database searches also indicated two classes of non-LTR retrotransposon, named Cili-1 and Cili-2. The Cili-1 sequences show matches to regions of the ORF2 product of mammalian L1 elements. The Cili-2 sequences possess similarity to the RNaseH domain of Lian-Aa1, a mosquito non-LTR retrotransposon. The most abundant element was a short interspersed nucleotide element named Cics-1 with a copy number estimated at 40,000. Cics-1 consists of two conserved domains separated by an A-rich stretch. The 172-bp 5' domain is related to tRNA sequences, whereas the 110-bp 3' domain is unique. Cics-1 is unusual, not just in its modular structure, but also in its lack of a 3' poly(A) tail or direct flanking repeats. A second abundant element, Cimi-1, has an A+T-rich 193-bp consensus sequence and 30-bp terminal inverted repeats (TIRs) and is usually flanked by A+T-rich 2-4-bp putative target site duplications-characteristics of miniature inverted-repeat transposable elements found in plants and insects. A single 2,444-bp foldback element was found, possessing long TIRs containing an A+T-rich internal domain, an array of subrepeats, and a flanking domain at the TIR ends; this is the first example of a chordate foldback element. This study provides the first systematic characterization of the families of transposable elements in a lower chordate.

Multiple lineages of the non-LTR retrotransposon Rex1 with varying success in invading fish genomes

Rex1, together with the related BABAR: elements, represents a new family of non-long-terminal-repeat (non-LTR) retrotransposons from fish, which might be related to the CR1 clade of LINE elements. Rex1/BABAR: retrotransposons encode a reverse transcriptase and an apurinic/apyrimidinic endonuclease, which is very frequently removed by incomplete reverse transcription. Different Rex1 elements show a conserved terminal 3' untranslated region followed by oligonucleotide tandem repeats of variable size and sequence. Phylogenetic analysis revealed that Rex1 retrotransposons were frequently active during fish evolution. They formed multiple ancient lineages, which underwent several independent and recent bursts of retrotransposition and invaded fish genomes with varying success (from <5 to 500 copies per haploid genome). At least three of these ancient Rex1 lineages were detected within the genome of poeciliids. One lineage is absent from some poeciliids but underwent successive rounds of retrotransposition in others, thereby increasing its copy number from <10 to about 200. At least three ancient Rex1 lineages were also detected in the genome project fish Fugu rubripes. Rex1 distribution within one of its major lineages is discontinuous: Rex1 was found in all Acanthopterygii (common ancestor in the main teleost lineage approximately 90 MYA) and in both European and Japanese eels (divergence from the main teleost lineage about 180 MYA) but not in trout, pike, carp, and zebrafish (divergence 100-120 MYA). This might either result from frequent loss or rapid divergence of Rex1 elements specifically in some fish lineages or represent one of the very rare examples of horizontal transfer of non-LTR retrotransposons. This analysis highlights the dynamics and complexity of retrotransposon evolution and the variability of the impact of retrotransposons on vertebrate genomes.

Phylogenetic analysis of a retroposon family as represented on the human X chromosome

SINE-R elements constitute a class of retroposons derived from the long terminal repeat (LTR) of the human endogenous retrovirus HERV-K family that are present in hominoid primates and active in the human genome. In an investigation of the X chromosome, we identified twenty-five SINE-R elements with between 89.6 and 97.7% homology with the SINE-R.C2 element that is human specific, originally identified in the gene for the C2 component of complement. SINE-R.C2 and a sequence HS307 that we previously identified in a region of Xq21.3 that has a recently created homology with a 4 Mb block in Yp11.2 are amongst the group of elements that have diverged furthest from the parent HERV-K10 sequence. The sequence on the X chromosome resemble those that we previously described on chromosomes 7 and 17 and the Y chromosome, with a similar range of variation. Phylogenetic analysis from the retroposon family including those of African great apes using the neighbor-joining method suggests that the SINE-R retroposon family have evolved independently during primate evolution. Further investigation of SINE-R elements on the sex chromosomes, particularly in recently created regions of X-Y homology, may cast light on the timing of the retroposition process and its possible relevance to recent evolutionary change.

The skeletal muscle alpha-actin gene of channel catfish (Ictalurus punctatus) and its association with piscine specific SINE elements

The alpha-actin gene of channel catfish (Ictalurus punctatus) was cloned and sequenced. The gene has a similar organization and exhibited a high level of sequence similarity to those from other vertebrate animals. The upstream region of the alpha-actin gene included a TATA box, a CAAT box, three E-boxes, and a CArG box. Nested deletion segments containing these transcriptional motifs were fused to the reporter gene chloramphenicol acetyl transferase (CAT). Transfection of the clones into C2C12 cells indicated that all these motifs are required for transcriptional activities. The channel catfish alpha-actin gene is associated with two distinct short interspersed repetitive elements (SINEs). The first SINE element showed high levels of sequence similarity to the zebrafish Mermaid element, while the second SINE element is not similar to the Mermaid element except for an 8bp sequence CCCCGTGC suggesting their evolutionary linkage. However, the second SINE element appeared to co-exist with the Mermaid element in most cases and therefore was designated as the Merman element. Approximately 9000 copies and 1200 copies of the Mermaid and Merman elements exist per haploid channel catfish genome, respectively. BLAST searches indicated that both the Mermaid and the Merman elements were frequently associated with gene sequences, mostly those of aquatic animals, suggesting their evolutionary origin in association with aquatic organisms and their function in shaping the evolution of genomes in aquatic animals.

Identification of seven genes in the major histocompatibility complex class I region of the zebrafish

Physical linkage of genes whose products are involved in similar physiological pathways may have functional significance. The identification of conserved gene linkage in distantly related organisms can therefore strengthen the hypothesis of selection acting towards keeping genes on a chromosome. We used the cDNA selection technique and the polymerase chain reaction (PCR) with generic primers for the identification of new genes on the genomic clones bearing the major histocompatibility complex (Mhc) class I genes of the zebrafish (Danio rerio). We found six new genes (BING1, DAXX, TAPBP, KNSL2, TAP2B and KE6) whose orthologues are known to be linked to the Mhc class II region in humans and mice. In addition, a new zebrafish Mhc class I gene, termed Dare-UFA, was detected. By contrast, a search for the human leucocyte antigen (HLA)-linked BING3, KE3 and SACM2L genes revealed that these loci are not located on the class I clones of the zebrafish. The zebrafish class I region contains repetitive elements with similarity to the DANA, SATA and LINE repeats, as well as Tc1 transposable elements. Our findings indicate a high degree of linkage conservation between the zebrafish class I and the mammalian class II regions.

Zebrafish genetic map with 2000 microsatellite markers

The zebrafish is the first vertebrate organism used for large-scale genetic screens seeking genes critical to development. These screens have been quite successful, with more than 1800 recessive mutations discovered that speak to morphogenesis of the vertebrate embryo. The cloning of the mutant genes depends on a dense genetic map. The 2000 markers we present here, using microsatellite (CA) repeats, provides 1.2-cM average resolution. One centimorgan in zebrafish is about 0. 74 megabase, so, for many mutations, these markers are close enough to begin positional cloning by YAC walks.

Analysis of the spermine synthase gene region in Fugu rubripes, Tetraodon fluviatilis, and Danio rerio

A prerequisite to understanding the evolution of the human X chromosome is the analysis of synteny of X-linked genes in different species. We have focused on the spermine synthase gene in human Xp22. 1. We show that whereas the human gene spans a genomic region of 54 kb, the Fugu rubripes gene is encompassed in a 4.7-kb region. However, we could not find conserved synteny between this region of human Xp22 and the equivalent F. rubripes region. A cosmid clone containing the F. rubripes gene does not contain other X-linked genes. Instead we identified homologs of human genes that are autosomally localized: the ryanodine receptor type I (RYRI), which is implicated in malignant hyperthermia and central core disease, and the HE6 gene. Comparison of the F. rubripes, Tetraodon fluviatilis, mouse, human, and Danio rerio 5'UTRs of spermine synthase highlights conserved sequences potentially involved in regulation. Interestingly, pseudogenes of this gene that are present in the human and mouse genomes seem to be absent in the compact F. rubripes genome. Analysis of a D. rerio PAC clone containing spermine synthase shows an intermediate genomic size in this fish. Sequence analysis of this PAC clone did not reveal other known genes: neither the RYRI gene, nor the HE6 gene, nor other human Xp22 genes were identified.

zfNLRR, a novel leucine-rich repeat protein is preferentially expressed during regeneration in zebrafish

zfNLRR is a novel transmembrane protein that is most prominently expressed during regeneration of the zebrafish central nervous system. Retinal ganglion cells and descending spinal cord neurons strongly increase zfNLRR mRNA levels after axotomy in the adult. In contrast, during development expression is hardly detectable and is restricted to a few sensory systems. In the adult brain, zfNLRR mRNA is found at low levels in several motor and premotor systems. Sequence analysis reveals that zfNLRR contains in its extracellular region 12 leucine-rich repeats, 1 immunoglobulin-like domain and 1 fibronectin type III-like domain. The same protein binding motifs were identified in transmembrane proteins from frog, mouse, and human. Together, they constitute a novel family of vertebrate neuronal leucine-rich repeat proteins. Three distinct isoforms are identified so far. On the basis of its structural features and expression pattern, we propose that zfNLRR functions as a neuronal-specific adhesion molecule or soluble ligand binding receptor, primarily during restoration of the nervous system after injury.

A role for DNA methylation in gastrulation and somite patterning

DNA methylation constitutes an important epigenetic factor in the control of genetic information. In this study, we analyzed expression of the DNA methyltransferase gene and examined DNA methylation patterns during early development of the zebrafish. Maternal transcripts of the zebrafish DNA methyltransferase gene (MTase) are ubiquitously present at high levels in early embryos with overall levels decreasing after the blastula stage. At 24 h, methyltransferase mRNA is predominantly found in the brain, neural tube, eyes, and differentiating somites. Expression of MTase in the somites is highest in the anterior cells of the somites. Despite the high levels of MTase mRNA in blastula-stage embryos, we observe DNA hypomethylation at the blastula and gastrula stages compared to sperm or older embryos. Zebrafish embryos treated with 5-azacytidine (5-azaC) and 5-aza-2-deoxycytidine (5-azadC), nucleotide analogs known to induce cellular differentiation and DNA hypomethylation in mammalian cells, exhibit DNA hypomethylation and developmental perturbations. These defects are specifically observed in embryos treated at the beginning of the blastula period, just prior to midblastula transition. The most common phenotype is the loss of tail and abnormal patterning of somites. Head development is also affected in some embryos. Histological and in situ hybridization analyses reveal whole or partial loss of a differentiated notochord and midline muscle in treated embryos. When examined during gastrulation, 5-azaC-treated embryos have a shortened and thickened axial mesoderm. We propose that DNA methylation is required for normal gastrulation and subsequent patterning of the dorsal mesoderm.

Zebrafish YAC, BAC, and PAC genomic libraries

Numerous positional cloning projects directed at isolating genes responsible for the myriads of observed developmental defects in the zebrafish are anticipated in the very near future. In this chapter, we have reviewed the YAC, BAC, and PAC large-insert genomic resources available to the zebrafish community. We have discussed how these resources are screened and used in a positional cloning scheme and have pointed out frequently formidable logistical considerations in the approach. Despite being extremely tedious, positional cloning projects in the zebrafish will be comparatively easier to accomplish than in human and mouse, because of unique biological advantages of the zebrafish system. Moreover, the ease and speed at which genes are identified and cloned should rapidly increase as more mapping reagents and information become available, thereby paving the way for meaningful biological studies.

Retroviruses, ascorbate, and mutations, in the evolution of Homo sapiens

Mutations, induced by free radicals, provide a rich molecular palette that other evolutionary forces can select for or against. A recent hypothesis proposed that large numbers of free radicals were produced when, millions of years ago, Anthropoidea lost the ability to produce endogenous ascorbate, increasing the frequency of mutations and accelerating the evolution of higher primates. Recognizing that retroviruses have been active throughout the period of primate evolution, we suggest that an endogenous retrovirus or other retroviral-like element may have been involved in mutating the gene coding for gulonolactone oxidase (GLO), the terminal step in ascorbate synthesis, approximately 45 million years ago. This possibility is supported by the presence of Alu elements (a common primate retroelement) adjacent to the site of a missing segment of the nonfunctional GLO gene. Although Homo sapiens and other higher primates produce other endogenous antioxidants, including superoxide dismutase and uric acid, they do not quench the same radicals as ascorbate and cannot fully compensate for a lack of endogenous ascorbate. As a consequence, a retrovirus may have played a pivotal role in primate and H. sapiens evolution, and the absence of endogenous ascorbate may be continuing to accelerate the rate of H. sapiens and primate evolution.

Use of the zebrafish (Danio rerio) to define hematopoiesis

Hematopoiesis in the vertebrate is characterized by the induction of ventral mesoderm to form hematopoietic stem cells and the eventual differentiation of these progenitors to form the peripheral blood lineages. Several genes have been implicated in the differentiation and development of hematopoietic and vascular progenitor cells, yet our understanding of the discrete steps involved in the induction of these cells from the ventral mesoderm is still incomplete. One method of delineating these processes is based on the use of lower vertebrates. The zebrafish (Danio rerio) is an especially robust vertebrate system for both isolating and characterizing genes involved in these processes. Hematopoietic mutants have been generated with defects in many of the steps of both the primitive and definitive hematopoietic programs. Cloning of the genes that underlie these mutations should yield valuable details of hematopoiesis and may have therapeutic implications for bone marrow transplantation and stem cell gene therapy.

Rapid identification and isolation of zebrafish cDNA clones

A fast and economical approach, referred to as cDNA clone tagging, was adapted to identify and isolate zebrafish cDNA clones. The basic approach was to partially sequence the coding region of size selected cDNA clones and the partial sequences were then used as tags for identifying the clones through homology search. To benefit maximally from the tagging approach, two cDNA libraries, derived from embryonic and adult fish poly(A)+ RNAs, respectively, were constructed by unidirectional cloning; conceptually, they have the potential to represent all expressed zebrafish genes. A total of 1084 clones were sequenced from the two libraries, and 511 clones were identified, based on sequence homology. These identified clones were derived from at least 261 genes, encoding 48 translational machinery proteins, 47 cytosolic proteins, 43 cytoskeletal proteins, 41 nuclear proteins, 32 membrane proteins, 22 secreted proteins, 20 mitochondrial proteins and 8 proteins with an unknown location. Of the 261 distinct cDNA clones identified, 254 were isolated for the first time in the zebrafish. These tagged cDNA clones, identified and unidentified, provide rich resources for developmental analysis as well as mapping of zebrafish genome. The long-term objective of this study is to establish a tagged zebrafish gene library that can be accessed both by hybridization screening against the plasmid DNAs and by electronic screening using the sequence information.

Repetitive elements and their genetic applications in zebrafish

Repetitive elements provide important clues about chromosome dynamics, evolutionary forces, and mechanisms for exchange of genetic information between organisms. Repetitive sequences, especially the mobile elements, have many potential applications in genetic research. DNA transposons and retroposons are routinely used for insertional mutagenesis, gene mapping, gene tagging, and gene transfer in several model systems. Once they are developed for the zebrafish, they will greatly facilitate the identification, mapping, and isolation of genes involved in development as well as the investigation of the evolutionary processes that have been shaping eukaryotic genomes. In this review repetitive elements are characterized in terms of their lengths and other physical properties, copy numbers, modes of amplification, and mobilities within a single genome and between genomes. Examples of how they can be used to screen genomes for species and individual strain differences are presented. This review does not cover repetitive gene families that encode well-studied products such as rRNAs, tRNAs, and the like.

The zebrafish BMP4 gene: sequence analysis and expression pattern during embryonic development

We have isolated zebrafish BMP4 gene from a zebrafish genomic DNA library. The size of the isolated BMP4 gene was approximately 14.9 kb. The isolated gene contained two exons which formed the complete coding region together with part of the 3'-noncoding region. The deduced BMP4 protein sequence contained 400 amino acids. Sequence comparison showed that it shared 73% amino acid sequence identity with that of human and mouse BMP4. An intron with a size of 8,963 bp was present between two coding exons. Danio retroposon A (DANA)-like retroposon was located in the intron. It contained four conserved boxes and was flanked by a pair of direct repeats of 9 nucleotide sequence (GTTTTAATA). During embryonic development of the zebrafish, a 3.8-kb BMP4 mRNA was detected from gastrula stage up to a month-old hatching larvae via Northern blot analysis. In addition, the use of reverse transcription polymerase chain reaction further demonstrated the presence of BMP4 mRNA in both the early developmental stages (i.e., cleavage and blastula) and in adult fish. Developmental expression of BMP4 protein was also analyzed. Trace amounts of an 18-kD protein were detected at pharyngula stage, while the production increased from hatching larvae to adult fish. In adult fish, the expression of BMP4 mRNA was observed in brain, heart, digestive tracts, testes, and jaw. The results suggest that the zebrafish BMP4 gene may play important roles during zebrafish development.