The effect of developmental variation on expression QTLs in a multi parental Caenorhabditis elegans population

Regulation of gene expression plays a crucial role in developmental processes and adaptation to changing environments. expression quantitative trait locus (eQTL) mapping is a technique used to study the genetic regulation of gene expression using the transcriptomes of recombinant inbred lines (RILs). Typically, the age of the inbred lines at the time of RNA sampling is carefully controlled. This is necessary because the developmental process causes changes in gene expression, complicating the interpretation of eQTL mapping experiments. However, due to genetics and variation in ambient micro-environments, organisms can differ in their "developmental age," even if they are of the same chronological age. As a result, eQTL patterns are affected by developmental variation in gene expression. The model organism Caenorhabditis elegans is particularly suited for studying the effect of developmental variation on eQTL mapping patterns. In a span of days, C. elegans transitions from embryo through 4 larval stages to adult while undergoing massive changes to its transcriptome. Here, we use C. elegans to investigate the effect of developmental age variation on eQTL patterns and present a normalization procedure. We used dynamical eQTL mapping, which includes the developmental age as a cofactor, to separate the variation in development from genotypic variation and explain variation in gene expression levels. We compare classical single marker eQTL mapping and dynamical eQTL mapping using RNA-seq data of ∼200 multi-parental RILs of C. elegans. The results show that (1) many eQTLs are caused by developmental variation, (2) most trans-bands are developmental QTLs, and (3) dynamical eQTL mapping detects additional eQTLs not found with classical eQTL mapping. We recommend that correction for variation in developmental age should be strongly considered in eQTL mapping studies given the large impact of processes like development on the transcriptome.

Natural allelic variation modifies acute ethanol response phenotypes in wild strains of C. elegans

BACKGROUND

Genetic variation contributes to the likelihood that an individual will develop an alcohol use disorder (AUD). Traditional laboratory studies in animal models have elucidated the molecular pharmacology of ethanol, but laboratory-derived genetic manipulations rarely model the naturally occurring genetic variation observed in wild populations. Rather, these manipulations are biased toward identifying genes of central importance in the phenotypes. Because changes in such genes can confer selective disadvantages, they are not ideal candidates for carrying AUD risk alleles in humans. We sought to exploit Caenorhabditis elegans to identify allelic variation existing in the wild that modulates ethanol response behaviors.

METHODS

We tested the acute ethanol responses of four strains recently isolated from the wild (JU1511, JU1926, JU1931, and JU1941) and 41 multiparental recombinant inbred lines (mpRILs) derived from them. We assessed locomotion at 10, 30, and 50 min on low and high ethanol concentrations. We performed principal component analyses (PCA) on the different phenotypes, tested for transgressive behavior, calculated heritability, and determined the correlations between behavioral responses.

RESULTS

We observed a range of responses to ethanol across the strains. We detected a low-concentration locomotor activation effect in some of the mpRILs not seen in the laboratory wild-type strain. PCA showed different ethanol response behaviors to be independent. We observed transgressive behavior for many of the measured phenotypes and found that multiple behaviors were uncorrelated. The average broad-sense heritability for all phenotypes was 23.2%.

CONCLUSIONS

Genetic variation significantly affects multiple acute ethanol response behaviors, many of which are independent of one another. This suggests that the genetic variation captured by these strains likely affects multiple biological mechanisms through which ethanol acts. Further study of these strains may allow these distinct mechanisms to be identified.

Cryptic genetic variation of eQTL architecture revealed by genetic perturbation in C. elegans

Genetic perturbation in different genetic backgrounds can cause a range of phenotypes within a species. These phenotypic differences can be the result of the interaction between the genetic background and the perturbation. Previously we reported that perturbation of gld-1, an important player in developmental control of C. elegans, released cryptic genetic variation affecting fitness in different genetic backgrounds. Here we investigated the change in transcriptional architecture. We found 414 genes with a cis-eQTL and 991 genes with a trans-eQTL that were specifically found in the gld-1 RNAi treatment. In total, we detected 16 eQTL-hotspots, of which 7 were only found in the gld-1 RNAi treatment. Enrichment analysis of those 7 hotspots showed that the regulated genes were associated with neurons and the pharynx. Furthermore, we found evidence of accelerated transcriptional aging in the gld-1 RNAi treated nematodes. Overall, our results illustrate that studying CGV leads to the discovery of hidden polymorphic regulators.

Differential expression of genes in C. elegans reveals transcriptional responses to indirect-acting xenobiotic compounds and insensitivity to 2,3,7,8-tetrachlorodibenzodioxin

Caenorhabditis elegans is a well-established model organism for toxicity testing of chemical substances. We recently demonstrated its potential for bioanalysis of the toxic potency of chemical contaminants in water. While many detoxification genes are homologues to those in mammalians, C. elegans is reported to be deficient in cytochrome CYP1-like P450 metabolism and that its aryl hydrocarbon receptor (AhR) homolog encoded by ahr-1 purportedly does not interact with dioxins or any other known xenobiotic ligand. This suggests that C. elegans is insensitive for compounds that require bioactivation (indirectly acting compounds) and for dioxins or dioxin-like compounds. This study analysed genome-wide gene expression of the nematode in response to 30 μM of aflatoxin B1 (AFB1), benzo(a)pyrene (B(a)P), Aroclor 1254 (PCB1254), and 10 μM of 2,3,7,8-tetrachlorodibenzodioxin (TCDD). After 24 h of exposure in the early L4 larval stage, microarray analysis revealed 182, 86, and 321 differentially expressed genes in the nematodes treated with 30 μM of AFB1, B(a)P, and PCB1254, respectively. Among these genes, many encode xenobiotic-metabolizing enzymes, and their transcription levels were among the highest-ranked fold-changed genes. Interestingly, only one gene (F59B1.8) was upregulated in the nematodes exposed to 10 μM TCDD. Genes related to metabolic processes and catalytic activity were the most induced by exposure to 30 μM of AFB1, B(a)P, and PCB1254. Despite the genotoxic nature of AFB1 and B(a)P, no differential expression was found in the genes encoding DNA repair and cell cycle checkpoint proteins. Analysis of concentration-response curves was performed to determine the Lowest Observed Transcriptomic Effect Levels (LOTEL) of AFB1, B(a)P, and PCB1254. The obtained LOTEL values showed that gene expression changes in C. elegans are more sensitive to toxicants than reproductive effects. Overall, transcriptional responses of metabolic enzymes suggest that the nematode does metabolize AFB1, B(a)P, and PCB1254. Our findings also support the assumption that the transcription factor AhR homolog in C. elegans does not bind typical xenobiotic ligands, rendering the nematode transcriptionally insensitive to TCDD effects.

Balancing Selection of the Intracellular Pathogen Response in Natural Caenorhabditis elegans Populations

Genetic variation in host populations may lead to differential viral susceptibilities. Here, we investigate the role of natural genetic variation in the Intracellular Pathogen Response (IPR), an important antiviral pathway in the model organism Caenorhabditis elegans against Orsay virus (OrV). The IPR involves transcriptional activity of 80 genes including the pals-genes. We examine the genetic variation in the pals-family for traces of selection and explore the molecular and phenotypic effects of having distinct pals-gene alleles. Genetic analysis of 330 global C. elegans strains reveals that genetic diversity within the IPR-related pals-genes can be categorized in a few haplotypes worldwide. Importantly, two key IPR regulators, pals-22 and pals-25, are in a genomic region carrying signatures of balancing selection, suggesting that different evolutionary strategies exist in IPR regulation. We infected eleven C. elegans strains that represent three distinct pals-22 pals-25 haplotypes with Orsay virus to determine their susceptibility. For two of these strains, N2 and CB4856, the transcriptional response to infection was also measured. The results indicate that pals-22 pals-25 haplotype shapes the defense against OrV and host genetic variation can result in constitutive activation of IPR genes. Our work presents evidence for balancing genetic selection of immunity genes in C. elegans and provides a novel perspective on the functional diversity that can develop within a main antiviral response in natural host populations.

