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.

Author Correction: Assembly and Analysis of Unmapped Genome Sequence Reads Reveal Novel Sequence and Variation in Dogs.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Assembly and Analysis of Unmapped Genome Sequence Reads Reveal Novel Sequence and Variation in Dogs.

Dogs are excellent animal models for human disease. They have extensive veterinary histories, pedigrees, and a unique genetic system due to breeding practices. Despite these advantages, one factor limiting their usefulness is the canine genome reference (CGR) which was assembled using a single purebred Boxer. Although a common practice, this results in many high-quality reads remaining unmapped. To address this whole-genome sequence data from three breeds, Border Collie (n = 26), Bearded Collie (n = 7), and Entlebucher Sennenhund (n = 8), were analyzed to identify novel, non-CGR genomic contigs using the previously validated pseudo-de novo assembly pipeline. We identified 256,957 novel contigs and paired-end relationships together with BLAT scores provided 126,555 (49%) high-quality contigs with genomic coordinates containing 4.6 Mb of novel sequence absent from the CGR. These contigs close 12,503 known gaps, including 2.4 Mb containing partially missing sequences for 11.5% of Ensembl, 16.4% of RefSeq and 12.2% of canFam3.1+ CGR annotated genes and 1,748 unmapped contigs containing 2,366 novel gene variants. Examples for six disease-associated genes (SCARF2, RD3, COL9A3, FAM161A, RASGRP1 and DLX6) containing gaps or alternate splice variants missing from the CGR are also presented. These findings from non-reference breeds support the need for improvement of the current Boxer-only CGR to avoid missing important biological information. The inclusion of the missing gene sequences into the CGR will facilitate identification of putative disease mutations across diverse breeds and phenotypes.

Baseline mRNA expression differs widely between common laboratory strains of zebrafish.

Common strains of wildtype zebrafish (Danio rerio) have unique genomic features including SNPs and CNV, but strain information often goes unreported in the literature. As a result, the confounding effects of interstrain variation makes repetition of studies in zebrafish challenging. Here we analyze hepatic mRNA expression patterns between three common zebrafish strains (AB, Tuebingen (TU), and WIK) using Agilent 4 × 44 K gene expression microarrays to establish baseline mRNA expression across strains and between sexes. We observed wide variation in sex-specific gene expression within AB and WIK strains (141 genes in AB and 67 genes in WIK), but no significant variation between sexes within TU. After partitioning the dataset into male and female subsets, we detected 421 unique mRNA transcripts with statistically significant differential expression; 269 mRNA transcripts varied between males, 212 mRNA transcripts varied between females, and 59 mRNA transcripts varied across the three strains, regardless of sex. It is not surprising that mRNA expression profiles differ between sexes and strains, but it is imperative to characterize the differences. These results highlight the complexity of variation within zebrafish and underscore the value of this model system as a valid representation of normal variation present in other species, including humans.

Anchored pseudo-de novo assembly of human genomes identifies extensive sequence variation from unmapped sequence reads.

The human genome reference (HGR) completion marked the genomics era beginning, yet despite its utility universal application is limited by the small number of individuals used in its development. This is highlighted by the presence of high-quality sequence reads failing to map within the HGR. Sequences failing to map generally represent 2-5 % of total reads, which may harbor regions that would enhance our understanding of population variation, evolution, and disease. Alternatively, complete de novo assemblies can be created, but these effectively ignore the groundwork of the HGR. In an effort to find a middle ground, we developed a bioinformatic pipeline that maps paired-end reads to the HGR as separate single reads, exports unmappable reads, de novo assembles these reads per individual and then combines assemblies into a secondary reference assembly used for comparative analysis. Using 45 diverse 1000 Genomes Project individuals, we identified 351,361 contigs covering 195.5 Mb of sequence unincorporated in GRCh38. 30,879 contigs are represented in multiple individuals with ~40 % showing high sequence complexity. Genomic coordinates were generated for 99.9 %, with 52.5 % exhibiting high-quality mapping scores. Comparative genomic analyses with archaic humans and primates revealed significant sequence alignments and comparisons with model organism RefSeq gene datasets identified novel human genes. If incorporated, these sequences will expand the HGR, but more importantly our data highlight that with this method low coverage (~10-20×) next-generation sequencing can still be used to identify novel unmapped sequences to explore biological functions contributing to human phenotypic variation, disease and functionality for personal genomic medicine.

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.

