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.

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.

Genetic Variants That Modify the Neuroendocrine Regulation of Foraging Behavior in C. elegans.

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. The neuroendocrine TGF-beta ligand, DAF-7, regulates diverse behavioral responses of Caenorhabditis elegans 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 and identified common variants in gap-2, encoding a GTPase-Activating Protein homologous to mammalian SynGAP proteins, which modify daf-7 expression cell-non-autonomously and promote exploratory foraging behavior in a 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.

Novel and improved Caenorhabditis briggsae gene models generated by community curation.

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.

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.

Novel and improved Caenorhabditis briggsae gene models generated by community curation.

Background: The nematode Caenorhabditis briggsae has been used as a model for genomics studies compared to Caenorhabditis elegans because of its striking morphological and behavioral similarities. These studies yielded numerous findings that have expanded our understanding of nematode development and evolution. However, the potential of C. briggsae to study nematode biology is limited by the quality of its genome resources. The reference genome and gene models 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 protein-coding gene predictions generated from short- and long-read transcriptomic data. Because of the limitations of gene prediction software, the existing gene models for QX1410 contain numerous errors in their structure and coding sequences. In this study, a team of researchers manually inspected over 21,000 software-derived gene models and underlying transcriptomic data to improve the protein-coding gene models of the C. briggsae QX1410 genome. Results: We designed a detailed workflow to train a team of nine students to manually curate genes using RNA read alignments and predicted gene models. We manually inspected the gene models using the genome annotation editor, Apollo, and proposed corrections to the coding sequences of over 8,000 genes. Additionally, we 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 before and after curation. Manual curation led to a substantial improvement in the protein sequence length accuracy of QX1410 genes. We also compared the curated QX1410 gene models against the existing AF16 gene models. The manual curation efforts yielded QX1410 gene models that are similar in quality to the extensively curated AF16 gene models in terms of protein-length accuracy and biological completeness scores. 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. Comparative genomic analysis using a related species with high-quality reference genome(s) and gene models can be used to quantify improvements in gene model quality in a newly sequenced genome. The detailed protocols provided in this work can be useful for future large-scale manual curation projects in other species. The chromosome-level reference genome for the C. briggsae strain QX1410 far surpasses the quality of the genome of the laboratory strain AF16, and our manual curation efforts have brought the QX1410 gene models to a comparable level of quality to the previous reference, AF16. The improved genome resources for C. briggsae provide reliable tools for the study of Caenorhabditis biology and other related nematodes.

Natural genetic variation in the pheromone production of C. elegans.

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.

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.

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.

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.

The turkey ascarid, Ascaridia dissimilis, as a model genetic system.

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.

Different structural variant prediction tools yield considerably different results in Caenorhabditis elegans.

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.

Characterization of larval growth in C. elegans cuticle mutants.

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.

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.

An anchored experimental design and meta-analysis approach to address batch effects in large-scale metabolomics.

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.

Natural variation in C. elegans short tandem repeats.

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.

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.