A sequence-based variation map of zebrafish

Immune gene variation associated with chromosome-scale differences among individual zebrafish genomes

Immune genes have evolved to maintain exceptional diversity, offering robust defense against pathogens. We performed genomic assembly to examine immune gene variation in zebrafish. Gene pathway analysis identified immune genes as significantly enriched among genes with evidence of positive selection. A large subset of genes was absent from analysis of coding sequences due to apparent lack of reads, prompting us to examine genes overlapping zero coverage regions (ZCRs), defined as 2 kb stretches without mapped reads. Immune genes were identified as highly enriched within ZCRs, including over 60% of major histocompatibility complex (MHC) genes and NOD-like receptor (NLR) genes, mediators of direct and indirect pathogen recognition. This variation was most highly concentrated throughout one arm of chromosome 4 carrying a large cluster of NLR genes, associated with large-scale structural variation covering more than half of the chromosome. Our genomic assemblies uncovered alternative haplotypes and distinct complements of immune genes among individual zebrafish, including the MHC Class II locus on chromosome 8 and the NLR gene cluster on chromosome 4. While previous studies have shown marked variation in NLR genes between vertebrate species, our study highlights extensive variation in NLR gene regions between individuals of the same species. Taken together, these findings provide evidence of immune gene variation on a scale previously unknown in other vertebrate species and raise questions about potential impact on immune function.

The Laboratory Domestication of Zebrafish: From Diverse Populations to Inbred Substrains

We know from human genetic studies that practically all aspects of biology are strongly influenced by the genetic background, as reflected in the advent of “personalized medicine.” Yet, with few exceptions, this is not taken into account when using laboratory populations as animal model systems for research in these fields. Laboratory strains of zebrafish (Danio rerio) are widely used for research in vertebrate developmental biology, behavior, and physiology, for modeling diseases, and for testing pharmaceutic compounds in vivo. However, all of these strains are derived from artificial bottleneck events and therefore are likely to represent only a fraction of the genetic diversity present within the species. Here, we use restriction site-associated DNA sequencing to genetically characterize wild populations of zebrafish from India, Nepal, and Bangladesh, and to compare them to previously published data on four common laboratory strains. We measured nucleotide diversity, heterozygosity, and allele frequency spectra, and find that wild zebrafish are much more diverse than laboratory strains. Further, in wild zebrafish, there is a clear signal of GC-biased gene conversion that is missing in laboratory strains. We also find that zebrafish populations in Nepal and Bangladesh are most distinct from all other strains studied, making them an attractive subject for future studies of zebrafish population genetics and molecular ecology. Finally, isolates of the same strains kept in different laboratories show a pattern of ongoing differentiation into genetically distinct substrains. Together, our findings broaden the basis for future genetic, physiological, pharmaceutic, and evolutionary studies in Danio rerio.

Development and characterization of a new DRCF cell line from Indian wild strain zebrafish Danio rerio (Hamilton 1822)

Danio rerio, zebrafish, has been widely used as a non-mammalian vertebrate model organism in various studies. The present research describes to develop and characterize a new cell line from a wild strain Indian zebrafish native to Brahmaputra River, Assam, India. The new cell line designated as DRCF was developed from the caudal fin of D. rerio. The cell line was successfully subcultured up to 31 passages. Growth studies revealed that cell growth of DRCF was optimal at 28 °C in L-15 medium supplemented with 20% FBS. Molecular characterization of the DRCF cell line using mitochondrial genes namely cytochrome oxidase subunit I gene (COI) and 16S rRNA authenticated the true origin of the cell line. The chromosome analysis of the DRCF cell line expressed its 50 diploid chromosome number of D. rerio. The immunocytochemical characterization of the cell line exhibited its fibroblastic morphology. The expression of the green fluorescent protein (GFP) following transfection revealed the suitability of the cell line for transfection studies.

Linking Virus Discovery to Immune Responses Visualized during Zebrafish Infections

