Consistent robustness analysis (CRA) identifies biologically relevant properties of regulatory network models

Background

A number of studies have previously demonstrated that “goodness of fit” is insufficient in reliably classifying the credibility of a biological model. Robustness and/or sensitivity analysis is commonly employed as a secondary method for evaluating the suitability of a particular model. The results of such analyses invariably depend on the particular parameter set tested, yet many parameter values for biological models are uncertain.

Results

Here, we propose a novel robustness analysis that aims to determine the “common robustness” of the model with multiple, biologically plausible parameter sets, rather than the local robustness for a particular parameter set. Our method is applied to two published models of the Arabidopsis circadian clock (the one-loop [1] and two-loop [2] models). The results reinforce current findings suggesting the greater reliability of the two-loop model and pinpoint the crucial role of TOC1 in the circadian network.

Conclusions

Consistent Robustness Analysis can indicate both the relative plausibility of different models and also the critical components and processes controlling each model.     

Regulation of BMAL1 Protein Stability and Circadian Function by GSK3β-Mediated Phosphorylation

Background

Circadian rhythms govern a large array of physiological and metabolic functions. To achieve plasticity in circadian regulation, proteins constituting the molecular clock machinery undergo various post-translational modifications (PTMs), which influence their activity and intracellular localization. The core clock protein BMAL1 undergoes several PTMs. Here we report that the Akt-GSK3β signaling pathway regulates BMAL1 protein stability and activity.

Principal Findings

GSK3β phosphorylates BMAL1 specifically on Ser 17 and Thr 21 and primes it for ubiquitylation. In the absence of GSK3β-mediated phosphorylation, BMAL1 becomes stabilized and BMAL1 dependent circadian gene expression is dampened. Dopamine D2 receptor mediated signaling, known to control the Akt-GSK3β pathway, influences BMAL1 stability and in vivo circadian gene expression in striatal neurons.

Conclusions

These findings uncover a previously unknown mechanism of circadian clock control. The GSK3β kinase phosphorylates BMAL1, an event that controls the stability of the protein and the amplitude of circadian oscillation. BMAL1 phosphorylation appears to be an important regulatory step in maintaining the robustness of the circadian clock.     

Circadian Timekeeping Is Disturbed in Rheumatoid Arthritis at Molecular Level

Introduction

Patients with rheumatoid arthritis (RA) have disturbances in the hypothalamic-pituitary-adrenal (HPA) axis. These are reflected in altered circadian rhythm of circulating serum cortisol, melatonin and IL-6 levels and in chronic fatigue. We hypothesized that the molecular machinery responsible for the circadian timekeeping is perturbed in RA. The aim of this study was to investigate the expression of circadian clock in RA.

Methods

Gene expression of thirteen clock genes was analyzed in the synovial membrane of RA and control osteoarthritis (OA) patients. BMAL1 protein was detected using immunohistochemistry. Cell autonomous clock oscillation was started in RA and OA synovial fibroblasts using serum shock. The effect of pro-inflammatory stimulus on clock gene expression in synovial fibroblasts was studied using IL-6 and TNF-α.

Results

Gene expression analysis disclosed disconcerted circadian timekeeping and immunohistochemistry revealed strong cytoplasmic localization of BMAL1 in RA patients. Perturbed circadian timekeeping is at least in part inflammation independent and cell autonomous, because RA synovial fibroblasts display altered circadian expression of several clock components, and perturbed circadian production of IL-6 and IL-1β after clock resetting. However, inflammatory stimulus disturbs the rhythm in cultured fibroblasts. Throughout the experiments ARNTL2 and NPAS2 appeared to be the most affected clock genes in human immune-inflammatory conditions.

Conclusion

We conclude that the molecular machinery controlling the circadian rhythm is disturbed in RA patients.     

The physiological period length of the human circadian clock in vivo is directly proportional to period in human fibroblasts

Background

Diurnal behavior in humans is governed by the period length of a circadian clock in the suprachiasmatic nuclei of the brain hypothalamus. Nevertheless, the cell-intrinsic mechanism of this clock is present in most cells of the body. We have shown previously that for individuals of extreme chronotype (“larks” and “owls”), clock properties measured in human fibroblasts correlated with extreme diurnal behavior.

Methodology/Principal Findings

