Lymphocyte Circadian Clocks Control Lymph Node Trafficking and Adaptive Immune Responses.

Lymphocytes circulate through lymph nodes (LN) in search for antigen in what is believed to be a continuous process. Here, we show that lymphocyte migration through lymph nodes and lymph occurred in a non-continuous, circadian manner. Lymphocyte homing to lymph nodes peaked at night onset, with cells leaving the tissue during the day. This resulted in strong oscillations in lymphocyte cellularity in lymph nodes and efferent lymphatic fluid. Using lineage-specific genetic ablation of circadian clock function, we demonstrated this to be dependent on rhythmic expression of promigratory factors on lymphocytes. Dendritic cell numbers peaked in phase with lymphocytes, with diurnal oscillations being present in disease severity after immunization to induce experimental autoimmune encephalomyelitis (EAE). These rhythms were abolished by genetic disruption of T cell clocks, demonstrating a circadian regulation of lymphocyte migration through lymph nodes with time-of-day of immunization being critical for adaptive immune responses weeks later.

Alternative polyadenylation factor CPSF6 regulates temperature compensation of the mammalian circadian clock.

A defining property of circadian clocks is temperature compensation, characterized by the resilience of their near 24-hour free-running periods against changes in environmental temperature within the physiological range. While temperature compensation is evolutionary conserved across different taxa of life and has been studied within many model organisms, its molecular underpinnings remain elusive. Posttranscriptional regulations such as temperature-sensitive alternative splicing or phosphorylation have been described as underlying reactions. Here, we show that knockdown of cleavage and polyadenylation specificity factor subunit 6 (CPSF6), a key regulator of 3'-end cleavage and polyadenylation, significantly alters circadian temperature compensation in human U-2 OS cells. We apply a combination of 3'-end-RNA-seq and mass spectrometry-based proteomics to globally quantify changes in 3' UTR length as well as gene and protein expression between wild-type and CPSF6 knockdown cells and their dependency on temperature. Since changes in temperature compensation behavior should be reflected in alterations of temperature responses within one or all of the 3 regulatory layers, we statistically assess differential responses upon changes in ambient temperature between wild-type and CPSF6 knockdown cells. By this means, we reveal candidate genes underlying circadian temperature compensation, including eukaryotic translation initiation factor 2 subunit 1 (EIF2S1).

Network switches and their role in circadian clocks.

Circadian rhythms are generated by complex interactions among genes and proteins. Self-sustained ∼24 h oscillations require negative feedback loops and sufficiently strong nonlinearities that are the product of molecular and network switches. Here, we review common mechanisms to obtain switch-like behavior, including cooperativity, antagonistic enzymes, multisite phosphorylation, positive feedback, and sequestration. We discuss how network switches play a crucial role as essential components in cellular circadian clocks, serving as integral parts of transcription-translation feedback loops that form the basis of circadian rhythm generation. The design principles of network switches and circadian clocks are illustrated by representative mathematical models that include bistable systems and negative feedback loops combined with Hill functions. This work underscores the importance of negative feedback loops and network switches as essential design principles for biological oscillations, emphasizing how an understanding of theoretical concepts can provide insights into the mechanisms generating biological rhythms.

Analyses of circRNA Expression throughout the Light-Dark Cycle Reveal a Strong Regulation of Cdr1as, Associated with Light Entrainment in the SCN.

Circular RNAs (circRNAs) are a large class of relatively stable RNA molecules that are highly expressed in animal brains. Many circRNAs have been associated with CNS disorders accompanied by an aberrant wake-sleep cycle. However, the regulation of circRNAs in brain homeostasis over daily light-dark (LD) cycles has not been characterized. Here, we aim to quantify the daily expression changes of circRNAs in physiological conditions in healthy adult animals. Using newly generated and public RNA-Seq data, we monitored circRNA expression throughout the 12:12 h LD cycle in various mouse brain regions. We identified that Cdr1as, a conserved circRNA that regulates synaptic transmission, is highly expressed in the suprachiasmatic nucleus (SCN), the master circadian pacemaker. Despite its high stability, Cdr1as has a very dynamic expression in the SCN throughout the LD cycle, as well as a significant regulation in the hippocampus following the entry into the dark phase. Computational integration of different public datasets predicted that Cdr1as is important for regulating light entrainment in the SCN. We hypothesize that the expression changes of Cdr1as in the SCN, particularly during the dark phase, are associated with light-induced phase shifts. Importantly, our work revises the current beliefs about natural circRNA stability and suggests that the time component must be considered when studying circRNA regulation.

