Search

Although many studies have examined quantitative trait variation across many species, only a small number of genes and thereby molecular mechanisms have been discovered. Without these data, we can only speculate about evolutionary processes that underlie trait variation. Here, we review how quantitative and molecular genetics in the nematode Caenorhabditis elegans led to the discovery and validation of 37 quantitative trait genes over the past 15 years. Using these data, we can start to make inferences about evolution from these quantitative trait genes, including the roles that coding versus noncoding variation, gene family expansion, common versus rare variants, pleiotropy, and epistasis play in trait variation across this species. Innovations in quantitative trait loci mapping and genome editing have led to the discovery and validation of 37 genes and variants underlying phenotypic variation in C. elegans. Numerous recombinant panels and a large collection of wild strains make C. elegans a formidable model to understand quantitative trait variation. Most of the identified quantitative trait genes have paralogs, providing evidence that gene duplication events are important for shaping quantitative traits. Pleiotropy is relatively common among C. elegans quantitative trait genes. [ABSTRACT FROM AUTHOR]

High-throughput imaging techniques have become widespread in many fields of biology. These powerful platforms generate large quantities of data that can be difficult to process and visualize efficiently using existing tools. We developed easyXpress to process and review C. elegans high-throughput microscopy data in the R environment. The package provides a logical workflow for the reading, analysis, and visualization of data generated using CellProfiler's WormToolbox. We equipped easyXpress with powerful functions to customize the filtering of noise in data, specifically by identifying and removing objects that deviate from expected animal measurements. This flexibility in data filtering allows users to optimize their analysis pipeline to match their needs. In addition, easyXpress includes tools for generating detailed visualizations, allowing the user to interactively compare summary statistics across wells and plates with ease. Researchers studying C. elegans benefit from this streamlined and extensible package as it is complementary to CellProfiler and leverages the R environment to rapidly process and analyze large high-throughput imaging datasets. [ABSTRACT FROM AUTHOR]

Phenotypic complexity is caused by the contributions of environmental factors and multiple genetic loci, interacting or acting independently. Studies of yeast and Arabidopsis often find that the majority of natural variation across phenotypes is attributable to independent additive quantitative trait loci (QTL). Detected loci in these organisms explain most of the estimated heritable variation. By contrast, many heritable components underlying phenotypic variation in metazoan models remain undetected. Before the relative impacts of additive and interactive variance components on metazoan phenotypic variation can be dissected, high replication and precise phenotypic measurements are required to obtain sufficient statistical power to detect loci contributing to this missing heritability. Here, we used a panel of 296 recombinant inbred advanced intercross lines of Caenorhabditis elegans and a highthroughput fitness assay to detect loci underlying responses to 16 different toxins, including heavy metals, chemotherapeutic drugs, pesticides, and neuropharmaceuticals. Using linkage mapping, we identified 82 QTL that underlie variation in responses to these toxins, and predicted the relative contributions of additive loci and genetic interactions across various growth parameters. Additionally, we identified three genomic regions that impact responses to multiple classes of toxins. These QTL hotspots could represent common factors impacting toxin responses. We went further to generate near-isogenic lines and chromosome substitution strains, and then experimentally validated these QTL hotspots, implicating additive and interactive loci that underlie toxin-response variation. [ABSTRACT FROM AUTHOR]

Telomeres are involved in the maintenance of chromosomes and the prevention of genome instability. Despite this central importance, significant variation in telomere length has been observed in a variety of organisms. The genetic determinants of telomerelength variation and their effects on organismal fitness are largely unexplored. Here, we describe natural variation in telomere length across the Caenorhabditis elegans species. We identify a large-effect variant that contributes to differences in telomere length. The variant alters the conserved oligonucleotide/oligosaccharide-binding fold of protection of telomeres 2 (POT-2), a homolog of a human telomere-capping shelterin complex subunit. Mutations within this domain likely reduce the ability of POT-2 to bind telomeric DNA, thereby increasing telomere length. We find that telomere-length variation does not correlate with offspring production or longevity in C. elegans wild isolates, suggesting that naturally long telomeres play a limited role in modifying fitness phenotypes in C. elegans. [ABSTRACT FROM AUTHOR]

The Hawaiian strain (CB4856) of Caenorhabditis elegans is one of the most divergent from the canonical laboratory strain N2 and has been widely used in developmental, population, and evolutionary studies. To enhance the utility of the strain, we have generated a draft sequence of the CB4856 genome, exploiting a variety of resources and strategies. When compared against the N2 reference, the CB4856 genome has 327,050 single nucleotide variants (SNVs) and 79,529 insertion--deletion events that result in a total of 3.3 Mb of N2 sequence missing from CB4856 and 1.4 Mb of sequence present in CB4856 but not present in N2. As previously reported, the density of SNVs varies along the chromosomes, with the arms of chromosomes showing greater average variation than the centers. In addition, we find 61 regions totaling 2.8 Mb, distributed across all six chromosomes, which have a greatly elevated SNV density, ranging from 2 to 16% SNVs. A survey of other wild isolates show that the two alternative haplotypes for each region are widely distributed, suggesting they have been maintained by balancing selection over long evolutionary times. These divergent regions contain an abundance of genes from large rapidly evolving families encoding F-box, MATH, BATH, seven-transmembrane G-coupled receptors, and nuclear hormone receptors, suggesting that they provide selective advantages in natural environments. The draft sequence makes available a comprehensive catalog of sequence differences between the CB4856 and N2 strains that will facilitate the molecular dissection of their phenotypic differences. Our work also emphasizes the importance of going beyond simple alignment of reads to a reference genome when assessing differences between genomes. [ABSTRACT FROM AUTHOR]

The genetic variability of toxicant responses among indisviduals in humans and mammalian models requires practically untenable sample sizes to create comprehensive chemical hazard risk evaluations. To address this need, tractable model systems enable reproducible and efficient experimental workflows to collect high-replication measurements of exposure cohorts. Caenorhabditis elegans is a premier toxicology model that has revolutionized our understanding of cellular responses to environmental pollutants and boasts robust genomic resources and high levels of genetic variation across the species. In this study, we performed dose-response analysis across 23 environmental toxicants using eight C. elegans strains representative of species-wide genetic diversity. We observed substantial variation in EC10 estimates and slope parameter estimates of dose-response curves of different strains, demonstrating that genetic background is a significant driver of differential toxicant susceptibility. We also showed that, across all toxicants, at least one C. elegans strain exhibited a significantly different EC10 or slope estimate compared to the reference strain, N2 (PD1074), indicating that population-wide differences among strains are necessary to understand responses to toxicants. Moreover, we quantified the heritability of responses (phenotypic variance attributable to genetic differences between individuals) to each toxicant exposure and observed a correlation between the exposure closest to the species-agnostic EC10 estimate and the exposure that exhibited the most heritable response. At least 20% of the variance in susceptibility to at least one exposure level of each compound was explained by genetic differences among the eight C. elegans strains. Taken together, these results provide robust evidence that heritable genetic variation explains differential susceptibility across an array of environmental pollutants and that genetically diverse C. elegans strains should be deployed to aid high-throughput toxicological screening efforts. [ABSTRACT FROM AUTHOR]

