A Post-Developmental Genetic Screen for Zebrafish Models of Inherited Liver Disease

Nonalcoholic fatty liver disease (NAFLD) is one of the most common causes of chronic liver disease such as simple steatosis, nonalcoholic steatohepatitis (NASH), cirrhosis and fibrosis. However, the molecular pathogenesis and genetic variations causing NAFLD are poorly understood. The high prevalence and incidence of NAFLD suggests that genetic variations on a large number of genes might be involved in NAFLD. To identify genetic variants causing inherited liver disease, we used zebrafish as a model system for a large-scale mutant screen, and adopted a whole genome sequencing approach for rapid identification of mutated genes found in our screen. Here, we report on a forward genetic screen of ENU mutagenized zebrafish. From 250 F2 lines of ENU mutagenized zebrafish during post-developmental stages (5 to 8 days post fertilization), we identified 19 unique mutant zebrafish lines displaying visual evidence of hepatomegaly and/or steatosis with no developmental defects. Histological analysis of mutants revealed several specific phenotypes, including common steatosis, micro/macrovesicular steatosis, hepatomegaly, ballooning, and acute hepatocellular necrosis. This work has identified multiple post-developmental mutants and establishes zebrafish as a novel animal model for post-developmental inherited liver disease.     

Phenotype sequencing: identifying the genes that cause a phenotype directly from pooled sequencing of independent mutants

Random mutagenesis and phenotype screening provide a powerful method for dissecting microbial functions, but their results can be laborious to analyze experimentally. Each mutant strain may contain 50-100 random mutations, necessitating extensive functional experiments to determine which one causes the selected phenotype. To solve this problem, we propose a "Phenotype Sequencing" approach in which genes causing the phenotype can be identified directly from sequencing of multiple independent mutants. We developed a new computational analysis method showing that 1. causal genes can be identified with high probability from even a modest number of mutant genomes; 2. costs can be cut many-fold compared with a conventional genome sequencing approach via an optimized strategy of library-pooling (multiple strains per library) and tag-pooling (multiple tagged libraries per sequencing lane). We have performed extensive validation experiments on a set of E. coli mutants with increased isobutanol biofuel tolerance. We generated a range of sequencing experiments varying from 3 to 32 mutant strains, with pooling on 1 to 3 sequencing lanes. Our statistical analysis of these data (4099 mutations from 32 mutant genomes) successfully identified 3 genes (acrB, marC, acrA) that have been independently validated as causing this experimental phenotype. It must be emphasized that our approach reduces mutant sequencing costs enormously. Whereas a conventional genome sequencing experiment would have cost $7,200 in reagents alone, our Phenotype Sequencing design yielded the same information value for only $1200. In fact, our smallest experiments reliably identified acrB and marC at a cost of only $110-$340.     

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.     

Fast homozygosity mapping and identification of a zebrafish ENU-induced mutation by whole-genome sequencing

Forward genetics using zebrafish is a powerful tool for studying vertebrate development through large-scale mutagenesis. Nonetheless, the identification of the molecular lesion is still laborious and involves time-consuming genetic mapping. Here, we show that high-throughput sequencing of the whole zebrafish genome can directly locate the interval carrying the causative mutation and at the same time pinpoint the molecular lesion. The feasibility of this approach was validated by sequencing the m1045 mutant line that displays a severe hypoplasia of the exocrine pancreas. We generated 13 Gb of sequence, equivalent to an eightfold genomic coverage, from a pool of 50 mutant embryos obtained from a map-cross between the AB mutant carrier and the WIK polymorphic strain. The chromosomal region carrying the causal mutation was localized based on its unique property to display high levels of homozygosity among sequence reads as it derives exclusively from the initial AB mutated allele. We developed an algorithm identifying such a region by calculating a homozygosity score along all chromosomes. This highlighted an 8-Mb window on chromosome 5 with a score close to 1 in the m1045 mutants. The sequence analysis of all genes within this interval revealed a nonsense mutation in the snapc4 gene. Knockdown experiments confirmed the assertion that snapc4 is the gene whose mutation leads to exocrine pancreas hypoplasia. In conclusion, this study constitutes a proof-of-concept that whole-genome sequencing is a fast and effective alternative to the classical positional cloning strategies in zebrafish.     

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.     

