Zebrafish Reference Genome Research Articles

Analysis of Differential Gene Expression under Acute Lead or Mercury Exposure in Larval Zebrafish Using RNA-Seq.

This study was first conducted to investigate the effects of acute lead exposure on developing zebrafish embryos or larvae from 24 to 120 h post-fertilization (hpf). Our data showed that treatment with 50-200 μM lead significantly affected larval survivability and morphology compared to the respective control. Second, we chose 120 hpf larvae treated with 12.5 μM lead for RNA sequencing due to its exposure level being sufficient to produce toxic effects with minimum death and lead bioaccumulation in developing zebrafish. A total of 137.45 million raw reads were obtained, and more than 86% of clean data were mapped to the zebrafish reference genome. Differential expression profiles generated 116 up- and 34 down-regulated genes upon lead exposure. The most enriched GO terms for representative DEGs were ion transport and lipid metabolism. Third, a comparison with the dataset of mercury-regulated gene expression identified 94 genes (64 up-regulated and 30 down-regulated) for exposure specific to lead, as well as 422 genes (338 up-regulated and 84 down-regulated) for exposure specific to mercury. In addition, 56 genes were co-regulated by micromolar mercury and lead treatment, and the expression of thirteen genes, including mt2, ctssb.1, prdx1, txn, sqrdl, tmprss13a, socs3a, trpv6, abcb6a, gsr, hbz, fads2, and zgc:92590 were validated by qRT-PCR. These genes were mainly associated with metal ion binding, proteolysis, antioxidant activity, signal transduction, calcium ion or oxygen transport, the fatty acid biosynthetic process, and protein metabolism. Taken together, these findings help better understand the genome-wide responses of developing zebrafish to lead or mercury and provide potential biomarkers for acute exposure to toxic metals.

Extensive sequence divergence between the reference genomes of two zebrafish strains, Tuebingen and AB.

The zebrafish (Danio rerio) is one of the most widely used model organisms for studying vertebrate gene function and human disease, given the 70% conserved protein-coding genes between zebrafish and human. Two of the most common laboratory zebrafish strains are Tuebingen and AB. Despite the fact that the zebrafish reference genome is derived from the Tuebingen strain, the AB strain is still widely used although a high-quality genome comparable to Tuebingen is lacking. Here, we report a 1.40-Gb representative de novo genome assembly of the AB strain (DrAB1), with contig N50length of 21Mb, by integrating Illumina short-read sequencing, Nanopore long-read sequencing and HiC-based chromatin mapping. Compared with the published zebrafish Zv11 reference genome (GRCz11), this genome assembly shows considerable improvements in both contiguity and completeness. In addition, substantial structural differences and extensive sequence divergence of unprecedented resolution have been uncovered, especially with respect to 9,029,929single nucleotide polymorphisms, 2,376,812 InDels, 32,623 insertions, 22,089 deletions and 220 inversions, which constitute ~2.6% of the DrAB1genome. Many of these variants may have potential functional effects on phenotype, which should be considered in further experimental designs. Consequently, our study provides additional genomic resources and a high-resolution structural variation map based on whole-genome alignment for the zebrafish community, which could also be an indispensable reference genome from a model species in future research on fish phylogenetic genomics, comparative genomics and adaptive evolution.

RNA-Seq reveals differential expression profiles and functional annotation of genes involved in retinal degeneration in Pde6c mutant Danio rerio

