Search

The heat shock response is the organized molecular response to stressors which disrupt proteostasis, potentially leading to protein misfolding and aggregation. While the regulation of the heat shock response is well-studied in single cells, its coordination at the cell, tissue, and systemic levels of a multicellular organism is poorly understood. To probe the interplay between systemic and cell-autonomous responses, we studied the upregulation of HSP-16.2, a molecular chaperone induced throughout the intestine ofCaenorhabditis elegansfollowing a heat shock, by taking longitudinal measurements in a microfluidic environment. Based on the dynamics of HSP-16.2 accumulation, we showed that a combination of heat shock temperature and duration define the intensity of stress inflicted on the worm and identified two regimes of low and high intensity stress. Modeling the underlying regulatory dynamics implicated the saturation of heat shock protein mRNA production in defining these two regimes and emphasized the importance of time separation between transcription and translation in establishing these dynamics. By applying a heat shock and measuring the response in separate parts of the animals, we implicated thermosensory neurons in accelerating the response and transducing information within the animal. We discuss possible implications of the systemic and cell level aspects and how they coordinate to facilitate the organismal response.

Exposure to high temperatures has an adverse effect on cellular processes and results in activation of the cellular heat shock response (HSR), a highly conserved program of inducible genes to maintain protein homeostasis1. The insulin/IGF-1 signaling (IIS) pathway, which has diverse roles from metabolism to stress response and longevity, is activated as part of the HSR2–4. Recent evidence suggest that the IIS pathway is able to affect proteostasis non-autonomously5,6, yet it is not known if it is activated autonomously in stressed cells or systemically as part of an organismic program. In Caenorhabditis elegans, the single forkhead box O (FOXO) homologue DAF-16 functions as the major target of the IIS pathway7 and, together with the heat-shock factor HSF-1, induce the expression of small heat shock proteins in response to heat shock8–10,3. Here we use a novel microfluidic device that allows precise control of the spatiotemporal temperature profile to show that cellular activation of DAF-16 integrates local temperature sensation with systemic signals. We demonstrate that DAF-16 activation in head sensory neurons is essential for DAF-16 activation in other tissues, but show that no known thermosensory neuron is individually required. Our findings demonstrate that systemic and cell-autonomous aspects of stress response act together to facilitate a coordinated cellular response at the organismic level.

Small non-coding RNAs can exert significant regulatory activity on gene expression in bacteria. In recent years, substantial progress has been made in understanding bacterial gene expression by sRNAs. However, recent findings that demonstrate that families of mRNAs show non-trivial sub-cellular distributions raise the question of how localization may affect the regulatory activity of sRNAs. Here we address this question within a simple mathematical model. We show that the non-uniform spatial distributions of mRNA can alter the threshold-linear response that characterizes sRNAs that act stoichiometrically, and modulate the hierarchy among targets co-regulated by the same sRNA. We also identify conditions where the sub-cellular organization of cofactors in the sRNA pathway can induce spatial heterogeneity on sRNA targets. Our results suggest that under certain conditions, interpretation and modeling of natural and synthetic gene regulatory circuits need to take into account the spatial organization of the transcripts of participating genes.

Two decades into the genomics era the question of mapping sequence to function has evolved from identifying functional elements to characterizing their quantitative properties including, in particular, their specificity and efficiency. Here, we use a large-scale approach to establish a quantitative map between the sequence of a bacterial regulatory RNA and its efficiency in modulating the expression of its targets. Our approach generalizes the sort-seq method, introduced recently to analyze promoter sequences, in order to accurately quantify the efficiency of a large library of sequence variants. We focus on two small RNAs (sRNAs) in E. coli, DsrA and RyhB, and their regulation of both repressed and activated targets. In addition to precisely identifying functional elements in the sRNAs, our data establish quantitative relationships between structural and energetic features of the sRNAs and their regulatory activity, and characterize a large set of direct and indirect interactions between nucleotides. A core of these interactions supports a model where specificity can be enhanced by a rigid molecular structure. Both sRNAs exhibit a modular design with limited cross-interactions, dividing the requirements for structural stability and target binding among modules.

