Whole-genome profiling of mutagenesis in Caenorhabditis elegans

Dopey-dependent regulation of extracellular vesicles maintains neuronal morphology

Mature neurons maintain their distinctive morphology for extended periods in adult life. Compared to developmental neurite outgrowth, axon guidance, and target selection, relatively little is known of mechanisms that maintain the morphology of mature neurons. Loss of function in C. elegans dip-2, a member of the conserved lipid metabolic regulator Dip2 family, results in progressive overgrowth of neurites in adults. We find that dip-2 mutants display specific genetic interactions with sax-2, the C. elegans ortholog of Drosophila Furry and mammalian FRY. Combined loss of dip-2 and sax-2 results in failure to maintain neuronal morphology and elevated release of neuronal extracellular vesicles (EVs). By screening for suppressors of dip-2(0) sax-2(0) double mutant defects, we identified gain-of-function (gf) mutations in the conserved Dopey family protein PAD-1 and its associated phospholipid flippase TAT-5/ATP9A that restore normal neuronal morphology and normal levels of EV release to dip-2(0) sax-2(0) double mutants. Neuron-specific knockdown suggests that PAD-1(gf) can act cell autonomously in neurons. PAD-1(gf) displays increased association with the plasma membrane in oocytes and inhibits EV release in multiple cell types. Our findings uncover a novel functional network of DIP-2, SAX-2, PAD-1, and TAT-5 that maintains neuronal morphology and modulates EV release.

Dopey-dependent regulation of extracellular vesicles maintains neuronal morphology

Mature neurons maintain their distinctive morphology for extended periods in adult life. Compared to developmental neurite outgrowth, axon guidance, and target selection, relatively little is known of mechanisms that maintain mature neuron morphology. Loss of function in C. elegans DIP-2, a member of the conserved lipid metabolic regulator Dip2 family, results in progressive overgrowth of neurites in adults. We find that dip-2 mutants display specific genetic interactions with sax-2, the C. elegans ortholog of Drosophila Furry and mammalian FRY. Combined loss of DIP-2 and SAX-2 results in severe disruption of neuronal morphology maintenance accompanied by increased release of neuronal extracellular vesicles (EVs). By screening for suppressors of dip-2 sax-2 double mutant defects we identified gain-of-function (gf) mutations in the conserved Dopey family protein PAD-1 and its associated phospholipid flippase TAT-5/ATP9A. In dip-2 sax-2 double mutants carrying either pad-1(gf) or tat-5(gf) mutation, EV release is reduced and neuronal morphology across multiple neuron types is restored to largely normal. PAD-1(gf) acts cell autonomously in neurons. The domain containing pad-1(gf) is essential for PAD-1 function, and PAD-1(gf) protein displays increased association with the plasma membrane and inhibits EV release. Our findings uncover a novel functional network of DIP-2, SAX-2, PAD-1, and TAT-5 that maintains morphology of neurons and other types of cells, shedding light on the mechanistic basis of neurological disorders involving human orthologs of these genes.

C. elegans: a prominent platform for modeling and drug screening in neurological disorders

INTRODUCTION

Human neurodevelopmental and neurodegenerative diseases (NDevDs and NDegDs, respectively) encompass a broad spectrum of disorders affecting the nervous system with an increasing incidence. In this context, the nematode C. elegans, has emerged as a benchmark model for biological research, especially in the field of neuroscience.

AREAS COVERED

The authors highlight the numerous advantages of this tiny worm as a model for exploring nervous system pathologies and as a platform for drug discovery. There is a particular focus given to describing the existing models of C. elegans for the study of NDevDs and NDegDs. Specifically, the authors underscore their strong applicability in preclinical drug development. Furthermore, they place particular emphasis on detailing the common techniques employed to explore the nervous system in both healthy and diseased states.

EXPERT OPINION

Drug discovery constitutes a long and expensive process. The incorporation of invertebrate models, such as C. elegans, stands as an exemplary strategy for mitigating costs and expediting timelines. The utilization of C. elegans as a platform to replicate nervous system pathologies and conduct high-throughput automated assays in the initial phases of drug discovery is pivotal for rendering therapeutic options more attainable and cost-effective.

Super-Resolution Microscopy Unveils Synergistic Structural Changes of Organelles Upon Point Mutation

Ethyl methanesulphonate (EMS), is a widely used chemical mutagen that causes high-frequency germline null mutation by inserting an alkyl group into the nucleotide guanine in eukaryotic cells. The effect of EMS on the dynamics of the aneuploid genome, increased cellular instability, and carcinogenicity in relation to benign and malignant tumors are reported, but the molecular level understanding of morphological changes of higher-order chromatin structure has poorly been understood. This is due to a lack of sufficient resolution in conventional microscopic techniques to see small structures below the diffraction limit. Here, using super-resolution radial fluctuation, a largely fragmented, decompaction, and less dense heterochromatin structure upon EMS treatment to HEK 293A cells without any change in nuclear DNA domains is observed. This result suggests an early stage of carcinogenicity happened due to the point mutation. In addition, the distinct structural changes with an elongated morphology of lysosomes are also observed. On the other hand, fragmented and increased heterogeneous populations with an increased cytoplasmic occupancy of mitochondria are observed.

Nuclear factor Y is a pervasive regulator of neuronal gene expression

Nervous system function relies on the establishment of complex gene expression programs that provide neuron-type-specific and core pan-neuronal features. These complementary regulatory paradigms are controlled by terminal selector and parallel-acting transcription factors (TFs), respectively. Here, we identify the nuclear factor Y (NF-Y) TF as a pervasive direct and indirect regulator of both neuron-type-specific and pan-neuronal gene expression. Mapping global NF-Y targets reveals direct binding to the cis-regulatory regions of pan-neuronal genes and terminal selector TFs. We show that NFYA-1 controls pan-neuronal gene expression directly through binding to CCAAT boxes in target gene promoters and indirectly by regulating the expression of terminal selector TFs. Further, we find that NFYA-1 regulation of neuronal gene expression is important for neuronal activity and motor function. Thus, our research sheds light on how global neuronal gene expression programs are buffered through direct and indirect regulatory mechanisms.

A suppressor screen in C. elegans identifies a multiprotein interaction that stabilizes the synaptonemal complex

Successful chromosome segregation into gametes depends on tightly regulated interactions between the parental chromosomes. During meiosis, chromosomes are aligned end-to-end by an interface called the synaptonemal complex, which also regulates exchanges between them. However, despite the functional and ultrastructural conservation of this essential interface, how protein-protein interactions within the synaptonemal complex regulate chromosomal interactions remains poorly understood. Here, we describe a genetic interaction in the C. elegans synaptonemal complex, comprised of short segments of three proteins, SYP-1, SYP-3, and SYP-4. We identified the interaction through a saturated suppressor screen of a mutant that destabilizes the synaptonemal complex. The specificity and tight distribution of suppressors suggest a charge-based interface that promotes interactions between synaptonemal complex subunits and, in turn, allows intimate interactions between chromosomes. Our work highlights the power of genetic studies to illuminate the mechanisms that underlie meiotic chromosome interactions.

Chitosan Nanoparticles Loaded with Capparis cartilaginea Decne Extract: Insights into Characterization and Antigenotoxicity In Vivo

Plant-based foods may enhance the prevention of cancer. The present investigation aimed to assess the antigenotoxic effects of chitosan nanoparticles (CNPs) when loaded with the ethanol extract of C. cartilaginea (CNPs/Cc). Synthesis of CNPs and CNPs/Cc and their characterization were carried out using TEM, EDS, DSC, and Zeta potential. For in vivo experiments, animal groups were treated in the following groups: negative control, ethyl methanesulfonate (EMS) (240 mg/kg), CNPs (350 mg/kg), high and low doses of CNPs/Cc, CNPs plus EMS, high dose of CNPs/Cc plus EMS, and low dose of CNPs/Cc plus EMS. Bone marrow chromosomal aberrations and sperm shape abnormalities were examined. TEM results showed that CNPs and CNPs/Cc are spherical particles. CNPs' physical stability was observed to be lower than that of CNPs/Cc due to the presence of more positive charges on CNPs/Cc. EMS significantly enhanced chromosomal abnormalities and sperm shape abnormalities. CNPs showed powerful antigenotoxic properties. For the first time, it could be concluded that loading chitosan nanoparticles with C. cartilaginea extract significantly promotes its protective properties.