Helminth Glycans at the Host-Parasite Interface and Their Potential for Developing Novel Therapeutics

Helminths are parasitic worms that have successfully co-evolved with their host immune system to sustain long-term infections. Their successful parasitism is mainly facilitated by modulation of the host immune system via the release of excretory-secretory (ES) products covered with glycan motifs such as Lewis X, fucosylated LDN, phosphorylcholine and tyvelose. Evidence is accumulating that these glycans play key roles in different aspects of helminth infection including interactions with immune cells for recognition and evasion of host defences. Moreover, antigenic properties of glycans can be exploited for improving the efficacy of anti-helminthic vaccines. Here, we illustrate that glycans have the potential to open new avenues for the development of novel biopharmaceuticals and effective vaccines based on helminth glycoproteins.

Development of a transcription-based bioanalytical tool to quantify the toxic potencies of hydrophilic compounds in water using the nematode Caenorhabditis elegans

Low concentrations of environmental contaminants can be difficult to detect with current analytical tools, yet they may pose a risk to human and environmental health. The development of bioanalytical tools can help to quantify toxic potencies of biologically active compounds even of hydrophilic contaminants that are hard to extract from water samples. In this study, we exposed the model organism Caenorhabditis elegans synchronized in larval stage L4 to hydrophilic compounds via the water phase and analyzed the effect on gene transcription abundance. The nematodes were exposed to three direct-acting genotoxicants (1 mM and 5 mM): N-ethyl-N-nitrosourea (ENU), formaldehyde (HCHO), and methyl methanesulfonate (MMS). Genome-wide gene expression analysis using microarrays revealed significantly altered transcription levels of 495 genes for HCHO, 285 genes for ENU, and 569 genes for MMS in a concentration-dependent manner. A relatively high number of differentially expressed genes was downregulated, suggesting a general stress in nematodes treated with toxicants. Gene ontology and Kyoto encyclopedia of genes and genomes analysis demonstrated that the upregulated genes were primarily associated with metabolism, xenobiotic detoxification, proteotoxic stress, and innate immune response. Interestingly, genes downregulated by MMS were linked to the inhibition of neurotransmission, and this is in accordance with the observed decreased locomotion in MMS-exposed nematodes. Unexpectedly, the expression level of DNA damage response genes such as cell-cycle checkpoints or DNA-repair proteins were not altered. Overall, the current study shows that gene expression profiling of nematodes can be used to identify the potential mechanisms underlying the toxicity of chemical compounds. C. elegans is a promising test organism to further develop into a bioanalytical tool for quantification of the toxic potency of a wide array of hydrophilic contaminants.

Virus infection modulates male sexual behaviour in Caenorhabditis elegans

Mating dynamics follow from natural selection on mate choice and individuals maximizing their reproductive success. Mate discrimination reveals itself by a plethora of behaviours and morphological characteristics, each of which can be affected by pathogens. A key question is how pathogens affect mate choice and outcrossing behaviour. Here we investigated the effect of Orsay virus on the mating dynamics of the androdiecious (male and hermaphrodite) nematode Caenorhabditis elegans. We tested genetically distinct strains and found that viral susceptibility differed between sexes in a genotype-dependent manner with males of reference strain N2 being more resistant than hermaphrodites. Males displayed a constitutively higher expression of intracellular pathogen response (IPR) genes, whereas the antiviral RNAi response did not have increased activity in males. Subsequent monitoring of sex ratios over 10 generations revealed that viral presence can change mating dynamics in isogenic populations. Sexual attraction assays showed that males preferred mating with uninfected rather than infected hermaphrodites. Together our results illustrate for the first time that viral infection can significantly affect male mating choice and suggest altered mating dynamics as a novel cause benefitting outcrossing under pathogenic stress conditions in C. elegans.

On the role of dauer in the adaptation of nematodes to a parasitic lifestyle

Nematodes are presumably the most abundant Metazoa on Earth, and can even be found in some of the most hostile environments of our planet. Various types of hypobiosis evolved to adapt their life cycles to such harsh environmental conditions. The five most distal major clades of the phylum Nematoda (Clades 8-12), formerly referred to as the Secernentea, contain many economically relevant parasitic nematodes. In this group, a special type of hypobiosis, dauer, has evolved. The dauer signalling pathway, which culminates in the biosynthesis of dafachronic acid (DA), is intensively studied in the free-living nematode Caenorhabditis elegans, and it has been hypothesized that the dauer stage may have been a prerequisite for the evolution of a wide range of parasitic lifestyles among other nematode species. Biosynthesis of DA is not specific for hypobiosis, but if it results in exit of the hypobiotic state, it is one of the main criteria to define certain behaviour as dauer. Within Clades 9 and 10, the involvement of DA has been validated experimentally, and dauer is therefore generally accepted to occur in those clades. However, for other clades, such as Clade 12, this has hardly been explored. In this review, we provide clarity on the nomenclature associated with hypobiosis and dauer across different nematological subfields. We discuss evidence for dauer-like stages in Clades 8 to 12 and support this with a meta-analysis of available genomic data. Furthermore, we discuss indications for a simplified dauer signalling pathway in parasitic nematodes. Finally, we zoom in on the host cues that induce exit from the hypobiotic stage and introduce two hypotheses on how these signals might feed into the dauer signalling pathway for plant-parasitic nematodes. With this work, we contribute to the deeper understanding of the molecular mechanisms underlying hypobiosis in parasitic nematodes. Based on this, novel strategies for the control of parasitic nematodes can be developed.

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.

Heat Stress Reduces the Susceptibility of Caenorhabditis elegans to Orsay Virus Infection

The nematode Caenorhabditis elegans has been a versatile model for understanding the molecular responses to abiotic stress and pathogens. In particular, the response to heat stress and virus infection has been studied in detail. The Orsay virus (OrV) is a natural virus of C. elegans and infection leads to intracellular infection and proteostatic stress, which activates the intracellular pathogen response (IPR). IPR related gene expression is regulated by the genes pals-22 and pals-25, which also control thermotolerance and immunity against other natural pathogens. So far, we have a limited understanding of the molecular responses upon the combined exposure to heat stress and virus infection. We test the hypothesis that the response of C. elegans to OrV infection and heat stress are co-regulated and may affect each other. We conducted a combined heat-stress-virus infection assay and found that after applying heat stress, the susceptibility of C. elegans to OrV was decreased. This difference was found across different wild types of C. elegans. Transcriptome analysis revealed a list of potential candidate genes associated with heat stress and OrV infection. Subsequent mutant screens suggest that pals-22 provides a link between viral response and heat stress, leading to enhanced OrV tolerance of C. elegans after heat stress.

Transcriptional analysis of the response of C. elegans to ethanol exposure

