Temperature-driven coordination of circadian transcriptional regulation.

The circadian clock is an evolutionarily-conserved molecular oscillator that enables species to anticipate rhythmic changes in their environment. At a molecular level, the core clock genes induce circadian oscillations in thousands of genes in a tissue-specific manner, orchestrating myriad biological processes. While previous studies have investigated how the core clock circuit responds to environmental perturbations such as temperature, the downstream effects of such perturbations on circadian regulation remain poorly understood. By analyzing bulk-RNA sequencing of Drosophila fat bodies harvested from flies subjected to different environmental conditions, we demonstrate a highly condition-specific circadian transcriptome: genes are cycling in a temperature-specific manner, and the distributions of their phases also differ between the two conditions. Further employing a reference-based gene regulatory network (Reactome), we find evidence of increased gene-gene coordination at low temperatures and synchronization of rhythmic genes that are network neighbors. We report that the phase differences between cycling genes increase as a function of geodesic distance in the low temperature condition, suggesting increased coordination of cycling on the gene regulatory network. Our results suggest a potential mechanism whereby the circadian clock mediates the fly's response to seasonal changes in temperature.

Temperature-driven coordination of circadian transcriptome regulation.

The circadian rhythm is an evolutionarily-conserved molecular oscillator that enables species to anticipate rhythmic changes in their environment. At a molecular level, the core clock genes induce a circadian oscillation in thousands of genes in a tissue-specific manner, orchestrating myriad biological processes. While studies have investigated how the core clock circuit responds to environmental perturbations such as temperature, the downstream effects of such perturbations on circadian regulation remain poorly understood. By analyzing bulk-RNA sequencing of Drosophila fat bodies harvested from flies subjected to different environmental conditions, we demonstrate a highly condition-specific circadian transcriptome. Further employing a reference-based gene regulatory network (Reactome), we find evidence of increased gene-gene coordination at low temperatures and synchronization of rhythmic genes that are network neighbors. Our results point to the mechanisms by which the circadian clock mediates the fly's response to seasonal changes in temperature.

Rapid threat assessment in the Drosophila thermosensory system.

Neurons that participate in sensory processing often display "ON" responses, i.e., fire transiently at the onset of a stimulus. ON transients are widespread, perhaps universal to sensory coding, yet their function is not always well-understood. Here, we show that ON responses in the Drosophila thermosensory system extrapolate the trajectory of temperature change, priming escape behavior if unsafe thermal conditions are imminent. First, we show that second-order thermosensory projection neurons (TPN-IIIs) and their Lateral Horn targets (TLHONs), display ON responses to thermal stimuli, independent of direction of change (heating or cooling) and of absolute temperature. Instead, they track the rate of temperature change, with TLHONs firing exclusively to rapid changes (>0.2 °C/s). Next, we use connectomics to track TLHONs' output to descending neurons that control walking and escape, and modeling and genetic silencing to demonstrate how ON transients can flexibly amplify aversive responses to small thermal change. Our results suggest that, across sensory systems, ON transients may represent a general mechanism to systematically anticipate and respond to salient or dangerous conditions.

Circadian programming of the ellipsoid body sleep homeostat in Drosophila.

Homeostatic and circadian processes collaborate to appropriately time and consolidate sleep and wake. To understand how these processes are integrated, we scheduled brief sleep deprivation at different times of day in Drosophila and find elevated morning rebound compared to evening. These effects depend on discrete morning and evening clock neurons, independent of their roles in circadian locomotor activity. In the R5 ellipsoid body sleep homeostat, we identified elevated morning expression of activity dependent and presynaptic gene expression as well as the presynaptic protein BRUCHPILOT consistent with regulation by clock circuits. These neurons also display elevated calcium levels in response to sleep loss in the morning, but not the evening consistent with the observed time-dependent sleep rebound. These studies reveal the circuit and molecular mechanisms by which discrete circadian clock neurons program a homeostatic sleep center.

EMBEDR: Distinguishing signal from noise in single-cell omics data.