Anthelmintic resistance (AR) in parasitic nematodes poses a global health problem in livestock and domestic animals and is an emerging problem in humans. Consequently, we must understand the mechanisms of AR, including target-site resistance (TSR), in which mutations affect drug binding, and non-target site resistance (NTSR), which involves alterations in drug metabolism and detoxification processes. Because much of the focus has been on TSR, NTSR has received less attention. Here, we describe how metabolomics approaches using Caenorhabditis elegans offer the ability to disentangle nematode drug metabolism, identify metabolic changes associated with resistance, uncover novel biomarkers, and enhance diagnostic methods., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Transposable elements (TEs) can alter host gene structure and expression, whereas host organisms develop mechanisms to repress TE activities. In the nematode Caenorhabditis elegans , a small interfering RNA pathway dependent on the helicase ERI-6/7 primarily silences retrotransposons and recent genes of likely viral origin. By studying gene expression variation among wild C. elegans strains, we found that structural variants and transposon remnants likely underlie expression variation in eri-6/7 and the pathway targets. We further found that multiple insertions of the DNA transposons, Polintons, reshuffled the eri-6/7 locus and induced inversion of eri-6 in some wild strains. In the inverted configuration, gene function was previously shown to be repaired by unusual trans-splicing mediated by direct repeats. We identified that these direct repeats originated from terminal inverted repeats of Polintons . Our findings highlight the role of host-transposon interactions in driving rapid host genome diversification among natural populations and shed light on evolutionary novelty in genes and splicing mechanisms.

Ascaridia galli and Ascaridia dissimilis are the most common and economically impactful nematode parasites of commercial poultry. These infections rarely cause clinical disease, but reduction in feed conversion efficiency is detected. To determine if feed conversion efficiency reductions correlate with any physiological measures independent of clinical disease, we determined if ascarid infections correlate with changes in the weights of the small intestine, liver, or total animal weight (quantitative measures of animal health). No correlation between parasite burden and these metrics were observed, supporting the concept that feed conversion is the only production metric impacted by ascarid infections., Competing Interests: The authors declare that there are no conflicts of interest present., (Copyright: © 2024 by the authors.)

Benzimidazole (BZ) anthelmintics are among the most important treatments for parasitic nematode infections in the developing world. Widespread BZ resistance in veterinary parasites and emerging resistance in human parasites raise major concerns for the continued use of BZs. Knowledge of the mechanisms of resistance is necessary to make informed treatment decisions and circumvent resistance. Benzimidazole resistance has traditionally been associated with mutations and natural variants in the C. elegans beta-tubulin gene ben-1 and orthologs in parasitic species. However, variants in ben-1 alone do not explain the differences in BZ responses across parasite populations. Here, we examined the roles of five C. elegans beta-tubulin genes (tbb-1, mec-7, tbb-4, ben-1, and tbb-6) in the BZ response as well as to determine if another beta-tubulin acts redundantly with ben-1. We generated C. elegans strains with a loss of each beta-tubulin gene, as well as strains with a loss of tbb-1, mec-7, tbb-4, or tbb-6 in a genetic background that also lacks ben-1. We found that the loss of ben-1 conferred the maximum level of resistance following exposure to a single concentration of albendazole, and the loss of a second beta-tubulin gene did not alter the level of resistance. However, additional traits other than larval development could be affected by the loss of additional beta-tubulins, and the roles of other beta-tubulin genes might be revealed at different albendazole concentrations. Therefore, further work is needed to fully define the possible roles of other beta-tubulin genes in the BZ response., Competing Interests: Declaration of competing interest None., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Phoresy is an interspecies interaction that facilitates spatial dispersal by attaching to a more mobile species. Hitchhiking species have evolved specific traits for physical contact and successful phoresy, but the regulatory mechanisms involved in such traits and their evolution are largely unexplored. The nematode Caenorhabditis elegans displays a hitchhiking behavior known as nictation during its stress-induced developmental stage. Dauer-specific nictation behavior has an important role in natural C. elegans populations, which experience boom-and-bust population dynamics. In this study, we investigated the nictation behavior of 137 wild C. elegans strains sampled throughout the world. We identified species-wide natural variation in nictation and performed a genome-wide association mapping. We show that the variants in the promoter of nta-1 , encoding a putative steroidogenic enzyme, underlie differences in nictation. This difference is due to the changes in nta-1 expression in glial cells, which implies that glial steroid metabolism regulates phoretic behavior. Population genetic analysis and geographic distribution patterns suggest that balancing selection maintained two nta-1 haplotypes that existed in ancestral C. elegans populations. Our findings contribute to further understanding of the molecular mechanism of species interaction and the maintenance of genetic diversity within natural populations., Competing Interests: Competing interests statement:The authors declare no competing interest.

Left unchecked, plant-parasitic nematodes have the potential to devastate crops globally. Highly effective but non-selective nematicides are justifiably being phased-out, leaving farmers with limited options for managing nematode infestation. Here, we report our discovery of a 1,3,4-oxadiazole thioether scaffold called Cyprocide that selectively kills nematodes including diverse species of plant-parasitic nematodes. Cyprocide is bioactivated into a lethal reactive electrophilic metabolite by specific nematode cytochrome P450 enzymes. Cyprocide fails to kill organisms beyond nematodes, suggesting that the targeted lethality of this pro-nematicide derives from P450 substrate selectivity. Our findings demonstrate that Cyprocide is a selective nematicidal scaffold with broad-spectrum activity that holds the potential to help safeguard our global food supply., (© 2024. The Author(s).)

