Cortical signatures of sleep are altered following effective deep brain stimulation for depression

Patient and trial details

Data reported here were collected as part of an experimental trial to derive biomarkers to optimize SCC DBS for TRD (ClinicalTrials.gov: NCT01984710). Ten subjects with treatment-resistant major depressive disorder were consecutively enrolled. All participants provided written informed consent, and the protocol adhered to the Declaration of Helsinki. The protocol was approved by the Institutional Review Boards of Emory University, Icahn School of Medicine at Mount Sinai and Georgia Institute of Technology, as well as the US Food and Drug Administration under a physician-sponsored Investigational Device Exemption (IDE G130107), with monitoring by the Emory University Department of Psychiatry and Behavioral Sciences Data and Safety Monitoring Board.

The study design included the following phases: 1) four weeks of presurgical baseline, during which only clinical assessments were performed, 2) at least four weeks of recovery following surgical implantation of the device, during which LFP was recorded daily with stimulation turned OFF, 3) twenty-four weeks of therapeutic DBS, and 4) one week when stimulation was discontinued in a single-blind manner (i.e. patient unaware) with daily LFP recordings (Figure 1A). All patients showed improvement on the Hamilton Depression Rating Scale (HDRS-17) following treatment (Table 1), though some improvement was already observed in the Pre-DBS post-operative period, consistent with placebo or lesion effects [53].

• Download figure
• Open in new tab

Figure 1: Quantifying slow-wave activity in bilateral SCC LFP segments collected at different times of day.

A: Overview of the data collection timeline. LFP segments from weeks 3 and 4 after surgery, but before initiation of therapeutic DBS comprised the Pre-DBS baseline phase (mean number of Pre-DBS samples 121±59; typical sampling frequency 1 sample / 3.06 hours), while LFP segments recorded 1 week after 24 weeks of therapeutic DBS (during which DBS was switched off) were used for the Post-DBS phase (mean number of Post-DBS samples 386±74; typical sampling frequency 1 sample / 0.375 hours). B, C: Left (B) and right (C) SCC LFP segments (low-pass filtered with a cut-off at 30Hz) collected at different times of a single day from one example patient during the Pre-DBS phase. Scale bar represents 100uV. D: SWA (0.5-4.5Hz) power z-scores corresponding to the left (black) and right (grey) hemisphere LFP segments shown in B&C. E: Histogram of log₁₀ SWA power (averaged between left and right) for all Pre-DBS LFP segments for one example patient. The red dotted line represents the 80^th percentile that was used as a cut-off for ‘high SWA’ putative sleep episodes. F: Average power spectrum (Welch’s method) of ‘high SWA’ (>80^th percentile) LFP segments compared to ‘low SWA’ (<50^th percentile) LFP segments. G: Scatter plot of the correlation between left and right hemisphere SWA power. Red data points represent ‘high SWA’ LFPs. H: Correlation coefficients of left and right SWA power for all 6 patients for which bilateral data was obtained, for the Pre-DBS and Post-DBS phase (Individual data points in red, mean and standard deviation provided in black; Pre-DBS: r = 0.805±0.0889, Post-DBS: r = 0.782±0.105, signed-rank = 17, p = 0.219, n = 6).

Clinical scoring

Clinical symptom severity was evaluated by an independent rater using the 17-item Hamilton Depression Rating Scale (HDRS-17), Montgomery-Åsberg Depression Rating Scale (MADRS), and self-reported Beck Depression Inventory (BDI) during weekly visits to the laboratory. In addition to DBS, participants patients were prescribed with different medications depending on their clinical history, which included antidepressants and benzodiazepines (Supplementary Table 1). The dosage of these medications was held constant throughout the study.

Surgical implantation, device details and DBS parameters