The Mouse MC13 Mutant Is a Novel ENU Mutation in Collagen Type II,Alpha 1

Phenotype-driven mutagenesis experiments are a powerful approach to identifying novel alleles in a variety of contexts. The traditional disadvantage of this approach has been the subsequent task of identifying the affected locus in the mutants of interest. Recent advances in bioinformatics and sequencing have reduced the burden of cloning these ENU mutants. Here we report our experience with an ENU mutagenesis experiment and the rapid identification of a mutation in a previously known gene. A combination of mapping the mutation with a high-density SNP panel and a candidate gene approach has identified a mutation in collagen type II, alpha I (Col2a1). Col2a1 has previously been studied in the mouse and our mutant phenotype closely resembles mutations made in the Col2a1 locus.     

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.     

Mutational and structural analysis of the nitrate reductase heme domain of Nicotiana plumbaginifolia.

We have analyzed four Nicotiana plumbaginifolia null mutants presumably affected in the heme domain of nitrate reductase. The DNA sequence of this domain has been determined for each mutant and for the wild type. Two mutations were identified as single base changes leading to, respectively, the substitution of a histidine residue by an asparagine (mutant E56) and to the appearance of an ochre stop codon (mutant E64). Based on the amino acid sequence homology between the nitrate reductase heme domain and mammalian cytochrome b5, we have predicted the three-dimensional structure of this domain. This showed that the nitrate reductase heme domain is structurally very similar to cytochrome b5 and it also confirmed that the residue involved in E56 mutation is one of the two heme-binding histidines. The two other mutations (mutants A1 and K21) were found to be, respectively, -1 and +1 frameshift mutations resulting in the appearance of an opal stop codon. These sequence data confirmed previous genetic and biochemical hypotheses on nitrate reductase-deficient mutants. Northern blot analysis of these mutants indicated that mutant E56 overexpressed the nitrate reductase mRNA, whereas the nonsense mutations present in the other mutants led to reduced levels of nitrate reductase mRNA.     

Genome-wide studies of copy number variation and exome sequencing identify rare variants in BAG3 as a cause of dilated cardiomyopathy

Dilated cardiomyopathy commonly causes heart failure and is the most frequent precipitating cause of heart transplantation. Familial dilated cardiomyopathy has been shown to be caused by rare variant mutations in more than 30 genes but only ~35% of its genetic cause has been identified, principally by using linkage-based or candidate gene discovery approaches. In a multigenerational family with autosomal dominant transmission, we employed whole-exome sequencing in a proband and three of his affected family members, and genome-wide copy number variation in the proband and his affected father and unaffected mother. Exome sequencing identified 428 single point variants resulting in missense, nonsense, or splice site changes. Genome-wide copy number analysis identified 51 insertion deletions and 440 copy number variants > 1 kb. Of these, a 8733 bp deletion, encompassing exon 4 of the heat shock protein cochaperone BCL2-associated athanogene 3 (BAG3), was found in seven affected family members and was absent in 355 controls. To establish the relevance of variants in this protein class in genetic DCM, we sequenced the coding exons in BAG3 in 311 other unrelated DCM probands and identified one frameshift, two nonsense, and four missense rare variants absent in 355 control DNAs, four of which were familial and segregated with disease. Knockdown of bag3 in a zebrafish model recapitulated DCM and heart failure. We conclude that new comprehensive genomic approaches have identified rare variants in BAG3 as causative of DCM.     

Mutation detection with next-generation resequencing through a mediator genome