In this study, we have measured circadian period in human primary fibroblasts taken from normal individuals and, for the first time, compared it directly with physiological period measured in vivo in the same subjects. Human physiological period length was estimated via the secretion pattern of the hormone melatonin in two different groups of sighted subjects and one group of totally blind subjects, each using different methods. Fibroblast period length was measured via cyclical expression of a lentivirally delivered circadian reporter. Within each group, a positive linear correlation was observed between circadian period length in physiology and in fibroblast gene expression. Interestingly, although blind individuals showed on average the same fibroblast clock properties as sighted ones, their physiological periods were significantly longer.

Conclusions/Significance

We conclude that the period of human circadian behaviour is mostly driven by cellular clock properties in normal individuals and can be approximated by measurement in peripheral cells such as fibroblasts. Based upon differences among sighted and blind subjects, we also speculate that period can be modified by prolonged unusual conditions such as the total light deprivation of blindness.     

Quantitative analysis of phase wave of gene expression in the mammalian central circadian clock network

Background

The suprachiasmatic nucleus (SCN), the master circadian clock, is a heterogeneous oscillator network, yet displays a robust synchronization dynamics. Recent single-cell bioluminescent imaging revealed temporal gradients in circadian clock gene expression in the SCN ex vivo. However, due to technical difficulty in biological approaches to elucidate the entire network structure of the SCN, characteristics of the gradient, which we refer to as phase wave, remain unknown.

Methodology/Principal Findings

We implemented new approaches, i.e., quantitative analysis and model simulation to characterize the phase waves in Per2::Luciferase clock reporter gene expression of the rat SCN slice. Our quantitative study demonstrated not only a high degree of synchronization between the neurons and regular occurrence of the phase wave propagation, but also a significant amount of phase fluctuations contained in the wave. In addition, our simulations based on local coupling model suggest that the intercellular coupling strength estimated by the model simulations is significantly higher than the critical value for generating the phase waves. Model simulations also suggest that heterogeneity of the SCN neurons is one of the main factors causing the phase wave fluctuations. Furthermore, robustness of the SCN network against dynamical noise and variation of the natural frequencies inherent in these neurons was quantitatively assessed.

Conclusions/Significance

To our knowledge, this is the first quantitative evaluation of the phase wave and further characterization of the SCN neuronal network features generating the wave i.e., intercellular synchrony, phase fluctuation, strong local coupling, heterogeneous periodicity and robustness. Our present study provides an approach, which will lead to a comprehensive understanding of mechanistic and/or biological significance of the phase wave in the central circadian oscillatory system.     

Expressions of Tight Junction Proteins Occludin and Claudin-1 Are under the Circadian Control in the Mouse Large Intestine: Implications in Intestinal Permeability and Susceptibility to Colitis

Background & Aims

The circadian clock drives daily rhythms in behavior and physiology. A recent study suggests that intestinal permeability is also under control of the circadian clock. However, the precise mechanisms remain largely unknown. Because intestinal permeability depends on tight junction (TJ) that regulates the epithelial paracellular pathway, this study investigated whether the circadian clock regulates the expression levels of TJ proteins in the intestine.

Methods

The expression levels of TJ proteins in the large intestinal epithelium and colonic permeability were analyzed every 4, 6, or 12 hours between wild-type mice and mice with a mutation of a key clock gene Period2 (Per2; mPer2^m/m). In addition, the susceptibility to dextran sodium sulfate (DSS)-induced colitis was compared between wild-type mice and mPer2^m/m mice.

Results

The mRNA and protein expression levels of Occludin and Claudin-1 exhibited daily variations in the colonic epithelium in wild-type mice, whereas they were constitutively high in mPer2^m/m mice. Colonic permeability in wild-type mice exhibited daily variations, which was inversely associated with the expression levels of Occludin and Claudin-1 proteins, whereas it was constitutively low in mPer2^m/m mice. mPer2^m/m mice were more resistant to the colonic injury induced by DSS than wild-type mice.

Conclusions

Occludin and Claudin-1 expressions in the large intestine are under the circadian control, which is associated with temporal regulation of colonic permeability and also susceptibility to colitis.     

Shift work in nurses: contribution of phenotypes and genotypes to adaptation

Background

