Taking the lag out of jet lag through model-based schedule design

Travel across multiple time zones results in desynchronization of environmental time cues and the sleep–wake schedule from their normal phase relationships with the endogenous circadian system. Circadian misalignment can result in poor neurobehavioral performance, decreased sleep efficiency, and inappropriately timed physiological signals including gastrointestinal activity and hormone release. Frequent and repeated transmeridian travel is associated with long-term cognitive deficits, and rodents experimentally exposed to repeated schedule shifts have increased death rates. One approach to reduce the short-term circadian, sleep–wake, and performance problems is to use mathematical models of the circadian pacemaker to design countermeasures that rapidly shift the circadian pacemaker to align with the new schedule. In this paper, the use of mathematical models to design sleep–wake and countermeasure schedules for improved performance is demonstrated. We present an approach to designing interventions that combines an algorithm for optimal placement of countermeasures with a novel mode of schedule representation. With these methods, rapid circadian resynchrony and the resulting improvement in neurobehavioral performance can be quickly achieved even after moderate to large shifts in the sleep–wake schedule. The key schedule design inputs are endogenous circadian period length, desired sleep–wake schedule, length of intervention, background light level, and countermeasure strength. The new schedule representation facilitates schedule design, simulation studies, and experiment design and significantly decreases the amount of time to design an appropriate intervention. The method presented in this paper has direct implications for designing jet lag, shift-work, and non-24-hour schedules, including scheduling for extreme environments, such as in space, undersea, or in polar regions

Can People Sleep Too Much? Effects of Extended Sleep Opportunity on Sleep Duration and Timing

Many people are concerned about whether they are getting “enough” sleep, and if they can “sleep too much.” These concerns can be approached scientifically using experiments probing long-term (i.e., multi-night) sleep homeostatic processes, since homeostatic processes move the system toward its physiological setpoint (i.e., between “not enough” and “too much”). We analyzed sleep data from two human studies with sleep opportunities much longer than people usually stay in bed (i.e., conditions in which sleep homeostatic responses could be documented): sleep opportunities were 14–16 h per day for 3–28 days. Across the nights of the extended sleep opportunities, total sleep duration, Rapid Eye Movement (REM) sleep duration and non-REM sleep durations decreased and sleep latency increased. Multiple nights were required to reach approximately steady-state values. These results suggest a multi-day homeostatic sleep process responding to self-selected insufficient sleep duration prior to the study. Once steady state-values were reached, there were large night-to-night variations in total sleep time and other sleep metrics. Our results therefore answer these concerns about sleep amount and are important for understanding the basic physiology of sleep and for two sleep-related topics: (i) the inter-individual and intra-individual variability are relevant to understanding “normal” sleep patterns and for people with insomnia and (ii) the multiple nights of sleep required for recovery from insufficient sleep from self-selected sleep loss is important for public health and other efforts for reducing the adverse effects of sleep loss on multiple areas of physiology

Biological Time Series Analysis Using a Context Free Language: Applicability to Pulsatile Hormone Data