Small animals such as the roundworm C. elegans are excellent models for studying bacterial infection and host response, as well as for genetic and chemical screens. A key methodology is the killing assay, in which the number of surviving animals is tracked as a function of the time post infection. This is a labor-intensive procedure, prone to human error and subjective choices, and often involves undesired perturbation to the animals and their environment. In addition, the survival of animals is just one aspect of a multi-dimensional complex biological process. Here we report a microfluidic-based approach for performing killing assays in worms, compatible with standard assays performed on solid media. In addition to providing accurate and reproducible survival curves at a considerably reduced labor, this approach allows acquisition of a multitude of quantitative data with minimal undesired perturbations. These measurements are obtained automatically at a worm-by-worm resolution using a custom image processing workflow. The proposed approach is simple, scalable, and extendable, and is significantly more economical than standard manual protocols.

Biological homeostasis invokes modulatory responses aimed at stabilizing internal conditions. Using tunable photo- and mechano-stimulation, we identified two distinct categories of homeostatic responses during the sleep-like state of Caenorhabditis elegans (lethargus). In the presence of weak or no stimuli, extended motion caused a subsequent extension of quiescence. The neuropeptide Y receptor homolog, NPR-1, and an inhibitory neuropeptide known to activate it, FLP-18, were required for this process. In the presence of strong stimuli, the correlations between motion and quiescence were disrupted for several minutes but homeostasis manifested as an overall elevation of the time spent in quiescence. This response to strong stimuli required the function of the DAF-16/FOXO transcription factor in neurons, but not that of NPR-1. Conversely, response to weak stimuli did not require the function of DAF-16/FOXO. These findings suggest that routine homeostatic stabilization of sleep may be distinct from homeostatic compensation following a strong disturbance. DOI: http://dx.doi.org/10.7554/eLife.04380.001, eLife digest The regenerative properties of sleep are required by all animals, with even the simplest animal, the nematode Caenorhabditis elegans, displaying a sleep-like state called lethargus. During development, nematodes must pass through four larval stages en route to adulthood, and the end of each stage is preceded by a period of lethargus lasting 2 to 3 hr. Human sleep is divided into distinct stages that recur in a prescribed order throughout the night. Nematodes, on the other hand, simply experience alternating periods of activity and stillness as they sleep. Nevertheless, in both species, any disruptions to sleep automatically lead to adjustments of the rest of the sleep cycle to compensate for the disturbance and to ensure that the organism gets an adequate amount of sleep overall. To date, it has been assumed that a single mechanism is responsible for adjusting the sleep cycle after any disturbance, regardless of its severity. However, Nagy, Tramm, Sanders et al. now show that this is not the case in C. elegans. Sleeping nematodes that were lightly disturbed by exposing them to light or to vibrations—causing them to briefly increase their activity levels—compensated for the disturbance by lengthening their next inactive period. By contrast, worms that were vigorously agitated by stronger vibrations showed a different response: the alternating pattern of stillness and activity was disrupted for several minutes, followed by an overall increase in the length of time spent in the stillness phase. Experiments using genetically modified worms revealed that these two responses involve distinct molecular pathways. A signaling molecule called neuropeptide Y affects the response to minor sleep disruptions, whereas a transcription factor called DAF-16/FOXO is involved in the corresponding role after major disruptions. Given that neuropeptide Y has already been implicated in sleep regulation in humans and flies, it is not implausible that similar mechanisms may occur in response to disturbances of our own sleep. DOI: http://dx.doi.org/10.7554/eLife.04380.002

