Analysis Method and Experimental Conditions Affect Computed Circadian Phase from Melatonin Data

Accurate determination of circadian phase is necessary for research and clinical purposes because of the influence of the master circadian pacemaker on multiple physiologic functions. Melatonin is presently the most accurate marker of the activity of the human circadian pacemaker. Current methods of analyzing the plasma melatonin rhythm can be grouped into three categories: curve-fitting, threshold-based and physiologically-based linear differential equations. To determine which method provides the most accurate assessment of circadian phase, we compared the ability to fit the data and the variability of phase estimates for seventeen different markers of melatonin phase derived from these methodological categories. We used data from three experimental conditions under which circadian rhythms - and therefore calculated melatonin phase - were expected to remain constant or progress uniformly. Melatonin profiles from older subjects and subjects with lower melatonin amplitude were less likely to be fit by all analysis methods. When circadian drift over multiple study days was algebraically removed, there were no significant differences between analysis methods of melatonin onsets (P = 0.57), but there were significant differences between those of melatonin offsets (P<0.0001). For a subset of phase assessment methods, we also examined the effects of data loss on variability of phase estimates by systematically removing data in 2-hour segments. Data loss near onset of melatonin secretion differentially affected phase estimates from the methods, with some methods incorrectly assigning phases too early while other methods assigning phases too late; missing data at other times did not affect analyses of the melatonin profile. We conclude that melatonin data set characteristics, including amplitude and completeness of data collection, differentially affect the results depending on the melatonin analysis method used

Deconvolution of Serum Cortisol Levels by Using Compressed Sensing