Ethanol-induced transcriptional changes underlie important physiological responses to ethanol that are likely to contribute to the addictive properties of the drug. We examined the transcriptional responses of Caenorhabditis elegans across a timecourse of ethanol exposure, between 30 min and 8 h, to determine what genes and genetic pathways are regulated in response to ethanol in this model. We found that short exposures to ethanol (up to 2 h) induced expression of metabolic enzymes involved in metabolizing ethanol and retinol, while longer exposure (8 h) had much more profound effects on the transcriptome. Several genes that are known to be involved in the physiological response to ethanol, including direct ethanol targets, were regulated at 8 h of exposure. This longer exposure to ethanol also resulted in the regulation of genes involved in cilia function, which is consistent with an important role for the effects of ethanol on cilia in the deleterious effects of chronic ethanol consumption in humans. Finally, we found that food deprivation for an 8-h period induced gene expression changes that were somewhat ameliorated by the presence of ethanol, supporting previous observations that worms can use ethanol as a calorie source.

Punctuated Loci on Chromosome IV Determine Natural Variation in Orsay Virus Susceptibility of Caenorhabditis elegans Strains Bristol N2 and Hawaiian CB4856

Host-pathogen interactions play a major role in evolutionary selection and shape natural genetic variation. The genetically distinct Caenorhabditis elegans strains, Bristol N2 and Hawaiian CB4856, are differentially susceptible to the Orsay virus (OrV). Here, we report the dissection of the genetic architecture of susceptibility to OrV infection. We compare OrV infection in the relatively resistant wild-type CB4856 strain to the more susceptible canonical N2 strain. To gain insight into the genetic architecture of viral susceptibility, 52 fully sequenced recombinant inbred lines (CB4856 × N2 RILs) were exposed to OrV. This led to the identification of two loci on chromosome IV associated with OrV resistance. To verify the two loci and gain additional insight into the genetic architecture controlling virus infection, introgression lines (ILs) that together cover chromosome IV, were exposed to OrV. Of the 27 ILs used, 17 had an CB4856 introgression in an N2 background, and 10 had an N2 introgression in a CB4856 background. Infection of the ILs confirmed and fine-mapped the locus underlying variation in OrV susceptibility, and we found that a single nucleotide polymorphism in cul-6 may contribute to the difference in OrV susceptibility between N2 and CB4856. An allele swap experiment showed the strain CB4856 became as susceptible as the N2 strain by having an N2 cul-6 allele, although having the CB4856 cul-6 allele did not increase resistance in N2. In addition, we found that multiple strains with nonoverlapping introgressions showed a distinct infection phenotype from the parental strain, indicating that there are punctuated locations on chromosome IV determining OrV susceptibility. Thus, our findings reveal the genetic complexity of OrV susceptibility in C. elegans and suggest that viral susceptibility is governed by multiple genes.IMPORTANCE Genetic variation determines the viral susceptibility of hosts. Yet, pinpointing which genetic variants determine viral susceptibility remains challenging. Here, we have exploited the genetic tractability of the model organism Caenorhabditis elegans to dissect the genetic architecture of Orsay virus infection. Our results provide novel insight into natural determinants of Orsay virus infection.

Metabolomics and lipidomics in Caenorhabditis elegans using a single-sample preparation

Comprehensive metabolomic and lipidomic mass spectrometry methods are in increasing demand; for instance, in research related to nutrition and aging. The nematode Caenorhabditis elegans is a key model organism in these fields, owing to the large repository of available C. elegans mutants and their convenient natural lifespan. Here, we describe a robust and sensitive analytical method for the semi-quantitative analysis of >100 polar (metabolomics) and >1000 apolar (lipidomics) metabolites in C. elegans, using a single-sample preparation. Our method is capable of reliably detecting a wide variety of biologically relevant metabolic aberrations in, for example, glycolysis and the tricarboxylic acid cycle, pyrimidine metabolism and complex lipid biosynthesis. In conclusion, we provide a powerful analytical tool that maximizes metabolic data yield from a single sample. This article has an associated First Person interview with the joint first authors of the paper.

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.

Genetic Variation in Caenorhabditis elegans Responses to Pathogenic Microbiota

The bacterivorous nematode Caenorhabditis elegans is an important model species for understanding genetic variation of complex traits. So far, most studies involve axenic laboratory settings using Escherichia coli as the sole bacterial species. Over the past decade, however, investigations into the genetic variation of responses to pathogenic microbiota have increasingly received attention. Quantitative genetic analyses have revealed detailed insight into loci, genetic variants, and pathways in C. elegans underlying interactions with bacteria, microsporidia, and viruses. As various quantitative genetic platforms and resources like C. elegans Natural Diversity Resource (CeNDR) and Worm Quantitative Trait Loci (WormQTL) have been developed, we anticipate that expanding C. elegans research along the lines of genetic variation will be a treasure trove for opening up new insights into genetic pathways and gene functionality of microbiota interactions.

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

Transcriptome resilience predicts thermotolerance in Caenorhabditis elegans

BACKGROUND

The detrimental effects of a short bout of stress can persist and potentially turn lethal, long after the return to normal conditions. Thermotolerance, which is the capacity of an organism to withstand relatively extreme temperatures, is influenced by the response during stress exposure, as well as the recovery process afterwards. While heat-shock response mechanisms have been studied intensively, predicting thermal tolerance remains a challenge.

RESULTS

Here, we use the nematode Caenorhabditis elegans to measure transcriptional resilience to heat stress and predict thermotolerance. Using principal component analysis in combination with genome-wide gene expression profiles collected in three high-resolution time series during control, heat stress, and recovery conditions, we infer a quantitative scale capturing the extent of stress-induced transcriptome dynamics in a single value. This scale provides a basis for evaluating transcriptome resilience, defined here as the ability to depart from stress-expression dynamics during recovery. Independent replication across multiple highly divergent genotypes reveals that the transcriptional resilience parameter measured after a spike in temperature is quantitatively linked to long-term survival after heat stress.

CONCLUSION

Our findings imply that thermotolerance is an intrinsic property that pre-determines long-term outcome of stress and can be predicted by the transcriptional resilience parameter. Inferring the transcriptional resilience parameters of higher organisms could aid in evaluating rehabilitation strategies after stresses such as disease and trauma.

Genetic background modifies phenotypic and transcriptional responses in a C. elegans model of α-synuclein toxicity

Background

Accumulation of protein aggregates are a major hallmark of progressive neurodegenerative disorders such as Parkinson’s disease and Alzheimer’s disease. Transgenic Caenorhabditis elegans nematodes expressing the human synaptic protein α-synuclein in body wall muscle show inclusions of aggregated protein, which affects similar genetic pathways as in humans. It is not however known how the effects of α-synuclein expression in C. elegans differs among genetic backgrounds. Here, we compared gene expression patterns and investigated the phenotypic consequences of transgenic α-synuclein expression in five different C. elegans genetic backgrounds.

Results

Transcriptome analysis indicates that α-synuclein expression effects pathways associated with nutrient storage, lipid transportation and ion exchange and that effects vary depending on the genetic background. These gene expression changes predict that a range of phenotypes will be affected by α-synuclein expression. We confirm this, showing that α-synuclein expression delayed development, reduced lifespan, increased rate of matricidal hatching, and slows pharyngeal pumping. Critically, these phenotypic effects depend on the genetic background and coincide with the core changes in gene expression.

Conclusions

Together, our results show genotype-specific effects and core alterations in both gene expression and in phenotype in response to α-synuclein expression. We conclude that the effects of α-synuclein expression are substantially modified by the genetic background, illustrating that genetic background needs to be considered in C. elegans models of neurodegenerative disease.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5597-1) contains supplementary material, which is available to authorized users.

