Natural variation in a single amino acid substitution underlies physiological responses to topoisomerase II poisons

Biosafety evaluation of etoposide lipid nanomedicines in C. elegans

Neuroblastoma is a pediatric tumor that originates during embryonic development and progresses into aggressive tumors, primarily affecting children under two years old. Many patients are diagnosed as high-risk and undergo chemotherapy, often leading to short- and long-term toxicities. Nanomedicine offers a promising solution to enhance drug efficacy and improve physical properties. In this study, lipid-based nanomedicines were developed with an average size of 140 nm, achieving a high encapsulation efficiency of over 90% for the anticancer drug etoposide. Then, cytotoxicity and apoptosis-inducing effects of these etoposide nanomedicines were assessed in vitro using human cell lines, both cancerous and non-cancerous. The results demonstrated that etoposide nanomedicines exhibited high toxicity and selectively induced apoptosis only in cancerous cells.Next, the biosafety of these nanomedicines in C. elegans, a model organism, was evaluated by measuring survival, body size, and the effect on dividing cells. The findings showed that the nanomedicines had a safer profile than the free etoposide in this model. Notably, nanomedicines exerted etoposide's antiproliferative effect only in highly proliferative germline cells. Therefore, the developed nanomedicines hold promise as safe drug delivery systems for etoposide, potentially leading to an improved therapeutic index for neuroblastoma treatment.

Natural variation in protein kinase D modifies alcohol sensitivity in Caenorhabditis elegans

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

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.

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

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

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

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

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

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

A Single Nucleotide Variant in the PPARγ-homolog Eip75B Affects Fecundity in Drosophila

Single nucleotide polymorphisms are the most common type of genetic variation, but how these variants contribute to the adaptation of complex phenotypes is largely unknown. Experimental evolution and genome-wide association studies have demonstrated that variation in the PPARγ-homolog Eip75B has associated with longevity and life-history differences in the fruit fly Drosophila melanogaster. Using RNAi knockdown, we first demonstrate that reduced expression of Eip75B in adult flies affects lifespan, egg-laying rate, and egg volume. We then tested the effects of a naturally occurring SNP within a cis-regulatory domain of Eip75B by applying two complementary approaches: a Mendelian randomization approach using lines of the Drosophila Genetic Reference Panel, and allelic replacement using precise CRISPR/Cas9-induced genome editing. Our experiments reveal that this natural polymorphism has a significant pleiotropic effect on fecundity and egg-to-adult viability, but not on longevity or other life-history traits. Our results provide a rare functional validation at the nucleotide level and identify a natural allelic variant affecting fitness and life-history adaptation.

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.

Linkage mapping reveals loci that underlie differences in Caenorhabditis elegans growth

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

High-throughput small molecule screen identifies inhibitors of microsporidia invasion and proliferation in C. elegans

Microsporidia are a diverse group of fungal-related obligate intracellular parasites that infect most animal phyla. Despite the emerging threat that microsporidia represent to humans and agricultural animals, few reliable treatment options exist. Here, we develop a high-throughput screening method for the identification of chemical inhibitors of microsporidia infection, using liquid cultures of Caenorhabditis elegans infected with the microsporidia species Nematocida parisii. We screen a collection of 2560 FDA-approved compounds and natural products, and identify 11 candidate microsporidia inhibitors. Five compounds prevent microsporidia infection by inhibiting spore firing, whereas one compound, dexrazoxane, slows infection progression. The compounds have in vitro activity against several other microsporidia species, including those known to infect humans. Together, our results highlight the effectiveness of C. elegans as a model host for drug discovery against intracellular pathogens, and provide a scalable high-throughput system for the identification and characterization of microsporidia inhibitors.

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

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

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

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

Natural genetic variation as a tool for discovery in Caenorhabditis nematodes

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

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

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

•

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

•

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

•

Five beta-tubulin alleles confer recessive benzimidazole resistance.

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

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

Natural variation in Caenorhabditis elegans responses to the anthelmintic emodepside

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

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

Although many studies have examined quantitative trait variation across many species, only a small number of genes and thereby molecular mechanisms have been discovered. Without these data, we can only speculate about evolutionary processes that underlie trait variation. Here, we review how quantitative and molecular genetics in the nematode Caenorhabditis elegans led to the discovery and validation of 37 quantitative trait genes over the past 15 years. Using these data, we can start to make inferences about evolution from these quantitative trait genes, including the roles that coding versus noncoding variation, gene family expansion, common versus rare variants, pleiotropy, and epistasis play in trait variation across this species.

