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.     

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.     

Targeted next-generation sequencing of a 12.5 Mb homozygous region reveals ANO10 mutations in patients with autosomal-recessive cerebellar ataxia

Autosomal-recessive cerebellar ataxias comprise a clinically and genetically heterogeneous group of neurodegenerative disorders. In contrast to their dominant counterparts, unraveling the molecular background of these ataxias has proven to be more complicated and the currently known mutations provide incomplete coverage for genotyping of patients. By combining SNP array-based linkage analysis and targeted resequencing of relevant sequences in the linkage interval with the use of next-generation sequencing technology, we identified a mutation in a gene and have shown its association with autosomal-recessive cerebellar ataxia. In a Dutch consanguineous family with three affected siblings a homozygous 12.5 Mb region on chromosome 3 was targeted by array-based sequence capture. Prioritization of all detected sequence variants led to four candidate genes, one of which contained a variant with a high base pair conservation score (phyloP score: 5.26). This variant was a leucine-to-arginine substitution in the DUF 590 domain of a 16K transmembrane protein, a putative calcium-activated chloride channel encoded by anoctamin 10 (ANO10). The analysis of ANO10 by Sanger sequencing revealed three additional mutations: a homozygous mutation (c.1150_1151del [p.Leu384fs]) in a Serbian family and a compound-heterozygous splice-site mutation (c.1476+1G>T) and a frameshift mutation (c.1604del [p.Leu535X]) in a French family. This illustrates the power of using initial homozygosity mapping with next-generation sequencing technology to identify genes involved in autosomal-recessive diseases. Moreover, identifying a putative calcium-dependent chloride channel involved in cerebellar ataxia adds another pathway to the list of pathophysiological mechanisms that may cause cerebellar ataxia.     

First TILLING Platform in Cucurbita pepo: A New Mutant Resource for Gene Function and Crop Improvement

Although the availability of genetic and genomic resources for Cucurbita pepo has increased significantly, functional genomic resources are still limited for this crop. In this direction, we have developed a high throughput reverse genetic tool: the first TILLING (Targeting Induced Local Lesions IN Genomes) resource for this species. Additionally, we have used this resource to demonstrate that the previous EMS mutant population we developed has the highest mutation density compared with other cucurbits mutant populations. The overall mutation density in this first C. pepo TILLING platform was estimated to be 1/133 Kb by screening five additional genes. In total, 58 mutations confirmed by sequencing were identified in the five targeted genes, thirteen of which were predicted to have an impact on the function of the protein. The genotype/phenotype correlation was studied in a peroxidase gene, revealing that the phenotype of seedling homozygous for one of the isolated mutant alleles was albino. These results indicate that the TILLING approach in this species was successful at providing new mutations and can address the major challenge of linking sequence information to biological function and also the identification of novel variation for crop breeding.     

Whole-Genome Sequencing of a Single Proband Together with Linkage Analysis Identifies a Mendelian Disease Gene

Although more than 2,400 genes have been shown to contain variants that cause Mendelian disease, there are still several thousand such diseases yet to be molecularly defined. The ability of new whole-genome sequencing technologies to rapidly indentify most of the genetic variants in any given genome opens an exciting opportunity to identify these disease genes. Here we sequenced the whole genome of a single patient with the dominant Mendelian disease, metachondromatosis (OMIM 156250), and used partial linkage data from her small family to focus our search for the responsible variant. In the proband, we identified an 11 bp deletion in exon four of PTPN11, which alters frame, results in premature translation termination, and co-segregates with the phenotype. In a second metachondromatosis family, we confirmed our result by identifying a nonsense mutation in exon 4 of PTPN11 that also co-segregates with the phenotype. Sequencing PTPN11 exon 4 in 469 controls showed no such protein truncating variants, supporting the pathogenicity of these two mutations. This combination of a new technology and a classical genetic approach provides a powerful strategy to discover the genes responsible for unexplained Mendelian disorders.     

A clonal analysis of the roles of somatic cells and germ line during oogenesis in Drosophila

