Extreme allelic heterogeneity at a Caenorhabditis elegans beta-tubulin locus explains natural resistance to benzimidazoles

Development of a size-separation technique to isolate Caenorhabditis elegans embryos using mesh filters

The free-living nematode Caenorhabditis elegans has been routinely used to study gene functions, genetic interactions, and conserved signaling pathways. Most experiments require that the animals are synchronized to be at the same specific developmental stage. Bleach synchronization is traditionally used to obtain a population of staged embryos, but the method can have harmful effects on the embryos. The physical separation of differently sized animals is preferred but often difficult to perform because some developmental stages are the same sizes as others. Microfluidic device filters have been used as alternatives, but they are expensive and require customization to scale up the preparation of staged animals. Here, we present a protocol for the synchronization of embryos using mesh filters. Using filtration, we obtained a higher yield of embryos per plate than using the standard bleach synchronization protocol and at a scale beyond microfluidic devices. Importantly, filtration has no deleterious effects on downstream larval development assays. In conclusion, we have exploited the differences in the sizes of C. elegans developmental stages to isolate embryo cultures suitable for use in high-throughput assays.

Independent mechanisms of benzimidazole resistance across Caenorhabditis nematodes

Benzimidazoles (BZs), a widely used class of anthelmintic drugs, target beta-tubulin proteins, disrupt microtubule formation, and cause nematode death. In parasitic nematode species, mutations in beta-tubulin genes (e.g., isotype-1 beta-tubulin) are predicted to inhibit BZ binding and are associated with BZ resistance. Similarly, in the free-living nematode Caenorhabditis elegans, mutations in an isotype-1 beta-tubulin ortholog, ben-1, are the primary drivers of BZ resistance. The recurrent association of BZ resistance with beta-tubulins suggests that BZ resistance is repeatedly caused by mutations in beta-tubulin genes, an example of repeated evolution of drug resistance across nematode species. To evaluate the hypothesis of repeated evolution of BZ resistance mediated by beta-tubulin, we identified predicted resistance alleles in beta-tubulin genes across wild strains from three Caenorhabditis species: C. elegans, Caenorhabditis briggsae, and Caenorhabditis tropicalis. We hypothesized that, if these species experienced similar selective pressures, they would evolve resistance to BZs by mutations in any of three beta-tubulin genes (ben-1, tbb-1, and tbb-2). Using high-throughput development assays, we tested the association of predicted beta-tubulin alleles with BZ resistance. We found that a heterogeneous set of variants identified in C. elegans ben-1 were associated with BZ resistance. In C. briggsae, only two variants in ben-1, predicted to encode a premature stop codon (W21stop) and a missense substitution (Q134H), were associated with BZ resistance. In C. tropicalis, two missense variants were identified in ben-1, but neither was associated with BZ resistance. C. briggsae and C. tropicalis might have evolved BZ resistance by mutations in other beta-tubulin genes, but we found that variants in tbb-1 or tbb-2 in these species were not associated with BZ resistance. Our findings reveal a lack of repeated evolution of BZ resistance across the three Caenorhabditis species and highlight the importance of defining BZ resistance mechanisms outside of beta-tubulins.

Naturally occurring variation in a cytochrome P450 modifies thiabendazole responses independently of beta-tubulin

Widespread anthelmintic resistance has complicated the management of parasitic nematodes. Resistance to the benzimidazole (BZ) drug class is nearly ubiquitous in many species and is associated with mutations in beta-tubulin genes. However, mutations in beta-tubulin alone do not fully explain all BZ resistance. We performed a genome-wide association study using a genetically diverse panel of Caenorhabditis elegans strains to identify loci that contribute to resistance to the BZ drug thiabendazole (TBZ). We identified a quantitative trait locus (QTL) on chromosome V independent of all beta-tubulin genes and overlapping with two promising candidate genes, the cytochrome P450 gene cyp-35D1 and the nuclear hormone receptor nhr-176. Both genes were previously demonstrated to play a role in TBZ metabolism. NHR-176 binds TBZ and induces the expression of CYP-35D1, which metabolizes TBZ. We generated single gene deletions of cyp-35D1 and nhr-176 and found that both genes play a role in TBZ response. A predicted high-impact lysine-to-glutamate substitution at position 267 (K267E) in CYP-35D1 was identified in a sensitive strain, and reciprocal allele replacement strains in different genetic backgrounds were used to show that the lysine allele conferred increased TBZ resistance. Using competitive fitness assays, we found that neither allele was deleterious, but the lysine allele was selected in the presence of TBZ. Additionally, we found that the lysine allele significantly increased the rate of TBZ metabolism compared to the glutamate allele. Moreover, yeast expression assays showed that the lysine version of CYP-35D1 had twice the enzymatic activity of the glutamate allele. To connect our results to parasitic nematodes, we analyzed four Haemonchus contortus cytochrome P450 orthologs but did not find variation at the 267 position in fenbendazole-resistant populations. Overall, we confirmed that variation in this cytochrome P450 gene is the first locus independent of beta-tubulin to play a role in BZ resistance.

Integrating metabolomics into the diagnosis and investigation of anthelmintic resistance

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.

Quantitative tests of albendazole resistance in Caenorhabditis elegans beta-tubulin mutants

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.

Identification of a family of species-selective complex I inhibitors as potential anthelmintics

Soil-transmitted helminths (STHs) are major pathogens infecting over a billion people. There are few classes of anthelmintics and there is an urgent need for new drugs. Many STHs use an unusual form of anaerobic metabolism to survive the hypoxic conditions of the host gut. This requires rhodoquinone (RQ), a quinone electron carrier. RQ is not made or used by vertebrate hosts making it an excellent therapeutic target. Here we screen 480 structural families of natural products to find compounds that kill Caenorhabditis elegans specifically when they require RQ-dependent metabolism. We identify several classes of compounds including a family of species-selective inhibitors of mitochondrial respiratory complex I. These identified complex I inhibitors have a benzimidazole core and we determine key structural requirements for activity by screening 1,280 related compounds. Finally, we show several of these compounds kill adult STHs. We suggest these species-selective complex I inhibitors are potential anthelmintics.

