University of Groningen Social Jetlag and Obesity Roenneberg, Till; Allebrandt, Karla V.; Merrow, Martha; Vetter, Celine Published in: Current Biology DOI: 10.1016/j.cub.2012.03.038 IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please check the document version below. Document Version Publisher's PDF, also known as Version of record Publication date: 2012 Link to publication in University of Groningen/UMCG research database Citation for published version (APA): Roenneberg, T., Allebrandt, K. V., Merrow, M., & Vetter, C. (2012). Social Jetlag and Obesity. Current Biology, 22(10), 939-943. DOI: 10.1016/j.cub.2012.03.038 Copyright Other than for strictly personal use, it is not permitted to download or to forward/distribute the text or part of it without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license (like Creative Commons). Take-down policy If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim. Downloaded from the University of Groningen/UMCG research database (Pure): http://www.rug.nl/research/portal. For technical reasons the number of authors shown on this cover page is limited to 10 maximum. Download date: 15-02-2018
Current Biology, Volume 22 Supplemental Information Social Jetlag and Obesity Till Roenneberg, Karla V. Allebrandt, Martha Merrow, and Céline Vetter Supplemental Inventory 1. Supplemental Figures and Tables Figure S1, related to Figure 1B Figure S2 Figure S3, related to Figure 2 Table S1, related to Figure 2 Table S2, related to Supplemental Experimental Procedures Table S3, related to Figure 3 Table S4, related to Figure 2 2. Supplemental Experimental Procedures 3. Supplemental References
Figure S1. Social Jetlag (A) Six-week long sleep-log of an extremely late chronotype (MSF 7), exemplifying the typical scalloping between sleep time on workdays and on free days (e.g., weekends). The top horizontal axis represents local time and the vertical axis represent the days of the sleep-log. The coloured bars sow the timing and duration of sleep on the respective days (red: workdays; green: free days). The difference between the mid-sleep point on free days, MSF) and that on workdays (MSW) is used to quantify social jetlag. Note that sleep on workdays in this late type is interrupted by the alarm clock (constant sleep end at around 7 a.m., corresponding to the MSF point of this subject). Although this is an extreme example of social jetlag (due to the late chronotype in combination with an early work start), the majority of the population shows similar patterns. (B) Distribution of social jetlag (MSF-MSW) in the population. Note that early chronotypes generally sleep outside of their circadian sleep window on evenings before free days due to the social pressure of a latechronotype society that makes them go to bed after their biological bedtime (hence negative social jetlag values). The majority of the population experiences the discrepancy between biological and social time on workdays (as shown in A), leading to positive values of social jetlag. (C) Age- and sex-dependencies of social jetlag (see Evaluation of the MCTQ below). Circles represent average values of age groups for both sexes (ages 10-65: 1 year bins; >65: 5 year bins; note that the age group 75 has no entries for social jetlag). A univariate ANOVAs with age and sex as covariates shows that these two variables influence social jetlag: age, F(1, 64,107) = 13,044.64, p < 0.0001, partial η2 = 0.17 and sex, F(1, 64,107)= 362,5, p < 0.0001, partial η2 = 0.01. Curves are polynomial fits (black: both sexes; red: females; blue: males). Vertical lines represent standard errors of the mean (±SEM; in most cases they are smaller than the respective symbols).
Figure S2. Age- and Sex-Dependencies of BMI Circles represent average values of age groups (1 year bins). A univariate ANOVAs with age and sex as covariates shows that these two variables influence BMI: age, F(1, 64,107)= 5,965.9, p < 0.0001, partial η2 = 0.09, and sex, F(1, 64.107)= 2,229.5, p < 0.0001, partial η2 = 0.03. Curves are polynomial fits (black: both sexes; red: females; blue: males). Vertical lines represent standard errors of the mean (±SEM; in most cases they are smaller than the respective symbols). Inset: Distribution of BMI in the dataset used in this study (normalised for age and sex, BMI corr ). The standardised categories indicating normal, overweight and obese body weight are shown in green, pink and red, respectively. Note, that we divided the population into two groups (normal: BMI < 25; and overweight/obese: BMI 25) for analysing the relationship between social jetlag and BMI.
Figure S3. The Relationship between Sleep Duration and BMI on Workdays (A) and on Free Days (B) Was Assessed by Partial Correlations Controlling for Age and Sex The association between sleep duration and BMI was almost twofold higher on work days (panel A: r = - 0.079; p < 0.0001) than on free days (panel B: r = -0.046, p < 0.0001, n = 64,106). Vertical lines are ±SEM of the 30-min data bins. For illustration purposes, sleep duration and BMI were corrected for age and sex (SDcorr; BMIcorr) according to the fits shown in Figures 1 and S1, see Table S3 for the fit parameters.
Table S1A. Predicting BMI in the Lower BMI Group (BMI < 25, n = 43.308), Related to Figure 3 Step Predictor variables B SE β 95% CI Step 1 Constant 1.302 0.002 1.298 to 1.305 Age 0.066 0.001 0.244*** 0.063 to 0.068 Sex 0.018 0.000 0.262*** 0.017 to 0.018 Step 2 Constant 1.320 0.005 1.310 to 1.331 Age 0.065 0.001 0.242*** 0.063 to 0.067 Sex 0.018 0.000 0.260*** 0.017 to 0.018 Average sleep duration -0.017 0.004-0.017*** -0.025 to -0.008 Step 3 Constant 1.325 0.006 1.310 to 1.331 Age 0.064 0.001 0.239*** 0.062 to 0.067 Sex 0.018 0.000 0.260*** 0.017 to 0.018 Average sleep duration -0.017 0.004-0.018*** -0.026 to -0.009 Social Jetlag -0.003 0.002-0.007-0.007 to 0.001 Step 4 Constant 1.333 0.007 1.320 to 1.346 Age 0.063 0.001 0.236*** 0.061 to 0.066 Sex 0.018 0.000 0.262*** 0.017 to 0.018 Average sleep duration -0.019 0.004-0.020*** -0.028 to -0.010 Social Jetlag -0.001 0.002 0.004-0.003 to 0.006 Chronotype -0.010 0.003-0.020*** -0.015 to -0.004 This table summarises the 4-step hierarchical multiple regression model (enter method). B: unstandardised regression coefficients; SE: standard error; ß: standardised regression coefficients; CI: confidence intervals (lower to upper bound). r 2 = 0.127 for Step 1, r 2 change = 0.000 (F (1,43304) = 14.36), p = 0.000) for Step 2, r 2 change = 0.000 (F (1,43303) = 1.88, p = 0.17) for Step 3, and r 2 change = 0.000 (F (1,43302) = 12.05, p < 0.000) for Step 4.
Table S1B. Predicting BMI in the Higher BMI Group (BMI 25, n = 20.735), Related to Figure 3 Step Predictor variables B SE β 95% CI Step 1 Constant 1.492 0.004 1.48 to 1.50 Age 0.018 0.003 0.046*** 0.013 to 0.024 Sex -0.010 0.001-0.105*** -0.012 to -0.009 Step 2 Constant 1.585 0.011 1.63 to 1.608 Age 0.014 0.003 0.036*** 0.009 to 0.020 Sex -0.011 0.001-0.110*** -0.012 to -0.009 Average sleep duration -0.083 0.009-0.063*** -0.101 to -0.065 Step 3 Constant 1.556 0.013 1.53 to 1.58 Age 0.020 0.003 0.050*** 0.014 to 0.026 Sex -0.011 0.001-0.112*** -0.012 to -0.010 Average sleep duration -0.078 0.009-0.059*** -0.096 to -0.060 Social Jetlag 0.021 0.004-0.035*** 0.012 to 0.029 Step 4 Constant 1.558 0.014 1.531 to 1.586 Age 0.020 0.003 0.050*** 0.014 to 0.026 Sex -0.011 0.001-0.111*** -0.012 to -0.010 Average sleep duration -0.078 0.009-0.059*** -0.097 to -0.060 Social Jetlag 0.021 0.006-0.037*** 0.011 to 0.031 Chronotype -0.002 0.006-0.003-0.014 t0 0.010 We used a 4-step hierarchical multiple regression model (enter method). Abbreviations as in Table 2. r 2 = 0.012 for Step 1, r 2 change = 0.004 (F (1,20.731) = 79,49, p < 0.000) for Step 2, r 2 change = 0.001 (F (1,20.730) = 21,95, p < 0.000) for Step 3, and r 2 change = 0.000 (F (1,20.730) = 0,103, p < 0.000) for Step 4. Note: in case the regression model was run with a stepwise method, chronotype was deleted from the model.
Supplemental Experimental Procedures The Munich ChronoType Questionnaire (MCTQ) is the central instrument of an on-going internet-based study about circadian behaviour (see chronotype study @ www.thewep.org). Users are either the general internet-using public or individuals who are invited to participate in a dedicated project. The online MCTQ database grows on average by 825 entries/month (not including dedicated projects). The MCTQ asks simple questions about peoples sleep habits when they go to bed, when they prepare for sleep, how long they take to fall asleep, when they wake up, whether wake-up is assisted by an alarm clock and when they get out of bed. These questions are asked separately for workdays and work-free days. The first assessment of chronotype is based on the midpoint of sleep on free days (MSF), which is then corrected for oversleep on free days (MSF sc ), a phenomenon reflecting compensation for accumulated sleep debt during the workweek. The detailed evaluation of the MCTQ and the computation of all derived variables is shown at the end of the SOM. BMI was computed by using the standard formula: weight (kg)/height (m) 2. We also asked how much time people spend outdoors without a roof above their head. This study is based on the database as of the end of July 2010. It contained 81,888 entries of the general internet-using public (i.e., dedicated projects were not used for the evaluations presented here). This dataset was inspected for and cleaned (see table S2) of extreme values, for individuals who have no regular week schedules (i.e., not working or having no free days during the week) as well as for those who use alarm clocks on free days (since chronotype cannot be calculated for individuals who never report sleep times without restrictions). After cleaning, the number of entries for each year are as shown in table S3. Except for the trend analyses shown in Fig. 4, the year 2002 was excluded from analyses due to the strong imbalance of entries compared to other years. Age- and sex corrections were performed by calculating the difference (in %) between the actual value and the value yielded by the sex-specific fit for the given age. This difference was then added to the value of the fit for the total population for the age of 30 years (see respective fit parameters in table S4). Table S2. Exclusion Criteria Applied to the Database and the Remaining Number of Participants in the Sample, Related to Supplemental Experimental Procedures n.i. = not indicated. Exclusion criteria Number of Exclusions Age < 16 and > 65 or n.i. 2,746 BMI <16 and > 60 or n.i. 1,813 Sleep offset on workdays > 12:00 or n.i. 1,763 Number of work days 0 or 7 2,681 Alarm on free days 8,661 MSF > 12:00 15 Average sleep duration < 3h or > 13h 6 Sleep duration on free days < 3h 18 Sleep duration on work days < 3h or > 13h 55 Age and sex corrected sleep duration on free days > 13h 20 Final sample size 64,110
Table S3. Number of Entries Per Year in the Final (Cleaned) Database, Related to Figure 3 2002 673 2003 13,357 2004 6,067 2005 6,207 2006 15,395 2007 4,270 2008 9,570 2009 4,876 2010 3,695
Table S4. Parameters of the Polynomial Fits (y = a + bx + cx 2 + dx 3 + ex 4 + fx 5 + gx 6 ), Related to Figure 2 cat a b c d e f g R 30 y SJL M 17 26.071-6.3573 0.51212-0.012828 0.976 M>17 11.58-1.0743 0.05204-0.0014569 2.39E-05-2.13E-07 7.85E-10 0.997 F 17-22.384 5.1481-0.3706 0.0091667 0.982 F>17 13.989-1.5627 0.084105-0.002481 4.16E-05-3.72E-07 1.36E-09 0.993 TP 17 6.8042-1.7169 0.15321-0.0038889 0.00E+00 0.00E+00 0.00E+00 0.989 TP>17 10.451-0.92698 0.042422-1.12E-03 1.76E-05-1.55E-07 5.74E-10 0.996 1.56 SD w M 19 1.5752 2.3268-0.21291 0.005564 0.990 M>19 6.0106 0.09936-0.002421 1.47E-05-5.73E-07 1.89E-08-1.42E-10 0.959 F 19 7.7934 1.0746-0.13395 0.004032 0.978 F>19-0.18031 1.1436-0.066409 0.0019379-3.06E-05 2.50E-07-8.24E-10 0.958 TP 19 5.1134 1.6146-0.16808 0.004697 0.984 TP>19 0.88507 0.93715-0.053674 0.0015756-2.57E-05 2.22E-07-7.91E-10 0.978 7.24 SD f M 6.3124 0.72941-0.058922 0.0020518-3.64E-05 3.23E-07-1.14E-09 0.977 F 7.0395 0.71279-0.063496 0.002372-4.43E-05 4.06E-07-1.45E-09 0.987 TP 7.3667 0.61357-0.054727 0.0020192-3.73E-05 3.41E-07-1.22E-09 0.992 8.18 ØSD M 19 6.2888 0.99604-0.096265 0.0024864 0.977 M>19 13.023-0.86933 0.056094-0.0018537 3.22E-05-2.79E-07 9.46E-10 0.962 F 19 4.885 1.4568-0.1393 0.0037374 0.967 F>19 3.9514 0.6445-0.03985 0.0011973-1.93E-05 1.61E-07-5.41E-10 0.965 TP 19 7.1109 0.87745-0.091677 0.0024747 0.973 TP>19 8.7488-0.13886 0.0091939-0.0003475 6.54E-06-5.78E-08 1.93E-10 0.976 7.50 BMI M 8.5475 1.1942-0.03188 0.00040281-1.98E-06 0.974 F 13.829 0.59412-0.01233 0.00012135-4.72E-07 0.923 TP 11.148 0.91046-0.023483 2.95E-04-1.45E-06 0.971 Abbreviations: cat = category (e.g., M 17, males younger than 18); M = males; F = females; TP = total population; SJL = social jetlag; SDw = sleep duration on workdays; SDf = sleep duration on free days; ØSD = average sleep duration over the work week; BMI = body mass index; R = regression coefficient of the fit; 30y = fit result for the age of 30 in the total population (serving as the basis for the normalisation).
Evaluation of the MCTQ The basic evaluation of the MCTQ comprises many different, interrelated variables. In the MCTQshift, these variables are calculated separately for each shift (e.g., morning, evening and night shift). Asterisks indicate variables that are computed from the direct answers in the MCTQ. 1. Basic Variables BTw Local time of going to bed on work days; SPrepw Local time of preparing to sleep on workdays; SLatw Sleep latency on workdays I need min to fall asleep ; * SOw Sleep onset on workdays = SPrepw + SLatw; SEw Sleep end on workdays; SIw Sleep inertia on workdays after min, I get up * GUw Local time of getting out of bed on workdays = SEw + SIw; * SDw Sleep duration on workdays = SEw SOw; * TBTw total time in bed on workdays = GUw BTw; BTf Local time of going to bed on free days; SPrepf Local time of preparing to sleep on free days; SLatf Sleep latency on free days I need min to fall asleep ; * SOf Sleep onset on free days = SPrepf + SLatf; SEf Sleep end on free days; SIf Sleep inertia on free days after min, I get up * GUf Local time of getting out of bed on free days; * SDf Sleep duration on free days = SEf SOf; * TBTf total time in bed on free days = GUf BTf; WD = number of workdays; * SDweek (Average SD across the week) = (SDw x WD + SDf x (7-WD))/7; * SLOSSweek (Sleep loss across the week): if SDweek > SDw: SLOSSweek = (SDweek SDw) x WD; if SDweek SDw: SLOSSweek = (SDweek SDf) x (7 WD); 2. Chronotype The basis for estimating chronotype is the Mid-Sleep Time on Free days (MSF) MSF = SOf + (SDf)/2; MSF is then corrected for oversleep on free days that subjects use to compensates the sleep debt accumulated during the workweek: MSFsc = MSF - (SDf - SDw)/2; This correction is only applied to individuals who sleep longer on free days than on workdays. For all others: MSFsc = MSF; MSFsc is the basic assessment for chronotype for an individual under the current circumstances; (MSFsc depends on developmental and environmental conditions, e.g., age and light exposure). For epidemiological and genetic studies, MSFsc is normalised for age and sex to make populations of different age and sex compositions comparable. 3. Social Jetlag The relative social jetlag (SJLrel) is the difference between the Mid-Sleep on work- and on free days: MSW = SOw + (SDw)/2; SJLrel = MSF MSW; The absolute social jetlag (SJL) is used for most assessments of the consequences of social jetlag: SJL = abs(sjlrel);