Efficient mapping and cloning of mutations in zebrafish by low-coverage whole-genome sequencing

Exploring hematopoiesis in zebrafish using forward genetic screening

Zebrafish have emerged as a powerful animal model for investigating the genetic basis of hematopoiesis. Owing to its close genetic and developmental similarities to humans, combined with its rapid reproduction and extensive genomic resources, zebrafish have become a versatile and efficient platform for genetic studies. In particular, the forward genetic screening approach has enabled the unbiased identification of novel genes and pathways related to blood development, from hematopoietic stem cell formation to terminal differentiation. Recent advances in mutant gene mapping have further expanded the scope of forward genetic screening, facilitating the identification of previously unknown genes and pathways relevant to hematopoiesis. In this review, we provide an overview of the zebrafish forward screening approach for hematopoietic gene discovery and highlight the key genes and pathways identified using this method. This review emphasizes the importance of zebrafish as a model system for understanding the genetic basis of hematopoiesis and its associated disorders.

The RNA-binding protein Adad1 is necessary for germ cell maintenance and meiosis in zebrafish

The double stranded RNA binding protein Adad1 (adenosine deaminase domain containing 1) is a member of the adenosine deaminase acting on RNAs (Adar) protein family with germ cell-specific expression. In mice, Adad1 is necessary for sperm differentiation, however its function outside of mammals has not been investigated. Here, through an N-ethyl-N-nitrosourea (ENU) based forward genetic screen, we identified an adad1 mutant zebrafish line that develops as sterile males. Further histological examination revealed complete lack of germ cells in adult mutant fish, however germ cells populated the gonad, proliferated, and entered meiosis in larval and juvenile fish. Although meiosis was initiated in adad1 mutant testes, the spermatocytes failed to progress beyond the zygotene stage. Thus, Adad1 is essential for meiosis and germline maintenance in zebrafish. We tested if spermatogonial stem cells were affected using nanos2 RNA FISH and a label retaining cell (LRC) assay, and found that the mutant testes had fewer LRCs and nanos2-expressing cells compared to wild-type siblings, suggesting that failure to maintain the spermatogonial stem cells resulted in germ cell loss by adulthood. To identify potential molecular processes regulated by Adad1, we sequenced bulk mRNA from mutants and wild-type testes and found mis-regulation of genes involved in RNA stability and modification, pointing to a potential broader role in post-transcriptional regulation. Our findings suggest that the RNA regulatory protein Adad1 is required for fertility through regulation of spermatogonial stem cell maintenance in zebrafish.

Genetically engineered zebrafish as models of skeletal development and regeneration

Zebrafish (Danio rerio) are aquatic vertebrates with significant homology to their terrestrial counterparts. While zebrafish have a centuries-long track record in developmental and regenerative biology, their utility has grown exponentially with the onset of modern genetics. This is exemplified in studies focused on skeletal development and repair. Herein, the numerous contributions of zebrafish to our understanding of the basic science of cartilage, bone, tendon/ligament, and other skeletal tissues are described, with a particular focus on applications to development and regeneration. We summarize the genetic strengths that have made the zebrafish a powerful model to understand skeletal biology. We also highlight the large body of existing tools and techniques available to understand skeletal development and repair in the zebrafish and introduce emerging methods that will aid in novel discoveries in skeletal biology. Finally, we review the unique contributions of zebrafish to our understanding of regeneration and highlight diverse routes of repair in different contexts of injury. We conclude that zebrafish will continue to fill a niche of increasing breadth and depth in the study of basic cellular mechanisms of skeletal biology.

Zebrafish Models of Cardiac Disease: From Fortuitous Mutants to Precision Medicine

Heart disease is the leading cause of death worldwide. Despite decades of research, most heart pathologies have limited treatments, and often the only curative approach is heart transplantation. Thus, there is an urgent need to develop new therapeutic approaches for treating cardiac diseases. Animal models that reproduce the human pathophysiology are essential to uncovering the biology of diseases and discovering therapies. Traditionally, mammals have been used as models of cardiac disease, but the cost of generating and maintaining new models is exorbitant, and the studies have very low throughput. In the last decade, the zebrafish has emerged as a tractable model for cardiac diseases, owing to several characteristics that made this animal popular among developmental biologists. Zebrafish fertilization and development are external; embryos can be obtained in high numbers, are cheap and easy to maintain, and can be manipulated to create new genetic models. Moreover, zebrafish exhibit an exceptional ability to regenerate their heart after injury. This review summarizes 25 years of research using the zebrafish to study the heart, from the classical forward screenings to the contemporary methods to model mutations found in patients with cardiac disease. We discuss the advantages and limitations of this model organism and introduce the experimental approaches exploited in zebrafish, including forward and reverse genetics and chemical screenings. Last, we review the models used to induce cardiac injury and essential ideas derived from studying natural regeneration. Studies using zebrafish have the potential to accelerate the discovery of new strategies to treat cardiac diseases.

A combinatorial strategy to identify various types of QTLs for quantitative traits using extreme phenotype individuals in an F2 population

Theoretical and applied studies demonstrate the difficulty of detecting extremely over-dominant and small-effect genes for quantitative traits via bulked segregant analysis (BSA) in an F₂ population. To address this issue, we proposed an integrated strategy for mapping various types of quantitative trait loci (QTLs) for quantitative traits via a combination of BSA and whole-genome sequencing. In this strategy, the numbers of read counts of marker alleles in two extreme pools were used to predict the numbers of read counts of marker genotypes. These observed and predicted numbers were used to construct a new statistic, G_w, for detecting quantitative trait genes (QTGs), and the method was named dQTG-seq1. This method was significantly better than existing BSA methods. If the goal was to identify extremely over-dominant and small-effect genes, another reserved DNA/RNA sample from each extreme phenotype F₂ plant was sequenced, and the observed numbers of marker alleles and genotypes were used to calculate G_w to detect QTGs; this method was named dQTG-seq2. In simulated and real rice dataset analyses, dQTG-seq2 could identify many more extremely over-dominant and small-effect genes than BSA and QTL mapping methods. dQTG-seq2 may be extended to other heterogeneous mapping populations. The significance threshold of G_w in this study was determined by permutation experiments. In addition, a handbook for the R software dQTG.seq, which is available at https://cran.r-project.org/web/packages/dQTG.seq/index.html, has been provided in the supplemental materials for the users' convenience. This study provides a new strategy for identifying all types of QTLs for quantitative traits in an F₂ population.

