Whole-genome profiling of mutagenesis in Caenorhabditis elegans

CloudMap: a cloud-based pipeline for analysis of mutant genome sequences

Whole genome sequencing (WGS) allows researchers to pinpoint genetic differences between individuals and significantly shortcuts the costly and time-consuming part of forward genetic analysis in model organism systems. Currently, the most effort-intensive part of WGS is the bioinformatic analysis of the relatively short reads generated by second generation sequencing platforms. We describe here a novel, easily accessible and cloud-based pipeline, called CloudMap, which greatly simplifies the analysis of mutant genome sequences. Available on the Galaxy web platform, CloudMap requires no software installation when run on the cloud, but it can also be run locally or via Amazon's Elastic Compute Cloud (EC2) service. CloudMap uses a series of predefined workflows to pinpoint sequence variations in animal genomes, such as those of premutagenized and mutagenized Caenorhabditis elegans strains. In combination with a variant-based mapping procedure, CloudMap allows users to sharply define genetic map intervals graphically and to retrieve very short lists of candidate variants with a few simple clicks. Automated workflows and extensive video user guides are available to detail the individual analysis steps performed (http://usegalaxy.org/cloudmap). We demonstrate the utility of CloudMap for WGS analysis of C. elegans and Arabidopsis genomes and describe how other organisms (e.g., Zebrafish and Drosophila) can easily be accommodated by this software platform. To accommodate rapid analysis of many mutants from large-scale genetic screens, CloudMap contains an in silico complementation testing tool that allows users to rapidly identify instances where multiple alleles of the same gene are present in the mutant collection. Lastly, we describe the application of a novel mapping/WGS method ("Variant Discovery Mapping") that does not rely on a defined polymorphic mapping strain, and we integrate the application of this method into CloudMap. CloudMap tools and documentation are continually updated at http://usegalaxy.org/cloudmap.

RNA-seq analysis of the C. briggsae transcriptome

Curation of a high-quality gene set is the critical first step in genome research, enabling subsequent analyses such as ortholog assignment, cis-regulatory element finding, and synteny detection. In this project, we have reannotated the genome of Caenorhabditis briggsae, the best studied sister species of the model organism Caenorhabditis elegans. First, we applied a homology-based gene predictor genBlastG to annotate the C. briggsae genome. We then validated and further improved the C. briggsae gene annotation through RNA-seq analysis of the C. briggsae transcriptome, which resulted in the first validated C. briggsae gene set (23,159 genes), among which 7347 genes (33.9% of all genes with introns) have all of their introns confirmed. Most genes (14,812, or 68.3%) have at least one intron validated, compared with only 3.9% in the most recent WormBase release (WS228). Of all introns in the revised gene set (103,083), 61,503 (60.1%) have been confirmed. Additionally, we have identified numerous trans-splicing leaders (SL1 and SL2 variants) in C. briggsae, leading to the first genome-wide annotation of operons in C. briggsae (1105 operons). The majority of the annotated operons (564, or 51.0%) are perfectly conserved in C. elegans, with an additional 345 operons (or 31.2%) somewhat divergent. Additionally, RNA-seq analysis revealed over 10 thousand small-size assembly errors in the current C. briggsae reference genome that can be readily corrected. The revised C. briggsae genome annotation represents a solid platform for comparative genomics analysis and evolutionary studies of Caenorhabditis species.

Oligoarray comparative genomic hybridization-mediated mapping of suppressor mutations generated in a deletion-biased mutagenesis screen

Suppressor screens are an invaluable method for identifying novel genetic interactions between genes in the model organism Caenorhabditis elegans. However, traditionally this approach has suffered from the laborious and protracted process of mapping mutations at the molecular level. Using a mutagen known to generate small deletions, coupled with oligoarray comparative genomic hybridization (aCGH), we have identified mutations in two genes that suppress the lethality associated with a mutation of the essential receptor tyrosine kinase rol-3. First, we find that deletion of the Bicaudal-C ortholog, bcc-1, suppresses rol-3-associated lethality. Second, we identify several duplications that also suppress rol-3-associated lethality. We establish that overexpression of srap-1, a single gene present in these duplications, mediates the suppression. This study demonstrates the suitability of deletion-biased mutagenesis screening in combination with aCGH characterization for the rapid identification of novel suppressor mutations. In addition to detecting small deletions, this approach is suitable for identifying copy number suppressor mutations, a class of suppressor not easily characterized using alternative approaches.

Whole-Genome Profiling of a Novel Mutagenesis Technique Using Proofreading-Deficient DNA Polymerase δ