Weak coupling between intracellular feedback loops explains dissociation of clock gene dynamics.

Circadian rhythms are generated by interlocked transcriptional-translational negative feedback loops (TTFLs), the molecular process implemented within a cell. The contributions, weighting and balancing between the multiple feedback loops remain debated. Dissociated, free-running dynamics in the expression of distinct clock genes has been described in recent experimental studies that applied various perturbations such as slice preparations, light pulses, jet-lag, and culture medium exchange. In this paper, we provide evidence that this "presumably transient" dissociation of circadian gene expression oscillations may occur at the single-cell level. Conceptual and detailed mechanistic mathematical modeling suggests that such dissociation is due to a weak interaction between multiple feedback loops present within a single cell. The dissociable loops provide insights into underlying mechanisms and general design principles of the molecular circadian clock.

Coherency of circadian rhythms in the SCN is governed by the interplay of two coupling factors.

Circadian clocks are autonomous oscillators driving daily rhythms in physiology and behavior. In mammals, a network of coupled neurons in the suprachiasmatic nucleus (SCN) is entrained to environmental light-dark cycles and orchestrates the timing of peripheral organs. In each neuron, transcriptional feedbacks generate noisy oscillations. Coupling mediated by neuropeptides such as VIP and AVP lends precision and robustness to circadian rhythms. The detailed coupling mechanisms between SCN neurons are debated. We analyze organotypic SCN slices from neonatal and adult mice in wild-type and multiple knockout conditions. Different degrees of rhythmicity are quantified by pixel-level analysis of bioluminescence data. We use empirical orthogonal functions (EOFs) to characterize spatio-temporal patterns. Simulations of coupled stochastic single cell oscillators can reproduce the diversity of observed patterns. Our combination of data analysis and modeling provides deeper insight into the enormous complexity of the data: (1) Neonatal slices are typically stronger oscillators than adult slices pointing to developmental changes of coupling. (2) Wild-type slices are completely synchronized and exhibit specific spatio-temporal patterns of phases. (3) Some slices of Cry double knockouts obey impaired synchrony that can lead to co-existing rhythms ("splitting"). (4) The loss of VIP-coupling leads to desynchronized rhythms with few residual local clusters. Additional information was extracted from co-culturing slices with rhythmic neonatal wild-type SCNs. These co-culturing experiments were simulated using external forcing terms representing VIP and AVP signaling. The rescue of rhythmicity via co-culturing lead to surprising results, since a cocktail of AVP-antagonists improved synchrony. Our modeling suggests that these counter-intuitive observations are pointing to an antagonistic action of VIP and AVP coupling. Our systematic theoretical and experimental study shows that dual coupling mechanisms can explain the astonishing complexity of spatio-temporal patterns in SCN slices.

A Robust Model for Circadian Redox Oscillations.

The circadian clock is an endogenous oscillator that controls daily rhythms in metabolism, physiology, and behavior. Although the timekeeping components differ among species, a common design principle is a transcription-translation negative feedback loop. However, it is becoming clear that other mechanisms can contribute to the generation of 24 h rhythms. Peroxiredoxins (Prxs) exhibit 24 h rhythms in their redox state in all kingdoms of life. In mammalian adrenal gland, heart and brown adipose tissue, such rhythms are generated as a result of an inactivating hyperoxidation reaction that is reduced by coordinated import of sulfiredoxin (Srx) into the mitochondria. However, a quantitative description of the Prx/Srx oscillating system is still missing. We investigate the basic principles that generate mitochondrial Prx/Srx rhythms using computational modeling. We observe that the previously described delay in mitochondrial Srx import, in combination with an appropriate separation of fast and slow reactions, is sufficient to generate robust self-sustained relaxation-like oscillations. We find that our conceptual model can be regarded as a series of three consecutive phases and two temporal switches, highlighting the importance of delayed negative feedback and switches in the generation of oscillations.