Participants underwent bilateral implantation of DBS electrodes (3387, Medtronic, Minneapolis, MN) at the intersection of four white matter pathways (forceps minor, cingulum bundle, uncinate fasciculus, and frontostriatal fibers) as described previously [52,54]. Stimulation was delivered using a voltage-controlled implanted pulse generator (Activa PC+S, Medtronic Inc., Minneapolis, MN) that allowed the collection of local field potential from contacts immediately adjacent of the stimulation contact on either side in a differential configuration. DBS therapy started four weeks after the implantation surgery to allow stable recordings of the depressed baseline, as per our past published protocols [3,4]. Therapy consisted of bilateral monopolar stimulation on a single contact per hemisphere at 130 Hz with 90 µs pulse width. Stimulation amplitude, initialized at 3.5 V, was varied by study psychiatrists based on clinical judgement throughout the 24 weeks. Final stimulation dose ranged from 3.5 – 5 V; without change in contact, pulse width or frequency.

Acquisition and analysis of LFP data

LFPs were acquired in 12 s segments by the Activa PC+S at 422 Hz at regular intervals. In the post-operative phase, ‘Pre-DBS’ LFPs were acquired approximately every 2-4 hours. Data were downloaded from the device weekly. In the Post-DBS phase, LFPs were collected approximately every 20-30 minutes with daily downloads; the precise sampling regime for individual patients was determined pragmatically by the data storage capacity of the device and the interval between successive clinical visits (For each patient’s exact sampling times, see Supplementary Figures S2-S10, panels A&B). We report mean typical sampling frequencies defined as follows: For each patient, we took the median across days of the median sampling interval for each individual day; typical sampling frequency for the study phase was then calculated as the mean across patients. The Pre-DBS typical sampling frequency was 1 LFP / 3.06 hours, while the Post-DBS typical sampling frequency was 1 LFP / 0.375 hours. During clinical visits, the LFP data were downloaded from the device and saved in text format. LFP data were compiled across days and stored in shareable formats (HDF-5, .mat).

We obtained SCC LFP data from 18 hemispheres (9 patients). Data from one patient was contaminated by amplifier artifact bilaterally (P005). In P001 and P003, data from only one hemisphere was of adequate quality as the other hemispheres had artifacts driven by electro-cardiographic activity. The remaining dataset represented 14 hemispheres (8 patients; See Table 1 for the included hemispheres for each participant). Baseline Pre-DBS recordings were collected between 14 and 28 days post-operatively, prior to DBS therapy (Figure 1A); for one patient (P008) data from only 4 days at the start of this period were available due to operational error. Post-DBS data were obtained after 24 weeks of therapeutic DBS (Figure 1A).

The first 2s of each LFP segment were rejected to remove the amplifier switching artifact. LFP segments with a voltage range that exceeded a Z-score of 6 (calculated across each individual patient’s LFP segments for the Pre-DBS and Post-DBS periods separately) were excluded from further analysis. This resulted in a mean of 121±59 useable LFP segments (10s each) per hemisphere in the Pre-DBS phase and a mean of 386±74 useable LFP segments in the Post-DBS phase.

Calculation of spectral power

SWA power was calculated as band power in the 0.5-4.5Hz range for each LFP segment and each hemisphere using the MATLAB bandpower() function. Spindle band power (for Figure 3) was similarly calculated as power in the 9-16Hz range. To determine the peak frequency of spindle-like activity during LFPs collected during putative NREM sleep, we used continuous wavelet transform (MATLAB’s cwt()) to calculate a global wavelet spectrum by taking the square of the absolute value of the wavelet transform and taking the mean across the time dimension.

Stability of SWA

The stability of SWA was quantified using the Temporal Variability Index (TVI), a measure of oscillatory stability independent of overall signal power [55]. Briefly, we took the Hilbert envelope of the signal filtered between 0.5 and 4.5 Hz, and divided this by its mean value to remove absolute amplitude differences. We next took the derivative of this signal – a timeseries with normally distributed values centred on zero. The TVI was defined as the standard deviation of these timeseries values.

Quantifying the diurnal pattern of SWA

To quantify the time towards which SWA was biased at a patient level, we averaged the SWA power of both hemispheres, unless one of the SCCs was excluded from analysis (in which case the SWA estimate was based solely on the included hemisphere). SWA power estimates were log-transformed (base 10) after adding a constant of 1 to avoid zero or negative values, and then normalised to their median. This transform reduced the right-tailedness of the distribution of SWA power estimates and made it more amenable to fitting and visualisation.

