Datasets description

Analysis.

A. Analysis flowchart. IRASA – Irregularly Resampled Auto-Spectral Analysis, sgolayfilt – Savitzky-Golay filter. B. Outputs of some of the analysis steps in an example healthy 26-year-old individual. From top to bottom: time-frequency representation of the total spectral power, raw and smoothed time series of the fractal slopes and hypnogram. Frontal spectral power and its slopes were calculated in the 0.3 – 30 Hz range for each 30 seconds of sleep.

Fractal cycles in healthy adults. A – B. Individual fractal and classical sleep cycles.

Time series of smoothed z-normalized fractal slopes (bottom) and corresponding hypnograms (top) observed in two participants. Time series of the fractal slopes and corresponding hypnograms for all participants are reported in Supplementary Material 1. The duration of the fractal cycle is a time interval between two successive peaks (blue diamonds). A: S15 from Dataset 3 shows a one-to-one match between fractal cycles defined by the algorithm and classical (non-REM – REM) cycles defined by the hypnogram. B: In S22 from dataset 5, the second part of night has many wake epochs, some of them are identified by the algorithm as local peaks. This results in a higher number of fractal cycles as compared to the classical ones and a poor match between the fractal cycles No. 3 – 7 and classical cycles No. 2 – 5. The algorithm does not distinguish between the wake and REM-related fractal slopes and can define both as local peaks. Since the duration of the fractal cycles is defined as an interval of time between two adjacent peaks, more awakenings/arousals during sleep (usually associated with aging, Fig.S5 B) are expected to result in more peaks and, consequently, more fractal cycles, i.e., a shorter cycle duration. This is one of the possible explanations for the correlation between the fractal cycle duration and age (shown in Fig. S5 A). C. Scatterplots: each dot represents the duration of the cycles averaged over one participant. The durations of the fractal and classical sleep cycles averaged over each participant correlate in all analyzed datasets, raw (non-ranked) values are shown. D. Cycle-to-cycle overnight dynamics show an inverted U shape of the duration of both fractal and classical cycles across a night and a gradual decrease in absolute amplitudes of the fractal descents and ascents from early to late cycles. SWS – slow-wave sleep, REM – rapid eye movement, r – Spearman’s correlation coefficient.

Demographic, sleep and fractal characteristics

Fractal cycles in children and adolescents.

A – B: Individual cycles: time series of smoothed z-normalized fractal slopes (bottom) and corresponding hypnograms (top). The duration of the fractal cycle is a time interval between two successive peaks (blue diamonds). A: In this 9.9-year-old participant, we split the first 150-minute-long classical cycle into two cycles according to the definitions of a “skipped” cycle presented in Methods. The fractal cycle algorithm successfully detected this skipped cycle. B: This 14.9-year-old participant has a 156-minute-long first classical cycle. Visual inspection shows that it should be divided into 3 skipped cycles, however, our a priori definition of skipped cycles did not include an option to subdivide a long cycle into three short cycles; hence, we split it into two short cycles. The fractal cycle algorithm was sensitive to these sleep lightenings and detected all three short cycles. Classical cycle 4 looks like a skipped cycle as it has two clear episodes of slow-wave sleep separated by non-REM stage 2. However, the length of this cycle is shorter than 110 min (the threshold defined a priori), therefore, we did not split the classical cycle 4 into two cycles. The fractal cycle algorithm was sensitive to this lightening of sleep and defined two fractal cycles during this period. C. Histograms: The frequency distribution of fractal (left) and classical (right) cycle durations in children and adolescents (mean age: 12.4 ± 3.1 years) compared to young adults (mean age: 24.8 ± 0.9 years). Kolmogorov-Smirnov’s test rejected the assumption that cycle duration comes from a standard normal distribution. D. Box plots: in each box, a vertical central line represents the median, the left and right edges of the box indicate the 25th and 75th percentiles, respectively, the whiskers extend to the most extreme data points not considered outliers, and a plus sign represents outliers. Children and adolescents show shorter fractal cycle duration compared to young adults. E. Overnight dynamics: cycle-to-cycle dynamics show that only the mean durations of the first two fractal cycles are significantly shorter in the pediatric compared to control group, * marks a statistically significant difference between the groups, SWS – slow-wave sleep, REM – rapid eye movement sleep.

Fractal cycles in MDD.

A. Individual fractal cycles: time series of smoothed z-normalized fractal slopes observed in a 22 y.o. MDD patient (Dataset B) in their unmedicated (top) and 7-day medicated (bottom) states. Fractal cycle durations (defined as an interval of time between two successive peaks) are longer in the medicated compared to unmedicated states, reflecting shallower fluctuations of fractal (aperiodic) activity. Two additional patients are shown in Fig.S9 (Supplemental Material 2). B. Box plots: the fractal cycle duration is longer in medicated MDD patients (red) compared to age and gender-matched healthy controls (black) in all datasets. In Dataset B, fractal cycles are longer in the medicated vs. patients’ own unmedicated state and in patients who took REM-suppressive vs. REM-non-suppressive antidepressants. A vertical central line represents the median in each box, the left and right edges of the box indicate the 25th and 75th percentiles, respectively, the whiskers extend to the most extreme data points not considered outliers, and a plus sign represents outliers (individual cycles). C. Frequency distribution: individual fractal and classical cycles pooled from three MDD datasets (A – C) are counted separately for medicated MDD patients and HC. D. Overnight dynamics: cycle-to-cycle dynamics of the duration of both fractal and classical cycles show a gradual decrease in medicated patients vs. an inverted U shape in controls. The between-group difference in cycle duration is the largest for the first cycle. Patients show flatter fractal descents of the second cycle and steeper fractal descents of the fourth cycle compared to controls. Contrary to controls, patients do not show a gradual decrease in absolute amplitudes of the fractal descents from the second to the fourth cycles. Patients and controls show comparable cycle-to-cycle dynamics of fractal ascents, * marks a statistically significant difference between the groups, MDD – major depressive disorder, HC – healthy controls, unmed. – unmedicated, med. – medicated, SWS – slow-wave sleep, REM – rapid eye movement sleep.

Fractal cycles in MDD

Hypothetical functional significance of fractal cycles