publications

2023

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

   Amanda O Shaver ,  Janneke Wit ,  Clayton M Dilks , and 4 more authors

   PLoS Pathog., Apr 2023

   Abs

   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.

2. CaeNDR, the Caenorhabditis Natural Diversity Resource

   Timothy A Crombie ,  Ryan McKeown ,  Nicolas D Moya , and 12 more authors

   Nucleic Acids Res., Oct 2023

   Abs

   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.

2022

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

   Timothy A Crombie ,  Paul Battlay ,  Robyn E Tanny , and 12 more authors

   Mol. Ecol., Apr 2022

   Abs

   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.

2. An automated approach to quantify chemotaxis index in C. elegans

   Timothy A Crombie ,  Chido Chikuturudzi ,  Daniel E Cook , and 1 more author

   MicroPubl Biol, May 2022

   Abs

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

3. A Highly Scalable Approach to Perform Ecological Surveys of Selfing Caenorhabditis Nematodes

   Timothy A Crombie ,  Robyn E Tanny ,  Claire M Buchanan , and 2 more authors

   J. Vis. Exp., Mar 2022

   Abs

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

4. C. elegans toxicant responses vary among genetically diverse individuals

   Samuel J Widmayer ,  Timothy A Crombie ,  Joy N Nyaanga , and 2 more authors

   Toxicology, Aug 2022

   Abs

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

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

   Bennett W Fox ,  Olga Ponomarova ,  Yong-Uk Lee , and 14 more authors

   Nature, Jul 2022

   Abs

   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.

6. A Highly Scalable Approach to Perform Ecological Surveys of Selfing Caenorhabditis Nematodes

   Timothy A Crombie ,  Robyn E Tanny ,  Claire M Buchanan , and 2 more authors

   J. Vis. Exp., Mar 2022

   Abs

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

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

   Timothy A Crombie ,  Paul Battlay ,  Robyn E Tanny , and 12 more authors

   Mol. Ecol., Apr 2022

   Abs

   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.

2021

1. Mutability of mononucleotide repeats, not oxidative stress, explains the discrepancy between laboratory-accumulated mutations and the natural allele-frequency spectrum in C. elegans

   Moein Rajaei ,  Ayush Shekhar Saxena ,  Lindsay M Johnson , and 6 more authors

   Genome Res., Sep 2021

   Abs

   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 laboratory as in nature. The base-substitution spectrum differs between C. elegans laboratory mutation accumulation (MA) experiments and the standing site-frequency spectrum, which has been argued to be in part owing to increased oxidative stress in the laboratory environment. Using genome sequence data from C. elegans MA lines carrying a mutation (mev-1) that increases the cellular titer of reactive oxygen species (ROS), leading to increased oxidative stress, we find the base-substitution spectrum is similar between mev-1, its wild-type progenitor (N2), and another set of MA lines derived from a different wild strain (PB306). Conversely, the rate of short insertions is greater in mev-1, consistent with studies in other organisms in which environmental stress increased the rate of insertion-deletion mutations. Further, the mutational properties of mononucleotide repeats in all strains are different from those of nonmononucleotide sequence, both for indels and base-substitutions, and whereas the nonmononucleotide spectra are fairly similar between MA lines and wild isolates, the mononucleotide spectra are very different, with a greater frequency of A:T →T:A transversions and an increased proportion of \pm1-bp indels. The discrepancy in mutational spectra between laboratory MA experiments and natural variation is likely owing to a consistent (but unknown) effect of the laboratory environment that manifests itself via different modes of mutability and/or repair at mononucleotide loci.

2. easyFulcrum: An R package to process and analyze ecological sampling data generated using the Fulcrum mobile application

   Matteo Di Bernardo ,  Timothy A Crombie ,  Daniel E Cook , and 1 more author

   PLoS One, Oct 2021

   Abs

   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 can be used with 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.