To generate a fit to estimate SWA density across the 24h of the diurnal cycle, we converted the time of day associated with each of these SWA power samples to a numeric x value between 0 (00:00) and 24 (24:00). X values and associated SWA power estimates were duplicated with an offset of −24 to enable smooth fitting around midnight (00:00, x =0). This allowed us to perform a smoothing spline fit to generate an SWA estimate for each time of day. The time-of-day fit was generated using MATLAB’s fit() function using a smoothing spline fit type with a smoothing parameter of 0.9. To quantify the size and the timing of the most concentrated night-time SWA activity, we found the time and the value of the largest peak in the estimate between 18:00 in the evening and 10:00 in the morning (the fit segment for −6 < x < 10).

Resampling control

LFPs were more sparsely sampled in the Pre-DBS phase (typical sampling frequency 1 LFP / 3.06 hours; Supplementary Figures S2-S10, panel A) compared to the Post-DBS phase (typical sampling frequency 1 LFP / 0.375 hours; Supplementary Figures S2-S10, panel B), which might have had an effect on the SWA power time-of-day fit and our estimates of the night-time SWA peak. We therefore investigated whether resampling the Post-DBS phase according to the Pre-DBS sample numbers and times of day would be expected to substantially alter our estimates of SWA timing. For 1000 iterations, we generated a new set of Post-DBS power estimates in the following way: For each unique time-of-day at which any Pre-DBS samples were taken, we determined the number of Pre-DBS samples taken at that exact time of day, and then randomly selected the same number of Post-DBS LFPs that were collected within half an hour on either side of this clock time. We computed an SWA time-of-day fit as described above for each of these resampled data sets, and determined its night-time peak and peak time. This resulted in a distribution of 1000 fits, 1000 peak values and 1000 peak times. For visual comparison, we compared the mean of these 1000 fits to the Pre-DBS fit (Figure S1G). For the comparison of SWA density and timing, we compared the median peak height and median peak time from the Post-DBS period with the peak height and time from the original Pre-DBS fit (Figure S1H,I). The resampled Post-DBS SWA fits for each patient and the resulting distribution of peak estimates are shown in Supplementary Figures S2-S10, panels G&H.

Spindle detection and analysis

Spindle detection was carried out using a simplified, frequency-specific version of the spindle detection approach described by Purcell and colleagues [50]. For each hemisphere, high-SWA LFPs (>80^th percentile of SWA power) during the night (21:00-09:00) were selected as putative night-time NREM sleep episodes. For each high-SWA LFP, a global wavelet spectrum was calculated as follows: a continuous wavelet transform was carried out using MATLAB’s cwt() function with a Morse wavelet. The resulting complex transformed data was squared and averaged across the time dimension to result in a wavelet-based estimate of power at each frequency. The resulting vector was multiplied by frequency squared, which brought out a clear peak in the spindle frequency range for most data sets. The frequency corresponding to the peak in the average vector across patients and hemispheres (11.66Hz) was then used to set the centre frequency of a 7-cycle complex Morlet wavelet.

To detect spindles, each high-SWA LFP segment was convolved with this complex Morlet wavelet. Spindle power was defined as the absolute value of the complex result. We used the median absolute deviation (MAD) as a measure of variability of this signal. Any peak in the spindle power signal that exceeded median power + 5 * MAD was considered a putative spindle. For each of these putative spindles, onset was calculated as the last crossing of the lower threshold value, median spindle power + 1 * MAD, before the peak. Similarly, the spindle offset was the first time the spindle power crossed this lower threshold in the negative direction. Onset and offset values were then used to calculate spindle duration – any putative spindle with a duration under 0.5s was rejected. For each spindle, amplitude was calculated as the maximum spindle power value attained between spindle onset and offset. Spindle density was calculated as the number of spindles observed divided by the cumulative length of time of all eligible (i.e. high-SWA) LFP segments, and expressed as number of spindles per minute (assessed across all relevant 10-second LFP segments).

Statistical analysis

Data were analysed in MATLAB R2021B (The Mathworks, Natick, MA, USA). Summary data are reported as mean ± standard deviation. For comparison of two means, both samples were tested for normality using the Shapiro-Wilk test (using A. BenSaida’s MATLAB implementation swtest(), revision 3.0, 2014). If no samples significantly (α < 0.05) deviated from a normal distribution, means were compared using two-tailed t-tests (paired or unpaired as appropriate). If normality was violated, medians were compared using the Wilcoxon signed-rank or Wilcoxon rank-sum tests.

