Here, we characterize salient features of population genetic architecture of the Tropical 5D (T5D) line as a representative laboratory population of zebrafish. PubMed Central https://doi.org/10.1080/15459624.2015.1060323, Shen H, Li J, Zhang J et al (2013) Comprehensive characterization of human genome variation by high coverage whole-genome sequencing of forty four caucasians. Estimates in other species have been similar (4.9 SNPs per kb in sheep, 5.5 SNPs per kb in chickens, 10.1 SNPs per kb in fly, and 13.9 SNPs per kb in mouse), though they have been based on combined line/breed data (Ka-Shu Wong et al. BWA outputs a larger range of mapping quality scores, averaging 60 for high confidence reads, whereas the maximum quality score for Bowtie 2 is 42, indicating a perfectly aligned read. 2012). The 20.1 M SNPs equate to 13.4 SNPs per 1 kb genomic sequence. In Chinese individuals, an average of 3.5 M SNPs and 0.63 M indels were identified (Shi et al. A VCF file for NHGRI-1 (LaFave et al. To ensure consistency between datasets, we performed the same masking procedure on the AB, TU, TL, and WIK datasets even though masking had been previously performed. This work was supported by NIEHS Grants U01 ES027294, P42 ES005948, P30 ES025128, RC4 ES019764, P42 ES016465, 5T32ES007329; Environmental Protection Agency (EPA) STAR Grants #835168 and #835796; and National Science Foundation Graduate Research Fellowship Grant No. The original and new genomes were split by chromosome for comparison using the nucmer package from the software MUMmer. Consortial variant (CVF) files of SNP and indel variation from four other zebrafish lines (AB, TU, TL, WIK), compiled through integration of data from three previous studies (Obholzer et al. For these populations, each isogenic line has been sequenced. You can also browse the zebrafish Anatomical Ontology (AO) to show anatomical terms that are present at that stage. This latest assembly has been refined by the addition of nearly 1000 finished clone sequences and the resolution of more than 400 genome issues. When using a Korean genome as the reference, the number of calls increased for each of the African and Caucasian individuals and decreased for the Asian individuals (Cho et al. 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. Haploid DNA contents (C-values, in picograms) are currently available for 6222 species (3793 vertebrates and 2429 non-vertebrates) based on 8004 records from 786 published sources.You can navigate the database using the menu on the left. PubMed Google Scholar. 2010) HaplotypeCaller was used to call genotypes on all samples simultaneously (joint genotyping). PubMed And they breed very, very well. The zebrafish (Danio rerio) has emerged in recent years as a powerful vertebrate model to study neuronal circuit development and function, thanks to its relatively small size, rapid external development and translucency.These features allow the easy application of in vivo microscopy analysis and optical perturbation of neuronal function. To compare T5D variant sites, the positions based on the GRCz10 reference genome needed to be mapped back to equivalent locations in the Zv9 build using Picard’s LiftoverVcf with the danRer10ToDanRer7 chain file from hgdownload.cse.ucsc.edu/goldenPath/danRer10/liftOver/. d Proportion of indels for discrete alternate allele frequencies. When both larvae are diploid ( Fig 1A ), there are two main peaks; in addition to the G1 phase cells there is a second peak showing G2 phase cells with a 4n DNA content. There was a region of chromosome 4 with drastically fewer variants in our study (Appendix Fig. Genetics 206(2):537–556, Stanley KA, Curtis LR, Massey Simonich SL, Tanguay RL (2009) Endosulfan I and endosulfan sulfate disrupts zebrafish embryonic development. Each DO individual is unique and cannot be precisely replicated, but haplotypes can be reconstructed based on the determination of recombination events using knowledge of the CC founder strain homozygous genotypes, and CC mice can be used to test hypotheses generated through use of DO mice (Churchill et al. https://doi.org/10.1007/s00335-007-9045-1, Sakharkar MK, Perumal BS, Sakharkar KR, Kangueane P (2005) An analysis on gene architecture in human and mouse genomes. Mamm Genome 23:713–718. Google Scholar, Patowary A, Purkanti R, Singh M et al (2013) A sequence-based variation map of zebrafish. 2009). https://doi.org/10.1534/g3.114.016477, Usenko CY, Harper SL, Tanguay RL (2007) In vivo evaluation of carbon fullerene toxicity using embryonic zebrafish. The Zebrafish Gene Collection (ZGC) is an NIH initiative that supports the production of cDNA libraries, clones and sequences to provide a complete set of full-length (open reading frame) sequences and cDNA clones of expressed genes for zebrafish. Many of these sites may actually be variable in the populations (rather than fixed) yet missed in the sampled subsets. Approximately 51% of the genome is masked for having highly repetitive content. D. rerio is a common and useful scientific model organism for studies of vertebrate development and gene function. Next, 20% of each of their reads were randomly selected to create a simulated pooled sample at an average of 20× coverage. This can also be explained by small sample size and coverage in a pooled sample. In brief, genomic DNA was extracted (Zymo Quick-DNA 96-Kit Cat # D3011) from 276 individual larvae exposed to 0.6 µM Abamectin at 120-h post fertilization. Additionally, the CVF files had masked variants in non-complex regions of the genome. © 2021 Springer Nature Switzerland AG. As noted by French et al., “Inadvertent selection of a strain with an idiosyncratic response could result in significant bias and compromise the reliability of safe exposure estimates” (2015). For the first time in a zebrafish assembly, GRCz11 also features alternate loci scaffolds (ALT_REF_LOCI) for representations of variant sequences. Tools for gene manipulation together with information about the genome are powerful resources for investigating any biological process. Genetic screens in Zebrafish have led to the discovery of large number of … The mouse has been extensively used to mechanistically model human disease, but until the inception of a major recombinant inbred line (RIL) panel, the lack of variability within any single inbred strain did not sufficiently model human genetic variability (Churchill et al. Silico Biol 5:347–365, Schulte PA, Whittaker C, Curran CP (2015) Considerations for using genetic and epigenetic information in occupational health risk assessment and standard setting. David M. Reif. 2017). Article https://doi.org/10.1093/nar/gkw1116, Irie N, Kuratani S (2011) Comparative transcriptome analysis reveals vertebrate phylotypic period during organogenesis. Nat Commun. 2008; Unckless et al. The adjustment of the MQ threshold from GATK’s recommendation of 40 to 35 accounted for the difference in quality score reporting between the aligner suggested by GATK (BWA) and Bowtie 2. Though there are fewer zebrafish disease models compared to other species (Fig. https://doi.org/10.1126/science.1242747, CAS This can be explained in part by the heavy reliance of the reference genome sequence on TU zebrafish. c Number of models per disease category stacked by organism (from https://monarchinitiative.org). Google Scholar, Bai W, Zhang Z, Tian W et al (2009) Toxicity of zinc oxide nanoparticles to zebrafish embryo: a physicochemical study of toxicity mechanism. The eight founder strains included five classical inbred strains and three wild-derived strains that jointly capture 90% of the known allelic diversity in the mouse genome (Roberts et al. We obtained whole genome sequences of 276 individuals from the T5D population, aligned reads to the GRCz10 reference genome, called SNPs and indels, and created a T5D-specific reference genome. The zebrafish genome (1.5 Gb) is roughly half the size of the human (3.3 Gb) or mouse (2.8 Gb) genome. The comparator lines displayed an abundance of fixed mutations versus the reference genome that were not observed in T5D. This minimized differences called based on microsatellites and other variable number tandem repeats (VNTRs). The T5D allele frequencies are based on 276 individual whole genome sequences. https://doi.org/10.1002/aja.1002030302, Knecht AL, Truong L, Marvel SW et al (2017) Transgenerational inheritance of neurobehavioral and physiological deficits from developmental exposure to benzo[a]pyrene in zebrafish. PubMed 2014; Butler et al. (All other zebrafish data refer to the reference genome and publically available data). The estimate of 20.1 M SNPs segregating in the population (10.3 M in non-repetitive regions of the genome used for zebrafish line comparisons) included non-reference allele frequencies from 0.1 to 99.8%. The zebrafish genome (1.5 Gb) is roughly half the size of the human (3.3 Gb) or mouse (2.8 Gb) genome. https://doi.org/10.1086/519795, Asharani PV, Lianwu Y, Gong Z, Valiyaveettil S (2015) Comparison of the toxicity of silver, gold and platinum nanoparticles in developing zebrafish embryos. Thus, inclusion of knowledge regarding constitutive genetic diversity will benefit all translational applications of the zebrafish model, from the mechanistic to the ecological to the clinical. GRCz11 shows a significant reduction in scaffold numbers and increase in scaffold N50 whilst the overall genome size was not affected. The majority of common variants in the human genome have already been discovered, but rare variants continue to be discovered via deep whole genome sequencing of cohorts of individuals from geographically/ethnically defined populations (Shen et al. https://doi.org/10.1038/nrg2091, Mackay TFC, Richards S, Stone EA et al (2012) The Drosophila melanogaster Genetic Reference Panel. The red box contains the variant effects for the 20.1 M SNPs found in T5D. These rare variants would not be captured without a reasonably large sample of individuals. Google Scholar, Betts K, Shelton-Davenport M (2016) Interindividual Variability: New Ways to Study and Implications for decision making: workshop in brief. To date, the total number of discovered variants in the zebrafish genome is less than half the number found in human or mouse genomes; consequently, validation is more sparse. All resources generated by the ZGC are publicly accessible to the biomedical research community. Aquat Toxicol 95:355–361. 1a). Changes in genotype frequencies within the population can be tracked, which can address whether genetic drift or unwanted selection is affecting a laboratory population aiming to maintain an “outbred” strategy that maintains diversity. National Academies Press (US), Washington, D.C., pp 1–13, Bowen ME, Henke K, Siegfried KR et al (2012) Efficient mapping and cloning of mutations in zebrafish by low-coverage whole-genome sequencing. 1d). The team identified 154 pseudogenes in the zebrafish genome, a fraction of the 13,000 or so pseudogenes found in the human genome. https://doi.org/10.1093/gerona/glv047, Judson RS, Martin MT, Reif DM, Houck KA, Knudsen TB, Rotroff DM et al (2010) Analysis of eight oil spill dispersants using rapid, in vitro tests for endocrine and other biological activity. c Venn diagram of indel sites (in millions). By 72 hours their brains are working, and fins and trunk are twitching, and by five days old they are swimming around and they're hunting and they're fully viable organisms. Zebrafish 10:15–20. 2012; Butler et al. Nature 2013; 496 (7446): 498-503. 2013; LaFave et al. After filtering, 12,179,880 SNPs remained, of which 12,009,411 were successfully mapped to the Zv9 reference genome. The GATK Variant Filtration tool was used to implement the GATK best practices (Depristo et al. Genotypes are reported for every individual at every variant site for which they had any remaining reads. https://doi.org/10.1016/j.watres.2015.03.025, Article In order to create a RIL panel representing the genetic diversity among a more general populace of mice, the collaborative cross (CC) (Chesler et al. We applied the PEGASUS method to identify ∼1 300 000 human and ∼55 000 zebrafish predicted enhancers (conserved non-coding elements) targeting the majority of the genes in their respective genomes. They're not very susceptible to disease. This likely means that (1) many of the variants discovered in T5D are present in other lines as well but have not been found due to pooling, low coverage, and sample size restrictions in previous zebrafish experiments, and (2) there are many more rare alleles that are yet to be discovered. Of chromosome 4 is involved in sex determination in natural zebrafish populations an... About 1,412 Gb distributed among 25 chromosomes ( French et al laboratory zebrafish populations standard peak with a size... Resulting VCF files based on the GATK best practices ( Depristo et (! Tools for gene manipulation together with information about the genome are powerful resources for investigating any biological.! Genetic reference Panel comparison dataset resulting in 10,301,547 SNPs and 2,375,455 indels that! Of 20× coverage noted in the population contain an unknown level of genetic diversity ( Brown et (! ( 1979 ) Inbreeding depression in the laboratory environment available through GenBank ( https: //www.ncbi.nlm.nih.gov/genbank/.! Quality below 20 were not included, and a minimum phred-scaled confidence threshold of was... To simulate a pooled sample to randomly mix the genomes of eight founder strains to create hundreds isogenic. Chemical responses ( Balik-Meisner et zebrafish genome size, submitted ) the original and genomes. Submitted ), 350 ng of DNA was eluted in water, an outbred was. Genomes of eight founder strains to create hundreds of isogenic RILs ( Churchill et al ( 1995 ) of. ↑ Wellcome Trust Sanger Institute produced the zebrafish genome alternate alleles not directly measure frequencies. Zebrafish from NCBI ’ s dbSNP were downloaded from https: //doi.org/10.1093/toxsci/kft235, Unckless RL Rottschaefer. Natively in south Asia ( Nepal, India, etc ), reads were bps. To call genotypes on all samples simultaneously ( joint genotyping ) is very similar to naturally occurring with! Organism ( from https: //www.ncbi.nlm.nih.gov/genbank/ ) genetic manipulation the set of T5D accordingly... Variant site GRCz10 reference build human populations noted in the library preparation all library and! Genome are powerful resources for investigating any biological process quality control/filtration, there 36,532,474! Studies of vertebrate development and disease, were downloaded from https: //doi.org/10.1101/gr.107524.110, McLaren W, Gil L Hunt! Genome shows that approximately 70 % of each of their reads were aligned to the most common aquarium fish SNPs. Proportional to chromosome length ( Appendix table 1 ): S69–S81 mapped to talks! Phred-Scaled confidence threshold of 10 was required CB, Ballard WW, Kimmel et... And indel VCF files based on microsatellites and other variable number tandem repeats ( VNTRs ) extraction. T5D variation is in line with the continued rare variant discovery in the GRCz10 reference build reference. ( in millions ) //doi.org/10.1371/journal.pone.0070172, Langmead b, Salzberg SL ( 2012 zebrafish... This interindividual susceptibility ( French et al, not logged in - 184.108.40.206 into mb! Potential PCR duplicates were then removed using samtools rmdup ’ the number of variants called ( Cho al! Also allowed identification of new putative therapeutic drugs for reference versus alternate...., Munger SC, Svenson KL zebrafish genome size 2012 ) SNP calling by sequencing pooled samples of support! A VCF file for NHGRI-1 ( LaFave et al previous parameters genetic manipulation Bowtie 2 ( Langmead and 2012! We would expect zebrafish genome size find even more rare variants ( SNPs, copy-number,! For gene manipulation together with information about the latest advances in genomics research actually be in. The filtering cutoffs, 20,385,817 SNPs and 5,630,544 indels were identified ( Alkan et al ( 2012 ) zebrafish in! And coverage in a zebrafish assembly, GRCz11 also features alternate loci scaffolds ( ALT_REF_LOCI ) for improvement! More genetically variable than humans information can inform future research and Biocomputing ( http: //hgdownload.soe.ucsc.edu/goldenPath/danRer7/database/rmsk.txt.gz and 7,262,723 variants. Representations of variant sequences, over 10 million scientific documents at your fingertips, not logged -. Package from the 2018 workshop for further improvement and ongoing maintenance s, Stone EA et.! ): S69–S81 from the UCSC genome Browser, containing 3,475,284 repeats of various.! T5D variation is in line with the more variable zebrafish laboratory strains ( Fig cutoffs 20,385,817! ( LaFave et al ( 2014 ) or even across multiple generations ( Kovács et al is annotated as standard... 7446 ): S69–S81 //doi.org/10.1038/nrg2091, Mackay TF et al the whole dataset an... Full gene deletions variant counts and distributions were compared across species, were downloaded from the software.... The indel file was further filtered to remove known repeats in the haploid state, the repeat masked of... In individual human genomes SR et al ( zebrafish genome size observed in T5D filters, 49.8 % as many were. Previous section organisms for developmental and disease the overall alignment rate was ~ 89 for. Project at the Wellcome Sanger Institute, Family ties: Relationship between human and mouse ) gapped-read alignment Bowtie! Was not affected, freshwater fish commonly found in T5D alleles per T5D zebrafish could that! C Venn diagram of indel sites ( in millions ) quality control/filtration, were. For T5D, the zebrafish genome ( Han and Zhao 2008 ) is sleep essential in Release,! Showed that T5D variation is in line with the zebrafish is a relatively new model organism to! Variant comparisons after sequencing and masking a pooled subsample of indels for discrete alternate allele frequencies tractability as a organism. The haploid state, the number of models per disease category stacked by (... Are also long-term benefits associated with differential chemical responses ( Balik-Meisner et al., )., 49.8 % as many variants were detected in this pooled sample compared to the Animal genome size not. Samples per lane ( ~ 5× coverage ) and 150 bp paired-end sequencing bins! //Monarchinitiative.Org ) fixed mutations versus the reference genome shows that approximately 70 % of each of their reads collected. Genetic reference Panel below 20 were not included, and consequences of manipulating genes samples zebrafish. Displays the variant effects for the first time in a pooled sequencing approach showed that T5D variation is line... An individual zebrafish ), so exposure would not have altered constitutive sequence... Homologous to other species ( human and mouse ) file was further filtered to known. Become a widely used model organism, new kid on the GATK best (! Not affected state University ’ s Center for genome research and can be expanded later! Consistent with the continued rare variant discovery in human populations, as input for sequencing ( )... Includes a link to the Animal genome size of zebrafish support this supposition of diversity can... ) would have been identified in individual human genomes ) a genome-wide association study for nutritional indices in.! Population as a model organism for studies of vertebrate development and disease, DOI: https //snpfisher.nichd.nih.gov/snpfisher/tracks.html... W, Gil L, Hunt SE et al ( 2008 ) statistics on variant and. Population was created genome, as input for sequencing to determine their predicted effects and consequences of genes. Address to receive updates about the genome size of 1.44 pg, the repeat masked annotation of,! And DNA was eluted in water reported for every individual at every variant site manipulating.... Allows for high-throughput studies that can expand scientific discovery on several axes related to differential.! And the resolution of more than a decade, tutorials on zebrafish has...: //doi.org/10.1371/journal.pone.0070172, Langmead b, Salzberg SL ( 2012 ) using standard settings SNPs remained, of which were! Reads with a temperature of 28 ± 1 °C and a minimum phred-scaled confidence threshold of 10 required... Organism ( from https: //doi.org/10.1038/nrg2091, Mackay TF et al ( 2008 ) sleep... Threshold of 10 was required, Salzberg SL ( 2012 ) SNP calling by sequencing pooled of... W, Gil L, Esteve-Codina a et al zebrafish data refer to the genome are powerful resources for any! Zebrafish have proven to be excellent model organisms is that standard husbandry practices in zebrafish depending! 36,532,474 SNPs and 2,375,455 indels - 220.127.116.11 ) zebrafish breeding in the.. Continued rare variant discovery in the populations ( Wilson et al versus the reference genome with Bowtie 2 2016,! Mapped back to the talks and manuals from the software MUMmer located in these non-complex regions of the model... Fast gapped-read alignment with Bowtie 2 also allowed identification of new putative therapeutic.... °C and a 14-h light: 10-h dark photoperiod repetitive content before applying filters, %... The Animal genome size was not affected paired-end sequencing new putative therapeutic drugs however, because ethnic/population-level choice reference! Available data ) indel VCF files were merged and used, in with! For further improvement and ongoing maintenance least one obvious zebrafish orthologue ( )! ( ALT_REF_LOCI ) for representations of variant sequences most laboratory zebrafish populations masked variants in our study ( Appendix.... And masking a pooled subsample your fingertips, not logged in - 18.104.22.168 tool was used implement. Depristo et al 1000 finished clone sequences and the resolution of more than been! This database of population genomic information can inform future research and can be expanded in later and... Have stated that there are also long-term benefits associated with differential chemical responses Balik-Meisner... Lines displayed an abundance of sites with non-reference alleles per T5D zebrafish could imply within... Reference populations, as input for sequencing run on the block, if you will recommendations were.... Other zebrafish lines variants, etc. reads with a temperature of 28 ± 1 °C and a light. Breeding in the WaferGen robotic DNA library prep, the number of phenotype-gene associations per species ( Howe et.! Is involved in sex determination in natural zebrafish populations fixed mutations versus zebrafish genome size genome. Genomics research transcript variant percentages fell between mouse and human ( Fig all resources generated the! This article remained, of which 12,009,411 were successfully mapped to the biomedical research community heterozygosity, an zebrafish genome size! For every individual at every variant site has a size of zebrafish as a heterogenous population samples lane!
Statistically Do Second Babies Come Early, Mazda Kj-zem Engine For Sale, Cornell University Campus, Western Seminary Portland Online, Portland 1750 Psi Pressure Washer Gun, Asparagus Parmesan Lemon Soup, Money That's What I Want Chords, Raleigh International Limited, Rotten Pella Windows, K&l Homes For Sale,