A multi-parent recombinant inbred line population of C. elegans allows identification of novel QTLs for complex life history traits

BACKGROUND

The nematode Caenorhabditis elegans has been extensively used to explore the relationships between complex traits, genotypes, and environments. Complex traits can vary across different genotypes of a species, and the genetic regulators of trait variation can be mapped on the genome using quantitative trait locus (QTL) analysis of recombinant inbred lines (RILs) derived from genetically and phenotypically divergent parents. Most RILs have been derived from crossing two parents from globally distant locations. However, the genetic diversity between local C. elegans populations can be as diverse as between global populations and could thus provide means of identifying genetic variation associated with complex traits relevant on a broader scale.

RESULTS

To investigate the effect of local genetic variation on heritable traits, we developed a new RIL population derived from 4 parental wild isolates collected from 2 closely located sites in France: Orsay and Santeuil. We crossed these 4 genetically diverse parental isolates to generate a population of 200 multi-parental RILs and used RNA-seq to obtain sequence polymorphisms identifying almost 9000 SNPs variable between the 4 genotypes with an average spacing of 11 kb, doubling the mapping resolution relative to currently available RIL panels for many loci. The SNPs were used to construct a genetic map to facilitate QTL analysis. We measured life history traits such as lifespan, stress resistance, developmental speed, and population growth in different environments, and found substantial variation for most traits. We detected multiple QTLs for most traits, including novel QTLs not found in previous QTL analysis, including those for lifespan and pathogen responses. This shows that recombining genetic variation across C. elegans populations that are in geographical close proximity provides ample variation for QTL mapping.

CONCLUSION

Taken together, we show that using more parents than the classical two parental genotypes to construct a RIL population facilitates the detection of QTLs and that the use of wild isolates facilitates the detection of QTLs. The use of multi-parent RIL populations can further enhance our understanding of local adaptation and life history trade-offs.

Natural genetic variation in C. elegans identified genomic loci controlling metabolite levels

Metabolic homeostasis is sustained by complex biological networks that respond to nutrient availability. Genetic and environmental factors may disrupt this equilibrium, leading to metabolic disorders, including obesity and type 2 diabetes. To identify the genetic factors controlling metabolism, we performed quantitative genetic analysis using a population of 199 recombinant inbred lines (RILs) in the nematode Caenorhabditis elegans We focused on the genomic regions that control metabolite levels by measuring fatty acid (FA) and amino acid (AA) composition in the RILs using targeted metabolomics. The genetically diverse RILs showed a large variation in their FA and AA levels with a heritability ranging from 32% to 82%. We detected strongly co-correlated metabolite clusters and 36 significant metabolite quantitative trait loci (mQTL). We focused on mQTL displaying highly significant linkage and heritability, including an mQTL for the FA C14:1 on Chromosome I, and another mQTL for the FA C18:2 on Chromosome IV. Using introgression lines (ILs), we were able to narrow down both mQTL to a 1.4-Mbp and a 3.6-Mbp region, respectively. RNAi-based screening focusing on the Chromosome I mQTL identified several candidate genes for the C14:1 mQTL, including lagr-1, Y87G2A.2, nhr-265, nhr-276, and nhr-81 Overall, this systems approach provides us with a powerful platform to study the genetic basis of C. elegans metabolism. Furthermore, it allows us to investigate interventions such as nutrients and stresses that maintain or disturb the regulatory network controlling metabolic homeostasis, and identify gene-by-environment interactions.

Genome-wide association mapping of the architecture of susceptibility to the root-knot nematode Meloidogyne incognita in Arabidopsis thaliana

Susceptibility to the root-knot nematode Meloidogyne incognita in plants is thought to be a complex trait based on multiple genes involved in cell differentiation, growth and defence. Previous genetic analyses of susceptibility to M. incognita have mainly focused on segregating dominant resistance genes in crops. It is not known if plants harbour significant genetic variation in susceptibility to M. incognita independent of dominant resistance. To study the genetic architecture of susceptibility to M. incognita, we analysed nematode reproduction on a highly diverse set of 340 natural inbred lines of Arabidopsis thaliana with genome-wide association mapping. We observed a surprisingly large variation in nematode reproduction among these lines. Genome-wide association mapping revealed four quantitative trait loci (QTLs) located on chromosomes 1 and 5 of A. thaliana significantly associated with reproductive success of M. incognita, none of which harbours typical resistance gene homologues. Mutant analysis of three genes located in two QTLs showed that the transcription factor BRASSINAZOLE RESISTANT1 and an F-box family protein may function as (co-)regulators of susceptibility to M. incognita in Arabidopsis. Our data suggest that breeding for loss-of-susceptibility, based on allelic variants critically involved in nematode feeding, could be used to make crops more resilient to root-knot nematodes.

Rictor/TORC2 mediates gut-to-brain signaling in the regulation of phenotypic plasticity in C. elegans

Animals integrate external cues with information about internal conditions such as metabolic state to execute the appropriate behavioral and developmental decisions. Information about food quality and quantity is assessed by the intestine and transmitted to modulate neuronal functions via mechanisms that are not fully understood. The conserved Target of Rapamycin complex 2 (TORC2) controls multiple processes in response to cellular stressors and growth factors. Here we show that TORC2 coordinates larval development and adult behaviors in response to environmental cues and feeding state in the bacterivorous nematode C. elegans. During development, pheromone, bacterial food, and temperature regulate expression of the daf-7 TGF-β and daf-28 insulin-like peptide in sensory neurons to promote a binary decision between reproductive growth and entry into the alternate dauer larval stage. We find that TORC2 acts in the intestine to regulate neuronal expression of both daf-7 and daf-28, which together reflect bacterial-diet dependent feeding status, thus providing a mechanism for integration of food signals with external cues in the regulation of neuroendocrine gene expression. In the adult, TORC2 similarly acts in the intestine to modulate food-regulated foraging behaviors via a PDF-2/PDFR-1 neuropeptide signaling-dependent pathway. We also demonstrate that genetic variation affects food-dependent larval and adult phenotypes, and identify quantitative trait loci (QTL) associated with these traits. Together, these results suggest that TORC2 acts as a hub for communication of feeding state information from the gut to the brain, thereby contributing to modulation of neuronal function by internal state.

Decision-making in all animals, including humans, involves weighing available information about the external environment as well as the animals’ internal conditions. Information about the environment is obtained via the sensory nervous system, whereas internal state can be assessed via cues such as levels of hormones or nutrients. How multiple external and internal inputs are processed in the nervous system to drive behavior or development is not fully understood. In this study, we examine how the nematode C. elegans integrates dietary information received by the gut with environmental signals to alter nervous system function. We have found that a signaling complex, called TORC2, acts in the gut to relay nutrition signals to alter hormonal signaling by the nervous system in C. elegans. Altered neuronal signaling in turn affects a food-dependent binary developmental decision in larvae, as well as food-dependent foraging behaviors in adults. Our results provide a mechanism by which animals prioritize specific signals such as feeding status to appropriately alter their development and/or behavior.

Temporal dynamics of gene expression in heat-stressed Caenorhabditis elegans