A suppressor screen in C. elegans identifies a multi-protein interaction interface that stabilizes the synaptonemal complex

Successful chromosome segregation into gametes depends on tightly-regulated interactions between the parental chromosomes. During meiosis, chromosomes are aligned end-to-end by an interface called the synaptonemal complex, which also regulates exchanges between them. However, despite the functional and ultrastructural conservation of this essential interface, how protein-protein interactions within the synaptonemal complex regulate chromosomal interactions remains poorly understood. Here we describe a novel interaction interface in the C. elegans synaptonemal complex, comprised of short segments of three proteins, SYP-1, SYP-3 and SYP-4. We identified the interface through a saturated suppressor screen of a mutant that destabilizes the synaptonemal complex. The specificity and tight distribution of suppressors point to a charge-based interface that promotes interactions between synaptonemal complex subunits and, in turn, allows intimate interactions between chromosomes. Our work highlights the power of genetic studies to illuminate the mechanisms that underly meiotic chromosome interactions.

Evaluation of mutagenic susceptibility of different stages in germ cell development of Caenorhabditis elegans using whole genome sequencing

In contrast to somatic mutations, mutations in germ cells affect every cell of any organism derived from the germ cell and therefore are related to numerous genetic diseases. However, there is no suitable assay to evaluate the mutagenic sensitivities of both male and female germ cells. The main type of Caenorhabditis elegans (C. elegans) is hermaphroditic, where spermatogenesis and oogenesis occur chronologically at specific stages, allowing induction of mutations in either sperm or eggs exclusively. In this study, we used the alkylating agent ethyl methanesulfonate and N-ethyl-N-nitrosourea to induce germline mutations in C. elegans at different developmental stages and analyzed mutation frequency and mutational spectrum from data gathered using next-generation sequencing (NGS) technology. Our results revealed low spontaneous mutation rates of C. elegans, along with distinct mutagenic effects elicited by the two mutagens. Our data show that the parental worms treated during germ cell mitosis, spermatogenesis, and oogenesis resulted in different mutation frequencies in their offspring, and female germ cells could be very susceptible to mutagen exposure during oogenesis. In summary, our study indicates that the use of C. elegans and its specific chronological hermaphroditism would be a promising way to explore the sensitivities of both male and female germ cells to mutagens.

The effects of mutations on gene expression and alternative splicing

Understanding the relationship between mutations and their genomic and phenotypic consequences has been a longstanding goal of evolutionary biology. However, few studies have investigated the impact of mutations on gene expression and alternative splicing on the genome-wide scale. In this study, we aim to bridge this knowledge gap by utilizing whole-genome sequencing data and RNA sequencing data from 16 obligately parthenogenetic Daphnia mutant lines to investigate the effects of ethyl methanesulfonate-induced mutations on gene expression and alternative splicing. Using rigorous analyses of mutations, expression changes and alternative splicing, we show that trans-effects are the major contributor to the variance in gene expression and alternative splicing between the wild-type and mutant lines, whereas cis mutations only affected a limited number of genes and do not always alter gene expression. Moreover, we show that there is a significant association between differentially expressed genes and exonic mutations, indicating that exonic mutations are an important driver of altered gene expression.

GALA: a computational framework for de novo chromosome-by-chromosome assembly with long reads

High-quality genome assembly has wide applications in genetics and medical studies. However, it is still very challenging to achieve gap-free chromosome-scale assemblies using current workflows for long-read platforms. Here we report on GALA (Gap-free long-read Assembly tool), a computational framework for chromosome-based sequencing data separation and de novo assembly implemented through a multi-layer graph that identifies discordances within preliminary assemblies and partitions the data into chromosome-scale scaffolding groups. The subsequent independent assembly of each scaffolding group generates a gap-free assembly likely free from the mis-assembly errors which usually hamper existing workflows. This flexible framework also allows us to integrate data from various technologies, such as Hi-C, genetic maps, and even motif analyses to generate gap-free chromosome-scale assemblies. As a proof of principle we de novo assemble the C. elegans genome using combined PacBio and Nanopore sequencing data and a rice cultivar genome using Nanopore sequencing data from publicly available datasets. We also demonstrate the proposed method's applicability with a gap-free assembly of the human genome using PacBio high-fidelity (HiFi) long reads. Thus, our method enables straightforward assembly of genomes with multiple data sources and overcomes barriers that at present restrict the application of de novo genome assembly technology.

Random Mutagenesis as a Promising Tool for Microalgal Strain Improvement towards Industrial Production

Microalgae have become a promising novel and sustainable feedstock for meeting the rising demand for food and feed. However, microalgae-based products are currently hindered by high production costs. One major reason for this is that commonly cultivated wildtype strains do not possess the robustness and productivity required for successful industrial production. Several strain improvement technologies have been developed towards creating more stress tolerant and productive strains. While classical methods of forward genetics have been extensively used to determine gene function of randomly generated mutants, reverse genetics has been explored to generate specific mutations and target phenotypes. Site-directed mutagenesis can be accomplished by employing different gene editing tools, which enable the generation of tailor-made genotypes. Nevertheless, strategies promoting the selection of randomly generated mutants avoid the introduction of foreign genetic material. In this paper, we review different microalgal strain improvement approaches and their applications, with a primary focus on random mutagenesis. Current challenges hampering strain improvement, selection, and commercialization will be discussed. The combination of these approaches with high-throughput technologies, such as fluorescence-activated cell sorting, as tools to select the most promising mutants, will also be discussed.

Systematic prediction of EMS-induced mutations in a sorghum mutant population

The precise detection of causal DNA mutations (deoxyribonucleic acid) is very crucial for forward genetic studies. Several sources of errors contribute to false-positive detections by current variant-calling algorithms, which impact associating phenotypes with genotypes. To improve the accuracy of mutation detection, we implemented a binning method for the accurate detection of likely ethyl methanesulfonate (EMS)-induced mutations in a sequenced mutant population. We also implemented a clustering algorithm for detecting likely false negatives with high accuracy. Sorghum bicolor is a very valuable crop species with tremendous potential for uncovering novel gene functions associated with highly desirable agronomical traits. We demonstrate the precision of the described approach in the detection of likely EMS-induced mutations from the publicly available low-cost sequencing of the M₃ generation from 600 sorghum BTx623 mutants. The approach detected 3,274,606 single nucleotide polymorphisms (SNPs), of which 96% (3,141,908) were G/C to A/T DNA substitutions, as expected by EMS-mutagenesis mode of action. We demonstrated the general applicability of the described method and showed a high concordance, 94% (3,074,759) SNPs overlap between SAMtools-based and GATK-based variant-calling algorithms. Our clustering algorithm uncovered evidence for an additional 223,048 likely false-negative shared EMS-induced mutations. The final 3,497,654 SNPs represent an 87% increase in SNPs detected from the previous analysis of the mutant population, with an average of one SNP per 125 kb in the sorghum genome. Annotation of the final SNPs revealed 10,263 high-impact and 136,639 moderate-impact SNPs, including 7217 stop-gained mutations, which averages 12 stop-gained mutations per mutant, and four high- or medium-impact SNPs per sorghum gene. We have implemented a public search database for this new genetic resource of 30,285 distinct sorghum genes containing medium- or high-impact EMS-induced mutations. Seedstock for a select 486 of the 600 described mutants are publicly available in the Germplasm Resources Information Network (GRIN) database.