Daily cycles of sleep/wake, hormones, and physiological processes are often misaligned with behavioral patterns during shift work, leading to an increased risk of developing cardiovascular/metabolic/gastrointestinal disorders, some types of cancer, and mental disorders including depression and anxiety. It is unclear how sleep timing, chronotype, and circadian clock gene variation contribute to adaptation to shift work.

Methods

Newly defined sleep strategies, chronotype, and genotype for polymorphisms in circadian clock genes were assessed in 388 hospital day- and night-shift nurses.

Results

Night-shift nurses who used sleep deprivation as a means to switch to and from diurnal sleep on work days (∼25%) were the most poorly adapted to their work schedule. Chronotype also influenced efficacy of adaptation. In addition, polymorphisms in CLOCK, NPAS2, PER2, and PER3 were significantly associated with outcomes such as alcohol/caffeine consumption and sleepiness, as well as sleep phase, inertia and duration in both single- and multi-locus models. Many of these results were specific to shift type suggesting an interaction between genotype and environment (in this case, shift work).

Conclusions

Sleep strategy, chronotype, and genotype contribute to the adaptation of the circadian system to an environment that switches frequently and/or irregularly between different schedules of the light-dark cycle and social/workplace time. This study of shift work nurses illustrates how an environmental “stress” to the temporal organization of physiology and metabolism can have behavioral and health-related consequences. Because nurses are a key component of health care, these findings could have important implications for health-care policy.     

PDFR and CRY signaling converge in a subset of clock neurons to modulate the amplitude and phase of circadian behavior in Drosophila

Background

To synchronize their molecular rhythms, circadian pacemaker neurons must input both external and internal timing cues and, therefore, signal integration between sensory information and internal clock status is fundamental to normal circadian physiology.

Methodology/Principal Findings

We demonstrate the specific convergence of clock-derived neuropeptide signaling with that of a deep brain photoreceptor. We report that the neuropeptide PDF receptor and the circadian photoreceptor CRYPTOCROME (CRY) are precisely co-expressed in a subset of pacemakers, and that these pathways together provide a requisite drive for circadian control of daily locomotor rhythms. These convergent signaling pathways influence the phase of rhythm generation, but also its amplitude. In the absence of both pathways, PER rhythms were greatly reduced in only those specific pacemakers that receive convergent inputs and PER levels remained high in the nucleus throughout the day. This suggested a large-scale dis-regulation of the pacemaking machinery. Behavioral rhythms were likewise disrupted: in light∶dark conditions they were aberrant, and under constant dark conditions, they were lost.

Conclusions/Significance

We speculate that the convergence of environmental and clock-derived signals may produce a coincident detection of light, synergistic responses to it, and thus more accurate and more efficient re-setting properties.     

Intramolecular Regulation of Phosphorylation Status of the Circadian Clock Protein KaiC

Background

KaiC, a central clock protein in cyanobacteria, undergoes circadian oscillations between hypophosphorylated and hyperphosphorylated forms in vivo and in vitro. Structural analyses of KaiC crystals have identified threonine and serine residues in KaiC at three residues (T426, S431, and T432) as potential sites at which KaiC is phosphorylated; mutation of any of these three sites to alanine abolishes rhythmicity, revealing an essential clock role for each residue separately and for KaiC phosphorylation in general. Mass spectrometry studies confirmed that the S431 and T432 residues are key phosphorylation sites, however, the role of the threonine residue at position 426 was not clear from the mass spectrometry measurements.

Methodology and Principal Findings

Mutational approaches and biochemical analyses of KaiC support a key role for T426 in control of the KaiC phosphorylation status in vivo and in vitro and demonstrates that alternative amino acids at residue 426 dramatically affect KaiC''s properties in vivo and in vitro, especially genetic dominance/recessive relationships, KaiC dephosphorylation, and the formation of complexes of KaiC with KaiA and KaiB. These mutations alter key circadian properties, including period, amplitude, robustness, and temperature compensation. Crystallographic analyses indicate that the T426 site is phosphorylatible under some conditions, and in vitro phosphorylation assays of KaiC demonstrate labile phosphorylation of KaiC when the primary S431 and T432 sites are blocked.

Conclusions and Significance

T426 is a crucial site that regulates KaiC phosphorylation status in vivo and in vitro and these studies underscore the importance of KaiC phosphorylation status in the essential cyanobacterial circadian functions. The regulatory roles of these phosphorylation sites–including T426–within KaiC enhance our understanding of the molecular mechanism underlying circadian rhythm generation in cyanobacteria.     