There is considerable insight into pathways and genes associated with heat-stress conditions. Most genes involved in stress response have been identified using mutant screens or gene knockdowns. Yet, there is limited understanding of the temporal dynamics of global gene expression in stressful environments. Here, we studied global gene expression profiles during 12 hours of heat stress in the nematode C. elegans. Using a high-resolution time series of increasing stress exposures, we found a distinct shift in gene expression patterns between 3–4 hours into the stress response, separating an initially highly dynamic phase from a later relatively stagnant phase. This turning point in expression dynamics coincided with a phenotypic turning point, as shown by a strong decrease in movement, survival and, progeny count in the days following the stress. Both detectable at transcriptional and phenotypic level, this study pin-points a relatively small time frame during heat stress at which enough damage is accumulated, making it impossible to recover the next few days.

Why do Individuals Differ in Viral Susceptibility? A Story Told by Model Organisms

Viral susceptibility and disease progression is determined by host genetic variation that underlies individual differences. Genetic polymorphisms that affect the phenotype upon infection have been well-studied for only a few viruses, such as HIV-1 and Hepatitis C virus. However, even for well-studied viruses the genetic basis of individual susceptibility differences remains elusive. Investigating the effect of causal polymorphisms in humans is complicated, because genetic methods to detect rare or small-effect polymorphisms are limited and genetic manipulation is not possible in human populations. Model organisms have proven a powerful experimental platform to identify and characterize polymorphisms that underlie natural variations in viral susceptibility using quantitative genetic tools. We summarize and compare the genetic tools available in three main model organisms, Mus musculus, Drosophila melanogaster, and Caenorhabditis elegans, and illustrate how these tools can be applied to detect polymorphisms that determine the viral susceptibility. Finally, we analyse how candidate polymorphisms from model organisms can be used to shed light on the underlying mechanism of individual variation. Insights in causal polymorphisms and mechanisms underlying individual differences in viral susceptibility in model organisms likely provide a better understanding in humans.

Forward and reverse genetics approaches to uncover metabolic aging pathways in Caenorhabditis elegans

The biological mechanisms of aging have been studied in depth and prominent findings in this field promote the development of new therapies for age-associated disorders. Various model organisms are used for research on aging; among these, the nematode Caenorhabditis elegans has been widely used and has provided valuable knowledge in determining the regulatory mechanisms driving the aging process. Many genes involved in lifespan regulation are associated with metabolic pathways and are influenced by genetic and environmental factors. In line with this, C. elegans provides a promising platform to study such gene by environment interactions, in either a reverse or forward genetics approach. In this review, we discuss longevity mechanisms related to metabolic networks that have been discovered in C. elegans. We also highlight the use of wild populations to study the complex genetic basis of natural variation for quantitative traits that mediate longevity.

Contribution of trans regulatory eQTL to cryptic genetic variation in C. elegans

Background

Cryptic genetic variation (CGV) is the hidden genetic variation that can be unlocked by perturbing normal conditions. CGV can drive the emergence of novel complex phenotypes through changes in gene expression. Although our theoretical understanding of CGV has thoroughly increased over the past decade, insight into polymorphic gene expression regulation underlying CGV is scarce. Here we investigated the transcriptional architecture of CGV in response to rapid temperature changes in the nematode Caenorhabditis elegans. We analyzed regulatory variation in gene expression (and mapped eQTL) across the course of a heat stress and recovery response in a recombinant inbred population.

Results

We measured gene expression over three temperature treatments: i) control, ii) heat stress, and iii) recovery from heat stress. Compared to control, exposure to heat stress affected the transcription of 3305 genes, whereas 942 were affected in recovering animals. These affected genes were mainly involved in metabolism and reproduction. The gene expression pattern in recovering animals resembled both the control and the heat-stress treatment. We mapped eQTL using the genetic variation of the recombinant inbred population and detected 2626 genes with an eQTL in the heat-stress treatment, 1797 in the control, and 1880 in the recovery. The cis-eQTL were highly shared across treatments. A considerable fraction of the trans-eQTL (40–57%) mapped to 19 treatment specific trans-bands. In contrast to cis-eQTL, trans-eQTL were highly environment specific and thus cryptic. Approximately 67% of the trans-eQTL were only induced in a single treatment, with heat-stress showing the most unique trans-eQTL.

Conclusions

These results illustrate the highly dynamic pattern of CGV across three different environmental conditions that can be evoked by a stress response over a relatively short time-span (2 h) and that CGV is mainly determined by response related trans regulatory eQTL.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3899-8) contains supplementary material, which is available to authorized users.

The background puzzle: how identical mutations in the same gene lead to different disease symptoms

Identical disease-causing mutations can lead to different symptoms in different people. The reason for this has been a puzzling problem for geneticists. Differential penetrance and expressivity of mutations has been observed within individuals with different and similar genetic backgrounds. Attempts have been made to uncover the underlying mechanisms that determine differential phenotypic effects of identical mutations through studies of model organisms. From these studies evidence is accumulating that to understand disease mechanism or predict disease prevalence, an understanding of the influence of genetic background is as important as the putative disease-causing mutations of relatively large effect. This review highlights current insights into phenotypic variation due to gene interactions, epigenetics and stochasticity in model organisms, and discusses their importance for understanding the mutational effect on disease symptoms.

AraQTL - workbench and archive for systems genetics in Arabidopsis thaliana

Genetical genomics studies uncover genome-wide genetic interactions between genes and their transcriptional regulators. High-throughput measurement of gene expression in recombinant inbred line populations has enabled investigation of the genetic architecture of variation in gene expression. This has the potential to enrich our understanding of the molecular mechanisms affected by and underlying natural variation. Moreover, it contributes to the systems biology of natural variation, as a substantial number of experiments have resulted in a valuable amount of interconnectable phenotypic, molecular and genotypic data. A number of genetical genomics studies have been published for Arabidopsis thaliana, uncovering many expression quantitative trait loci (eQTLs). However, these complex data are not easily accessible to the plant research community, leaving most of the valuable genetic interactions unexplored as cross-analysis of these studies is a major effort. We address this problem with AraQTL (http://www.bioinformatics.nl/Ara QTL/), an easily accessible workbench and database for comparative analysis and meta-analysis of all published Arabidopsis eQTL datasets. AraQTL provides a workbench for comparing, re-using and extending upon the results of these experiments. For example, one can easily screen a physical region for specific local eQTLs that could harbour candidate genes for phenotypic QTLs, or detect gene-by-environment interactions by comparing eQTLs under different conditions.

reGenotyper: Detecting mislabeled samples in genetic data

In high-throughput molecular profiling studies, genotype labels can be wrongly assigned at various experimental steps; the resulting mislabeled samples seriously reduce the power to detect the genetic basis of phenotypic variation. We have developed an approach to detect potential mislabeling, recover the “ideal” genotype and identify “best-matched” labels for mislabeled samples. On average, we identified 4% of samples as mislabeled in eight published datasets, highlighting the necessity of applying a “data cleaning” step before standard data analysis.

Contrasting invertebrate immune defense behaviors caused by a single gene, the Caenorhabditis elegans neuropeptide receptor gene npr-1

Background

The invertebrate immune system comprises physiological mechanisms, physical barriers and also behavioral responses. It is generally related to the vertebrate innate immune system and widely believed to provide nonspecific defense against pathogens, whereby the response to different pathogen types is usually mediated by distinct signalling cascades. Recent work suggests that invertebrate immune defense can be more specific at least at the phenotypic level. The underlying genetic mechanisms are as yet poorly understood.

Results

We demonstrate in the model invertebrate Caenorhabditis elegans that a single gene, a homolog of the mammalian neuropeptide Y receptor gene, npr-1, mediates contrasting defense phenotypes towards two distinct pathogens, the Gram-positive Bacillus thuringiensis and the Gram-negative Pseudomonas aeruginosa. Our findings are based on combining quantitative trait loci (QTLs) analysis with functional genetic analysis and RNAseq-based transcriptomics. The QTL analysis focused on behavioral immune defense against B. thuringiensis, using recombinant inbred lines (RILs) and introgression lines (ILs). It revealed several defense QTLs, including one on chromosome X comprising the npr-1 gene. The wildtype N2 allele for the latter QTL was associated with reduced defense against B. thuringiensis and thus produced an opposite phenotype to that previously reported for the N2 npr-1 allele against P. aeruginosa. Analysis of npr-1 mutants confirmed these contrasting immune phenotypes for both avoidance behavior and nematode survival. Subsequent transcriptional profiling of C. elegans wildtype and npr-1 mutant suggested that npr-1 mediates defense against both pathogens through p38 MAPK signaling, insulin-like signaling, and C-type lectins. Importantly, increased defense towards P. aeruginosa seems to be additionally influenced through the induction of oxidative stress genes and activation of GATA transcription factors, while the repression of oxidative stress genes combined with activation of Ebox transcription factors appears to enhance susceptibility to B. thuringiensis.

Conclusions

Our findings highlight the role of a single gene, npr-1, in fine-tuning nematode immune defense, showing the ability of the invertebrate immune system to produce highly specialized and potentially opposing immune responses via single regulatory genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2603-8) contains supplementary material, which is available to authorized users.

