Coupled mutagenesis screens and genetic mapping in zebrafish

Forward genetic analysis is one of the principal advantages of the zebrafish model system. However, managing zebrafish mutant lines derived from mutagenesis screens and mapping the corresponding mutations and integrating them into the larger collection of mutations remain arduous tasks. To simplify and focus these endeavors, we developed an approach that facilitates the rapid mapping of new zebrafish mutations as they are generated through mutagenesis screens. We selected a minimal panel of 149 simple sequence length polymorphism markers for a first-pass genome scan in crosses involving C32 and SJD inbred lines. We also conducted a small chemical mutagenesis screen that identified several new mutations affecting zebrafish embryonic melanocyte development. Using our first-pass marker panel in bulked-segregant analysis, we were able to identify the genetic map positions of these mutations as they were isolated in our screen. Rapid mapping of the mutations facilitated stock management, helped direct allelism tests, and should accelerate identification of the affected genes. These results demonstrate the efficacy of coupling mutagenesis screens with genetic mapping.     

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

The generation and analysis of mutants in zebrafish has been instrumental in defining the genetic regulation of vertebrate development, physiology, and disease. However, identifying the genetic changes that underlie mutant phenotypes remains a significant bottleneck in the analysis of mutants. Whole-genome sequencing has recently emerged as a fast and efficient approach for identifying mutations in nonvertebrate model organisms. However, this approach has not been applied to zebrafish due to the complicating factors of having a large genome and lack of fully inbred lines. Here we provide a method for efficiently mapping and detecting mutations in zebrafish using these new parallel sequencing technologies. This method utilizes an extensive reference SNP database to define regions of homozygosity-by-descent by low coverage, whole-genome sequencing of pooled DNA from only a limited number of mutant F(2) fish. With this approach we mapped each of the five different zebrafish mutants we sequenced and identified likely causative nonsense mutations in two and candidate mutations in the remainder. Furthermore, we provide evidence that one of the identified mutations, a nonsense mutation in bmp1a, underlies the welded mutant phenotype.     

High-density detection of restriction-site-associated DNA markers for rapid mapping of mutated loci in Neurospora

The wealth of sequence information available for Neurospora crassa and other fungi has greatly facilitated evolutionary and molecular analyses of this group. Although "reverse" genetics, in which genes are first identified by their sequence rather than by their mutant phenotypes, serves as a valuable new approach for elucidating biological processes, classical "forward" genetic analysis is still extremely useful. Unfortunately, mapping mutations and identifying the corresponding genes has typically been slow and laborious. To facilitate forward genetics in Neurospora, we have adapted microarray-based restriction-site-associated DNA (RAD) mapping for use with N. crassa oligonucleotide microarrays. This technique was used to simultaneously detect an unprecedented number of genomewide restriction site polymorphisms from two N. crassa strains: Mauriceville and Oak Ridge. Furthermore, RAD mapping was used to quickly map a previously unknown gene, defective in methylation-7 (dim-7).     

Rapid SNP discovery and genetic mapping using sequenced RAD markers

Single nucleotide polymorphism (SNP) discovery and genotyping are essential to genetic mapping. There remains a need for a simple, inexpensive platform that allows high-density SNP discovery and genotyping in large populations. Here we describe the sequencing of restriction-site associated DNA (RAD) tags, which identified more than 13,000 SNPs, and mapped three traits in two model organisms, using less than half the capacity of one Illumina sequencing run. We demonstrated that different marker densities can be attained by choice of restriction enzyme. Furthermore, we developed a barcoding system for sample multiplexing and fine mapped the genetic basis of lateral plate armor loss in threespine stickleback by identifying recombinant breakpoints in F(2) individuals. Barcoding also facilitated mapping of a second trait, a reduction of pelvic structure, by in silico re-sorting of individuals. To further demonstrate the ease of the RAD sequencing approach we identified polymorphic markers and mapped an induced mutation in Neurospora crassa. Sequencing of RAD markers is an integrated platform for SNP discovery and genotyping. This approach should be widely applicable to genetic mapping in a variety of organisms.     

Rapid mapping and identification of mutations in Caenorhabditis elegans by restriction site-associated DNA mapping and genomic interval pull-down sequencing

