Utilization of a whole genome SNP panel for efficient genetic mapping in the mouse

Phenotype-driven genetics can be used to create mouse models of human disease and birth defects. However, the utility of these mutant models is limited without identification of the causal gene. To facilitate genetic mapping, we developed a fixed single nucleotide polymorphism (SNP) panel of 394 SNPs as an alternative to analyses using simple sequence length polymorphism (SSLP) marker mapping. With the SNP panel, chromosomal locations for 22 monogenic mutants were identified. The average number of affected progeny genotyped for mapped monogenic mutations is nine. Map locations for several mutants have been obtained with as few as four affected progeny. The average size of genetic intervals obtained for these mutants is 43 Mb, with a range of 17-83 Mb. Thus, our SNP panel allows for identification of moderate resolution map position with small numbers of mice in a high-throughput manner. Importantly, the panel is suitable for mapping crosses from many inbred and wild-derived inbred strain combinations. The chromosomal localizations obtained with the SNP panel allow one to quickly distinguish between potentially novel loci or remutations in known genes, and facilitates fine mapping and positional cloning. By using this approach, we identified DNA sequence changes in two ethylnitrosourea-induced mutants.     

Single nucleotide polymorphisms in wild isolates of Caenorhabditis elegans

Caenorhabditis elegans (isolate N2 from Bristol, UK) is the first animal of which the complete genome sequence was available. We sampled genomic DNA of natural isolates of C. elegans from four different locations (Australia, Germany, California, and Wisconsin) and found single nucleotide polymorphisms (SNPs) by comparing with the Bristol strain. SNPs are under-represented in coding regions, and many were found to be third base silent codon mutations. We tested 19 additional natural isolates for the presence and distribution of SNPs originally found in one of the four strains. Most SNPs are present in isolates from around the globe and thus are older than the latest contact between these strains. An exception is formed by an isolate from an island (Hawaii) that contains many unique SNPs, absent in the tested isolates from the rest of the world. It has been noticed previously that conserved genes (as defined by homology to genes in Saccharomyces cerevisiae) cluster in the chromosome centers. We found that the SNP frequency outside these regions is 4.5 times higher, supporting the notion of a higher rate of evolution of genes on the chromosome arms.     

Efficient target-selected mutagenesis in Caenorhabditis elegans: toward a knockout for every gene

Reverse genetic or gene-driven knockout approaches have contributed significantly to the success of model organisms for fundamental and biomedical research. Although various technologies are available for C. elegans, none of them scale very well for genome-wide application. To address this, we implemented a target-selected knockout approach that is based on random chemical mutagenesis and detection of single nucleotide mutations in genes of interest using high-throughput resequencing. A clonal library of 6144 EMS-mutagenized worms was established and screened, resulting in the identification of 1044 induced mutations in 109 Mbp, which translates into an average spacing between exonic mutations in the library of only 17 bp. We covered 25% of the open reading frames of 32 genes and identified one or more inactivating mutations (nonsense or splice site) in 84% of them. Extrapolation of our results indicates that nonsense mutations for >90% of all C. elegans genes are present in the library. To identify all of these mutations, one only needs to inspect those positions that--given the known specificity of the mutagen--can result in the introduction of a stop codon. We define these positions as nonsense introducing mutations (NIMs). The genome-wide collection of possible NIMs can be calculated for any organism with a sequenced genome and reduces the screening complexity by 200- to 2000-fold, depending on the organism and mutagen. For EMS-mutagenized C. elegans, there are only approximately 500,000 NIMs. We show that a NIM genotyping approach employing high-density microarrays can, in principle, be used for the genome-wide identification of C. elegans knockouts.     

A genetic dissociation of learning and recall in Caenorhabditis elegans

A learning event can be dissociated into 3 components: acquisition, storage, and recall. When the laboratory wild-type strain of Caenorhabditis elegans (N2 strain) is exposed to benzaldehyde in the absence of food, the worms display a reduction of their attractive response to this volatile odorant. This results from the association between benzaldehyde and a nutrient-deficient environment. Another wild-type isolate, the CB4856 strain, fails to display this decreased response to benzaldehyde after exposure to benzaldehyde in the absence of food. However, like the N2 strain, when tested to isoamyl alcohol after benzaldehyde conditioning, the CB4856 strain displays a decreased isoamyl alcohol response. Therefore, the CB4856 strain does not have an acquisition deficit, but it suffers from a recall deficit specific to benzaldehyde.     