Single-cell "omics"-based measurements are often high dimensional so that dimensionality reduction (DR) algorithms are necessary for data visualization and analysis. The lack of methods for separating signal from noise in DR outputs has limited their utility in generating data-driven discoveries in single-cell data. In this work we present EMBEDR, which assesses the output of any DR algorithm to distinguish evidence of structure from algorithm-induced noise in DR outputs. We apply EMBEDR to DR-generated representations of single-cell omics data of several modalities to show where they visually show real-not spurious-structure. EMBEDR generates a "p" value for each sample, allowing for direct comparisons of DR algorithms and facilitating optimization of algorithm hyperparameters. We show that the scale of a sample's neighborhood can thus be determined and used to generate a novel "cell-wise optimal" embedding. EMBEDR is available as a Python package for immediate use.

The E3 ubiquitin ligase adaptor Tango10 links the core circadian clock to neuropeptide and behavioral rhythms.

Circadian transcriptional timekeepers in pacemaker neurons drive profound daily rhythms in sleep and wake. Here we reveal a molecular pathway that links core transcriptional oscillators to neuronal and behavioral rhythms. Using two independent genetic screens, we identified mutants of Transport and Golgi organization 10 (Tango10) with poor behavioral rhythmicity. Tango10 expression in pacemaker neurons expressing the neuropeptide PIGMENT-DISPERSING FACTOR (PDF) is required for robust rhythms. Loss of Tango10 results in elevated PDF accumulation in nerve terminals even in mutants lacking a functional core clock. TANGO10 protein itself is rhythmically expressed in PDF terminals. Mass spectrometry of TANGO10 complexes reveals interactions with the E3 ubiquitin ligase CULLIN 3 (CUL3). CUL3 depletion phenocopies Tango10 mutant effects on PDF even in the absence of the core clock gene timeless Patch clamp electrophysiology in Tango10 mutant neurons demonstrates elevated spontaneous firing potentially due to reduced voltage-gated Shaker-like potassium currents. We propose that Tango10/Cul3 transduces molecular oscillations from the core clock to neuropeptide release important for behavioral rhythms.

Robustness and plasticity in Drosophila heat avoidance.

Simple innate behavior is often described as hard-wired and largely inflexible. Here, we show that the avoidance of hot temperature, a simple innate behavior, contains unexpected plasticity in Drosophila. First, we demonstrate that hot receptor neurons of the antenna and their molecular heat sensor, Gr28B.d, are essential for flies to produce escape turns away from heat. High-resolution fly tracking combined with a 3D simulation of the thermal environment shows that, in steep thermal gradients, the direction of escape turns is determined by minute temperature differences between the antennae (0.1°-1 °C). In parallel, live calcium imaging confirms that such small stimuli reliably activate both peripheral thermosensory neurons and central circuits. Next, based on our measurements, we evolve a fly/vehicle model with two symmetrical sensors and motors (a "Braitenberg vehicle") which closely approximates basic fly thermotaxis. Critical differences between real flies and the hard-wired vehicle reveal that fly heat avoidance involves decision-making, relies on rapid learning, and is robust to new conditions, features generally associated with more complex behavior.

TimeTrial: An Interactive Application for Optimizing the Design and Analysis of Transcriptomic Time-Series Data in Circadian Biology Research.

The circadian rhythm drives the oscillatory expression of thousands of genes across all tissues, coordinating physiological processes. The effect of this rhythm on health has generated increasing interest in discovering genes under circadian control by searching for periodic patterns in transcriptomic time-series experiments. While algorithms for detecting cycling transcripts have advanced, there remains little guidance quantifying the effect of experimental design and analysis choices on cycling detection accuracy. We present TimeTrial, a user-friendly benchmarking framework using both real and synthetic data to investigate cycle detection algorithms' performance and improve circadian experimental design. Results show that the optimal choice of analysis method depends on the sampling scheme, noise level, and shape of the waveform of interest and provides guidance on the impact of sampling frequency and duration on cycling detection accuracy. The TimeTrial software is freely available for download and may also be accessed through a web interface. By supplying a tool to vary and optimize experimental design considerations, TimeTrial will enhance circadian transcriptomics studies.