Genome-wide detection of ultralow-frequency substitution mutations in cultures of mouse lymphoma L5178Y cells and Caenorhabditis elegans worms by PacBio sequencing

Many conventional genetic toxicology assays require specialized cell cultures or animals and can only detect mutations that inactivate the function of a reporter gene. These limitations make such assays incompatible with many toxicological models but could be overcome by the development of techniques capable of directly detecting genome-wide somatic mutations through DNA sequencing. PacBio sequencing can generate almost error-free consensus reads by repeatedly inspecting both DNA strands from circularized molecules (a method known as PacBio HiFi). In this study, we show that PacBio HiFi can detect genome-wide ultralow-frequency substitution mutations in cultures of mouse lymphoma L5178Y cells and Caenorhabditis elegans worms. The mutation frequencies (MFs) of unexposed samples in both models were ~1 × 10^-7 mutations per base pair. Compared to these controls, PacBio HiFi detected MF increases of 23-fold in cultures of L5178Y cells exposed to 5 mM ethyl methanosulfonate (EMS) for 4 h, and 5-, 12-, and 29-fold in cultures of C. elegans worms exposed to 12.5, 25, and 50 mM EMS for 4 h, respectively. In both models, the mutation spectra of controls were diverse, while those derived from EMS-exposed samples were dominated by C:G → T:A transitions. To validate these results, clone sequencing analyses were performed on the same cultures of L5178Y cells. The results obtained by clone sequencing and PacBio HiFi were almost identical. Our results suggest that PacBio sequencing could be used for the detection, quantitation, and characterization of mutations in any DNA-containing sample, including those that are not compatible with conventional mutation detection approaches.

Identification of Essential Genes in Caenorhabditis elegans with Lethal Mutations Maintained by Genetic Balancers

Genetic balancer systems, which allow effective capture and maintenance of lethal mutations stably, play an important role in identifying essential genes. Whole-genome sequencing (WGS) followed by bioinformatics analysis, combined with genetic mapping data analysis, allows for an efficient and economical means of identifying genomic mutations in essential genes. Using this approach, we successfully identified 104 essential genes on ChrI, ChrIII, and ChrV in C. elegans. In this report, we described a protocol that sequences the genome of prebalanced Caenorhabditis elegans (C. elegans) strains to carry lethal mutations and identifies candidate causal mutations and candidate essential genes using a robust bioinformatics procedure.

Whole genome sequencing facilitates intragenic variant interpretation following modifier screening in C. elegans

Background

Intragenic modifiers (in-phase, second-site variants) are known to have dramatic effects on clinical outcomes, affecting disease attributes such as severity or age of onset. However, despite their clinical importance, the focus of many genetic screens in model systems is on the discovery of extragenic variants, with many labs still relying upon more traditional methods to identify modifiers. However, traditional methods such as PCR and Sanger sequencing can be time-intensive and do not permit a thorough understanding of the intragenic modifier effects in the context of non-isogenic genomic backgrounds.

Results

Here, we apply high throughput approaches to identify and understand intragenic modifiers using Caenorhabditis elegans. Specifically, we applied whole genome sequencing (WGS) to a mutagen-induced forward genetic screen to identify intragenic suppressors of a temperature-sensitive zyg-1(it25) allele in C. elegans. ZYG-1 is a polo kinase that is important for centriole function and cell divisions, and mutations that truncate its human orthologue, PLK4, have been associated with microcephaly. Combining WGS and CRISPR/Cas9, we rapidly identify intragenic modifiers, show that these variants are distributed non-randomly throughout zyg-1 and that genomic context plays an important role on phenotypic outcomes.

Conclusions

Ultimately, our work shows that WGS facilitates high-throughput identification of intragenic modifiers in clinically relevant genes by reducing hands-on research time and overall costs and by allowing thorough understanding of the intragenic phenotypic effects in the context of different genetic backgrounds.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08142-8.

The genome-wide rate and spectrum of EMS-induced heritable mutations in the microcrustacean Daphnia: on the prospect of forward genetics

Forward genetic screening using the alkylating mutagen ethyl methanesulfonate (EMS) is an effective method for identifying phenotypic mutants of interest, which can be further genetically dissected to pinpoint the causal genetic mutations. An accurate estimate of the rate of EMS-induced heritable mutations is fundamental for determining the mutant sample size of a screening experiment that aims to saturate all the genes in a genome with mutations. This study examines the genome-wide EMS-induced heritable base-substitutions in three species of the freshwater microcrustacean Daphnia to help guide screening experiments. Our results show that the 10 mM EMS treatment induces base substitutions at an average rate of 1.17 × 10^-6/site/generation across the three species, whereas a significantly higher average mutation rate of 1.75 × 10^-6 occurs at 25 mM. The mutation spectrum of EMS-induced base substitutions at both concentration is dominated by G:C to A:T transitions. Furthermore, we find that female Daphnia exposed to EMS (F₀ individuals) can asexually produce unique mutant offspring (F₁) for at least 3 consecutive broods, suggestive of multiple broods as F₁ mutants. Lastly, we estimate that about 750 F₁s are needed for all genes in the Daphnia genome to be mutated at least once with a 95% probability. We also recommend 4-5 F₂s should be collected from each F₁ mutant through sibling crossing so that all induced mutations could appear in the homozygous state in the F₂ population at 70-80% probability.

Identification of essential genes in Caenorhabditis elegans through whole genome sequencing of legacy mutant collections

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 apparent sperm-specific functions. In total, 58 putative essential genes were identified on chromosomes III–V, of which 52 genes are represented by novel alleles in this collection. Of these 52 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 catalog 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.

Harnessing model organism genomics to underpin the machine learning-based prediction of essential genes in eukaryotes - Biotechnological implications

The availability of high-quality genomes and advances in functional genomics have enabled large-scale studies of essential genes in model eukaryotes, including the 'elegant worm' (Caenorhabditis elegans; Nematoda) and the 'vinegar fly' (Drosophila melanogaster; Arthropoda). However, this is not the case for other, much less-studied organisms, such as socioeconomically important parasites, for which functional genomic platforms usually do not exist. Thus, there is a need to develop innovative techniques or approaches for the prediction, identification and investigation of essential genes. A key approach that could enable the prediction of such genes is machine learning (ML). Here, we undertake an historical review of experimental and computational approaches employed for the characterisation of essential genes in eukaryotes, with a particular focus on model ecdysozoans (C. elegans and D. melanogaster), and discuss the possible applicability of ML-approaches to organisms such as socioeconomically important parasites. We highlight some recent results showing that high-performance ML, combined with feature engineering, allows a reliable prediction of essential genes from extensive, publicly available 'omic data sets, with major potential to prioritise such genes (with statistical confidence) for subsequent functional genomic validation. These findings could 'open the door' to fundamental and applied research areas. Evidence of some commonality in the essential gene-complement between these two organisms indicates that an ML-engineering approach could find broader applicability to ecdysozoans such as parasitic nematodes or arthropods, provided that suitably large and informative data sets become/are available for proper feature engineering, and for the robust training and validation of algorithms. This area warrants detailed exploration to, for example, facilitate the identification and characterisation of essential molecules as novel targets for drugs and vaccines against parasitic diseases. This focus is particularly important, given the substantial impact that such diseases have worldwide, and the current challenges associated with their prevention and control and with drug resistance in parasite populations.

Preservation of lagging strand integrity at sites of stalled replication by Pol α-primase and 9-1-1 complex

During genome duplication, the replication fork encounters a plethora of obstacles in the form of damaged bases, DNA-cross-linked proteins, and secondary structures. How cells protect DNA integrity at sites of stalled replication is currently unknown. Here, by engineering "primase deserts" into the Caenorhabditis elegans genome close to replication-impeding G-quadruplexes, we show that de novo DNA synthesis downstream of the blocked fork suppresses DNA loss. We next identify the pol α-primase complex to limit deletion mutagenesis, a conclusion substantiated by whole-genome analysis of animals carrying mutated POLA2/DIV-1. We subsequently identify a new role for the 9-1-1 checkpoint clamp in protecting Okazaki fragments from resection by EXO1. Together, our results provide a mechanistic model for controlling the fate of replication intermediates at sites of stalled replication.