Comparative Analysis of the Shared Sex-Determination Region (SDR) among Salmonid Fishes.

Salmonids present an excellent model for studying evolution of young sex-chromosomes. Within the genus, Oncorhynchus, at least six independent sex-chromosome pairs have evolved, many unique to individual species. This variation results from the movement of the sex-determining gene, sdY, throughout the salmonid genome. While sdY is known to define sexual differentiation in salmonids, the mechanism of its movement throughout the genome has remained elusive due to high frequencies of repetitive elements, rDNA sequences, and transposons surrounding the sex-determining regions (SDR). Despite these difficulties, bacterial artificial chromosome (BAC) library clones from both rainbow trout and Atlantic salmon containing the sdY region have been reported. Here, we report the sequences for these BACs as well as the extended sequence for the known SDR in Chinook gained through genome walking methods. Comparative analysis allowed us to study the overlapping SDRs from three unique salmonid Y chromosomes to define the specific content, size, and variation present between the species. We found approximately 4.1 kb of orthologous sequence common to all three species, which contains the genetic content necessary for masculinization. The regions contain transposable elements that may be responsible for the translocations of the SDR throughout salmonid genomes and we examine potential mechanistic roles of each one.

Cross-species array comparative genomic hybridization identifies novel oncogenic events in zebrafish and human embryonal rhabdomyosarcoma.

Human cancer genomes are highly complex, making it challenging to identify specific drivers of cancer growth, progression, and tumor maintenance. To bypass this obstacle, we have applied array comparative genomic hybridization (array CGH) to zebrafish embryonal rhabdomyosaroma (ERMS) and utilized cross-species comparison to rapidly identify genomic copy number aberrations and novel candidate oncogenes in human disease. Zebrafish ERMS contain small, focal regions of low-copy amplification. These same regions were commonly amplified in human disease. For example, 16 of 19 chromosomal gains identified in zebrafish ERMS also exhibited focal, low-copy gains in human disease. Genes found in amplified genomic regions were assessed for functional roles in promoting continued tumor growth in human and zebrafish ERMS--identifying critical genes associated with tumor maintenance. Knockdown studies identified important roles for Cyclin D2 (CCND2), Homeobox Protein C6 (HOXC6) and PlexinA1 (PLXNA1) in human ERMS cell proliferation. PLXNA1 knockdown also enhanced differentiation, reduced migration, and altered anchorage-independent growth. By contrast, chemical inhibition of vascular endothelial growth factor (VEGF) signaling reduced angiogenesis and tumor size in ERMS-bearing zebrafish. Importantly, VEGFA expression correlated with poor clinical outcome in patients with ERMS, implicating inhibitors of the VEGF pathway as a promising therapy for improving patient survival. Our results demonstrate the utility of array CGH and cross-species comparisons to identify candidate oncogenes essential for the pathogenesis of human cancer.

Genomic amplification of an endogenous retrovirus in zebrafish T-cell malignancies.

Genomic instability plays a crucial role in oncogenesis. Somatically acquired mutations can disable some genes and inappropriately activate others. In addition, chromosomal rearrangements can amplify, delete, or even fuse genes, altering their functions and contributing to malignant phenotypes. Using array comparative genomic hybridization (aCGH), a technique to detect numeric variations between different DNA samples, we examined genomes from zebrafish (Danio rerio) T-cell leukemias of three cancer-prone lines. In all malignancies tested, we identified recurring amplifications of a zebrafish endogenous retrovirus. This retrovirus, ZFERV, was first identified due to high expression of proviral transcripts in thymic tissue from larval and adult fish. We confirmed ZFERV amplifications by quantitative PCR analyses of DNA from wild-type fish tissue and normal and malignant D. rerio T cells. We also quantified ZFERV RNA expression and found that normal and neoplastic T cells both produce retrovirally encoded transcripts, but most cancers show dramatically increased transcription. In aggregate, these data imply that ZFERV amplification and transcription may be related to T-cell leukemogenesis. Based on these data and ZFERV's phylogenetic relation to viruses of the murine-leukemia-related virus class of gammaretroviridae, we posit that ZFERV may be oncogenic via an insertional mutagenesis mechanism.

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.

Molecular cytogenetic methodologies and a BAC probe panel resource for genomic analyses in the zebrafish.