The molecular mechanisms underlying diversity in animal behavior are not well understood. A major experimental challenge is determining the contribution of genetic variants that affect neuronal gene expression to differences in behavioral traits. In Caenorhabditis elegans , the neuroendocrine transforming growth factor-β ligand, DAF-7, regulates diverse behavioral responses to bacterial food and pathogens. The dynamic neuron-specific expression of daf-7 is modulated by environmental and endogenous bacteria-derived cues. Here, we investigated natural variation in the expression of daf-7 from the ASJ pair of chemosensory neurons. We identified common genetic variants in gap-2 , encoding a Ras guanosine triphosphatase (GTPase)-activating protein homologous to mammalian synaptic Ras GTPase-activating protein, which modify daf-7 expression cell nonautonomously and promote exploratory foraging behavior in a partially DAF-7-dependent manner. Our data connect natural variation in neuron-specific gene expression to differences in behavior and suggest that genetic variation in neuroendocrine signaling pathways mediating host-microbe interactions may give rise to diversity in animal behavior.

Differences in naïve alcohol sensitivity between individuals are a strong predictor of later life alcohol use disorders (AUD). However, the genetic bases for alcohol sensitivity (beyond ethanol metabolism) and pharmacological approaches to modulate alcohol sensitivity remain poorly understood. We used a high-throughput behavioral screen to measure acute behavioral sensitivity to alcohol, a model of intoxication, in a genetically diverse set of over 150 wild strains of the nematode Caenorhabditis elegans . We performed a genome-wide association study to identify loci that underlie natural variation in alcohol sensitivity. We identified five quantitative trait loci (QTL) and further show that variants in the C. elegans ortholog of protein kinase D, dkf-2 , likely underlie the chromosome V QTL. We found that resistance to intoxication was conferred by dkf-2 loss-of-function mutations as well as partly by a PKD inhibitor in a dkf-2 -dependent manner. Protein kinase D might represent a conserved, druggable target to modify alcohol sensitivity with application towards AUD., Competing Interests: Competing interests: The authors declare no competing interest.

Albendazole (a benzimidazole) and ivermectin (a macrocyclic lactone) are the two most commonly co-administered anthelmintic drugs in mass-drug administration programs worldwide. Despite emerging resistance, we do not fully understand the mechanisms of resistance to these drugs nor the consequences of delivering them in combination. Albendazole resistance has primarily been attributed to variation in the drug target, a beta-tubulin gene. Ivermectin targets glutamate-gated chloride channels (GluCls), but it is unknown whether GluCl genes are involved in ivermectin resistance in nature. Using Caenorhabditis elegans, we defined the fitness costs associated with loss of the drug target genes singly or in combinations of the genes that encode GluCl subunits. We quantified the loss-of-function effects on three traits: (i) multi-generational competitive fitness, (ii) fecundity, and (iii) development. In competitive fitness and development assays, we found that a deletion of the beta-tubulin gene ben-1 conferred albendazole resistance, but ivermectin resistance required the loss of two GluCl genes (avr-14 and avr-15). The fecundity assays revealed that loss of ben-1 did not provide any fitness benefit in albendazole conditions and that no GluCl deletion mutants were resistant to ivermectin. Next, we searched for evidence of multi-drug resistance across the three traits. Loss of ben-1 did not confer resistance to ivermectin, nor did loss of any single GluCl subunit or combination confer resistance to albendazole. Finally, we assessed the development of 124 C. elegans wild strains across six benzimidazoles and seven macrocyclic lactones to identify evidence of multi-drug resistance between the two drug classes and found a strong phenotypic correlation within a drug class but not across drug classes. Because each gene affects various aspects of nematode physiology, these results suggest that it is necessary to assess multiple fitness traits to evaluate how each gene contributes to anthelmintic resistance., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Shaver et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Studies of model organisms have provided important insights into how natural genetic differences shape trait variation. These discoveries are driven by the growing availability of genomes and the expansive experimental toolkits afforded to researchers using these species. For example, Caenorhabditis elegans is increasingly being used to identify and measure the effects of natural genetic variants on traits using quantitative genetics. Since 2016, the C. elegans Natural Diversity Resource (CeNDR) has facilitated many of these studies by providing an archive of wild strains, genome-wide sequence and variant data for each strain, and a genome-wide association (GWA) mapping portal for the C. elegans community. Here, we present an updated platform, the Caenorhabditis Natural Diversity Resource (CaeNDR), that enables quantitative genetics and genomics studies across the three Caenorhabditis species: C. elegans, C. briggsae and C. tropicalis. The CaeNDR platform hosts several databases that are continually updated by the addition of new strains, whole-genome sequence data and annotated variants. Additionally, CaeNDR provides new interactive tools to explore natural variation and enable GWA mappings. All CaeNDR data and tools are accessible through a freely available web portal located at caendr.org., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The ascarids are a large group of parasitic nematodes that infect a wide range of animal species. In humans, they cause neglected diseases of poverty; many animal parasites also cause zoonotic infections in people. Control measures include hygiene and anthelmintic treatments, but they are not always appropriate or effective and this creates a continuing need to search for better ways to reduce the human, welfare and economic costs of these infections. To this end, Le Studium Institute of Advanced Studies organized a two-day conference to identify major gaps in our understanding of ascarid parasites with a view to setting research priorities that would allow for improved control. The participants identified several key areas for future focus, comprising of advances in genomic analysis and the use of model organisms, especially Caenorhabditis elegans, a more thorough appreciation of the complexity of host-parasite (and parasite-parasite) communications, a search for novel anthelmintic drugs and the development of effective vaccines. The participants agreed to try and maintain informal links in the future that could form the basis for collaborative projects, and to co-operate to organize future meetings and workshops to promote ascarid research., (Copyright © 2024. Published by Elsevier Ltd.)

Background: The nematode Caenorhabditis briggsae has been used as a model in comparative genomics studies with Caenorhabditis elegans because of their striking morphological and behavioral similarities. However, the potential of C. briggsae for comparative studies is limited by the quality of its genome resources. The genome resources for the C. briggsae laboratory strain AF16 have not been developed to the same extent as C. elegans. The recent publication of a new chromosome-level reference genome for QX1410, a C. briggsae wild strain closely related to AF16, has provided the first step to bridge the gap between C. elegans and C. briggsae genome resources. Currently, the QX1410 gene models consist of software-derived gene predictions that contain numerous errors in their structure and coding sequences. In this study, a team of researchers manually inspected over 21,000 gene models and underlying transcriptomic data to repair software-derived errors., Results: We designed a detailed workflow to train a team of nine students to manually curate gene models using RNA read alignments. We manually inspected the gene models, proposed corrections to the coding sequences of over 8,000 genes, and modeled thousands of putative isoforms and untranslated regions. We exploited the conservation of protein sequence length between C. briggsae and C. elegans to quantify the improvement in protein-coding gene model quality and showed that manual curation led to substantial improvements in the protein sequence length accuracy of QX1410 genes. Additionally, collinear alignment analysis between the QX1410 and AF16 genomes revealed over 1,800 genes affected by spurious duplications and inversions in the AF16 genome that are now resolved in the QX1410 genome., Conclusions: Community-based, manual curation using transcriptome data is an effective approach to improve the quality of software-derived protein-coding genes. The detailed protocols provided in this work can be useful for future large-scale manual curation projects in other species. Our manual curation efforts have brought the QX1410 gene models to a comparable level of quality as the extensively curated AF16 gene models. The improved genome resources for C. briggsae provide reliable tools for the study of Caenorhabditis biology and other related nematodes., (© 2023. BioMed Central Ltd., part of Springer Nature.)