Translational relevance of forward genetic screens in animal models for the study of psychiatric disease

Psychiatric disorders represent a significant burden in our societies. Despite the convincing evidence pointing at gene and gene-environment interaction contributions, the role of genetics in the etiology of psychiatric disease is still poorly understood. Forward genetic screens in animal models have helped elucidate causal links. Here we discuss the application of mutagenesis-based forward genetic approaches in common animal model species: two invertebrates, nematodes (Caenorhabditis elegans) and fruit flies (Drosophila sp.); and two vertebrates, zebrafish (Danio rerio) and mice (Mus musculus), in relation to psychiatric disease. We also discuss the use of large scale genomic studies in human populations. Despite the advances using data from human populations, animal models coupled with next-generation sequencing strategies are still needed. Although with its own limitations, zebrafish possess characteristics that make them especially well-suited to forward genetic studies exploring the etiology of psychiatric disorders.

Insights Into Central Nervous System Glial Cell Formation and Function From Zebrafish

The term glia describes a heterogenous collection of distinct cell types that make up a large proportion of our nervous system. Although once considered the glue of the nervous system, the study of glial cells has evolved significantly in recent years, with a large body of literature now highlighting their complex and diverse roles in development and throughout life. This progress is due, in part, to advances in animal models in which the molecular and cellular mechanisms of glial cell development and function as well as neuron-glial cell interactions can be directly studied in vivo in real time, in intact neural circuits. In this review we highlight the instrumental role that zebrafish have played as a vertebrate model system for the study of glial cells, and discuss how the experimental advantages of the zebrafish lend themselves to investigate glial cell interactions and diversity. We focus in particular on recent studies that have provided insight into the formation and function of the major glial cell types in the central nervous system in zebrafish.

Modulation of bioelectric cues in the evolution of flying fishes

Changes to allometry, or the relative proportions of organs and tissues within organisms, is a common means for adaptive character change in evolution. However, little is understood about how relative size is specified during development and shaped during evolution. Here, through a phylogenomic analysis of genome-wide variation in 35 species of flying fishes and relatives, we identify genetic signatures in both coding and regulatory regions underlying the convergent evolution of increased paired fin size and aerial gliding behaviors. To refine our analysis, we intersected convergent phylogenomic signatures with mutants with altered fin size identified in distantly related zebrafish. Through these paired approaches, we identify a surprising role for an L-type amino acid transporter, lat4a, and the potassium channel, kcnh2a, in the regulation of fin proportion. We show that interaction between these genetic loci in zebrafish closely phenocopies the observed fin proportions of flying fishes. The congruence of experimental and phylogenomic findings point to conserved, non-canonical signaling integrating bioelectric cues and amino acid transport in the establishment of relative size in development and evolution.

Proteolytic Restriction of Chordin Range Underlies BMP Gradient Formation

A fundamental question in developmental biology is how morphogens, such as bone morphogenetic protein (BMP), form precise signaling gradients to impart positional and functional identity to the cells of the early embryo. We combine rigorous mutant analyses with quantitative immunofluorescence to determine that the proteases Bmp1a and Tolloid spatially restrict the BMP antagonist Chordin in dorsoventral (DV) axial patterning of the early zebrafish gastrula. We show that maternally deposited Bmp1a plays an unexpected and non-redundant role in establishing the BMP signaling gradient, while the Bmp1a/Tolloid antagonist Sizzled is surprisingly dispensable. Combining computational modeling and in vivo analyses with an immobile Chordin construct, we demonstrate that long-range Chordin diffusion is not necessary for BMP gradient formation and DV patterning. Our data do not support a counter-gradient of Chordin and instead favor a Chordin sink, established by Bmp1a and Tolloid, as the primary mechanism that drives BMP gradient formation.

The BMP morphogen generates a precise signaling gradient during axial patterning. In the zebrafish embryo, Tuazon et al. find that proteases Bmp1a/Tolloid are key to this process, preventing the long-range diffusion of the BMP antagonist, Chordin. By regionally restricting Chordin, Bmp1a/Tolloid establish the signaling sink that drives BMP gradient formation.

Latent developmental potential to form limb-like skeletal structures in zebrafish

Changes in appendage structure underlie key transitions in vertebrate evolution. Addition of skeletal elements along the proximal-distal axis facilitated critical transformations, including the fin-to-limb transition that permitted generation of diverse modes of locomotion. Here, we identify zebrafish mutants that form supernumerary long bones in their pectoral fins. These new bones integrate into musculature, form joints, and articulate with neighboring elements. This phenotype is caused by activating mutations in previously unrecognized regulators of appendage patterning, vav2 and waslb, that function in a common pathway. This pathway is required for appendage development across vertebrates, and loss of Wasl in mice causes defects similar to those seen in murine Hox mutants. Concordantly, formation of supernumerary bones requires Hox11 function, and mutations in the vav2/wasl pathway drive enhanced expression of hoxa11b, indicating developmental homology with the forearm. Our findings reveal a latent, limb-like pattern ability in fins that is activated by simple genetic perturbation.

Evolution of the potassium channel gene Kcnj13 underlies colour pattern diversification in Danio fish

The genetic basis of morphological variation provides a major topic in evolutionary developmental biology. Fish of the genus Danio display colour patterns ranging from horizontal stripes, to vertical bars or spots. Stripe formation in zebrafish, Danio rerio, is a self-organizing process based on cell-contact mediated interactions between three types of chromatophores with a leading role of iridophores. Here we investigate genes known to regulate chromatophore interactions in zebrafish that might have evolved to produce a pattern of vertical bars in its sibling species, Danio aesculapii. Mutant D. aesculapii indicate a lower complexity in chromatophore interactions and a minor role of iridophores in patterning. Reciprocal hemizygosity tests identify the potassium channel gene obelix/Kcnj13 as evolved between the two species. Complementation tests suggest evolutionary change through divergence in Kcnj13 function in two additional Danio species. Thus, our results point towards repeated and independent evolution of this gene during colour pattern diversification.