The affordability of next generation sequencing (NGS) is transforming the field of mutation analysis in bacteria. The genetic basis for phenotype alteration can be identified directly by sequencing the entire genome of the mutant and comparing it to the wild-type (WT) genome, thus identifying acquired mutations. A major limitation for this approach is the need for an a-priori sequenced reference genome for the WT organism, as the short reads of most current NGS approaches usually prohibit de-novo genome assembly. To overcome this limitation we propose a general framework that utilizes the genome of relative organisms as mediators for comparing WT and mutant bacteria. Under this framework, both mutant and WT genomes are sequenced with NGS, and the short sequencing reads are mapped to the mediator genome. Variations between the mutant and the mediator that recur in the WT are ignored, thus pinpointing the differences between the mutant and the WT. To validate this approach we sequenced the genome of Bdellovibrio bacteriovorus 109J, an obligatory bacterial predator, and its prey-independent mutant, and compared both to the mediator species Bdellovibrio bacteriovorus HD100. Although the mutant and the mediator sequences differed in more than 28,000 nucleotide positions, our approach enabled pinpointing the single causative mutation. Experimental validation in 53 additional mutants further established the implicated gene. Our approach extends the applicability of NGS-based mutant analyses beyond the domain of available reference genomes.     

Five gametophytic mutations affecting pollen development and pollen tube growth in Arabidopsis thaliana

Mutant analysis represents one of the most reliable approaches to identifying genes involved in plant development. The screening of the Versailles collection of Arabidopsis thaliana T-DNA insertion transformants has allowed us to isolate different mutations affecting male gametophytic functions and viability. Among several mutated lines, five have been extensively studied at the genetic, molecular, and cytological levels. For each mutant, several generations of selfing and outcrossing have been carried out, leading to the conclusion that all these mutations are tagged and affect only the male gametophyte. However, only one out of the five mutations is completely penetrant. A variable number of T-DNA copies has integrated in the mutant lines, although all segregate at one mutated locus. Two mutants could be defined as "early mutants": the mutated genes are presumably expressed during pollen grain maturation and their alteration leads to the production of nonfunctional pollen grains. Two other mutants could be defined as "late mutant" since their pollen is able to germinate but pollen tube growth is highly disturbed. Screening for segregation ratio distortions followed by thorough genetic analysis proved to be a powerful tool for identifying gametophytic mutations of all phases of pollen development.     

Identification of EMS-Induced Mutations in Drosophila melanogaster by Whole-Genome Sequencing

Next-generation methods for rapid whole-genome sequencing enable the identification of single-base-pair mutations in Drosophila by comparing a chromosome bearing a new mutation to the unmutagenized sequence. To validate this approach, we sought to identify the molecular lesion responsible for a recessive EMS-induced mutation affecting egg shell morphology by using Illumina next-generation sequencing. After obtaining sufficient sequence from larvae that were homozygous for either wild-type or mutant chromosomes, we obtained high-quality reads for base pairs composing ~70% of the third chromosome of both DNA samples. We verified 103 single-base-pair changes between the two chromosomes. Nine changes were nonsynonymous mutations and two were nonsense mutations. One nonsense mutation was in a gene, encore, whose mutations produce an egg shell phenotype also observed in progeny of homozygous mutant mothers. Complementation analysis revealed that the chromosome carried a new functional allele of encore, demonstrating that one round of next-generation sequencing can identify the causative lesion for a phenotype of interest. This new method of whole-genome sequencing represents great promise for mutant mapping in flies, potentially replacing conventional methods.     

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.     

Mutations in prickle orthologs cause seizures in flies, mice, and humans

Epilepsy is heritable, yet few causative gene mutations have been identified, and thus far no human epilepsy gene mutations have been found to produce seizures in invertebrates. Here we show that mutations in prickle genes are associated with seizures in humans, mice, and flies. We identified human epilepsy patients with heterozygous mutations in either PRICKLE1 or PRICKLE2. In overexpression assays in zebrafish, prickle mutations resulted in aberrant prickle function. A seizure phenotype was present in the Prickle1-null mutant mouse, two Prickle1 point mutant (missense and nonsense) mice, and a Prickle2-null mutant mouse. Drosophila with prickle mutations displayed seizures that were responsive to anti-epileptic medication, and homozygous mutant embryos showed neuronal defects. These results suggest that prickle mutations have caused seizures throughout evolution.     

Maternal control of vertebrate development before the midblastula transition: mutants from the zebrafish I