Slow-wave activity can be quantified using on-device LFP recordings from an SCC DBS implant

10-second LFP segments were obtained at regular intervals from the SCC of 9 patients implanted with bilateral DBS leads before and after a 6-month period of continuous high frequency DBS (Figure 1A-C). Slow wave activity (SWA, 0.5-4.5Hz) power was calculated for each LFP segment (Figure 1D). LFPs with SWA power above the 80^th percentile were classed as ‘high SWA power’ episodes that might reflect NREM sleep (Figure 1E). In these LFPs, the increase in low-frequency LFP power was specific to the lower frequencies and in particular the delta band (e.g. Figure 1F), i.e. it did not reflect a broad-band or otherwise non-specific increase in overall signal power.

While sleep has local components on shorter time scales, it is a global brain state with high levels of brain-wide synchronization, and SWA power is therefore expected to correlate strongly between 10-second LFP samples across hemispheres. To confirm this, we quantified the correlation of left and right SWA power (Figure 1G; individual participant correlations provided in Figures S2-10C,D). In the patients where both hemispheres met inclusion criteria (n = 6), correlation of left and right log SWA power (Pearson’s) was r = 0.80±0.09 (Pre-DBS) and r = 0.78±0.11 (Post-DBS; Figure 1H). The difference between the Pre-DBS and Post-DBS phase was not significant (signed-rank = 17, p = 0.219, n = 6). SWA power correlations between hemispheres for all patients are visualised as scatter plots in Supplementary Figures S2-S10, C&D.

Slow-wave activity is more temporally focused and occurs earlier in the night after successful SCC DBS for depression

We characterised SWA power across the 24 hours of the day during the Pre- and Post-DBS periods using a smoothing spline fit to SWA power in the LFP segments corresponding to each time of day (Figure 2A-F). This revealed that, across patients, there was an earlier and sharper increase in SWA power at night in the Post-DBS period (Figure 2G). We quantified the size as well as the timing of the greatest night-time (18:00-10:00) peak of the SWA fit, and found that across patients, the SWA fit peak was significantly greater in the Post-DBS phase (Figure 2H; Pre-DBS: 2.00±0.53, Post-DBS: 2.41±0.40, t(7) = −3.81, p = 0.0066, n = 8) and occurred significantly earlier in the night (Figure 2I; 03:34±02:08, Post-DBS: 01:19±01:33, t(7) = 3.26, p = 0.0139, n = 8). The Pre-DBS and Post-DBS SWA data points and fits are provided for all patients in Supplementary Figures S2-S10, E&F. As LFPs were less frequently and less densely sampled in the Pre-DBS phase compared to the Post-DBS phase, we replicated this analysis using a random resampling approach during which Post-DBS was only sampled according to the sampling clock times and frequency from the Pre-DBS period. We generated 1000 resampled data sets per patient, and computed 1000 SWA power time-of-day fits with corresponding fit peaks (Figure S1A-G). Comparing the median of each patient’s Post-DBS fit peak distribution with the Pre-DBS fit peak still yielded a significant difference in both SWA peak size and timing (Figure S1 H,I; Peak size: Pre-DBS: 2.00±0.53, Post-DBS: 2.19±0.50, t(7) = −2.86, p = 0.024, n = 8; Peak time in 24h clock time: Pre-DBS: 03:34±02:08, Post-DBS: 01:49±01:13, t(7) = 2.69, p = 0.031, n = 8). Each patient’s distribution of Post-DBS fits and fit peak estimates is visualised in Supplementary Figures S2-S10 G&H.

• Download figure
• Open in new tab

Figure 2: Slow-wave activity is more temporally consistent and earlier in the night after SCC DBS.

