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

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.

Nematode gene annotation by machine-learning-assisted proteotranscriptomics enables proteome-wide evolutionary analysis

Nematodes encompass more than 24,000 described species, which were discovered in almost every ecological habitat, and make up >80% of metazoan taxonomic diversity in soils. The last common ancestor of nematodes is believed to date back to ∼650-750 million years, generating a large and phylogenetically diverse group to be explored. However, for most species high-quality gene annotations are incomprehensive or missing. Combining short-read RNA sequencing with mass spectrometry-based proteomics and machine-learning quality control in an approach called proteotranscriptomics, we improve gene annotations for nine genome-sequenced nematode species and provide new gene annotations for three additional species without genome assemblies. Emphasizing the sensitivity of our methodology, we provide evidence for two hitherto undescribed genes in the model organism Caenorhabditis elegans Extensive phylogenetic systems analysis using this comprehensive proteome annotation provides new insights into evolutionary processes of this metazoan group.

An automated approach to quantify chemotaxis index in C. elegans

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

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.

Preconditioning With Natural Microbiota Strain Ochrobactrum vermis MYb71 Influences Caenorhabditis elegans Behavior

In comparison with the standard monoxenic maintenance in the laboratory, rearing the nematode Caenorhabditis elegans on its natural microbiota improves its fitness and immunity against pathogens. Although C. elegans is known to exhibit choice behavior and pathogen avoidance behavior, little is known about whether C. elegans actively chooses its (beneficial) microbiota and whether the microbiota influences worm behavior. We examined eleven natural C. elegans isolates in a multiple-choice experiment for their choice behavior toward four natural microbiota bacteria and found that microbiota choice varied among C. elegans isolates. The natural C. elegans isolate MY2079 changed its choice behavior toward microbiota isolate Ochrobactrum vermis MYb71 in both multiple-choice and binary-choice experiments, in particular on proliferating bacteria: O. vermis MYb71 was chosen less than other microbiota bacteria or OP50, but only after preconditioning with MYb71. Examining escape behavior and worm fitness on MYb71, we ruled out pathogenicity of MYb71 and consequently learned pathogen avoidance behavior as the main driver of the behavioral change toward MYb71. The change in behavior of C. elegans MY2079 toward microbiota bacterium MYb71 demonstrates how the microbiota influences the worm's choice. These results might give a baseline for future research on host-microbiota interaction in the C. elegans model.

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 fecundity is correlated with species-wide levels of divergence in Caenorhabditis elegans

Life history traits underlie the fitness of organisms and are under strong natural selection. A new mutation that positively impacts a life history trait will likely increase in frequency and become fixed in a population (e.g., a selective sweep). The identification of the beneficial alleles that underlie selective sweeps provides insights into the mechanisms that occurred during the evolution of a species. In the global population of Caenorhabditis elegans, we previously identified selective sweeps that have drastically reduced chromosomal-scale genetic diversity in the species. Here, we measured the fecundity of 121 wild C. elegans strains, including many recently isolated divergent strains from the Hawaiian islands and found that strains with larger swept genomic regions have significantly higher fecundity than strains without evidence of the recent selective sweeps. We used genome-wide association (GWA) mapping to identify three quantitative trait loci (QTL) underlying the fecundity variation. In addition, we mapped previous fecundity data from wild C. elegans strains and C. elegans recombinant inbred advanced intercross lines that were grown in various conditions and detected eight QTL using GWA and linkage mappings. These QTL show the genetic complexity of fecundity across this species. Moreover, the haplotype structure in each GWA QTL region revealed correlations with recent selective sweeps in the C. elegans population. North American and European strains had significantly higher fecundity than most strains from Hawaii, a hypothesized origin of the C. elegans species, suggesting that beneficial alleles that caused increased fecundity could underlie the selective sweeps during the worldwide expansion of C. elegans.

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.

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.

A three-dimensional habitat for C. elegans environmental enrichment