Using Zebrafish to Analyze the Genetic and Environmental Etiologies of Congenital Heart Defects

Congenital heart defects (CHDs) are among the most common human birth defects. However, the etiology of a large proportion of CHDs remains undefined. Studies identifying the molecular and cellular mechanisms that underlie cardiac development have been critical to elucidating the origin of CHDs. Building upon this knowledge to understand the pathogenesis of CHDs requires examining how genetic or environmental stress changes normal cardiac development. Due to strong molecular conservation to humans and unique technical advantages, studies using zebrafish have elucidated both fundamental principles of cardiac development and have been used to create cardiac disease models. In this chapter we examine the unique toolset available to zebrafish researchers and how those tools are used to interrogate the genetic and environmental contributions to CHDs.

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.

celsr1a is essential for tissue homeostasis and onset of aging phenotypes in the zebrafish

The use of genetics has been invaluable in defining the complex mechanisms of aging and longevity. Zebrafish, while a prominent model for vertebrate development, have not been used systematically to address questions of how and why we age. In a mutagenesis screen focusing on late developmental phenotypes, we identified a new mutant that displays aging phenotypes at young adult stages. We find that the phenotypes are due to loss-of-function in the non-classical cadherin celsr1a. The premature aging is not associated with increased cellular senescence or telomere length but is a result of a failure to maintain progenitor cell populations. We show that celsr1a is essential for maintenance of stem cell progenitors in late stages. Caloric restriction can ameliorate celsr1a aging phenotypes. These data suggest that celsr1a function helps to mediate stem cell maintenance during maturation and homeostasis of tissues and thus regulates the onset or expressivity of aging phenotypes.

Whole-Genome Resequencing of Twenty Branchiostoma belcheri Individuals Provides a Brand-New Variant Dataset for Branchiostoma

As the extant representatives of the basal chordate lineage, amphioxi (including the genera Branchiostoma, Asymmetron and Epigonichthys) play important roles in tracing the state of chordate ancestry. Previous studies have reported that members of the Branchiostoma species have similar morphological phenotypic characteristics, but in contrast, there are high levels of genetic polymorphisms in the populations. Here, we resequenced 20 Branchiostomabelcheri genomes to an average depth of approximately 12.5X using the Illumina HiSeq 2000 platform. In this study, over 52 million variations (~12% of the total genome) were detected in the B. belcheri population, and an average of 12.8 million variations (~3% of the total genome) were detected in each individual, confirming that Branchiostoma is one of the most genetically diverse species sequenced to date. Demographic inference analysis highlighted the role of historical global temperature in the long-term population dynamics of Branchiostoma, and revealed a population expansion at the Greenlandian stage of the current geological epoch. We detected 594 Single nucleotide polymorphism and 148 Indels in the Branchiostoma mitochondrial genome, and further analyzed their genetic mutations. A recent study found that the epithelial cells of the digestive tract in Branchiostoma can directly phagocytize food particles and convert them into absorbable nontoxic nutrients using powerful digestive and immune gene groups. In this study, we predicted all potential mutations in intracellular digestion-associated genes. The results showed that most “probably damaging” mutations were related to rare variants (MAF<0.05) involved in strengthening or weakening the intracellular digestive capacity of Branchiostoma. Due to the extremely high number of polymorphisms in the Branchiostoma genome, our analysis with a depth of approximately 12.5X can only be considered a preliminary analysis. However, the novel variant dataset provided here is a valuable resource for further investigation of phagocytic intracellular digestion in Branchiostoma and determination of the phenotypic and genotypic features of Branchiostoma.

Integrated K+ channel and K+Cl- cotransporter functions are required for the coordination of size and proportion during development

The coordination of growth during development establishes proportionality within and among the different anatomic structures of organisms. Innate memory of this proportionality is preserved, as shown in the ability of regenerating structures to return to their original size. Although the regulation of this coordination is incompletely understood, mutant analyses of zebrafish with long-finned phenotypes have uncovered important roles for bioelectric signaling in modulating growth and size of the fins and barbs. To date, long-finned mutants identified are caused by hypermorphic mutations, leaving unresolved whether such signaling is required for normal development. We isolated a new zebrafish mutant, schleier, with proportional overgrowth phenotypes caused by a missense mutation and loss of function in the K+-Cl- cotransporter Kcc4a. Creation of dominant negative Kcc4a in wild-type fish leads to loss of growth restriction in fins and barbs, supporting a requirement for Kcc4a in regulation of proportion. Epistasis experiments suggest that Kcc4a and the two-pore potassium channel Kcnk5b both contribute to a common bioelectrical signaling response in the fin. These data suggest that an integrated bioelectric signaling pathway is required for the coordination of size and proportion during development.

Trait mapping in diverse arthropods by bulked segregant analysis

Bulked segregant analysis (BSA) is a cross-based method for genetic mapping in sexually reproducing organisms. The method's use of bulked (pooled) samples markedly reduces the genotyping effort associated with traditional linkage mapping studies. Further, it can be applied to species with life histories or physical attributes (as for micro-insects) that render genetic mapping with other methods impractical. Recent studies in both insects and mites have revealed that advanced BSA experimental designs can resolve causal loci to narrow genomic intervals, facilitating follow-up investigations. As high-quality genomes become more widely available, BSA methods are poised to become an increasingly important tool for the rapid mapping of both monogenic and polygenic traits in diverse arthropod species.

Cyclin-dependent kinase 21 is a novel regulator of proliferation and meiosis in the male germline of zebrafish