A novel mutagenesis technique using error-prone DNA polymerase δ (polδ), the disparity mutagenesis model of evolution, has been successfully employed to generate novel microorganism strains with desired traits. However, little else is known about the spectra of mutagenic effects caused by disparity mutagenesis. We evaluated and compared the performance of the polδMKII mutator, which expresses the proofreading-deficient and low-fidelity polδ, in Saccharomyces cerevisiae haploid strain with that of the commonly used chemical mutagen ethyl methanesulfonate (EMS). This mutator strain possesses exogenous mutant polδ supplied from a plasmid, tthereby leaving the genomic one intact. We measured the mutation rate achieved by each mutagen and performed high-throughput next generation sequencing to analyze the genome-wide mutation spectra produced by the 2 mutagenesis methods. The mutation frequency of the mutator was approximately 7 times higher than that of EMS. Our analysis confirmed the strong G/C to A/T transition bias of EMS, whereas we found that the mutator mainly produces transversions, giving rise to more diverse amino acid substitution patterns. Our present study demonstrated that the polδMKII mutator is a useful and efficient method for rapid strain improvement based on in vivo mutagenesis.

Establishing developmental genetics in a self-fertilizing fish (Krytolebias marmoratus)

Kryptolebias marmoratus is a synchronous hermaphroditic vertebrate that utilizes an ovotestis for reproduction. This fish develops externally, is easy to maintain, and has about a 100-day life cycle, making it a desirable developmental genetic model organism. Here, we present a pilot zygotic mutant screen utilizing the common chemical mutagen, N-ethyl-N-nitrosourea (ENU) to establish genetics in this model species. Selection of clonal stocks and optimal conditions for mutagenizing this fish are presented and the types and frequencies of zygotic mutants are documented in comparison to other fish models. Kryptolebias marmoratus is an exemplar model organism that will complement future developmental genetic screens in vertebrates.

Whole genome sequencing for quantifying germline mutation frequency in humans and model species: cautious optimism

Factors affecting the type and frequency of germline mutations in animals are of significant interest from health and toxicology perspectives. However, studies in this field have been limited by the use of markers with low detection power or uncertain relevance to phenotype. Whole genome sequencing (WGS) is now a potential option to directly determine germline mutation type and frequency in family groups at all loci simultaneously. Medical studies have already capitalized on WGS to identify novel mutations in human families for clinical purposes, such as identifying candidate genes contributing to inherited conditions. However, WGS has not yet been used in any studies of vertebrates that aim to quantify changes in germline mutation frequency as a result of environmental factors. WGS is a promising tool for detecting mutation induction, but it is currently limited by several technical challenges. Perhaps the most pressing issue is sequencing error rates that are currently high in comparison to the intergenerational mutation frequency. Different platforms and depths of coverage currently result in a range of 10-10(3) false positives for every true mutation. In addition, the cost of WGS is still relatively high, particularly when comparing mutation frequencies among treatment groups with even moderate sample sizes. Despite these challenges, WGS offers the potential for unprecedented insight into germline mutation processes. Refinement of available tools and emergence of new technologies may be able to provide the improved accuracy and reduced costs necessary to make WGS viable in germline mutation studies in the very near future. To streamline studies, researchers may use multiple family triads per treatment group and sequence a targeted (reduced) portion of each genome with high (20-40 ×) depth of coverage. We are optimistic about the application of WGS for quantifying germline mutations, but caution researchers regarding the resource-intensive nature of the work using existing technology.

Whole-Genome Sequencing of Sordaria macrospora Mutants Identifies Developmental Genes

The study of mutants to elucidate gene functions has a long and successful history; however, to discover causative mutations in mutants that were generated by random mutagenesis often takes years of laboratory work and requires previously generated genetic and/or physical markers, or resources like DNA libraries for complementation. Here, we present an alternative method to identify defective genes in developmental mutants of the filamentous fungus Sordaria macrospora through Illumina/Solexa whole-genome sequencing. We sequenced pooled DNA from progeny of crosses of three mutants and the wild type and were able to pinpoint the causative mutations in the mutant strains through bioinformatics analysis. One mutant is a spore color mutant, and the mutated gene encodes a melanin biosynthesis enzyme. The causative mutation is a G to A change in the first base of an intron, leading to a splice defect. The second mutant carries an allelic mutation in the pro41 gene encoding a protein essential for sexual development. In the mutant, we detected a complex pattern of deletion/rearrangements at the pro41 locus. In the third mutant, a point mutation in the stop codon of a transcription factor-encoding gene leads to the production of immature fruiting bodies. For all mutants, transformation with a wild type-copy of the affected gene restored the wild-type phenotype. Our data demonstrate that whole-genome sequencing of mutant strains is a rapid method to identify developmental genes in an organism that can be genetically crossed and where a reference genome sequence is available, even without prior mapping information.