As we learn more about the importance of gene-environment interactions and the effects of environmental enrichment, it becomes evident that minimalistic laboratory conditions can affect gene expression patterns and behaviors of model organisms. In the laboratory, Caenorhabditis elegans is generally cultured on two-dimensional, homogeneous agar plates abundantly covered with axenic bacteria culture as a food source. However, in the wild, this nematode thrives in rotting fruits and plant stems feeding on bacteria and small eukaryotes. This contrast in habitat complexity suggests that studying C. elegans in enriched laboratory conditions can deepen our understanding of its fundamental traits and behaviors. Here, we developed a protocol to create three-dimensional habitable scaffolds for trans-generational culture of C. elegans in the laboratory. Using decellularization and sterilization of fruit tissue, we created an axenic environment that can be navigated throughout and where the microbial environment can be strictly controlled. C. elegans were maintained over generations on this habitat, and showed a clear behavioral bias for the enriched environment. As an initial assessment of behavioral variations, we found that dauer populations in scaffolds exhibit high-frequency, complex nictation behavior including group towering and jumping behavior.

Bacterial diets differentially alter lifespan and healthspan trajectories in C. elegans

Diet is one of the more variable aspects in life due to the variety of options that organisms are exposed to in their natural habitats. In the laboratory, C. elegans are raised on bacterial monocultures, traditionally the E. coli B strain OP50, and spontaneously occurring microbial contaminants are removed to limit experimental variability because diet-including the presence of contaminants-can exert a potent influence over animal physiology. In order to diversify the menu available to culture C. elegans in the lab, we have isolated and cultured three such microbes: Methylobacterium, Xanthomonas, and Sphingomonas. The nutritional composition of these bacterial foods is unique, and when fed to C. elegans, can differentially alter multiple life history traits including development, reproduction, and metabolism. In light of the influence each food source has on specific physiological attributes, we comprehensively assessed the impact of these bacteria on animal health and devised a blueprint for utilizing different food combinations over the lifespan, in order to promote longevity. The expansion of the bacterial food options to use in the laboratory will provide a critical tool to better understand the complexities of bacterial diets and subsequent changes in physiology and gene expression.

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

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.

The Local Coexistence Pattern of Selfing Genotypes in Caenorhabditis elegans Natural Metapopulations

To study the interplay of rare outcrossing and metapopulation structure, we focus on the nematode Caenorhabditis elegans Its remarkably low outcrossing rate is at the extreme end of the spectrum for facultative selfing organisms. At the demographic level, C. elegans natural populations undergo boom and bust dynamics on ephemeral resources, with the dauer diapause larva acting as the dispersal form. Here we investigate the small-scale genetic structure of C. elegans populations in two localities over several years, using 2b restriction-associated DNA sequencing of nearly 1000 individuals. We find a remarkably small number of genome-wide haplotypes, almost exclusively in the homozygous state, confirming the low effective outcrossing rate. Most strikingly, the major haplotypes in a locality remain intact and do not effectively recombine over several years. From the spatial pattern of diversity, we estimate that each subpopulation or deme is seeded by a mean of 3-10 immigrating individuals. Populations are thus formed by clones that compete at two levels, within a subpopulation and at the metapopulation level. We test for the presence of local phenotypic variation in pathogen resistance and dauer larva nictation, which could possibly explain the maintenance of different genotypes by heterogeneous selection in different local environments or lifecycles. This study is the first to address the local spatiotemporal genetic structure of C. elegans on feeding substrates. We conclude that these animals coexist as competing homozygous clones at the smallest population scale as well as in the metapopulation.

Inverse correlation between longevity and developmental rate among wild C. elegans strains

Genetic studies using model organisms have shown that many long-lived mutants display impaired fitness, such as reduced fecundity and delayed development. However, in several wild animals, the association between longevity and fitness does not seem to be inevitable. Thus, the relationship between longevity and fitness in wild organisms remains inconclusive. Here, we determined the correlation between lifespan and fitness, developmental rate and brood size, by using 16 wild-derived C. elegans strains originated from various geographic areas. We found a negative correlation between lifespan and developmental rate. In contrast, we did not find such negative correlation between longevity and developmental rate among the individuals of C. elegans strains. These data imply that polymorphic genetic variants among wild isolates determine resource allocation to longevity and developmental rate.