Forward genetic screens provide a powerful approach for inferring gene function on the basis of the phenotypes associated with mutated genes. However, determining the causal mutation by traditional mapping and candidate gene sequencing is often the rate-limiting step, especially when analyzing many mutants. We report two genomic approaches for more rapidly determining the identity of the affected genes in Caenorhabditis elegans mutants. First, we report our use of restriction site-associated DNA (RAD) polymorphism markers for rapidly mapping mutations after chemical mutagenesis and mutant isolation. Second, we describe our use of genomic interval pull-down sequencing (GIPS) to selectively capture and sequence megabase-sized portions of a mutant genome. Together, these two methods provide a rapid and cost-effective approach for positional cloning of C. elegans mutant loci, and are also applicable to other genetic model systems.     

Centromere-linked microsatellite markers for linkage groups 3, 4, 6, 7, 13, and 20 of zebrafish (Danio rerio)

A large number of interesting mutations affecting development and organogenesis have been identified through genetic screens in zebrafish. Mapping of these mutations to a chromosomal region can be rapidly accomplished using half-tetrad analysis. However, knowledge of centromere-linked markers on every chromosome is essential to this mapping method. Centromeres on all 25 linkage groups have been mapped on the RAPD zebrafish genetic map. However, species specificity and the lack of codominance make RAPD markers less practical for mapping than microsatellite-based markers. On the microsatellite-based genetic map, centromere-linked markers have been identified for 19 linkage groups. No direct evidence has been published linking microsatellite markers to the centromeres of linkage groups 3, 4, 6, 7, 13, and 20. Therefore, we compared the microsatellite-based genetic map with the RAPD map to identify markers most likely linked to the centromeres of these 6 linkage groups. These candidate markers were tested for potential centromere linkage using four panels of half-tetrad embryos derived by early-pressure treatment of eggs from four different female zebrafish. We have identified microsatellite markers for linkage groups 3, 4, 6, 7, 13, and 20 to within 1.7 cM of their centromeres. These markers will greatly facilitate the rapid mapping of mutations in zebrafish by half-tetrad analysis.     

Forward genetic analysis of visual behavior in zebrafish

The visual system converts the distribution and wavelengths of photons entering the eye into patterns of neuronal activity, which then drive motor and endocrine behavioral responses. The gene products important for visual processing by a living and behaving vertebrate animal have not been identified in an unbiased fashion. Likewise, the genes that affect development of the nervous system to shape visual function later in life are largely unknown. Here we have set out to close this gap in our understanding by using a forward genetic approach in zebrafish. Moving stimuli evoke two innate reflexes in zebrafish larvae, the optomotor and the optokinetic response, providing two rapid and quantitative tests to assess visual function in wild-type (WT) and mutant animals. These behavioral assays were used in a high-throughput screen, encompassing over half a million fish. In almost 2,000 F2 families mutagenized with ethylnitrosourea, we discovered 53 recessive mutations in 41 genes. These new mutations have generated a broad spectrum of phenotypes, which vary in specificity and severity, but can be placed into only a handful of classes. Developmental phenotypes include complete absence or abnormal morphogenesis of photoreceptors, and deficits in ganglion cell differentiation or axon targeting. Other mutations evidently leave neuronal circuits intact, but disrupt phototransduction, light adaptation, or behavior-specific responses. Almost all of the mutants are morphologically indistinguishable from WT, and many survive to adulthood. Genetic linkage mapping and initial molecular analyses show that our approach was effective in identifying genes with functions specific to the visual system. This collection of zebrafish behavioral mutants provides a novel resource for the study of normal vision and its genetic disorders.     

Spontaneous Mitotic Recombination in Yeast: The Hyper-Recombinational Rem1 Mutations Are Alleles of the Rad3 Gene

The RAD3 gene of Saccharomyces cerevisiae is required for UV excision-repair and is essential for cell viability. We have identified the rem1 mutations (enhanced spontaneous mitotic recombination and mutation) of Saccharomyces cerevisiae as alleles of RAD3 by genetic mapping, complementation with the cloned wild-type gene, and DNA hybridization. The high levels of spontaneous mitotic gene conversion, crossing over, and mutation conferred upon cells by the rem1 mutations are distinct from the effects of all other alleles of RAD3. We present preliminary data on the localization of the rem1 mutations within the RAD3 gene. The interaction of the rem1 mutant alleles with a number of radiation-sensitive mutations is also different than the interactions reported for previously described (UV-sensitive) alleles of RAD3. Double mutants of rem1 and a defect in the recombination-repair pathway are inviable, while double mutants containing UV-sensitive alleles of RAD3 are viable. The data presented here demonstrate that: (1) rem1 strains containing additional mutations in other excision-repair genes do not exhibit elevated gene conversion; (2) triple mutants containing rem1 and mutations in both excision-repair and recombination-repair are viable; (3) such triple mutants containing rad52 have reduced levels of gene conversion but wild-type frequencies of crossing over. We have interpreted these observations in a model to explain the effects of rem1. Consistent with the predictions of the model, we find that the size of DNA from rem1 strains, as measured by neutral sucrose gradients, is smaller than wild type.     