Allele-specific suppression in C. elegans reveals details of EMS mutagenesis and a possible moonlighting interaction between the vesicular acetylcholine transporter and ERD2 receptors

A missense mutant, unc-17(e245), which affects the Caenorhabditis elegans vesicular acetylcholine transporter UNC-17, has a severe uncoordinated phenotype, allowing efficient selection of dominant suppressors that revert this phenotype to wild-type. Such selections permitted isolation of numerous suppressors after EMS (ethyl methanesulfonate) mutagenesis, leading to demonstration of delays in mutation fixation after initial EMS treatment, as has been shown in T4 bacteriophage but not previously in eukaryotes. Three strong dominant extragenic suppressor loci have been defined, all of which act specifically on allele e245, which causes a G347R mutation in UNC-17. Two of the suppressors (sup-1 and sup-8/snb-1) have previously been shown to encode synaptic proteins able to interact directly with UNC-17. We found that the remaining suppressor, sup-2, corresponds to a mutation in erd-2.1, which encodes an endoplasmic reticulum retention protein; sup-2 causes a V186E missense mutation in transmembrane helix 7 of ERD-2.1. The same missense change introduced into the redundant paralogous gene erd-2.2 also suppressed unc-17(e245). Suppression presumably occurred by compensatory charge interactions between transmembrane helices of UNC-17 and ERD-2.1 or ERD-2.2, as previously proposed in work on suppression by SUP-1(G84E) or SUP-8(I97D)/synaptobrevin. erd-2.1(V186E) homozygotes were fully viable, but erd-2.1(V186E); erd-2.2(RNAi) exhibited synthetic lethality (like erd-2.1(RNAi); erd-2.2(RNAi)), indicating that the missense change in ERD-2.1 impairs its normal function in the secretory pathway but may allow it to adopt a novel moonlighting function as an unc-17 suppressor.

Efficient collection of a large number of mutations by mutagenesis of DNA damage response defective animals

With the development of massive parallel sequencing technology, it has become easier to establish new model organisms that are ideally suited to the specific biological phenomena of interest. Considering the history of research using classical model organisms, we believe that the efficient construction and sharing of gene mutation libraries will facilitate the progress of studies using these new model organisms. Using C. elegans, we applied the TMP/UV mutagenesis method to animals lacking function in the DNA damage response genes atm-1 and xpc-1. This method produces genetic mutations three times more efficiently than mutagenesis of wild-type animals. Furthermore, we confirmed that the use of next-generation sequencing and the elimination of false positives through machine learning could automate the process of mutation identification with an accuracy of over 95%. Eventually, we sequenced the whole genomes of 488 strains and isolated 981 novel mutations generated by the present method; these strains have been made available to anyone who wants to use them. Since the targeted DNA damage response genes are well conserved and the mutagens used in this study are also effective in a variety of species, we believe that our method is generally applicable to a wide range of animal species.

DNA repair, recombination, and damage signaling

DNA must be accurately copied and propagated from one cell division to the next, and from one generation to the next. To ensure the faithful transmission of the genome, a plethora of distinct as well as overlapping DNA repair and recombination pathways have evolved. These pathways repair a large variety of lesions, including alterations to single nucleotides and DNA single and double-strand breaks, that are generated as a consequence of normal cellular function or by external DNA damaging agents. In addition to the proteins that mediate DNA repair, checkpoint pathways have also evolved to monitor the genome and coordinate the action of various repair pathways. Checkpoints facilitate repair by mediating a transient cell cycle arrest, or through initiation of cell suicide if DNA damage has overwhelmed repair capacity. In this chapter, we describe the attributes of Caenorhabditis elegans that facilitate analyses of DNA repair, recombination, and checkpoint signaling in the context of a whole animal. We review the current knowledge of C. elegans DNA repair, recombination, and DNA damage response pathways, and their role during development, growth, and in the germ line. We also discuss how the analysis of mutational signatures in C. elegans is helping to inform cancer mutational signatures in humans.

Bedrock radioactivity influences the rate and spectrum of mutation

All organisms on Earth are exposed to low doses of natural radioactivity but some habitats are more radioactive than others. Yet, documenting the influence of natural radioactivity on the evolution of biodiversity is challenging. Here, we addressed whether organisms living in naturally more radioactive habitats accumulate more mutations across generations using 14 species of waterlice living in subterranean habitats with contrasted levels of radioactivity. We found that the mitochondrial and nuclear mutation rates across a waterlouse species’ genome increased on average by 60% and 30%, respectively, when radioactivity increased by a factor of three. We also found a positive correlation between the level of radioactivity and the probability of G to T (and complementary C to A) mutations, a hallmark of oxidative stress. We conclude that even low doses of natural bedrock radioactivity influence the mutation rate possibly through the accumulation of oxidative damage, in particular in the mitochondrial genome.

Terminal web and vesicle trafficking proteins mediate nematode single-cell tubulogenesis

Single-celled tubules represent a complicated structure that forms during development, requiring extension of a narrow cytoplasm surrounding a lumen exerting osmotic pressure that can burst the luminal membrane. Genetic studies on the excretory canal cell of Caenorhabditis elegans have revealed many proteins that regulate the cytoskeleton, vesicular transport, and physiology of the narrow canals. Here, we show that βH-spectrin regulates the placement of intermediate filament proteins forming a terminal web around the lumen, and that the terminal web in turn retains a highly conserved protein (EXC-9/CRIP1) that regulates apical endosomal trafficking. EXC-1/IRG, the binding partner of EXC-9, is also localized to the apical membrane and affects apical actin placement and RAB-8-mediated vesicular transport. The results suggest that an intermediate filament protein acts in a novel pathway to direct the traffic of vesicles to locations of lengthening apical surface during single-celled tubule development.

Studies of an insecticidal inhibitor of acetyl-CoA carboxylase in the nematode C. elegans

We have studied the mode of action of the insecticide spirotetramat in the nematode Caenorhabditis elegans. A combination of symptomology, forward genetics and genome editing show that spirotetramat acts on acetyl-CoA carboxylase (ACC) in C. elegans, as it does in insects. We found C. elegans embryos exposed to spirotetramat show a cell division defect which closely resembles the phenotype of loss-of-function mutations in the gene pod-2, which encodes ACC. We then identified two mutations in the carboxyl transferase domain of pod-2 (ACC) which confer resistance and were confirmed using CRISPR/Cas9. One of these mutations substitutes an invertebrate-specific amino acid with one ubiquitous in other taxa; this residue may, therefore, be a determinant of the selectivity of spirotetramat for invertebrates. Such a mutation may also be the target of selection for resistance in the field. Our study is a further demonstration of the utility of C. elegans in studying bioactive chemicals.

A Rapid, SuperSelective Method for Detection of Single Nucleotide Variants in Caenorhabditis elegans

With the widespread use of single nucleotide variants generated through mutagenesis screens and genome editing technologies, there is pressing need for an efficient and low-cost strategy to genotype single nucleotide substitutions. We have developed a rapid and inexpensive method for detection of point mutants through optimization of SuperSelective (SS) primers for end-point PCR in Caenorhabditis elegans Each SS primer consists of a 5' "anchor" that hybridizes to the template, followed by a noncomplementary "bridge," and a "foot" corresponding to the target allele. The foot sequence is short, such that a single mismatch at the terminal 3' nucleotide destabilizes primer binding and prevents extension, enabling discrimination of different alleles. We explored how length and sequence composition of each SS primer segment affected selectivity and efficiency in various genetic contexts in order to develop simple rules for primer design that allow for differentiation between alleles over a broad range of annealing temperatures. Manipulating bridge length affected amplification efficiency, while modifying the foot sequence altered discriminatory power. Changing the anchor position enabled SS primers to be used for genotyping in regions with sequences that are challenging for standard primer design. After defining primer design parameters, we demonstrated the utility of SS primers for genotyping crude C. elegans lysates, suggesting that this approach could also be used for SNP mapping and screening of CRISPR mutants. Further, since SS primers reliably detect point mutations, this method has potential for broad application in all genetic systems.