Ras/MAPK Modifier Loci Revealed by eQTL in Caenorhabditis elegans

The oncogenic Ras/MAPK pathway is evolutionarily conserved across metazoans. Yet, almost all our knowledge on this pathway comes from studies using single genetic backgrounds, whereas mutational effects can be highly background dependent. Therefore, we lack insight in the interplay between genetic backgrounds and the Ras/MAPK-signaling pathway. Here, we used a Caenorhabditis elegans RIL population containing a gain-of-function mutation in the Ras/MAPK-pathway gene let-60 and measured how gene expression regulation is affected by this mutation. We mapped eQTL and found that the majority (∼73%) of the 1516 detected cis-eQTL were not specific for the let-60 mutation, whereas most (∼76%) of the 898 detected trans-eQTL were associated with the let-60 mutation. We detected six eQTL trans-bands specific for the interaction between the genetic background and the mutation, one of which colocalized with the polymorphic Ras/MAPK modifier amx-2. Comparison between transgenic lines expressing allelic variants of amx-2 showed the involvement of amx-2 in 79% of the trans-eQTL for genes mapping to this trans-band. Together, our results have revealed hidden loci affecting Ras/MAPK signaling using sensitized backgrounds in C. elegans. These loci harbor putative polymorphic modifier genes that would not have been detected using mutant screens in single genetic backgrounds.

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.

Genetic variation in neurodegenerative diseases and its accessibility in the model organism Caenorhabditis elegans

BACKGROUND

Neurodegenerative diseases (NGDs) such as Alzheimer's and Parkinson's are debilitating and largely untreatable conditions strongly linked to age. The clinical, neuropathological, and genetic components of NGDs indicate that neurodegeneration is a complex trait determined by multiple genes and by the environment.

MAIN BODY

The symptoms of NGDs differ among individuals due to their genetic background, and this variation affects the onset and progression of NGD and NGD-like states. Such genetic variation affects the molecular and cellular processes underlying NGDs, leading to differential clinical phenotypes. So far, we have a limited understanding of the mechanisms of individual background variation. Here, we consider how variation between genetic backgrounds affects the mechanisms of aging and proteostasis in NGD phenotypes. We discuss how the nematode Caenorhabditis elegans can be used to identify the role of variation between genetic backgrounds. Additionally, we review advances in C. elegans methods that can facilitate the identification of NGD regulators and/or networks.

CONCLUSION

Genetic variation both in disease genes and in regulatory factors that modulate onset and progression of NGDs are incompletely understood. The nematode C. elegans represents a valuable system in which to address such questions.

The Natural Biotic Environment of Caenorhabditis elegans

Organisms evolve in response to their natural environment. Consideration of natural ecological parameters are thus of key importance for our understanding of an organism's biology. Curiously, the natural ecology of the model species Caenorhabditis elegans has long been neglected, even though this nematode has become one of the most intensively studied models in biological research. This lack of interest changed ∼10 yr ago. Since then, an increasing number of studies have focused on the nematode's natural ecology. Yet many unknowns still remain. Here, we provide an overview of the currently available information on the natural environment of C. elegans We focus on the biotic environment, which is usually less predictable and thus can create high selective constraints that are likely to have had a strong impact on C. elegans evolution. This nematode is particularly abundant in microbe-rich environments, especially rotting plant matter such as decomposing fruits and stems. In this environment, it is part of a complex interaction network, which is particularly shaped by a species-rich microbial community. These microbes can be food, part of a beneficial gut microbiome, parasites and pathogens, and possibly competitors. C. elegans is additionally confronted with predators; it interacts with vector organisms that facilitate dispersal to new habitats, and also with competitors for similar food environments, including competitors from congeneric and also the same species. Full appreciation of this nematode's biology warrants further exploration of its natural environment and subsequent integration of this information into the well-established laboratory-based research approaches.

Correlations of Genotype with Climate Parameters Suggest Caenorhabditis elegans Niche Adaptations