An Inactivation Switch Enables Rhythms in a Neurospora Clock Model.

Autonomous endogenous time-keeping is ubiquitous across many living organisms, known as the circadian clock when it has a period of about 24 h. Interestingly, the fundamental design principle with a network of interconnected negative and positive feedback loops is conserved through evolution, although the molecular components differ. Filamentous fungus Neurospora crassa is a well-established chrono-genetics model organism to investigate the underlying mechanisms. The core negative feedback loop of the clock of Neurospora is composed of the transcription activator White Collar Complex (WCC) (heterodimer of WC1 and WC2) and the inhibitory element called FFC complex, which is made of FRQ (Frequency protein), FRH (Frequency interacting RNA Helicase) and CK1a (Casein kinase 1a). While exploring their temporal dynamics, we investigate how limit cycle oscillations arise and how molecular switches support self-sustained rhythms. We develop a mathematical model of 10 variables with 26 parameters to understand the interactions and feedback among WC1 and FFC elements in nuclear and cytoplasmic compartments. We performed control and bifurcation analysis to show that our novel model produces robust oscillations with a wild-type period of 22.5 h. Our model reveals a switch between WC1-induced transcription and FFC-assisted inactivation of WC1. Using the new model, we also study the possible mechanisms of glucose compensation. A fairly simple model with just three nonlinearities helps to elucidate clock dynamics, revealing a mechanism of rhythms' production. The model can further be utilized to study entrainment and temperature compensation.

Moran's I quantifies spatio-temporal pattern formation in neural imaging data.

MOTIVATION: Neural activities of the brain occur through the formation of spatio-temporal patterns. In recent years, macroscopic neural imaging techniques have produced a large body of data on these patterned activities, yet a numerical measure of spatio-temporal coherence has often been reduced to the global order parameter, which does not uncover the degree of spatial correlation. Here, we propose to use the spatial autocorrelation measure Moran's I, which can be applied to capture dynamic signatures of spatial organization. We demonstrate the application of this technique to collective cellular circadian clock activities measured in the small network of the suprachiasmatic nucleus (SCN) in the hypothalamus. RESULTS: We found that Moran's I is a practical quantitative measure of the degree of spatial coherence in neural imaging data. Initially developed with a geographical context in mind, Moran's I accounts for the spatial organization of any interacting units. Moran's I can be modified in accordance with the characteristic length scale of a neural activity pattern. It allows a quantification of statistical significance levels for the observed patterns. We describe the technique applied to synthetic datasets and various experimental imaging time-series from cultured SCN explants. It is demonstrated that major characteristics of the collective state can be described by Moran's I and the traditional Kuramoto order parameter R in a complementary fashion. AVAILABILITY AND IMPLEMENTATION: Python 2.7 code of illustrative examples can be found in the Supplementary Material. CONTACT: christoph.schmal@charite.de or grigory.bordyugov@hu-berlin.de. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Identification of Novel Nuclear Factor of Activated T Cell (NFAT)-associated Proteins in T Cells.

Transcription factors of the nuclear factor of activated T cell (NFAT) family are essential for antigen-specific T cell activation and differentiation. Their cooperative DNA binding with other transcription factors, such as AP1 proteins (FOS, JUN, and JUNB), FOXP3, IRFs, and EGR1, dictates the gene regulatory action of NFATs. To identify as yet unknown interaction partners of NFAT, we purified biotin-tagged NFATc1/αA, NFATc1/ßC, and NFATc2/C protein complexes and analyzed their components by stable isotope labeling by amino acids in cell culture-based mass spectrometry. We revealed more than 170 NFAT-associated proteins, half of which are involved in transcriptional regulation. Among them are many hitherto unknown interaction partners of NFATc1 and NFATc2 in T cells, such as Raptor, CHEK1, CREB1, RUNX1, SATB1, Ikaros, and Helios. The association of NFATc2 with several other transcription factors is DNA-dependent, indicating cooperative DNA binding. Moreover, our computational analysis discovered that binding motifs for RUNX and CREB1 are found preferentially in the direct vicinity of NFAT-binding motifs and in a distinct orientation to them. Furthermore, we provide evidence that mTOR and CHEK1 kinase activity influence NFAT's transcriptional potency. Finally, our dataset of NFAT-associated proteins provides a good basis to further study NFAT's diverse functions and how these are modulated due to the interplay of multiple interaction partners.