Maternal factors control development prior to the activation of the embryonic genome. In vertebrates, little is known about the molecular mechanisms by which maternal factors regulate embryonic development. To understand the processes controlled by maternal factors and identify key genes involved, we embarked on a maternal-effect mutant screen in the zebrafish. We identified 68 maternal-effect mutants. Here we describe 15 mutations in genes controlling processes prior to the midblastula transition, including egg development, blastodisc formation, embryonic polarity, initiation of cell cleavage, and cell division. These mutants exhibit phenotypes not previously observed in zygotic mutant screens. This collection of maternal-effect mutants provides the basis for a molecular genetic analysis of the maternal control of embryogenesis in vertebrates.     

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.     

Efficient mutation identification in zebrafish by microarray capturing and next generation sequencing

Fish models like medaka, stickleback or zebrafish provide a valuable resource to study vertebrate genes. However, finding genetic variants e.g. mutations in the genome is still arduous. Here we used a combination of microarray capturing and next generation sequencing to identify the affected gene in the mozartkugel^p11cv (mzl^p11cv) mutant zebrafish. We discovered a 31-bp deletion in macf1 demonstrating the potential of this technique to efficiently isolate mutations in a vertebrate genome.     

Analysis of a zebrafish VEGF receptor mutant reveals specific disruption of angiogenesis

Blood vessels form either by the assembly and differentiation of mesodermal precursor cells (vasculogenesis) or by sprouting from preexisting vessels (angiogenesis). Endothelial-specific receptor tyrosine kinases and their ligands are known to be essential for these processes. Targeted disruption of vascular endothelial growth factor (VEGF) or its receptor kdr (flk1, VEGFR2) in mouse embryos results in a severe reduction of all blood vessels, while the complete loss of flt1 (VEGFR1) leads to an increased number of hemangioblasts and a disorganized vasculature. In a large-scale forward genetic screen, we identified two allelic zebrafish mutants in which the sprouting of blood vessels is specifically disrupted without affecting the assembly and differentiation of angioblasts. Molecular cloning revealed nonsense mutations in flk1. Analysis of mRNA expression in flk1 mutant embryos showed that flk1 expression was severely downregulated, while the expression of other genes (scl, gata1, and fli1) involved in vasculogenesis or hematopoiesis was unchanged. Overexpression of vegf(121+165) led to the formation of additional vessels only in sibling larvae, not in flk1 mutants. We demonstrate that flk1 is not required for proper vasculogenesis and hematopoiesis in zebrafish embryos. However, the disruption of flk1 impairs the formation or function of vessels generated by sprouting angiogenesis.     

New tools for mitochondrial genetics of Chlamydomonas reinhardtii: manganese mutagenesis and cytoduction

Summary A novel and efficient genetic procedure is described for generating mitochondrial mutants of the green alga Chlamydomonas reinhardtii. The development of a mutagenesis procedure using manganese cations and the application of cytoduction techniques resulted in a combined approach for the generation and analysis of mitochondrial mutants. Although mitochondrial mutations are inherited in sexual crosses from the minus mating type parent, the cytoduction technique can be used to transfer mitochondrial mutations into recipient strains with different genetic backgrounds, irrespective of their mating type. Cytoduction allows the transfer of mitochondrial markers from diploid to haploid cells also, which is of great benefit since diploid cells do not germinate in C. reinhardtii. We report here the isolation and characterisation of eight mutants, which are resistant to the antibiotics myxothiazol and mucidin. The mutants all have point mutations in the mitochondrial gene for apocytochrome b. Using in vitro-amplified cytb gene fragments as probes for direct DNA sequencing, three different types of single base pair substitutions were revealed in all mutants tested. In particular, amino acid substitutions in the mutant apocytochrome b polypeptide have been identified at residues 129, 132 and 137, which have been implicated in forming part of an antibiotic-binding niche. The amino acid substitution at position 132 has not been so far described for mutant apocytochrome b in any other organism, prokaryotic or eukaryotic. The genetic approach presented here confirms C. reinhardtii as a model system that is unique among plant cells.