During nervous system development, axons navigate complex environments to reach synaptic targets. Early extending axons must interact with guidance cues in the surrounding tissue, while later extending axons can interact directly with earlier "pioneering" axons, "following" their path. In Caenorhabditis elegans, the AVG neuron pioneers the right axon tract of the ventral nerve cord. We previously found that aex-3, a rab-3 guanine nucleotide exchange factor, is essential for AVG axon navigation in a nid-1 mutant background and that aex-3 might be involved in trafficking of UNC-5, a receptor for the guidance cue UNC-6/netrin. Here, we describe a new gene in this pathway: ccd-5, a putative cdk-5 binding partner. ccd-5 mutants exhibit increased navigation defects of AVG pioneer as well as interneuron and motor neuron follower axons in a nid-1 mutant background. We show that ccd-5 acts in a pathway with cdk-5, aex-3, and unc-5. Navigation defects of follower interneuron and motoneuron axons correlate with AVG pioneer axon defects. This suggests that ccd-5 mostly affects pioneer axon navigation and that follower axon defects are largely a secondary consequence of pioneer navigation defects. To determine the consequences for nervous system function, we assessed various behavioral and movement parameters. ccd-5 single mutants have no significant movement defects, and nid-1 ccd-5 double mutants are less responsive to mechanosensory stimuli compared with nid-1 single mutants. These surprisingly minor defects indicate either a high tolerance for axon guidance defects within the motor circuit and/or an ability to maintain synaptic connections among commonly misguided axons.

It has been estimated that 15-30% of the ∼20,000 genes in C. elegans are essential, yet many of these genes remain to be identified or characterized. With the goal of identifying unknown essential genes, we performed whole genome sequencing on complementation pairs from legacy collections of maternal-effect lethal and sterile mutants. This approach uncovered maternal genes required for embryonic development and genes with putative sperm-specific functions. In total, 58 essential genes were identified on chromosomes III, IV, and V, of which 49 genes are represented by novel alleles in this collection. Of these 49 genes, 19 (40 alleles) were selected for further functional characterization. The terminal phenotypes of embryos were examined, revealing defects in cell division, morphogenesis, and osmotic integrity of the eggshell. Mating assays with wild-type males revealed previously unknown male-expressed genes required for fertilization and embryonic development. The result of this study is a catalogue of mutant alleles in essential genes that will serve as a resource to guide further study toward a more complete understanding of this important model organism. As many genes and developmental pathways in C. elegans are conserved and essential genes are often linked to human disease, uncovering the function of these genes may also provide insight to further our understanding of human biology.

We used whole-genome sequencing (WGS) data from a Caenorhabditis elegans strain homozygous for the reciprocal translocation hT2(I;III) to identify its breakpoints molecularly. The translocation structure is fairly straightforward, with only minor secondary rearrangement in addition to the primary breakpoints. The graphical representation below depicts the two hT2 half-translocations for ease of conceptualization.