Genetic analysis of chromosomal rearrangements in the cyclops region of the zebrafish genome.

Genetic screens in zebrafish have provided mutations in hundreds of genes with essential functions in the developing embryo. To investigate the possible uses of chromosomal rearrangements in the analysis of these mutations, we genetically characterized three gamma-ray induced alleles of cyclops (cyc), a gene required for development of midline structures. We show that cyc maps near one end of Linkage Group 12 (LG 12) and that this region is involved in a reciprocal translocation with LG 2 in one gamma-ray induced mutation, cyc(b213). The translocated segments together cover approximately 5% of the genetic map, and we show that this rearrangement is useful for mapping cloned genes that reside in the affected chromosomal regions. The other two alleles, cyc(b16) and cyc(b229), have deletions in the distal region of LG 12. Interestingly, both of these mutations suppress recombination between genetic markers in LG 12, including markers at a distance from the deletion. This observation raises the possibility that these deletions affect a site required for meiotic recombination on the LG 12 chromosome. The cyc(b16) and cyc(b229) mutations may be useful for balancing other lethal mutations located in the distal region of LG 12. These results show that chromosomal rearrangements can provide useful resources for mapping and genetic analyses in zebrafish.     

Linkage mapping and comparative genomics using next-generation RAD sequencing of a non-model organism

Restriction-site associated DNA (RAD) sequencing is a powerful new method for targeted sequencing across the genomes of many individuals. This approach has broad potential for genetic analysis of non-model organisms including genotype-phenotype association mapping, phylogeography, population genetics and scaffolding genome assemblies through linkage mapping. We constructed a RAD library using genomic DNA from a Plutella xylostella (diamondback moth) backcross that segregated for resistance to the insecticide spinosad. Sequencing of 24 individuals was performed on a single Illumina GAIIx lane (51 base paired-end reads). Taking advantage of the lack of crossing over in homologous chromosomes in female Lepidoptera, 3,177 maternally inherited RAD alleles were assigned to the 31 chromosomes, enabling identification of the spinosad resistance and W/Z sex chromosomes. Paired-end reads for each RAD allele were assembled into contigs and compared to the genome of Bombyx mori (n = 28) using BLAST, revealing 28 homologous matches plus 3 expected fusion/breakage events which account for the difference in chromosome number. A genome-wide linkage map (1292 cM) was inferred with 2,878 segregating RAD alleles inherited from the backcross father, producing chromosome and location specific sequenced RAD markers. Here we have used RAD sequencing to construct a genetic linkage map de novo for an organism that has no previous genome data. Comparative analysis of P. xyloxtella linkage groups with B. mori chromosomes shows for the first time, genetic synteny appears common beyond the Macrolepidoptera. RAD sequencing is a powerful system capable of rapidly generating chromosome specific data for non-model organisms.     

Screening mosaic F1 females for mutations affecting zebrafish heart induction and patterning

The genetic pathways underlying the induction and anterior-posterior patterning of the heart are poorly understood. The recent emergence of the zebrafish model system now allows a classical genetic approach to such challenging problems in vertebrate development. Two large-scale screens for mutations affecting zebrafish embryonic development have recently been completed; among the hundreds of mutations identified were several that affect specific aspects of cardiac morphogenesis, differentiation, and function. However, very few mutations affecting induction and/or anterior-posterior patterning of the heart were identified. We hypothesize that a directed approach utilizing molecular markers to examine these particular steps of heart development will uncover additional such mutations. To test this hypothesis, we are conducting two parallel screens for mutations that affect either the induction or the anterior-posterior patterning of the zebrafish heart. As an indicator of cardiac induction, we examine expression of nkx2.5, the earliest known marker of precardiac mesoderm; to assess anterior-posterior patterning, we distinguish ventricle from atrium with antibodies that recognize different myosin heavy chain isoforms. In order to expedite the examination of a large number of mutations, we are screening the haploid progeny of mosaic F1 females. In these ongoing screens, we have identified four mutations that affect nkx2.5 expression as well as 21 that disrupt either ventricular or atrial development and thus far have recovered several of these mutations, demonstrating the value of our approach. Future analysis of these and other cardiac mutations will provide further insight into the processes of induction and anterior-posterior patterning of the heart. Dev. Genet. 22:288–299, 1998. © 1998 Wiley-Liss, Inc.     