The discovery of new viruses currently outpaces our capacity for experimental examination of infection biology. To better couple virus discovery with immunology, we genetically modified zebrafish to visually report on virus infections. After generating a strain that expresses green fluorescent protein (GFP) under an interferon-stimulated gene promoter, we repeatedly observed transgenic larvae spontaneously expressing GFP days after hatching. RNA sequencing comparisons of co-housed GFP-positive and GFP-negative zebrafish revealed a naturally occurring picornavirus that induced a canonical interferon-mediated response and hundreds of antiviral defense genes not observed following immunostimulatory treatments or experimental infections with other viruses. Among the many genes induced by picornavirus infection was a large set encoding guanosine triphosphatase (GTPase) of immunity-associated proteins (GIMAPs). The GIMAP gene family is massively expanded in fish genomes and may also play a crucial role in antiviral responses in mammals, including humans. We subsequently detected zebrafish picornavirus in publicly available sequencing data from seemingly asymptomatic zebrafish in many research institutes and found that it altered gene expression in a previous study of zebrafish development. Experiments revealed a horizontal mode of virus transmission, highlighting a system for studying the spread of picornavirus infections within and between individuals. Our study describes a naturally occurring picornavirus that elicits strong antiviral responses in zebrafish and provides new strategies for simultaneously discovering viruses and their impact on vertebrate hosts.

Genome and Transcriptome Sequencing of casper and roy Zebrafish Mutants Provides Novel Genetic Clues for Iridophore Loss

casper has been a widely used transparent mutant of zebrafish. It possesses a combined loss of reflective iridophores and light-absorbing melanophores, which gives rise to its almost transparent trunk throughout larval and adult stages. Nevertheless, genomic causal mutations of this transparent phenotype are poorly defined. To identify the potential genetic basis of this fascinating morphological phenotype, we constructed genome maps by performing genome sequencing of 28 zebrafish individuals including wild-type AB strain, roy orbison (roy), and casper mutants. A total of 4.3 million high-quality and high-confidence homozygous single nucleotide polymorphisms (SNPs) were detected in the present study. We also identified a 6.0-Mb linkage disequilibrium block specifically in both roy and casper that was composed of 39 functional genes, of which the mpv17 gene was potentially involved in the regulation of iridophore formation and maintenance. This is the first report of high-confidence genomic mutations in the mpv17 gene of roy and casper that potentially leads to defective splicing as one major molecular clue for the iridophore loss. Additionally, comparative transcriptomic analyses of skin tissues from the AB, roy and casper groups revealed detailed transcriptional changes of several core genes that may be involved in melanophore and iridophore degeneration. In summary, our updated genome and transcriptome sequencing of the casper and roy mutants provides novel genetic clues for the iridophore loss. These new genomic variation maps will offer a solid genetic basis for expanding the zebrafish mutant database and in-depth investigation into pigmentation of animals.

Pattern of fin rays along the antero-posterior axis based on their connection to distal radials

BACKGROUND

Teleost paired fins are composed of two endoskeletal domains, proximal and distal radials, and an exoskeletal domain, the fin ray. The zebrafish pectoral fin displays elaborately patterned radials along the anteroposterior (AP) axis. Radials are considered homologous to tetrapod limb skeletons, and their patterning mechanisms in embryonic development are similar to those of limb development. Nevertheless, the pattern along the AP axis in fin rays has not been well described in the zebrafish pectoral fin, although several recent reports have revealed that fin ray development shares some cellular and genetic properties with fin/limb endoskeleton development. Thus, fin ray morphogenesis may involve developmental mechanisms for AP patterning in the fin/limb endoskeleton, and may have a specific pattern along the AP axis.

RESULTS

We conducted detailed morphological observations on fin rays and their connection to distal radials by comparing intra- and inter-strain zebrafish specimens. Although the number of fin rays varied, pectoral fin rays could be categorized into three domains along the AP axis, according to the connection between the fin rays and distal radials; additionally, the number of fin rays varied in the posterior part of the three domains. This result was confirmed by observation of the morphogenesis process of fin rays and distal radials, which showed altered localization of distal radials in the middle domain. We also evaluated the expression pattern of lhx genes, which have AP patterning activity in limb development, in fin rays and during distal radial development and found these genes to be expressed during morphogenesis in both fin rays and distal radials.