In order to test whether particular female sterile mutations block functions which normally occur in somatic or germ line derivatives, clones homozygous for each mutation were X-ray induced in heterozygous females. Using the germ line-dependent egg marker, fs(1)K10, it was possible to identify the eggs derived from clones which had been induced in the germ line. Mutations were classified as germ line dependent when these eggs also showed the phenotype associated with the female sterile mutation. Two mutations which caused early abnormalities in oogenesis (fs(1)116, fs(1)1304) were shown to affect germ cells, whereas two mutations which caused egg retention (fs(1)462, fs(1)1001) were somatically dependent. A mutation altering egg dimensions without affecting egg volume (short egg) was also shown to depend on somatic cells in the ovary. With one exception (fs(1)K4), mutations which caused production of fragile, collapsed eggs (fs(1)180, fs(1)473, fs(1)384, and fs(1)1163) were somatically dependent. Patches of mutant fs(1)384 morphology were found in the chorions of the eggs not derived from germ line clones. These patches are interpreted as being caused by homozygous clones in the somatically derived follicle cell epithelium and suggest that fs(1)384 affects processes occurring in these cells during the synthesis of the egg coverings.     

Genetic mapping of glutaric aciduria, type 3, to chromosome 7 and identification of mutations in c7orf10

While screening Old Order Amish children for glutaric aciduria type 1 (GA1) between 1989 and 1993, we found three healthy children who excreted abnormal quantities of glutaric acid but low 3-hydroxyglutaric acid, a pattern consistent with glutaric aciduria type 3 (GA3). None of these children had the GCDH c.1262C→T mutation that causes GA1 among the Amish. Using single-nucleotide polymorphism (SNP) genotypes, we identified a shared homozygous 4.7 Mb region on chromosome 7. This region contained 25 genes including C7orf10, an open reading frame with a putative mitochondrial targeting sequence and coenzyme-A transferase domain. Direct sequencing of C7orf10 revealed that the three Amish individuals were homozygous for a nonsynonymous sequence variant (c.895C→T, Arg299Trp). We then sequenced three non-Amish children with GA3 and discovered two nonsense mutations (c.322C→T, Arg108Ter, and c.424C→T, Arg142Ter) in addition to the Amish mutation. Two pathogenic alleles were identified in each of the six patients. There was no consistent clinical phenotype associated with GA3. In affected individuals, urine molar ratios of glutarate to its derivatives (3-hydroxyglutarate, glutarylcarnitine, and glutarylglycine) were elevated, suggesting impaired formation of glutaryl-CoA. These observations refine our understanding of the lysine-tryptophan degradation pathway and have important implications for the pathophysiology of GA1.     

Identification of a Nonsense Mutation in CWC15 Associated with Decreased Reproductive Efficiency in Jersey Cattle

With the recent advent of genomic tools for cattle, several recessive conditions affecting fertility have been identified and selected against, such as deficiency of uridine monophosphate synthase, complex vertebral malformation, and brachyspina. The current report refines the location of a recessive haplotype affecting fertility in Jersey cattle using crossover haplotypes, discovers the causative mutation using whole genome sequencing, and examines the gene’s role in embryo loss. In an attempt to identify unknown recessive lethal alleles in the current dairy population, a search using deep Mendelian sampling of 5,288 Jersey cattle was conducted for high-frequency haplotypes that have a deficit of homozygotes at the population level. This search led to the discovery of a putative recessive lethal in Jersey cattle on Bos taurus autosome 15. The haplotype, denoted JH1, was associated with reduced fertility, and further investigation identified one highly-influential Jersey bull as the putative source ancestor. By combining SNP analysis of whole-genome sequences aligned to the JH1 interval and subsequent SNP validation a nonsense mutation in CWC15 was identified as the likely causative mutation underlying the fertility phenotype. No homozygous recessive individuals were found in 749 genotyped animals, whereas all known carriers and carrier haplotypes possessed one copy of the mutant allele. This newly identified lethal has been responsible for a substantial number of spontaneous abortions in Jersey dairy cattle throughout the past half-century. With the mutation identified, selection against the deleterious allele in breeding schemes will aid in reducing the incidence of this defect in the population. These results also show that carrier status can be imputed with high accuracy. Whole-genome resequencing proved to be a powerful strategy to rapidly identify a previously mapped deleterious mutation in a known carrier of a recessive lethal allele.     

Next-generation sequencing as a tool to quickly identify causative EMS-generated mutations

The advent of next generation sequencing has influenced every aspect of biological research. Many labs are now using whole genome sequencing in Arabidopsis thaliana as a means to quickly identify EMS-generated mutations present in isolated mutants. Following identification of these mutations, examination of T-DNA insertional alleles defective in candidate genes or complementation of the mutant phenotype with a wild type copy of candidate genes can be used to verify which mutation is causative for the phenotype of interest. Here, we discuss the benefits and pitfalls of using this method to identify mutations underlying phenotypes.     