Two novel loci underlie natural differences in Caenorhabditis elegans abamectin responses

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

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

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

Natural variation in the sequestosome-related gene, sqst-5, underlies zinc homeostasis in Caenorhabditis elegans

Zinc is an essential trace element that acts as a co-factor for many enzymes and transcription factors required for cellular growth and development. Altering intracellular zinc levels can produce dramatic effects ranging from cell proliferation to cell death. To avoid such fates, cells have evolved mechanisms to handle both an excess and a deficiency of zinc. Zinc homeostasis is largely maintained via zinc transporters, permeable channels, and other zinc-binding proteins. Variation in these proteins might affect their ability to interact with zinc, leading to either increased sensitivity or resistance to natural zinc fluctuations in the environment. We can leverage the power of the roundworm nematode Caenorhabditis elegans as a tractable metazoan model for quantitative genetics to identify genes that could underlie variation in responses to zinc. We found that the laboratory-adapted strain (N2) is resistant and a natural isolate from Hawaii (CB4856) is sensitive to micromolar amounts of exogenous zinc supplementation. Using a panel of recombinant inbred lines, we identified two large-effect quantitative trait loci (QTL) on the left arm of chromosome III and the center of chromosome V that are associated with zinc responses. We validated and refined both QTL using near-isogenic lines (NILs) and identified a naturally occurring deletion in sqst-5, a sequestosome-related gene, that is associated with resistance to high exogenous zinc. We found that this deletion is relatively common across strains within the species and that variation in sqst-5 is associated with zinc resistance. Our results offer a possible mechanism for how organisms can respond to naturally high levels of zinc in the environment and how zinc homeostasis varies among individuals.

Dissecting the eQTL Micro-Architecture in Caenorhabditis elegans

The study of expression quantitative trait loci (eQTL) using natural variation in inbred populations has yielded detailed information about the transcriptional regulation of complex traits. Studies on eQTL using recombinant inbred lines (RILs) led to insights on cis and trans regulatory loci of transcript abundance. However, determining the underlying causal polymorphic genes or variants is difficult, but ultimately essential for the understanding of regulatory networks of complex traits. This requires insight into whether associated loci are single eQTL or a combination of closely linked eQTL, and how this QTL micro-architecture depends on the environment. We addressed these questions by testing for independent replication of previously mapped eQTL in Caenorhabditis elegans using new data from introgression lines (ILs). Both populations indicate that the overall heritability of gene expression, number, and position of eQTL differed among environments. Across environments we were able to replicate 70% of the cis- and 40% of the trans-eQTL using the ILs. Testing eight different simulation models, we suggest that additive effects explain up to 60-93% of RIL/IL heritability for all three environments. Closely linked eQTL explained up to 40% of RIL/IL heritability in the control environment whereas only 7% in the heat-stress and recovery environments. In conclusion, we show that reproducibility of eQTL was higher for cis vs. trans eQTL and that the environment affects the eQTL micro-architecture.

Caenorhabditis elegans in anthelmintic research - Old model, new perspectives

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

Natural variation in a glucuronosyltransferase modulates propionate sensitivity in a C. elegans propionic acidemia model

Mutations in human metabolic genes can lead to rare diseases known as inborn errors of human metabolism. For instance, patients with loss-of-function mutations in either subunit of propionyl-CoA carboxylase suffer from propionic acidemia because they cannot catabolize propionate, leading to its harmful accumulation. Both the penetrance and expressivity of metabolic disorders can be modulated by genetic background. However, modifiers of these diseases are difficult to identify because of the lack of statistical power for rare diseases in human genetics. Here, we use a model of propionic acidemia in the nematode Caenorhabditis elegans to identify genetic modifiers of propionate sensitivity. Using genome-wide association (GWA) mapping across wild strains, we identify several genomic regions correlated with reduced propionate sensitivity. We find that natural variation in the putative glucuronosyltransferase GLCT-3, a homolog of human B3GAT, partly explains differences in propionate sensitivity in one of these genomic intervals. We demonstrate that loss-of-function alleles in glct-3 render the animals less sensitive to propionate. Additionally, we find that C. elegans has an expansion of the glct gene family, suggesting that the number of members of this family could influence sensitivity to excess propionate. Our findings demonstrate that natural variation in genes that are not directly associated with propionate breakdown can modulate propionate sensitivity. Our study provides a framework for using C. elegans to characterize the contributions of genetic background in models of human inborn errors in metabolism.

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

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

•