Universal method for robust detection of circadian state from gene expression.

Circadian clocks play a key role in regulating a vast array of biological processes, with significant implications for human health. Accurate assessment of physiological time using transcriptional biomarkers found in human blood can significantly improve diagnosis of circadian disorders and optimize the delivery time of therapeutic treatments. To be useful, such a test must be accurate, minimally burdensome to the patient, and readily generalizable to new data. A major obstacle in development of gene expression biomarker tests is the diversity of measurement platforms and the inherent variability of the data, often resulting in predictors that perform well in the original datasets but cannot be universally applied to new samples collected in other settings. Here, we introduce TimeSignature, an algorithm that robustly infers circadian time from gene expression. We demonstrate its application in data from three independent studies using distinct microarrays and further validate it against a new set of samples profiled by RNA-sequencing. Our results show that TimeSignature is more accurate and efficient than competing methods, estimating circadian time to within 2 h for the majority of samples. Importantly, we demonstrate that once trained on data from a single study, the resulting predictor can be universally applied to yield highly accurate results in new data from other studies independent of differences in study population, patient protocol, or assay platform without renormalizing the data or retraining. This feature is unique among expression-based predictors and addresses a major challenge in the development of generalizable, clinically useful tests.

Control of Stochastic and Induced Switching in Biophysical Networks.

Noise caused by fluctuations at the molecular level is a fundamental part of intracellular processes. While the response of biological systems to noise has been studied extensively, there has been limited understanding of how to exploit it to induce a desired cell state. Here we present a scalable, quantitative method based on the Freidlin-Wentzell action to predict and control noise-induced switching between different states in genetic networks that, conveniently, can also control transitions between stable states in the absence of noise. We apply this methodology to models of cell differentiation and show how predicted manipulations of tunable factors can induce lineage changes, and further utilize it to identify new candidate strategies for cancer therapy in a cell death pathway model. This framework offers a systems approach to identifying the key factors for rationally manipulating biophysical dynamics, and should also find use in controlling other classes of noisy complex networks.

Spatial and functional heterogeneities shape collective behavior of tumor-immune networks.

Tumor growth involves a dynamic interplay between cancer cells and host cells, which collectively form a tumor microenvironmental network that either suppresses or promotes tumor growth under different conditions. The transition from tumor suppression to tumor promotion is mediated by a tumor-induced shift in the local immune state, and despite the clinical challenge this shift poses, little is known about how such dysfunctional immune states are initiated. Clinical and experimental observations have indicated that differences in both the composition and spatial distribution of different cell types and/or signaling molecules within the tumor microenvironment can strongly impact tumor pathogenesis and ultimately patient prognosis. How such "functional" and "spatial" heterogeneities confer such effects, however, is not known. To investigate these phenomena at a level currently inaccessible by direct observation, we developed a computational model of a nascent metastatic tumor capturing salient features of known tumor-immune interactions that faithfully recapitulates key features of existing experimental observations. Surprisingly, over a wide range of model formulations, we observed that heterogeneity in both spatial organization and cell phenotype drove the emergence of immunosuppressive network states. We determined that this observation is general and robust to parameter choice by developing a systems-level sensitivity analysis technique, and we extended this analysis to generate other parameter-independent, experimentally testable hypotheses. Lastly, we leveraged this model as an in silico test bed to evaluate potential strategies for engineering cell-based therapies to overcome tumor associated immune dysfunction and thereby identified modes of immune modulation predicted to be most effective. Collectively, this work establishes a new integrated framework for investigating and modulating tumor-immune networks and provides insights into how such interactions may shape early stages of tumor formation.

Histone methyltransferase MMSET/NSD2 alters EZH2 binding and reprograms the myeloma epigenome through global and focal changes in H3K36 and H3K27 methylation.