Allelic ratios and the mutational landscape reveal biologically significant heterozygous SNVs

The issue of heterozygosity continues to be a challenge in the analysis of genome sequences. In this article, we describe the use of allele ratios to distinguish biologically significant single-nucleotide variants from background noise. An application of this approach is the identification of lethal mutations in Caenorhabditis elegans essential genes, which must be maintained by the presence of a wild-type allele on a balancer. The h448 allele of let-504 is rescued by the duplication balancer sDp2. We readily identified the extent of the duplication when the percentage of read support for the lesion was between 70 and 80%. Examination of the EMS-induced changes throughout the genome revealed that these mutations exist in contiguous blocks. During early embryonic division in self-fertilizing C. elegans, alkylated guanines pair with thymines. As a result, EMS-induced changes become fixed as either G→A or C→T changes along the length of the chromosome. Thus, examination of the distribution of EMS-induced changes revealed the mutational and recombinational history of the chromosome, even generations later. We identified the mutational change responsible for the h448 mutation and sequenced PCR products for an additional four alleles, correlating let-504 with the DNA-coding region for an ortholog of a NFκB-activating protein, NKAP. Our results confirm that whole-genome sequencing is an efficient and inexpensive way of identifying nucleotide alterations responsible for lethal phenotypes and can be applied on a large scale to identify the molecular basis of essential genes.

ENU mutagenesis screen to establish motor phenotypes in wild-type mice and modifiers of a pre-existing motor phenotype in tau mutant mice

Modifier screening is a powerful genetic tool. While not widely used in the vertebrate system, we applied these tools to transgenic mouse strains that recapitulate key aspects of Alzheimer's disease (AD), such as tau-expressing mice. These are characterized by a robust pathology including both motor and memory impairment. The phenotype can be modulated by ENU mutagenesis, which results in novel mutant mouse strains and allows identifying the underlying gene/mutation. Here we discuss this strategy in detail. We firstly obtained pedigrees that modify the tau-related motor phenotype, with mapping ongoing. We further obtained transgene-independent motor pedigrees: (i) hyperactive, circling ENU 37 mice with a causal mutation in the Tbx1 gene—the complete knock-out of Tbx1 models DiGeorge Syndrome; (ii) ENU12/301 mice that show sudden jerky movements and tremor constantly; they have a causal mutation in the Kcnq1 gene, modelling aspects of the Romano-Ward and Jervell and Lange-Nielsen syndromes; and (iii) ENU16/069 mice with tremor and hypermetric gait that have a causal mutation in the Mpz (Myelin Protein Zero) gene, modelling Charcot-Marie-Tooth disease type 1 (CMT1B). Together, we provide evidence for a real potential of an ENU mutagenesis to dissect motor functions in wild-type and tau mutant mice.

Mapping mutations in C. elegans

At present, the principal goal of mapping is to establish correspondence between a mutation identified via a change in phenotype and an alteration in the DNA sequence of the genome. Recent advances in molecular biology and bioinformatics have greatly facilitated this procedure, but certain standard methods, such as the three-factor cross, continue to be extremely useful for high-resolution mapping and separation of tightly linked mutations. This chapter provides both general guidelines and specific procedures for the characterization and mapping of newly isolated mutations in C. elegans. Procedures are included for dealing with mutations that cannot be propagated as homozygotes, as well as mutations that can only be scored in specialized genetic backgrounds, for example, suppressor, enhancer, and modifier mutations.

From genes to function: the C. elegans genetic toolbox

This review aims to provide an overview of the technologies which make the nematode Caenorhabditis elegans an attractive genetic model system. We describe transgenesis techniques and forward and reverse genetic approaches to isolate mutants and clone genes. In addition, we discuss the new possibilities offered by genome engineering strategies and next-generation genome analysis tools.

WormBase 2012: more genomes, more data, new website