Harnessing the power of genetics: fast forward genetics in Caenorhabditis elegans

Forward genetics is a powerful tool to unravel molecular mechanisms of diverse biological processes. The success of genetic screens primarily relies on the ease of genetic manipulation of an organism and the availability of a plethora of genetic tools. The roundworm Caenorhabditis elegans has been one of the favorite models for genetic studies due to its hermaphroditic lifestyle, ease of maintenance, and availability of various genetic manipulation tools. The strength of C. elegans genetics is highlighted by the leading role of this organism in the discovery of several conserved biological processes. In this review, the principles and strategies for forward genetics in C. elegans are discussed. Further, the recent advancements that have drastically accelerated the otherwise time-consuming process of mutation identification, making forward genetic screens a method of choice for understanding biological functions, are discussed. The emphasis of the review has been on providing practical and conceptual pointers for designing genetic screens that will identify mutations, specifically disrupting the biological processes of interest.

Analysis of mutational signatures in C. elegans: Implications for cancer genome analysis

Genome integrity is constantly challenged by exogenous and endogenous insults, and mutations are associated with inherited disease and cancer. Here we summarize recent studies that utilized C. elegans whole genome next generation sequencing to experimentally determine mutational signatures associated with mutagen exposure, DNA repair deficiency or a combination of both and discuss the implications of these results for the understanding of cancer genome evolution. The experimental analysis of wild-type and DNA repair deficient nematodes propagated under unchallenged conditions over many generations revealed increased mutagenesis in approximately half of all DNA repair deficient strains, its rate, except for DNA mismatch repair, only being moderately increased. The exposure of wild-type and DNA repair defective strains to selected genotoxins, including UV-B and ionizing radiation, alkylating compounds, aristolochic acid, aflatoxin-B1, and cisplatin enabled the systematic analysis of the relative contributions of redundant repair modalities that mend DNA damage. Combining genotoxin exposure with DNA repair deficiency can manifest as altered mutation rates and/or as a change in mutational profiles, and reveals how different DNA alterations induced by one genotoxin are repaired by separate DNA repair pathways, often in a highly redundant way. Cancer genomes provide a snapshot of all mutational events that happened prior to cancer detection and sequencing, necessitating computational models to deduce mutational signatures using mathematical best fit approaches. While computationally deducing signatures from cancer genomes has been tremendously successful in associating some signatures to known mutagenic causes, many inferred signatures lack a clear link to a known mutagenic process. Moreover, analytical signatures might not reflect any distinct mutagenic processes. Nonetheless, combined effects of mutagen exposure and DNA damage-repair deficiency are also present in cancer genomes, but cannot be as easily detected owing to the unknown histories of genotoxic exposures and because biallelic in contrast to monoallelic DNA repair deficiency is rare. The impact of damage-repair interactions also manifests through selective pressure for DNA repair gene inactivation during cancer evolution. Using these considerations, we discuss a theoretical framework that explains why minute mutagenic changes, possibly too small to manifest as change in a signature, can have major effects in oncogenesis. Overall, the experimental analysis of mutational processes underscores that the interpretation of mutational signatures requires considering both the primary DNA lesion and repair status and imply that mutational signatures derived from cancer genomes may be more variable than currently anticipated.

NHJ-1 Is Required for Canonical Nonhomologous End Joining in Caenorhabditis elegans

DNA double-strand breaks (DSBs) are a particularly lethal form of DNA damage that must be repaired to restore genomic integrity. Canonical nonhomologous end joining (NHEJ), is a widely conserved pathway that detects and directly ligates the broken ends to repair the DSB. These events globally require the two proteins that form the Ku ring complex, Ku70 and Ku80, and the terminal ligase LIG4. While the NHEJ pathway in vertebrates is elaborated by more than a dozen factors of varying conservation and is similarly complex in other eukaryotes, the entire known NHEJ toolkit in Caenorhabditis elegans consists only of the core components CKU-70, CKU-80, and LIG-4 Here, we report the discovery of the first accessory NHEJ factor in C. elegans Our analysis of the DNA damage response in young larvae revealed that the canonical wild-type N2 strain consisted of two lines that exhibited a differential phenotypic response to ionizing radiation (IR). Following the mapping of the causative locus to a candidate on chromosome V and clustered regularly interspaced short palindromic repeats-Cas9 mutagenesis, we show that disruption of the nhj-1 sequence induces IR sensitivity in the N2 line that previously exhibited IR resistance. Using genetic and cytological analyses, we demonstrate that nhj-1 functions in the NHEJ pathway to repair DSBs. Double mutants of nhj-1 and lig-4 or cku-80 do not exhibit additive IR sensitivity, and the post-IR somatic and fertility phenotypes of nhj-1 mimic those of the other NHEJ factors. Furthermore, in com-1 mutants that permit repair of meiotic DSBs by NHEJ instead of restricting their repair to the homologous recombination pathway, loss of nhj-1 mimics the consequences of loss of lig-4 Diakinesis-stage nuclei in nhj-1; com-1 and nhj-1; lig-4 mutant germlines exhibit increased numbers of DAPI-staining bodies, consistent with increased chromosome fragmentation in the absence of NHEJ-mediated meiotic DSB repair. Finally, we show that NHJ-1 and LIG-4 localize to somatic nuclei in larvae, but are excluded from the germline progenitor cells, consistent with NHEJ being the dominant DNA repair pathway in the soma. nhj-1 shares no sequence homology with other known eukaryotic NHEJ factors and is taxonomically restricted to the Rhabditid family, underscoring the evolutionary plasticity of even highly conserved pathways.

Whole-genome sequencing of Puccinia striiformis f. sp. tritici mutant isolates identifies avirulence gene candidates

Background

The stripe rust pathogen, Puccinia striiformis f. sp. tritici (Pst), threats world wheat production. Resistance to Pst is often overcome by pathogen virulence changes, but the mechanisms of variation are not clearly understood. To determine the role of mutation in Pst virulence changes, in previous studies 30 mutant isolates were developed from a least virulent isolate using ethyl methanesulfonate (EMS) mutagenesis and phenotyped for virulence changes. The progenitor isolate was sequenced, assembled and annotated for establishing a high-quality reference genome. In the present study, the 30 mutant isolates were sequenced and compared to the wide-type isolate to determine the genomic variation and identify candidates for avirulence (Avr) genes.

Results

The sequence reads of the 30 mutant isolates were mapped to the wild-type reference genome to identify genomic changes. After selecting EMS preferred mutations, 264,630 and 118,913 single nucleotide polymorphism (SNP) sites and 89,078 and 72,513 Indels (Insertion/deletion) were detected among the 30 mutant isolates compared to the primary scaffolds and haplotigs of the wild-type isolate, respectively. Deleterious variants including SNPs and Indels occurred in 1866 genes. Genome wide association analysis identified 754 genes associated with avirulence phenotypes. A total of 62 genes were found significantly associated to 16 avirulence genes after selection through six criteria for putative effectors and degree of association, including 48 genes encoding secreted proteins (SPs) and 14 non-SP genes but with high levels of association (P ≤ 0.001) to avirulence phenotypes. Eight of the SP genes were identified as avirulence-associated effectors with high-confidence as they met five or six criteria used to determine effectors.

Conclusions

Genome sequence comparison of the mutant isolates with the progenitor isolate unraveled a large number of mutation sites along the genome and identified high-confidence effector genes as candidates for avirulence genes in Pst. Since the avirulence gene candidates were identified from associated SNPs and Indels caused by artificial mutagenesis, these avirulence gene candidates are valuable resources for elucidating the mechanisms of the pathogen pathogenicity, and will be studied to determine their functions in the interactions between the wheat host and the Pst pathogen.

