The rhythmic transcriptional landscape in Caenorhabditis elegans: daily, circadian and novel 16-hour cycling gene expression revealed by RNA-sequencing

The nematode Caenorhabditis elegans is an unconventional model in chronobiology, reported to exhibit physiological and behavioural circadian rhythms while lacking a defined transcriptional oscillator. The extent and importance of circadian rhythms in C. elegans remain uncertain, given a probable lack of functional conservation of key circadian proteins, relatively non-robust reported rhythms and an ambiguous diel ecology. Here, we investigated the temporal coordination of gene expression in C. elegans post-development using RNA sequencing. Over a circadian time series, in which we synchronised nematodes to combined light and temperature cycles, we found clear evidence of daily oscillations in 343 genes using JTK_Cycle. However, rhythms were not well-sustained in constant conditions, with only 13 genes remaining significantly rhythmic. Reanalysis of previous transcriptomic data echoed this finding in identifying far fewer rhythmic genes in constant conditions, while also identifying a greater number of rhythmic genes overall. Weighted gene co-expression network analysis (WGCNA), a hierarchical clustering approach, further confirmed prevalent environmentally driven daily oscillations in the RNA-seq data. This analysis additionally revealed a novel co-expression trend in which over 1000 genes exhibited hitherto unreported 16-hour oscillations, highlighting a new facet of temporal gene expression coordination in C. elegans.