Beta-tubulin alleles from parasitic nematodes were introduced into drug-susceptible C. elegans using CRISPR-Cas9.

•

Beta-tubulin variants are strongly selected in albendazole competitive fitness assays.

•

The E198V allele confers a fitness cost in control conditions as compared to other alleles.

The Gene scb-1 Underlies Variation in Caenorhabditis elegans Chemotherapeutic Responses

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

WormQTL2: an interactive platform for systems genetics in Caenorhabditis elegans

Quantitative genetics provides the tools for linking polymorphic loci to trait variation. Linkage analysis of gene expression is an established and widely applied method, leading to the identification of expression quantitative trait loci (eQTLs). (e)QTL detection facilitates the identification and understanding of the underlying molecular components and pathways, yet (e)QTL data access and mining often is a bottleneck. Here, we present WormQTL2, a database and platform for comparative investigations and meta-analyses of published (e)QTL data sets in the model nematode worm C. elegans. WormQTL2 integrates six eQTL studies spanning 11 conditions as well as over 1000 traits from 32 studies and allows experimental results to be compared, reused and extended upon to guide further experiments and conduct systems-genetic analyses. For example, one can easily screen a locus for specific cis-eQTLs that could be linked to variation in other traits, detect gene-by-environment interactions by comparing eQTLs under different conditions, or find correlations between QTL profiles of classical traits and gene expression. WormQTL2 makes data on natural variation in C. elegans and the identified QTLs interactively accessible, allowing studies beyond the original publications. Database URL: www.bioinformatics.nl/WormQTL2/.

Natural Variation and Genetic Determinants of Caenorhabditis elegans Sperm Size

The diversity in sperm shape and size represents a powerful paradigm to understand how selection drives the evolutionary diversification of cell morphology. Experimental work on the sperm biology of the male-hermaphrodite nematode Caenorhabditis elegans has elucidated diverse factors important for sperm fertilization success, including the competitive superiority of larger sperm. Yet despite extensive research, the molecular mechanisms regulating C. elegans sperm size and the genetic basis underlying natural variation in sperm size remain unknown. To address these questions, we quantified male sperm size variation of a worldwide panel of 97 genetically distinct C. elegans strains, allowing us to uncover significant genetic variation in male sperm size. Aiming to characterize the molecular genetic basis of C. elegans male sperm size variation using a genome-wide association study, we did not detect any significant quantitative trait loci. We therefore focused on the genetic analysis of pronounced sperm size differences observed between recently diverged laboratory strains (N2 vs. LSJ1/2). Using mutants and quantitative complementation tests, we demonstrate that variation in the gene nurf-1 underlies the evolution of small sperm in the LSJ lineage. Given the previous discovery that this same nurf-1 variation was central for hermaphrodite laboratory adaptation, the evolution of reduced male sperm size in LSJ strains likely reflects a pleiotropic consequence. Together, our results provide a comprehensive quantification of natural variation in C. elegans sperm size and first insights into the genetic determinants of Caenorhabditis sperm size, pointing at an involvement of the NURF chromatin remodeling complex.

Tightly linked antagonistic-effect loci underlie polygenic phenotypic variation in C. elegans

Recent work has provided strong empirical support for the classic polygenic model for trait variation. Population-based findings suggest that most regions of genome harbor variation affecting most traits. Here, we use the approach of experimental genetics to show that, indeed, most genomic regions carry variants with detectable effects on growth and reproduction in Caenorhabditis elegans populations sensitized by nickel stress. Nine of 15 adjacent intervals on the X chromosome, each encompassing ∼0.001 of the genome, have significant effects when tested individually in near-isogenic lines (NILs). These intervals have effects that are similar in magnitude to those of genome-wide significant loci that we mapped in a panel of recombinant inbred advanced intercross lines (RIAILs). If NIL-like effects were randomly distributed across the genome, the RIAILs would exhibit phenotypic variance that far exceeds the observed variance. However, the NIL intervals are arranged in a pattern that significantly reduces phenotypic variance relative to a random arrangement; adjacent intervals antagonize one another, cancelling each other's effects. Contrary to the expectation of small additive effects, our findings point to large-effect variants whose effects are masked by epistasis or linkage disequilibrium between alleles of opposing effect.

A Novel Gene Underlies Bleomycin-Response Variation in Caenorhabditis elegans