Species inhabit a variety of environmental niches, and the adaptation to a particular niche is often controlled by genetic factors, including gene-by-environment interactions. The genes that vary in order to regulate the ability to colonize a niche are often difficult to identify, especially in the context of complex ecological systems and in experimentally uncontrolled natural environments. Quantitative genetic approaches provide an opportunity to investigate correlations between genetic factors and environmental parameters that might define a niche. Previously, we have shown how a collection of 208 whole-genome sequenced wild Caenorhabditis elegans can facilitate association mapping approaches. To correlate climate parameters with the variation found in this collection of wild strains, we used geographic data to exhaustively curate daily weather measurements in short-term (3 month), middle-term (one year), and long-term (three year) durations surrounding the date of strain isolation. These climate parameters were used as quantitative traits in association mapping approaches, where we identified 11 quantitative trait loci (QTL) for three climatic variables: elevation, relative humidity, and average temperature. We then narrowed the genomic interval of interest to identify gene candidates with variants potentially underlying phenotypic differences. Additionally, we performed two-strain competition assays at high and low temperatures to validate a QTL that could underlie adaptation to temperature and found suggestive evidence supporting that hypothesis.

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

Enhancing metaproteomics--The value of models and defined environmental microbial systems

Metaproteomics--the large-scale characterization of the entire protein complement of environmental microbiota at a given point in time--has provided new features to study complex microbial communities in order to unravel these "black boxes." New technical challenges arose that were not an issue for classical proteome analytics before that could be tackled by the application of different model systems. Here, we review different current and future model systems for metaproteome analysis. Following a short introduction to microbial communities and metaproteomics, we introduce model systems for clinical and biotechnological research questions including acid mine drainage, anaerobic digesters, and activated sludge. Model systems are useful to evaluate the challenges encountered within (but not limited to) metaproteomics, including species complexity and coverage, biomass availability, or reliable protein extraction. The implementation of model systems can be considered as a step forward to better understand microbial community responses and ecological functions of single member organisms. In the future, improvements are necessary to fully explore complex environmental systems by metaproteomics.

Collective effects of common SNPs in foraging decisions in Caenorhabditis elegans and an integrative method of identification of candidate genes

Optimal foraging decision is a quantitative flexible behavior, which describes the time at which animals choose to abandon a depleting food supply. The total minor allele content (MAC) in an individual has been shown to correlate with quantitative variations in complex traits. We have studied the role of MAC in the decision to leave a food lawn in recombinant inbred advanced intercross lines (RIAILs) of Caenorhabditis elegans. We found a strong link between MAC and the food lawn leaving rates (Spearman r = 0.4, P = 0.005). We identified 28 genes of unknown functions whose expression levels correlated with both MAC and leaving rates. When examined by RNAi experiments, 8 of 10 tested among the 28 affected leaving rates, whereas only 2 of 9 did among genes that were only associated with leaving rates but not MAC (8/10 vs 2/9, P < 0.05). The results establish a link between MAC and the foraging behavior and identify 8 genes that may play a role in linking MAC with the quantitative nature of the trait. The method of correlations with both MAC and traits may find broad applications in high efficiency identification of target genes for other complex traits in model organisms and humans.

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.

Genotype-dependent lifespan effects in peptone deprived Caenorhabditis elegans

Dietary restriction appears to act as a general non-genetic mechanism that can robustly prolong lifespan. There have however been reports in many systems of cases where restricted food intake either shortens, or does not affect, lifespan. Here we analyze lifespan and the effect of food restriction via deprived peptone levels on lifespan in wild isolates and introgression lines (ILs) of the nematode Caenorhabditis elegans. These analyses identify genetic variation in lifespan, in the effect of this variation in diet on lifespan and also in the likelihood of maternal, matricidal, hatching. Importantly, in the wild isolates and the ILs, we identify genotypes in which peptone deprivation mediated dietary restriction reduces lifespan. We also identify, in recombinant inbred lines, a locus that affects maternal hatching, a phenotype closely linked to dietary restriction in C. elegans. These results indicate that peptone deprivation mediated dietary restriction affects lifespan in C. elegans in a genotype-dependent manner, reducing lifespan in some genotypes. This may operate by a mechanism similar to dietary restriction.

Ten years of life in compost: temporal and spatial variation of North German Caenorhabditis elegans populations