Anthelmintic drugs are used to treat parasitic roundworm and flatworm infections in humans and other animals. Caenorhabditis elegans is an established model to investigate anthelmintics used to treat roundworms. In this study, we use C. elegans to examine the mode of action and the mechanisms of resistance against the flatworm anthelmintic drug praziquantel (PZQ), used to treat trematode and cestode infections. We found that PZQ inhibited development and that this developmental delay varies by genetic background. Interestingly, both enantiomers of PZQ are equally effective against C. elegans, but the right-handed PZQ (R-PZQ) is most effective against schistosome infections. We conducted a genome-wide association mapping with 74 wild C. elegans strains to identify a region on chromosome IV that is correlated with differential PZQ susceptibility. Five candidate genes in this region: cct-8, znf-782, Y104H12D.4, Y104H12D.2, and cox-18, might underlie this variation. The gene cct-8, a subunit of the protein folding complex TRiC, has variation that causes a putative protein coding change (G226V), which is correlated with reduced developmental delay. Gene expression analysis suggests that this variant correlates with slightly increased expression of both cct-8 and hsp-70. Acute exposure to PZQ caused increased expression of hsp-70, indicating that altered TRiC function might be involved in PZQ responses. To test if this variant affects development upon exposure to PZQ, we used CRISPR-Cas9 genome editing to introduce the V226 allele into the N2 genetic background (G226) and the G226 allele into the JU775 genetic background (V226). These experiments revealed that this variant was not sufficient to explain the effects of PZQ on development. Nevertheless, this study shows that C. elegans can be used to study PZQ mode of action and resistance mechanisms. Additionally, we show that the TRiC complex requires further evaluation for PZQ responses in C. elegans., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Wit et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Parasitic nematodes cause significant effects on humans each year, with the most prevalent being Ascaris lumbricoides. Benzimidazoles (BZ) are the most widely used anthelmintic drug in humans, and although the biology of resistance to this drug class is understood in some species, resistance is poorly characterized in ascarids. Models such as Caenorhabditis elegans were essential in developing our current understanding of BZ resistance, but more closely related model nematodes are needed to understand resistance in ascarids. Here, we propose a new ascarid model species that infects turkeys, Ascaridia dissimilis, to develop a better understanding of BZ resistance., (Copyright © 2022 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

From bacterial quorum sensing to human language, communication is essential for social interactions. Nematodes produce and sense pheromones to communicate among individuals and respond to environmental changes. These signals are encoded by different types and mixtures of ascarosides, whose modular structures further enhance the diversity of this nematode pheromone language. Interspecific and intraspecific differences in this ascaroside pheromone language have been described previously, but the genetic basis and molecular mechanisms underlying the variation remain largely unknown. Here, we analyzed natural variation in the production of 44 ascarosides across 95 wild Caenorhabditis elegans strains using high-performance liquid chromatography coupled to high-resolution mass spectrometry. We discovered wild strains defective in the production of specific subsets of ascarosides ( e.g. , the aggregation pheromone icas#9) or short- and medium-chain ascarosides, as well as inversely correlated patterns between the production of two major classes of ascarosides. We investigated genetic variants that are significantly associated with the natural differences in the composition of the pheromone bouquet, including rare genetic variants in key enzymes participating in ascaroside biosynthesis, such as the peroxisomal 3-ketoacyl-CoA thiolase, daf-22 , and the carboxylesterase cest-3 . Genome-wide association mappings revealed genomic loci harboring common variants that affect ascaroside profiles. Our study yields a valuable dataset for investigating the genetic mechanisms underlying the evolution of chemical communication.

Short tandem repeats (STRs) have orders of magnitude higher mutation rates than single nucleotide variants (SNVs) and have been proposed to accelerate evolution in many organisms. However, only few studies have addressed the impact of STR variation on phenotypic variation at both the organismal and molecular levels. Potential driving forces underlying the high mutation rates of STRs also remain largely unknown. Here, we leverage the recently generated expression and STR variation data among wild Caenorhabditis elegans strains to conduct a genome-wide analysis of how STRs affect gene expression variation. We identify thousands of expression STRs (eSTRs) showing regulatory effects and demonstrate that they explain missing heritability beyond SNV-based expression quantitative trait loci. We illustrate specific regulatory mechanisms such as how eSTRs affect splicing sites and alternative splicing efficiency. We also show that differential expression of antioxidant genes and oxidative stresses might affect STR mutations systematically using both wild strains and mutation accumulation lines. Overall, we reveal the interplay between STRs and gene expression variation by providing novel insights into regulatory mechanisms of STRs and highlighting that oxidative stress could lead to higher STR mutation rates., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Treatment of parasitic nematode infections in humans and livestock relies on a limited arsenal of anthelmintic drugs that have historically reduced parasite burdens. However, anthelmintic resistance (AR) is increasing, and little is known about the molecular and genetic causes of resistance for most drugs. The free-living roundworm Caenorhabditis elegans has proven to be a tractable model to understand AR, where studies have led to the identification of molecular targets of all major anthelmintic drug classes. Here, we used genetically diverse C. elegans strains to perform dose-response analyses across 26 anthelmintic drugs that represent the three major anthelmintic drug classes (benzimidazoles, macrocyclic lactones, and nicotinic acetylcholine receptor agonists) in addition to seven other anthelmintic classes. First, we found that C. elegans strains displayed similar anthelmintic responses within drug classes and significant variation across drug classes. Next, we compared the effective concentration estimates to induce a 10% maximal response (EC10) and slope estimates of each dose-response curve of each strain to the laboratory reference strain, which enabled the identification of anthelmintics with population-wide differences to understand how genetics contribute to AR. Because genetically diverse strains displayed differential susceptibilities within and across anthelmintics, we show that C. elegans is a useful model for screening potential nematicides before applications to helminths. Third, we quantified the levels of anthelmintic response variation caused by genetic differences among individuals (heritability) to each drug and observed a significant correlation between exposure closest to the EC10 and the exposure that exhibited the most heritable responses. These results suggest drugs to prioritize in genome-wide association studies, which will enable the identification of AR genes., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Shaver et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Ancylostoma caninum is an important zoonotic gastrointestinal nematode of dogs worldwide and a close relative of human hookworms. We recently reported that racing greyhound dogs in the USA are infected with A. caninum that are commonly resistant to multiple anthelmintics. Benzimidazole resistance in A. caninum in greyhounds was associated with a high frequency of the canonical F167Y(TTC>TAC) isotype-1 β-tubulin mutation. In this work, we show that benzimidazole resistance is remarkably widespread in A. caninum from domestic dogs across the USA. First, we identified and showed the functional significance of a novel benzimidazole isotype-1 β-tubulin resistance mutation, Q134H(CAA>CAT). Several benzimidazole resistant A. caninum isolates from greyhounds with a low frequency of the F167Y(TTC>TAC) mutation had a high frequency of a Q134H(CAA>CAT) mutation not previously reported from any eukaryotic pathogen in the field. Structural modeling predicted that the Q134 residue is directly involved in benzimidazole drug binding and that the 134H substitution would significantly reduce binding affinity. Introduction of the Q134H substitution into the C. elegans β-tubulin gene ben-1, by CRISPR-Cas9 editing, conferred similar levels of resistance as a ben-1 null allele. Deep amplicon sequencing on A. caninum eggs from 685 hookworm positive pet dog fecal samples revealed that both mutations were widespread across the USA, with prevalences of 49.7% (overall mean frequency 54.0%) and 31.1% (overall mean frequency 16.4%) for F167Y(TTC>TAC) and Q134H(CAA>CAT), respectively. Canonical codon 198 and 200 benzimidazole resistance mutations were absent. The F167Y(TTC>TAC) mutation had a significantly higher prevalence and frequency in Western USA than in other regions, which we hypothesize is due to differences in refugia. This work has important implications for companion animal parasite control and the potential emergence of drug resistance in human hookworms., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Venkatesan et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The accurate characterization of structural variation is crucial for our understanding of how large chromosomal alterations affect phenotypic differences and contribute to genome evolution. Whole-genome sequencing is a popular approach for identifying structural variants, but the accuracy of popular tools remains unclear due to the limitations of existing benchmarks. Moreover, the performance of these tools for predicting variants in non-human genomes is less certain, as most tools were developed and benchmarked using data from the human genome. To evaluate the use of long-read data for the validation of short-read structural variant calls, the agreement between predictions from a short-read ensemble learning method and long-read tools were compared using real and simulated data from Caenorhabditis elegans. The results obtained from simulated data indicate that the best performing tool is contingent on the type and size of the variant, as well as the sequencing depth of coverage. These results also highlight the need for reference datasets generated from real data that can be used as 'ground truth' in benchmarks., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2022 Lesack et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Parasitic nematode infections cause an enormous global burden to both humans and livestock. Resistance to the limited arsenal of anthelmintic drugs used to combat these infections is widespread, including benzimidazole (BZ) compounds. Previous studies using the free-living nematode Caenorhabditis elegans to model parasitic nematode resistance have shown that loss-of-function mutations in the beta-tubulin gene ben-1 confer resistance to BZ drugs. However, the mechanism of resistance and the tissue-specific susceptibility are not well known in any nematode species. To identify in which tissue(s) ben-1 function underlies BZ susceptibility, transgenic strains that express ben-1 in different tissues, including hypodermis, muscles, neurons, intestine, and ubiquitous expression were generated. High-throughput fitness assays were performed to measure and compare the quantitative responses to BZ compounds among different transgenic lines. Significant BZ susceptibility was observed in animals expressing ben-1 in neurons, comparable to expression using the ben-1 promoter. This result suggests that ben-1 function in neurons underlies susceptibility to BZ. Subsetting neuronal expression of ben-1 based on the neurotransmitter system further restricted ben-1 function in cholinergic neurons to cause BZ susceptibility. These results better inform our current understanding of the cellular mode of action of BZs and also suggest additional treatments that might potentiate the effects of BZs in neurons., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Untargeted metabolomics studies are unbiased but identifying the same feature across studies is complicated by environmental variation, batch effects, and instrument variability. Ideally, several studies that assay the same set of metabolic features would be used to select recurring features to pursue for identification. Here, we developed an anchored experimental design. This generalizable approach enabled us to integrate three genetic studies consisting of 14 test strains of Caenorhabditis elegans prior to the compound identification process . An anchor strain, PD1074, was included in every sample collection, resulting in a large set of biological replicates of a genetically identical strain that anchored each study. This enables us to estimate treatment effects within each batch and apply straightforward meta-analytic approaches to combine treatment effects across batches without the need for estimation of batch effects and complex normalization strategies. We collected 104 test samples for three genetic studies across six batches to produce five analytical datasets from two complementary technologies commonly used in untargeted metabolomics. Here, we use the model system C. elegans to demonstrate that an augmented design combined with experimental blocks and other metabolomic QC approaches can be used to anchor studies and enable comparisons of stable spectral features across time without the need for compound identification. This approach is generalizable to systems where the same genotype can be assayed in multiple environments and provides biologically relevant features for downstream compound identification efforts. All methods are included in the newest release of the publicly available SECIMTools based on the open-source Galaxy platform., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Shaver, Garcia, Gouveia, Morse, Liu, Asef, Borges, Leach, Andersen, Amster, Fernández, Edison and McIntyre.)

In Caenorhabditis elegans, many genes involved in the formation of the cuticle are also known to influence body size and shape. We assessed post-embryonic growth of both long and short C. elegans body size mutants from the L1 to L4 stage. We found similar developmental trajectories of N2 and lon-3 animals. By contrast, we observed overall decreases in body length and increases in body width of tested dpy mutants compared to N2, consistent with the Dpy phenotype. We further show that the dynamics of animal shape in the mutant strains are consistent with a previously proposed "Stretcher" growth model., (Copyright: © 2022 by the authors.)

Short tandem repeats (STRs) represent an important class of genetic variation that can contribute to phenotypic differences. Although millions of single nucleotide variants (SNVs) and short indels have been identified among wild Caenorhabditis elegans strains, the natural diversity in STRs remains unknown. Here, we characterized the distribution of 31,991 STRs with motif lengths of 1-6 bp in the reference genome of C. elegans Of these STRs, 27,667 harbored polymorphisms across 540 wild strains and only 9691 polymorphic STRs (pSTRs) had complete genotype data for more than 90% of the strains. Compared with the reference genome, the pSTRs showed more contraction than expansion. We found that STRs with different motif lengths were enriched in different genomic features, among which coding regions showed the lowest STR diversity and constrained STR mutations. STR diversity also showed similar genetic divergence and selection signatures among wild strains as in previous studies using SNVs. We further identified STR variation in two mutation accumulation line panels that were derived from two wild strains and found background-dependent and fitness-dependent STR mutations. We also performed the first genome-wide association analyses between natural variation in STRs and organismal phenotypic variation among wild C. elegans strains. Overall, our results delineate the first large-scale characterization of STR variation in wild C. elegans strains and highlight the effects of selection on STR mutations., (© 2022 Zhang et al.; Published by Cold Spring Harbor Laboratory Press.)

Growth rate and body size are complex traits that contribute to the fitness of organisms. The identification of loci that underlie differences in these traits provides insights into the genetic contributions to development. Leveraging Caenorhabditis elegans as a tractable metazoan model for quantitative genetics, we can identify genomic regions that underlie differences in growth. We measured postembryonic growth of the laboratory-adapted wild-type strain (N2) and a wild strain from Hawaii (CB4856) and found differences in body size. Using linkage mapping, we identified three distinct quantitative trait loci (QTL) on chromosomes IV, V, and X that are associated with variation in body growth. We further examined these growth-associated quantitative trait loci using chromosome substitution strains and near-isogenic lines and validated the chromosome X quantitative trait loci. In addition, we generated a list of candidate genes for the chromosome X quantitative trait loci. These genes could potentially contribute to differences in animal growth and should be evaluated in subsequent studies. Our work reveals the genetic architecture underlying animal growth variation and highlights the genetic complexity of growth in Caenorhabditis elegans natural populations., (© The Author(s) 2022. Published by Oxford University Press on behalf of Genetics Society of America.)

Individuals can exhibit differences in metabolism that are caused by the interplay of genetic background, nutritional input, microbiota and other environmental factors 1-4 . It is difficult to connect differences in metabolism to genomic variation and derive underlying molecular mechanisms in humans, owing to differences in diet and lifestyle, among others. Here we use the nematode Caenorhabditis elegans as a model to study inter-individual variation in metabolism. By comparing three wild strains and the commonly used N2 laboratory strain, we find differences in the abundances of both known metabolites and those that have not to our knowledge been previously described. The latter metabolites include conjugates between 3-hydroxypropionate (3HP) and several amino acids (3HP-AAs), which are much higher in abundance in one of the wild strains. 3HP is an intermediate in the propionate shunt pathway, which is activated when flux through the canonical, vitamin-B 12 -dependent propionate breakdown pathway is perturbed 5 . We show that increased accumulation of 3HP-AAs is caused by genetic variation in HPHD-1, for which 3HP is a substrate. Our results suggest that the production of 3HP-AAs represents a 'shunt-within-a-shunt' pathway to accommodate a reduction-of-function allele in hphd-1. This study provides a step towards the development of metabolic network models that capture individual-specific differences of metabolism and more closely represent the diversity that is found in entire species., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

Phenotypic variation in organism-level traits has been studied in Caenorhabditis elegans wild strains, but the impacts of differences in gene expression and the underlying regulatory mechanisms are largely unknown. Here, we use natural variation in gene expression to connect genetic variants to differences in organismal-level traits, including drug and toxicant responses. We perform transcriptomic analyses on 207 genetically distinct C. elegans wild strains to study natural regulatory variation of gene expression. Using this massive dataset, we perform genome-wide association mappings to investigate the genetic basis underlying gene expression variation and reveal complex genetic architectures. We find a large collection of hotspots enriched for expression quantitative trait loci across the genome. We further use mediation analysis to understand how gene expression variation could underlie organism-level phenotypic variation for a variety of complex traits. These results reveal the natural diversity in gene expression and possible regulatory mechanisms in this keystone model organism, highlighting the promise of using gene expression variation to understand how phenotypic diversity is generated., (© 2022. The Author(s).)

Growth control establishes organism size, requiring mechanisms to sense and adjust growth during development. Studies of single cells revealed that size homeostasis uses distinct control methods. In multicellular organisms, mechanisms that regulate single cell growth must integrate control across organs and tissues during development to generate adult size and shape. We leveraged the roundworm Caenorhabditis elegans as a scalable and tractable model to collect precise growth measurements of thousands of individuals, measure feeding behavior, and quantify changes in animal size and shape during a densely sampled developmental time course. As animals transitioned from one developmental stage to the next, we observed changes in body aspect ratio while body volume remained constant. Then, we modeled a physical mechanism by which constraints on cuticle stretch could cause changes in C. elegans body shape. The model-predicted shape changes are consistent with those observed in the data. Theoretically, cuticle stretch could be sensed by the animal to initiate larval-stage transitions, providing a means for physical constraints to influence developmental timing and growth rate in C. elegans., (Copyright © 2022 Elsevier B.V. All rights reserved.)

Chemotaxis assays are used extensively to study behavioral responses of Caenorhabditis nematodes to environmental cues. These assays result in a chemotaxis index (CI) that denotes the behavioral response of a population of nematodes to a particular compound and can range from 1 (maximum attraction) to -1 (maximum avoidance). Traditional chemotaxis assays have low throughput because researchers must manually setup experimental populations and score CIs. Here, we describe an automated methodology that increases throughput by using liquid-handling robots to setup experimental populations and a custom image analysis package, ct, to automate the scoring of CIs from plate images., (Copyright: © 2022 by the authors.)

Caenorhabditis elegans is one of the major model organisms in biology, but only recently have researchers focused on its natural ecology. The relative sparsity of information about C. elegans in its natural context comes from the challenges involved in the identification of the small nematode in nature. Despite these challenges, an increasing focus on the ecology of C. elegans has caused a wealth of new information regarding its life outside of the laboratory. The intensified search for C. elegans in nature has contributed to the discovery of many new Caenorhabditis species and revealed that congeneric nematodes frequently cohabitate in the wild, where they feed on microbial blooms associated with rotting plant material. The identification of new species has also revealed that the androdioecious mating system of males and self-fertilizing hermaphrodites has evolved three times independently within Caenorhabditis. The other two selfing species, C. briggsae and C. tropicalis, share the experimental advantages of C. elegans and have enabled comparative studies into the mechanistic basis of important traits, including self-fertilization. Despite these advances, much remains to be learned about the ecology and natural diversity of these important species. For example, we still lack functional information for many of their genes, which might only be attained through an understanding of their natural ecology. To facilitate ecological research of selfing Caenorhabditis nematodes, we developed a highly scalable method to collect nematodes from the wild. Our method makes use of mobile data collection platforms, cloud-based databases, and the R software environment to enhance researchers' ability to collect nematodes from the wild, record associated ecological data, and identify wild nematodes using molecular barcodes.

Over the last 20 years, studies of Caenorhabditis elegans natural diversity have demonstrated the power of quantitative genetic approaches to reveal the evolutionary, ecological, and genetic factors that shape traits. These studies complement the use of the laboratory-adapted strain N2 and enable additional discoveries not possible using only one genetic background. In this chapter, we describe how to perform quantitative genetic studies in Caenorhabditis, with an emphasis on C. elegans. These approaches use correlations between genotype and phenotype across populations of genetically diverse individuals to discover the genetic causes of phenotypic variation. We present methods that use linkage, near-isogenic lines, association, and bulk-segregant mapping, and we describe the advantages and disadvantages of each approach. The power of C. elegans quantitative genetic mapping is best shown in the ability to connect phenotypic differences to specific genes and variants. We will present methods to narrow genomic regions to candidate genes and then tests to identify the gene or variant involved in a quantitative trait. The same features that make C. elegans a preeminent experimental model animal contribute to its exceptional value as a tool to understand natural phenotypic variation., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America.)