Selective sweeps and parallel mutation in the adaptive recovery from deleterious mutation in Caenorhabditis elegans

Deleterious mutation poses a serious threat to human health and the persistence of small populations. Although adaptive recovery from deleterious mutation has been well-characterized in prokaryotes, the evolutionary mechanisms by which multicellular eukaryotes recover from deleterious mutation remain unknown. We applied high-throughput DNA sequencing to characterize genomic divergence patterns associated with the adaptive recovery from deleterious mutation using a Caenorhabditis elegans recovery-line system. The C. elegans recovery lines were initiated from a low-fitness mutation-accumulation (MA) line progenitor and allowed to independently evolve in large populations (N ～ 1000) for 60 generations. All lines rapidly regained levels of fitness similar to the wild-type (N2) MA line progenitor. Although there was a near-zero probability of a single mutation fixing due to genetic drift during the recovery experiment, we observed 28 fixed mutations. Cross-generational analysis showed that all mutations went from undetectable population-level frequencies to a fixed state in 10-20 generations. Many recovery-line mutations fixed at identical timepoints, suggesting that the mutations, if not beneficial, hitchhiked to fixation during selective sweep events observed in the recovery lines. No MA line mutation reversions were detected. Parallel mutation fixation was observed for two sites in two independent recovery lines. Analysis using a C. elegans interactome map revealed many predicted interactions between genes with recovery line-specific mutations and genes with previously accumulated MA line mutations. Our study suggests that recovery-line mutations identified in both coding and noncoding genomic regions might have beneficial effects associated with compensatory epistatic interactions.     

Efficient high-resolution deletion discovery in Caenorhabditis elegans by array comparative genomic hybridization

We have developed array Comparative Genomic Hybridization for Caenorhabditis elegans as a means of screening for novel induced deletions in this organism. We designed three microarrays consisting of overlapping 50-mer probes to annotated exons and micro-RNAs, the first with probes to chromosomes X and II, the second with probes to chromosome II alone, and a third to the entire genome. These arrays were used to reliably detect both a large (50 kb) multigene deletion and a small (1 kb) single-gene deletion in homozygous and heterozygous samples. In one case, a deletion breakpoint was resolved to fewer than 50 bp. In an experiment designed to identify new mutations we used the X:II and II arrays to detect deletions associated with lethal mutants on chromosome II. One is an 8-kb deletion targeting the ast-1 gene on chromosome II and another is a 141-bp deletion in the gene C06A8.1. Others span large sections of the chromosome, up to >750 kb. As a further application of array Comparative Genomic Hybridization in C. elegans we used the whole-genome array to detect the extensive natural gene content variation (almost 2%) between the N2 Bristol strain and the strain CB4856, a strain isolated in Hawaii and JU258, a strain isolated in Madeira.     

A cSNP map and database for human chromosome 21

Single nucleotide polymorphisms (SNPs) are likely to contribute to the study of complex genetic diseases. The genomic sequence of human chromosome 21q was recently completed with 225 annotated genes, thus permitting efficient identification and precise mapping of potential cSNPs by bioinformatics approaches. Here we present a human chromosome 21 (HC21) cSNP database and the first chromosome-specific cSNP map. Potential cSNPs were generated using three approaches: (1) Alignment of the complete HC21 genomic sequence to cognate ESTs and mRNAs. Candidate cSNPs were automatically extracted using a novel program for context-dependent SNP identification that efficiently discriminates between true variation, poor quality sequencing, and paralogous gene alignments. (2) Multiple alignment of all known HC21 genes to all other human database entries. (3) Gene-targeted cSNP discovery. To date we have identified 377 cSNPs averaging ~1 SNP per 1.5 kb of transcribed sequence, covering 65% of known genes in the chromosome. Validation of our bioinformatics approach was demonstrated by a confirmation rate of 78% for the predicted cSNPs, and in total 32% of the cSNPs in our database have been confirmed. The database is publicly available at http://csnp.unige.ch or http://csnp.isb-sib.ch. These SNPs provide a tool to study the contribution of HC21 loci to complex diseases such as bipolar affective disorder and allele-specific contributions to Down syndrome phenotypes.     