Mutational Analysis of N-Ethyl-N-Nitrosourea (ENU) in the Fission Yeast Schizosaccharomyces pombe

Forward genetics in model organisms has boosted our knowledge of the genetic bases of development, aging, and human diseases. In this experimental pipeline, it is crucial to start by inducing a large number of random mutations in the genome of the model organism to search for phenotypes of interest. Many chemical mutagens are used to this end because most of them display particular reactivity properties and act differently over DNA. Here we report the use of N-ethyl-N-nitrosourea (ENU) as a mutagen in the fission yeast Schizosaccharomyces pombe As opposed to many other alkylating agents, ENU only induces an S N 1-type reaction with a low s constant (s = 0.26), attacking preferentially O2 and O4 in thymine and O6 deoxyguanosine, leading to base substitutions rather than indels, which are extremely rare in its resulting mutagenic repertoire. Using ENU, we gathered a collection of 13 temperature-sensitive mutants and 80 auxotrophic mutants including two deleterious alleles of the human ortholog ATIC. Defective alleles of this gene cause AICA-ribosiduria, a severe genetic disease. In this screen, we also identified 13 aminoglycoside-resistance inactivating mutations in APH genes. Mutations reported here may be of interest for metabolism related diseases and antibiotic resistance research fields.

Tyrosine aminotransferase is involved in the oxidative stress response by metabolizing meta-tyrosine in Caenorhabditis elegans

Under oxidative stress conditions, hydroxyl radicals can oxidize the phenyl ring of phenylalanine, producing the abnormal tyrosine isomer meta-tyrosine (m-tyrosine). m-Tyrosine levels are commonly used as a biomarker of oxidative stress, and its accumulation has recently been reported to adversely affect cells, suggesting a direct role for m-tyrosine in oxidative stress effects. We found that the Caenorhabditis elegans ortholog of tyrosine aminotransferase (TATN-1)-the first enzyme involved in the metabolic degradation of tyrosine-is up-regulated in response to oxidative stress and directly activated by the oxidative stress-responsive transcription factor SKN-1. Worms deficient in tyrosine aminotransferase activity displayed increased sensitivity to multiple sources of oxidative stress. Biochemical assays revealed that m-tyrosine is a substrate for TATN-1-mediated deamination, suggesting that TATN-1 also metabolizes m-tyrosine. Consistent with a toxic effect of m-tyrosine and a protective function of TATN-1, tatn-1 mutant worms exhibited delayed development, marked reduction in fertility, and shortened lifespan when exposed to m-tyrosine. A forward genetic screen identified a mutation in the previously uncharacterized gene F01D4.5-homologous with human transcription factor 20 (TCF20) and retinoic acid-induced 1 (RAI1)-that suppresses the adverse phenotypes observed in m-tyrosine-treated tatn-1 mutant worms. RNA-Seq analysis of F01D4.5 mutant worms disclosed a significant reduction in the expression of specific isoforms of genes encoding ribosomal proteins, suggesting that alterations in protein synthesis or ribosome structure could diminish the adverse effects of m-tyrosine. Our findings uncover a critical role for tyrosine aminotransferase in the oxidative stress response via m-tyrosine metabolism.

Biology is the root of variability: cautionary tales in Caenorhabditis elegans biology

Reproducibility is critical for the standardization, interpretation, and progression of research. However, many factors increase variability and reduce reproducibility. In Caenorhabditis elegans research, there are many possible causes of variability that may explain why experimental outcomes sometimes differ between laboratories and between experiments. Factors contributing to experimental variability include the genetic background of both C. elegans and its bacterial diet, differences in media composition, intergenerational and transgenerational effects that may be carried over for generations, and the use of chemicals or reagents that may have unexpected consequences. This review summarizes sources of variability in C. elegans research and serves to identify laboratory practices that could influence reproducibility.

Recompleting the Caenorhabditis elegans genome

Caenorhabditis elegans was the first multicellular eukaryotic genome sequenced to apparent completion. Although this assembly employed a standard C. elegans strain (N2), it used sequence data from several laboratories, with DNA propagated in bacteria and yeast. Thus, the N2 assembly has many differences from any C. elegans available today. To provide a more accurate C. elegans genome, we performed long-read assembly of VC2010, a modern strain derived from N2. Our VC2010 assembly has 99.98% identity to N2 but with an additional 1.8 Mb including tandem repeat expansions and genome duplications. For 116 structural discrepancies between N2 and VC2010, 97 structures matching VC2010 (84%) were also found in two outgroup strains, implying deficiencies in N2. Over 98% of N2 genes encoded unchanged products in VC2010; moreover, we predicted ≥53 new genes in VC2010. The recompleted genome of C. elegans should be a valuable resource for genetics, genomics, and systems biology.

Epidermal Remodeling in Caenorhabditis elegans Dauers Requires the Nidogen Domain Protein DEX-1

Phenotypic plasticity is a critical component of an organism's ability to thrive in a changing environment. The free-living nematode Caenorhabditis elegans adapts to unfavorable environmental conditions by pausing reproductive development and entering a stress-resistant larval stage known as dauer. The transition into dauer is marked by vast morphological changes, including remodeling of epidermis, neurons, and muscle. Although many of these dauer-specific traits have been described, the molecular basis of dauer-specific remodeling is still poorly understood. Here we show that the nidogen domain-containing protein DEX-1 facilitates stage-specific tissue remodeling during dauer morphogenesis. DEX-1 was previously shown to regulate sensory dendrite formation during embryogenesis. We find that DEX-1 is also required for proper remodeling of the stem cell-like epidermal seam cells. dex-1 mutant dauers lack distinct lateral cuticular alae during dauer and have increased sensitivity to sodium dodecyl sulfate. Furthermore, we find that DEX-1 is required for proper dauer mobility. We show that DEX-1 is secreted from the seam cells during dauer, but acts locally in a cell-autonomous manner. We find that dex-1 expression during dauer is regulated through DAF-16/FOXO-mediated transcriptional activation. Finally, we show that dex-1 acts with a family of zona pellucida domain-encoding genes to regulate dauer-specific epidermal remodeling. Taken together, our data indicate that DEX-1 is an extracellular matrix component that plays a central role in C. elegans epidermal remodeling during dauer.

Tubular Excretory Canal Structure Depends on Intermediate Filaments EXC-2 and IFA-4 in Caenorhabditis elegans

The excretory canals of Caenorhabditis elegans are a model for understanding the maintenance of apical morphology in narrow single-celled tubes. Light and electron microscopy shows that mutants in exc-2 start to form canals normally, but these swell to develop large fluid-filled cysts that lack a complete terminal web at the apical surface, and accumulate filamentous material in the canal lumen. Here, whole-genome sequencing and gene rescue show that exc-2 encodes intermediate filament protein IFC-2 EXC-2/IFC-2 protein, fluorescently tagged via clustered regularly interspaced short palindromic repeats/Cas9, is located at the apical surface of the canals independently of other intermediate filament proteins. EXC-2 is also located in several other tissues, though the tagged isoforms are not seen in the larger intestinal tube. Tagged EXC-2 binds via pulldown to intermediate filament protein IFA-4, which is also shown to line the canal apical surface. Overexpression of either protein results in narrow but shortened canals. These results are consistent with a model whereby three intermediate filaments in the canals-EXC-2, IFA-4, and IFB-1-restrain swelling of narrow tubules in concert with actin filaments that guide the extension and direction of tubule outgrowth, while allowing the tube to bend as the animal moves.

Whole-Genome Sequence Accuracy Is Improved by Replication in a Population of Mutagenized Sorghum