BackgroundRetinal degenerative diseases affect millions of people and represent the leading cause of vision loss around the world. Retinal degeneration has been attributed to a wide variety of causes, such as disruption of genes involved in phototransduction, biosynthesis, folding of the rhodopsin molecule, and the structural support of the retina. The molecular pathogenesis of the biological events in retinal degeneration is unclear; however, the molecular basis of the retinal pathological defect can be potentially determined by gene-expression profiling of the whole retina. In the present study, we analyzed the differential gene expression profile of the retina from a wild-type zebrafish and phosphodiesterase 6c (pde6c) mutant.ResultsThe datasets were downloaded from the Sequence Read Archive (SRA), and adaptors and unbiased bases were removed, and sequences were checked to ensure the quality. The reads were further aligned to the reference genome of zebrafish, and the gene expression was calculated. The differentially expressed genes (DEGs) were filtered based on the log fold change (logFC) (±4) and p-values (p < 0.001). We performed gene annotation (molecular function [MF], biological process [BP], cellular component [CC]), and determined the functional pathways Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway for the DEGs. Our result showed 216 upregulated and 3527 downregulated genes between normal and pde6c mutant zebrafish. These DEGs are involved in various KEGG pathways, such as the phototransduction (12 genes), mRNA surveillance (17 genes), phagosome (25 genes), glycolysis/gluconeogenesis (15 genes), adrenergic signaling in cardiomyocytes (29 genes), ribosome (20 genes), the citrate cycle (TCA cycle; 8 genes), insulin signaling (24 genes), oxidative phosphorylation (20 genes), and RNA transport (22 genes) pathways. Many more of all the pathway genes were down-regulated, while fewer were up-regulated in the retina of pde6c mutant zebrafish.ConclusionsOur data strongly indicate that, among these genes, the above-mentioned pathways’ genes as well as calcium-binding, neural damage, peptidase, immunological, and apoptosis proteins are mostly involved in the retinal and neural degeneration that cause abnormalities in photoreceptors or retinal pigment epithelium (RPE) cells.

Tissue-Specific eQTL in Zebrafish.

Copy number variants (CNVs) refer to the loss or gain of copies of a genomic DNA region. While some CNVs may play a role in species evolution by enriching the diversity of an organism, CNVs may also be linked to certain diseases such as neurological disorders, early onset obesity, and cancer. CNVs may affect gene expression by direct overlap of the genic region or by an indirect effect where the CNV is located outside the gene location. These indirect CNV regions may contain regulatory elements such as transcription enhancers or repressors as well as regulators such as miRNAs which may work at the level of transcription or translation. Danio rerio (zebrafish) is an excellent model organism for CNV studies. Zebrafish genomes contain a large amount of variation with 14.6% of the zebrafish reference genome found to be copy number variable. This level of variation is more than four times the percentage of reference genome sequence covered by similarly common CNVs in humans. It is this high level of variation that makes zebrafish interesting to investigate the effects of CNV on gene expression. Additionally, zebrafish share 70% of genetic similarities with humans, and 84% of genes associated with human disease are also found in zebrafish. Expressive quantitative trait loci (eQTL) analysis may be used in zebrafish to explore how CNVs may affect gene expression in both a direct and indirect manner. eQTL analysis may be performed for cis associations with a 1-Mb (megabase) window upstream and downstream from the transcription probe midpoint to CGH midpoint. Trans associations (variants that are located beyond the 1-Mb window of the gene either on the same chromosome as the gene or on a different chromosome) may be investigated as well through eQTL analysis; however, trans associations require more tests to be performed than cis associations, which limits power to detect associations. Pairwise associations between each pair of copy number variant and gene will be investigated separately from the same individual using Spearman rank correlations with significant associations found being followed with a multi-test correction technique to assess significance of those CNV gene expression associations. An association between a CNV to a gene expression phenotype should be considered significant only if the p value from the analysis of the observed data is lower than the 0.001 tail threshold from a distribution of the minimal p values (which are found from all comparisons for a given gene from 10,000 permutations of the expression phenotypes). Associations between CNVs and genes may be found to be direct or indirect as well as positive (increased copy number-increased expression) or negative (increased copy number-decreased expression, decreased copy number-increased expression). Ongoing analyses with these associations will investigate the impact of CNVs on gene functionality including immune function and potential disease susceptibility.

Protein Network Analysis of the Fifth Chromosome of Zebrafish.