The nematode Caenorhabditis elegans is a central laboratory model system in almost all biological disciplines, yet its natural life history and population biology are largely unexplored. Such information is essential for in-depth understanding of the nematode's biology because its natural ecology provides the context, in which its traits and the underlying molecular mechanisms evolved. We characterized natural phenotypic and genetic variation among North German C. elegans isolates. We used the unique opportunity to compare samples collected 10 years apart from the same compost heap and additionally included recent samples for this and a second site, collected across a 1.5-year period. Our analysis revealed significant population genetic differentiation between locations, across the 10-year time period, but for only one location a trend across the shorter time frame. Significant variation was similarly found for phenotypic traits of likely importance in nature, such as choice behavior and population growth in the presence of pathogens or naturally associated bacteria. Phenotypic variation was significantly influenced by C. elegans genotype, time of isolation, and sampling site. The here studied C. elegans isolates may provide a valuable, genetically variable resource for future dissection of naturally relevant gene functions.

Caenorhabditis evolution in the wild

Recent research has filled many gaps about Caenorhabditis natural history, simultaneously exposing how much remains to be discovered. This awareness now provides means of connecting ecological and evolutionary theory with diverse biological patterns within and among species in terms of adaptation, sexual selection, breeding systems, speciation, and other phenomena. Moreover, the heralded laboratory tractability of C. elegans, and Caenorhabditis species generally, provides a powerful case study for experimental hypothesis testing about evolutionary and ecological processes to levels of detail unparalleled by most study systems. Here, I synthesize pertinent theory with what we know and suspect about Caenorhabditis natural history for salient features of biodiversity, phenotypes, population dynamics, and interactions within and between species. I identify topics of pressing concern to advance Caenorhabditis biology and to study general evolutionary processes, including the key opportunities to tackle problems in dispersal dynamics, competition, and the dimensionality of niche space.

Remarkably Divergent Regions Punctuate the Genome Assembly of the Caenorhabditis elegans Hawaiian Strain CB4856

The Hawaiian strain (CB4856) of Caenorhabditis elegans is one of the most divergent from the canonical laboratory strain N2 and has been widely used in developmental, population, and evolutionary studies. To enhance the utility of the strain, we have generated a draft sequence of the CB4856 genome, exploiting a variety of resources and strategies. When compared against the N2 reference, the CB4856 genome has 327,050 single nucleotide variants (SNVs) and 79,529 insertion–deletion events that result in a total of 3.3 Mb of N2 sequence missing from CB4856 and 1.4 Mb of sequence present in CB4856 but not present in N2. As previously reported, the density of SNVs varies along the chromosomes, with the arms of chromosomes showing greater average variation than the centers. In addition, we find 61 regions totaling 2.8 Mb, distributed across all six chromosomes, which have a greatly elevated SNV density, ranging from 2 to 16% SNVs. A survey of other wild isolates show that the two alternative haplotypes for each region are widely distributed, suggesting they have been maintained by balancing selection over long evolutionary times. These divergent regions contain an abundance of genes from large rapidly evolving families encoding F-box, MATH, BATH, seven-transmembrane G-coupled receptors, and nuclear hormone receptors, suggesting that they provide selective advantages in natural environments. The draft sequence makes available a comprehensive catalog of sequence differences between the CB4856 and N2 strains that will facilitate the molecular dissection of their phenotypic differences. Our work also emphasizes the importance of going beyond simple alignment of reads to a reference genome when assessing differences between genomes.

The laboratory domestication of Caenorhabditis elegans

Model organisms are of great importance to our understanding of basic biology and to making advances in biomedical research. However, the influence of laboratory cultivation on these organisms is underappreciated, and especially how that environment can affect research outcomes. Recent experiments led to insights into how the widely used laboratory reference strain of the nematode Caenorhabditis elegans compares with natural strains. Here we describe potential selective pressures that led to the fixation of laboratory-derived alleles for the genes npr-1, glb-5, and nath-10. These alleles influence a large number of traits, resulting in behaviors that affect experimental interpretations. Furthermore, strong phenotypic effects caused by these laboratory-derived alleles hinder the discovery of natural alleles. We highlight strategies to reduce the influence of laboratory-derived alleles and to harness the full power of C. elegans.