Since its release in 2000, WormBase (http://www.wormbase.org) has grown from a small resource focusing on a single species and serving a dedicated research community, to one now spanning 15 species essential to the broader biomedical and agricultural research fields. To enhance the rate of curation, we have automated the identification of key data in the scientific literature and use similar methodology for data extraction. To ease access to the data, we are collaborating with journals to link entities in research publications to their report pages at WormBase. To facilitate discovery, we have added new views of the data, integrated large-scale datasets and expanded descriptions of models for human disease. Finally, we have introduced a dramatic overhaul of the WormBase website for public beta testing. Designed to balance complexity and usability, the new site is species-agnostic, highly customizable, and interactive. Casual users and developers alike will be able to leverage the public RESTful application programming interface (API) to generate custom data mining solutions and extensions to the site. We report on the growth of our database and on our work in keeping pace with the growing demand for data, efforts to anticipate the requirements of users and new collaborations with the larger science community.

The early bird catches the worm: new technologies for the Caenorhabditis elegans toolkit

The inherent simplicity of Caenorhabditis elegans and its extensive genetic toolkit make it ideal for studying complex biological processes. Recent developments further increase the usefulness of the worm, including new methods for: altering gene expression, altering physiology using optogenetics, manipulating large numbers of worms, automating laborious processes and processing high-resolution images. These developments both enhance the worm as a model for studying processes such as development and ageing and make it an attractive model in areas such as neurobiology and behaviour.

Damage-induced localized hypermutability

Genome instability continuously presents perils of cancer, genetic disease and death of a cell or an organism. At the same time, it provides for genome plasticity that is essential for development and evolution. We address here the genome instability confined to a small fraction of DNA adjacent to free DNA ends at uncapped telomeres and double-strand breaks. We found that budding yeast cells can tolerate nearly 20 kilobase regions of subtelomeric single-strand DNA that contain multiple UV-damaged nucleotides. During restoration to the double-strand state, multiple mutations are generated by error-prone translesion synthesis. Genome-wide sequencing demonstrated that multiple regions of damage-induced localized hypermutability can be tolerated, which leads to the simultaneous appearance of multiple mutation clusters in the genomes of UV- irradiated cells. High multiplicity and density of mutations suggest that this novel form of genome instability may play significant roles in generating new alleles for evolutionary selection as well as in the incidence of cancer and genetic disease.

Purging deleterious mutations under self fertilization: paradoxical recovery in fitness with increasing mutation rate in Caenorhabditis elegans

BACKGROUND

The accumulation of deleterious mutations can drastically reduce population mean fitness. Self-fertilization is thought to be an effective means of purging deleterious mutations. However, widespread linkage disequilibrium generated and maintained by self-fertilization is predicted to reduce the efficacy of purging when mutations are present at multiple loci.

METHODOLOGY/PRINCIPAL FINDINGS

We tested the ability of self-fertilizing populations to purge deleterious mutations at multiple loci by exposing obligately self-fertilizing populations of Caenorhabditis elegans to a range of elevated mutation rates and found that mutations accumulated, as evidenced by a reduction in mean fitness, in each population. Therefore, purging in obligate selfing populations is overwhelmed by an increase in mutation rate. Surprisingly, we also found that obligate and predominantly self-fertilizing populations exposed to very high mutation rates exhibited consistently greater fitness than those subject to lesser increases in mutation rate, which contradicts the assumption that increases in mutation rate are negatively correlated with fitness. The high levels of genetic linkage inherent in self-fertilization could drive this fitness increase.

CONCLUSIONS

Compensatory mutations can be more frequent under high mutation rates and may alleviate a portion of the fitness lost due to the accumulation of deleterious mutations through epistatic interactions with deleterious mutations. The prolonged maintenance of tightly linked compensatory and deleterious mutations facilitated by self-fertilization may be responsible for the fitness increase as linkage disequilibrium between the compensatory and deleterious mutations preserves their epistatic interaction.

Analysis of multiple ethyl methanesulfonate-mutagenized Caenorhabditis elegans strains by whole-genome sequencing

Whole-genome sequencing (WGS) of organisms displaying a specific mutant phenotype is a powerful approach to identify the genetic determinants of a plethora of biological processes. We have previously validated the feasibility of this approach by identifying a point-mutated locus responsible for a specific phenotype, observed in an ethyl methanesulfonate (EMS)-mutagenized Caenorhabditis elegans strain. Here we describe the genome-wide mutational profile of 17 EMS-mutagenized genomes as assessed with a bioinformatic pipeline, called MAQGene. Surprisingly, we find that while outcrossing mutagenized strains does reduce the total number of mutations, a striking mutational load is still observed even in outcrossed strains. Such genetic complexity has to be taken into account when establishing a causative relationship between genotype and phenotype. Even though unintentional, the 17 sequenced strains described here provide a resource of allelic variants in almost 1000 genes, including 62 premature stop codons, which represent candidate knockout alleles that will be of further use for the C. elegans community to study gene function.