3. High-throughput imaging of Caenorhabditis elegans aging using collective activity monitoring

   Anthony D Fouad ,  Matthew A Churgin ,  Julia Hayden , and 15 more authors

   Oct 2021

   Abs

   The genetic manipulability and short lifespan of C. elegans make it an important model for aging research. Widely applied methods for measurements of worm aging based on manual observation are labor intensive and low-throughput. Here, we describe the Worm Collective Activity Monitoring Platform (WormCamp), a system for assaying aging in C. elegans by monitoring activity of populations of worms in standard 24-well plates. We show that metrics based on the rate of decline in collective activity can be used to estimate the average lifespan and locomotor healthspan in the population. Using the WormCamp, we assay a panel of highly divergent natural isolates of C. elegans and show that both lifespan and locomotor healthspan display substantial heritability. To facilitate analysis of large numbers of worms, we developed a robotic imaging system capable of simultaneous automated monitoring of activity, lifespan, and locomotor healthspan in up to 2,304 populations containing a total of  90,000 animals. We applied the automated system to conduct a large-scale RNA interference screen for genes that affect lifespan and locomotor healthspan. The WormCamp system is complementary to other current automated methods for assessing C. elegans aging and is well suited for efficiently screening large numbers of conditions. ### Competing Interest Statement The authors have declared no competing interest.

4. easyXpress: An R package to analyze and visualize high-throughput C. elegans microscopy data generated using CellProfiler

   Joy Nyaanga ,  Timothy A Crombie ,  Samuel J Widmayer , and 1 more author

   PLoS One, Aug 2021

   Abs

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

5. Balancing selection maintains hyper-divergent haplotypes in Caenorhabditis elegans

   Daehan Lee ,  Stefan Zdraljevic ,  Lewis Stevens , and 12 more authors

   Nat Ecol Evol, Apr 2021

   Abs

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

2019

1. Natural variation in C. elegans arsenic toxicity is explained by differences in branched chain amino acid metabolism

   Stefan Zdraljevic ,  Bennett William Fox ,  Christine Strand , and 7 more authors

   Elife, Apr 2019

   Abs

   We find that variation in the dbt-1 gene underlies natural differences in Caenorhabditis elegans responses to the toxin arsenic. This gene encodes the E2 subunit of the branched-chain α-keto acid dehydrogenase (BCKDH) complex, a core component of branched-chain amino acid (BCAA) metabolism. We causally linked a non-synonymous variant in the conserved lipoyl domain of DBT-1 to differential arsenic responses. Using targeted metabolomics and chemical supplementation, we demonstrate that differences in responses to arsenic are caused by variation in iso-branched chain fatty acids. Additionally, we show that levels of branched chain fatty acids in human cells are perturbed by arsenic treatment. This finding has broad implications for arsenic toxicity and for arsenic-focused chemotherapeutics across human populations. Our study implicates the BCKDH complex and BCAA metabolism in arsenic responses, demonstrating the power of C. elegans natural genetic diversity to identify novel mechanisms by which environmental toxins affect organismal physiology. Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor’s assessment is that all the issues have been addressed (see decision letter).

2. Deep sampling of Hawaiian Caenorhabditis elegans reveals high genetic diversity and admixture with global populations

   Timothy A Crombie ,  Stefan Zdraljevic ,  Daniel E Cook , and 11 more authors

   Elife, Dec 2019

   Abs

   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.

3. Natural variation in C. elegans arsenic toxicity is explained by differences in branched chain amino acid metabolism

   Stefan Zdraljevic ,  Bennett William Fox ,  Christine Strand , and 7 more authors

   Elife, Apr 2019

   Abs

   We find that variation in the dbt-1 gene underlies natural differences in Caenorhabditis elegans responses to the toxin arsenic. This gene encodes the E2 subunit of the branched-chain α-keto acid dehydrogenase (BCKDH) complex, a core component of branched-chain amino acid (BCAA) metabolism. We causally linked a non-synonymous variant in the conserved lipoyl domain of DBT-1 to differential arsenic responses. Using targeted metabolomics and chemical supplementation, we demonstrate that differences in responses to arsenic are caused by variation in iso-branched chain fatty acids. Additionally, we show that levels of branched chain fatty acids in human cells are perturbed by arsenic treatment. This finding has broad implications for arsenic toxicity and for arsenic-focused chemotherapeutics across human populations. Our study implicates the BCKDH complex and BCAA metabolism in arsenic responses, demonstrating the power of C. elegans natural genetic diversity to identify novel mechanisms by which environmental toxins affect organismal physiology. Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor’s assessment is that all the issues have been addressed (see decision letter).