National Space Biomedical Research Institute (NASA NCC 9-58 HFP01603)National Space Biomedical Research Institute (NASA NCC 9-58 HPF00405)National Institutes of Health (U.S.) (NIH NCRR-GCRC-M01-RR-02635)United States. Air Force Office of Scientific Research (AFOSR F49620-95-1-0388)United States. Air Force Office of Scientific Research (AFOSR FA9550-06-0080)National Institutes of Health (U.S.) (NIH P01-AG09975)National Institutes of Health (U.S.) (NIH T32 HL07901-10)National Institutes of Health (U.S.) (NIH F31-GM095340-01)National Institutes of Health (U.S.) (NIH K24-HL105664)National Institutes of Health (U.S.) (NIH K02-HD045459)National Institutes of Health (U.S.) (NIH RC2-HL101340)National Institutes of Health (U.S.) (NIH R01-AR43130)National Institutes of Health (U.S.) (NIH K24-HL103845)National Institutes of Health (U.S.) (NIH R01-MH071847)National Institutes of Health (U.S.) (NIH R01 HL098433)National Institutes of Health (U.S.) (NIH R01 HL098433-02S1

Biological Time Series Analysis Using a Context Free Language: Applicability to Pulsatile Hormone Data

We present a novel approach for analyzing biological time-series data using a context-free language (CFL) representation that allows the extraction and quantification of important features from the time-series. This representation results in Hierarchically AdaPtive (HAP) analysis, a suite of multiple complementary techniques that enable rapid analysis of data and does not require the user to set parameters. HAP analysis generates hierarchically organized parameter distributions that allow multi-scale components of the time-series to be quantified and includes a data analysis pipeline that applies recursive analyses to generate hierarchically organized results that extend traditional outcome measures such as pharmacokinetics and inter-pulse interval. Pulsicons, a novel text-based time-series representation also derived from the CFL approach, are introduced as an objective qualitative comparison nomenclature. We apply HAP to the analysis of 24 hours of frequently sampled pulsatile cortisol hormone data, which has known analysis challenges, from 14 healthy women. HAP analysis generated results in seconds and produced dozens of figures for each participant. The results quantify the observed qualitative features of cortisol data as a series of pulse clusters, each consisting of one or more embedded pulses, and identify two ultradian phenotypes in this dataset. HAP analysis is designed to be robust to individual differences and to missing data and may be applied to other pulsatile hormones. Future work can extend HAP analysis to other time-series data types, including oscillatory and other periodic physiological signals

Stability of the timing of food intake at daily and monthly timescales in young adults

Cross-sectional observations have shown that the timing of eating may be important for health-related outcomes. Here we examined the stability of eating timing, using both clock hour and relative circadian time, across one semester (n = 14) at daily and monthly time-scales. At three time points ~ 1 month apart, circadian phase was determined during an overnight in-laboratory visit and eating was photographically recorded for one week to assess timing and composition. Day-to-day stability was measured using the Composite Phase Deviation (deviation from a perfectly regular pattern) and intraclass correlation coefficients (ICC) were used to determine individual stability across months (weekly average compared across months). Day-to-day clock timing of caloric events had poor stability within individuals (~ 3-h variation; ICC = 0.12–0.34). The timing of eating was stable across months (~ 1-h variation, ICCs ranging from 0.54–0.63), but less stable across months when measured relative to circadian timing (ICC = 0.33–0.41). Our findings suggest that though day-to-day variability in the timing of eating has poor stability, the timing of eating measured for a week is stable across months within individuals. This indicates two relevant timescales: a monthly timescale with more stability in eating timing than a daily timescale. Thus, a single day’s food documentation may not represent habitual (longer timescale) patterns

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

Impact of Common Diabetes Risk Variant in MTNR1B

The risk of type 2 diabetes (T2D) is increased by abnormalities in sleep quantity and quality, circadian alignment, and melatonin regulation. A common genetic variant in a receptor for the circadian-regulated hormone melatonin (MTNR1B) is associated with increased fasting blood glucose and risk of T2D, but whether sleep or circadian disruption mediates this risk is unknown. We aimed to test if MTNR1B diabetes risk variant rs10830963 associates with measures of sleep or circadian physiology in intensive in-laboratory protocols (n = 58–96) or cross-sectional studies with sleep quantity and quality and timing measures from self-report (n = 4,307–10,332), actigraphy (n = 1,513), or polysomnography (n = 3,021). In the in-laboratory studies, we found a significant association with a substantially longer duration of elevated melatonin levels (41 min) and delayed circadian phase of dim-light melatonin offset (1.37 h), partially mediated through delayed offset of melatonin synthesis. Furthermore, increased T2D risk in MTNR1B risk allele carriers was more pronounced in early risers versus late risers as determined by 7 days of actigraphy. Our results provide the surprising insight that the MTNR1B risk allele influences dynamics of melatonin secretion, generating a novel hypothesis that the MTNR1B risk allele may extend the duration of endogenous melatonin production later into the morning and that early waking may magnify the diabetes risk conferred by the risk allele

Taking the lag out of jet lag through model-based schedule design.

Travel across multiple time zones results in desynchronization of environmental time cues and the sleep-wake schedule from their normal phase relationships with the endogenous circadian system. Circadian misalignment can result in poor neurobehavioral performance, decreased sleep efficiency, and inappropriately timed physiological signals including gastrointestinal activity and hormone release. Frequent and repeated transmeridian travel is associated with long-term cognitive deficits, and rodents experimentally exposed to repeated schedule shifts have increased death rates. One approach to reduce the short-term circadian, sleep-wake, and performance problems is to use mathematical models of the circadian pacemaker to design countermeasures that rapidly shift the circadian pacemaker to align with the new schedule. In this paper, the use of mathematical models to design sleep-wake and countermeasure schedules for improved performance is demonstrated. We present an approach to designing interventions that combines an algorithm for optimal placement of countermeasures with a novel mode of schedule representation. With these methods, rapid circadian resynchrony and the resulting improvement in neurobehavioral performance can be quickly achieved even after moderate to large shifts in the sleep-wake schedule. The key schedule design inputs are endogenous circadian period length, desired sleep-wake schedule, length of intervention, background light level, and countermeasure strength. The new schedule representation facilitates schedule design, simulation studies, and experiment design and significantly decreases the amount of time to design an appropriate intervention. The method presented in this paper has direct implications for designing jet lag, shift-work, and non-24-hour schedules, including scheduling for extreme environments, such as in space, undersea, or in polar regions