Germ cell differentiation and maintenance relies on complex regulation of mitotic and meiotic progression. Cyclin-dependent kinases (CDKs) and their activating cyclin partners are known to have specialized roles in regulating cell cycle progression across tissues, including germ cells. Very little is known about CDK/cyclin function in zebrafish or the regulation of germ cell maintenance and differentiation. In a forward genetic screen for gonadogenesis defects in zebrafish, a mutation disrupting cdk21 (cyclin-dependent kinase 21) was identified, which caused gonad hypoplasia, reduced fertility and failure of female sex specification. The cdk21 gene is unique to fishes, though the encoded protein is related to the D-cyclin partners Cdk4 and Cdk6, which are known G1 cell cycle regulators. In the testis, cdk21 mutant germ cells exhibited cell cycle defects such as diminished proliferation, prolonged meiosis and delayed sperm differentiation. Furthermore, cdk21 mutants failed to maintain germ cells following breeding. Based on these findings, we propose that cdk21 regulates spermatogonial proliferation, progression through meiosis and germline stem cell activation in the testis. In addition, we investigated cdk4 and cdk6 in zebrafish development and found that each has distinct expression patterns in the gonads. Mutant analysis demonstrated that cdk6 was necessary for viability beyond larval stages. In contrast, cdk4 mutants were viable but were all male with low breeding success and sperm overabundance. Our analysis demonstrated that zebrafish harbor three genes of the cdk4/6 family, cdk4, cdk6 and cdk21, with cdk21 having an essential role in germ cell development in the testis.

Germline Structural Variations Are Preferential Sites of DNA Replication Timing Plasticity during Development

The DNA replication timing program is modulated throughout development and is also one of the main factors influencing the distribution of mutation rates across the genome. However, the relationship between the mutagenic influence of replication timing and its developmental plasticity remains unexplored. Here, we studied the distribution of copy number variations (CNVs) and single nucleotide polymorphisms across the zebrafish genome in relation to changes in DNA replication timing during embryonic development in this model vertebrate species. We show that CNV sites exhibit strong replication timing plasticity during development, replicating significantly early during early development but significantly late during more advanced developmental stages. Reciprocally, genomic regions that changed their replication timing during development contained a higher proportion of CNVs than developmentally constant regions. Developmentally plastic CNV sites, in particular those that become delayed in their replication timing, were enriched for the clustered protocadherins, a set of genes important for neuronal development that have undergone extensive genetic and epigenetic diversification during zebrafish evolution. In contrast, single nucleotide polymorphism sites replicated consistently early throughout embryonic development, highlighting a unique aspect of the zebrafish genome. Our results uncover a hitherto unrecognized interface between development and evolution.

Conserved but flexible modularity in the zebrafish skull: implications for craniofacial evolvability

Morphological variation is the outward manifestation of development and provides fodder for adaptive evolution. Because of this contingency, evolution is often thought to be biased by developmental processes and functional interactions among structures, which are statistically detectable through forms of covariance among traits. This can take the form of substructures of integrated traits, termed modules, which together comprise patterns of variational modularity. While modularity is essential to an understanding of evolutionary potential, biologists currently have little understanding of its genetic basis and its temporal dynamics over generations. To address these open questions, we compared patterns of craniofacial modularity among laboratory strains, defined mutant lines and a wild population of zebrafish (Danio rerio). Our findings suggest that relatively simple genetic changes can have profound effects on covariance, without greatly affecting craniofacial shape. Moreover, we show that instead of completely deconstructing the covariance structure among sets of traits, mutations cause shifts among seemingly latent patterns of modularity suggesting that the skull may be predisposed towards a limited number of phenotypes. This new insight may serve to greatly increase the evolvability of a population by providing a range of 'preset' patterns of modularity that can appear readily and allow for rapid evolution.

The role of a class III gibberellin 2-oxidase in tomato internode elongation

A network of environmental inputs and internal signaling controls plant growth, development and organ elongation. In particular, the growth-promoting hormone gibberellin (GA) has been shown to play a significant role in organ elongation. The use of tomato as a model organism to study elongation presents an opportunity to study the genetic control of internode-specific elongation in a eudicot species with a sympodial growth habit and substantial internodes that can and do respond to external stimuli. To investigate internode elongation, a mutant with an elongated hypocotyl and internodes but wild-type petioles was identified through a forward genetic screen. In addition to stem-specific elongation, this mutant, named tomato internode elongated -1 (tie-1) is more sensitive to the GA biosynthetic inhibitor paclobutrazol and has altered levels of intermediate and bioactive GAs compared with wild-type plants. The mutation responsible for the internode elongation phenotype was mapped to GA2oxidase 7, a class III GA 2-oxidase in the GA biosynthetic pathway, through a bulked segregant analysis and bioinformatic pipeline, and confirmed by transgenic complementation. Furthermore, bacterially expressed recombinant TIE protein was shown to have bona fide GA 2-oxidase activity. These results define a critical role for this gene in internode elongation and are significant because they further the understanding of the role of GA biosynthetic genes in organ-specific elongation.

Zebrafish as an Emerging Model for Osteoporosis: A Primary Testing Platform for Screening New Osteo-Active Compounds

Osteoporosis is metabolic bone disease caused by an altered balance between bone anabolism and catabolism. This dysregulated balance is responsible for fragile bones that fracture easily after minor falls. With an aging population, the incidence is rising and as yet pharmaceutical options to restore this imbalance is limited, especially stimulating osteoblast bone-building activity. Excitingly, output from large genetic studies on people with high bone mass (HBM) cases and genome wide association studies (GWAS) on the population, yielded new insights into pathways containing osteo-anabolic players that have potential for drug target development. However, a bottleneck in development of new treatments targeting these putative osteo-anabolic genes is the lack of animal models for rapid and affordable testing to generate functional data and that simultaneously can be used as a compound testing platform. Zebrafish, a small teleost fish, are increasingly used in functional genomics and drug screening assays which resulted in new treatments in the clinic for other diseases. In this review we outline the zebrafish as a powerful model for osteoporosis research to validate potential therapeutic candidates, describe the tools and assays that can be used to study bone homeostasis, and affordable (semi-)high-throughput compound testing.

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.

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.

Different mutational function of low- and high-linear energy transfer heavy-ion irradiation demonstrated by whole-genome resequencing of Arabidopsis mutants