CONCLUSION

The fin ray and its connection to the endoskeleton are patterned along the AP axis, and the pattern along the AP axis in the fin ray and the radial connection is constructed by the developmental mechanism related to AP patterning in the limb/fin bud. Our results indicate the possibility that the developmental mechanisms of fin rays and their connection are comparable to those of the distal element of the limb skeleton.

A genome-wide map of circular RNAs in adult zebrafish

Circular RNAs (circRNAs) are transcript isoforms generated by back-splicing of exons and circularisation of the transcript. Recent genome-wide maps created for circular RNAs in humans and other model organisms have motivated us to explore the repertoire of circular RNAs in zebrafish, a popular model organism. We generated RNA-seq data for five major zebrafish tissues - Blood, Brain, Heart, Gills and Muscle. The repertoire RNA sequence reads left over after reference mapping to linear transcripts were used to identify unique back-spliced exons utilizing a split-mapping algorithm. Our analysis revealed 3,428 novel circRNAs in zebrafish. Further in-depth analysis suggested that majority of the circRNAs were derived from previously well-annotated protein-coding and long noncoding RNA gene loci. In addition, many of the circular RNAs showed extensive tissue specificity. We independently validated a subset of circRNAs using polymerase chain reaction (PCR) and divergent set of primers. Expression analysis using quantitative real time PCR recapitulate selected tissue specificity in the candidates studied. This study provides a comprehensive genome-wide map of circular RNAs in zebrafish tissues.

An Interrogation of Shared and Unique Copy Number Variants Across Genetically Distinct Zebrafish Strains

Zebrafish (Danio rerio) are a widely utilized model system for human disorders, but common laboratory strains have distinct behavioral and physiological differences. Accompanying these known strain differences, commonly used "wildtype" zebrafish strains have both shared and unique suites of single nucleotide polymorphisms and copy number variants (CNVs). Despite this, genomic variation is often ignored in study design, and the actual strain used is often not adequately reported. The goal of this study was to assess CNVs across three common laboratory strains of zebrafish-AB, Tubingen (TU), and WIK-and provide these data as a tool for the zebrafish community. Herein we identified 1351 CNV regions within the most recent genome assembly (GRCz11) covering 1.9% of the zebrafish genome (31.7 Mb). CNVs were found across all chromosomes, and 2200 genes (5121 transcripts) lie within ±5 kb of identified CNVs, pointing to likely cis regulatory actions of CNVs on nearby gene neighbors. We have created a Public Session accessible on the UCSC Genome Browser to view CNVs from this study titled "danRer11 zebrafish CNV across strains" as a tool for the zebrafish community.

Elucidating Gene-by-Environment Interactions Associated with Differential Susceptibility to Chemical Exposure

BACKGROUND

Modern societies are exposed to vast numbers of potentially hazardous chemicals. Despite demonstrated linkages between chemical exposure and severe health effects, there are limited, often conflicting, data on how adverse health effects of exposure differ across individuals.

OBJECTIVES

We tested the hypothesis that population variability in response to certain chemicals could elucidate a role for gene-environment interactions (GxE) in differential susceptibility.

METHODS

High-throughput screening (HTS) data on thousands of chemicals in genetically heterogeneous zebrafish were leveraged to identify a candidate chemical (Abamectin) with response patterns indicative of population susceptibility differences. We tested the prediction by generating genome-wide sequence data for 276 individual zebrafish displaying susceptible (Affected) vs. resistant (Unaffected) phenotypes following identical chemical exposure.

RESULTS

We found GxE associated with differential susceptibility in the sox7 promoter region and then confirmed gene expression differences between phenotypic response classes.

CONCLUSIONS

The results for Abamectin in zebrafish demonstrate that GxE associated with naturally occurring, population genetic variation play a significant role in mediating individual response to chemical exposure. https://doi.org/10.1289/EHP2662.