Genetic adaptation of the human circadian clock to day-length latitudinal variations and relevance for affective disorders

Background

The temporal coordination of biological processes into daily cycles is a common feature of most living organisms. In humans, disruption of circadian rhythms is commonly observed in psychiatric diseases, including schizophrenia, bipolar disorder, depression and autism. Light therapy is the most effective treatment for seasonal affective disorder and circadian-related treatments sustain antidepressant response in bipolar disorder patients. Day/night cycles represent a major circadian synchronizing signal and vary widely with latitude.

Results

We apply a geographically explicit model to show that out-of-Africa migration, which led humans to occupy a wide latitudinal area, affected the evolutionary history of circadian regulatory genes. The SNPs we identify using this model display consistent signals of natural selection using tests based on population genetic differentiation and haplotype homozygosity. Signals of natural selection driven by annual photoperiod variation are detected for schizophrenia, bipolar disorder, and restless leg syndrome risk variants, in line with the circadian component of these conditions.

Conclusions

Our results suggest that human populations adapted to life at different latitudes by tuning their circadian clock systems. This process also involves risk variants for neuropsychiatric conditions, suggesting possible genetic modulators for chronotherapies and candidates for interaction analysis with photoperiod-related environmental variables, such as season of birth, country of residence, shift-work or lifestyle habits.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0499-7) contains supplementary material, which is available to authorized users.     

Differential contribution of rod and cone circadian clocks in driving retinal melatonin rhythms in Xenopus

Background

Although an endogenous circadian clock located in the retinal photoreceptor layer governs various physiological events including melatonin rhythms in Xenopus laevis, it remains unknown which of the photoreceptors, rod and/or cone, is responsible for the circadian regulation of melatonin release.

Methodology/Principal Findings

We selectively disrupted circadian clock function in either the rod or cone photoreceptor cells by generating transgenic Xenopus tadpoles expressing a dominant-negative CLOCK (XCLΔQ) under the control of a rod or cone-specific promoter. Eyecup culture and continuous melatonin measurement revealed that circadian rhythms of melatonin release were abolished in a majority of the rod-specific XCLΔQ transgenic tadpoles, although the percentage of arrhythmia was lower than that of transgenic tadpole eyes expressing XCLΔQ in both rods and cones. In contrast, whereas a higher percentage of arrhythmia was observed in the eyes of the cone-specific XCLΔQ transgenic tadpoles compare to wild-type counterparts, the rate was significantly lower than in rod-specific transgenics. The levels of the transgene expression were comparable between these two different types of transgenics. In addition, the average overall melatonin levels were not changed in the arrhythmic eyes, suggesting that CLOCK does not affect absolute levels of melatonin, only its temporal expression pattern.

Conclusions/Significance

These results suggest that although the Xenopus retina is made up of approximately equal numbers of rods and cones, the circadian clocks in the rod cells play a dominant role in driving circadian melatonin rhythmicity in the Xenopus retina, although some contribution of the clock in cone cells cannot be excluded.     

Phase computations and phase models for discrete molecular oscillators

Background

Biochemical oscillators perform crucial functions in cells, e.g., they set up circadian clocks. The dynamical behavior of oscillators is best described and analyzed in terms of the scalar quantity, phase. A rigorous and useful definition for phase is based on the so-called isochrons of oscillators. Phase computation techniques for continuous oscillators that are based on isochrons have been used for characterizing the behavior of various types of oscillators under the influence of perturbations such as noise.

Results

In this article, we extend the applicability of these phase computation methods to biochemical oscillators as discrete molecular systems, upon the information obtained from a continuous-state approximation of such oscillators. In particular, we describe techniques for computing the instantaneous phase of discrete, molecular oscillators for stochastic simulation algorithm generated sample paths. We comment on the accuracies and derive certain measures for assessing the feasibilities of the proposed phase computation methods. Phase computation experiments on the sample paths of well-known biological oscillators validate our analyses.

Conclusions

The impact of noise that arises from the discrete and random nature of the mechanisms that make up molecular oscillators can be characterized based on the phase computation techniques proposed in this article. The concept of isochrons is the natural choice upon which the phase notion of oscillators can be founded. The isochron-theoretic phase computation methods that we propose can be applied to discrete molecular oscillators of any dimension, provided that the oscillatory behavior observed in discrete-state does not vanish in a continuous-state approximation. Analysis of the full versatility of phase noise phenomena in molecular oscillators will be possible if a proper phase model theory is developed, without resorting to such approximations.     