Infections by parasitic nematodes cause large health and economic burdens worldwide. We use anthelmintic drugs to reduce these infections. However, resistance to anthelmintic drugs is extremely common and increasing worldwide. It is essential to understand the mechanisms of resistance to slow its spread. Recently, four new parasitic nematode beta-tubulin alleles have been identified in benzimidazole (BZ) resistant parasite populations: E198I, E198K, E198T, and E198stop. These alleles have not been tested for the ability to confer resistance or for any effects that they might have on organismal fitness. We introduced these four new alleles into the sensitive C. elegans laboratory-adapted N2 strain and exposed these genome-edited strains to both albendazole and fenbendazole. We found that all four alleles conferred resistance to both BZ drugs. Additionally, we tested for fitness consequences in both control and albendazole conditions over seven generations in competitive fitness assays. We found that none of the edited alleles had deleterious effects on fitness in control conditions and that all four alleles conferred strong and equivalent fitness benefits in BZ drug conditions. Because it is unknown if previously validated alleles confer a dominant or recessive BZ resistance phenotype, we tested the phenotypes caused by five of these alleles and found that none of them conferred a dominant BZ resistance phenotype. Accurate measurements of resistance, fitness effects, and dominance caused by the resistance alleles allow for the generation of better models of population dynamics and facilitate control practices that maximize the efficacy of this critical anthelmintic drug class., (Copyright © 2021 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Life history traits underlie the fitness of organisms and are under strong natural selection. A new mutation that positively impacts a life history trait will likely increase in frequency and become fixed in a population (e.g., a selective sweep). The identification of the beneficial alleles that underlie selective sweeps provides insights into the mechanisms that occurred during the evolution of a species. In the global population of Caenorhabditis elegans, we previously identified selective sweeps that have drastically reduced chromosomal-scale genetic diversity in the species. Here, we measured the fecundity of 121 wild C. elegans strains, including many recently isolated divergent strains from the Hawaiian islands and found that strains with larger swept genomic regions have significantly higher fecundity than strains without evidence of the recent selective sweeps. We used genome-wide association (GWA) mapping to identify three quantitative trait loci (QTL) underlying the fecundity variation. In addition, we mapped previous fecundity data from wild C. elegans strains and C. elegans recombinant inbred advanced intercross lines that were grown in various conditions and detected eight QTL using GWA and linkage mappings. These QTL show the genetic complexity of fecundity across this species. Moreover, the haplotype structure in each GWA QTL region revealed correlations with recent selective sweeps in the C. elegans population. North American and European strains had significantly higher fecundity than most strains from Hawaii, a hypothesized origin of the C. elegans species, suggesting that beneficial alleles that caused increased fecundity could underlie the selective sweeps during the worldwide expansion of C. elegans., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America.)