The accurate detection of induced mutations is critical for both forward and reverse genetics studies. Experimental chemical mutagenesis induces relatively few single base changes per individual. In a complex eukaryotic genome, false positive detection of mutations can occur at or above this mutagenesis rate. We demonstrate here, using a population of ethyl methanesulfonate (EMS)-treated Sorghum bicolor BTx623 individuals, that using replication to detect false positive-induced variants in next-generation sequencing (NGS) data permits higher throughput variant detection with greater accuracy. We used a lower sequence coverage depth (average of 7×) from 586 independently mutagenized individuals and detected 5,399,493 homozygous single nucleotide polymorphisms (SNPs). Of these, 76% originated from only 57,872 genomic positions prone to false positive variant calling. These positions are characterized by high copy number paralogs where the error-prone SNP positions are at copies containing a variant at the SNP position. The ability of short stretches of homology to generate these error-prone positions suggests that incompletely assembled or poorly mapped repeated sequences are one driver of these error-prone positions. Removal of these false positives left 1,275,872 homozygous and 477,531 heterozygous EMS-induced SNPs, which, congruent with the mutagenic mechanism of EMS, were >98% G:C to A:T transitions. Through this analysis, we generated a collection of sequence indexed mutants of sorghum. This collection contains 4035 high-impact homozygous mutations in 3637 genes and 56,514 homozygous missense mutations in 23,227 genes. Each line contains, on average, 2177 annotated homozygous SNPs per genome, including seven likely gene knockouts and 96 missense mutations. The number of mutations in a transcript was linearly correlated with the transcript length and also the G+C count, but not with the GC/AT ratio. Analysis of the detected mutagenized positions identified CG-rich patches, and flanking sequences strongly influenced EMS-induced mutation rates. This method for detecting false positive-induced mutations is generally applicable to any organism, is independent of the choice of in silico variant-calling algorithm, and is most valuable when the true mutation rate is likely to be low, such as in laboratory-induced mutations or somatic mutation detection in medicine.

Genomic Identification and Functional Characterization of Essential Genes in Caenorhabditis elegans

Using combined genetic mapping, Illumina sequencing, bioinformatics analyses, and experimental validation, we identified 60 essential genes from 104 lethal mutations in two genomic regions of Caenorhabditis elegans totaling ∼14 Mb on chromosome III(mid) and chromosome V(left). Five of the 60 genes had not previously been shown to have lethal phenotypes by RNA interference depletion. By analyzing the regions around the lethal missense mutations, we identified four putative new protein functional domains. Furthermore, functional characterization of the identified essential genes shows that most are enzymes, including helicases, tRNA synthetases, and kinases in addition to ribosomal proteins. Gene Ontology analysis indicated that essential genes often encode for enzymes that conduct nucleic acid binding activities during fundamental processes, such as intracellular DNA replication, transcription, and translation. Analysis of essential gene shows that they have fewer paralogs, encode proteins that are in protein interaction hubs, and are highly expressed relative to nonessential genes. All these essential gene traits in C. elegans are consistent with those of human disease genes. Most human orthologs (90%) of the essential genes in this study are related to human diseases. Therefore, functional characterization of essential genes underlines their importance as proxies for understanding the biological functions of human disease genes.

Use of a Sibling Subtraction Method for Identifying Causal Mutations in Caenorhabditis elegans by Whole-Genome Sequencing

Whole-genome sequencing (WGS) is an indispensable tool for identifying causal mutations obtained from genetic screens. To reduce the number of causal mutation candidates typically uncovered by WGS, Caenorhabditis elegans researchers have developed several strategies. One involves crossing N2-background mutants to the polymorphic Hawaiian (HA) strain, which can be used to simultaneously identify mutant strain variants and obtain high-density mapping information. This approach, however, is not well suited for uncovering mutations in complex genetic backgrounds, and HA polymorphisms can alter phenotypes. Other approaches make use of DNA variants present in the initial background or introduced by mutagenesis. This information is used to implicate genomic regions with high densities of DNA lesions that persist after backcrossing, but these methods can provide lower resolution than HA mapping. To identify suppressor mutations using WGS, we developed an approach termed the sibling subtraction method (SSM). This method works by eliminating variants present in both mutants and their nonmutant siblings, thus greatly reducing the number of candidates. We used this method with two members of the C. elegans NimA-related kinase family, nekl-2 and nekl-3. Combining weak aphenotypic alleles of nekl-2 and nekl-3 leads to penetrant molting defects and larval arrest. We isolated ∼50 suppressors of nekl-2; nekl-3 synthetic lethality using F1 clonal screening methods and a peel-1–based counterselection strategy. When applied to five of the suppressors, SSM led to only one to four suppressor candidates per strain. Thus SSM is a powerful approach for identifying causal mutations in any genetic background and provides an alternative to current methods.

MinION-based long-read sequencing and assembly extends the Caenorhabditis elegans reference genome

Advances in long-read single molecule sequencing have opened new possibilities for 'benchtop' whole-genome sequencing. The Oxford Nanopore Technologies MinION is a portable device that uses nanopore technology that can directly sequence DNA molecules. MinION single molecule long sequence reads are well suited for de novo assembly of complex genomes as they facilitate the construction of highly contiguous physical genome maps obviating the need for labor-intensive physical genome mapping. Long sequence reads can also be used to delineate complex chromosomal rearrangements, such as those that occur in tumor cells, that can confound analysis using short reads. Here, we assessed MinION long-read-derived sequences for feasibility concerning: (1) the de novo assembly of a large complex genome, and (2) the elucidation of complex rearrangements. The genomes of two Caenorhabditis elegans strains, a wild-type strain and a strain containing two complex rearrangements, were sequenced with MinION. Up to 42-fold coverage was obtained from a single flow cell, and the best pooled data assembly produced a highly contiguous wild-type C. elegans genome containing 48 contigs (N50 contig length = 3.99 Mb) covering >99% of the 100,286,401-base reference genome. Further, the MinION-derived genome assembly expanded the C. elegans reference genome by >2 Mb due to a more accurate determination of repetitive sequence elements and assembled the complete genomes of two co-extracted bacteria. MinION long-read sequence data also facilitated the elucidation of complex rearrangements in a mutagenized strain. The sequence accuracy of the MinION long-read contigs (∼98%) was improved using Illumina-derived sequence data to polish the final genome assembly to 99.8% nucleotide accuracy when compared to the reference assembly.

Different mutational function of low- and high-linear energy transfer heavy-ion irradiation demonstrated by whole-genome resequencing of Arabidopsis mutants

Heavy-ion irradiation is a powerful mutagen that possesses high linear energy transfer (LET). Several studies have indicated that the value of LET affects DNA lesion formation in several ways, including the efficiency and the density of double-stranded break induction along the particle path. We assumed that the mutation type can be altered by selecting an appropriate LET value. Here, we quantitatively demonstrate differences in the mutation type induced by irradiation with two representative ions, Ar ions (LET: 290 keV μm^-1 ) and C ions (LET: 30.0 keV μm^-1 ), by whole-genome resequencing of the Arabidopsis mutants produced by these irradiations. Ar ions caused chromosomal rearrangements or large deletions (≥100 bp) more frequently than C ions, with 10.2 and 2.3 per mutant genome under Ar- and C-ion irradiation, respectively. Conversely, C ions induced more single-base substitutions and small indels (<100 bp) than Ar ions, with 28.1 and 56.9 per mutant genome under Ar- and C-ion irradiation, respectively. Moreover, the rearrangements induced by Ar-ion irradiation were more complex than those induced by C-ion irradiation, and tended to accompany single base substitutions or small indels located close by. In conjunction with the detection of causative genes through high-throughput sequencing, selective irradiation by beams with different effects will be a powerful tool for forward genetics as well as studies on chromosomal rearrangements.

Applying inbreeding, hybridization and mutagenesis to improve oxidative stress tolerance and longevity of the entomopathogenic nematode Heterorhabditis bacteriophora