We used whole-genome sequencing (WGS) data from a number of balanced lethal strains in Caenorhabditis elegans to show that the crossover suppressor qC1 is an inversion. The rearrangement is complex, with a large primary inversion that contains several other smaller inverted regions. The graphical representation below depicts these various qC1 rearrangements for ease of conceptualization. It is the simplest chromosomal structure compatible with the data currently available, but even then it is worth noting that the complexity of the qC1 chromosome can make the graphical reconstruction difficult to understand, and it may seem a bit like relativity theory or artwork from M.C. Escher (https://moa.byu.edu/m-c-eschers-relativity/).

Improved genome engineering methods that enable automation of large and precise edits are essential for systematic investigations of genome function. We adapted peel-1 negative selection to an optimized Dual-Marker Selection (DMS) cassette protocol for CRISPR-Cas9 genome engineering in Caenorhabditis elegans and observed robust increases in multiple measures of efficiency that were consistent across injectors and four genomic loci. The use of Peel-1-DMS selection killed animals harboring transgenes as extrachromosomal arrays and spared genome-edited integrants, often circumventing the need for visual screening to identify genome-edited animals. To demonstrate the applicability of the approach, we created deletion alleles in the putative proteasomal subunit pbs-1 and the uncharacterized gene K04F10.3 and used machine vision to automatically characterize their phenotypic profiles, revealing homozygous essential and heterozygous behavioral phenotypes. These results provide a robust and scalable approach to rapidly generate and phenotype genome-edited animals without the need for screening or scoring by eye.

John Sulston changed the way we do science, not once, but three times - initially with the complete cell lineage of the nematode Caenorhabditis elegans, next with completion of the genome sequences of the worm and human genomes and finally with his strong and active advocacy for open data sharing. His contributions were widely recognized and in 2002 he received the Nobel Prize in Physiology and Medicine.

Whether generated within a lab setting or isolated from the wild, variant alleles continue to be an important resource for decoding gene function in model organisms such as Caenorhabditis elegans. With advances in massively parallel sequencing, multiple whole-genome sequenced (WGS) strain collections are now available to the research community. The Million Mutation Project (MMP) for instance, analyzed 2007 N2-derived, mutagenized strains. Individually, each strain averages ∼400 single nucleotide variants amounting to ∼80 protein-coding variants. The effects of these variants, however, remain largely uncharacterized and querying the breadth of these strains for phenotypic changes requires a method amenable to rapid and sensitive high-throughput analysis. Here we present a pooled competitive fitness approach to quantitatively phenotype subpopulations of sequenced collections via molecular inversion probes (PhenoMIP). We phenotyped the relative fitness of 217 mutant strains on multiple food sources and classified these into five categories. We also demonstrate on a subset of these strains, that their fitness defects can be genetically mapped. Overall, our results suggest that approximately 80% of MMP mutant strains may have a decreased fitness relative to the lab reference, N2. The costs of generating this form of analysis through WGS methods would be prohibitive while PhenoMIP analysis in this manner is accomplished at less than one-tenth of projected WGS costs. We propose methods for applying PhenoMIP to a broad range of population selection experiments in a cost-efficient manner that would be useful to the community at large.

Germ cell immortality, or transgenerational maintenance of the germ line, could be promoted by mechanisms that could occur in either mitotic or meiotic germ cells. Here we report for the first time that the GSP-2 PP1/Glc7 phosphatase promotes germ cell immortality. Small RNA-induced genome silencing is known to promote germ cell immortality, and we identified a separation-of-function allele of C. elegans gsp-2 that is compromised for germ cell immortality and is also defective for small RNA-induced genome silencing and meiotic but not mitotic chromosome segregation. Previous work has shown that GSP-2 is recruited to meiotic chromosomes by LAB-1, which also promoted germ cell immortality. At the generation of sterility, gsp-2 and lab-1 mutant adults displayed germline degeneration, univalents, histone methylation and histone phosphorylation defects in oocytes, phenotypes that mirror those observed in sterile small RNA-mediated genome silencing mutants. Our data suggest that a meiosis-specific function of GSP-2 ties small RNA-mediated silencing of the epigenome to germ cell immortality. We also show that transgenerational epigenomic silencing at hemizygous genetic elements requires the GSP-2 phosphatase, suggesting a functional link to small RNAs. Given that LAB-1 localizes to the interface between homologous chromosomes during pachytene, we hypothesize that small localized discontinuities at this interface could promote genomic silencing in a manner that depends on small RNAs and the GSP-2 phosphatase.

In a genetic screen to identify genes that promote GLP-1/Notch signaling in Caenorhabditis elegans germline stem cells, we found a single mutation, om40, defining a gene called ego-3. ego-3(om40) causes several defects in the soma and the germline, including paralysis during larval development, sterility, delayed proliferation of germline stem cells, and ectopic germline stem cell proliferation. Whole genome sequencing identified om40 as an allele of hsp-90, previously known as daf-21, which encodes the C. elegans ortholog of the cytosolic form of HSP90. This protein is a molecular chaperone with a central position in the protein homeostasis network, which is responsible for proper folding, structural maintenance, and degradation of proteins. In addition to its essential role in cellular function, HSP90 plays an important role in stem cell maintenance and renewal. Complementation analysis using a deletion allele of hsp-90 confirmed that ego-3 is the same gene. hsp-90(om40) is an I→N conservative missense mutation of a highly conserved residue in the middle domain of HSP-90. RNA interference-mediated knockdown of hsp-90 expression partially phenocopied hsp-90(om40), confirming the loss-of-function nature of hsp-90(om40). Furthermore, reduced HSP-90 activity enhanced the effect of reduced function of both the GLP-1 receptor and the downstream LAG-1 transcription factor. Taken together, our results provide the first experimental evidence of an essential role for HSP90 in Notch signaling in development.

The Caenorhabditis elegans Gene Knockout Consortium is tasked with obtaining null mutations in each of the more than 20,000 open reading frames (ORFs) of this organism. To date, approximately 15,000 ORFs have associated putative null alleles. As there has been substantial success in using CRISPR/Cas9 in C. elegans, this appears to be the most promising technique to complete the task. To enhance the efficiency of using CRISPR/Cas9 to generate gene deletions in C. elegans we provide a web-based interface to access our database of guide RNAs (http://genome.sfu.ca/crispr). When coupled with previously developed selection vectors, optimization for homology arm length, and the use of purified Cas9 protein, we demonstrate a robust and effective protocol for generating deletions for this large-scale project. Debate and speculation in the larger scientific community concerning off-target effects due to non-specific Cas9 cutting has prompted us to investigate through whole genome sequencing the occurrence of single nucleotide variants and indels accompanying targeted deletions. We did not detect any off-site variants above the natural spontaneous mutation rate and therefore conclude that this modified protocol does not generate off-target events to any significant degree in C. elegans. We did, however, observe a number of non-specific alterations at the target site itself following the Cas9-induced double-strand break and offer a protocol for best practice quality control for such events.

The Caenorhabditis elegans Gene Knockout (KO) Consortium is tasked with obtaining null mutations in each of the more than 20,000 open reading frames (ORFs) of this organism. To date, approximately15,000 ORFs have associated putative null alleles. A directed approach using CRISPR/Cas9 methodology is the most promising technique to complete the task. While there has been substantial success in using CRISPR/Cas9 in C. elegans , there has been little emphasis on optimizing the method for generating large insertions/deletions in this organism. To enhance the efficiency of using CRISPR/Cas9 to generate gene knockouts in C. elegans we have developed an online species-specific guide RNA selection tool (<http://genome.sfu.ca/crispr>). When coupled with previously developed selection vectors, optimization for homology arm length, and the use of purified Cas9 protein, we demonstrate a robust, efficient and effective protocol for generating deletions. Debate and speculation in the larger scientific community about off- target effects due to non-specific Cas9 cutting has prompted us to investigate through whole genome sequencing the occurrence of single nucleotide variants and indels accompanying targeted deletions. We did not detect any off-site variants above the natural spontaneous mutation rate and therefore conclude this modified protocol does not generate off-target events to any significant degree in C. elegans.

Genomic silencing can promote germ cell immortality, or transgenerational maintenance of the germ line, via mechanisms that may occur during mitosis or meiosis. Here we report that the gsp-2 PP1/Glc7 phosphatase promotes germ cell immortality. We identified a separation-of-function allele of C. elegans GSP-2 that caused a meiosis-specific chromosome segregation defect and defects in transgenerational small RNA-induced genome silencing. GSP-2 is recruited to meiotic chromosomes by LAB-1, which also promoted germ cell immortality. Sterile gsp-2 and lab-1 mutant adults displayed germline degeneration, univalents and histone phosphorylation defects in oocytes, similar to small RNA genome silencing mutants. Epistasis and RNA analysis suggested that GSP-2 functions downstream of small RNAs. We conclude that a meiosis-specific function of GSP-2/LAB-1 ties small RNA-mediated silencing of the epigenome to germ cell immortality. Given that hemizygous genetic elements can drive transgenerational epigenomic silencing, and given that LAB-1 promotes pairing of homologous chromosomes and localizes to the interface between homologous chromosomes during pachytene, we suggest that discontinuities at this interface could promote nuclear silencing in a manner that depends on GSP-2. Author Summary The germ line of an organism is considered immortal in its capacity to give rise to an unlimited number of future generations. To protect the integrity of the germ line, mechanisms act to suppress the accumulation of transgenerational damage to the genome or epigenome. Loss of germ cell immortality can result from mutations that disrupt the small RNA-mediated silencing pathway that helps to protect the integrity of the epigenome. Here we report for the first time that the C. elegans protein phosphatase GSP-2 that promotes core chromosome biology functions during meiosis is also required for germ cell immortality. Specifically, we identified a partial loss of function allele of gsp-2 that exhibits defects in meiotic chromosome segregation and is also dysfunctional for transgenerational small RNA-mediated genome silencing. Our results are consistent with a known role of Drosophila Protein Phosphatase 1 in heterochromatin silencing, and point to a meiotic phosphatase function that is relevant to germ cell immortality, conceivably related to its roles in chromosome pairing or sister chromatid cohesion.

Our ability to sequence genomes has vastly surpassed our ability to interpret the genetic variation we discover. This presents a major challenge in the clinical setting, where the recent application of whole-exome and whole-genome sequencing has uncovered thousands of genetic variants of uncertain significance. Here, we present a strategy for targeted human gene replacement and phenomic characterization, based on CRISPR-Cas9 genome engineering in the genetic model organism Caenorhabditis elegans, that will facilitate assessment of the functional conservation of human genes and structure-function analysis of disease-associated variants with unprecedented precision. We validate our strategy by demonstrating that direct single-copy replacement of the C. elegans ortholog (daf-18) with the critical human disease-associated gene phosphatase and tensin homolog (PTEN) is sufficient to rescue multiple phenotypic abnormalities caused by complete deletion of daf-18, including complex chemosensory and mechanosensory impairments. In addition, we used our strategy to generate animals harboring a single copy of the known pathogenic lipid phosphatase inactive PTEN variant (PTEN-G129E), and showed that our automated in vivo phenotypic assays could accurately and efficiently classify this missense variant as loss of function. The integrated nature of the human transgenes allows for analysis of both homozygous and heterozygous variants and greatly facilitates high-throughput precision medicine drug screens. By combining genome engineering with rapid and automated phenotypic characterization, our strategy streamlines the identification of novel conserved gene functions in complex sensory and learning phenotypes that can be used as in vivo functional assays to decipher variants of uncertain significance.

Significance Copy-number variants are ubiquitous in nature, yet their immediate functional consequences are obscure. We conducted a spontaneous mutation accumulation experiment at varying sizes in Caenorhabditis elegans , thereby enabling the simultaneous investigation of the mutational input and strength of selection on the evolution of copy-number changes. Whole-genome sequencing reveals the highest genome-wide rate of gene duplication for any species thus far. Our transcriptome analysis further demonstrates that gene duplication frequently results in a greater than two-fold change in transcription. Despite the adaptive role of duplication as the primary source of novel genes, we find duplications and deletions of highly transcribed genes to be more detrimental to fitness and evidence for selection against increase in transcript abundance.

Mitochondrial genomes of metazoans, given their elevated rates of evolution, have served as pivotal markers for phylogeographic studies and recent phylogenetic events. In order to determine the dynamics of spontaneous mitochondrial mutations in small populations in the absence and presence of selection, we evolved mutation accumulation (MA) lines of Caenorhabditis elegans in parallel over 409 consecutive generations at three varying population sizes of N = 1, 10, and 100 hermaphrodites. The N =1 populations should have a minimal influence of natural selection to provide the spontaneous mutation rate and the expected rate of neutral evolution, whereas larger population sizes should experience increasing intensity of selection. New mutations were identified by Illumina paired-end sequencing of 86 mtDNA genomes across 35 experimental lines and compared with published genomes of natural isolates. The spontaneous mitochondrial mutation rate was estimated at 1.05 × 10-7/site/generation. A strong G/C→A/T mutational bias was observed in both the MA lines and the natural isolates. This suggests that the low G + C content at synonymous sites is the product of mutation bias rather than selection as previously proposed. The mitochondrial effective population size per worm generation was estimated to be 62. Although it was previously concluded that heteroplasmy was rare in C. elegans, the vast majority of mutations in this study were heteroplasmic despite an experimental regime exceeding 400 generations. The frequencies of frameshift and nonsynonymous mutations were negatively correlated with population size, which suggests their deleterious effects on fitness and a potent role for selection in their eradication.

Mutants remain a powerful means for dissecting gene function in model organisms such as Caenorhabditis elegans. Massively parallel sequencing has simplified the detection of variants after mutagenesis but determining precisely which change is responsible for phenotypic perturbation remains a key step. Genetic mapping paradigms in C. elegans rely on bulk segregant populations produced by crosses with the problematic Hawaiian wild isolate and an excess of redundant information from whole-genome sequencing (WGS). To increase the repertoire of available mutants and to simplify identification of the causal change, we performed WGS on 173 temperature-sensitive (TS) lethal mutants and devised a novel mapping method. The mapping method uses molecular inversion probes (MIP-MAP) in a targeted sequencing approach to genetic mapping, and replaces the Hawaiian strain with a Million Mutation Project strain with high genomic and phenotypic similarity to the laboratory wild-type strain N2. We validated MIP-MAP on a subset of the TS mutants using a competitive selection approach to produce TS candidate mapping intervals with a mean size < 3 Mb. MIP-MAP successfully uses a non-Hawaiian mapping strain and multiplexed libraries are sequenced at a fraction of the cost of WGS mapping approaches. Our mapping results suggest that the collection of TS mutants contains a diverse library of TS alleles for genes essential to development and reproduction. MIP-MAP is a robust method to genetically map mutations in both viable and essential genes and should be adaptable to other organisms. It may also simplify tracking of individual genotypes within population mixtures.

Forward genetic screens represent powerful, unbiased approaches to uncover novel components in any biological process. Such screens suffer from a major bottleneck, however, namely the cloning of corresponding genes causing the phenotypic variation. Reverse genetic screens have been employed as a way to circumvent this issue, but can often be limited in scope. Here we demonstrate an innovative approach to gene discovery. Using C. elegans as a model system, we used a whole-genome sequenced multi-mutation library, from the Million Mutation Project, together with the Sequence Kernel Association Test (SKAT), to rapidly screen for and identify genes associated with a phenotype of interest, namely defects in dye-filling of ciliated sensory neurons. Such anomalies in dye-filling are often associated with the disruption of cilia, organelles which in humans are implicated in sensory physiology (including vision, smell and hearing), development and disease. Beyond identifying several well characterised dye-filling genes, our approach uncovered three genes not previously linked to ciliated sensory neuron development or function. From these putative novel dye-filling genes, we confirmed the involvement of BGNT-1.1 in ciliated sensory neuron function and morphogenesis. BGNT-1.1 functions at the trans-Golgi network of sheath cells (glia) to influence dye-filling and cilium length, in a cell non-autonomous manner. Notably, BGNT-1.1 is the orthologue of human B3GNT1/B4GAT1, a glycosyltransferase associated with Walker-Warburg syndrome (WWS). WWS is a multigenic disorder characterised by muscular dystrophy as well as brain and eye anomalies. Together, our work unveils an effective and innovative approach to gene discovery, and provides the first evidence that B3GNT1-associated Walker-Warburg syndrome may be considered a ciliopathy.