Why we need more ecology for genetic models such as C. elegans

Functional information about the large majority of the genes is still lacking in the classical eukaryotic model species Drosophila melanogaster, Caenorhabditis elegans, and Mus musculus. Because many of these genes are likely to be important in natural settings, considering explicit ecological information should increase our knowledge of gene function. Using C. elegans as an example, we discuss the importance of biotic factors as a driving force in shaping the composition and structure of the nematode genome. We highlight examples for which consideration of ecological information and natural variation have been key to the identification of novel, unexpected gene functions, and use these examples to define future research avenues for the classical genetic model taxa.

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.

Spread and transmission of bacterial pathogens in experimental populations of the nematode Caenorhabditis elegans

Caenorhabditis elegans is frequently used as a model species for the study of bacterial virulence and innate immunity. In recent years, diverse mechanisms contributing to the nematode's immune response to bacterial infection have been discovered. Yet despite growing interest in the biochemical and molecular basis of nematode-bacterium associations, many questions remain about their ecology. Although recent studies have demonstrated that free-living nematodes could act as vectors of opportunistic pathogens in soil, the extent to which worms may contribute to the persistence and spread of these bacteria has not been quantified. We conducted a series of experiments to test whether colonization of and transmission between C. elegans nematodes could enable two opportunistic pathogens (Salmonella enterica and Pseudomonas aeruginosa) to spread on agar plates occupied by Escherichia coli. We monitored the transmission of S. enterica and P. aeruginosa from single infected nematodes to their progeny and measured bacterial loads both within worms and on the plates. In particular, we analyzed three factors affecting the dynamics of bacteria: (i) initial source of the bacteria, (ii) bacterial species, and (iii) feeding behavior of the host. Results demonstrate that worms increased the spread of bacteria through shedding and transmission. Furthermore, we found that despite P. aeruginosa's relatively high transmission rate among worms, its pathogenic effects reduced the overall number of worms colonized. This study opens new avenues to understand the role of nematodes in the epidemiology and evolution of pathogenic bacteria in the environment.

Loss-of-function of β-catenin bar-1 slows development and activates the Wnt pathway in Caenorhabditis elegans

C. elegans is extensively used to study the Wnt-pathway and most of the core-signalling components are known. Four β-catenins are important gene expression regulators in Wnt-signalling. One of these, bar-1, is part of the canonical Wnt-pathway. Together with Wnt effector pop-1, bar-1 forms a transcription activation complex which regulates the transcription of downstream genes. The effects of bar-1 loss-of-function mutations on many phenotypes have been studied well. However, the effects on global gene expression are unknown. Here we report the effects of a loss-of-function mutation bar-1(ga80). By analysing the transcriptome and developmental phenotyping we show that bar-1(ga80) impairs developmental timing. This developmental difference confounds the comparison of the gene expression profile between the mutant and the reference strain. When corrected for this difference it was possible to identify genes that were directly affected by the bar-1 mutation. We show that the Wnt-pathway itself is activated, as well as transcription factors elt-3, pqm-1, mdl-1 and pha-4 and their associated genes. The outcomes imply that this response compensates for the loss of functional bar-1. Altogether we show that bar-1 loss-of function leads to delayed development possibly caused by an induction of a stress response, reflected by daf-16 activated genes.

Genome-wide variations in a natural isolate of the nematode Caenorhabditis elegans

Background

Increasing genetic and phenotypic differences found among natural isolates of C. elegans have encouraged researchers to explore the natural variation of this nematode species.

Results

Here we report on the identification of genomic differences between the reference strain N2 and the Hawaiian strain CB4856, one of the most genetically distant strains from N2. To identify both small- and large-scale genomic variations (GVs), we have sequenced the CB4856 genome using both Roche 454 (~400 bps single reads) and Illumina GA DNA sequencing methods (101 bps paired-end reads). Compared to previously described variants (available in WormBase), our effort uncovered twice as many single nucleotide variants (SNVs) and increased the number of small InDels almost 20-fold. Moreover, we identified and validated large insertions, most of which range from 150 bps to 1.2 kb in length in the CB4856 strain. Identified GVs had a widespread impact on protein-coding sequences, including 585 single-copy genes that have associated severe phenotypes of reduced viability in RNAi and genetics studies. Sixty of these genes are homologs of human genes associated with diseases. Furthermore, our work confirms previously identified GVs associated with differences in behavioural and biological traits between the N2 and CB4856 strains.