Feedback Loops of the Mammalian Circadian Clock Constitute Repressilator.

Mammals evolved an endogenous timing system to coordinate their physiology and behaviour to the 24h period of the solar day. While it is well accepted that circadian rhythms are generated by intracellular transcriptional feedback loops, it is still debated which network motifs are necessary and sufficient for generating self-sustained oscillations. Here, we systematically explore a data-based circadian oscillator model with multiple negative and positive feedback loops and identify a series of three subsequent inhibitions known as "repressilator" as a core element of the mammalian circadian oscillator. The central role of the repressilator motif is consistent with time-resolved ChIP-seq experiments of circadian clock transcription factors and loss of rhythmicity in core clock gene knockouts.

Ras-mediated deregulation of the circadian clock in cancer.

Circadian rhythms are essential to the temporal regulation of molecular processes in living systems and as such to life itself. Deregulation of these rhythms leads to failures in biological processes and eventually to the manifestation of pathological phenotypes including cancer. To address the questions as to what are the elicitors of a disrupted clock in cancer, we applied a systems biology approach to correlate experimental, bioinformatics and modelling data from several cell line models for colorectal and skin cancer. We found strong and weak circadian oscillators within the same type of cancer and identified a set of genes, which allows the discrimination between the two oscillator-types. Among those genes are IFNGR2, PITX2, RFWD2, PPARγ, LOXL2, Rab6 and SPARC, all involved in cancer-related pathways. Using a bioinformatics approach, we extended the core-clock network and present its interconnection to the discriminative set of genes. Interestingly, such gene signatures link the clock to oncogenic pathways like the RAS/MAPK pathway. To investigate the potential impact of the RAS/MAPK pathway - a major driver of colorectal carcinogenesis - on the circadian clock, we used a computational model which predicted that perturbation of BMAL1-mediated transcription can generate the circadian phenotypes similar to those observed in metastatic cell lines. Using an inducible RAS expression system, we show that overexpression of RAS disrupts the circadian clock and leads to an increase of the circadian period while RAS inhibition causes a shortening of period length, as predicted by our mathematical simulations. Together, our data demonstrate that perturbations induced by a single oncogene are sufficient to deregulate the mammalian circadian clock.

Transcription factor co-occupied regions in the murine genome constitute T-helper-cell subtype-specific enhancers.

Transcription factors (TFs) regulate cell-type-specific gene expression programs by combinatorial binding to cis-genomic elements, particularly enhancers, subsequently leading to the recruitment of cofactors, and the general transcriptional machinery to target genes. Using data integration of genome-wide TF binding profiles, we defined regions with combinatorial binding of lineage-specific master TFs (T-BET, GATA3, and ROR-γt) and STATs (STAT1 and STAT4, STAT6, and STAT3) in murine T helper (Th) 1, Th2, and Th17 cells, respectively. Stringently excluding promoter regions, we revealed precise genomic elements which were preferentially associated with the enhancer marks p300 and H3K4me1. Furthermore, closely adjacent TF co-occupied regions constituted larger enhancer domains in the respective Th-cell subset (177 in Th1, 141 in Th2, and 266 in Th17 cells) with characteristics of so-called super-enhancers. Importantly, 89% of these super-enhancer regions were Th-cell subtype-specific. Genes associated with super-enhancers, including relevant Th-cell genes (such as Ifng in Th1, Il13 in Th2, and Il17a in Th17 cells), showed strong transcriptional activity. Altogether, the discovered catalog of enhancers provides information about crucial Th-cell subtype-specific regulatory hubs, which will be useful for revealing cell-type-specific gene regulation processes.

The structural code of cyanobacterial genomes.