Overexpression of the histone methyltransferase MMSET in t(4;14)+ multiple myeloma patients is believed to be the driving factor in the pathogenesis of this subtype of myeloma. MMSET catalyzes dimethylation of lysine 36 on histone H3 (H3K36me2), and its overexpression causes a global increase in H3K36me2, redistributing this mark in a broad, elevated level across the genome. Here, we demonstrate that an increased level of MMSET also induces a global reduction of lysine 27 trimethylation on histone H3 (H3K27me3). Despite the net decrease in H3K27 methylation, specific genomic loci exhibit enhanced recruitment of the EZH2 histone methyltransferase and become hypermethylated on this residue. These effects likely contribute to the myeloma phenotype since MMSET-overexpressing cells displayed increased sensitivity to EZH2 inhibition. Furthermore, we demonstrate that such MMSET-mediated epigenetic changes require a number of functional domains within the protein, including PHD domains that mediate MMSET recruitment to chromatin. In vivo, targeting of MMSET by an inducible shRNA reversed histone methylation changes and led to regression of established tumors in athymic mice. Together, our work elucidates previously unrecognized interplay between MMSET and EZH2 in myeloma oncogenesis and identifies domains to be considered when designing inhibitors of MMSET function.

Intrinsic bursting of AII amacrine cells underlies oscillations in the rd1 mouse retina.

In many forms of retinal degeneration, photoreceptors die but inner retinal circuits remain intact. In the rd1 mouse, an established model for blinding retinal diseases, spontaneous activity in the coupled network of AII amacrine and ON cone bipolar cells leads to rhythmic bursting of ganglion cells. Since such activity could impair retinal and/or cortical responses to restored photoreceptor function, understanding its nature is important for developing treatments of retinal pathologies. Here we analyzed a compartmental model of the wild-type mouse AII amacrine cell to predict that the cell's intrinsic membrane properties, specifically, interacting fast Na and slow, M-type K conductances, would allow its membrane potential to oscillate when light-evoked excitatory synaptic inputs were withdrawn following photoreceptor degeneration. We tested and confirmed this hypothesis experimentally by recording from AIIs in a slice preparation of rd1 retina. Additionally, recordings from ganglion cells in a whole mount preparation of rd1 retina demonstrated that activity in AIIs was propagated unchanged to elicit bursts of action potentials in ganglion cells. We conclude that oscillations are not an emergent property of a degenerated retinal network. Rather, they arise largely from the intrinsic properties of a single retinal interneuron, the AII amacrine cell.

Adaptation to background light enables contrast coding at rod bipolar cell synapses.

Rod photoreceptors contribute to vision over an ∼ 6-log-unit range of light intensities. The wide dynamic range of rod vision is thought to depend upon light intensity-dependent switching between two parallel pathways linking rods to ganglion cells: a rod → rod bipolar (RB) cell pathway that operates at dim backgrounds and a rod → cone → cone bipolar cell pathway that operates at brighter backgrounds. We evaluated this conventional model of rod vision by recording rod-mediated light responses from ganglion and AII amacrine cells and by recording RB-mediated synaptic currents from AII amacrine cells in mouse retina. Contrary to the conventional model, we found that the RB pathway functioned at backgrounds sufficient to activate the rod → cone pathway. As background light intensity increased, the RB's role changed from encoding the absorption of single photons to encoding contrast modulations around mean luminance. This transition is explained by the intrinsic dynamics of transmission from RB synapses.

Balanced synaptic impact via distance-dependent synapse distribution and complementary expression of AMPARs and NMDARs in hippocampal dendrites.

Neuronal computation involves the integration of synaptic inputs that are often distributed over expansive dendritic trees, suggesting the need for compensatory mechanisms that enable spatially disparate synapses to influence neuronal output. In hippocampal CA1 pyramidal neurons, such mechanisms have indeed been reported, which normalize either the ability of distributed synapses to drive action potential initiation in the axon or their ability to drive dendritic spiking locally. Here we report that these mechanisms can coexist, through an elegant combination of distance-dependent regulation of synapse number and synaptic expression of AMPA and NMDA receptors. Together, these complementary gradients allow individual dendrites in both the apical and basal dendritic trees of hippocampal neurons to operate as facile computational subunits capable of supporting both global integration in the soma/axon and local integration in the dendrite.

