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 (EC 10 ) 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 EC 10 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.
Abstract Recent efforts to understand the natural niche of the keystone model organism Caenorhabditis elegans have suggested that this species is cosmopolitan and associated with rotting vegetation and fruits. However, most of the strains isolated from nature have low genetic diversity likely because recent chromosome-scale selective sweeps contain alleles that increase fitness in human-associated habitats. Strains from the Hawaii Islands are highly divergent from non-Hawaiian strains. This result suggests that Hawaiian strains might contain ancestral genetic diversity that was purged from most non-Hawaiian strains by the selective sweeps. To characterize the genetic diversity and niche of Hawaiian C. elegans , we sampled across the Hawaiian Islands and isolated 100 new C. elegans strains. We found that C. elegans strains are not associated with any one substrate but are found in cooler climates at high elevations. These Hawaiian strains are highly diverged compared to the rest of the global population. Admixture analysis identified 11 global populations, four of which are from Hawaii. Surprisingly, one of the Hawaiian populations shares recent ancestry with non-Hawaiian populations, including portions of globally swept haplotypes. This discovery provides the first evidence of gene flow between Hawaiian and non-Hawaiian populations. Most importantly, the high levels of diversity observed in Hawaiian strains might represent the complex patterns of ancestral genetic diversity in the C. elegans species before human influence.
Hawaiian isolates of the nematode species Caenorhabditis elegans have long been known to harbor genetic diversity greater than the rest of the worldwide population, but this observation was supported by only a small number of wild strains. To better characterize the niche and genetic diversity of Hawaiian C. elegans and other Caenorhabditis species, we sampled different substrates and niches across the Hawaiian islands. We identified hundreds of new Caenorhabditis strains from known species and a new species, Caenorhabditis oiwi. Hawaiian C. elegans are found in cooler climates at high elevations but are not associated with any specific substrate, as compared to other Caenorhabditis species. Surprisingly, admixture analysis revealed evidence of shared ancestry between some Hawaiian and non-Hawaiian C. elegans strains. We suggest that the deep diversity we observed in Hawaii might represent patterns of ancestral genetic diversity in the C. elegans species before human influence.
Abstract Large-scale ecological sampling can be difficult and costly, especially for organisms that are too small to be easily identified in a natural environment by eye. Typically, these microscopic floral and fauna are sampled by collecting substrates from nature and then separating organisms from substrates in the laboratory. In many cases, diverse organisms can be identified to the species-level using molecular barcodes. To facilitate large-scale ecological sampling of microscopic organisms, we used a geographic data-collection platform for mobile devices called Fulcrum that streamlines the organization of geospatial sampling data, substrate photographs, and environmental data at natural sampling sites. These sampling data are then linked to organism isolation data from the laboratory. Here, we describe the easyFulcrum R package, which can be used to clean, process, and visualize ecological field sampling and isolation data exported from the Fulcrum mobile application. We developed this package for wild nematode sampling, but it is extensible to other organisms. The advantages of using Fulcrum combined with easyFulcrum are (1) the elimination of transcription errors by replacing manual data entry and/or spreadsheets with a mobile application, (2) the ability to clean, process, and visualize sampling data using a standardized set of functions in the R software environment, and (3) the ability to join disparate data to each other, including environmental data from the field and the molecularly defined identities of individual specimens isolated from samples.
Abstract Important clues about natural selection can be gleaned from discrepancies between the properties of segregating genetic variants and of mutations accumulated experimentally under minimal selection, provided the mutational process is the same in the lab as in nature. The ratio of transitions to transversions (Ts/Tv) is consistently lower in C. elegans mutation accumulation (MA) experiments than in nature, which has been argued to be in part due to increased oxidative stress in the lab environment. Using whole-genome sequence data from a set of C. elegans MA lines carrying a mutation ( mev-1 ) that increases the cellular titer of reactive oxygen species (ROS), leading to increased endogenous oxidative stress, we find that the base-substitution spectrum is similar between mev-1 lines, its wild-type progenitor (N2), and another set of MA lines derived from a different wild strain (PB306). By contrast, the rate of short insertions is greater in the mev-1 lines, consistent with studies in other organisms in which environmental stress led to an increase in the rate of insertion-deletion mutations. Further, the mutational properties of mononucleotide repeats in all strains are qualitatively different from those of non-mononucleotide sequence, both for indels and base-substitutions, and whereas the non-mononucleotide spectra are fairly similar between MA lines and wild isolates, the mononucleotide spectra are very different. The discrepancy in mutational spectra between lab MA experiments and natural variation is likely due to a consistent (but unknown) effect of the lab environment that manifests itself via different modes of mutability and/or repair at mononucleotide loci.
Abstract Organismal fitness is relevant in many contexts in biology. The most meaningful experimental measure of fitness is competitive fitness, when two or more entities (e.g., genotypes) are allowed to compete directly. In theory, competitive fitness is simple to measure: an experimental population is initiated with the different types in known proportions and allowed to evolve under experimental conditions to a predefined endpoint. In practice, there are several obstacles to obtaining robust estimates of competitive fitness in multicellular organisms, the most pervasive of which is simply the time it takes to count many individuals of different types from many replicate populations. Methods by which counting can be automated in high throughput are desirable, but for automated methods to be useful, the bias and technical variance associated with the method must be (a) known, and (b) sufficiently small relative to other sources of bias and variance to make the effort worthwhile. The nematode Caenorhabditis elegans is an important model organism, and the fitness effects of genotype and environmental conditions are often of interest. We report a comparison of three experimental methods of quantifying competitive fitness, in which wild-type strains are competed against GFP-marked competitors under standard laboratory conditions. Population samples were split into three replicates and counted (1) “by eye” from a saved image, (2) from the same image using CellProfiler image analysis software, and (3) with a large particle flow cytometer (a “worm sorter”). From 720 replicate samples, neither the frequency of wild-type worms nor the among-sample variance differed significantly between the three methods. CellProfiler and the worm sorter provide at least a tenfold increase in sample handling speed with little (if any) bias or increase in variance.
Abstract Across diverse taxa, selfing species have evolved independently from outcrossing species thousands of times. The transition from outcrossing to selfing significantly decreases the effective population size, effective recombination rate, and heterozygosity within a species. These changes lead to a reduction in genetic diversity, and therefore adaptive potential, by intensifying the effects of random genetic drift and linked selection. Within the nematode genus Caenorhabditis , selfing has evolved at least three times and all three species, including in the model organism Caenorhabditis elegans , show substantially reduced genetic diversity relative to outcrossing species. Selfing and outcrossing Caenorhabditis species are often found in the same niches, but we still do not know how selfing species with limited genetic diversity can adapt to these environments. Here, we examine the whole-genome sequences from 609 wild C. elegans strains isolated worldwide and show that genetic variation is concentrated in punctuated hyper-divergent regions that cover 20% of the C. elegans reference genome. These regions are enriched in environmental response genes that mediate sensory perception, pathogen response, and xenobiotic stress response. Population genomic evidence suggests that genetic diversity in these regions has been maintained by long-term balancing selection. Using long-read genome assemblies for 15 wild strains, we show that hyper-divergent haplotypes contain unique sets of genes and show levels of divergence comparable to levels found between Caenorhabditis species that diverged millions of years ago. These results provide an example for how species can avoid the evolutionary “dead end” associated with selfing.