Quantifying the fitness effects of resistance alleles with and without anthelmintic selection pressure using Caenorhabditis elegans

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.

Quantitative tests of albendazole resistance in beta-tubulin mutants

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 examine the roles of five C. elegans beta-tubulin genes (tbb-1, mec-7, tbb-4, ben-1, and tbb-6) to identify the role each gene plays in BZ response. 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 to test beta-tubulin redundancy in BZ response. We found that only the individual loss of ben-1 conferred a substantial level of BZ resistance, although the loss of tbb-1 was found to confer a small benefit in the presence of albendazole (ABZ). The loss of ben-1 was found to confer an almost complete rescue of animal development in the presence of 30 μM ABZ, likely explaining why no additive effects caused by the loss of a second beta-tubulin were observed. We demonstrate that ben-1 is the only beta-tubulin gene in C. elegans where loss confers substantial BZ resistance.

Altered larval activation response associated with multidrug resistance in the canine hookworm Ancylostoma caninum

Parasitic gastrointestinal nematodes pose significant health risks to humans, livestock, and companion animals, and their control relies heavily on the use of anthelmintic drugs. Overuse of these drugs has led to the emergence of resistant nematode populations. Herein, a naturally occurring isolate (referred to as BCR) of the dog hookworm, Ancylostoma caninum, that is resistant to 3 major classes of anthelmintics is characterized. Various drug assays were used to determine the resistance of BCR to thiabendazole, ivermectin, moxidectin and pyrantel pamoate. When compared to a drug-susceptible isolate of A. caninum, BCR was shown to be significantly resistant to all 4 of the drugs tested. Multiple single nucleotide polymorphisms have been shown to impart benzimidazole resistance, including the F167Y mutation in the β-tubulin isotype 1 gene, which was confirmed to be present in BCR through molecular analysis. The frequency of the resistant allele in BCR was 76.3% following its first passage in the lab, which represented an increase from approximately 50% in the founding hookworm population. A second, recently described mutation in codon 134 (Q134H) was also detected at lower frequency in the BCR population. Additionally, BCR exhibits an altered larval activation phenotype compared to the susceptible isolate, suggesting differences in the signalling pathways involved in the activation process which may be associated with resistance. Further characterization of this isolate will provide insights into the mechanisms of resistance to macrocyclic lactones and tetrahydropyrimidine anthelmintics.

Getting around the roundworms: Identifying knowledge gaps and research priorities for the ascarids

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.

Quantifying the fitness effects of resistance alleles with and without anthelmintic selection pressure using Caenorhabditis elegans

Albendazole and ivermectin 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 channel (GluCl) genes, but it is unknown whether these 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 loss of two GluCl genes (avr-14 and avr-15) or loss of three GluCl genes (avr-14, avr-15, and glc-1). The fecundity assays revealed that loss of ben-1 did not provide any fitness benefit in albendazole 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.

CaeNDR, the Caenorhabditis Natural Diversity Resource

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.

Understanding anthelmintic resistance in livestock using "omics" approaches

Widespread and improper use of various anthelmintics, genetic, and epidemiological factors has resulted in anthelmintic-resistant (AR) helminth populations in livestock. This is currently quite common globally in different livestock animals including sheep, goats, and cattle to gastrointestinal nematode (GIN) infections. Therefore, the mechanisms underlying AR in parasitic worm species have been the subject of ample research to tackle this challenge. Current and emerging technologies in the disciplines of genomics, transcriptomics, metabolomics, and proteomics in livestock species have advanced the understanding of the intricate molecular AR mechanisms in many major parasites. The technologies have improved the identification of possible biomarkers of resistant parasites, the ability to find actual causative genes, regulatory networks, and pathways of parasites governing the AR development including the dynamics of helminth infection and host-parasite infections. In this review, various "omics"-driven technologies including genome scan, candidate gene, quantitative trait loci, transcriptomic, proteomic, and metabolomic approaches have been described to understand AR of parasites of veterinary importance. Also, challenges and future prospects of these "omics" approaches are also discussed.

Screening of the Pandemic Response Box identifies anti-microsporidia compounds

Microsporidia are fungal obligate intracellular pathogens, which infect most animals and cause microsporidiosis. Despite the serious threat that microsporidia pose to humans and agricultural animals, few drugs are available for the treatment and control of microsporidia. To identify novel inhibitors, we took advantage of the model organism Caenorhabditis elegans infected with its natural microsporidian Nematocida parisii. We used this system to screen the Pandemic Response Box, a collection of 400 diverse compounds with known antimicrobial activity. After testing these compounds in a 96-well format at high (100 μM) and low (40 μM) concentrations, we identified four inhibitors that restored the ability of C. elegans to produce progeny in the presence of N. parisii. All four compounds reduced the pathogen load of both N. parisii and Pancytospora epiphaga, a C. elegans-infecting microsporidia related to human-infecting species. One of these compounds, a known inhibitor of a viral protease, MMV1006203, inhibited invasion and prevented the firing of spores. A bis-indole derivative, MMV1593539, decreased spore viability. An albendazole analog, MMV1782387, inhibited proliferation of N. parisii. We tested albendazole as well as 5 other analogs and observed that MMV1782387 was amongst the strongest inhibitors of N. parisii and displayed the least host toxicity. Our study further demonstrates the effectiveness of the C. elegans-N. parisii system for discovering microsporidia inhibitors and the compounds we identified provide potential scaffolds for anti-microsporidia drug development.