Post-transcriptional regulation by miRNAs is a widespread and highly conserved phenomenon in metazoans, with several hundreds to thousands of conserved binding sites for each miRNA, and up to two thirds of all genes under miRNA regulation. At the same time, the effect of miRNA regulation on mRNA and protein levels is usually quite modest and associated phenotypes are often weak or subtle. This has given rise to the notion that the highly interconnected miRNA regulatory network exerts its function less through any individual link and more via collective effects that lead to a functional interdependence of network links. We present a Bayesian framework to quantify conservation of miRNA target sites using vertebrate whole-genome alignments. The increased statistical power of our phylogenetic model allows detection of evolutionary correlation in the conservation patterns of site pairs. Such correlations could result from collective functions in the regulatory network. For instance, co-conservation of target site pairs supports a selective benefit of combinatorial regulation by multiple miRNAs. We find that some miRNA families are under pronounced co-targeting constraints, indicating a high connectivity in the regulatory network, while others appear to function in a more isolated way. By analyzing coordinated targeting of different curated gene sets, we observe distinct evolutionary signatures for protein complexes and signaling pathways that could reflect differences in control strategies. Our method is easily scalable to analyze upcoming larger data sets, and readily adaptable to detect high-level selective constraints between other genomic loci. We thus provide a proof-of-principle method to understand regulatory networks from an evolutionary perspective., Author Summary Sequence conservation patterns can be used to assess the functional importance of a genomic locus, e.g. a binding site of a regulatory factor. If one locus is functionally coupled with another, they do not evolve independently and their conservation patterns are correlated. We used evolutionary correlations to study the vertebrate miRNA regulatory network, which is densely connected and shows rapid evolutionary rewiring. Using a Bayesian framework, we identify topological features of this network that are under weak but significant selection, such as combinatorial regulation of one gene by multiple miRNAs, or the coordinated targeting of protein complexes and signaling pathways. Our method is easily scalable to include more genomes as they are being sequenced and can be used to study correlations between other kinds of loci.

Small regulatory RNAs are central players in the regulation of many cellular processes across all kingdoms of life. Experiments in mouse and human have shown that a typical small RNA may regulate the expression of many different genes, suggesting that small RNAs act as global regulators. It is noted though that most targets respond only weakly to the presence of the small RNA. At the same time, evidence in bacteria and animals suggest that the phenotypes associated with small RNA mutants are only due to a few of their targets. Here we assume that targets regulated by a small RNA to control function is in fact small, and propose that the role of the many other weak targets is to confer robustness to the regulation of these few principal targets. Through mathematical modeling we show that auxiliary targets may significantly buffer both number and kinetic fluctuations of the principal targets, with only minor slowdown in the kinetics of response. Analysis of genomic data suggests that auxiliary targets experience a nonspecific evolutionary pressure, playing a role at the system level. Our work is of importance for studies on small RNA functions, and impacts on the understanding of small RNA evolution.

High-density DNA arrays, used to monitor gene expression at a genomic scale, have produced vast amounts of information which require the development of efficient computational methods to analyze them. The important first step is to extract the fundamental patterns of gene expression inherent in the data. This paper describes the application of a novel clustering algorithm, Super-Paramagnetic Clustering (SPC) to analysis of gene expression profiles that were generated recently during a study of the yeast cell cycle. SPC was used to organize genes into biologically relevant clusters that are suggestive for their co-regulation. Some of the advantages of SPC are its robustness against noise and initialization, a clear signature of cluster formation and splitting, and an unsupervised self-organized determination of the number of clusters at each resolution. Our analysis revealed interesting correlated behavior of several groups of genes which has not been previously identified.

Background Sort-seq is an effective approach for simultaneous activity measurements in a large-scale library, combining flow cytometry, deep sequencing, and statistical inference. Such assays enable the characterization of functional landscapes at unprecedented scale for a wide-reaching array of biological molecules and functionalities in vivo. Applications of sort-seq range from footprinting to establishing quantitative models of biological systems and rational design of synthetic genetic elements. Nearly as diverse are implementations of this technique, reflecting key design choices with extensive impact on the scope and accuracy the results. Yet how to make these choices remains unclear. Here we investigate the effects of alternative sort-seq designs and inference methods on the information output using mathematical formulation and simulations. Results We identify key intrinsic properties of any system of interest with practical implications for sort-seq assays, depending on the experimental goals. The fluorescence range and cell-to-cell variability specify the number of sorted populations needed for quantitative measurements that are precise and unbiased. These factors also indicate cases where an enrichment-based approach that uses a single sorted population can offer satisfactory results. These predications of our model are corroborated using re-analysis of published data. We explore implications of these results for quantitative modeling and library design. Conclusions Sort-seq assays can be streamlined by reducing the number of sorted populations, saving considerable resources. Simple preliminary experiments can guide optimal experiment design, minimizing cost while maintaining the maximal information output and avoiding latent biases. These insights can facilitate future applications of this highly adaptable technique. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-2533-5) contains supplementary material, which is available to authorized users.