The pulsatile release of cortisol from the adrenal glands is controlled by a hierarchical system that involves corticotropin releasing hormone (CRH) from the hypothalamus, adrenocorticotropin hormone (ACTH) from the pituitary, and cortisol from the adrenal glands. Determining the number, timing, and amplitude of the cortisol secretory events and recovering the infusion and clearance rates from serial measurements of serum cortisol levels is a challenging problem. Despite many years of work on this problem, a complete satisfactory solution has been elusive. We formulate this question as a non-convex optimization problem, and solve it using a coordinate descent algorithm that has a principled combination of (i) compressed sensing for recovering the amplitude and timing of the secretory events, and (ii) generalized cross validation for choosing the regularization parameter. Using only the observed serum cortisol levels, we model cortisol secretion from the adrenal glands using a second-order linear differential equation with pulsatile inputs that represent cortisol pulses released in response to pulses of ACTH. Using our algorithm and the assumption that the number of pulses is between 15 to 22 pulses over 24 hours, we successfully deconvolve both simulated datasets and actual 24-hr serum cortisol datasets sampled every 10 minutes from 10 healthy women. Assuming a one-minute resolution for the secretory events, we obtain physiologically plausible timings and amplitudes of each cortisol secretory event with R[superscript 2] above 0.92. Identification of the amplitude and timing of pulsatile hormone release allows (i) quantifying of normal and abnormal secretion patterns towards the goal of understanding pathological neuroendocrine states, and (ii) potentially designing optimal approaches for treating hormonal disorders.National Science Foundation (U.S.). Graduate Research Fellowship ProgramNational Institutes of Health (U.S.) (NIH DP1 OD003646)National Science Foundation (U.S.) (0836720)National Science Foundation (U.S.). Office of Emerging Frontiers in Research and Innovation (EFRI-0735956

Amplitude Reduction and Phase Shifts of Melatonin, Cortisol and Other Circadian Rhythms after a Gradual Advance of Sleep and Light Exposure in Humans

Background: The phase and amplitude of rhythms in physiology and behavior are generated by circadian oscillators and entrained to the 24-h day by exposure to the light-dark cycle and feedback from the sleep-wake cycle. The extent to which the phase and amplitude of multiple rhythms are similarly affected during altered timing of light exposure and the sleepwake cycle has not been fully characterized. Methodology/Principal Findings: We assessed the phase and amplitude of the rhythms of melatonin, core body temperature, cortisol, alertness, performance and sleep after a perturbation of entrainment by a gradual advance of the sleep-wake schedule (10 h in 5 days) and associated light-dark cycle in 14 healthy men. The light-dark cycle consisted either of moderate intensity ‘room ’ light (,90–150 lux) or moderate light supplemented with bright light (,10,000 lux) for 5 to 8 hours following sleep. After the advance of the sleep-wake schedule in moderate light, no significant advance of the melatonin rhythm was observed whereas, after bright light supplementation the phase advance was 8.1 h (SEM 0.7 h). Individual differences in phase shifts correlated across variables. The amplitude of the melatonin rhythm assessed under constant conditions was reduced after moderate light by 54 % (17–94%) and after bright light by 52 % (range 12–84%), as compared to the amplitude at baseline in the presence of a sleep-wake cycle. Individual differences in amplitude reduction of the melatonin rhythm correlated with the amplitude of body temperature, cortisol and alertness

Power Law versus Exponential State Transition Dynamics: Application to Sleep-Wake Architecture

BACKGROUND: Despite the common experience that interrupted sleep has a negative impact on waking function, the features of human sleep-wake architecture that best distinguish sleep continuity versus fragmentation remain elusive. In this regard, there is growing interest in characterizing sleep architecture using models of the temporal dynamics of sleep-wake stage transitions. In humans and other mammals, the state transitions defining sleep and wake bout durations have been described with exponential and power law models, respectively. However, sleep-wake stage distributions are often complex, and distinguishing between exponential and power law processes is not always straightforward. Although mono-exponential distributions are distinct from power law distributions, multi-exponential distributions may in fact resemble power laws by appearing linear on a log-log plot. METHODOLOGY/PRINCIPAL FINDINGS: To characterize the parameters that may allow these distributions to mimic one another, we systematically fitted multi-exponential-generated distributions with a power law model, and power law-generated distributions with multi-exponential models. We used the Kolmogorov-Smirnov method to investigate goodness of fit for the "incorrect" model over a range of parameters. The "zone of mimicry" of parameters that increased the risk of mistakenly accepting power law fitting resembled empiric time constants obtained in human sleep and wake bout distributions. CONCLUSIONS/SIGNIFICANCE: Recognizing this uncertainty in model distinction impacts interpretation of transition dynamics (self-organizing versus probabilistic), and the generation of predictive models for clinical classification of normal and pathological sleep architecture

Obstructive Sleep Apnea Alters Sleep Stage Transition Dynamics

Enhanced characterization of sleep architecture, compared with routine polysomnographic metrics such as stage percentages and sleep efficiency, may improve the predictive phenotyping of fragmented sleep. One approach involves using stage transition analysis to characterize sleep continuity.We analyzed hypnograms from Sleep Heart Health Study (SHHS) participants using the following stage designations: wake after sleep onset (WASO), non-rapid eye movement (NREM) sleep, and REM sleep. We show that individual patient hypnograms contain insufficient number of bouts to adequately describe the transition kinetics, necessitating pooling of data. We compared a control group of individuals free of medications, obstructive sleep apnea (OSA), medical co-morbidities, or sleepiness (n = 374) with mild (n = 496) or severe OSA (n = 338). WASO, REM sleep, and NREM sleep bout durations exhibited multi-exponential temporal dynamics. The presence of OSA accelerated the "decay" rate of NREM and REM sleep bouts, resulting in instability manifesting as shorter bouts and increased number of stage transitions. For WASO bouts, previously attributed to a power law process, a multi-exponential decay described the data well. Simulations demonstrated that a multi-exponential process can mimic a power law distribution.OSA alters sleep architecture dynamics by decreasing the temporal stability of NREM and REM sleep bouts. Multi-exponential fitting is superior to routine mono-exponential fitting, and may thus provide improved predictive metrics of sleep continuity. However, because a single night of sleep contains insufficient transitions to characterize these dynamics, extended monitoring of sleep, probably at home, would be necessary for individualized clinical application

Harmonics of Circadian Gene Transcription in Mammals

The circadian clock is a molecular and cellular oscillator found in most mammalian tissues that regulates rhythmic physiology and behavior. Numerous investigations have addressed the contribution of circadian rhythmicity to cellular, organ, and organismal physiology. We recently developed a method to look at transcriptional oscillations with unprecedented precision and accuracy using high-density time sampling. Here, we report a comparison of oscillating transcription from mouse liver, NIH3T3, and U2OS cells. Several surprising observations resulted from this study, including a 100-fold difference in the number of cycling transcripts in autonomous cellular models of the oscillator versus tissues harvested from intact mice. Strikingly, we found two clusters of genes that cycle at the second and third harmonic of circadian rhythmicity in liver, but not cultured cells. Validation experiments show that 12-hour oscillatory transcripts occur in several other peripheral tissues as well including heart, kidney, and lungs. These harmonics are lost ex vivo, as well as under restricted feeding conditions. Taken in sum, these studies illustrate the importance of time sampling with respect to multiple testing, suggest caution in use of autonomous cellular models to study clock output, and demonstrate the existence of harmonics of circadian gene expression in the mouse

Cerebrospinal fluid sodium rhythms

Background: Cerebrospinal fluid (CSF) sodium levels have been reported to rise during episodic migraine. Since migraine frequently starts in early morning or late afternoon, we hypothesized that natural sodium chronobiology may predispose susceptible persons when extracellular CSF sodium increases. Since no mammalian brain sodium rhythms are known, we designed a study of healthy humans to test if cation rhythms exist in CSF. Methods: Lumbar CSF was collected every ten minutes at 0.1 mL/min for 24 h from six healthy participants. CSF sodium and potassium concentrations were measured by ion chromatography, total protein by fluorescent spectrometry, and osmolarity by freezing point depression. We analyzed cation and protein distributions over the 24 h period and spectral and permutation tests to identify significant rhythms. We applied the False Discovery Rate method to adjust significance levels for multiple tests and Spearman correlations to compare sodium fluctuations with potassium, protein, and osmolarity. Results: The distribution of sodium varied much more than potassium, and there were statistically significant rhythms at 12 and 1.65 h periods. Curve fitting to the average time course of the mean sodium of all six subjects revealed the lowest sodium levels at 03.20 h and highest at 08.00 h, a second nadir at 09.50 h and a second peak at 18.10 h. Sodium levels were not correlated with potassium or protein concentration, or with osmolarity. Conclusion: These CSF rhythms are the first reports of sodium chronobiology in the human nervous system. The results are consistent with our hypothesis that rising levels of extracellular sodium may contribute to the timing of migraine onset. The physiological importance of sodium in the nervous system suggests that these rhythms may have additional repercussions on ultradian functions

Changes in the Circadian Rhythm in Patients with Primary Glaucoma

Purpose The current study was undertaken to investigate whether glaucoma affects the sleep quality and whether there is any difference between patients with primary glaucoma (primary open angle glaucoma, POAG and primary angle-closure glaucoma, PACG) and healthy subjects, using a validated self-rated questionnaire, the Pittsburgh Sleep Quality Index (PSQI). Methods The sleep quality of patients with POAG and PACG was tested against normal controls. Subjects were divided into three sub-groups according to age. Differences in the frequency of sleep disturbances (PSQI score >7) were assessed. The differences of sleep quality within the three groups and within the POAG group depending on the patients’ intraocular pressure (IOP) and impairment of visual field (VF) were also studied. Results 92 POAG patients, 48 PACG patients and 199 controls were included. Sleep quality declined with age in control and POAG group (tendency chi-square, P0.05). No significant differences were found in POAG group between patients with a highest IOP in daytime and at nighttime (χ2-test, P>0.05). Conclusions The prevalence of sleep disorders was higher in patients with POAG and PACG than in controls. PACG patients seemed to have a more serious problem of sleep disorders than POAG patients between 61 to 80 years old. No correlation was found between the prevalence of sleep disorders and impairment of VF or the time when POAG patients showed a highest IOP

Dose finding of melatonin for chronic idiopathic childhood sleep onset insomnia: an RCT

Contains fulltext : 86695.pdf (publisher's version ) (Open Access)Rationale Pharmacokinetics of melatonin in children might differ from that in adults. Objectives This study aims to establish a dose–response relationship for melatonin in advancing dim light melatonin onset (DLMO), sleep onset (SO), and reducing sleep onset latency (SOL) in children between 6 and 12 years with chronic sleep onset insomnia (CSOI). Methods The method used for this study is the randomized, placebo-controlled double-blind trial. Children with CSOI (n=72) received either melatonin 0.05, 0.1, and 0.15 mg/kg or placebo during 1 week. Sleep was assessed with log and actigraphy during this week and the week before. Outcomes were the shifts in DLMO, SO, and SOL. Results Treatment with melatonin significantly advanced SO and DLMO by approximately 1 h and decreased SOL by 35 min. Within the three melatonin groups, effect size was not different, but the circadian time of administration (TOA) correlated significantly with treatment effect on DLMO (rs=-0.33, p=0.022) and SO (rs=-0.38, p=0.004), whereas clock TOA was correlated with SO shift (r=-0.35, p=0.006) and not with DLMO shift. Conclusions No dose–response relationship of melatonin with SO, SOL, and DLMO is found within a dosage range of 0.05–0.15 mg/kg. The effect of exogenous melatonin on SO, SOL, and DLMO increases with an earlier circadian TOA. The soporific effects of melatonin enhance the SO shift. This study demonstrates that melatonin for treatment of CSOI in children is effective in a dosage of 0.05 mg/kg given at least 1 to 2 h before DLMO and before desired bedtime.13 p

Plasticity of the Intrinsic Period of the Human Circadian Timing System

Human expeditions to Mars will require adaptation to the 24.65-h Martian solar day-night cycle (sol), which is outside the range of entrainment of the human circadian pacemaker under lighting intensities to which astronauts are typically exposed. Failure to entrain the circadian time-keeping system to the desired rest-activity cycle disturbs sleep and impairs cognitive function. Furthermore, differences between the intrinsic circadian period and Earth's 24-h light-dark cycle underlie human circadian rhythm sleep disorders, such as advanced sleep phase disorder and non-24-hour sleep-wake disorders. Therefore, first, we tested whether exposure to a model-based lighting regimen would entrain the human circadian pacemaker at a normal phase angle to the 24.65-h Martian sol and to the 23.5-h day length often required of astronauts during short duration space exploration. Second, we tested here whether such prior entrainment to non-24-h light-dark cycles would lead to subsequent modification of the intrinsic period of the human circadian timing system. Here we show that exposure to moderately bright light (∼450 lux; ∼1.2 W/m2) for the second or first half of the scheduled wake episode is effective for entraining individuals to the 24.65-h Martian sol and a 23.5-h day length, respectively. Estimations of the circadian periods of plasma melatonin, plasma cortisol, and core body temperature rhythms collected under forced desynchrony protocols revealed that the intrinsic circadian period of the human circadian pacemaker was significantly longer following entrainment to the Martian sol as compared to following entrainment to the 23.5-h day. The latter finding of after-effects of entrainment reveals for the first time plasticity of the period of the human circadian timing system. Both findings have important implications for the treatment of circadian rhythm sleep disorders and human space exploration