Natural Genetic Variation Influences Protein Abundances in C. elegans Developmental Signalling Pathways

Complex traits, including common disease-related traits, are affected by many different genes that function in multiple pathways and networks. The apoptosis, MAPK, Notch, and Wnt signalling pathways play important roles in development and disease progression. At the moment we have a poor understanding of how allelic variation affects gene expression in these pathways at the level of translation. Here we report the effect of natural genetic variation on transcript and protein abundance involved in developmental signalling pathways in Caenorhabditis elegans. We used selected reaction monitoring to analyse proteins from the abovementioned four pathways in a set of recombinant inbred lines (RILs) generated from the wild-type strains N2 (Bristol) and CB4856 (Hawaii) to enable quantitative trait locus (QTL) mapping. About half of the cases from the 44 genes tested showed a statistically significant change in protein abundance between various strains, most of these were however very weak (below 1.3-fold change). We detected a distant QTL on the left arm of chromosome II that affected protein abundance of the phosphatidylserine receptor protein PSR-1, and two separate QTLs that influenced embryonic and ionizing radiation-induced apoptosis on chromosome IV. Our results demonstrate that natural variation in C. elegans is sufficient to cause significant changes in signalling pathways both at the gene expression (transcript and protein abundance) and phenotypic levels.

Natural Genetic Variation Differentially Affects the Proteome and Transcriptome in Caenorhabditis elegans

Natural genetic variation is the raw material of evolution and influences disease development and progression. An important question is how this genetic variation translates into variation in protein abundance. To analyze the effects of the genetic background on gene and protein expression in the nematode Caenorhabditis elegans, we quantitatively compared the two genetically highly divergent wild-type strains N2 and CB4856. Gene expression was analyzed by microarray assays, and proteins were quantified using stable isotope labeling by amino acids in cell culture. Among all transcribed genes, we found 1,532 genes to be differentially transcribed between the two wild types. Of the total 3,238 quantified proteins, 129 proteins were significantly differentially expressed between N2 and CB4856. The differentially expressed proteins were enriched for genes that function in insulin-signaling and stress-response pathways, underlining strong divergence of these pathways in nematodes. The protein abundance of the two wild-type strains correlates more strongly than protein abundance versus transcript abundance within each wild type. Our findings indicate that in C. elegans only a fraction of the changes in protein abundance can be explained by the changes in mRNA abundance. These findings corroborate with the observations made across species.

On predicting regulatory genes by analysis of functional networks in C. elegans

Background

Connectivity networks, which reflect multiple interactions between genes and proteins, possess not only a descriptive but also a predictive value, as new connections can be extrapolated and tested by means of computational analysis. Integration of different types of connectivity data (such as co-expression and genetic interactions) in one network has proven to benefit ‘guilt by association’ analysis. However predictive values of connectives of different types, that had their specific functional meaning and topological characteristics were not obvious, and have been addressed in this analysis.

Methods

eQTL data for 3 experimental C.elegans age groups were retrieved from WormQTL. WormNet has been used to obtain pair-wise gene interactions. The Shortest Path Function (SPF) has been adopted for statistical validation of the co-expressed gene clusters and for computational prediction of their potential gene expression regulators from a network context. A new SPF-based algorithm has been applied to genetic interactions sub-networks adjacent to the clusters of co-expressed genes for ranking the most likely gene expression regulators causal to eQTLs.

Results

We have demonstrated that known co-expression and genetic interactions between C. elegans genes can be complementary in predicting gene expression regulators. Several algorithms were compared in respect to their predictive potential in different network connectivity contexts. We found that genes associated with eQTLs are highly clustered in a C. elegans co-expression sub-network, and their adjacent genetic interactions provide the optimal functional connectivity environment for application of the new SPF-based algorithm. It was successfully tested in the reverse-prediction analysis on groups of genes with known regulators and applied to co-expressed genes and experimentally observed expression quantitative trait loci (eQTLs).

Conclusions

This analysis demonstrates differences in topology and connectivity of co-expression and genetic interactions sub-networks in WormNet. The modularity of less continuous genetic interaction network does not correspond to modularity of the dense network comprised by gene co-expression interactions. However the genetic interaction network can be used much more efficiently with the SPF method in prediction of potential regulators of gene expression. The developed method can be used for validation of functional significance of suggested eQTLs and a discovery of new regulatory modules.

Widespread genomic incompatibilities in Caenorhabditis elegans

In the Bateson-Dobzhansky-Muller (BDM) model of speciation, incompatibilities emerge from the deleterious interactions between alleles that are neutral or advantageous in the original genetic backgrounds, i.e., negative epistatic effects. Within species such interactions are responsible for outbreeding depression and F2 (hybrid) breakdown. We sought to identify BDM incompatibilities in the nematode Caenorhabditis elegans by looking for genomic regions that disrupt egg laying; a complex, highly regulated, and coordinated phenotype. Investigation of introgression lines and recombinant inbred lines derived from the isolates CB4856 and N2 uncovered multiple incompatibility quantitative trait loci (QTL). These QTL produce a synthetic egg-laying defective phenotype not seen in CB4856 and N2 nor in other wild isolates. For two of the QTL regions, results are inconsistent with a model of pairwise interaction between two loci, suggesting that the incompatibilities are a consequence of complex interactions between multiple loci. Analysis of additional life history traits indicates that the QTL regions identified in these screens are associated with effects on other traits such as lifespan and reproduction, suggesting that the incompatibilities are likely to be deleterious. Taken together, these results indicate that numerous BDM incompatibilities that could contribute to reproductive isolation can be detected and mapped within C. elegans.

Worm variation made accessible: Take your shopping cart to store, link, and investigate!