2018

1. Head-to-head comparison of three experimental methods of quantifying competitive fitness in C. elegans

   Timothy A Crombie ,  Sayran Saber ,  Ayush Shekhar Saxena , and 2 more authors

   PLoS One, Oct 2018

   Abs

   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.

2. The mutational decay of male-male and hermaphrodite-hermaphrodite competitive fitness in the androdioecious nematode C. elegans

   Shu-Dan Yeh ,  Ayush Shekhar Saxena ,  Timothy A Crombie , and 6 more authors

   Heredity, Jan 2018

   Abs

   Androdioecious Caenorhabditis have a high frequency of self-compatible hermaphrodites and a low frequency of males. The effects of mutations on male fitness are of interest for two reasons. First, when males are rare, selection on male-specific mutations is less efficient than in hermaphrodites. Second, males may present a larger mutational target than hermaphrodites because of the different ways in which fitness accrues in the two sexes. We report the first estimates of male-specific mutational effects in an androdioecious organism. The rate of male-specific inviable or sterile mutations is ⩽5 \times 10-4/generation, below the rate at which males would be lost solely due to those kinds of mutations. The rate of mutational decay of male competitive fitness is   0.17%/generation; that of hermaphrodite competitive fitness is   0.11%/generation. The point estimate of   1.5X faster rate of mutational decay of male fitness is nearly identical to the same ratio in Drosophila. Estimates of mutational variance (VM) for male mating success and competitive fitness are not significantly different from zero, whereas VM for hermaphrodite competitive fitness is similar to that of non-competitive fitness. Two independent estimates of the average selection coefficient against mutations affecting hermaphrodite competitive fitness agree to within two-fold, 0.33-0.5%.

3. Head-to-head comparison of three experimental methods of quantifying competitive fitness in C. elegans

   Timothy A Crombie ,  Sayran Saber ,  Ayush Shekhar Saxena , and 2 more authors

   PLoS One, Oct 2018

   Abs

   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.

4. Head-to-Head Comparison of Three Methods of Quantifying Competitive Fitness in C. elegans

   Timothy A Crombie ,  Sayran Saber ,  Ayush Shekhar Saxena , and 2 more authors

   Jul 2018

   Abs

   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.

5. The mutational decay of male-male and hermaphrodite-hermaphrodite competitive fitness in the androdioecious nematode C. elegans

   Shu-Dan Yeh ,  Ayush Shekhar Saxena ,  Timothy A Crombie , and 6 more authors

   Heredity, Jan 2018

   Abs

   Androdioecious Caenorhabditis have a high frequency of self-compatible hermaphrodites and a low frequency of males. The effects of mutations on male fitness are of interest for two reasons. First, when males are rare, selection on male-specific mutations is less efficient than in hermaphrodites. Second, males may present a larger mutational target than hermaphrodites because of the different ways in which fitness accrues in the two sexes. We report the first estimates of male-specific mutational effects in an androdioecious organism. The rate of male-specific inviable or sterile mutations is ⩽5 \times 10-4/generation, below the rate at which males would be lost solely due to those kinds of mutations. The rate of mutational decay of male competitive fitness is   0.17%/generation; that of hermaphrodite competitive fitness is   0.11%/generation. The point estimate of   1.5X faster rate of mutational decay of male fitness is nearly identical to the same ratio in Drosophila. Estimates of mutational variance (VM) for male mating success and competitive fitness are not significantly different from zero, whereas VMfor hermaphrodite competitive fitness is similar to that of non-competitive fitness. Two independent estimates of the average selection coefficient against mutations affecting hermaphrodite competitive fitness agree to within two-fold, 0.33-0.5%.

2017

1. The effect of temperature adaptation on the ubiquitin–proteasome pathway in notothenioid fishes

   A E Todgham ,  T A Crombie ,  and  others

   Journal of Experimental, Jan 2017

   Abs

   There is an accumulating body of evidence suggesting that the sub-zero Antarctic marine environment places physiological constraints on protein homeostasis. Levels of ubiquitin (Ub)-conjugated proteins, 20S proteasome activity and mRNA expression of many proteins …