线虫基因组功能的高通量研究

线虫是研究细胞生物学的一个理想模型,RNA i现象也是最早在线虫中发现的。由于RNA i技术不仅具有高度特异性靶定基因的特点,而且简便、易于操作,从而成为研究基因组功能的高通量方法。很多研究小组已经使用大规模RNA i系统地研究线虫各个基因的功能缺失表型。线虫独特的RNA i方法和以此为基础的基因组大规模分析有利于对线虫的各种生物学过程产生新的认识。     

Pattern of Sequence Variation Across 213 Environmental Response Genes

To promote the clinical and epidemiological studies that improve our understanding of human genetic susceptibility to environmental exposure, the Environmental Genome Project (EGP) has scanned 213 environmental response genes involved in DNA repair, cell cycle regulation, apoptosis, and metabolism for single nucleotide polymorphisms (SNPs). Many of these genes have been implicated by loss-of-function mutations associated with severe diseases attributable to decreased protection of genomic integrity. Therefore, the hypothesis for these studies is that individuals with functionally significant polymorphisms within these genes may be particularly susceptible to genotoxic environmental agents. On average, 20.4 kb of baseline genomic sequence or 86% of each gene, including a substantial amount of introns, all exons, and 1.3 kb upstream and downstream, were scanned for variations in the 90 samples of the Polymorphism Discovery Resource panel. The average nucleotide diversity across the 4.2 MB of these 213 genes is 6.7 × 10^-4, or one SNP every 1500 bp, when two random chromosomes are compared. The average candidate environmental response gene contains 26 PHASE inferred haplotypes, 34 common SNPs, 6.2 coding SNPs (cSNPs), and 2.5 nonsynonymous cSNPs. SIFT and Polyphen analysis of 541 nonsynonymous cSNPs identified 57 potentially deleterious SNPs. An additional eight polymorphisms predict altered protein translation. Because these genes represent 1% of all known human genes, extrapolation from these data predicts the total genomic set of cSNPs, nonsynonymous cSNPs, and potentially deleterious nonsynonymous cSNPs. The implications for the use of these data in direct and indirect association studies of environmentally induced diseases are discussed.     

novoSNP, a novel computational tool for sequence variation discovery

Technological improvements shifted sequencing from low-throughput, work-intensive, gel-based systems to high-throughput capillary systems. This resulted in a broad use of genomic resequencing to identify sequence variations in genes and regulatory, as well as extended genomic regions. We describe a software package, novoSNP, that conscientiously discovers single nucleotide polymorphisms (SNPs) and insertion-deletion polymorphisms (INDELs) in sequence trace files in a fast, reliable, and user-friendly way. We compared the performance of novoSNP with that of PolyPhred and PolyBayes on two data sets. The first data set comprised 1028 sequence trace files obtained from diagnostic mutation analyses of SCN1A (neuronal voltage-gated sodium channel alpha-subunit type I gene). The second data set comprised 9062 sequence trace files from a genomic resequencing project aiming at the construction of a high-density SNP map of MAPT (microtubule-associated protein tau gene). Visual inspection of these data sets had identified 38 sequence variations for SCN1A and 488 for MAPT. novoSNP automatically identified all 38 SCN1A variations including five INDELs, while for MAPT only 15 of the 488 variations were not correctly marked. PolyPhred detected far fewer SNPs as compared to novoSNP and missed nearly all INDELs. PolyBayes, designed for the sequence analysis of cloned templates, detected only a limited number of the variations present in the data set. Besides the significant improvement in the automated detection of sequence variations both in diagnostic mutation analyses and in SNP discovery projects, novoSNP also offers a user-friendly interface for inspecting possible genetic variations.     

Novel antigen identification method for discovery of protective malaria antigens by rapid testing of DNA vaccines encoding exons from the parasite genome