Bleomycin is a powerful chemotherapeutic drug used to treat a variety of cancers. However, individual patients vary in their responses to bleomycin. The identification of genetic differences that underlie this response variation could improve treatment outcomes by tailoring bleomycin dosages to each patient. We used the model organism Caenorhabditis elegans to identify genetic determinants of bleomycin-response differences by performing linkage mapping on recombinants derived from a cross between the laboratory strain (N2) and a wild strain (CB4856). This approach identified a small genomic region on chromosome V that underlies bleomycin-response variation. Using near-isogenic lines, and strains with CRISPR-Cas9 mediated deletions and allele replacements, we discovered that a novel nematode-specific gene (scb-1) is required for bleomycin resistance. Although the mechanism by which this gene causes variation in bleomycin responses is unknown, we suggest that a rare variant present in the CB4856 strain might cause differences in the potential stress-response function of scb-1 between the N2 and CB4856 strains, thereby leading to differences in bleomycin resistance.

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

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).

Analysis of Epistasis in Natural Traits Using Model Organisms

The ability to detect and understand epistasis in natural populations is important for understanding how biological traits are influenced by genetic variation. However, identification and characterization of epistasis in natural populations remains difficult due to statistical issues that arise as a result of multiple comparisons, and the fact that most genetic variants segregate at low allele frequencies. In this review, we discuss how model organisms may be used to manipulate genotypic combinations to power the detection of epistasis as well as test interactions between specific genes. Findings from a number of species indicate that statistical epistasis is pervasive between natural genetic variants. However, the properties of experimental systems that enable analysis of epistasis also constrain extrapolation of these results back into natural populations.

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

Benzimidazoles (BZ) are essential components of the limited chemotherapeutic arsenal available to control the global burden of parasitic nematodes. The emerging threat of BZ resistance among multiple nematode species necessitates the development of novel strategies to identify genetic and molecular mechanisms underlying this resistance. All detection of parasitic helminth resistance to BZ is focused on the genotyping of three variant sites in the orthologs of the β-tubulin gene found to confer resistance in the free-living nematode Caenorhabditis elegans. Because of the limitations of laboratory and field experiments in parasitic nematodes, it is difficult to look beyond these three sites to identify additional mechanisms that might contribute to BZ resistance in the field. Here, we took an unbiased genome-wide mapping approach in the free-living nematode species C. elegans to identify the genetic underpinnings of natural resistance to the commonly used BZ, albendazole (ABZ). We found a wide range of natural variation in ABZ resistance in natural C. elegans populations. In agreement with known mechanisms of BZ resistance in parasites, we found that a majority of the variation in ABZ resistance among wild C. elegans strains is caused by variation in the β-tubulin gene ben-1. This result shows empirically that resistance to ABZ naturally exists and segregates within the C. elegans population, suggesting that selection in natural niches could enrich for resistant alleles. We identified 25 distinct ben-1 alleles that are segregating at low frequencies within the C. elegans population, including many novel molecular variants. Population genetic analyses indicate that ben-1 variation arose multiple times during the evolutionary history of C. elegans and provide evidence that these alleles likely occurred recently because of local selective pressures. Additionally, we find purifying selection at all five β-tubulin genes, despite predicted loss-of-function variants in ben-1, indicating that BZ resistance in natural niches is a stronger selective pressure than loss of one β-tubulin gene. Furthermore, we used genome-editing to show that the most common parasitic nematode β-tubulin allele that confers BZ resistance, F200Y, confers resistance in C. elegans. Importantly, we identified a novel genomic region that is correlated with ABZ resistance in the C. elegans population but independent of ben-1 and the other β-tubulin loci, suggesting that there are multiple mechanisms underlying BZ resistance. Taken together, our results establish a population-level resource of nematode natural diversity as an important model for the study of mechanisms that give rise to BZ resistance.

Shared Genomic Regions Underlie Natural Variation in Diverse Toxin Responses

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

From "the Worm" to "the Worms" and Back Again: The Evolutionary Developmental Biology of Nematodes

Since the earliest days of research on nematodes, scientists have noted the developmental and morphological variation that exists within and between species. As various cellular and developmental processes were revealed through intense focus on Caenorhabditis elegans, these comparative studies have expanded. Within the genus Caenorhabditis, they include characterization of intraspecific polymorphisms and comparisons of distinct species, all generally amenable to the same laboratory culture methods and supported by robust genomic and experimental tools. The C. elegans paradigm has also motivated studies with more distantly related nematodes and animals. Combined with improved phylogenies, this work has led to important insights about the evolution of nematode development. First, while many aspects of C. elegans development are representative of Caenorhabditis, and of terrestrial nematodes more generally, others vary in ways both obvious and cryptic. Second, the system has revealed several clear examples of developmental flexibility in achieving a particular trait. This includes developmental system drift, in which the developmental control of homologous traits has diverged in different lineages, and cases of convergent evolution. Overall, the wealth of information and experimental techniques developed in C. elegans is being leveraged to make nematodes a powerful system for evolutionary cellular and developmental biology.