In Caenorhabditis elegans, the recent advances in high-throughput quantitative analyses of natural genetic and phenotypic variation have led to a wealth of data on genotype phenotype relations. This data has resulted in the discovery of genes with major allelic effects and insights in the effect of natural genetic variation on a whole range of complex traits as well as how this variation is distributed across the genome. Regardless of the advances presented in specific studies, the majority of the data generated in these studies had yet to be made easily accessible, allowing for meta-analysis. Not only data in figures or tables but meta-data should be accessible for further investigation and comparison between studies. A platform was created where all the data, phenotypic measurements, genotypes, and mappings can be stored, compared, and new linkages within and between published studies can be discovered. WormQTL focuses on quantitative genetics in Caenorhabditis and other nematode species, whereas WormQTL^HD quantitatively links gene expression quantitative trait loci (eQTL) in C. elegans to gene–disease associations in humans.

Nematode endogenous small RNA pathways

The discovery of small RNA silencing pathways has greatly extended our knowledge of gene regulation. Small RNAs have been presumed to play a role in every field of biology because they affect many biological processes via regulation of gene expression and chromatin remodeling. Most well-known examples of affected processes are development, fertility, and maintenance of genome stability. Here we review the role of the three main endogenous small RNA silencing pathways in Caenorhabditis elegans: microRNAs, endogenous small interfering RNAs, and PIWI-interacting RNAs. After providing an entry-level overview on how these pathways function, we discuss research on other nematode species providing insight into the evolution of these small RNA pathways. In understanding the differences between the endogenous small RNA pathways and their evolution, a more comprehensive picture is formed of the functions and effects of small RNAs.

A heritable antiviral RNAi response limits Orsay virus infection in Caenorhabditis elegans N2

Orsay virus (OrV) is the first virus known to be able to complete a full infection cycle in the model nematode species Caenorhabditis elegans. OrV is transmitted horizontally and its infection is limited by antiviral RNA interference (RNAi). However, we have no insight into the kinetics of OrV replication in C. elegans. We developed an assay that infects worms in liquid, allowing precise monitoring of the infection. The assay revealed a dual role for the RNAi response in limiting Orsay virus infection in C. elegans. Firstly, it limits the progression of the initial infection at the step of recognition of dsRNA. Secondly, it provides an inherited protection against infection in the offspring. This establishes the heritable RNAi response as anti-viral mechanism during OrV infections in C. elegans. Our results further illustrate that the inheritance of the anti-viral response is important in controlling the infection in the canonical wild type Bristol N2. The OrV replication kinetics were established throughout the worm life-cycle, setting a standard for further quantitative assays with the OrV-C. elegans infection model.

WormQTLHD--a web database for linking human disease to natural variation data in C. elegans

Interactions between proteins are highly conserved across species. As a result, the molecular basis of multiple diseases affecting humans can be studied in model organisms that offer many alternative experimental opportunities. One such organism—Caenorhabditis elegans—has been used to produce much molecular quantitative genetics and systems biology data over the past decade. We present WormQTL^HD (Human Disease), a database that quantitatively and systematically links expression Quantitative Trait Loci (eQTL) findings in C. elegans to gene–disease associations in man. WormQTL^HD, available online at http://www.wormqtl-hd.org, is a user-friendly set of tools to reveal functionally coherent, evolutionary conserved gene networks. These can be used to predict novel gene-to-gene associations and the functions of genes underlying the disease of interest. We created a new database that links C. elegans eQTL data sets to human diseases (34 337 gene–disease associations from OMIM, DGA, GWAS Central and NHGRI GWAS Catalogue) based on overlapping sets of orthologous genes associated to phenotypes in these two species. We utilized QTL results, high-throughput molecular phenotypes, classical phenotypes and genotype data covering different developmental stages and environments from WormQTL database. All software is available as open source, built on MOLGENIS and xQTL workbench.

Genetic mapping of variation in dauer larvae development in growing populations of Caenorhabditis elegans

In the nematode Caenorhabditis elegans, the appropriate induction of dauer larvae development within growing populations is likely to be a primary determinant of genotypic fitness. The underlying genetic architecture of natural genetic variation in dauer formation has, however, not been thoroughly investigated. Here, we report extensive natural genetic variation in dauer larvae development within growing populations across multiple wild isolates. Moreover, bin mapping of introgression lines (ILs) derived from the genetically divergent isolates N2 and CB4856 reveals 10 quantitative trait loci (QTLs) affecting dauer formation. Comparison of individual ILs to N2 identifies an additional eight QTLs, and sequential IL analysis reveals six more QTLs. Our results also show that a behavioural, laboratory-derived, mutation controlled by the neuropeptide Y receptor homolog npr-1 can affect dauer larvae development in growing populations. These findings illustrate the complex genetic architecture of variation in dauer larvae formation in C. elegans and may help to understand how the control of variation in dauer larvae development has evolved.

Gene-environment and protein-degradation signatures characterize genomic and phenotypic diversity in wild Caenorhabditis elegans populations

Background

Analyzing and understanding the relationship between genotypes and phenotypes is at the heart of genetics. Research on the nematode Caenorhabditis elegans has been instrumental for unraveling genotype-phenotype relations, and has important implications for understanding the biology of mammals, but almost all studies, including forward and reverse genetic screens, are limited by investigations in only one canonical genotype. This hampers the detection and functional analysis of allelic variants, which play a key role in controlling many complex traits. It is therefore essential to explore the full potential of the natural genetic variation and evolutionary context of the genotype-phenotype map in wild C. elegans populations.

Results

We used multiple wild C. elegans populations freshly isolated from local sites to investigate gene sequence polymorphisms and a multitude of phenotypes including the transcriptome, fitness, and behavioral traits. The genotype, transcriptome, and a number of fitness traits showed a direct link with the original site of the strains. The separation between the isolation sites was prevalent on all chromosomes, but chromosome V was the largest contributor to this variation. These results were supported by a differential food preference of the wild isolates for naturally co-existing bacterial species. Comparing polymorphic genes between the populations with a set of genes extracted from 19 different studies on gene expression in C. elegans exposed to biotic and abiotic factors, such as bacteria, osmotic pressure, and temperature, revealed a significant enrichment for genes involved in gene-environment interactions and protein degradation.

Conclusions

We found that wild C. elegans populations are characterized by gene-environment signatures, and we have unlocked a wealth of genotype-phenotype relations for the first time. Studying natural isolates provides a treasure trove of evidence compared with that unearthed by the current research in C. elegans, which covers only a diminutive part of the myriad of genotype-phenotype relations that are present in the wild.

Worms under stress: C. elegans stress response and its relevance to complex human disease and aging

Many organisms have stress response pathways, components of which share homology with players in complex human disease pathways. Research on stress response in the nematode worm Caenorhabditis elegans has provided detailed insights into the genetic and molecular mechanisms underlying complex human diseases. In this review we focus on four different types of environmental stress responses - heat shock, oxidative stress, hypoxia, and osmotic stress - and on how these can be used to study the genetics of complex human diseases. All four types of responses involve the genetic machineries that underlie a number of complex human diseases such as cancer and neurodegenerative diseases, including Alzheimer's and Parkinson's. We highlight the types of stress response experiments required to detect the genes and pathways underlying human disease and suggest that studying stress biology in worms can be translated to understanding human disease and provide potential targets for drug discovery.

WormQTL--public archive and analysis web portal for natural variation data in Caenorhabditis spp