We describe a novel approach for identifying target antigens for preerythrocytic malaria vaccines. Our strategy is to rapidly test hundreds of DNA vaccines encoding exons from the Plasmodium yoelii yoelii genomic sequence. In this antigen identification method, we measure reduction in parasite burden in the liver after sporozoite challenge in mice. Orthologs of protective P. y. yoelii genes can then be identified in the genomic databases of Plasmodium falciparum and Plasmodium vivax and investigated as candidate antigens for a human vaccine. A pilot study to develop the antigen identification method approach used 192 P. y. yoelii exons from genes expressed during the sporozoite stage of the life cycle. A total of 182 (94%) exons were successfully cloned into a DNA immunization vector with the Gateway cloning technology. To assess immunization strategies, mice were vaccinated with 19 of the new DNA plasmids in addition to the well-characterized protective plasmid encoding P. y. yoelii circumsporozoite protein. Single plasmid immunization by gene gun identified a novel vaccine target antigen which decreased liver parasite burden by 95% and which has orthologs in P. vivax and P. knowlesi but not P. falciparum. Intramuscular injection of DNA plasmids produced a different pattern of protective responses from those seen with gene gun immunization. Intramuscular immunization with plasmid pools could reduce liver parasite burden in mice despite the fact that none of the plasmids was protective when given individually. We conclude that high-throughput cloning of exons into DNA vaccines and their screening is feasible and can rapidly identify new malaria vaccine candidate antigens.     

New developmental insights from high-throughput biological analysis in Caenorhabditis elegans

The use of Caenorhabditis elegans as a model system for understanding animal development and human disease has long been recognized as an efficient tool of discovery. Recent developments, particularly in our understanding of RNA-mediated interference and its ability to modify gene activity, have facilitated the use of C. elegans in determining gene function via high-throughput analysis. These new strategies have provided a framework that allows investigators to analyse gene function globally at the genomic level and will likely become a prototypic model for biological analysis in the post-genome era.     

Novel genomic techniques open new avenues in the analysis of monogenic disorders

The molecular genetic cause of over 3,000 monogenic disorders is currently unknown. This review discusses how novel genomic techniques like Next‐Generation DNA Sequencing (NGS) and genotyping arrays open new avenues in the elucidation of genetic defects causing monogenic disorders. They will not only speed up disease gene identification but will enable us to systematically tackle previously intractable monogenic disorders. These are mainly disorders not amenable to classic linkage analysis, for example, due to insufficient family size. Most monogenic diseases are caused by exonic mutations or splice‐site mutations changing the amino acid sequence of the affected gene. These mutations can be identified by sequencing of all exons in the human genome (exome sequencing) rendering whole genome sequencing unnecessary in most cases. Genotyping arrays containing 10⁵–2×10⁶ single nucleotide polymorphisms (SNPs) and nonpolymorphic markers allow highly accurate mapping of genomic deletions and duplications not detectable by exome sequencing, which are the second most common cause of monogenic disorders. However, several hundred rare, previously unknown sequence variants affecting the amino acid sequence of the encoded protein are found in the exome of every human individual. Therefore, the main challenge will be the differentiation between the many rare benign variants detected by novel genomic techniques and disease causing mutations. Hum Mutat 32:144–151, 2011. © 2011 Wiley‐Liss, Inc.     

Scaffolding a Caenorhabditis nematode genome with RNA-seq

Efficient sequencing of animal and plant genomes by next-generation technology should allow many neglected organisms of biological and medical importance to be better understood. As a test case, we have assembled a draft genome of Caenorhabditis sp. 3 PS1010 through a combination of direct sequencing and scaffolding with RNA-seq. We first sequenced genomic DNA and mixed-stage cDNA using paired 75-nt reads from an Illumina GAII. A set of 230 million genomic reads yielded an 80-Mb assembly, with a supercontig N50 of 5.0 kb, covering 90% of 429 kb from previously published genomic contigs. Mixed-stage poly(A)(+) cDNA gave 47.3 million mappable 75-mers (including 5.1 million spliced reads), which separately assembled into 17.8 Mb of cDNA, with an N50 of 1.06 kb. By further scaffolding our genomic supercontigs with cDNA, we increased their N50 to 9.4 kb, nearly double the average gene size in C. elegans. We predicted 22,851 protein-coding genes, and detected expression in 78% of them. Multigenome alignment and data filtering identified 2672 DNA elements conserved between PS1010 and C. elegans that are likely to encode regulatory sequences or previously unknown ncRNAs. Genomic and cDNA sequencing followed by joint assembly is a rapid and useful strategy for biological analysis.     