In the past 20 years, zebrafish has gradually become an important animal model for studying the function of human genes. At the same time, identification of zebrafish reference genome sequences and >10,000 protein-coding genes indicates that they are at least 75% homologous to human genes, further validating its utility as a research model for growth and developments (GDs). However, the molecular mechanism of zebrafish GDs has almost no molecular interactions, so a protein-protein interaction (PPI) network is highly desirable. The project extracted 58 genes on the fifth chromosome of the model organism, and studied and modeled its GDs mechanisms, which encompassed 296 interactions between the 58 proteins involved in GDs. This includes not only accurate predictions of PPI but also those molecular interactions collected from the literature and experimentally derived. These molecules then interact, modularize, look for central genes, analyze, and predict their GDs processes, and hope to help scholars study the process of GDs, providing hypotheses and help.

The Zebrafish Information Network: new support for non-coding genes, richer Gene Ontology annotations and the Alliance of Genome Resources.

The Zebrafish Information Network (ZFIN) (https://zfin.org/) is the database for the model organism, zebrafish (Danio rerio). ZFIN expertly curates, organizes and provides a wide array of zebrafish genetic and genomic data, including genes, alleles, transgenic lines, gene expression, gene function, mutant phenotypes, orthology, human disease models, nomenclature and reagents. New features at ZFIN include increased support for genomic regions and for non-coding genes, and support for more expressive Gene Ontology annotations. ZFIN has recently taken over maintenance of the zebrafish reference genome sequence as part of the Genome Reference Consortium. ZFIN is also a founding member of the Alliance of Genome Resources, a collaboration of six model organism databases (MODs) and the Gene Ontology Consortium (GO). The recently launched Alliance portal (https://alliancegenome.org) provides a unified, comparative view of MOD, GO, and human data, and facilitates foundational and translational biomedical research.

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.

Transcriptomic characterization of zebrafish larvae in response to mercury exposure

Mercury is a widespread toxicant in aquatic environment that can cause deleterious effects on fish. Although a number of mercury-regulated genes have been investigated in adult fish, the transcriptional responses of fish larvae to acute mercury exposure are not well understood. In this study, RNA sequencing was used to examine the transcriptional changes in developing zebrafish larvae under a low concentration of mercuric chloride exposure from 24 to 120hpf. Our initial results showed that a total of 142.59 million raw reads were obtained from sequencing libraries and about 86% of the processed reads were mapped to the reference genome of zebrafish. Differential expression analysis identified 391 up- and 87 down-regulated genes. Gene ontology enrichment analysis revealed that most of the differential expressed genes are closely related to the regulation of cellular process, metabolic process, multicellular organismal process, biological regulation, pigmentation, and response to stimulus. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis demonstrated that antigen processing and presentation was the most significantly enriched pathway. Moreover, we characterized a novel and sensitive mercury-induced ABCB (ATP- binding cassette B subfamily) transporter gene - abcb5. This gene is localized on zebrafish chromosome 16 and contains a 4014bp open-reading frame. The deduced polypeptide is composed of 1337 amino acids and possesses most of functional domains and critical residues defined in human and mouse ABCB5/Abcb5. Functional analysis in vitro demonstrated that overexpression of zebrafish abcb5 gene can significantly decrease the cytotoxicity of mercury in LLC-PK1 cells, implying it is a potential efflux transporter of mercury. Thus, these findings provide useful insights to help further understand the transcriptional response and detoxification ability of zebrafish larvae following acute exposure to mercury.

Alternative haplotypes of antigen processing genes in zebrafish diverged early in vertebrate evolution

Antigen processing and presentation genes found within the MHC are among the most highly polymorphic genes of vertebrate genomes, providing populations with diverse immune responses to a wide array of pathogens. Here, we describe transcriptome, exome, and whole-genome sequencing of clonal zebrafish, uncovering the most extensive diversity within the antigen processing and presentation genes of any species yet examined. Our CG2 clonal zebrafish assembly provides genomic context within a remarkably divergent haplotype of the core MHC region on chromosome 19 for six expressed genes not found in the zebrafish reference genome: mhc1uga, proteasome-β 9b (psmb9b), psmb8f, and previously unknown genes psmb13b, tap2d, and tap2e We identify ancient lineages for Psmb13 within a proteasome branch previously thought to be monomorphic and provide evidence of substantial lineage diversity within each of three major trifurcations of catalytic-type proteasome subunits in vertebrates: Psmb5/Psmb8/Psmb11, Psmb6/Psmb9/Psmb12, and Psmb7/Psmb10/Psmb13. Strikingly, nearby tap2 and MHC class I genes also retain ancient sequence lineages, indicating that alternative lineages may have been preserved throughout the entire MHC pathway since early diversification of the adaptive immune system ∼500 Mya. Furthermore, polymorphisms within the three MHC pathway steps (antigen cleavage, transport, and presentation) are each predicted to alter peptide specificity. Lastly, comparative analysis shows that antigen processing gene diversity is far more extensive than previously realized (with ancient coelacanth psmb8 lineages, shark psmb13, and tap2t and psmb10 outside the teleost MHC), implying distinct immune functions and conserved roles in shaping MHC pathway evolution throughout vertebrates.

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.

Approaches to Inactivate Genes in Zebrafish.

Animal models of tumor initiation and tumor progression are essential components toward understanding cancer and designing/validating future therapies. Zebrafish is a powerful model for studying tumorigenesis and has been successfully exploited in drug discovery. According to the zebrafish reference genome, 82 % of disease-associated genes in the Online Mendelian Inheritance in Man (OMIM) database have clear zebrafish orthologues. Using a variety of large-scale random mutagenesis methods developed to date, zebrafish can provide a unique opportunity to identify gene mutations that may be associated with cancer predisposition. On the other hand, newer technologies enabling targeted mutagenesis can facilitate reverse cancer genetic studies and open the door for complex genetic analysis of tumorigenesis. In this chapter, we will describe the various technologies for conducting genome editing in zebrafish with special emphasis on the approaches to inactivate genes.

Transcriptomic characterization of the dorsal lobes after hepatectomy of the ventral lobe in zebrafish.

BackgroundThe liver possesses an ability of compensatory growth after removing three of five lobes in mammals or one of three lobes in zebrafish. The reenter of hepatocytes into the cell cycle is one of the hallmarks for the initiation of liver compensatory growth, but cellular and molecular mechanisms underlying the activation of hepatocytes remain largely unknown.ResultsTo better understand the process, transcriptional profiles of the remaining liver dorsal lobes in female zebrafish were generated with RNA-seq. About 44 million raw reads were obtained from three sequencing libraries and 71 % of raw reads were mapped to the reference genome of zebrafish. A total number of 5652 genes were differentially expressed in at least one of two time points during the compensatory growth of liver dorsal lobes and classified into different functional categories. A number of genes encoding angiogenesis-related growth factors/ligands and apoptosis-associated cytokines were strongly expressed at 6-h time point after the removal of the ventral lobe. Gene ontology enrichment analysis of genes up-regulated during early stages of liver compensatory growth revealed that small GTPase-mediated signal transduction, RNA processing and intracellular protein transport were the most highly overrepresented biological processes and SNARE interactions in vesicular transport, proteasome and basal transcription factors were the most highly enriched pathways. Moreover, 477 genes differently expressed during liver compensatory growth of both female zebrafish and mice were involved in the response to stimulus, DNA replication, metabolic processes of fatty acid, lipid and steroid, multicellular organismal homeostasis and extracellular matrix constituent secretion.ConclusionsMultiple biological processes and signaling pathways are immediately activated in remaining dorsal lobes of female zebrafish right after removal of the ventral lobe and these findings provide crucial clues for further identification of cis-elements and trans-factors that are extensively involved in the initiation of liver compensatory growth.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-2145-5) contains supplementary material, which is available to authorized users.

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.