Maternal feeding controls fetal biological clock

Background

It is widely accepted that circadian physiological rhythms of the fetus are affected by oscillators in the maternal brain that are coupled to the environmental light-dark (LD) cycle.

Methodology/Principal Findings

To study the link between fetal and maternal biological clocks, we investigated the effects of cycles of maternal food availability on the rhythms of Per1 gene expression in the fetal suprachiasmatic nucleus (SCN) and liver using a transgenic rat model whose tissues express luciferase in vitro. Although the maternal SCN remained phase-locked to the LD cycle, maternal restricted feeding phase-advanced the fetal SCN and liver by 5 and 7 hours respectively within the 22-day pregnancy.

Conclusions/Significance

Our results demonstrate that maternal feeding entrains the fetal SCN and liver independently of both the maternal SCN and the LD cycle. This indicates that maternal-feeding signals can be more influential for the fetal SCN and particular organ oscillators than hormonal signals controlled by the maternal SCN, suggesting the importance of a regular maternal feeding schedule for appropriate fetal molecular clockwork during pregnancy.     

Characterization of Circadian Behavior in the Starlet Sea Anemone,Nematostella vectensis

Background

Although much is known about how circadian systems control daily cycles in the physiology and behavior of Drosophila and several vertebrate models, marine invertebrates have often been overlooked in circadian rhythms research. This study focuses on the starlet sea anemone, Nematostella vectensis, a species that has received increasing attention within the scientific community for its potential as a model research organism. The recently sequenced genome of N. vectensis makes it an especially attractive model for exploring the molecular evolution of circadian behavior. Critical behavioral data needed to correlate gene expression patterns to specific behaviors are currently lacking in N. vectensis.

Methodology/Principal Findings

To detect the presence of behavioral oscillations in N. vectensis, locomotor activity was evaluated using an automated system in an environmentally controlled chamber. Animals exposed to a 24 hr photoperiod (12 hr light: 12 hr dark) exhibited locomotor behavior that was both rhythmic and predominantly nocturnal. The activity peak occurred in the early half of the night with a 2-fold increase in locomotion. Upon transfer to constant lighting conditions (constant light or constant dark), an approximately 24 hr rhythm persisted in most animals, suggesting that the rhythm is controlled by an endogenous circadian mechanism. Fourier analysis revealed the presence of multiple peaks in some animals suggesting additional rhythmic components could be present. In particular, an approximately 12 hr oscillation was often observed. The nocturnal increase in generalized locomotion corresponded to a 24 hr oscillation in animal elongation.

Conclusions/Significance

These data confirm the presence of a light-entrainable circadian clock in Nematostella vectensis. Additional components observed in some individuals indicate that an endogenous clock of approximately 12 hr frequency may also be present. By describing rhythmic locomotor behavior in N. vectensis, we have made important progress in developing the sea anemone as a model organism for circadian rhythm research.     

RNA-Interference Knockdown of Drosophila Pigment Dispersing Factor in Neuronal Subsets: The Anatomical Basis of a Neuropeptide's Circadian Functions

Background

In animals, neuropeptide signaling is an important component of circadian timekeeping. The neuropeptide pigment dispersing factor (PDF) is required for several aspects of circadian activity rhythms in Drosophila.

Methodology/Principal Findings

Here we investigate the anatomical basis for PDF''s various circadian functions by targeted PDF RNA-interference in specific classes of Drosophila neuron. We demonstrate that PDF is required in the ventro-lateral neurons (vLNs) of the central brain and not in the abdominal ganglion for normal activity rhythms. Differential knockdown of PDF in the large or small vLNs indicates that PDF from the small vLNs is likely responsible for the maintenance of free-running activity rhythms and that PDF is not required in the large vLNs for normal behavior. PDF''s role in setting the period of free-running activity rhythms and the proper timing of evening activity under light:dark cycles emanates from both subtypes of vLN, since PDF in either class of vLN was sufficient for these aspects of behavior.

Conclusions/Significance

These results reveal the neuroanatomical basis PDF''s various circadian functions and refine our understanding of the clock neuron circuitry of Drosophila.