L.B. Snoek

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.

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

BackgroundAnalyzing 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.ResultsWe 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.ConclusionsWe 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.

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.

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.

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.

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.

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.

Maintenance of muscle myosin levels in adult C. elegans requires both the double bromodomain protein BET-1 and sumoylation.

Summary Attenuation of RAS-mediated signalling is a conserved process essential to control cell proliferation, differentiation, and apoptosis. Cooperative interactions between histone modifications such as acetylation, methylation and sumoylation are crucial for proper attenuation in C. elegans, implying that the proteins recognising these histone modifications could also play an important role in attenuation of RAS-mediated signalling. We sought to systematically identify these proteins and found BET-1. BET-1 is a conserved double bromodomain protein that recognises acetyl-lysines on histone tails and maintains the stable fate of various lineages. Unexpectedly, adults lacking both BET-1 and SUMO-1 are depleted of muscle myosin, an essential component of myofibrils. We also show that this muscle myosin depletion does not occur in all animals at a specific time, but rather that the penetrance of the phenotype increases with age. To gain mechanistic insights into this process, we sought to delay the occurrence of the muscle myosin depletion phenotype and found that it requires caspase activity and MEK-dependent signalling. We also performed transcription profiling on these mutants and found an up-regulation of the FGF receptor, egl-15, a tyrosine kinase receptor acting upstream of MEK. Consistent with a MEK requirement, we could delay the muscle phenotype by systemic or hypodermal knock down of egl-15. Thus, this work uncovered a caspase- and MEK-dependent mechanism that acts specifically on ageing adults to maintain the appropriate net level of muscle myosin.

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

BackgroundConnectivity 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.MethodseQTL 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.ResultsWe 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).ConclusionsThis 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.