Realistic control of network dynamics.

The control of complex networks is of paramount importance in areas as diverse as ecosystem management, emergency response and cell reprogramming. A fundamental property of networks is that perturbations to one node can affect other nodes, potentially causing the entire system to change behaviour or fail. Here we show that it is possible to exploit the same principle to control network behaviour. Our approach accounts for the nonlinear dynamics inherent to real systems, and allows bringing the system to a desired target state even when this state is not directly accessible due to constraints that limit the allowed interventions. Applications show that this framework permits reprogramming a network to a desired task, as well as rescuing networks from the brink of failure-which we illustrate through the mitigation of cascading failures in a power-grid network and the identification of potential drug targets in a signalling network of human cancer.

Hippocampal pyramidal neurons comprise two distinct cell types that are countermodulated by metabotropic receptors.

Relating the function of neuronal cell types to information processing and behavior is a central goal of neuroscience. In the hippocampus, pyramidal cells in CA1 and the subiculum process sensory and motor cues to form a cognitive map encoding spatial, contextual, and emotional information, which they transmit throughout the brain. Do these cells constitute a single class or are there multiple cell types with specialized functions? Using unbiased cluster analysis, we show that there are two morphologically and electrophysiologically distinct principal cell types that carry hippocampal output. We show further that these two cell types are inversely modulated by the synergistic action of glutamate and acetylcholine acting on metabotropic receptors that are central to hippocampal function. Combined with prior connectivity studies, our results support a model of hippocampal processing in which the two pyramidal cell types are predominantly segregated into two parallel pathways that process distinct modalities of information.

Synaptic amplification by dendritic spines enhances input cooperativity.

Dendritic spines are the nearly ubiquitous site of excitatory synaptic input onto neurons and as such are critically positioned to influence diverse aspects of neuronal signalling. Decades of theoretical studies have proposed that spines may function as highly effective and modifiable chemical and electrical compartments that regulate synaptic efficacy, integration and plasticity. Experimental studies have confirmed activity-dependent structural dynamics and biochemical compartmentalization by spines. However, there is a longstanding debate over the influence of spines on the electrical aspects of synaptic transmission and dendritic operation. Here we measure the amplitude ratio of spine head to parent dendrite voltage across a range of dendritic compartments and calculate the associated spine neck resistance (R(neck)) for spines at apical trunk dendrites in rat hippocampal CA1 pyramidal neurons. We find that R(neck) is large enough (~500 MΩ) to amplify substantially the spine head depolarization associated with a unitary synaptic input by ~1.5- to ~45-fold, depending on parent dendritic impedance. A morphologically realistic compartmental model capable of reproducing the observed spatial profile of the amplitude ratio indicates that spines provide a consistently high-impedance input structure throughout the dendritic arborization. Finally, we demonstrate that the amplification produced by spines encourages electrical interaction among coactive inputs through an R(neck)-dependent increase in spine head voltage-gated conductance activation. We conclude that the electrical properties of spines promote nonlinear dendritic processing and associated forms of plasticity and storage, thus fundamentally enhancing the computational capabilities of neurons.

The mechanisms of repetitive spike generation in an axonless retinal interneuron.

Several types of retinal interneurons exhibit spikes but lack axons. One such neuron is the AII amacrine cell, in which spikes recorded at the soma exhibit small amplitudes (<10 mV) and broad time courses (>5 ms). Here, we used electrophysiological recordings and computational analysis to examine the mechanisms underlying this atypical spiking. We found that somatic spikes likely represent large, brief action potential-like events initiated in a single, electrotonically distal dendritic compartment. In this same compartment, spiking undergoes slow modulation, likely by an M-type K conductance. The structural correlate of this compartment is a thin neurite that extends from the primary dendritic tree: local application of TTX to this neurite, or excision of it, eliminates spiking. Thus, the physiology of the axonless AII is much more complex than would be anticipated from morphological descriptions and somatic recordings; in particular, the AII possesses a single dendritic structure that controls its firing pattern.