2. The effect of temperature adaptation on the ubiquitin–proteasome pathway in notothenioid fishes

   Anne E Todgham ,  Timothy A Crombie ,  and  Gretchen E Hofmann

   J. Exp. Biol., Feb 2017

   Abs

   There is an accumulating body of evidence suggesting that the sub-zero Antarctic marine environment places physiological constraints on protein homeostasis. Levels of ubiquitin (Ub)-conjugated proteins, 20S proteasome activity and mRNA expression of many proteins involved in both the Ub tagging of damaged proteins as well as the different complexes of the 26S proteasome were measured to examine whether there is thermal compensation of the Ub–proteasome pathway in Antarctic fishes to better understand the efficiency of the protein degradation machinery in polar species. Both Antarctic ( Trematomus bernacchii , Pagothenia borchgrevinki ) and non-Antarctic ( Notothenia angustata , Bovichtus variegatus ) notothenioids were included in this study to investigate the mechanisms of cold adaptation of this pathway in polar species. Overall, there were significant differences in the levels of Ub-conjugated proteins between the Antarctic notothenioids and B. variegatus , with N. angustata possessing levels very similar to those of the Antarctic fishes. Proteasome activity in the gills of Antarctic fishes demonstrated a high degree of temperature compensation such that activity levels were similar to activities measured in their temperate relatives at ecologically relevant temperatures. A similar level of thermal compensation of proteasome activity was not present in the liver of two Antarctic fishes. Higher gill proteasome activity is likely due in part to higher cellular levels of proteins involved in the Ub–proteasome pathway, as evidenced by high mRNA expression of relevant genes. Reduced activity of the Ub–proteasome pathway does not appear to be the mechanism responsible for elevated levels of denatured proteins in Antarctic fishes, at least in the gills.

2016

1. Inhibition of the oxidative stress response by heat stress in Caenorhabditis elegans

   Timothy A Crombie ,  Lanlan Tang ,  Keith P Choe , and 1 more author

   J. Exp. Biol., Jul 2016

   Abs

   It has long been recognized that simultaneous exposure to heat stress and oxidative stress shows a synergistic interaction that reduces organismal fitness, but relatively little is known about the mechanisms underlying this interaction. We investigated the role of molecular stress responses in driving this synergistic interaction using the nematode Caenorhabditis elegans To induce oxidative stress, we used the pro-oxidant compounds acrylamide, paraquat and juglone. As expected, we found that heat stress and oxidative stress interact synergistically to reduce survival. Compared with exposure to each stressor alone, during simultaneous sublethal exposure to heat stress and oxidative stress the normal induction of key oxidative-stress response (OxSR) genes was generally inhibited, whereas the induction of key heat-shock response (HSR) genes was not. Genetically activating the SKN-1-dependent OxSR increased a marker for protein aggregation and decreased whole-worm survival during heat stress alone, with the latter being independent of HSF-1. In contrast, compared with wild-type worms, inactivating the HSR by HSF-1 knockdown, which would be expected to decrease basal heat shock protein expression, increased survival during oxidative stress alone. Taken together, these data suggest that, in C. elegans, the HSR and OxSR cannot be simultaneously activated to the same extent that each can be activated during a single stressor exposure. We conclude that the observed synergistic reduction in survival during combined exposure to heat stress and oxidative stress is due, at least in part, to inhibition of the OxSR during activation of the HSR.