Treatment of parasitic nematode infections depends primarily on the use of anthelmintics. However, this drug arsenal is limited, and resistance against most anthelmintics is widespread. Emodepside is a new anthelmintic drug effective against gastrointestinal and filarial nematodes. Nematodes that are resistant to other anthelmintic drug classes are susceptible to emodepside, indicating that the emodepside mode of action is distinct from previous anthelmintics. The laboratory-adapted Caenorhabditis elegans strain N2 is sensitive to emodepside, and genetic selection and in vitro experiments implicated slo-1, a large K + conductance (BK) channel gene, in emodepside mode of action. In an effort to understand how natural populations will respond to emodepside, we measured brood sizes and developmental rates of wild C. elegans strains after exposure to the drug and found natural variation across the species. Some of the observed variation in C. elegans emodepside responses correlates with amino acid substitutions in slo-1, but genetic mechanisms other than slo-1 coding variants likely underlie emodepside resistance in wild C. elegans strains. Additionally, the assayed strains have higher offspring production in low concentrations of emodepside (a hormetic effect). We find that natural variation affects emodepside sensitivity, supporting the suitability of C. elegans as a model system to study emodepside responses across natural nematode populations., (Copyright © 2021 The Authors. Published by Elsevier Ltd.. All rights reserved.)