Acute Effects of Drugs on Caenorhabditis elegans Movement Reveal Complex Responses and Plasticity

Many drugs act very rapidly - they can turn on or off their targets within minutes in a whole animal. What are the acute effects of drug treatment and how does an animal respond to these? We developed a simple assay to measure the acute effects of drugs on C. elegans movement and examined the effects of a range of compounds including neuroactive drugs, toxins, environmental stresses and novel compounds on worm movement over a time period of 3 hr. We found a wide variety of acute responses. Many compounds cause rapid paralysis which may be permanent or followed by one or more recovery phases. The recoveries are not the result of some generic stress response but are specific to the drug e.g., recovery from paralysis due to a neuroactive drug requires neurotransmitter pathways whereas recovery from a metabolic inhibitor requires metabolic changes. Finally, we also find that acute responses can vary greatly across development and that there is extensive natural variation in acute responses. In summary, acute responses are sensitive probes of the ability of biological networks to respond to drug treatment and these responses can reveal the action of unexplored pathways.

Caenorhabditis elegans as an emerging model system in environmental epigenetics

The roundworm Caenorhabitis elegans has been an established model organism for the study of genetics and developmental biology, including studies of transcriptional regulation, since the 1970s. This model organism has continued to be used as a classical model system as the field of transcriptional regulation has expanded to include scientific advances in epigenetics and chromatin biology. In the last several decades, C. elegans has emerged as a powerful model for environmental toxicology, particularly for the study of chemical genotoxicity. Here, we outline the utility and applicability of C. elegans as a powerful model organism for mechanistic studies of environmental influences on the epigenome. Our goal in this article is to inform the field of environmental epigenetics of the strengths and limitations of the well-established C. elegans model organism as an emerging model for medium-throughput, in vivo exploration of the role of exogenous chemical stimuli in transcriptional regulation, developmental epigenetic reprogramming, and epigenetic memory and inheritance. As the field of environmental epigenetics matures, and research begins to map mechanisms underlying observed associations, new toolkits and model systems, particularly manipulable, scalable in vivo systems that accurately model human transcriptional regulatory circuits, will provide an essential experimental bridge between in vitro biochemical experiments and mammalian model systems. Environ. Mol. Mutagen. 59:560-575, 2018. © 2018 Wiley Periodicals, Inc.

Polyamine-containing etoposide derivatives as poisons of human type II topoisomerases: Differential effects on topoisomerase IIα and IIβ

Etoposide is an anticancer drug that acts by inducing topoisomerase II-mediated DNA cleavage. Despite its wide use, etoposide is associated with some very serious side-effects including the development of treatment-related acute myelogenous leukemias. Etoposide targets both human topoisomerase IIα and IIβ. However, the contributions of the two enzyme isoforms to the therapeutic vs. leukemogenic properties of the drug are unclear. In order to develop an etoposide-based drug with specificity for cancer cells that express an active polyamine transport system, the sugar moiety of the drug has been replaced with a polyamine tail. To analyze the effects of this substitution on the specificity of hybrid molecules toward the two enzyme isoforms, we analyzed the activity of a series of etoposide-polyamine hybrids toward human topoisomerase IIα and IIβ. All of the compounds displayed an ability to induce enzyme-mediated DNA cleavage that was comparable to or higher than that of etoposide. Relative to the parent drug, the hybrid compounds displayed substantially higher activity toward topoisomerase IIβ than IIα. Modeling studies suggest that the enhanced specificity may result from interactions with Gln778 in topoisomerase IIβ. The corresponding residue in the α isoform is a methionine.

Topoisomerase II Poisons for Glioblastoma; Existing Challenges and Opportunities to Personalize Therapy