Poor shelf-life and sensitivity to environmental stress of entomopathogenic nematodes (EPNs) are traits, which deserve attention for improvement. Recently, a strong positive correlation between oxidative stress tolerance and longevity of Heterorhabditis bacteriophora dauer juveniles (DJs) has been reported. In this study, the improvement of H. bacteriophora DJ longevity was achieved by hybridization and mutagenesis. A hybrid pool deriving from two oxidative stress tolerant and long-living parental strains was generated. This hybrid AU1 × HU2 survived 2.6 days and 18 days longer than its best parent under oxidative stress and control conditions, respectively. In addition to the natural genetic variability, an EMS-mutant pool (M-OXI) with high longevity was generated and one of the derived mutagenized inbred lines (MOX-IL6) survived 5.8 days and 28.4 days longer than its donor line (IL3) under oxidative stress and control conditions, respectively. A genetic cross between the mutagenized inbred line and its donor line (MOX-IL × IL3) still survived 2.5 days and 18.5 days longer than the donor line under oxidative stress and control conditions, respectively. Concerning virulence and reproductive potential, trade-off effects were not observed as a result of hybridization and mutagenesis. These results underline the potential of classical genetic approaches for trait improvement in the nematode H. bacteriophora.

Mitochondrial Mutation Rate, Spectrum and Heteroplasmy in Caenorhabditis elegans Spontaneous Mutation Accumulation Lines of Differing Population Size

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.

Targeted mutagenesis in a human-parasitic nematode

Parasitic nematodes infect over 1 billion people worldwide and cause some of the most common neglected tropical diseases. Despite their prevalence, our understanding of the biology of parasitic nematodes has been limited by the lack of tools for genetic intervention. In particular, it has not yet been possible to generate targeted gene disruptions and mutant phenotypes in any parasitic nematode. Here, we report the development of a method for introducing CRISPR-Cas9-mediated gene disruptions in the human-parasitic threadworm Strongyloides stercoralis. We disrupted the S. stercoralis twitchin gene unc-22, resulting in nematodes with severe motility defects. Ss-unc-22 mutations were resolved by homology-directed repair when a repair template was provided. Omission of a repair template resulted in deletions at the target locus. Ss-unc-22 mutations were heritable; we passed Ss-unc-22 mutants through a host and successfully recovered mutant progeny. Using a similar approach, we also disrupted the unc-22 gene of the rat-parasitic nematode Strongyloides ratti. Our results demonstrate the applicability of CRISPR-Cas9 to parasitic nematodes, and thereby enable future studies of gene function in these medically relevant but previously genetically intractable parasites.

Parasitic worms are a widespread public health burden, yet very little is known about the cellular and molecular mechanisms that contribute to their parasitic lifestyle. One of the major barriers to better understanding these mechanisms is that there are currently no available methods for making targeted gene knockouts in any parasitic worm species. Here, we describe the first mutant phenotype in a parasitic worm resulting from a targeted gene disruption. We applied CRISPR-Cas9-mediated mutagenesis to parasitic worms in the genus Strongyloides and developed a method that overcomes many of the challenges that have previously inhibited generating mutant parasitic worms. We characterize heritable mutant phenotypes and outline a toolkit that will be applicable to many other genes with potential roles in parasitism. Importantly, we developed our method for gene knockouts in a human-parasitic worm. By directly investigating the genes and molecular pathways that enable worms to parasitize humans, we may be able to develop novel anthelmintic therapies or other measures for preventing nematode infections.

PQN-75 is expressed in the pharyngeal gland cells of Caenorhabditiselegans and is dispensable for germline development

In Caenorhabditis elegans, five pharyngeal gland cells reside in the terminal bulb of the pharynx and extend anterior processes to five contact points in the pharyngeal lumen. Pharyngeal gland cells secrete mucin-like proteins thought to facilitate digestion, hatching, molting and assembly of the surface coat of the cuticle, but supporting evidence has been sparse. Here we show pharyngeal gland cell expression of PQN-75, a unique protein containing an N-terminal signal peptide, nucleoporin (Nup)-like phenylalanine/glycine (FG) repeats, and an extensive polyproline repeat domain with similarities to human basic salivary proline-rich pre-protein PRB2. Imaging of C-terminal tagged PQN-75 shows localization throughout pharyngeal gland cell processes but not the pharyngeal lumen; instead, aggregates of PQN-75 are occasionally found throughout the pharynx, suggesting secretion from pharyngeal gland cells into the surrounding pharyngeal muscle. PQN-75 does not affect fertility and brood size in C. elegans but confers some degree of stress resistance and thermotolerance through unknown mechanisms.

Summary: PQN-75 is expressed in pharyngeal gland cells and shares similarity with human basic salivary proline-rich protein PBR2, suggesting evolutionary conservation between gland cells in the upper digestive tract.

Revitalization of a Forward Genetic Screen Identifies Three New Regulators of Fungal Secondary Metabolism in the Genus Aspergillus

The study of aflatoxin in Aspergillus spp. has garnered the attention of many researchers due to aflatoxin's carcinogenic properties and frequency as a food and feed contaminant. Significant progress has been made by utilizing the model organism Aspergillus nidulans to characterize the regulation of sterigmatocystin (ST), the penultimate precursor of aflatoxin. A previous forward genetic screen identified 23 A. nidulans mutants involved in regulating ST production. Six mutants were characterized from this screen using classical mapping (five mutations in mcsA) and complementation with a cosmid library (one mutation in laeA). The remaining mutants were backcrossed and sequenced using Illumina and Ion Torrent sequencing platforms. All but one mutant contained one or more sequence variants in predicted open reading frames. Deletion of these genes resulted in identification of mutant alleles responsible for the loss of ST production in 12 of the 17 remaining mutants. Eight of these mutations were in genes already known to affect ST synthesis (laeA, mcsA, fluG, and stcA), while the remaining four mutations (in laeB, sntB, and hamI) were in previously uncharacterized genes not known to be involved in ST production. Deletion of laeB, sntB, and hamI in A. flavus results in loss of aflatoxin production, confirming that these regulators are conserved in the aflatoxigenic aspergilli. This report highlights the multifaceted regulatory mechanisms governing secondary metabolism in Aspergillus Additionally, these data contribute to the increasing number of studies showing that forward genetic screens of fungi coupled with whole-genome resequencing is a robust and cost-effective technique.IMPORTANCE In a postgenomic world, reverse genetic approaches have displaced their forward genetic counterparts. The techniques used in forward genetics to identify loci of interest were typically very cumbersome and time-consuming, relying on Mendelian traits in model organisms. The current work was pursued not only to identify alleles involved in regulation of secondary metabolism but also to demonstrate a return to forward genetics to track phenotypes and to discover genetic pathways that could not be predicted through a reverse genetics approach. While identification of mutant alleles from whole-genome sequencing has been done before, here we illustrate the possibility of coupling this strategy with a genetic screen to identify multiple alleles of interest. Sequencing of classically derived mutants revealed several uncharacterized genes, which represent novel pathways to regulate and control the biosynthesis of sterigmatocystin and of aflatoxin, a societally and medically important mycotoxin.

MIP-MAP: High-Throughput Mapping of Caenorhabditis elegans Temperature-Sensitive Mutants via Molecular Inversion Probes

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.

Toward Universal Forward Genetics: Using a Draft Genome Sequence of the Nematode Oscheius tipulae To Identify Mutations Affecting Vulva Development

Mapping-by-sequencing has become a standard method to map and identify phenotype-causing mutations in model species. Here, we show that a fragmented draft assembly is sufficient to perform mapping-by-sequencing in nonmodel species. We generated a draft assembly and annotation of the genome of the free-living nematode Oscheius tipulae, a distant relative of the model Caenorhabditis elegans We used this draft to identify the likely causative mutations at the O. tipulae cov-3 locus, which affect vulval development. The cov-3 locus encodes the O. tipulae ortholog of C. elegans mig-13, and we further show that Cel-mig-13 mutants also have an unsuspected vulval-development phenotype. In a virtuous circle, we were able to use the linkage information collected during mutant mapping to improve the genome assembly. These results showcase the promise of genome-enabled forward genetics in nonmodel species.