RNA secondary structure profiling in zebrafish reveals unique regulatory features

Background

RNA is known to play diverse roles in gene regulation. The clues for this regulatory function of RNA are embedded in its ability to fold into intricate secondary and tertiary structure.

Results

We report the transcriptome-wide RNA secondary structure in zebrafish at single nucleotide resolution using Parallel Analysis of RNA Structure (PARS). This study provides the secondary structure map of zebrafish coding and non-coding RNAs. The single nucleotide pairing probabilities of 54,083 distinct transcripts in the zebrafish genome were documented. We identified RNA secondary structural features embedded in functional units of zebrafish mRNAs. Translation start and stop sites were demarcated by weak structural signals. The coding regions were characterized by the three-nucleotide periodicity of secondary structure and display a codon base specific structural constrain. The splice sites of transcripts were also delineated by distinct signature signals. Relatively higher structural signals were observed at 3’ Untranslated Regions (UTRs) compared to Coding DNA Sequence (CDS) and 5’ UTRs. The 3′ ends of transcripts were also marked by unique structure signals. Secondary structural signals in long non-coding RNAs were also explored to better understand their molecular function.

Conclusions

Our study presents the first PARS-enabled transcriptome-wide secondary structure map of zebrafish, which documents pairing probability of RNA at single nucleotide precision. Our findings open avenues for exploring structural features in zebrafish RNAs and their influence on gene expression.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4497-0) contains supplementary material, which is available to authorized users.

Population genetic diversity in zebrafish lines

Toxicological and pharmacological researchers have seized upon the many benefits of zebrafish, including the short generation time, well-characterized development, and early maturation as clear embryos. A major difference from many model organisms is that standard husbandry practices in zebrafish are designed to maintain population diversity. While this diversity is attractive for translational applications in human and ecological health, it raises critical questions on how interindividual genetic variation might contribute to chemical exposure or disease susceptibility differences. Findings from pooled samples of zebrafish support this supposition of diversity yet cannot directly measure allele frequencies for reference versus alternate alleles. Using the Tanguay lab Tropical 5D zebrafish line (T5D), we performed whole genome sequencing on a large group (n = 276) of individual zebrafish embryos. Paired-end reads were collected on an Illumina 3000HT, then aligned to the most recent zebrafish reference genome (GRCz10). These data were used to compare observed population genetic variation across species (humans, mice, zebrafish), then across lines within zebrafish. We found more single nucleotide polymorphisms (SNPs) in T5D than have been reported in SNP databases for any of the WIK, TU, TL, or AB lines. We theorize that some subset of the novel SNPs may be shared with other zebrafish lines but have not been identified in other studies due to the limitations of capturing population diversity in pooled sequencing strategies. We establish T5D as a model that is representative of diversity levels within laboratory zebrafish lines and demonstrate that experimental design and analysis can exert major effects when characterizing genetic diversity in heterogeneous populations.

Zebrafish: A Premier Vertebrate Model for Biomedical Research in Indian Scenario

The zebrafish (Danio rerio) is a versatile model organism that has been used in biomedical research for several decades to study a wide range of biological phenomena. There are many technical advantages of using zebrafish over other vertebrate models. They are readily available, hardy, easy, and inexpensive to maintain in the laboratory, have a short life cycle, and have excellent fecundity. Due to its optical clarity and reproducible capabilities, it has become one of the predominant models of human genetic diseases. Zebrafish research has made rapid strides in the United States and Europe, but in India the field is at an early stage and many researchers still remain unaware of the full research potential of this tiny fish. The zebrafish model system was introduced into India in the early 2000s. Up to now, more than 200 scientific referred articles have been published by Indian researchers. This review gives an overview of the current state of knowledge for zebrafish research in India, with the aim of promoting wider utilization of zebrafish for high level biological studies.

Loss of the chromatin modifier Kdm2aa causes BrafV600E-independent spontaneous melanoma in zebrafish