Wild sex in zebrafish: loss of the natural sex determinant in domesticated strains.

Sex determination can be robustly genetic, strongly environmental, or genetic subject to environmental perturbation. The genetic basis of sex determination is unknown for zebrafish (Danio rerio), a model for development and human health. We used RAD-tag population genomics to identify sex-linked polymorphisms. After verifying this "RAD-sex" method on medaka (Oryzias latipes), we studied two domesticated zebrafish strains (AB and TU), two natural laboratory strains (WIK and EKW), and two recent isolates from nature (NA and CB). All four natural strains had a single sex-linked region at the right tip of chromosome 4, enabling sex genotyping by PCR. Genotypes for the single nucleotide polymorphism (SNP) with the strongest statistical association to sex suggested that wild zebrafish have WZ/ZZ sex chromosomes. In natural strains, "male genotypes" became males and some "female genotypes" also became males, suggesting that the environment or genetic background can cause female-to-male sex reversal. Surprisingly, TU and AB lacked detectable sex-linked loci. Phylogenomics rooted on D. nigrofasciatus verified that all strains are monophyletic. Because AB and TU branched as a monophyletic clade, we could not rule out shared loss of the wild sex locus in a common ancestor despite their independent domestication. Mitochondrial DNA sequences showed that investigated strains represent only one of the three identified zebrafish haplogroups. Results suggest that zebrafish in nature possess a WZ/ZZ sex-determination mechanism with a major determinant lying near the right telomere of chromosome 4 that was modified during domestication. Strains providing the zebrafish reference genome lack key components of the natural sex-determination system but may have evolved variant sex-determining mechanisms during two decades in laboratory culture.