Molecular cytogenetics is a field that emerged in the 1980s, based on a technique referred to as fluorescence in situ hybridization, (FISH). Using FISH methodologies, a specific DNA sequence or collection of DNA fragments may be selectively labeled with a hapten molecule or fluorescent dye and hybridized to denatured chromosomes, interphase cells, or even chromatin fibers. DNA hybridization kinetics permit these labeled probes to anneal to their complementary sequences on such chromosomal DNA preparations allowing for direct visualization of the sequence of interest in the genome being interrogated. If present, the relative chromosomal position of the sequence can sometimes also be ascertained. Progress in molecular cytogenetic research has advanced the genetic characterization of zebrafish models of human diseases as well as assisted with accurate annotation of the zebrafish reference genome by anchoring large DNA fragments to specific chromosome regions. Using the procedures described in this chapter, hundreds of ambiguous zebrafish bacterial artificial chromosome (BAC) clones have already been assigned to individual genetic linkage groups. Molecular cytogenetic techniques can also be used to study gene duplication events and study the molecular mechanisms by which they arise. Moreover, the availability of a new molecular cytogenetic technique, array-based comparative genomic hybridization (aCGH), is now able to identify gains and losses of DNA segments in zebrafish DNA samples in a genome-wide manner and in a single assay.

Lack of a heritable reproductive defect in the offspring of male rainbow trout exposed to the environmental estrogen 17alpha-ethynylestradiol.

Endocrine disruptors, including environmental estrogens, have been shown to induce heritable effects through both genetic and epigenetic mechanisms in mammals. Despite this information and the wealth of knowledge regarding the significant reproductive impacts endocrine disruptors impose on fishes, no studies have reported whether the observed effects are heritable. Without this information it is difficult to establish the long-term consequences for exposed populations. To determine potential consequences of long-term effects we must consider the possibility that induced reproductive defects in fishes may be heritable. Using rainbow trout (Oncorhynchus mykiss) as a model this study aims to determine whether a specific reproductive defect observed in 17alpha-ethynylestradiol exposed male parents, diminished progeny survival, is heritable in the unexposed surviving F1 males. Semen was collected from anesthetized males of the F1 generation upon sexual maturation at two time-points, one year old precocious males and two years old males. In vitro fertilization was used to produce an F2 generation. F2 embryos were then analyzed for survival at 19 days post-fertilization (eye pigmentation) and the different treatment groups statistically compared to the controls. Analysis indicated that F2 offspring survival from F1 males propagated from both exposed and unexposed parents survive normally and no heritable effect was observed in males from the F1 generation for this specific reproductive defect. These results provide scope for the recovery of fish populations exposed to environmental estrogens should the contaminant be removed.

Aneuploid sperm formation in rainbow trout exposed to the environmental estrogen 17{alpha}-ethynylestradiol.

Environmental contaminants that mimic native estrogens (i.e., environmental estrogens) are known to significantly impact a wide range of vertebrate species and have been implicated as a source for increasing human male reproductive deficiencies and diseases. Despite the widespread occurrence of environmental estrogens and recognized detrimental effects on male vertebrate reproduction, no specific mechanism has been determined indicating how reduced fertility and/or fecundity is achieved. Previous studies show that male rainbow trout, Oncorhynchus mykiss, exposed to the environmental estrogen 17alpha-ethynylestradiol (EE2) before gamete formation and fertilization produce progeny with significantly reduced embryonic survival. To determine whether this observed decrease results from sperm chromosome alterations during spermatogenesis, male rainbow trout were exposed to 10 ng of EE2/l for 50 days. After exposure, semen was collected and sperm aneuploidy levels analyzed with two chromosome markers by fluorescent in situ hybridization. In vitro fertilizations were also conducted by using control and exposed sperm crossed to eggs from an unexposed female for offspring analysis. Evaluations for nucleolar organizer region number and karyotype were performed on developing embryos to determine whether sperm aneuploidy translated into embryonic aneuploidy. Results conclusively show increased aneuploid sperm formation due to EE2 exposure. Additionally, embryonic cells from propagated progeny of individuals possessing elevated sperm aneuploidy display high levels of embryonic aneuploidy. This study concludes that EE2 exposure in sexually developing male rainbow trout increases levels of aneuploid sperm, providing a mechanism for decreased embryonic survival and ultimately diminished reproductive success in EE2 exposed males.

Effect of parental exposure to trenbolone and the brominated flame retardant BDE-47 on fertility in rainbow trout (Oncorhynchus mykiss).