A, B: Log-transformed and median-normalized SWA power (mean across both hemispheres; black) and a smoothing spline fit (blue) around the 24 hours of the day in the Pre-DBS phase for two example patients. C: Mean (solid line) plus and minus standard deviation (dotted lines) across patients of Pre-DBS time-of-day fits as illustrated in A and B. D, E: Log-transformed and median-normalized SWA power (mean across both hemispheres; black) and a smoothing spline fit (orange) around the 24h of the day in the Post-DBS phase for the same example patients shown in A and B. F: Mean (solid line) plus and minus standard deviation (dotted lines) across patients of Post-DBS time-of-day fits as illustrated in D and E. G: Mean Pre-DBS (blue line) and Post-DBS (orange line) time-of-day fits across patients highlighting increased SWA in the late evening and early night in the Post-DBS phase. H: Pre-DBS vs. Post-DBS maximum value of the SWA time-of-day fit (normalised to median fit value; Pre-DBS: 2.00±0.53, Post-DBS: 2.41±0.40, t(7) = −3.81, p = 0.0066, n = 8). Gray lines represent female patients, black lines represent male patients, red dots indicate differences. I: Pre-DBS vs. Post-DBS peak time (i.e. the clock time associated with the maximum value of the SWA time-of-day fit; Pre-DBS: 03:34±02:08, Post-DBS: 01:19±01:33, t(7) = 3.26, p = 0.0139, n = 8). Gray lines represent female patients, black lines represent male patients, red dots indicate differences.

• Download figure
• Open in new tab

Figure 3: Power and stability of night-time SCC slow wave activity Pre- and Post-DBS.

A: SCC night-time sleep LFP segment from an example patient during the Pre-DBS phase (low-pass filtered with a cutoff of 30Hz). B: SCC night-time sleep LFP segment from the same example patient during the Post-DBS phase (low-pass filtered with a cutoff of 30Hz). C: Mean night-time LFP power spectra across patients Pre-DBS (blue) and Post-DBS (orange). D: Mean (black line) and standard deviation (grey shaded region) night-time sleep LFP power spectral differences (multiplied by frequency for visualization) between the Pre-DBS and Post-DBS phases. E: Relative SWA power of night-time SCC LFPs before and after DBS (Pre-DBS: 12.9±2.49%, Post-DBS: 12.3±1.88%, t(7) = 0.895, p = 0.4, n = 8). Gray lines represent female patients, black lines represent male patients, red dots indicate differences. F: Temporal variability index (TVI) of night-time SCC LFPs before and after DBS (Pre-DBS: 6.8±0.489, Post-DBS: 6.98±0.4, t(7) = −1.75, p = 0.123, n = 8). Gray lines represent female patients, black lines represent male patients, red dots indicate differences. G: Relative spindle power of night time SCC LFPs before and after DBS (Pre-DBS: 1.79±0.415, Post-DBS: 2.15±0.837, t(7) = −1.89, p = 0.1, n = 8). Gray lines represent female patients, black lines represent male patients, red dots indicate differences.

We next compared these SWA sleep metrics to patients self-reported sleep ratings from the HDRS. Interestingly, patients already scored relatively low on the individual sleep-related items in the Pre-DBS phase, and as such there was no clear difference in scores on these items Post-DBS (see Table 2). In addition, initial SWA peak height and peak timing did not significantly correlate with Pre-DBS Hamilton score or sleep sub-scores, and changes in SWA peak height and peak timing did not significantly correlate with changes in the Hamilton overall score or sleep sub-scores (Supplementary Tables 2&3).

Power and stability of slow-wave activity are not altered after SCC DBS for depression