Heavy-ion irradiation is a powerful mutagen that possesses high linear energy transfer (LET). Several studies have indicated that the value of LET affects DNA lesion formation in several ways, including the efficiency and the density of double-stranded break induction along the particle path. We assumed that the mutation type can be altered by selecting an appropriate LET value. Here, we quantitatively demonstrate differences in the mutation type induced by irradiation with two representative ions, Ar ions (LET: 290 keV μm^-1 ) and C ions (LET: 30.0 keV μm^-1 ), by whole-genome resequencing of the Arabidopsis mutants produced by these irradiations. Ar ions caused chromosomal rearrangements or large deletions (≥100 bp) more frequently than C ions, with 10.2 and 2.3 per mutant genome under Ar- and C-ion irradiation, respectively. Conversely, C ions induced more single-base substitutions and small indels (<100 bp) than Ar ions, with 28.1 and 56.9 per mutant genome under Ar- and C-ion irradiation, respectively. Moreover, the rearrangements induced by Ar-ion irradiation were more complex than those induced by C-ion irradiation, and tended to accompany single base substitutions or small indels located close by. In conjunction with the detection of causative genes through high-throughput sequencing, selective irradiation by beams with different effects will be a powerful tool for forward genetics as well as studies on chromosomal rearrangements.

Whole Genome Sequencing-Based Mapping and Candidate Identification of Mutations from Fixed Zebrafish Tissue

As forward genetic screens in zebrafish become more common, the number of mutants that cannot be identified by gross morphology or through transgenic approaches, such as many nervous system defects, has also increased. Screening for these difficult-to-visualize phenotypes demands techniques such as whole-mount in situ hybridization (WISH) or antibody staining, which require tissue fixation. To date, fixed tissue has not been amenable for generating libraries for whole genome sequencing (WGS). Here, we describe a method for using genomic DNA from fixed tissue and a bioinformatics suite for WGS-based mapping of zebrafish mutants. We tested our protocol using two known zebrafish mutant alleles, gpr126^st49 and egr2b^fh227, both of which cause myelin defects. As further proof of concept we mapped a novel mutation, stl64, identified in a zebrafish WISH screen for myelination defects. We linked stl64 to chromosome 1 and identified a candidate nonsense mutation in the F-box and WD repeat domain containing 7 (fbxw7) gene. Importantly, stl64 mutants phenocopy previously described fbxw7^vu56 mutants, and knockdown of fbxw7 in wild-type animals produced similar defects, demonstrating that stl64 disrupts fbxw7. Together, these data show that our mapping protocol can map and identify causative lesions in mutant screens that require tissue fixation for phenotypic analysis.

Sequence-Based Mapping and Genome Editing Reveal Mutations in Stickleback Hps5 Cause Oculocutaneous Albinism and the casper Phenotype

Here, we present and characterize the spontaneous X-linked recessive mutation casper, which causes oculocutaneous albinism in threespine sticklebacks (Gasterosteus aculeatus). In humans, Hermansky-Pudlak syndrome results in pigmentation defects due to disrupted formation of the melanin-containing lysosomal-related organelle (LRO), the melanosome. casper mutants display not only reduced pigmentation of melanosomes in melanophores, but also reductions in the iridescent silver color from iridophores, while the yellow pigmentation from xanthophores appears unaffected. We mapped casper using high-throughput sequencing of genomic DNA from bulked casper mutants to a region of the stickleback X chromosome (chromosome 19) near the stickleback ortholog of Hermansky-Pudlak syndrome 5 (Hps5). casper mutants have an insertion of a single nucleotide in the sixth exon of Hps5, predicted to generate an early frameshift. Genome editing using CRISPR/Cas9 induced lesions in Hps5 and phenocopied the casper mutation. Injecting single or paired Hps5 guide RNAs revealed higher incidences of genomic deletions from paired guide RNAs compared to single gRNAs. Stickleback Hps5 provides a genetic system where a hemizygous locus in XY males and a diploid locus in XX females can be used to generate an easily scored visible phenotype, facilitating quantitative studies of different genome editing approaches. Lastly, we show the ability to better visualize patterns of fluorescent transgenic reporters in Hps5 mutant fish. Thus, Hps5 mutations present an opportunity to study pigmented LROs in the emerging stickleback model system, as well as a tool to aid in assaying genome editing and visualizing enhancer activity in transgenic fish.

Genetic Screen for Postembryonic Development in the Zebrafish (Danio rerio): Dominant Mutations Affecting Adult Form

Large-scale forward genetic screens have been instrumental for identifying genes that regulate development, homeostasis, and regeneration, as well as the mechanisms of disease. The zebrafish, Danio rerio, is an established genetic and developmental model used in genetic screens to uncover genes necessary for early development. However, the regulation of postembryonic development has received less attention as these screens are more labor intensive and require extensive resources. The lack of systematic interrogation of late development leaves large aspects of the genetic regulation of adult form and physiology unresolved. To understand the genetic control of postembryonic development, we performed a dominant screen for phenotypes affecting the adult zebrafish. In our screen, we identified 72 adult viable mutants showing changes in the shape of the skeleton as well as defects in pigmentation. For efficient mapping of these mutants and mutation identification, we devised a new mapping strategy based on identification of mutant-specific haplotypes. Using this method in combination with a candidate gene approach, we were able to identify linked mutations for 22 out of 25 mutants analyzed. Broadly, our mutational analysis suggests that there are key genes and pathways associated with late development. Many of these pathways are shared with humans and are affected in various disease conditions, suggesting constraint in the genetic pathways that can lead to change in adult form. Taken together, these results show that dominant screens are a feasible and productive means to identify mutations that can further our understanding of gene function during postembryonic development and in disease.

Utility of quantitative micro-computed tomographic analysis in zebrafish to define gene function during skeletogenesis