KDM2A is a histone demethylase associated with transcriptional silencing, however very little is known about its in vivo role in development and disease. Here we demonstrate that loss of the orthologue kdm2aa in zebrafish causes widespread transcriptional disruption and leads to spontaneous melanomas at a high frequency. Fish homozygous for two independent premature stop codon alleles show reduced growth and survival, a strong male sex bias, and homozygous females exhibit a progressive oogenesis defect. kdm2aa mutant fish also develop melanomas from early adulthood onwards which are independent from mutations in braf and other common oncogenes and tumour suppressors as revealed by deep whole exome sequencing. In addition to effects on translation and DNA replication gene expression, high-replicate RNA-seq in morphologically normal individuals demonstrates a stable regulatory response of epigenetic modifiers and the specific de-repression of a group of zinc finger genes residing in constitutive heterochromatin. Together our data reveal a complex role for Kdm2aa in regulating normal mRNA levels and carcinogenesis. These findings establish kdm2aa mutants as the first single gene knockout model of melanoma biology.

India's Computational Biology Growth and Challenges

India's computational science is growing swiftly due to the outburst of internet and information technology services. The bioinformatics sector of India has been transforming rapidly by creating a competitive position in global bioinformatics market. Bioinformatics is widely used across India to address a wide range of biological issues. Recently, computational researchers and biologists are collaborating in projects such as database development, sequence analysis, genomic prospects and algorithm generations. In this paper, we have presented the Indian computational biology scenario highlighting bioinformatics-related educational activities, manpower development, internet boom, service industry, research activities, conferences and trainings undertaken by the corporate and government sectors. Nonetheless, this new field of science faces lots of challenges.

Working with zebrafish at postembryonic stages

As the processes of embryogenesis become increasingly well understood, there is growing interest in the development that occurs at later, postembryonic stages. Postembryonic development holds tremendous potential for discoveries of both fundamental and translational importance. Zebrafish, which are small, rapidly and externally developing, and which boast a wealth of genetic resources, are an outstanding model of vertebrate postembryonic development. Nonetheless, there are specific challenges posed by working with zebrafish at these stages, and this chapter is meant to serve as a primer for those working with larval and juvenile zebrafish. Since accurate staging is critical for high-quality results and experimental reproducibility, we outline best practices for reporting postembryonic developmental progress. Emphasizing the importance of accurate staging, we present new data showing that rates of growth and size-stage relationships can differ even between wild-type strains. Finally, since rapid and uniform development is particularly critical when working at postembryonic stages, we briefly describe methods that we use to achieve high rates of growth and developmental uniformity through postembryonic stages in both wild-type and growth-compromised zebrafish.

Pseudo-De Novo Assembly and Analysis of Unmapped Genome Sequence Reads in Wild Zebrafish Reveal Novel Gene Content

Zebrafish represents the third vertebrate with an officially completed genome, yet it remains incomplete with additions and corrections continuing with the current release, GRCz10, having 13% of zebrafish cDNA sequences unmapped. This disparity may result from population differences, given that the genome reference was generated from clonal individuals with limited genetic diversity. This is supported by the recent analysis of a single wild zebrafish, which identified over 5.2 million SNPs and 1.6 million in/dels in the previous genome build, zv9. Re-examination of this sequence data set indicated that 13.8% of quality sequence reads failed to align to GRCz10. Using a novel bioinformatics de novo assembly pipeline on these unmappable reads, we identified 1,514,491 novel contigs covering ∼224 Mb of genomic sequence. Among these, 1083 contigs were found to contain a potential gene coding sequence. RNA-seq data comparison confirmed that 362 contigs contained a transcribed DNA sequence, suggesting that a large amount of functional genomic sequence remains unannotated in the zebrafish reference genome. By utilizing the bioinformatics pipeline developed in this study, the zebrafish genome will be bolstered as a model for human disease research. Adaptation of the pipeline described here also offers a cost-efficient and effective method to identify and map novel genetic content across any genome and will ultimately aid in the completion of additional genomes for a broad range of species.

Chamber Specific Gene Expression Landscape of the Zebrafish Heart