Praziquantel inhibits Caenorhabditis elegans development and species-wide differences might be cct-8-dependent

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.

Beyond the reference: gene expression variation and transcriptional response to RNA interference in Caenorhabditis elegans

Though natural systems harbor genetic and phenotypic variation, research in model organisms is often restricted to a reference strain. Focusing on a reference strain yields a great depth of knowledge but potentially at the cost of breadth of understanding. Furthermore, tools developed in the reference context may introduce bias when applied to other strains, posing challenges to defining the scope of variation within model systems. Here, we evaluate how genetic differences among 5 wild Caenorhabditis elegans strains affect gene expression and its quantification, in general and after induction of the RNA interference (RNAi) response. Across strains, 34% of genes were differentially expressed in the control condition, including 411 genes that were not expressed at all in at least 1 strain; 49 of these were unexpressed in reference strain N2. Reference genome mapping bias caused limited concern: despite hyperdiverse hotspots throughout the genome, 92% of variably expressed genes were robust to mapping issues. The transcriptional response to RNAi was highly strain- and target-gene-specific and did not correlate with RNAi efficiency, as the 2 RNAi-insensitive strains showed more differentially expressed genes following RNAi treatment than the RNAi-sensitive reference strain. We conclude that gene expression, generally and in response to RNAi, differs across C. elegans strains such that the choice of strain may meaningfully influence scientific inferences. Finally, we introduce a resource for querying gene expression variation in this dataset at https://wildworm.biosci.gatech.edu/rnai/.

Interplay Between Polymorphic Short Tandem Repeats and Gene Expression Variation in Caenorhabditis elegans

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.

Variation in anthelmintic responses are driven by genetic differences among diverse C. elegans wild strains

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.

Beyond the reference: gene expression variation and transcriptional response to RNAi in C. elegans

A universal feature of living systems is that natural variation in genotype underpins variation in phenotype. Yet, research in model organisms is often constrained to a single genetic background, the reference strain. Further, genomic studies that do evaluate wild strains typically rely on the reference strain genome for read alignment, leading to the possibility of biased inferences based on incomplete or inaccurate mapping; the extent of reference bias can be difficult to quantify. As an intermediary between genome and organismal traits, gene expression is well positioned to describe natural variability across genotypes generally and in the context of environmental responses, which can represent complex adaptive phenotypes. C. elegans sits at the forefront of investigation into small-RNA gene regulatory mechanisms, or RNA interference (RNAi), and wild strains exhibit natural variation in RNAi competency following environmental triggers. Here, we examine how genetic differences among five wild strains affect the C. elegans transcriptome in general and after inducing RNAi responses to two germline target genes. Approximately 34% of genes were differentially expressed across strains; 411 genes were not expressed at all in at least one strain despite robust expression in others, including 49 genes not expressed in reference strain N2. Despite the presence of hyper-diverse hotspots throughout the C. elegans genome, reference mapping bias was of limited concern: over 92% of variably expressed genes were robust to mapping issues. Overall, the transcriptional response to RNAi was strongly strain-specific and highly specific to the target gene, and the laboratory strain N2 was not representative of the other strains. Moreover, the transcriptional response to RNAi was not correlated with RNAi phenotypic penetrance; the two germline RNAi incompetent strains exhibited substantial differential gene expression following RNAi treatment, indicating an RNAi response despite failure to reduce expression of the target gene. We conclude that gene expression, both generally and in response to RNAi, differs across C. elegans strains such that choice of strain may meaningfully influence scientific conclusions. To provide a public, easily accessible resource for querying gene expression variation in this dataset, we introduce an interactive website at https://wildworm.biosci.gatech.edu/rnai/ .

Molecular evidence of widespread benzimidazole drug resistance in Ancylostoma caninum from domestic dogs throughout the USA and discovery of a novel β-tubulin benzimidazole resistance mutation

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.

Functions of the tubulin code in the C. elegans nervous system

Due to their elongated and polarized morphology, neurons rely on the microtubule (MT) cytoskeleton for their shape, as well as for efficient intracellular transport that maintains neuronal function, survival, and connectivity. Although all MTs are constructed from α- and β-tubulins that are highly conserved throughout eukaryotes, different MT networks within neurons exhibit different dynamics and functions. For example, molecular motors must be able to differentially recognize the axonal and dendritic MTs to deliver appropriate cargos to sensory endings and synaptic regions. The Tubulin Code hypothesis proposes that MTs can be specialized in form and function by chemical differences in their composition by inclusion of different α- and β-tubulins into the MT lattice, as well as differences in post-translational enzymatic modifications. The chemical differences encode information that allow MTs to regulate interactions with various microtubule-based molecular motors such as kinesins and dyneins as well as with structural microtubule-associated proteins (MAPs), which can, in turn, modify the function or stability of MTs. Here, we review studies involving C. elegans, a model organism with a relatively simple nervous system that is amenable to genetic analysis, that have contributed to our understanding of how the Tubulin Code can specialize neuronal MT networks to establish differences in neuronal morphology and function. Such studies have revealed molecules and mechanisms that are conserved in vertebrates and have the potential to inform our understanding of neurological diseases involving defects in the cytoskeleton and intracellular transport.

Benzimidazoles cause lethality by inhibiting the function of Caenorhabditis elegans neuronal beta-tubulin

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.

Identification of Novel Genes and Proteoforms in Angiostrongylus costaricensis through a Proteogenomic Approach