Conclusions

The identified GVs provide a rich resource for future studies that aim to explain the genetic basis for other trait differences between the N2 and CB4856 strains.

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.

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.

Multigenic natural variation underlies Caenorhabditis elegans olfactory preference for the bacterial pathogen Serratia marcescens

The nematode Caenorhabditis elegans can use olfaction to discriminate among different kinds of bacteria, its major food source. We asked how natural genetic variation contributes to choice behavior, focusing on differences in olfactory preference behavior between two wild-type C. elegans strains. The laboratory strain N2 strongly prefers the odor of Serratia marcescens, a soil bacterium that is pathogenic to C. elegans, to the odor of Escherichia coli, a commonly used laboratory food source. The divergent Hawaiian strain CB4856 has a weaker attraction to Serratia than the N2 strain, and this behavioral difference has a complex genetic basis. At least three quantitative trait loci (QTLs) from the CB4856 Hawaii strain (HW) with large effect sizes lead to reduced Serratia preference when introgressed into an N2 genetic background. These loci interact and have epistatic interactions with at least two antagonistic QTLs from HW that increase Serratia preference. The complex genetic architecture of this C. elegans trait is reminiscent of the architecture of mammalian metabolic and behavioral traits.

The prevalence of Caenorhabditis elegans across 1.5 years in selected North German locations: the importance of substrate type, abiotic parameters, and Caenorhabditis competitors

BACKGROUND

Although the nematode Caenorhabditis elegans is a major model organism in diverse biological areas and well studied under laboratory conditions, little is known about its ecology. Therefore, characterization of the species' natural habitats should provide a new perspective on its otherwise well-studied biology. The currently best characterized populations are in France, demonstrating that C. elegans prefers nutrient- and microorganism-rich substrates such as rotting fruits and decomposing plant matter. In order to extend these findings, we sampled C. elegans continuously across 1.5 years from rotting apples and compost heaps in three North German locations.

RESULTS

C. elegans was found throughout summer and autumn in both years. It shares its habitat with the related nematode species C. remanei, which could thus represent an important competitor for a similar ecological niche. The two species were isolated from the same site, but rarely the same substrate sample. In fact, C. elegans was mainly found on compost and C. remanei on rotten apples, possibly suggesting niche separation. The occurrence of C. elegans itself was related to environmental humidity and rain, although the correlation was significant for only one sampling site each. Additional associations between nematode prevalence and abiotic parameters could not be established.

CONCLUSIONS

Taken together, our findings vary from the previous results for French C. elegans populations in that the considered German populations always coexisted with the congeneric species C. remanei (rather than C. briggsae as in France) and that C. elegans prevalence can associate with humidity and rain (rather than temperature, as suggested for French populations). Consideration of additional locations and time points is thus essential for full appreciation of the nematode's natural ecology.

A rapid and massive gene expression shift marking adolescent transition in C. elegans

Organismal development is the most dynamic period of the life cycle, yet we have only a rough understanding of the dynamics of gene expression during adolescent transition. Here we show that adolescence in Caenorhabditis elegans is characterized by a spectacular expression shift of conserved and highly polymorphic genes. Using a high resolution time series we found that in adolescent worms over 10,000 genes changed their expression. These genes were clustered according to their expression patterns. One cluster involved in chromatin remodelling showed a brief up-regulation around 50 h post-hatch. At the same time a spectacular shift in expression was observed. Sequence comparisons for this cluster across many genotypes revealed diversifying selection. Strongly up-regulated genes showed signs of purifying selection in non-coding regions, indicating that adolescence-active genes are constrained on their regulatory properties. Our findings improve our understanding of adolescent transition and help to eliminate experimental artefacts due to incorrect developmental timing.

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.