The organization of structure and function of cardiac chambers in vertebrates is defined by chamber-specific distinct gene expression. This peculiarity and uniqueness of the genetic signatures demonstrates functional resolution attributed to the different chambers of the heart. Altered expression of the cardiac chamber genes can lead to individual chamber related dysfunctions and disease patho-physiologies. Information on transcriptional repertoire of cardiac compartments is important to understand the spectrum of chamber specific anomalies. We have carried out a genome wide transcriptome profiling study of the three cardiac chambers in the zebrafish heart using RNA sequencing. We have captured the gene expression patterns of 13,396 protein coding genes in the three cardiac chambers—atrium, ventricle and bulbus arteriosus. Of these, 7,260 known protein coding genes are highly expressed (≥10 FPKM) in the zebrafish heart. Thus, this study represents nearly an all-inclusive information on the zebrafish cardiac transcriptome. In this study, a total of 96 differentially expressed genes across the three cardiac chambers in zebrafish were identified. The atrium, ventricle and bulbus arteriosus displayed 20, 32 and 44 uniquely expressing genes respectively. We validated the expression of predicted chamber-restricted genes using independent semi-quantitative and qualitative experimental techniques. In addition, we identified 23 putative novel protein coding genes that are specifically restricted to the ventricle and not in the atrium or bulbus arteriosus. In our knowledge, these 23 novel genes have either not been investigated in detail or are sparsely studied. The transcriptome identified in this study includes 68 differentially expressing zebrafish cardiac chamber genes that have a human ortholog. We also carried out spatiotemporal gene expression profiling of the 96 differentially expressed genes throughout the three cardiac chambers in 11 developmental stages and 6 tissue types of zebrafish. We hypothesize that clustering the differentially expressed genes with both known and unknown functions will deliver detailed insights on fundamental gene networks that are important for the development and specification of the cardiac chambers. It is also postulated that this transcriptome atlas will help utilize zebrafish in a better way as a model for studying cardiac development and to explore functional role of gene networks in cardiac disease pathogenesis.

The identification of additional zebrafish DICP genes reveals haplotype variation and linkage to MHC class I genes

Bony fish encode multiple multi-gene families of membrane receptors that are comprised of immunoglobulin (Ig) domains and are predicted to function in innate immunity. One of these families, the diverse immunoglobulin (Ig) domain-containing protein (DICP) genes, maps to three chromosomal loci in zebrafish. Most DICPs possess one or two Ig ectodomains and include membrane-bound and secreted forms. Membrane-bound DICPs include putative inhibitory and activating receptors. Recombinant DICP Ig domains bind lipids with varying specificity, a characteristic shared with mammalian CD300 and TREM family members. Numerous DICP transcripts amplified from different lines of zebrafish did not match the zebrafish reference genome sequence suggesting polymorphic and haplotypic variation. The expression of DICPs in three different lines of zebrafish has been characterized employing PCR-based strategies. Certain DICPs exhibit restricted expression in adult tissues whereas others are expressed ubiquitously. Transcripts of a subset of DICPs can be detected during embryonic development suggesting roles in embryonic immunity or other developmental processes. Transcripts representing 11 previously uncharacterized DICP sequences were identified. The assignment of two of these sequences to an unplaced genomic scaffold resulted in the identification of an alternative DICP haplotype that is linked to a MHC class I Z lineage haplotype on zebrafish chromosome 3. The linkage of DICP and MHC class I genes also is observable in the genomes of the related grass carp (Ctenopharyngodon idellus) and common carp (Cyprinus carpio) suggesting that this is a shared character with the last common Cyprinidae ancestor.

Short stories on zebrafish long noncoding RNAs