A genetic linkage map and comparative genome analysis of common carp (Cyprinus carpio L.) using microsatellites and SNPs

A genetic linkage map is a powerful research tool for mapping traits of interest and is essential to understanding genome evolution. The aim of this study is to provide an expanded genetic linkage map of common carp to effectively carry out quantitative trait loci analysis and conduct comparative mapping analysis between lineages. Here, we constructed a genetic linkage map of common carp (Cyprinus carpio L.) using microsatellite and single-nucleotide polymorphism (SNP) markers in a 159 sibling family. A total of 246 microsatellites and 306 SNP polymorphic markers were genotyped in this family. Linkage analysis using JoinMap 4.0 organized 427 markers (186 microsatellites and 241 SNPs) to 50 linkage groups, ranging in size from 1.4 to 130.1 cM. Each group contained 2-30 markers. The linkage map covered a genetic distance of 2,039.2 cM and the average interval for markers within the linkage groups was approximately 6.4 cM. In addition, comparative genome analysis within five model teleost fish revealed a high percentage (74.7%) of conserved loci corresponding to zebrafish chromosomes. In most cases, each zebrafish chromosome comprised two common carp linkage groups. The comparative analysis also revealed independent chromosome rearrangements in common carp and zebrafish. The linkage map will be of great assistance in mapping genes of interest and serve as a reference to approach comparative mapping and enable further insights into the comprehensive investigations of genome evolution of common carp.     

A meta-analysis of kidney microarray datasets: investigation of cytokine gene detection and correlation with rt-PCR and detection thresholds

Background

Large-scale mutagenesis screens in the zebrafish employing the mutagen ENU have isolated several hundred mutant loci that represent putative developmental control genes. In order to realize the potential of such screens, systematic genetic mapping of the mutations is necessary. Here we report on a large-scale effort to map the mutations generated in mutagenesis screening at the Max Planck Institute for Developmental Biology by genome scanning with microsatellite markers.

Results

We have selected a set of microsatellite markers and developed methods and scoring criteria suitable for efficient, high-throughput genome scanning. We have used these methods to successfully obtain a rough map position for 319 mutant loci from the Tübingen I mutagenesis screen and subsequent screening of the mutant collection. For 277 of these the corresponding gene is not yet identified. Mapping was successful for 80 % of the tested loci. By comparing 21 mutation and gene positions of cloned mutations we have validated the correctness of our linkage group assignments and estimated the standard error of our map positions to be approximately 6 cM.

Conclusion

By obtaining rough map positions for over 300 zebrafish loci with developmental phenotypes, we have generated a dataset that will be useful not only for cloning of the affected genes, but also to suggest allelism of mutations with similar phenotypes that will be identified in future screens. Furthermore this work validates the usefulness of our methodology for rapid, systematic and inexpensive microsatellite mapping of zebrafish mutations.     

Patterning during organogenesis: genetic analysis of cardiac chamber formation.

A classical genetic approach, in which mutagenized organisms are screened for phenotypes of interest, is appealing for the analysis of developmental processes. Here, we describe the advantages of zebrafish genetics for the study of heart development. As an example of the utility of this strategy, we discuss its potential to illuminate the molecular mechanisms of cardiac chamber formation. The signals that specify ventricular and atrial lineages and the differentiation pathways that produce distinct chambers are poorly understood. Recently identified zebrafish mutations that disrupt ventricular or atrial development promise to reveal genes essential for these processes.     

Acridine mutagenesis of zebrafish (Danio rerio)

Mutagenesis screening, in which heritable traits are isolated following damage to the genome, is a powerful approach for investigating gene function. Among vertebrate model organisms, the zebrafish (Danio rerio) is ideally suited to mutagenesis screens. The success of large-scale screens is dependent on the way in which changes are identified. The type of damage induced is also pivotal. Single base coding region deletions and insertions are suited to abolition of gene function whilst inducing a small physical alteration to the genome. Such mutations are not commonly found following mutagenesis schemes reported to date. Here, we show that an acridine mutagen, ICR191, which in other model organisms frequently induces single base deletions and insertions, is mutagenic in zebrafish. ICR191 induces hallmark phenotypes associated with genetic damage in treated embryos. Alterations are heritable. Offspring of mutagenised fish had mutations in a marker gene and were found to produce offspring with abnormal development. Using an adaptation of a molecular mutation detection method, fluorescent arbitrary primed PCR, we identified an induced alteration directly. The estimated frequency of induced mutations was sufficiently high to make it feasible to employ this approach for mutagenesis screening.     