RNA sequencing (RNA-Seq) and mass-spectrometry-based proteomics data are often integrated in proteogenomic studies to assist in the prediction of eukaryote genome features, such as genes, splicing, single-nucleotide (SNVs), and single-amino-acid variants (SAAVs). Most genomes of parasite nematodes are draft versions that lack transcript- and protein-level information and whose gene annotations rely only on computational predictions. Angiostrongylus costaricensis is a roundworm species that causes an intestinal inflammatory disease, known as abdominal angiostrongyliasis (AA). Currently, there is no drug available that acts directly on this parasite, mostly due to the sparse understanding of its molecular characteristics. The available genome of A. costaricensis, specific to the Costa Rica strain, is a draft version that is not supported by transcript- or protein-level evidence. This study used RNA-Seq and MS/MS data to perform an in-depth annotation of the A. costaricensis genome. Our prediction improved the reference annotation with (a) novel coding and non-coding genes; (b) pieces of evidence of alternative splicing generating new proteoforms; and (c) a list of SNVs between the Brazilian (Crissiumal) and the Costa Rica strain. To the best of our knowledge, this is the first time that a multi-omics approach has been used to improve the genome annotation of A. costaricensis. We hope this improved genome annotation can assist in the future development of drugs, kits, and vaccines to treat, diagnose, and prevent AA caused by either the Brazil strain (Crissiumal) or the Costa Rica strain.

Genomic landscape of drug response reveals mediators of anthelmintic resistance

Like other pathogens, parasitic helminths can rapidly evolve resistance to drug treatment. Understanding the genetic basis of anthelmintic drug resistance in parasitic nematodes is key to tracking its spread and improving the efficacy and sustainability of parasite control. Here, we use an in vivo genetic cross between drug-susceptible and multi-drug-resistant strains of Haemonchus contortus in a natural host-parasite system to simultaneously map resistance loci for the three major classes of anthelmintics. This approach identifies new alleles for resistance to benzimidazoles and levamisole and implicates the transcription factor cky-1 in ivermectin resistance. This gene is within a locus under selection in ivermectin-resistant populations worldwide; expression analyses and functional validation using knockdown experiments support that cky-1 is associated with ivermectin survival. Our work demonstrates the feasibility of high-resolution forward genetics in a parasitic nematode and identifies variants for the development of molecular diagnostics to combat drug resistance in the field.

Genomic signatures of selection associated with benzimidazole drug treatments in Haemonchus contortus field populations

Genome-wide methods offer a powerful approach to detect signatures of drug selection. However, limited availability of suitable reference genomes and the difficulty of obtaining field populations with well-defined, distinct drug treatment histories mean there is little information on the signatures of selection in parasitic nematodes and on how best to detect them. This study addresses these knowledge gaps by using field populations of Haemonchus contortus with well-defined benzimidazole treatment histories, leveraging a recently completed chromosomal-scale reference genome assembly. We generated a panel of 49,393 genomic markers to genotype 20 individual adult worms from each of four H. contortus populations: two from closed sheep flocks with an approximate 20 year history of frequent benzimidazole treatment, and two populations with a history of little or no treatment. Sampling occurred in the same geographical region to limit genetic differentiation and maximise the detection sensitivity. A clear signature of selection was detected on chromosome I, centred on the isotype-1 β-tubulin gene. Two additional, but weaker, signatures of selection were detected; one near the middle of chromosome I spanning 3.75 Mbp and 259 annotated genes, and one on chromosome II spanning a region of 3.3 Mbp and 206 annotated genes, including the isotype-2 β-tubulin locus. We also assessed how sensitivity was impacted by sequencing depth, worm number, and pooled versus individual worm sequence data. This study provides the first known direct genome-wide evidence for any parasitic nematode, that the isotype-1 β-tubulin gene is quantitatively the single most important benzimidazole resistance locus. It also identified two additional genomic regions that likely contain benzimidazole resistance loci of secondary importance. This study provides an experimental framework to maximise the power of genome-wide approaches to detect signatures of selection driven by anthelmintic drug treatments in field populations of parasitic nematodes.

Transgenic Expression of Haemonchus contortus Cytochrome P450 Hco-cyp-13A11 Decreases Susceptibility to Particular but Not All Macrocyclic Lactones in the Model Organism Caenorhabditis elegans

The number of reported macrocyclic lactones (ML) resistance cases across all livestock hosts is steadily increasing. Different studies in the parasitic nematode Haemonchus contortus assume the participation of cytochrome P450s (Cyps) enzymes in ML resistance. Still, functional data about their individual contribution to resistance or substrate specificity is missing. Via microinjection, transgenic Caenorhabditis elegans expressing HCON_00141052 (transgene-Hco-cyp-13A11) from extrachromosomal arrays were generated. After 24 h of exposure to different concentrations of ivermectin (IVM), ivermectin aglycone (IVMa), selamectin (SEL), doramectin (DRM), eprinomectin (EPR), and moxidectin (MOX), motility assays were performed to determine the impact of the H. contortus Cyp to the susceptibility of the worms against each ML. While transgene-Hco-cyp-13A11 significantly decreased susceptibility to IVM (four-fold), IVMa (2-fold), and SEL (3-fold), a slight effect for DRM and no effect for MOX, and EPR was observed. This substrate specificity of Hco-cyp-13A11 could not be explained by molecular modeling and docking studies. Hco-Cyp-13A11 molecular models were obtained for alleles from isolates with different resistance statuses. Although 14 amino acid polymorphisms were detected, none was resistance specific. In conclusion, Hco-cyp-13A11 decreased IVM, IVMa, and SEL susceptibility to a different extent, but its potential impact on ML resistance is not driven by polymorphisms.

Evaluating the power and limitations of genome-wide association studies in Caenorhabditis elegans