The recent re-annotation of the transcriptome of human and other model organisms, using next-generation sequencing approaches, has unravelled a hitherto unknown repertoire of transcripts that do not have a potential to code for proteins. These transcripts have been largely classified into an amorphous class popularly known as long noncoding RNAs (lncRNA). This discovery of lncRNAs in human and other model systems have added a new layer to the understanding of gene regulation at the transcriptional and post-transcriptional levels. In recent years, three independent studies have discovered a number of lncRNAs expressed in different stages of zebrafish development and adult tissues using a high-throughput RNA sequencing approach, significantly adding to the repertoire of genes known in zebrafish. A subset of these transcripts also shows distinct and specific spatiotemporal patterns of gene expression, pointing to a tight regulatory control and potential functional roles in development, organogenesis, and/ or homeostasis. This review provides an overview of the lncRNAs in zebrafish and discusses how their discovery could provide new insights into understanding biology, explaining mutant phenotypes, and helping in potentially modeling disease processes.

SNPfisher: tools for probing genetic variation in laboratory-reared zebrafish

Single nucleotide polymorphisms (SNPs) are the benchmark molecular markers for modern genomics. Until recently, relatively few SNPs were known in the zebrafish genome. The use of next-generation sequencing for the positional cloning of zebrafish mutations has increased the number of known SNP positions dramatically. Still, the identified SNP variants remain under-utilized, owing to scant annotation of strain specificity and allele frequency. To address these limitations, we surveyed SNP variation in three common laboratory zebrafish strains using whole-genome sequencing. This survey identified an average of 5.04 million SNPs per strain compared with the Zv9 reference genome sequence. By comparing the three strains, 2.7 million variants were found to be strain specific, whereas the remaining variants were shared among all (2.3 million) or some of the strains. We also demonstrate the broad usefulness of our identified variants by validating most in independent populations of the same laboratory strains. We have made all of the identified SNPs accessible through 'SNPfisher', a searchable online database (snpfisher.nichd.nih.gov). The SNPfisher website includes the SNPfisher Variant Reporter tool, which provides the genomic position, alternate allele read frequency, strain specificity, restriction enzyme recognition site changes and flanking primers for all SNPs and Indels in a user-defined gene or region of the zebrafish genome. The SNPfisher site also contains links to display our SNP data in the UCSC genome browser. The SNPfisher tools will facilitate the use of SNP variation in zebrafish research as well as vertebrate genome evolution.

Advancing biology through a deeper understanding of zebrafish ecology and evolution

Over the last two decades, the zebrafish has joined the ranks of premier model organisms for biomedical research, with a full suite of tools and genomic resources. Yet we still know comparatively little about its natural history. Here I review what is known about the natural history of the zebrafish, where significant gaps in our knowledge remain, and how a fuller appreciation of this organism's ecology and behavior, population genetics, and phylogeny can inform a variety of research endeavors.

Multigene families of immunoglobulin domain-containing innate immune receptors in zebrafish: deciphering the differences

Five large multigene families encoding innate-type immune receptors that are comprised of immunoglobulin domains have been identified in bony fish, of which four do not possess definable mammalian orthologs. The members of some of the multigene families exhibit unusually extensive patterns of divergence and the individual family members demonstrate marked variation in interspecific comparisons. As a group, the gene families reveal striking differences in domain type and content, mechanisms of intracellular signaling, basic structural features, haplotype and allelic variation and ligand binding. The potential functional roles of these innate immune receptors, their relationships to immune genes in higher vertebrate species and the basis for their adaptive evolution are of broad interest. Ongoing investigations are expected to provide new insight into alternative mechanisms of immunity.

The MHC class I genes of zebrafish

Major histocompatibility complex (MHC) molecules play a central role in the immune response and in the recognition of non-self. Found in all jawed vertebrate species, including zebrafish and other teleosts, MHC genes are considered the most polymorphic of all genes. In this review we focus on the multi-faceted diversity of zebrafish MHC class I genes, which are classified into three sequence lineages: U, Z, and L. We examine the polygenic, polymorphic, and haplotypic diversity of the zebrafish MHC class I genes, discussing known and postulated functional differences between the different class I lineages. In addition, we provide the first comprehensive nomenclature for the L lineage genes in zebrafish, encompassing at least 15 genes, and characterize their sequence properties. Finally, we discuss how recent findings have shed new light on the remarkably diverse MHC loci of this species.