Despite advances in surgery, radiotherapy, and chemotherapy, glioblastoma (GBM) remains a malignancy with poor prognosis. The molecular profile of GBM is diverse across patients, and individual responses to therapy are highly variable. Yet, patients diagnosed with GBM are treated with a rather uniform paradigm. Exploiting these molecular differences and inter-individual responses to therapy may present an opportunity to improve the otherwise bleak prognosis of patients with GBM. This review aims to examine one group of chemotherapeutics: Topoisomerase 2 (TOP2) poisons, a class of drugs that enables TOP2 to induce DNA damage, but interferes with its ability to repair it. These potent chemotherapeutic agents are currently used for a number of malignancies and have shown promise in the treatment of GBM. Despite their robust efficacy in vitro, some of these agents have fallen short of achieving similar results in clinical trials for this tumor. In this review, we explore reasons for this discrepancy, focusing on drug delivery and individual susceptibility differences as challenges for effective TOP2-targeting for GBM. We critically review the evidence implicating genes in susceptibility to TOP2 poisons and categorize this evidence as experimental, correlative or both. This is important as mere experimental evidence does not necessarily lead to identification of genes that serve as good biomarkers of susceptibility for personalizing the use of these drugs.

Discovery of genomic intervals that underlie nematode responses to benzimidazoles

Parasitic nematodes impose a debilitating health and economic burden across much of the world. Nematode resistance to anthelmintic drugs threatens parasite control efforts in both human and veterinary medicine. Despite this threat, the genetic landscape of potential resistance mechanisms to these critical drugs remains largely unexplored. Here, we exploit natural variation in the model nematodes Caenorhabditis elegans and Caenorhabditis briggsae to discover quantitative trait loci (QTL) that control sensitivity to benzimidazoles widely used in human and animal medicine. High-throughput phenotyping of albendazole, fenbendazole, mebendazole, and thiabendazole responses in panels of recombinant lines led to the discovery of over 15 QTL in C. elegans and four QTL in C. briggsae associated with divergent responses to these anthelmintics. Many of these QTL are conserved across benzimidazole derivatives, but others show drug and dose specificity. We used near-isogenic lines to recapitulate and narrow the C. elegans albendazole QTL of largest effect and identified candidate variants correlated with the resistance phenotype. These QTL do not overlap with known benzimidazole target resistance genes from parasitic nematodes and present specific new leads for the discovery of novel mechanisms of nematode benzimidazole resistance. Analyses of orthologous genes reveal conservation of candidate benzimidazole resistance genes in medically important parasitic nematodes. These data provide a basis for extending these approaches to other anthelmintic drug classes and a pathway towards validating new markers for anthelmintic resistance that can be deployed to improve parasite disease control.

The treatment of roundworm (nematode) infections in both humans and animals relies on a small number of anti-parasitic drugs. Resistance to these drugs has appeared in veterinary parasite populations and is a growing concern in human medicine. A better understanding of the genetic basis for parasite drug resistance can be used to help maintain the effectiveness of anti-parasitic drugs and to slow or to prevent the spread of drug resistance in parasite populations. This goal is hampered by the experimental intractability of nematode parasites. Here, we use non-parasitic model nematodes to systematically explore responses to the critical benzimidazole class of anti-parasitic compounds. Using a quantitative genetics approach, we discovered unique genomic intervals that control drug effects, and we identified differences in the effects of these intervals across compounds and doses. We were able to narrow a major-effect genomic region associated with albendazole resistance and to establish that candidate genes discovered in our genetic mappings are largely conserved in important human and animal parasites. This work provides new leads for understanding parasite drug resistance and contributes a powerful template that can be extended to other anti-parasitic drug classes.

Natural diversity facilitates the discovery of conserved chemotherapeutic response mechanisms

Organismal fitness depends on adaptation to complex niches where chemical compounds and pathogens are omnipresent. These stresses can lead to the fixation of alleles in both xenobiotic responses and proliferative signaling pathways that promote survival in these niches. However, both xenobiotic responses and proliferative pathways vary within and among species. For example, genetic differences can accumulate within populations because xenobiotic exposures are not constant and selection is variable. Additionally, neutral genetic variation can accumulate in conserved proliferative pathway genes because these systems are robust to genetic perturbations given their essential roles in normal cell-fate specification. For these reasons, sensitizing mutations or chemical perturbations can disrupt pathways and reveal cryptic variation. With this fundamental view of how organisms respond to cytotoxic compounds and cryptic variation in conserved signaling pathways, it is not surprising that human patients have highly variable responses to chemotherapeutic compounds. These different responses result in the low FDA-approval rates for chemotherapeutics and underscore the need for new approaches to understand these diseases and therapeutic interventions. Model organisms, especially the classic invertebrate systems of Caenorhabditis elegans and Drosophila melanogaster, can be used to combine studies of natural variation across populations with responses to both xenobiotic compounds and chemotherapeutics targeted to conserved proliferative signaling pathways.