Quantitative genetics in Caenorhabditis elegans seeks to identify naturally segregating genetic variants that underlie complex traits. Genome-wide association studies scan the genome for individual genetic variants that are significantly correlated with phenotypic variation in a population, or quantitative trait loci. Genome-wide association studies are a popular choice for quantitative genetic analyses because the quantitative trait loci that are discovered segregate in natural populations. Despite numerous successful mapping experiments, the empirical performance of genome-wide association study has not, to date, been formally evaluated in C. elegans. We developed an open-source genome-wide association study pipeline called NemaScan and used a simulation-based approach to provide benchmarks of mapping performance in collections of wild C. elegans strains. Simulated trait heritability and complexity determined the spectrum of quantitative trait loci detected by genome-wide association studies. Power to detect smaller-effect quantitative trait loci increased with the number of strains sampled from the C. elegans Natural Diversity Resource. Population structure was a major driver of variation in mapping performance, with populations shaped by recent selection exhibiting significantly lower false discovery rates than populations composed of more divergent strains. We also recapitulated previous genome-wide association studies of experimentally validated quantitative trait variants. Our simulation-based evaluation of performance provides the community with critical context to pursue quantitative genetic studies using the C. elegans Natural Diversity Resource to elucidate the genetic basis of complex traits in C. elegans natural populations.

C. elegans as a model for inter-individual variation in metabolism

Individuals can exhibit differences in metabolism that are caused by the interplay of genetic background, nutritional input, microbiota and other environmental factors1-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-B12-dependent propionate breakdown pathway is perturbed5. 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.

Interactions of C. elegans β-tubulins with the microtubule inhibitor and anthelmintic drug albendazole

Parasitic nematodes are major human and agricultural pests, and benzimidazoles are amongst the most important broad spectrum anthelmintic drug class used for their control. Benzimidazole resistance is now widespread in many species of parasitic nematodes in livestock globally and an emerging concern for the sustainable control of human soil transmitted helminths. β-tubulin is the major benzimidazole target, although other genes may influence resistance. Among the six C. elegans β-tubulin genes, loss of ben-1 causes resistance without other apparent defects. Here, we explored the genetics of C. elegans β-tubulin genes in relation to the response to the benzimidazole derivative albendazole. The most highly expressed β-tubulin isotypes, encoded by tbb-1 and tbb-2, were known to be redundant with each other for viability, and their products are predicted not to bind benzimidazoles. We found that tbb-2 mutants, and to a lesser extent tbb-1 mutants, were hypersensitive to albendazole. The double mutant tbb-2 ben-1 is uncoordinated and short, resembling the wild type exposed to albendazole, but the tbb-1 ben-1 double mutant did not show the same phenotypes. These results suggest that tbb-2 is a modifier of ABZ sensitivity. To better understand how BEN-1 mutates to cause benzimidazole resistance, we isolated mutants resistant to albendazole and found that 15 of 16 mutations occurred in the ben-1 coding region. Mutations ranged from likely nulls to hypomorphs, and several corresponded to residues that cause resistance in other organisms. Null alleles of ben-1 are albendazole-resistant and BEN-1 shows high sequence identity with tubulins from other organisms, suggesting that many amino acid changes could cause resistance. However, our results suggest that missense mutations conferring resistance are not evenly distributed across all possible conserved sites. Independent of their roles in benzimidazole resistance, tbb-1 and tbb-2 may have specialized functions as null mutants of tbb-1 or tbb-2 were cold or heat sensitive, respectively.

The impact of species-wide gene expression variation on Caenorhabditis elegans complex traits

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.

Absence of Polymorphisms in Codons 167, 198 and 200 of All Seven β-tubulin Isotypes of Benzimidazole Susceptible and Resistant Parascaris spp. Specimens from Australia

Benzimidazoles resistance is widespread in strongyle parasitic nematodes and associated with polym orphisms in the codons 167, 198 and 200 of isotype 1 β-tubulin (tbb-1). In ascarids, benzimidazole (BZ) resistance has rarely been reported and in none of these cases were any of these polymorphisms detected. Here, available genome and transcriptome data from WormBase ParaSite were used to compare the complete β-tubulin reservoirs of Parascaris univalens, Ascaris suum and Ascaris lumbricoides. Adult Parascaris spp. specimens collected in Australia from horses after BZ treatment (susceptible, n = 13) or surviving BZ treatment and collected after ivermectin treatment (resistant, n = 10) were genotyped regarding codons 167, 198 and 200 using Sanger sequencing. Phylogenetic analyses clearly showed that there are no one-to-one ascarid orthologs of strongyle tbb-1 genes. In the reference genomes, as well as phenotypically susceptible and resistant Parascaris spp. from Australia, six out of seven β-tubulin genes showed a BZ-susceptible genotype (F167, E198, F200). The only exception were the testis-specific β-tubulin D genes from all three ascarid species that encode tyrosine at codon 200. This was observed independently of the BZ-susceptibility phenotype of Parascaris spp. These data suggest that different mechanisms lead to BZ resistance in ascarid and strongyle nematodes.

Local adaptation and spatiotemporal patterns of genetic diversity revealed by repeated sampling of Caenorhabditis elegans across the Hawaiian Islands

The nematode Caenorhabditis elegans is among the most widely studied organisms, but relatively little is known about its natural ecology. Genetic diversity is low across much of the globe but high in the Hawaiian Islands and across the Pacific Rim. To characterize the niche and genetic diversity of C. elegans on the Hawaiian Islands and to explore how genetic diversity might be influenced by local adaptation, we repeatedly sampled nematodes over a three-year period, measured various environmental parameters at each sampling site, and whole-genome sequenced the C. elegans isolates that we identified. We found that the typical Hawaiian C. elegans niche comprises moderately moist native forests at high elevations (500-1,500 m) where ambient air temperatures are cool (15-20°C). Compared to other Caenorhabditis species found on the Hawaiian Islands (e.g., Caenorhabditis briggsae and Caenorhabditis tropicalis), we found that C. elegans were enriched in native habitats. We measured levels of genetic diversity and differentiation among Hawaiian C. elegans and found evidence of seven genetically distinct groups distributed across the islands. Then, we scanned these genomes for signatures of local adaptation and identified 18 distinct regions that overlap with hyper-divergent regions, which may be maintained by balancing selection and are enriched for genes related to environmental sensing, xenobiotic detoxification, and pathogen resistance. These results provide strong evidence of local adaptation among Hawaiian C. elegans and contribute to our understanding of the forces that shape genetic diversity on the most remote volcanic archipelago in the world.