A periodic bias in nucleotide frequency with a period of about 11 bp is characteristic for bacterial genomes. This signal is commonly interpreted to relate to the helical pitch of negatively supercoiled DNA. Functions in supercoiling-dependent RNA transcription or as a 'structural code' for DNA packaging have been suggested. Cyanobacterial genomes showed especially strong periodic signals and, on the other hand, DNA supercoiling and supercoiling-dependent transcription are highly dynamic and underlie circadian rhythms of these phototrophic bacteria. Focusing on this phylum and dinucleotides, we find that a minimal motif of AT-tracts (AT2) yields the strongest signal. Strong genome-wide periodicity is ancestral to a clade of unicellular and polyploid species but lost upon morphological transitions into two baeocyte-forming and a symbiotic species. The signal is intermediate in heterocystous species and weak in monoploid picocyanobacteria. A pronounced 'structural code' may support efficient nucleoid condensation and segregation in polyploid cells. The major source of the AT2 signal are protein-coding regions, where it is encoded preferentially in the first and third codon positions. The signal shows only few relations to supercoiling-dependent and diurnal RNA transcription in Synechocystis sp. PCC 6803. Strong and specific signals in two distinct transposons suggest roles in transposase transcription and transpososome formation.

Coupling Controls the Synchrony of Clock Cells in Development and Knockouts.

In mammals, a network of coupled neurons within the hypothalamus coordinates physiological rhythms with daily changes in the environment. In each neuron, delayed negative transcriptional feedbacks generate oscillations, albeit noisy and unreliable ones. Coupling mediated by diffusible neuropeptides lends precision and robustness to circadian rhythms. The double knockout of Cryptochrome Cry turns adult mice arrhythmic. But, remarkably, double knockout neonates continue to show robust oscillation much like wild-type neonates and appear to lose rhythmicity with development. We study quantitatively dispersed neurons and brain slices from wild-type and Cry double knockout mice to understand the links between single cell rhythmicity and intercellular coupling. We quantify oscillator properties of dispersed cells using nonlinear regression and study bifurcations diagrams of network models. We find that varying just three parameters-oscillator strength, strength of coupling, and timing of coupling-can reproduce experimentally observed features. In particular, modeling reveals that minor changes in timing of coupling can destroy synchronization as observed in adult slices from knockout mice.

Timing of neuropeptide coupling determines synchrony and entrainment in the mammalian circadian clock.

Robust synchronization is a critical feature of several systems including the mammalian circadian clock. The master circadian clock in mammals consists of about 20000 'sloppy' neuronal oscillators within the hypothalamus that keep robust time by synchronization driven by inter-neuronal coupling. The complete understanding of this synchronization in the mammalian circadian clock and the mechanisms underlying it remain an open question. Experiments and computational studies have shown that coupling individual oscillators can achieve robust synchrony, despite heterogeneity and different network topologies. But, much less is known regarding the mechanisms and circuits involved in achieving this coupling, due to both system complexity and experimental limitations. Here, we computationally study the coupling mediated by the primary coupling neuropeptide, vasoactive intestinal peptide (VIP) and its canonical receptor, VPAC2R, using the transcriptional elements and generic mode of VIP-VPAC2R signaling. We find that synchrony is only possible if VIP (an inducer of Per expression) is released in-phase with activators of Per expression. Moreover, anti-phasic VIP release suppresses coherent rhythms by moving the network into a desynchronous state. Importantly, experimentally observed rhythms in VPAC2R have little effect on network synchronization, but can improve the amplitude of the SCN network rhythms while narrowing the network entrainment range. We further show that these findings are valid across several computational network models. Thus, we identified a general design principle to achieve robust synchronization: An activating coupling agent, such as VIP, must act in-phase with the activity of core-clock promoters. More generally, the phase of coupling is as critical as the strength of coupling from the viewpoint of synchrony and entrainment.

Coupling allows robust mammalian redox circadian rhythms despite heterogeneity and noise.