We next investigated whether the nature of SWA during putative sleep episodes was different. To exclude circadian and light cycle effects, only night-time LFPs with relatively high SWA amplitudes (LFPs < 80^th percentile of SWA power collected after 9pm and before 9am) were included in this analysis. Figure 3A&B show representative segments from the Pre-DBS and Post-DBS phases. We compared spectral power in these LFPs, which showed a relatively broad-band increase in the low frequencies between the Pre-DBS and Post-DBS phases (Figure 3C,D). To exclude broad-band effects resulting from any changes to the electrode-brain interface during the 6 months of DBS, we expressed power as a percentage of total signal power determined from the power spectrum. This revealed no difference in relative SWA power after SCC DBS (Figure 3e; Pre-DBS: 12.9±2.49%, Post-DBS: 12.3±1.88%, t(7) = 0.895, p = 0.4, n = 8 patients), though there was an increase in absolute power (Pre-DBS: 179±82µV²/Hz, Post-DBS: 227±111µV²/Hz, t(7) = −2.92, p = 0.0223, n = 8 patients – data not shown). We quantified the stability of slow-wave activity using a temporal variability index (TVI) that quantifies fluctuations in power envelope but is insensitive to overall changes in power [55]; see Supplementary Methods). The TVI showed no consistent difference in SWA stability between Pre-DBS and Post-DBS phases (Figure 3F; Pre-DBS: 6.8±0.489, Post-DBS: 6.98±0.4, t(7) = −1.75, p = 0.123, n = 8 patients). Finally, we compared relative power in the sleep spindle frequency range (9-16Hz), which revealed no significant increase in relative spindle power after DBS (Figure 3G; Pre-DBS: 1.79±0.415, Post-DBS: 2.15±0.837, t(7) = - 1.89, p = 0.1, n = 8 patients). However, an increase in spindle power was present in absolute terms (Pre-DBS: 25.5±16.4 µV²/Hz, Post-DBS: 37.4±22.5µV²/Hz, signed-rank = 0, p = 0.00781, n = 8 patients– data not shown).

Sleep spindle activity is increased after successful SCC DBS treatment for depression

To investigate treatment effects on sleep spindles in more detail, we used a wavelet-based spindle detection approach (see Methods). Figure 4A shows a single night-time LFP with high SWA power from a single patient during the Post-DBS phase. Night-time SWA LFPs for this patient Post-DBS showed a peak in the frequency-adjusted wavelet spectrum in the spindle range (Figure 4B). A clear peak was present for most patients, with a mean spindle peak value of 11.66Hz, consistent with frequencies observed in cingulate cortex from previous human intracranial recordings [51]. This mean frequency was used to set the frequency of a complex 7-cycle Morlet wavelet for a wavelet transform to detect spindle events (see Methods), which successfully identified events consistent with sleep spindles (Figure 4C,D). The timing of detected spindles showed strong overlap with SWA power (e.g. Figure 4E,F), and for most patients, there was a strong correlation between times of high SWA power and times of high spindle activity (Figure 4G; Pre DBS: r = 0.414±0.385, Post DBS: r = 0.662±0.165). Moreover, we observed a shift in spindle activity between the Pre-DBS and Post-DBS phase (Figure 4H), analogous to the shift in SWA activity illustrated in Figure 2. Figure 4I shows the correlation between SWA power and percentage of spindle activity across all LFP snippets from all hemispheres (r = 0.43, p = 4.1*10^-241). There was an increased density of spindles in the Post-DBS phase (Figure 4J; Pre DBS: 5.51±3.82 spindles/min, Post DBS: 7.26±4 spindles/min, t(7) = - 3.77, p = 0.00702, n = 8). A comparison of the properties of individual spindles revealed no change in spindle duration (Figure 4K; Pre DBS: 0.773±0.165s, Post DBS: 0.78±0.0639s, t(7) = −0.13, p = 0.901, n = 8) or amplitude (Figure 4L, Pre DBS: 4.25±0.587, Post DBS: 4.51±0.542, t(7) = −2.17, p = 0.0668, n = 8). Changes in spindle density did not significantly correlate with changes in the Hamilton overall score or sleep items (Supplementary Table 4).

• Download figure
• Open in new tab

Figure 4: Sleep spindle activity is increased in the post-DBS phase.