We exposed sexually maturing male rainbow trout (Oncorhynchus mykiss) to BDE-47 (a polybrominated diphenyl ether) and female rainbow trout to trenbolone (an anabolic steroid). Male trout were orally exposed for 17 days to 55 microg/kg/day BDE-47 and female trout continuously exposed for 60-77 days to a measured trenbolone water concentration of 35 ng/L. After the exposure, eggs and semen were collected and in vitro fertilization trials performed using a sperm:egg ratio of 300,000:1. In the BDE-47 study, eggs from control females were fertilized with semen from exposed males, while in the trenbolone study, eggs from exposed females were fertilized with semen from control males. All treatments were evaluated at two-three early developmental time-points representing first cleavage (0.5 day), embryonic keel (9 days), and eyed stages (19 days), respectively. The results indicated that BDE-47 exposure did not alter fertility as embryonic survival was similar between control and exposed groups. Trenbolone exposure also did not alter embryo survival. However, in the embryos fertilized with eggs from trenbolone exposed females, a noticeable delay in developmental progress was observed. On day 19 when eye development is normally complete, the majority of the embryos either lacked eyes or displayed under-developed eyes, in contrast to control embryos. This finding suggests steroidal androgen exposure in sexually maturing female rainbow trout can impact developmental timing of F1 offspring.

Reduced embryonic survival in rainbow trout resulting from paternal exposure to the environmental estrogen 17alpha-ethynylestradiol during late sexual maturation.

Exposure of fishes to environmental estrogens is known to affect sexual development and spawning, but little information exists regarding effects on gametes. This study evaluated embryonic survival of offspring from male rainbow trout (Oncorhynchus mykiss) exposed to 17alpha-ethynylestradiol (EE(2)) using an in vitro fertilization protocol. Males were exposed at either 1800 or 6700 degree days ( degrees d) (i.e. 161 or 587 days post-fertilization (dpf)) to test for effects on testes linked to reproductive ontogeny. At 1800 degrees d, fish were beginning testicular differentiation and were exposed to 109 ng EE(2)/l for 21 days. At 6700 degrees d, fish have testes containing spermatocytes and spermatids and were exposed for 56 days to either 0.8, 8.3, or 65 ng EE(2)/l. Semen was collected at full sexual maturity in each group and used to fertilize eggs pooled from several non-exposed females. Significant decreases in embryonic survival were observed only with the 6700 degrees d exposure. In 0.8 and 8.3 ng EE(2)/l treatments, embryo survival was significantly reduced at 19 dpf when compared with the control. In contrast, an immediate decrease in embryonic survival at 0.5 dpf was observed in the 65 ng EE(2)/l treatment. Blood samples collected at spawning from 6700 degrees d exposed males revealed a significant decrease in 11-ketotestosterone and a significant increase in luteinizing hormone levels for the 65 ng EE(2)/l treatment when compared with the other treatment groups. Results indicate that sexually maturing male rainbow trout are susceptible to EE(2) exposure with these fish exhibiting two possible mechanisms of reduced embryonic survival through sperm varying dependant on EE(2) exposure concentrations experienced.

Intraspecific variation in the rainbow trout mitochondrial DNA genome.

Four complete mitochondrial DNA genomes for rainbow trout, Oncorhynchus mykiss, were sequenced and compared to the previously reported complete mitochondrial sequence. Analyses revealed 13 additional bp and 31 insertion/deletion (indels) sites between the four new sequences and the previously published sequence. Twenty-eight indels were single base pairs while the remaining were multiple base pair insertions, two in the ATPase6 gene, one of three bases and one of nine bases, and a two base insert in the control region. Indels in the new genomes were compared to a wide range of salmonid species and corresponded with sequences in closely related species, indicating that errors were likely incorporated into the previously published sequence. Although variable nucleotide positions were observed throughout the newly sequenced genomes, non-synonymous amino acid substitutions were observed in only six of seven ND genes. While there were only 14 amino acid substitutions between genomes, K(a)/K(s) ratios ranged between 0 and 1. A single amino acid deletion in the ND1 coding region unique to rainbow trout was also detected. Using complete mitochondrial genomes for all sequenced Salmonidae species, a phylogenetic analysis was performed to establish a salmonid mitochondrial phylogeny revealing support for Salvelinus, rather than Salmo, as the sister taxa to Oncorhynchus.