The use of quality control samples in metabolomics ensures data quality, reproducibility, and comparability between studies, analytical platforms, and laboratories. Long-term, stable, and sustainable reference materials (RMs) are a critical component of the quality assurance/quality control (QA/QC) system; however, the limited selection of currently available matrix-matched RMs reduces their applicability for widespread use. To produce an RM in any context, for any matrix that is robust to changes over the course of time, we developed iterative batch averaging method (IBAT). To illustrate this method, we generated 11 independently grown Escherichia coli batches and made an RM over the course of 10 IBAT iterations. We measured the variance of these materials by nuclear magnetic resonance (NMR) and showed that IBAT produces a stable and sustainable RM over time. This E. coli RM was then used as a food source to produce a Caenorhabditis elegans RM for a metabolomics experiment. The metabolite extraction of this material, alongside 41 independently grown individual C. elegans samples of the same genotype, allowed us to estimate the proportion of sample variation in preanalytical steps. From the NMR data, we found that 40% of the metabolite variance is due to the metabolite extraction process and analysis and 60% is due to sample-to-sample variance. The availability of RMs in untargeted metabolomics is one of the predominant needs of the metabolomics community that reach beyond quality control practices. IBAT addresses this need by facilitating the production of biologically relevant RMs and increasing their widespread use.