Use of Targeted Exome Sequencing in Genetic Diagnosis of Chinese Familial Hypercholesterolemia

Familial hypercholesterolemia is an autosomal dominant inherited disease characterized by elevated plasma low-density lipoprotein cholesterol (LDL-C). It is mainly caused by mutations of the low-density lipoprotein receptor (LDLR) gene. Currently, the methods of whole genome sequencing or whole exome sequencing for screening mutations in familial hypercholesterolemia are not applicable in China due to high cost. We performed targeted exome sequencing of 167 genes implicated in the homozygous phenotype of a proband pedigree to identify candidate mutations, validated them in the family of the proband, studied the functions of the mutant protein, and followed up serum lipid levels after treatment. We discovered that exon 9 c.1268 T>C and exon 8 c.1129 T>G compound heterozygous mutations in the LDLR gene in the proband derived from the mother and father, respectively, in which the mutation of c.1129 T>G has not been reported previously. The mutant LDL-R protein had 57% and 52% binding and internalization functions, respectively, compared with that of the wild type. After 6 months of therapy, the LDL-C level of the proband decreased by more than 50% and the LDL-C of the other family members with heterozygous mutation also reduced to normal. Targeted exome sequencing is an effective method for screening mutation genes in familial hypercholesterolemia. The exon 8 and 9 mutations of the LDLR gene were pedigree mutations. The functions of the mutant LDL-R protein were decreased significantly compared with that of the wild type. Simvastatin plus ezetimibe was proven safe and effective in this preschool-age child.     

Discovery of Candidate Disease Genes in ENU–Induced Mouse Mutants by Large-Scale Sequencing,Including a Splice-Site Mutation in Nucleoredoxin

An accurate and precisely annotated genome assembly is a fundamental requirement for functional genomic analysis. Here, the complete DNA sequence and gene annotation of mouse Chromosome 11 was used to test the efficacy of large-scale sequencing for mutation identification. We re-sequenced the 14,000 annotated exons and boundaries from over 900 genes in 41 recessive mutant mouse lines that were isolated in an N-ethyl-N-nitrosourea (ENU) mutation screen targeted to mouse Chromosome 11. Fifty-nine sequence variants were identified in 55 genes from 31 mutant lines. 39% of the lesions lie in coding sequences and create primarily missense mutations. The other 61% lie in noncoding regions, many of them in highly conserved sequences. A lesion in the perinatal lethal line l11Jus13 alters a consensus splice site of nucleoredoxin (Nxn), inserting 10 amino acids into the resulting protein. We conclude that point mutations can be accurately and sensitively recovered by large-scale sequencing, and that conserved noncoding regions should be included for disease mutation identification. Only seven of the candidate genes we report have been previously targeted by mutation in mice or rats, showing that despite ongoing efforts to functionally annotate genes in the mammalian genome, an enormous gap remains between phenotype and function. Our data show that the classical positional mapping approach of disease mutation identification can be extended to large target regions using high-throughput sequencing.     

Novel splice, missense, and nonsense mutations in the fumarylacetoacetase gene causing tyrosinemia type 1.

In six unrelated patients with hereditary tyrosinemia type 1 (HT1), three different disease-causing mutations were found by DNA sequencing. Two Pakistani patients, with acute and intermediate forms of HT1, were homozygous for a G192-->T mutation in the last nucleotide of exon 2. This caused aberrant splicing with partial intron 2 retention and premature termination. Three Turkish patients with chronic and intermediate forms of HT1 were homozygous for an A698-->T mutation substituting aspartic acid 233 with valine. A Norwegian patient with an intermediate clinical phenotype was heterozygous for G786-->A, introducing a TGA stop codon for Trp262 (W262X). Site-directed mutagenesis and expression in a rabbit reticulocyte lysate system demonstrated that the nonsense and missense mutations abolished fumarylacetoacetase activity and gave reduced amounts of a truncated and a full-length protein, respectively. Simple tests were established to identify the three mutations by restriction digestion of PCR-amplified genomic DNA. Among 30 additional HT1 patients investigated, 2 were found to be homozygous and 1 heterozygous for G192-->T. Two other patients were homozygous and one was heterozygous for W262X.     