High-throughput semiquantitative analysis of insertional mutations in heterogeneous tumors

Retroviral and transposon-based insertional mutagenesis (IM) screens are widely used for cancer gene discovery in mice. Exploiting the full potential of IM screens requires methods for high-throughput sequencing and mapping of transposon and retroviral insertion sites. Current protocols are based on ligation-mediated PCR amplification of junction fragments from restriction endonuclease-digested genomic DNA, resulting in amplification biases due to uneven genomic distribution of restriction enzyme recognition sites. Consequently, sequence coverage cannot be used to assess the clonality of individual insertions. We have developed a novel method, called shear-splink, for the semiquantitative high-throughput analysis of insertional mutations. Shear-splink employs random fragmentation of genomic DNA, which reduces unwanted amplification biases. Additionally, shear-splink enables us to assess clonality of individual insertions by determining the number of unique ligation points (LPs) between the adapter and genomic DNA. This parameter serves as a semiquantitative measure of the relative clonality of individual insertions within heterogeneous tumors. Mixing experiments with clonal cell lines derived from mouse mammary tumor virus (MMTV)-induced tumors showed that shear-splink enables the semiquantitative assessment of the clonality of MMTV insertions. Further, shear-splink analysis of 16 MMTV- and 127 Sleeping Beauty (SB)-induced tumors showed enrichment for cancer-relevant insertions by exclusion of irrelevant background insertions marked by single LPs, thereby facilitating the discovery of candidate cancer genes. To fully exploit the use of the shear-splink method, we set up the Insertional Mutagenesis Database (iMDB), offering a publicly available web-based application to analyze both retroviral- and transposon-based insertional mutagenesis data.     

Genomics of the future: identification of quantitative trait loci in the mouse

Positional cloning of quantitative trait loci in rodents is a common approach to identify genes involved in complex phenotypes, including genes important to human disease. However, cloning the causative genes has proved to be more difficult than determining their positions. New tools such as genomic sequence, clone libraries, and new genomic-based methods offer new approaches to identify these genes. Here we review how these new tools and approaches will improve our ability to discover the genes important in complex traits.     

Haplotype structure and association to Crohn's disease of CARD15 mutations in two ethnically divergent populations

Current debate focuses on the relevance of linkage disequilibrium (LD), ethnicity and underlying haplotype structure to the search for genes involved in complex disorders. The recently described association between single nucleotide polymorphisms (SNPs) of the CARD15 (NOD2) gene and Crohn's disease (CD) in populations of north-European descent provides a test case that we have subjected to detailed SNP haplotype based analyses. We examined 23 SNPs spanning 290 kb, including CARD15, in large North-European and Korean samples of patients with Crohn's disease and normal controls. In Europeans we confirmed that the three disease-associated SNPs occur independently but share a common background haplotype. This suggests a common origin and the possibility of an undiscovered more strongly predisposing mutation. Korean CD patients present a phenotype identical to the European patients and have not previously been screened for CARD15. The three disease-associated SNPs were absent and there was no evidence of association between CARD15 and CD. Consequently, the disease-associated mutations in the Europeans, which are rare, have arisen recently (after the Asian-European split). Our results highlight important issues relevant to mapping the genes that predispose to complex disorders. First, although ethnically divergent populations may present identical phenotypes they do not necessarily share the same set of predisposing genes. Second, although single-locus tests of association showed consistent association with markers throughout the gene, pair-wise LD between markers (r(2) and D') yielded very little information about actual disease-association. Third, a population comparative approach allowed refining of the marker set through the examination of shared polymorphisms and common LD-groups. This approach, in conjunction with the examination of the mutational steps in a haplotype network, allows unambiguous identification of the potentially causative mutations.