Erratum: Corrigendum: The zebrafish reference genome sequence and its relationship to the human genome

Nature 496, 498–503 (2013); doi:10.1038/nature12111 In this Letter, five authors were inadvertently omitted: Sharmin Begum and Christine Lloyd from the Wellcome Trust Sanger Institute, and Christa Lanz, Gunter Raddatz and Stephan C. Schuster from the Max Planck Institute for Developmental Biology. David Elliot was incorrectly listed as David Eliot, Beverley Mortimore was incorrectly listed as Beverly Mortimer, and James D.

Transcriptomic characterization of cold acclimation in larval zebrafish

BackgroundTemperature is one of key environmental parameters that affect the whole life of fishes and an increasing number of studies have been directed towards understanding the mechanisms of cold acclimation in fish. However, the adaptation of larvae to cold stress and the cold-specific transcriptional alterations in fish larvae remain largely unknown. In this study, we characterized the development of cold-tolerance in zebrafish larvae and investigated the transcriptional profiles under cold stress using RNA-seq.ResultsPre-exposure of 96 hpf zebrafish larvae to cold stress (16°C) for 24 h significantly increased their survival rates under severe cold stress (12°C). RNA-seq generated 272 million raw reads from six sequencing libraries and about 92% of the processed reads were mapped to the reference genome of zebrafish. Differential expression analysis identified 1,431 up- and 399 down-regulated genes. Gene ontology enrichment analysis of cold-induced genes revealed that RNA splicing, ribosome biogenesis and protein catabolic process were the most highly overrepresented biological processes. Spliceosome, proteasome, eukaryotic ribosome biogenesis and RNA transport were the most highly enriched pathways for genes up-regulated by cold stress. Moreover, alternative splicing of 197 genes and promoter switching of 64 genes were found to be regulated by cold stress. A shorter isoform of stk16 that lacks 67 amino acids at the N-terminus was specifically generated by skipping the second exon in cold-treated larvae. Alternative promoter usage was detected for per3 gene under cold stress, which leading to a highly up-regulated transcript encoding a truncated protein lacking the C-terminal domains.ConclusionsThese findings indicate that zebrafish larvae possess the ability to build cold-tolerance under mild low temperature and transcriptional and post-transcriptional regulations are extensively involved in this acclimation process.

A systematic genome-wide analysis of zebrafish protein-coding gene function