Circadian clocks are endogenous oscillators present in almost all cells that drive daily rhythms in physiology and behavior. There are two mechanisms that have been proposed to explain how circadian rhythms are generated in mammalian cells: through a transcription-translation feedback loop (TTFL) and based on oxidation/reduction reactions, both of which are intrinsically stochastic and heterogeneous at the single cell level. In order to explore the emerging properties of stochastic and heterogeneous redox oscillators, we simplify a recently developed kinetic model of redox oscillations to an amplitude-phase oscillator with 'twist' (period-amplitude correlation) and subject to Gaussian noise. We show that noise and heterogeneity alone lead to fast desynchronization, and that coupling between noisy oscillators can establish robust and synchronized rhythms with amplitude expansions and tuning of the period due to twist. Coupling a network of redox oscillators to a simple model of the TTFL also contributes to synchronization, large amplitudes and fine-tuning of the period for appropriate interaction strengths. These results provide insights into how the circadian clock compensates randomness from intracellular sources and highlight the importance of noise, heterogeneity and coupling in the context of circadian oscillators.

Stable IL-2 decision making by endogenous c-Fos amounts in peripheral memory T-helper cells.

The cytokine IL-2 performs opposite functions supporting efficient immune responses and playing a key role in peripheral tolerance. Therefore, precise fine-tuning of IL-2 expression is crucial for adjusting the immune response. Combining transcription factor analysis at the single cell and the single nucleus level using flow cytometry with statistical analysis, we showed that physiological differences in the expression levels of c-Fos and NFATc2, but not of c-Jun and NF-κBp65, are limiting for the decision whether IL-2 is expressed in a strongly activated human memory T-helper (Th) cell. Variation in the expression of c-Fos leads to substantial diversity of IL-2 expression in ∼40% of the memory Th cells. The remaining cells exhibit an equally high c-Fos expression level, thereby ensuring robustness in IL-2 response within the population. These findings reveal how memory Th cells benefit from regulated variation in transcription factor expression to achieve a certain stability and variability of cytokine expression in a controlled manner.

Reverse engineering a hierarchical regulatory network downstream of oncogenic KRAS.

RAS mutations are highly relevant for progression and therapy response of human tumours, but the genetic network that ultimately executes the oncogenic effects is poorly understood. Here, we used a reverse-engineering approach in an ovarian cancer model to reconstruct KRAS oncogene-dependent cytoplasmic and transcriptional networks from perturbation experiments based on gene silencing and pathway inhibitor treatments. We measured mRNA and protein levels in manipulated cells by microarray, RT-PCR and western blot analysis, respectively. The reconstructed model revealed complex interactions among the transcriptional and cytoplasmic components, some of which were confirmed by double pertubation experiments. Interestingly, the transcription factors decomposed into two hierarchically arranged groups. To validate the model predictions, we analysed growth parameters and transcriptional deregulation in the KRAS-transformed epithelial cells. As predicted by the model, we found two functional groups among the selected transcription factors. The experiments thus confirmed the predicted hierarchical transcription factor regulation and showed that the hierarchy manifests itself in downstream gene expression patterns and phenotype.

3D-Image analysis platform monitoring relocation of pluripotency genes during reprogramming.

Nuclear organization of chromatin is an important level of genome regulation with positional changes of genes occurring during reprogramming. Inherent variability of biological specimens, wide variety of sample preparation and imaging conditions, though pose significant challenges to data analysis and comparison. Here, we describe the development of a computational image analysis toolbox overcoming biological variability hurdles by a novel single cell randomizing normalization. We performed a comparative analysis of the relationship between spatial positioning of pluripotency genes with their genomic activity and determined the degree of similarity between fibroblasts, induced pluripotent stem cells and embryonic stem cells. Our analysis revealed a preferred positioning of actively transcribed Sox2, Oct4 and Nanog away from the nuclear periphery, but not from pericentric heterochromatin. Moreover, in the silent state, we found no common nuclear localization for any of the genes. Our results suggest that the surrounding gene density hinders relocation from an internal nuclear position. Altogether, our data do not support the hypothesis that the nuclear periphery acts as a general transcriptional silencer, rather suggesting that internal nuclear localization is compatible with expression in pluripotent cells but not sufficient for expression in mouse embryonic fibroblasts. Thus, our computational approach enables comparative analysis of topological relationships in spite of stark morphological variability typical of biological data sets.