Advances in our understanding of nematode ion channels as potential anthelmintic targets

Ion channels are specialized multimeric proteins that underlie cell excitability. These channels integrate with a variety of neuromuscular and biological functions. In nematodes, the physiological behaviors including locomotion, navigation, feeding and reproduction, are regulated by these protein entities. Majority of the antinematodal chemotherapeutics target the ion channels to disrupt essential biological functions. Here, we have summarized current advances in our understanding of nematode ion channel pharmacology. We review cys-loop ligand gated ion channels (LGICs), including nicotinic acetylcholine receptors (nAChRs), acetylcholine-chloride gated ion channels (ACCs), glutamate-gated chloride channels (GluCls), and GABA (γ-aminobutyric acid) receptors, and other ionotropic receptors (transient receptor potential (TRP) channels and potassium ion channels). We have provided an update on the pharmacological properties of these channels from various nematodes. This article catalogs the differences in ion channel composition and resulting pharmacology in the phylum Nematoda. This diversity in ion channel subunit repertoire and pharmacology emphasizes the importance of pursuing species-specific drug target research. In this review, we have provided an overview of recent advances in techniques and functional assays available for screening ion channel properties and their application.

ConVarT: a search engine for matching human genetic variants with variants from non-human species

The availability of genetic variants, together with phenotypic annotations from model organisms, facilitates comparing these variants with equivalent variants in humans. However, existing databases and search tools do not make it easy to scan for equivalent variants, namely 'matching variants' (MatchVars) between humans and other organisms. Therefore, we developed an integrated search engine called ConVarT (http://www.convart.org/) for matching variants between humans, mice, and Caenorhabditis elegans. ConVarT incorporates annotations (including phenotypic and pathogenic) into variants, and these previously unexploited phenotypic MatchVars from mice and C. elegans can give clues about the functional consequence of human genetic variants. Our analysis shows that many phenotypic variants in different genes from mice and C. elegans, so far, have no counterparts in humans, and thus, can be useful resources when evaluating a relationship between a new human mutation and a disease.

Natural genetic variation as a tool for discovery in Caenorhabditis nematodes

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.

Identification of enzymes that have helminth-specific active sites and are required for Rhodoquinone-dependent metabolism as targets for new anthelmintics

Soil transmitted helminths (STHs) are major human pathogens that infect over a billion people. Resistance to current anthelmintics is rising and new drugs are needed. Here we combine multiple approaches to find druggable targets in the anaerobic metabolic pathways STHs need to survive in their mammalian host. These require rhodoquinone (RQ), an electron carrier used by STHs and not their hosts. We identified 25 genes predicted to act in RQ-dependent metabolism including sensing hypoxia and RQ synthesis and found 9 are required. Since all 9 have mammalian orthologues, we used comparative genomics and structural modeling to identify those with active sites that differ between host and parasite. Together, we found 4 genes that are required for RQ-dependent metabolism and have different active sites. Finding these high confidence targets can open up in silico screens to identify species selective inhibitors of these enzymes as new anthelmintics.

Newly identified parasitic nematode beta-tubulin alleles confer resistance to benzimidazoles

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.

•

Four newly identified parasitic nematode beta-tubulin alleles confer benzimidazole resistance.

•

The four newly identified alleles do not confer deleterious fitness consequences.

•

Five beta-tubulin alleles confer recessive benzimidazole resistance.

Characterization of the β-tubulin gene family in Ascaris lumbricoides and Ascaris suum and its implication for the molecular detection of benzimidazole resistance

Background

The treatment coverage of control programs providing benzimidazole (BZ) drugs to eliminate the morbidity caused by soil-transmitted helminths (STHs) is unprecedently high. This high drug pressure may result in the development of BZ resistance in STHs and so there is an urgent need for surveillance systems detecting molecular markers associated with BZ resistance. A critical prerequisite to develop such systems is an understanding of the gene family encoding β-tubulin proteins, the principal targets of BZ drugs.

Methodology and principal findings

First, the β-tubulin gene families of Ascaris lumbricoides and Ascaris suum were characterized through the analysis of published genomes. Second, RNA-seq and RT-PCR analyses on cDNA were applied to determine the transcription profiles of the different gene family members. The results revealed that Ascaris species have at least seven different β-tubulin genes of which two are highly expressed during the entire lifecycle. Third, deep amplicon sequencing was performed on these two genes in more than 200 adult A. lumbricoides (Ethiopia and Tanzania) and A. suum (Belgium) worms, to investigate the intra- and inter-species genetic diversity and the presence of single nucleotide polymorphisms (SNPs) that are associated with BZ resistance in other helminth species; F167Y (TTC>TAC or TTT>TAT), E198A (GAA>GCA or GAG>GCG), E198L (GAA>TTA) and F200Y (TTC>TAC or TTT>TAT). These particular SNPs were absent in the two investigated genes in all three Ascaris populations.

Significance

This study demonstrated the presence of at least seven β-tubulin genes in Ascaris worms. A new nomenclature was proposed and prioritization of genes for future BZ resistance research was discussed. This is the first comprehensive description of the β-tubulin gene family in Ascaris and provides a framework to investigate the prevalence and potential role of β-tubulin sequence polymorphisms in BZ resistance in a more systematic manner than previously possible.

Benzimidazole (BZ) drugs remain the standard of treatment in large-scale deworming programs that aim to control the morbidity caused by intestinal worms. As these deworming programs are expanding world-wide, there is an increasing risk of worms becoming resistant to BZ drugs, highlighting the necessity for tools to detect gene mutations associated with drug resistance. However, the development of such tools is impeded by a lack of insights into the genes that are coding for β-tubulin proteins, which are the principal targets of BZ drugs. The aim of this study was to comprehensively characterize these genes in the worm species Ascaris lumbricoides and Ascaris suum. The findings highlight that these species have at least seven β-tubulin genes. Only two genes are highly expressed throughout the different life stages of the worm, and hence are more likely to be involved in the development of BZ resistance. No mutations that have previously been associated with BZ resistance in other intestinal worms were found. This study provides a baseline towards more efficient and accurate monitoring of drug resistance in large-scale deworming programs.

The genetics of gene expression in a Caenorhabditis elegans multiparental recombinant inbred line population

Studying genetic variation of gene expression provides a powerful way to unravel the molecular components underlying complex traits. Expression quantitative trait locus (eQTL) studies have been performed in several different model species, yet most of these linkage studies have been based on the genetic segregation of two parental alleles. Recently, we developed a multiparental segregating population of 200 recombinant inbred lines (mpRILs) derived from four wild isolates (JU1511, JU1926, JU1931, and JU1941) in the nematode Caenorhabditis elegans. We used RNA-seq to investigate how multiple alleles affect gene expression in these mpRILs. We found 1789 genes differentially expressed between the parental lines. Transgression, expression beyond any of the parental lines in the mpRILs, was found for 7896 genes. For expression QTL mapping almost 9000 SNPs were available. By combining these SNPs and the RNA-seq profiles of the mpRILs, we detected almost 6800 eQTLs. Most trans-eQTLs (63%) co-locate in six newly identified trans-bands. The trans-eQTLs found in previous two-parental allele eQTL experiments and this study showed some overlap (17.5-46.8%), highlighting on the one hand that a large group of genes is affected by polymorphic regulators across populations and conditions, on the other hand, it shows that the mpRIL population allows identification of novel gene expression regulatory loci. Taken together, the analysis of our mpRIL population provides a more refined insight into C. elegans complex trait genetics and eQTLs in general, as well as a starting point to further test and develop advanced statistical models for detection of multiallelic eQTLs and systems genetics studying the genotype-phenotype relationship.

Multiple drug resistance in hookworms infecting greyhound dogs in the USA

Ancylostoma caninum is the most prevalent nematode parasite of dogs. We confirmed multiple-drug resistance (MDR) in several A. caninum isolates to all anthelmintic drug classes approved for the treatment of hookworms in dogs in the USA. Cases of MDR hookworms appear to be highly overrepresented in greyhounds. The aims of this study were to evaluate the drug-resistant phenotypes and genotypes of the A. caninum infecting greyhounds. Fecal samples from greyhounds of the USA were acquired from two greyhound adoption kennels, one active greyhound racing kennel, and three veterinary practices. Fecal egg counts (FECs) were performed on fecal samples from 219 greyhounds, and despite treatment with anthelmintics, the mean FEC was 822.4 eggs per gram (EPG). Resistance to benzimidazoles and macrocyclic lactones were measured using the egg hatch assay (EHA) and the larval development assay (LDA), respectively. We performed 23 EHA and 22 LDA on either individual or pooled feces, representing 54 animals. Mean and median IC₅₀ and IC₉₅ values for the EHA were 5.3 μM, 3.6 μM, and 24.5 μM, 23.4 μM, respectively. For the LDA, the median IC₅₀ value was >1000 nM. These values ranged 62–81 times higher than our susceptible laboratory isolate. Only post-treatment samples were available. For samples collected <10 days post-treatment with albendazole, moxidectin, or a combination of febantel-pyrantel-moxidectin, the mean FEC were 349, 333, and 835 EPG, respectively. We obtained DNA from hookworm eggs isolated from 70 fecal samples, comprised of 60 individual dogs and 10 pools. Deep sequencing of the isotype 1 β-tubulin gene only revealed the presence of the F167Y (TTC>TAC) resistance polymorphism in 99% of these samples. These clinical, in vitro, and genetic data provide strong evidence that greyhound dogs in the USA are infected with MDR A. caninum at very high levels in prevalence and infection intensity.

•

Conclusive evidence that most racing greyhounds in the USA are infected with MDR A. caninum.

•

79% of the samples from racing or retired greyhounds were positive.

•

IC₅₀ values for BZs and MLs were 62–81 times higher in greyhounds.

•

The F167Y SNP was detected in 99% of samples, and in more than 2/3 the frequency was ≥ 75%.

Discovery of new vascular disrupting agents based on evolutionarily conserved drug action, pesticide resistance mutations, and humanized yeast

Thiabendazole (TBZ) is an FDA-approved benzimidazole widely used for its antifungal and antihelminthic properties. We showed previously that TBZ is also a potent vascular disrupting agent and inhibits angiogenesis at the tissue level by dissociating vascular endothelial cells in newly formed blood vessels. Here, we uncover TBZ's molecular target and mechanism of action. Using human cell culture, molecular modeling, and humanized yeast, we find that TBZ selectively targets only 1 of 9 human β-tubulin isotypes (TUBB8) to specifically disrupt endothelial cell microtubules. By leveraging epidemiological pesticide resistance data and mining chemical features of commercially used benzimidazoles, we discover that a broader class of benzimidazole compounds, in extensive use for 50 years, also potently disrupt immature blood vessels and inhibit angiogenesis. Thus, besides identifying the molecular mechanism of benzimidazole-mediated vascular disruption, this study presents evidence relevant to the widespread use of these compounds while offering potential new clinical applications.

Natural variation in fecundity is correlated with species-wide levels of divergence in Caenorhabditis elegans

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.

Natural variation in Caenorhabditis elegans responses to the anthelmintic emodepside

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.

From QTL to gene: C. elegans facilitates discoveries of the genetic mechanisms underlying natural variation

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.

The Caenorhabditis elegans and Haemonchus contortus beta-tubulin genes cannot substitute for loss of the Saccharomyces cerevisiae beta-tubulin gene

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.

Two novel loci underlie natural differences in Caenorhabditis elegans abamectin responses

Parasitic nematodes cause a massive worldwide burden on human health along with a loss of livestock and agriculture productivity. Anthelmintics have been widely successful in treating parasitic nematodes. However, resistance is increasing, and little is known about the molecular and genetic causes of resistance for most of these drugs. The free-living roundworm Caenorhabditis elegans provides a tractable model to identify genes that underlie resistance. Unlike parasitic nematodes, C. elegans is easy to maintain in the laboratory, has a complete and well annotated genome, and has many genetic tools. Using a combination of wild isolates and a panel of recombinant inbred lines constructed from crosses of two genetically and phenotypically divergent strains, we identified three genomic regions on chromosome V that underlie natural differences in response to the macrocyclic lactone (ML) abamectin. One locus was identified previously and encodes an alpha subunit of a glutamate-gated chloride channel (glc-1). Here, we validate and narrow two novel loci using near-isogenic lines. Additionally, we generate a list of prioritized candidate genes identified in C. elegans and in the parasite Haemonchus contortus by comparison of ML resistance loci. These genes could represent previously unidentified resistance genes shared across nematode species and should be evaluated in the future. Our work highlights the advantages of using C. elegans as a model to better understand ML resistance in parasitic nematodes.

Susceptible trichostrongyloid species mask presence of benzimidazole-resistant Haemonchus contortus in cattle

Background

Benzimidazole (BZ) anthelmintics are widely used to control infections with parasitic nematodes, but BZ resistance is an emerging threat among several nematode species infecting humans and animals. In Sudan, BZ-resistant Haemonchus contortus populations were recently reported in goats in South Darfur State. The objective of this study was to collect data regarding the situation of BZ resistance in cattle parasitic nematodes in South Darfur using phenotypic and molecular approaches, besides providing some epidemiological data on nematodes in cattle.

Methods

The faecal egg count reduction test and the egg hatch test (EHT) were used to evaluate benzimidazole efficacy in cattle nematodes in five South Darfur study areas: Beleil, Kass, Nyala, Rehed Al-Birdi and Tulus. Genomic DNA was extracted from pools of third-stage larvae (L3) (n = 40) during trials, before and after treatment, and pools of adult male Haemonchus spp. (n = 18) from abattoirs. The polymorphisms F167Y, E198A and F200Y in isotype 1 β-tubulin genes of H. contortus and H. placei were analysed using Sanger and pyrosequencing.

Results

Prevalence of gastro-intestinal helminths in cattle was 71% (313/443). Reduced albendazole faecal egg count reduction efficacy was detected in three study areas: Nyala (93.7%), Rehed Al-Birdi (89.7%) and Tulus (88.2%). In the EHT, EC₅₀ values of these study areas ranged between 0.032 and 0.037 µg/ml thiabendazole. Genus-specific PCRs detected the genera Haemonchus, Trichostrongylus and Cooperia in L3 samples collected after albendazole treatment. Sanger sequencing followed by pyrosequencing assays did not detect elevated frequencies of known BZ resistance-associated alleles in codon F167Y, E198A and F200Y in isotype 1 β-tubulin gene of H. placei (≤ 11.38%). However, polymorphisms were detected in H. contortus and in samples with mixed infections with H. contortus and H. placei at codon 198, including E198L (16/58), E198V (2/58) and potentially E198Stop (1/58). All pooled L3 samples post-albendazole treatment (n = 13) were identified as H. contortus with an E198L substitution at codon 198.

Conclusions

To the knowledge of the authors, this is the first report of reduced albendazole efficacy in cattle in Sudan and is the first study describing an E198L substitution in phenotypically BZ-resistant nematodes collected from cattle.

Considerations for anthelmintic resistance emergence in hookworm at a single locus

Over 800 million people are infected with hookworms around the world. Hookworms of the genus Ancylostoma and Necator are examples of nematodes that harbor the ability to enter a host by penetrating the skin, and after entry the infective larvae migrate to the small intestine where they encounter host-specific signals that initiate developmental pathways and culminate in maturation to the adult stage. Currently no vaccine is available for the treatment of hookworm infection. The control strategy is limited to anthelmintic drugs, which run the risk of losing efficacy as resistance grows. Genetic resistance has developed against all classes of anthelmintic drugs against livestock parasites, and recently markers of anthelmintic resistance in human hookworm populations have been reported. As anthelmintic resistance develops in human populations of hookworm, new drugs and novel control methods like vaccines will be required in the future to control hookworm transmission. This review outlines how population genetics and anthelmintic resistance could interact at a single locus to influence current control strategies.

•

This review outlines how population genetics and anthelmintic resistance may interact to influence current control strategies.

•

Altering mutation rate modifies rate of emergence of anthelmintic resistance.

•

Genotypic fitness changes the efficacy of anthelmintic administration.

Complementary Approaches with Free-living and Parasitic Nematodes to Understanding Anthelmintic Resistance

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.

Quantitative benzimidazole resistance and fitness effects of parasitic nematode beta-tubulin alleles

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.