The zebrafish is a powerful experimental model to investigate the genetic and morphologic basis of vertebrate development. Analysis of skeletogenesis in this fish is challenging as a result of the small size of the developing and adult zebrafish. Many of the bones of small fishes such as the zebrafish and medaka are quite thin, precluding many standard assays of bone quality and morphometrics commonly used on bones of larger animals. Microcomputed tomography (microCT) is a common imaging technique used for detailed analysis of the skeleton of the zebrafish and determination of mutant phenotypes. However, the utility of this modality for analysis of the zebrafish skeleton, and the effect of inherent variation among individual zebrafish, including variables such as sex, age and strain, is not well understood. Given the increased use and accessibility of microCT, we set out to define the sensitivity of microCT methods in developing and adult zebrafish. We assessed skeletal shape and density measures in the developing vertebrae and parasphenoid of the skull base. We found most skeletal variables are tightly correlated to standard length, but that at later growth stages (>3months) there are age dependent effects on some skeletal measures. Further we find modest strain but not sex differences in skeletal measures. These data suggest that the appropriate control for assessing mutant phenotypes should be age and strain matched, ideally a wild-type sibling. By analyzing two mutants exhibiting skeletal dysplasia, we show that microCT imaging can be a sensitive method to quantify distinct skeletal parameters of adults. Finally, as developing zebrafish skeletons remain difficult to resolve by radiographic means, we define a contrast agent specific for bone that enhances resolution at early stages, permitting detailed morphometric analysis of the forming skeleton. This increased capability for detection extends the use of this imaging modality to leverage the zebrafish model to understand the development causes of skeletal dysplasias.

LITTLE FISH, BIG DATA: ZEBRAFISH AS A MODEL FOR CARDIOVASCULAR AND METABOLIC DISEASE

The burden of cardiovascular and metabolic diseases worldwide is staggering. The emergence of systems approaches in biology promises new therapies, faster and cheaper diagnostics, and personalized medicine. However, a profound understanding of pathogenic mechanisms at the cellular and molecular levels remains a fundamental requirement for discovery and therapeutics. Animal models of human disease are cornerstones of drug discovery as they allow identification of novel pharmacological targets by linking gene function with pathogenesis. The zebrafish model has been used for decades to study development and pathophysiology. More than ever, the specific strengths of the zebrafish model make it a prime partner in an age of discovery transformed by big-data approaches to genomics and disease. Zebrafish share a largely conserved physiology and anatomy with mammals. They allow a wide range of genetic manipulations, including the latest genome engineering approaches. They can be bred and studied with remarkable speed, enabling a range of large-scale phenotypic screens. Finally, zebrafish demonstrate an impressive regenerative capacity scientists hope to unlock in humans. Here, we provide a comprehensive guide on applications of zebrafish to investigate cardiovascular and metabolic diseases. We delineate advantages and limitations of zebrafish models of human disease and summarize their most significant contributions to understanding disease progression to date.

AMAP: A pipeline for whole-genome mutation detection in Arabidopsis thaliana

Detection of mutations at the whole-genome level is now possible by the use of high-throughput sequencing. However, determining mutations is a time-consuming process due to the number of false positives provided by mutation-detecting programs. AMAP (automated mutation analysis pipeline) was developed to overcome this issue. AMAP integrates a set of well-validated programs for mapping (BWA), removal of potential PCR duplicates (Picard), realignment (GATK) and detection of mutations (SAMtools, GATK, Pindel, BreakDancer and CNVnator). Thus, all types of mutations such as base substitution, deletion, insertion, translocation and chromosomal rearrangement can be detected by AMAP. In addition, AMAP automatically distinguishes false positives by comparing lists of candidate mutations in sequenced mutants. We tested AMAP by inputting already analyzed read data derived from three individual Arabidopsis thaliana mutants and confirmed that all true mutations were included in the list of candidate mutations. The result showed that the number of false positives was reduced to 12% of that obtained in a previous analysis that lacked a process of reducing false positives. Thus, AMAP will accelerate not only the analysis of mutation induction by individual mutagens but also the process of forward genetics.

Genetic screens for mutations affecting adult traits and parental-effect genes

Forward genetics remains an important approach for the unbiased identification of factors involved in biological pathways. Forward genetic analysis in the zebrafish has until now largely been restricted to the developmental period from zygotic genome activation through the end of embryogenesis. However, the use of the zebrafish as a model system for the analysis of late larval, juvenile and adult traits, including fertility and maternal and paternal effects, continues to gain momentum. Here, we describe two approaches, based on an F3-extended family and gynogenetic methods, that allow genetic screening for, and recovery of mutations affecting post-embryonic stages, including adult traits, fertility, and parental effects. For each approach, we also describe strategies to maintain, map, and molecularly clone the identified mutations.

Mapping Challenging Mutations by Whole-Genome Sequencing

Whole-genome sequencing provides a rapid and powerful method for identifying mutations on a global scale, and has spurred a renewed enthusiasm for classical genetic screens in model organisms. The most commonly characterized category of mutation consists of monogenic, recessive traits, due to their genetic tractability. Therefore, most of the mapping methods for mutation identification by whole-genome sequencing are directed toward alleles that fulfill those criteria (i.e., single-gene, homozygous variants). However, such approaches are not entirely suitable for the characterization of a variety of more challenging mutations, such as dominant and semidominant alleles or multigenic traits. Therefore, we have developed strategies for the identification of those classes of mutations, using polymorphism mapping in Caenorhabditis elegans as our model for validation. We also report an alternative approach for mutation identification from traditional recombinant crosses, and a solution to the technical challenge of sequencing sterile or terminally arrested strains where population size is limiting. The methods described herein extend the applicability of whole-genome sequencing to a broader spectrum of mutations, including classes that are difficult to map by traditional means.

Oligodendrocyte differentiation

In the nervous system, axons transmit information in the form of electrical impulses over long distances. The speed of impulse conduction is enhanced by myelin, a lipid-rich membrane that wraps around axons. Myelin also is required for the long-term health of axons by providing metabolic support. Accordingly, myelin deficiencies are implicated in a wide range of neurodevelopmental and neuropsychiatric disorders, intellectual disabilities, and neurodegenerative conditions. Central nervous system myelin is formed by glial cells called oligodendrocytes. During development, oligodendrocyte precursor cells migrate from their origins to their target axons, extend long membrane processes that wrap axons, and produce the proteins and lipids that provide myelin membrane with its unique characteristics. Myelination is a dynamic process that involves intricate interactions between multiple cell types. Therefore, an in vivo myelination model, such as the zebrafish, which allows for live observation of cell dynamics and cell-to-cell interactions, is well suited for investigating oligodendrocyte development. Zebrafish offer several advantages to investigating myelination, including the use of transgenic reporter lines, live imaging, forward genetic screens, chemical screens, and reverse genetic approaches. This chapter will describe how these tools and approaches have provided new insights into the regulatory mechanisms that guide myelination.