Interaction between checkpoint genes RAD9, RAD17, RAD24, and RAD53 involved in the determination of yeast Saccharomyces cerevisiae sensitivity to ionizing radiation

Mechanisms for genetic control of cell division cycle (checkpoint control) have been studied in most detail in yeast Saccharomyces cerevisiae. To clarify the role of checkpoint genes RAD9, RAD17, RAD24, and RAD53 in cell radioresistance, double mutants were analyzed for cell sensitivity to ionizing radiation. Double mutants carrying mutations in combination with mutation rad9delta were shown to manifest the epistatic type of interaction. Our results suggest that checkpoint genes RAD9, RAD17, RAD24, and RAD53 belong to a single epistatic group designated RAD9 and govern the same pathway. Genes RAD9 and RAD53 have a positive effect on sensitivity to gamma-radiation, whereas RAD17 and RAD24 have a negative effect. Interactions between mutations may differ when considering their sensitivity to gamma-radiation and UV light; mutations rad9delta and rad24delta were shown to manifest the additive effect in the first case and epistatic effect in the second.     

RAD1, an excision repair gene of Saccharomyces cerevisiae, is also involved in recombination.

The RAD1 gene of Saccharomyces cerevisiae is required for the incision step of excision repair of damaged DNA. In this paper, we report our observations on the effect of the RAD1 gene on genetic recombination. Mitotic intrachromosomal and interchromosomal recombination in RAD+, rad1, rad52, and other rad mutant strains was examined. The rad1 deletion mutation and some rad1 point mutations reduced the frequency of intrachromosomal recombination of a his3 duplication, in which one his3 allele is deleted at the 3' end while the other his3 allele is deleted at the 5' end. Mutations in the other excision repair genes, RAD2, RAD3, and RAD4, did not lower recombination frequencies in the his3 duplication. As expected, recombination between the his3 deletion alleles in the duplication was reduced in the rad52 mutant. The frequency of HIS3+ recombinants fell synergistically in the rad1 rad52 double mutant, indicating that the RAD1 and RAD52 genes affect this recombination via different pathways. In contrast to the effect of mutations in the RAD52 gene, mutations in the RAD1 gene did not lower intrachromosomal and interchromosomal recombination between heteroalleles that carry point mutations rather than partial deletions; however, the rad1 delta mutation did lower the frequency of integration of linear plasmids and DNA fragments into homologous genomic sequences. We suggest that RAD1 plays a role in recombination after the formation of the recombinogenic substrate.     

Interaction between checkpoint genes RAD9, RAD17, RAD24, and RAD53 involved in the determination of yeast Saccharomyces cerevisiae sensitivity to ionizing radiation

Mechanisms for genetic control of cell division cycle (checkpoint control) have been studied in most detail in yeast Saccharomyces cerevisiae. To clarify the role of checkpoint genes RAD9, RAD17, RAD24, and RAD53 in cell radioresistance, double mutants were analyzed for cell sensitivity to ionizing radiation. Double mutants carrying mutations in combination with mutation rad9Delta were shown to manifest the epistatic type of interaction. Our results suggest that checkpoint genes RAD9, RAD17, RAD24, and RAD53 belong to a single epistatic group designated RAD9 and govern the same pathway. Genes RAD9 and RAD53 have a positive effect on sensitivity to gamma-radiation, whereas RAD17 and RAD24 have a negative effect. Interactions between mutations may differ when considering their sensitivity to gamma-radiation and UV light; mutations rad9Delta and rad24Delta were shown to manifest the additive effect in the first case and epistatic effect in the second.     

Zebrafish heart failure models: opportunities and challenges

Heart failure is a complex pathophysiological syndrome of pumping functional failure that results from injury, infection or toxin-induced damage on the myocardium, as well as genetic influence. Gene mutations associated with cardiomyopathies can lead to various pathologies of heart failure. In recent years, zebrafish, Danio rerio, has emerged as an excellent model to study human cardiovascular diseases such as congenital heart defects, cardiomyopathy, and preclinical development of drugs targeting these diseases. In this review, we will first summarize zebrafish genetic models of heart failure arose from cardiomyopathy, which is caused by mutations in sarcomere, calcium or mitochondrial-associated genes. Moreover, we outline zebrafish heart failure models triggered by chemical compounds. Elucidation of these models will improve the understanding of the mechanism of pathogenesis and provide potential targets for novel therapies.