Here, we present WormQTL (http://www.wormqtl.org), an easily accessible database enabling search, comparative analysis and meta-analysis of all data on variation in Caenorhabditis spp. Over the past decade, Caenorhabditis elegans has become instrumental for molecular quantitative genetics and the systems biology of natural variation. These efforts have resulted in a valuable amount of phenotypic, high-throughput molecular and genotypic data across different developmental worm stages and environments in hundreds of C. elegans strains. WormQTL provides a workbench of analysis tools for genotype-phenotype linkage and association mapping based on but not limited to R/qtl (http://www.rqtl.org). All data can be uploaded and downloaded using simple delimited text or Excel formats and are accessible via a public web user interface for biologists and R statistic and web service interfaces for bioinformaticians, based on open source MOLGENIS and xQTL workbench software. WormQTL welcomes data submissions from other worm researchers.

Genetic variation for stress-response hormesis in C. elegans lifespan

Increased lifespan can be associated with greater resistance to many different stressors, most notably thermal stress. Such hormetic effects have also been found in C. elegans where short-term exposure to heat lengthens the lifespan. Genetic investigations have been carried out using mutation perturbations in a single genotype, the wild type Bristol N2. Yet, induced mutations do not yield insight regarding the natural genetic variation of thermal tolerance and lifespan. We investigated the genetic variation of heat-shock recovery, i.e. hormetic effects on lifespan and associated quantitative trait loci (QTL) in C. elegans. Heat-shock resulted in an 18% lifespan increase in wild type CB4856 whereas N2 did not show a lifespan elongation. Using recombinant inbred lines (RILs) derived from a cross between wild types N2 and CB4856 we found natural variation in stress-response hormesis in lifespan. Approx. 28% of the RILs displayed a hormesis effect in lifespan. We did not find any hormesis effects for total offspring. Across the RILs there was no relation between lifespan and offspring. The ability to recover from heat-shock mapped to a significant QTL on chromosome II which overlapped with a QTL for offspring under heat-shock conditions. The QTL was confirmed by introgressing relatively small CB4856 regions into chromosome II of N2. Our observations show that there is natural variation in hormetic effects on C. elegans lifespan for heat-shock and that this variation is genetically determined.

A fitness assay for comparing RNAi effects across multiple C. elegans genotypes

Background

RNAi technology by feeding of E. coli containing dsRNA in C. elegans has significantly contributed to further our understanding of many different fields, including genetics, molecular biology, developmental biology and functional genomics. Most of this research has been carried out in a single genotype or genetic background. However, RNAi effects in one genotype do not reveal the allelic effects that segregate in natural populations and contribute to phenotypic variation.

Results

Here we present a method that allows for rapidly comparing RNAi effects among diverse genotypes at an improved high throughput rate. It is based on assessing the fitness of a population of worms by measuring the rate at which E. coli is consumed. Critically, we demonstrate the analytical power of this method by QTL mapping the loss of RNAi sensitivity (in the germline) in a recombinant inbred population derived from a cross between Bristol and a natural isolate from Hawaii. Hawaii has lost RNAi sensitivity in the germline. We found that polymorphisms in ppw-1 contribute to this loss of RNAi sensitivity, but that other loci are also likely to be important.

Conclusions

In summary, we have established a fast method that improves the throughput of RNAi in liquid, that generates quantitative data, that is easy to implement in most laboratories, and importantly that enables QTL mapping using RNAi.

Gene expression modifications by temperature-toxicants interactions in Caenorhabditis elegans

Although organophosphorus pesticides (OP) share a common mode of action, there is increased awareness that they elicit a diverse range of gene expression responses. As yet however, there is no clear understanding of these responses and how they interact with ambient environmental conditions. In the present study, we investigated genome-wide gene expression profiles in the nematode Caenorhabditis elegans exposed to two OP, chlorpyrifos and diazinon, in single and combined treatments at different temperatures. Our results show that chlorpyrifos and diazinon induced expression of different genes and that temperature affected the response of detoxification genes to the pesticides. The analysis of transcriptional responses to a combination of chlorpyrifos and diazinon shows interactions between toxicants that affect gene expression. Furthermore, our combined analysis of the transcriptional responses to OP at different temperatures suggests that the combination of OP and high temperatures affect detoxification genes and modified the toxic levels of the pesticides.

Genome-wide gene expression analysis in response to organophosphorus pesticide chlorpyrifos and diazinon in C. elegans

Organophosphorus pesticides (OPs) were originally designed to affect the nervous system by inhibiting the enzyme acetylcholinesterase, an important regulator of the neurotransmitter acetylcholine. Over the past years evidence is mounting that these compounds affect many other processes. Little is known, however, about gene expression responses against OPs in the nematode Caenorhabditis elegans. This is surprising because C. elegans is extensively used as a model species in toxicity studies. To address this question we performed a microarray study in C. elegans which was exposed for 72 hrs to two widely used Ops, chlorpyrifos and diazinon, and a low dose mixture of these two compounds. Our analysis revealed transcriptional responses related to detoxification, stress, innate immunity, and transport and metabolism of lipids in all treatments. We found that for both compounds as well as in the mixture, these processes were regulated by different gene transcripts. Our results illustrate intense, and unexpected crosstalk between gene pathways in response to chlorpyrifos and diazinon in C. elegans.

Genome-wide gene expression regulation as a function of genotype and age in C. elegans

Gene expression becomes more variable with age, and it is widely assumed that this is due to a decrease in expression regulation. But currently there is no understanding how gene expression regulatory patterns progress with age. Here we explored genome-wide gene expression variation and regulatory loci (eQTL) in a population of developing and aging C. elegans recombinant inbred worms. We found almost 900 genes with an eQTL, of which almost half were found to have a genotype-by-age effect ((gxa)eQTL). The total number of eQTL decreased with age, whereas the variation in expression increased. In developing worms, the number of genes with increased expression variation (1282) was similar to the ones with decreased expression variation (1328). In aging worms, the number of genes with increased variation (1772) was nearly five times higher than the number of genes with a decreased expression variation (373). The number of cis-acting eQTL in juveniles decreased by almost 50% in old worms, whereas the number of trans-acting loci decreased by approximately 27%, indicating that cis-regulation becomes relatively less frequent than trans-regulation in aging worms. Of the 373 genes with decreased expression level variation in aging worms, approximately 39% had an eQTL compared with approximately 14% in developing worms. (gxa)eQTL were found for approximately 21% of these genes in aging worms compared with only approximately 6% in developing worms. We highlight three examples of linkages: in young worms (pgp-6), in old worms (daf-16), and throughout life (lips-16). Our findings demonstrate that eQTL patterns are strongly affected by age, and suggest that gene network integrity declines with age.

Rapid divergence of genetic variance-covariance matrix within a natural population

The matrix of genetic variances and covariances (G matrix) represents the genetic architecture of multiple traits sharing developmental and genetic processes and is central for predicting phenotypic evolution. These predictions require that the G matrix be stable. Yet the timescale and conditions promoting G matrix stability in natural populations remain unclear. We studied stability of the G matrix in a 20-year evolution field experiment, where a population of the cosmopolitan parthenogenetic soil nematode Acrobeloides nanus was subjected to drift and divergent selection (benign and stress environments). Selection regime did not influence the level of absolute genetic constraints: under both regimes, two genetic dimensions for three life-history traits were identified. A substantial response to selection in principal components structure and in general matrix pattern was indicated by three statistical methods. G structure was also influenced by drift, with higher divergence under benign conditions. These results show that the G matrix might evolve rapidly in natural populations. The observed high dynamics of G structure probably represents the general feature of asexual species and limits the predictive power of G in phenotypic evolution analyses.