Insights into mechanisms of central nervous system myelination using zebrafish

Myelin is the multi-layered membrane that surrounds most axons and is produced by oligodendrocytes in the central nervous system (CNS). In addition to its important role in enabling rapid nerve conduction, it has become clear in recent years that myelin plays additional vital roles in CNS function. Myelinating oligodendrocytes provide metabolic support to axons and active myelination is even involved in regulating forms of learning and memory formation. However, there are still large gaps in our understanding of how myelination by oligodendrocytes is regulated. The small tropical zebrafish has become an increasingly popular model organism to investigate many aspects of nervous system formation, function, and regeneration. This is mainly due to two approaches for which the zebrafish is an ideally suited vertebrate model--(1) in vivo live cell imaging using vital dyes and genetically encoded reporters, and (2) gene and target discovery using unbiased screens. This review summarizes how the use of zebrafish has helped understand mechanisms of oligodendrocyte behavior and myelination in vivo and discusses the potential use of zebrafish to shed light on important future questions relating to myelination in the context of CNS development, function and repair.

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.

Genetic mapping of the nulliplex-branch gene (gb_nb1) in cotton using next-generation sequencing

Using bulked segregant analysis based on next-generation sequencing, the recessive nulliplex-branch gene was mapped between two SNP markers ~600 kb apart. In a "nulliplex-branch" cotton mutant, most of the flowers arise directly from leaf axils on the main shoot, which usually does not have a fruiting branch. A nulliplex-branch is a useful trait by which to study cotton architecture; however, the genetic basis of this mutant has remained elusive. In this study, bulked segregant analysis combined with next-generation sequencing technology was used to finely map the underlying genes that result in a nulliplex-branch plant. The nulliplex-branch Pima cotton variety, Xinhai-18, was crossed with the normal branch upland cotton line, TM-1, resulting in an F2 population. The nulliplex-branch trait was found to be controlled by the recessive gene gb_nb1. Allelic single-nucleotide polymorphisms (SNPs) were discovered by reduced-representation sequencing between the parents, and their profiles were also characterized in the nulliplex-branch and normal branch bulks constructed using the F2 plants. A candidate ~9.0 Mb-long region comprising 42 SNP markers was found to be associated with gb_nb1, which helped localize it at the ~600-kb interval on Chr 16 by segregation analysis in the F2 population. The closely linked markers with gb_nb1 developed in this study will facilitate the marker-assisted selection of the nulliplex-branch trait, and the fine map constructed will accelerate map-based cloning of gb_nb1.

Localization of causal locus in the genome of the brown macroalga Ectocarpus: NGS-based mapping and positional cloning approaches

Mutagenesis is the only process by which unpredicted biological gene function can be identified. Despite that several macroalgal developmental mutants have been generated, their causal mutation was never identified, because experimental conditions were not gathered at that time. Today, progresses in macroalgal genomics and judicious choices of suitable genetic models make mutated gene identification possible. This article presents a comparative study of two methods aiming at identifying a genetic locus in the brown alga Ectocarpus siliculosus: positional cloning and Next-Generation Sequencing (NGS)-based mapping. Once necessary preliminary experimental tools were gathered, we tested both analyses on an Ectocarpus morphogenetic mutant. We show how a narrower localization results from the combination of the two methods. Advantages and drawbacks of these two approaches as well as potential transfer to other macroalgae are discussed.

Fast forward genetics to identify mutations causing a high light tolerant phenotype in Chlamydomonas reinhardtii by whole-genome-sequencing

BACKGROUND

High light tolerance of microalgae is a desired phenotype for efficient cultivation in large scale production systems under fluctuating outdoor conditions. Outdoor cultivation requires the use of either wild-type or non-GMO derived mutant strains due to safety concerns. The identification and molecular characterization of such mutants derived from untagged forward genetics approaches was limited previously by the tedious and time-consuming methods involving techniques such as classical meiotic mapping. The combination of mapping with next generation sequencing technologies offers alternative strategies to identify genes involved in high light adaptation in untagged mutants.

RESULTS

We used the model alga Chlamydomonas reinhardtii in a non-GMO mutation strategy without any preceding crossing step or pooled progeny to identify genes involved in the regulatory processes of high light adaptation. To generate high light tolerant mutants, wildtype cells were mutagenized only to a low extent, followed by a stringent selection. We performed whole-genome sequencing of two independent mutants hit1 and hit2 and the parental wildtype. The availability of a reference genome sequence and the removal of shared bakground variants between the wildtype strain and each mutant, enabled us to identify two single nucleotide polymorphisms within the same gene Cre02.g085050, hereafter called LRS1 (putative Light Response Signaling protein 1). These two independent single amino acid exchanges are both located in the putative WD40 propeller domain of the corresponding protein LRS1. Both mutants exhibited an increased rate of non-photochemical-quenching (NPQ) and an improved resistance against chemically induced reactive oxygen species. In silico analyses revealed homology of LRS1 to the photoregulatory protein COP1 in plants.

CONCLUSIONS

In this work we identified the nuclear encoded gene LRS1 as an essential factor for high light adaptation in C. reinhardtii. The causative random mutation within this gene was identified by a rapid and efficient method, avoiding any preceding crossing step, meiotic mapping, or pooled progeny. Our results open up new insights into mechanisms of high light adaptation in microalgae and at the same time provide a simplified strategy for non-GMO forward genetics, a crucial precondition that could result in the identification of key factors for economically relevant biological processes within algae.

Kinesin family member 6 (kif6) is necessary for spine development in zebrafish

BACKGROUND

Idiopathic scoliosis is a form of spinal deformity that affects 2-3% of children and results in curvature of the spine without structural defects of the vertebral units. The pathogenesis of idiopathic scoliosis remains poorly understood, in part due to the lack of a relevant animal model.

RESULTS