2. Inhibition of the oxidative stress response by heat stress in Caenorhabditis elegans

   Timothy A Crombie ,  Lanlan Tang ,  Keith P Choe , and 1 more author

   J. Exp. Biol., Jul 2016

   Abs

   It has long been recognized that simultaneous exposure to heat stress and oxidative stress shows a synergistic interaction that reduces organismal fitness, but relatively little is known about the mechanisms underlying this interaction. We investigated the role of molecular stress responses in driving this synergistic interaction using the nematode Caenorhabditis elegans To induce oxidative stress, we used the pro-oxidant compounds acrylamide, paraquat and juglone. As expected, we found that heat stress and oxidative stress interact synergistically to reduce survival. Compared with exposure to each stressor alone, during simultaneous sublethal exposure to heat stress and oxidative stress the normal induction of key oxidative-stress response (OxSR) genes was generally inhibited, whereas the induction of key heat-shock response (HSR) genes was not. Genetically activating the SKN-1-dependent OxSR increased a marker for protein aggregation and decreased whole-worm survival during heat stress alone, with the latter being independent of HSF-1. In contrast, compared with wild-type worms, inactivating the HSR by HSF-1 knockdown, which would be expected to decrease basal heat shock protein expression, increased survival during oxidative stress alone. Taken together, these data suggest that, in C. elegans, the HSR and OxSR cannot be simultaneously activated to the same extent that each can be activated during a single stressor exposure. We conclude that the observed synergistic reduction in survival during combined exposure to heat stress and oxidative stress is due, at least in part, to inhibition of the OxSR during activation of the HSR.

2008

1. Aggression frequency and intensity, independent of testosterone levels, relate to neural activation within the dorsolateral subdivision of the ventromedial …

   D Kabelik ,  T Crombie ,  and  M C Moore

   Horm. Behav., Jul 2008

   Abs

   The mechanisms by which testosterone regulates aggression are unclear and may involve changes that alter the activity levels of one or more brain nuclei. We estimate neural activity by counting immunopositive cells against phosphorylated cyclic AMP response element …

2. Aggression frequency and intensity, independent of testosterone levels, relate to neural activation within the dorsolateral subdivision of the ventromedial hypothalamus in the tree lizard Urosaurus ornatus

   David Kabelik ,  Tim Crombie ,  and  Michael C Moore

   Horm. Behav., Jun 2008

   Abs

   The mechanisms by which testosterone regulates aggression are unclear and may involve changes that alter the activity levels of one or more brain nuclei. We estimate neural activity by counting immunopositive cells against phosphorylated cyclic AMP response element binding protein (pCREB). We demonstrate increased pCREB immunoreactivity within the dorsolateral subdivision of the ventromedial hypothalamus (VMHdl) following an aggressive encounter in male tree lizards Urosaurus ornatus. This immunoreactivity is induced both by exposure to and performance of aggressive behaviors. This dual activation of the VMHdl suggests its possible role as an integration center for assessment and expression of aggressive behavior. Furthermore, pCREB induction was greater in encounters involving higher frequency and intensity of aggressive display, demonstrating a direct relationship between neural activation and behavior. The VMHdl is also rich in steroid receptors. In a second experiment involving hormone manipulations, testosterone treatment increased aggression levels, though it did not increase the number of pCREB positive cells within the VMHdl. This lack of an effect of testosterone on pCREB induction within the VMHdl may be due to induction arising from the behaviors of conspecifics (especially in low-testosterone, low-aggression individuals), variation in aggression mediated by other variables, or regulation of aggression by circuits outside of the VMHdl. Together, these findings support a notion of the VMHdl as a nucleus involved in integrating afferent and efferent information within the neural aggression-control circuit.

2006

1. Systematic mapping of genetic interactions in Caenorhabditis elegans identifies common modifiers of diverse signaling pathways

   Ben Lehner ,  Catriona Crombie ,  Julia Tischler , and 2 more authors

   Nat. Genet., Aug 2006

   Abs

   Most heritable traits, including disease susceptibility, are affected by interactions between multiple genes. However, we understand little about how genes interact because very few possible genetic interactions have been explored experimentally. We have used RNA interference in Caenorhabditis elegans to systematically test approximately 65,000 pairs of genes for their ability to interact genetically. We identify approximately 350 genetic interactions between genes functioning in signaling pathways that are mutated in human diseases, including components of the EGF/Ras, Notch and Wnt pathways. Most notably, we identify a class of highly connected ’hub’ genes: inactivation of these genes can enhance the phenotypic consequences of mutation of many different genes. These hub genes all encode chromatin regulators, and their activity as genetic hubs seems to be conserved across animals. We propose that these genes function as general buffers of genetic variation and that these hub genes may act as modifier genes in multiple, mechanistically unrelated genetic diseases in humans.