A: A night-time sleep LFP segment from the right SCC of an example patient, collected during the Post-DBS period (low-pass filtered with a cutoff frequency of 30Hz). B: The mean (black line) and standard deviation (grey area) of the global wavelet spectra of all Post-DBS night-time sleep LFP segments for this example hemisphere, multiplied by frequency² to bring out the spindle peak. The point and dotted line indicate the spindle peak for this particular patient. C: All spindle events detected in night-time sleep LFP segments from the Post-DBS phase for this hemisphere, aligned to the spindle peak (transparent orange lines) and their mean (black line). D: The same LFP segment as in A, narrow-band filtered around the spindle peak frequency (11-12 Hz pass band), with orange sections indicating detected spindle events. E: Median night-time SWA power for this example hemisphere at different times of day for this example hemisphere. F: Percentage of spindle activity across LFP segments collected at different times of day for this example hemisphere. Note the close association with SWA power in E. G: Pearson’s correlation coefficients of SWA power and spindle percentage for the Pre-DBS and Post-DBS periods Pre DBS: r = 0.414±0.385, Post DBS: r = 0.662±0.165. H: Mean percentage of spindle activity across patients in the Pre-DBS (blue) and Post-DBS (orange) period across the 24h diurnal cycle. Note the shift towards earlier spindle times in the Post-DBS period. I: Scatter plot of log10 slow wave power vs. spindle percentage for all individual LFP snippets in the Post-DBS data set across participants, with regression line (purple) showing a clear positive correlation. J: Paired line plot representing spindle density before and after SCC DBS (Pre DBS: 5.51±3.82 spindles/min, Post DBS: 7.26±4 spindles/min, t(7) = −3.77, p = 0.00702, n = 8). Gray lines represent female patients, black lines represent male patients, red dots indicate differences. K: Paired line plot representing spindle duration before and after SCC DBS (Pre DBS: 0.773±0.165s, Post DBS: 0.78±0.0639s, t(7) = −0.13, p = 0.901, n = 8). Gray lines represent female patients, black lines represent male patients, red dots indicate differences. L: Paired line plot representing spindle amplitude before and after SCC DBS (Pre DBS: 4.25±0.587, Post DBS: 4.51±0.542, t(7) = −2.17, p = 0.0668, n = 8). Gray lines represent female patients, black lines represent male patients, red dots indicate differences.

Quality control and limitations

Recording LFPs regularly over long periods provides a unique opportunity characterise brain activity associated with healthy and disease-related activity and their modulation by therapies. However, as we and others have previously shown, obtaining such recordings while patients go about their daily lives provides several sources of possible artifact [43,70]. While the majority of the data sets were of sufficient quality for analysing sleep-related activity, data from a number of hemispheres had to be excluded on the basis of ECG artifact or amplifier saturation. One simple quality control measure we used was to look at the correlation of SWA across hemispheres, as during NREM sleep, SWA is found across the entire forebrain [58]. We found SWA in SCC LFPs to be consistent with global synchronisation, as the power of SWA was highly correlated across hemispheres. Both the temporal profile and synchronisation of SWA power were consistent with accurate detection of these activities, though we acknowledge that ground truth polysomnographic sleep scoring would be necessary to conclusively verify this. The sampling density of LFP data in the Pre- and Post-DBS phases was determined based on the storage capacity of the device and the time period between in-clinic visits that would allow for the downloading of the data. Moreover, the timing of some clinical visits themselves led to some gaps in the sampling of the data (as shown in Supplementary Figures S2-S10 for individual patients). As this might have particularly influenced our estimates of SWA timing, we controlled for this by using a resampling approach, which tested the consistency of the Post-DBS SWA fits and metrics with a sparser sampling regime. While it would be ideal for sampling to be consistent between patients and phases, datasets recorded over weeks away from the clinic will often require sparse sampling to preserve battery life. Importantly, we demonstrate that sparse and/or irregular sampling does not prevent meaningful analysis and such datasets remain highly valuable. Moreover, this approach allows data to be recorded with minimal disruption to the patient’s routine by visits to the clinic. Next-generation DBS sensing devices[71–73] may allow more controlled night recordings to replicate and extend these initial observations.

In the Pre-DBS phase, we only included data acquired at least 14 days post-surgery, to minimize the impact of surgical recovery on our analysis. However, it is likely that after-effects of neurosurgery (whether behavioural or neurophysiological) still had an impact on our Pre-DBS results. Specifically, during this time, there was already a significant improvement in depressive symptoms, which could indicate the presence of a lesion effect at this time [53]. Secondly, the requirement to come into the clinic daily during the Post-DBS phase, as noted, could have had an effect on patient behavioural routine and therefore have influenced bed times and their consistency across days. It would be less likely to affect sleep architecture or sleep spindle density. Finally, it should be noted that many antidepressants [66] and benzodiazepines [57] have known effects on sleep architecture. However, despite the variance of medications used across the subjects, all subjects remained on stable doses of their pre-surgical medications throughout the study, so such medication-specific effects are not expected to underlie the consistent findings across subjects.