Traffic Lines: New Tools for Genetic Analysis in Arabidopsis thaliana

Genetic analysis requires the ability to identify the genotypes of individuals in a segregating population. This task is straightforward if each genotype has a distinctive phenotype, but is difficult if these genotypes are phenotypically similar or identical. We show that Arabidopsis seeds homozygous or heterozygous for a mutation of interest can be identified in a segregating family by placing the mutation in trans to a chromosome carrying a pair of seed-expressed green and red fluorescent transgenes (a “traffic line”) that flank the mutation. Nonfluorescent seeds in the self-pollinated progeny of such a heterozygous plant are usually homozygous for the mutation, whereas seeds with intermediate green and red fluorescence are typically heterozygous for the mutation. This makes it possible to identify seedlings homozygous for mutations that lack an obvious seedling phenotype, and also facilitates the analysis of lethal or sterile mutations, which must be propagated in heterozygous condition. Traffic lines can also be used to identify progeny that have undergone recombination within a defined region of the genome, facilitating genetic mapping and the production of near-isogenic lines. We produced 488 transgenic lines containing single genome-mapped insertions of NAP:dsRED and NAP:eGFP in Columbia (330 lines) and Landsberg erecta (158 lines) and generated sets of traffic lines that span most regions of the Arabidopsis genome. We demonstrated the utility of these lines for identifying seeds of a specific genotype and for generating near-isogenic lines using mutations of WUSCHEL and SHOOTMERISTEMLESS. This new resource significantly decreases the effort and cost of genotyping segregating families and increases the efficiency of experiments that rely on the ability to detect recombination in a defined chromosomal segment.     

The phenotype of a germline mutation in PIGA: the gene somatically mutated in paroxysmal nocturnal hemoglobinuria

Phosphatidylinositol glycan class A (PIGA) is involved in the first step of glycosylphosphatidylinositol (GPI) biosynthesis. Many proteins, including CD55 and CD59, are anchored to the cell by GPI. Loss of CD55 and CD59 on erythrocytes causes complement-mediated lysis in paroxysmal nocturnal hemoglobinuria (PNH), a disease that manifests after clonal expansion of hematopoietic cells with somatic PIGA mutations. Although somatic PIGA mutations have been identified in many PNH patients, it has been proposed that germline mutations are lethal. We report a family with an X-linked lethal disorder involving cleft palate, neonatal seizures, contractures, central nervous system (CNS) structural malformations, and other anomalies. An X chromosome exome next-generation sequencing screen identified a single nonsense PIGA mutation, c.1234C>T, which predicts p.Arg412^∗. This variant segregated with disease and carrier status in the family, is similar to mutations known to cause PNH as a result of PIGA dysfunction, and was absent in 409 controls. PIGA-null mutations are thought to be embryonic lethal, suggesting that p.Arg412^∗ PIGA has residual function. Transfection of a mutant p.Arg412^∗ PIGA construct into PIGA-null cells showed partial restoration of GPI-anchored proteins. The genetic data show that the c.1234C>T (p.Arg412^∗) mutation is present in an affected child, is linked to the affected chromosome in this family, is rare in the population, and results in reduced, but not absent, biosynthesis of GPI anchors. We conclude that c.1234C>T in PIGA results in the lethal X-linked phenotype recognized in the reported family.     

Comparing Platforms for C. elegans Mutant Identification Using High-Throughput Whole-Genome Sequencing

Background

Whole-genome sequencing represents a promising approach to pinpoint chemically induced mutations in genetic model organisms, thereby short-cutting time-consuming genetic mapping efforts.

Principal Findings

We compare here the ability of two leading high-throughput platforms for paired-end deep sequencing, SOLiD (ABI) and Genome Analyzer (Illumina; “Solexa”), to achieve the goal of mutant detection. As a test case we used a mutant C. elegans strain that harbors a mutation in the lsy-12 locus which we compare to the reference wild-type genome sequence. We analyzed the accuracy, sensitivity, and depth-coverage characteristics of the two platforms. Both platforms were able to identify the mutation that causes the phenotype of the mutant C. elegans strain, lsy-12. Based on a 4 MB genomic region in which individual variants were validated by Sanger sequencing, we observe tradeoffs between rates of false positives and false negatives when using both platforms under similar coverage and mapping criteria.

Significance