We performed a forward mutagenesis screen in zebrafish to identify new models for idiopathic scoliosis. We isolated a recessive zebrafish mutant, called skolios, which develops isolated spinal curvature that arises independent of vertebral malformations. Using meiotic mapping and whole genome sequencing, we identified a nonsense mutation in kinesin family member 6 (kif6(gw326) ) unique to skolios mutants. Three additional kif6 frameshift alleles (gw327, gw328, gw329) were generated with transcription activator-like effector nucleases (TALENs). Zebrafish homozygous or compound heterozygous for kif6 frameshift mutations developed a scoliosis phenotype indistinguishable from skolios mutants, confirming that skolios is caused by the loss of kif6. Although kif6 may play a role in cilia, no evidence for cilia dysfunction was seen in kif6(gw326) mutants.

CONCLUSIONS

Overall, these findings demonstrate a novel role for kif6 in spinal development and identify a new candidate gene for human idiopathic scoliosis.

Using next-generation sequencing to isolate mutant genes from forward genetic screens

The long-lasting success of forward genetic screens relies on the simple molecular basis of the characterized phenotypes, which are typically caused by mutations in single genes. Mapping the location of causal mutations using genetic crosses has traditionally been a complex, multistep procedure, but next-generation sequencing now allows the rapid identification of causal mutations at single-nucleotide resolution even in complex genetic backgrounds. Recent advances of this mapping-by-sequencing approach include methods that are independent of reference genome sequences, genetic crosses and any kind of linkage information, which make forward genetics amenable for species that have not been considered for forward genetic screens so far.

Fish is Fish: the use of experimental model species to reveal causes of skeletal diversity in evolution and disease

Fishes are wonderfully diverse. This variety is a result of the ability of ray-finned fishes to adapt to a wide range of environments, and has made them more specious than the rest of vertebrates combined. With such diversity it is easy to dismiss comparisons between distantly related fishes in efforts to understand the biology of a particular fish species. However, shared ancestry and the conservation of developmental mechanisms, morphological features and physiology provide the ability to use comparative analyses between different organisms to understand mechanisms of development and physiology. The use of species that are amenable to experimental investigation provides tools to approach questions that would not be feasible in other 'non-model' organisms. For example, the use of small teleost fishes such as zebrafish and medaka has been powerful for analysis of gene function and mechanisms of disease in humans, including skeletal diseases. However, use of these fish to aid in understanding variation and disease in other fishes has been largely unexplored. This is especially evident in aquaculture research. Here we highlight the utility of these small laboratory fishes to study genetic and developmental factors that underlie skeletal malformations that occur under farming conditions. We highlight several areas in which model species can serve as a resource for identifying the causes of variation in economically important fish species as well as to assess strategies to alleviate the expression of the variant phenotypes in farmed fish. We focus on genetic causes of skeletal deformities in the zebrafish and medaka that closely resemble phenotypes observed both in farmed as well as natural populations of fishes.

Production of haploid zebrafish embryos by in vitro fertilization

The zebrafish has become a mainstream vertebrate model that is relevant for many disciplines of scientific study. Zebrafish are especially well suited for forward genetic analysis of developmental processes due to their external fertilization, embryonic size, rapid ontogeny, and optical clarity – a constellation of traits that enable the direct observation of events ranging from gastrulation to organogenesis with a basic stereomicroscope. Further, zebrafish embryos can survive for several days in the haploid state. The production of haploid embryos in vitro is a powerful tool for mutational analysis, as it enables the identification of recessive mutant alleles present in first generation (F1) female carriers following mutagenesis in the parental (P) generation. This approach eliminates the necessity to raise multiple generations (F2, F3, etc.) which involves breeding of mutant families, thus saving the researcher time along with reducing the needs for zebrafish colony space, labor, and the husbandry costs. Although zebrafish have been used to conduct forward screens for the past several decades, there has been a steady expansion of transgenic and genome editing tools. These tools now offer a plethora of ways to create nuanced assays for next generation screens that can be used to further dissect the gene regulatory networks that drive vertebrate ontogeny. Here, we describe how to prepare haploid zebrafish embryos. This protocol can be implemented for novel future haploid screens, such as in enhancer and suppressor screens, to address the mechanisms of development for a broad number of processes and tissues that form during early embryonic stages.

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.

Using zebrafish to study podocyte genesis during kidney development and regeneration

During development, vertebrates form a progression of up to three different kidneys that are comprised of functional units termed nephrons. Nephron composition is highly conserved across species, and an increasing appreciation of the similarities between zebrafish and mammalian nephron cell types has positioned the zebrafish as a relevant genetic system for nephrogenesis studies. A key component of the nephron blood filter is a specialized epithelial cell known as the podocyte. Podocyte research is of the utmost importance as a vast majority of renal diseases initiate with the dysfunction or loss of podocytes, resulting in a condition known as proteinuria that causes nephron degeneration and eventually leads to kidney failure. Understanding how podocytes develop during organogenesis may elucidate new ways to promote nephron health by stimulating podocyte replacement in kidney disease patients. In this review, we discuss how the zebrafish model can be used to study kidney development, and how zebrafish research has provided new insights into podocyte lineage specification and differentiation. Further, we discuss the recent discovery of podocyte regeneration in adult zebrafish, and explore how continued basic research using zebrafish can provide important knowledge about podocyte genesis in embryonic and adult environments. genesis 52:771-792, 2014. © 2014 Wiley Periodicals, Inc.

Rapid identification of kidney cyst mutations by whole exome sequencing in zebrafish

Forward genetic approaches in zebrafish have provided invaluable information about developmental processes. However, the relative difficulty of mapping and isolating mutations has limited the number of new genetic screens. Recent improvements in the annotation of the zebrafish genome coupled to a reduction in sequencing costs prompted the development of whole genome and RNA sequencing approaches for gene discovery. Here we describe a whole exome sequencing (WES) approach that allows rapid and cost-effective identification of mutations. We used our WES methodology to isolate four mutations that cause kidney cysts; we identified novel alleles in two ciliary genes as well as two novel mutants. The WES approach described here does not require specialized infrastructure or training and is therefore widely accessible. This methodology should thus help facilitate genetic screens and expedite the identification of mutants that can inform basic biological processes and the causality of genetic disorders in humans.