Since the publication of the human reference genome, the identities of specific genes associated with human diseases are being discovered at an enormous rate. A central problem is that the biological activity of these genes is often unclear. Detailed investigations in vertebrate model organisms, typically mice, have been essential for understanding the activities of many orthologues of these disease-associated genes. Although gene-targeting approaches1-3 and phenotype analysis have led to a detailed understanding of nearly 6,000 protein-coding genes3,4, this number falls significantly short of all >22,000 mouse protein-coding genes5. Similarly, in zebrafish genetics, one-by-one gene studies using positional cloning6, insertional mutagenesis7-9, antisense morpholino oligonucleotides10, targeted re-sequencing11-13 and zinc finger and TAL endonucleases14-17 have made significant contributions to our understanding of the biological activity of vertebrate genes, but the number of genes studied again falls well short of the >26,000 zebrafish protein-coding genes18. Importantly, for both mice and zebrafish, none of these strategies is particularly suited to the rapid generation of knockouts in thousands of genes and the assessment of their biological activity. Enabled by a well-annotated zebrafish reference genome sequence18,19, high-throughput sequencing and efficient chemical mutagenesis, we describe an active project that aims to identify and phenotype disruptive mutations in every zebrafish protein-coding gene. Thus far we have identified potentially disruptive mutations in more than 38% of all known protein coding genes. We have developed a multi-allelic phenotyping scheme to efficiently assess the effects of each allele during embryogenesis and have analysed the phenotypic consequences of over 1000 alleles. All mutant alleles and data are available to the community and our phenotyping scheme is adaptable to phenotypic analysis beyond embryogenesis.

A Sequence-Based Variation Map of Zebrafish

Zebrafish (Danio rerio) is a popular vertebrate model organism largely deployed using outbred laboratory animals. The nonisogenic nature of the zebrafish as a model system offers the opportunity to understand natural variations and their effect in modulating phenotype. In an effort to better characterize the range of natural variation in this model system and to complement the zebrafish reference genome project, the whole genome sequence of a wild zebrafish at 39-fold genome coverage was determined. Comparative analysis with the zebrafish reference genome revealed approximately 5.2 million single nucleotide variations and over 1.6 million insertion-deletion variations. This dataset thus represents a new catalog of genetic variations in the zebrafish genome. Further analysis revealed selective enrichment for variations in genes involved in immune function and response to the environment, suggesting genome-level adaptations to environmental niches. We also show that human disease gene orthologs in the sequenced wild zebrafish genome show a lower ratio of nonsynonymous to synonymous single nucleotide variations.

Extensive genetic diversity and substructuring among zebrafish strains revealed through copy number variant analysis

Copy number variants (CNVs) represent a substantial source of genomic variation in vertebrates and have been associated with numerous human diseases. Despite this, the extent of CNVs in the zebrafish, an important model for human disease, remains unknown. Using 80 zebrafish genomes, representing three commonly used laboratory strains and one native population, we constructed a genome-wide, high-resolution CNV map for the zebrafish comprising 6,080 CNV elements and encompassing 14.6% of the zebrafish reference genome. This amount of copy number variation is four times that previously observed in other vertebrates, including humans. Moreover, 69% of the CNV elements exhibited strain specificity, with the highest number observed for Tubingen. This variation likely arose, in part, from Tubingen's large founding size and composite population origin. Additional population genetic studies also provided important insight into the origins and substructure of these commonly used laboratory strains. This extensive variation among and within zebrafish strains may have functional effects that impact phenotype and, if not properly addressed, such extensive levels of germ-line variation and population substructure in this commonly used model organism can potentially confound studies intended for translation to human diseases.

When Technology Meets Conservation: Increased Microsatellite Marker Production Using 454 Genome Sequencing on the Endangered Okaloosa Darter (Etheostoma okaloosae)

High-throughput sequencing affords a cost and time effective means of obtaining large numbers of genetic markers for conservation studies. Here, we present thousands of novel microsatellite loci developed for the Okaloosa darter, Etheostoma okaloosae, an endangered percid fish. We sequenced more than 29 million bp using 454 whole genome shotgun sequencing and employed free user-friendly bioinformatics tools to screen for microsatellite loci and design appropriate primers. We tested 39 primer sets for polymorphism and ran population-level analyses on a population of Okaloosa darters. Of these, 30 markers were variable with an observed and expected heterozygosity of 0.382 and 0.430, respectively, and allele numbers ranging from 2 to 13. Comparisons against the zebra fish reference genome, Danio rerio, revealed that these loci represent an adequate chromosomal coverage of the darter genome, although total genomic coverage was only 2.4-3.3%. We also tested these loci on the brown darter, E. edwini, and identified loci that will be useful for hybridization studies between these taxa.