A defined zebrafish line for high-throughput genetics and genomics: NHGRI-1

Substantial intrastrain variation at the nucleotide level complicates molecular and genetic studies in zebrafish, such as the use of CRISPRs or morpholinos to inactivate genes. In the absence of robust inbred zebrafish lines, we generated NHGRI-1, a healthy and fecund strain derived from founder parents we sequenced to a depth of ∼50×. Within this strain, we have identified the majority of the genome that matches the reference sequence and documented most of the variants. This strain has utility for many reasons, but in particular it will be useful for any researcher who needs to know the exact sequence (with all variants) of a particular genomic region or who wants to be able to robustly map sequences back to a genome with all possible variants defined.

The Zebrafish GenomeWiki: a crowdsourcing approach to connect the long tail for zebrafish gene annotation

A large repertoire of gene-centric data has been generated in the field of zebrafish biology. Although the bulk of these data are available in the public domain, most of them are not readily accessible or available in nonstandard formats. One major challenge is to unify and integrate these widely scattered data sources. We tested the hypothesis that active community participation could be a viable option to address this challenge. We present here our approach to create standards for assimilation and sharing of information and a system of open standards for database intercommunication. We have attempted to address this challenge by creating a community-centric solution for zebrafish gene annotation. The Zebrafish GenomeWiki is a 'wiki'-based resource, which aims to provide an altruistic shared environment for collective annotation of the zebrafish genes. The Zebrafish GenomeWiki has features that enable users to comment, annotate, edit and rate this gene-centric information. The credits for contributions can be tracked through a transparent microattribution system. In contrast to other wikis, the Zebrafish GenomeWiki is a 'structured wiki' or rather a 'semantic wiki'. The Zebrafish GenomeWiki implements a semantically linked data structure, which in the future would be amenable to semantic search. Database URL: http://genome.igib.res.in/twiki.

Zebrafish sex: a complicated affair

In this review, we provide a detailed overview of studies on the elusive sex determination (SD) and gonad differentiation mechanisms of zebrafish (Danio rerio). We show that the data obtained from most studies are compatible with polygenic sex determination (PSD), where the decision is made by the allelic combinations of several loci. These loci are typically dispersed throughout the genome, but in some teleost species a few of them might be located on a preferential pair of (sex) chromosomes. The PSD system has a much higher level of variation of SD genotypes both at the level of gametes and the sexual genotype of individuals, than that of the chromosomal sex determination systems. The early sexual development of zebrafish males is a complicated process, as they first develop a ‘juvenile ovary’, that later undergoes a transformation to give way to a testis. To date, three major developmental pathways were shown to be involved with gonad differentiation through the modulation of programmed cell death. In our opinion, there are more pathways participating in the regulation of zebrafish gonad differentiation/transformation. Introduction of additional powerful large-scale genomic approaches into the analysis of zebrafish reproduction will result in further deepening of our knowledge as well as identification of additional pathways and genes associated with these processes in the near future.

The CRISPR system--keeping zebrafish gene targeting fresh

We are entering a new era in our ability to modify and edit the genomes of model organisms. Zinc finger nucleases (ZFNs) opened the door to the first custom nuclease-targeted genome engineering in the late 1990s. However, ZFNs remained out of reach for most research labs because of the difficulty of production, high costs, and modest efficacy in many applications. Transcription activator-like effector nucleases (TALENs) were built upon a DNA binding system discovered in a group of plant bacterial pathogens and broadened custom nuclease technology, showing significant improvements in both targeting flexibility and efficiency. Perhaps most importantly, TALENs are open source and easy to produce, providing zebrafish laboratories around the world with affordable tools that can be made in-house rapidly, at low cost, and with reliably high activity. Now a new system for targeted genome engineering derived from the CRISPR/Cas system in eubacteria and archaea promises to simplify this process further. Together, these tools will help overcome many of the bottlenecks that have constrained gene targeting in zebrafish, paving the way for advanced genome engineering applications in this model teleost.