In conclusion, whole-genome sequencing conducted by either platform is a viable approach for the identification of single-nucleotide variations in the C. elegans genome.     

High-precision, whole-genome sequencing of laboratory strains facilitates genetic studies

Whole-genome sequencing is a powerful technique for obtaining the reference sequence information of multiple organisms. Its use can be dramatically expanded to rapidly identify genomic variations, which can be linked with phenotypes to obtain biological insights. We explored these potential applications using the emerging next-generation sequencing platform Solexa Genome Analyzer, and the well-characterized model bacterium Bacillus subtilis. Combining sequencing with experimental verification, we first improved the accuracy of the published sequence of the B. subtilis reference strain 168, then obtained sequences of multiple related laboratory strains and different isolates of each strain. This provides a framework for comparing the divergence between different laboratory strains and between their individual isolates. We also demonstrated the power of Solexa sequencing by using its results to predict a defect in the citrate signal transduction pathway of a common laboratory strain, which we verified experimentally. Finally, we examined the molecular nature of spontaneously generated mutations that suppress the growth defect caused by deletion of the stringent response mediator relA. Using whole-genome sequencing, we rapidly mapped these suppressor mutations to two small homologs of relA. Interestingly, stable suppressor strains had mutations in both genes, with each mutation alone partially relieving the relA growth defect. This supports an intriguing three-locus interaction module that is not easily identifiable through traditional suppressor mapping. We conclude that whole-genome sequencing can drastically accelerate the identification of suppressor mutations and complex genetic interactions, and it can be applied as a standard tool to investigate the genetic traits of model organisms.     

Mutations in IMPG2, Encoding Interphotoreceptor Matrix Proteoglycan 2, Cause Autosomal-Recessive Retinitis Pigmentosa

Retinitis pigmentosa (RP) is a heterogeneous group of inherited retinal diseases caused by progressive degeneration of the photoreceptor cells. Using autozygosity mapping, we identified two families, each with three affected siblings sharing large overlapping homozygous regions that harbored the IMPG2 gene on chromosome 3. Sequence analysis of IMPG2 in the two index cases revealed homozygous mutations cosegregating with the disease in the respective families: three affected siblings of Iraqi Jewish ancestry displayed a nonsense mutation, and a Dutch family displayed a 1.8 kb genomic deletion that removes exon 9 and results in the absence of seven amino acids in a conserved SEA domain of the IMPG2 protein. Transient transfection of COS-1 cells showed that a construct expressing the wild-type SEA domain is properly targeted to the plasma membrane, whereas the mutant lacking the seven amino acids appears to be retained in the endoplasmic reticulum. Mutation analysis in ten additional index cases that were of Dutch, Israeli, Italian, and Pakistani origin and had homozygous regions encompassing IMPG2 revealed five additional mutations; four nonsense mutations and one missense mutation affecting a highly conserved phenylalanine residue. Most patients with IMPG2 mutations showed an early-onset form of RP with progressive visual-field loss and deterioration of visual acuity. The patient with the missense mutation, however, was diagnosed with maculopathy. The IMPG2 gene encodes the interphotoreceptor matrix proteoglycan IMPG2, which is a constituent of the interphotoreceptor matrix. Our data therefore show that mutations in a structural component of the interphotoreceptor matrix can cause arRP.     

Mutations in CSTA, encoding Cystatin A, underlie exfoliative ichthyosis and reveal a role for this protease inhibitor in cell-cell adhesion

Autosomal-recessive exfoliative ichthyosis presents shortly after birth as dry, scaly skin over most of the body with coarse peeling of nonerythematous skin on the palms and soles, which is exacerbated by excessive moisture and minor trauma. Using whole-genome homozygosity mapping, candidate-gene analysis and deep sequencing, we have identified loss-of-function mutations in the gene for protease inhibitor cystatin A (CSTA) as the underlying genetic cause of exfoliative ichthyosis. We found two homozygous mutations, a splice-site and a nonsense mutation, in two consanguineous families of Bedouin and Turkish origin. Electron microscopy of skin biopsies from affected individuals revealed that the level of detachment occurs in the basal and lower suprabasal layers. In addition, in vitro modeling suggests that in the absence of cystatin A protein, there is a cell-cell adhesion defect in human keratinocytes that is particularly prominent when cells are subject to mechanical stress. We show here evidence of a key role for a protease inhibitor in epidermal adhesion within the lower layers of the human epidermis.