To better understand the mechanism of resistance caused by putative interactions between beta-tubulin and benzimidazole compounds, we sought to purify nematode-specific beta-tubulins using heterologous expression after replacement of the single Saccharomyces cerevisiae beta-tubulin gene. However, we found that haploid yeast cells containing nematode-specific beta-tubulin genes were not viable, suggesting that nematode beta-tubulin cannot substitute for the loss of the yeast beta-tubulin gene., (Copyright: © 2021 by the authors.)

Across diverse taxa, selfing species have evolved independently from outcrossing species thousands of times. The transition from outcrossing to selfing decreases the effective population size, effective recombination rate and heterozygosity within a species. These changes lead to a reduction in genetic diversity, and therefore adaptive potential, by intensifying the effects of random genetic drift and linked selection. Within the nematode genus Caenorhabditis, selfing has evolved at least three times, and all three species, including the model organism Caenorhabditis elegans, show substantially reduced genetic diversity relative to outcrossing species. Selfing and outcrossing Caenorhabditis species are often found in the same niches, but we still do not know how selfing species with limited genetic diversity can adapt to these environments. Here, we examine the whole-genome sequences from 609 wild C. elegans strains isolated worldwide and show that genetic variation is concentrated in punctuated hyper-divergent regions that cover 20% of the C. elegans reference genome. These regions are enriched in environmental response genes that mediate sensory perception, pathogen response and xenobiotic stress response. Population genomic evidence suggests that genetic diversity in these regions has been maintained by long-term balancing selection. Using long-read genome assemblies for 15 wild strains, we show that hyper-divergent haplotypes contain unique sets of genes and show levels of divergence comparable to levels found between Caenorhabditis species that diverged millions of years ago. These results provide an example of how species can avoid the evolutionary dead end associated with selfing.

Anthelmintic drugs are the major line of defense against parasitic nematode infections, but the arsenal is limited and resistance threatens sustained efficacy of the available drugs. Discoveries of the modes of action of these drugs and mechanisms of resistance have predominantly come from studies of a related nonparasitic nematode species, Caenorhabditis elegans, and the parasitic nematode Haemonchus contortus. Here, we discuss how our understanding of anthelmintic resistance and modes of action came from the interplay of results from each of these species. We argue that this 'cycle of discovery', where results from one species inform the design of experiments in the other, can use the complementary strengths of both to understand anthelmintic modes of action and mechanisms of resistance., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

For more than four decades, the free-living nematode Caenorhabditis elegans has been extensively used in anthelmintic research. Classic genetic screens and heterologous expression in the C. elegans model enormously contributed to the identification and characterization of molecular targets of all major anthelmintic drug classes. Although these findings provided substantial insights into common anthelmintic mechanisms, a breakthrough in the treatment and control of parasitic nematodes is still not in sight. Instead, we are facing increasing evidence that the enormous diversity within the phylum Nematoda cannot be recapitulated by any single free-living or parasitic species and the development of novel broad-spectrum anthelmintics is not be a simple goal. In the present review, we summarize certain milestones and challenges of the C. elegans model with focus on drug target identification, anthelmintic drug discovery and identification of resistance mechanisms. Furthermore, we present new perspectives and strategies on how current progress in C. elegans research will support future anthelmintic research., (Copyright © 2020 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Infections by parasitic nematodes inflict a huge burden on the health of humans and livestock throughout the world. Anthelmintic drugs are the first line of defense against these infections. Unfortunately, resistance to these drugs is rampant and continues to spread. To improve treatment strategies, we must understand the genetics and molecular mechanisms that underlie resistance. Studies of the fungus Aspergillus nidulans and the free-living nematode Caenorhabditis elegans discovered that a beta-tubulin gene is mutated in benzimidazole (BZ) resistant strains. In parasitic nematode populations, three beta-tubulin alleles, F167Y, E198A, and F200Y, have long been correlated with resistance. Additionally, improvements in sequencing technologies have identified new alleles - E198V, E198L, E198K, E198I, and E198Stop - also correlated with BZ resistance. However, none of these alleles have been proven to cause resistance. To empirically demonstrate this point, we independently introduced the F167Y, E198A, and F200Y alleles as well as two of the newly identified alleles, E198V and E198L, into the BZ susceptible C. elegans N2 genetic background using the CRISPR-Cas9 system. These genome-edited strains were exposed to both albendazole and fenbendazole to quantitatively measure animal responses to BZs. We used a range of concentrations for each BZ compound to define response curves and found that all five of the alleles conferred resistance to BZ compounds equal to a loss of the entire beta-tubulin gene. These results prove that the parasite beta-tubulin alleles cause resistance. The E198V allele is found at low frequencies along with the E198L allele in natural parasite populations, suggesting that it could affect fitness. We performed competitive fitness assays and demonstrated that the E198V allele reduces animal health, supporting the hypothesis that this allele might be less fit in field populations. Overall, we present a powerful platform to quantitatively assess anthelmintic resistance and effects of specific resistance alleles on organismal fitness in the presence or absence of the drug., (Copyright © 2020 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Pleiotropy, the concept that a single gene controls multiple distinct traits, is prevalent in most organisms and has broad implications for medicine and agriculture. The identification of the molecular mechanisms underlying pleiotropy has the power to reveal previously unknown biological connections between seemingly unrelated traits. Additionally, the discovery of pleiotropic genes increases our understanding of both genetic and phenotypic complexity by characterizing novel gene functions. Quantitative trait locus (QTL) mapping has been used to identify several pleiotropic regions in many organisms. However, gene knockout studies are needed to eliminate the possibility of tightly linked, non-pleiotropic loci. Here, we use a panel of 296 recombinant inbred advanced intercross lines of Caenorhabditis elegans and a high-throughput fitness assay to identify a single large-effect QTL on the center of chromosome V associated with variation in responses to eight chemotherapeutics. We validate this QTL with near-isogenic lines and pair genome-wide gene expression data with drug response traits to perform mediation analysis, leading to the identification of a pleiotropic candidate gene, scb-1 , for some of the eight chemotherapeutics. Using deletion strains created by genome editing, we show that scb-1 , which was previously implicated in response to bleomycin, also underlies responses to other double-strand DNA break-inducing chemotherapeutics. This finding provides new evidence for the role of scb-1 in the nematode drug response and highlights the power of mediation analysis to identify causal genes., (Copyright © 2020 Evans, Andersen.)