- 1Child Development Center, University Children's Hospital Zurich, Zurich, Switzerland
- 2Children's Research Center, University Children's Hospital Zurich, Zurich, Switzerland
- 3Department of Pulmonology, University Hospital Zurich, Zurich, Switzerland
- 4University Hospital of Child and Adolescent Psychiatry and Psychotherapy, University of Bern, Bern, Switzerland
- 5Department of Psychology, University of Fribourg, Fribourg, Switzerland
- 6Stichting Epilepsie Instellingen Nederland, Heemstede, Netherlands
- 7Department of Neurology, Leiden University Medical Center, Leiden, Netherlands
- 8Pediatrics and Neurodevelopmental Neuroscience, National Institute of Mental Health, Bethesda, MD, United States
- 9Child Psychiatry Branch, National Institute of Mental Health, Bethesda, MD, United States
- 10Paediatric Neurology and Neurophysiology Unit, Department of Women's and Children's Health, University Hospital of Padova, Padova, Italy
- 11Neuroimmunology Group, Paediatric Research Institute Città della Speranza, Padova, Italy
- 12Department of Pediatric Neurology, Children's Hospital Datteln, Witten/Herdecke University, Datteln, Germany
- 13Department of Medical Sciences, Pediatric Section, University of Ferrara, Ferrara, Italy
- 14Department of Child and Adolescent Psychiatry and Psychotherapy, Psychiatric University Hospital Zurich, University of Zurich, Zurich, Switzerland
- 15Neuroscience Center Zurich, University of Zurich and Swiss Federal Institute of Technology Zurich, Zurich, Switzerland
- 16Zurich Center for Integrative Human Physiology, University of Zurich, Zurich, Switzerland
- 17Translational Research Center, University Hospital of Psychiatry and Psychotherapy, University of Bern, Bern, Switzerland
- 18Department of Pediatric Neurology, University Children's Hospital Zurich, Zurich, Switzerland
- 19Department of Pediatric Neurology, Children's Hospital of Eastern Switzerland, St. Gallen, Switzerland
Background: Sleep disturbances are intertwined with the progression and pathophysiology of psychotic symptoms in schizophrenia. Reductions in sleep spindles, a major electrophysiological oscillation during non-rapid eye movement sleep, have been identified in patients with schizophrenia as a potential biomarker representing the impaired integrity of the thalamocortical network. Altered glutamatergic neurotransmission within this network via a hypofunction of the N-methyl-D-aspartate receptor (NMDAR) is one of the hypotheses at the heart of schizophrenia. This pathomechanism and the symptomatology are shared by anti-NMDAR encephalitis (NMDARE), where antibodies specific to the NMDAR induce a reduction of functional NMDAR. However, sleep spindle parameters have yet to be investigated in NMDARE and a comparison of these rare patients with young individuals with schizophrenia and healthy controls (HC) is lacking. This study aims to assess and compare sleep spindles across young patients affected by Childhood-Onset Schizophrenia (COS), Early-Onset Schizophrenia, (EOS), or NMDARE and HC. Further, the potential relationship between sleep spindle parameters in COS and EOS and the duration of the disease is examined.
Methods: Sleep EEG data of patients with COS (N = 17), EOS (N = 11), NMDARE (N = 8) aged 7–21 years old, and age- and sex-matched HC (N = 36) were assessed in 17 (COS, EOS) or 5 (NMDARE) electrodes. Sleep spindle parameters (sleep spindle density, maximum amplitude, and sigma power) were analyzed.
Results: Central sleep spindle density, maximum amplitude, and sigma power were reduced when comparing all patients with psychosis to all HC. Between patient group comparisons showed no differences in central spindle density but lower central maximum amplitude and sigma power in patients with COS compared to patients with EOS or NMDARE. Assessing the topography of spindle density, it was significantly reduced over 15/17 electrodes in COS, 3/17 in EOS, and 0/5 in NMDARE compared to HC. In the pooled sample of COS and EOS, a longer duration of illness was associated with lower central sigma power.
Conclusions: Patients with COS demonstrated more pronounced impairments of sleep spindles compared to patients with EOS and NMDARE. In this sample, there is no strong evidence that changes in NMDAR activity are related to spindle deficits.
1. Introduction
Once thought to be merely a phenotype of psychosis, sleep dysfunction is now considered to be intricately intertwined in the pathophysiology of disorders along the psychosis spectrum (1–4). Schizophrenia (SZ) is one of several conditions that falls within the psychosis spectrum and is characterized by psychotic symptoms, marking a disconnection from reality, and considerable changes in affective state, perception, and cognition. Disturbed sleep often precedes the onset of SZ and predicts relapses in remitted patients (5–7). Thus, increased research efforts are made to study sleep quality, structure, and electroencephalographic characteristics in order to identify specific parameters potentially contributing to or resulting from the emergence of the disorder. Specifically sleep spindles, a sleep electroencephalogram (EEG) characteristic, caught much attention as a potential endophenotype as well as a target for novel therapeutic approaches (8). Sleep spindles are major oscillations occurring in the sigma frequency range (12–15 Hz) and are most prominent during stage 2 of Non-Rapid Eye Movement (NREM) sleep (3). The reduction of sleep spindle density (SSD; number per minute) in patients affected by SZ when compared to healthy controls (HC) has been the most consistent finding across studies and diverse patient populations with varying age, medication status, or chronicity of the disorder (3, 8–10). This finding is most pronounced in central regions of the cortex in the fast frequency range of spindles (13–15 Hz) (3, 9, 11). Reductions in sleep spindle maximum amplitude (SSMA) and sigma power (SP) are also frequently reported, but tend to be less consistently reduced across studies. This may be in part due to differences in electrode density montages, detection methods, and chronicity of the disorder (1, 3, 9–11).
Sleep spindles are generated by inhibitory GABAergic (gamma-aminobutyric acid) neurons in the thalamic reticular nucleus (12). Thalamic reticular nucleus neurons inhibit excitatory, glutamatergic neurons that propagate signals from the thalamus to the cortex. Cortical projection neurons further communicate with thalamic reticular nucleus neurons via the activation of specific receptors, including the glutamate-specific ionotrophic receptor, the N-methyl-D-aspartate receptor (NMDAR) (13, 14). This NMDAR-mediated glutamatergic feedback loop is necessary to modulate and maintain sleep spindles (15, 16). The NMDAR is also highly implicated in the emergence of SZ. According to the glutamate hypothesis, glutamatergic activity is disrupted in patients with SZ due to a hypofunction of the NMDAR (17, 18). This provokes a cascade of altered neurotransmission that affects brain development and results in the diverse symptomatology of SZ, including cognitive deficits. Sleep spindles play a role in sleep-dependent memory consolidation, and the decrease in sleep spindles in schizophrenia has been hypothesized to contribute to the cognitive deficits that are integral to the disorder (8, 19). It is currently unclear if the hypofunction of the NMDAR is also a key element contributing to the detected sleep spindle reduction in SZ.
A rare form of autoimmune encephalitis, anti-NMDAR encephalitis (NMDARE), is induced by autoantibodies that target the NMDAR leading to a functional deficit (20). The symptomatology of NMDARE and SZ also greatly overlap, often presenting with sleep and cognitive disturbances and culminating with a range of psychotic symptoms (21, 22). This has encouraged research in the context of the glutamate hypothesis of SZ (23–28). Although a previous study found that SSD and SP were only reduced in patients with SZ when comparing them to other patients with mental health disorders who presented with psychosis, sleep spindles have yet to be compared to patient populations with somatic illnesses who present with psychosis (29). The shared phenotype and key pathomechanistic element make SZ and NMDARE very promising to study in concert. Further insight may be gained into the role of the NMDAR within the thalamo-cortical network and the detectable sleep spindle phenotype of both disorders.
Highly dynamic phases of brain development could potentially influence the relationship between sleep and psychosis spectrum disorders. During childhood and adolescence, individuals undergo drastic brain maturation that affects sleep structure and characteristics (25). Additionally, these periods introduce vulnerabilities in the developing cortex which may coincide with the onset of mental health disorders (30, 31). SZ rarely emerges during these young, malleable periods, being classified as Childhood-Onset (COS) if the onset age precedes the age of 13 and Early-Onset Schizophrenia (EOS) if the onset age occurs between age 13 and 18. Previous studies indicate that these populations have reductions in SSD (32, 33), although the two have never been directly compared. In COS, a longer duration of illness was associated with lower SSMA and SP and a recent meta-analysis in adults with SZ points to an association of a longer duration of illness with SSD deficits (33, 34). Combining data from young patients with SZ will allow to further assess the relationship between the duration of illness and sleep spindle parameters.
Pooling sleep data from young patients affected by COS, EOS, and NMDARE, this study aims to assess sleep spindle parameters (SSD, SSMA, SP) between the patient groups and to compare them to age- and sex-matched HC. Due to evidence demonstrating reduced sleep spindles in SZ and the overlapping phenotype and pathomechanism of SZ and NMDARE, it is hypothesized that the pooled patient group will have reduced sleep spindle parameters when compared to HC. Between-patient group comparisons will be performed to elucidate the contribution of each patient population to this presumed finding. Further, to analyze potential regional differences across the cortex, topographical distribution of spindle parameters will be compared between each patient group and matched HC. Lastly, the impact of the duration of illness on spindle alterations will be assessed in an exploratory correlation analysis in young patients with SZ.
2. Materials and methods
2.1. Participants
Sleep data from three patient groups were gathered for the purpose of this study. Screening procedures and characteristics of patient groups have previously been published (22, 32, 33). All-night EEG, demographic, and clinical data from individuals with COS (N = 17, 16.0 ± 3.6 years, 70.6% female), EOS (N = 11, 16.4 ± 1.4 years, 36.4% female), NMDARE (N = 8, 13.3 ± 4.5 years, 75.0% female), and 36 age- and sex-matched HC were pooled in this project. Patients with SZ (COS, EOS) were either recruited at the National Institute of Mental Health in the United States (COS) or at the Department of Child and Adolescent Psychiatry of the Psychiatric University Hospital Zurich in Switzerland (EOS). The diagnosis for patients with COS was assessed by DSM-III-R or DSM-IV criteria and confirmed by two child psychiatrists after an in-patient stay where medication was discontinued (35, 36). For patients with EOS, mental disorders were assessed by the Mini International Neuropsychiatric Interview for Children and Adolescents (MINI-Kid), and criteria for EOS were further verified according to the DSM-IV (37). Age of onset, duration of illness, and antipsychotic medication measured as chlorpromazine equivalent were provided for each patient with SZ. IQ scores were available from 12 patients with COS and eight with EOS. One individual with EOS declined further use of non-genetic personal health data and was excluded in this study.
In the NMDARE group, three individuals had two EEG recordings available. Sleep structure and spindle parameters were calculated individually for each of these EEG and averaged to have one value for each participant. Patient data for this group was provided by three clinics, specifically the Department of Neuropediatrics, University Children's Hospital Zurich (Switzerland) (N = 4; EEG = 6); the Pediatric Neurology and Neurophysiology Unit in Padua (Italy) (N = 2; EEG = 3); the Department of Pediatric Neurology, Children's Hospital in Datteln, University Witten/Herdecke (Germany) (N = 2; EEG = 2). Clinically suspected diagnosis of NMDARE was confirmed by detection of ANMDAR antibodies in serum and/or cerebrospinal fluid.
The three patient groups had different age and sex characteristics (Table 1). Thirty-six HC were used in total; 17 were age- and sex-matched to the COS group (HC-C), accordingly 11 were matched to the EOS group (HC-E), and eight to the NMDARE group (HC-N). HC were recruited at the University Children's Hospital Zurich and underwent a screening process, excluding individuals with a history of mental health disorders and who took psychotropic medication. Written informed consent was obtained from all participants and/or their legal guardian. This research project was approved by swissethics (BASEC 2021-00215). All study procedures were performed according to the Declaration of Helsinki.
2.2. EEG recordings
All-night polysomnography recordings from patients with COS were measured using a 21-channel 10-20 system [TWin or Nihon Kohden; sampling rate (SR): 200 Hz]. EEG data from patients with EOS and HC were recorded using a 128-channel system (Electrical Geodesic Sensory Net; SR: 500 Hz). Nineteen electrodes were shared between patients with SZ and associated HC: two prefrontal (Fp1, Fp2), three frontal (F3, Fz, F4), three central (C3, C4, Cz; Cz available as a recording electrode following re-referencing), six temporal (F7, F8, T3, T4, T5, T6), three parietal (P3, Pz, P4), and two occipital (O1, O2). The patients with NMDARE had heterogeneous EEG recording protocols: (a) Zurich, 21-channel 10-20 system (NeurofileXP EE; SR: 256 Hz), (b) Padua, 8-10 channels (Galileo NT Software; SR: 256 Hz), and (c) Datteln, 11 channels (Neurofax EEG-9210 System; SR: 200 Hz). Seven electrodes were shared across the EEG montages for all three groups (COS, EOS and NMDARE): two prefrontal (Fp1, Fp2), two temporal (T3, T4), two occipital (O1, O2), and one central (Cz).
Epoch-specific sleep staging and artifact rejection were available from previous analyses (22, 32, 33). Sleep stages were scored based on standard criteria by a sleep expert and either validated by a second sleep expert (COS, EOS, HC) or by one rater only to reduce inter-rater variability (NMDARE) (38). Standard sleep stages were assigned to each epoch, composed of 20-s (EOS, HC, NMDARE) or 30-s (COS) time windows. Specific epochs were excluded on a channel to channel basis from the analyses if high levels of artifacts were found. A semi-automatic procedure was used to detect artifacts based on power thresholds in the muscle frequency band (20–40 Hz) and slow-wave frequency band (0.8–4.6 Hz, COS; 0.75–4.5 Hz, EOS, NMDARE) (39, 40). Channels with artifacts in < 97% (COS, EOS, HC) or 95% (NMDARE) of their epochs were considered for subsequent analyses.
2.3. EEG data analyses
EEG data were analyzed in MATLAB 2017b. Data were low-pass filtered (−6 dB cut-off = 40.5 Hz; Finite Impulse Response (FIR) filter), down-sampled to 125 Hz, and high-pass filtered (−6 dB cut-off = 0.38 Hz, Kaiser Window FIR filter). Electrodes were re-referenced to the mean of the occipital channels (O1, O2). The occipital channels were chosen as they are shared by all participants. Further, all other possible channel locations shared by all EEG montages (frontal, central, temporal) have previously shown disturbances in spindles in SZ (8, 32, 41–46).
Sleep spindles were analyzed during the first hour of artifact-free NREM Stage 2 and 3 for each participant. This is considered the most robust part of sleep, with little or no fragmentation. Sleep spindles were automatically detected based on an established amplitude-threshold algorithm (10). EEG signal was band-pass filtered (−6 dB cut-offs: 12.01 and 16.97 Hz, Chebyshev Type II filter) and the average signal amplitude was calculated for each channel. As signal amplitude varies for each electrode, a threshold was set relative to the average signal amplitude for each channel. Therefore, sleep spindles were detected if fluctuations in amplitude for a given channel surpassed an upper threshold of five times the average signal amplitude. The upper threshold value was determined based on previous studies examining young adults and was confirmed as the best representation of spindle detection after visual inspection (47, 48). The beginning and end of each spindle were defined as the instances where the signal surrounding the peak amplitude dipped below a lower threshold of two times the average signal. The lower duration threshold of a spindle was determined to be 0.3 s (49).
Sleep spindle density (SSD), maximum amplitude (SSMA), and sigma power (SP) were investigated in this project. SSD refers to the number of spindles per minute (13.75–14.75 Hz). SSMA refers to the largest amplitude value (positive or negative; μV) of a given spindle and was calculated for every spindle within the frequency range of interest and averaged together. Spectral analysis of EEG signals were performed for each epoch (Fast Fourier transform, 4-s Hanning window, overlapping by 2-s), resulting in a 0.25 Hz frequency resolution. Sigma power (13.75–14.75 Hz; μV2) was calculated for each epoch.
2.4. Statistical analyses
Sleep structure parameters and demographic and clinical characteristics were compared between patient groups and their associated HC (HC-C, HC-E, HC-N) using non-parametric Mann-Whitney U-tests. The same variables were compared between patient groups using non-parametric Kruskal Wallis tests. Multiple comparisons and determining the direction of interaction terms was carried out with Tukey-Kramer post-hoc tests for pairwise comparisons. Significant differences were considered with a p-value < 0.05. Sleep structure parameter values were calculated for either (1) the first hour of artifact-free NREM Stage 2 and 3 for each participant (the same time frame used in the sleep spindle analysis), or (2) the longest sleep recording period for each participant. The comparison of sleep structure considers the five recording electrodes shared by all groups. It was not possible to calculate sleep latency or efficiency for the NMDARE group as the time when lights were turned off was not documented.
In order to isolate the frequency band and area of interest, the a priori hypothesis regarding SSD in psychosis was investigated. This hypothesis claims that central SSD in the fast spindle frequency range is reduced in patients with psychosis compared to HC. Unpaired t-tests assessed SSD differences between the collective patient group (PSY; N = 36) and the entire HC sample (N = 36). A test was conducted for each recording electrode shared by all participants (Fp1, Fp2, T3, T4, Cz) at each frequency bin of interest (12–15 Hz, discretized into 0.25 steps). Values from the 12–12.5 and 15 Hz frequency bins were excluded from this analysis as the filter attenuated < 20% of the signal amplitude or there were low levels of spindles, respectively. Multiple comparisons were corrected for with Benjamini and Hochberg False Discovery Rate statistical approach. Furthermore, the subsequent spindle analyses focused on values pertaining to the 13.75–14.75 Hz frequency band in the central electrode (Figure 1).
Figure 1. Sleep spindle density differences between PSY and HC for all channels and frequency bins. Sleep spindle density (SSD) is presented for each recording channel for all relevant frequency bins of interest for the collective patient group [PSY; (A)] and collective healthy control group [HC; (B)]. In (C), the percent change in SSD is directly compared between the two groups; an increase in graph (C) represents an increase in SSD in PSY compared to HC. Frequency bin-channel combinations that were significantly different from unpaired t-tests between the two groups are outlined in black.
In addition to SSD, SSMA and SP were compared between the PSY and HC groups with unpaired t-tests. Differences in spindle parameters between the three patient groups were analyzed with one-way analysis of the variance (ANOVA); Tukey-Kramer post-hoc tests were applied subsequently. These two analyses were performed using spindle values specific to the central electrode (Cz). Due to large amounts of artifacts in the central electrode of two patients with NMDARE, they and their age- and sex-matched HC were not included in two above mentioned analyses. Topographical analyses of sleep spindles between each patient group and associated HC were performed with electrode-wise Student's unpaired t-tests; multiple comparisons were corrected for with Benjamini and Hochberg False Discovery Rate statistical approach. The comparison considers 17 electrodes for the COS and EOS groups and five electrodes for the NMDARE group. Normal distribution of the data was visualized with quantile-quantile plots and objectively confirmed with Shapiro-Wilk tests.
Furthermore, to assess potential relationships between sleep spindle parameters and illness duration, partial Pearson correlation analysis were performed controlling for age and chlorpromazine equivalent. This analysis was conducted in patients with SZ (COS, EOS). As only patients with SZ were included in this analysis, the value of sleep spindle parameters were taken from the mean of all available central electrodes (C3, Cz, C4).
3. Results
3.1. Demographic and clinical characteristics and sleep parameters
3.1.1. Demographic and clinical characteristics
No significant difference was found for age and sex between patient groups and associated HC and inter-patient groups (see Table 1). The COS group had a significantly longer duration of illness and a lower age of disease onset when compared to the EOS group. Additionally, the COS group had significantly higher levels of medication, measured in chlorpromazine equivalent, than the EOS group. No significant differences were found in IQ between the COS and EOS groups.
3.1.2. Sleep parameters
Differences in sleep structure and quality for the entire night are reported in Table 2. Patients with COS took longer to fall asleep and spent more time asleep compared to HC-C. After initial sleep onset, they spent less time awake during the entire night. Sleep characteristics were comparable between the patients with EOS and HC-E. Patients with NMDARE slept for significantly longer than HC-N. Differences in sleep parameters were also detected when comparing patient groups. Patients with EOS slept for a shorter amount of time than those with COS. Patients with COS spent less time awake during the entire night following initial sleep onset than those with NMDARE. A statistical difference in percent of all-night time spent in NREM sleep stage 1 was found in the omnibus test, yet no significant differences between the means of each group (COS, EOS, NMDARE) were detected when assessing with post-hoc multiple comparisons.
Sleep parameters were also calculated for the time frame used in the sleep spindle analysis (Table 3). This refers to the amount of time it took to isolate 60 minu of artifact-free NREM stage 2 and 3 sleep from each participant. This amount of time was shorter for patients with COS compared to HC-C. Patients with COS also spent less time awake following initial descent into sleep and time in NREM sleep stage 1 during the first hour. The amount of time needed to isolate 60 min of artifact-free sleep in the NMDARE group was significantly longer as when compared to HC-N. Additionally, they spent more time awake and in REM sleep during the first hour compared to HC-N. When comparing patient groups, it was found that the absolute length of the first hour was less in COS when compared to NMDARE. Additionally, the COS group spent less time awake and in NREM sleep stage 1 compared to the EOS and NMDARE group. Finally, both the COS and EOS groups spent less time than the NMDARE group in REM sleep.
3.2. Sleep spindle parameters
3.2.1. Sleep spindle differences between PSY and HC and patient groups
Based on results from the unpaired t-test, sleep spindle parameter values from the frequency range of 13.75–14.5 Hz pertaining to the Cz electrode were used to compare the PSY and HC groups and the three patient groups. Results assessing SSD, SSMA, and SP differences between the PSY and HC groups are presented in the left column of Figure 2. For all parameters, the PSY group showed significant reductions in spindle characteristics when compared to the HC group [SSD: t(66) = −5.24, p < 0.00001; SSMA: t(66) = −3.43, p = 0.0011; SP: t(66) = −4.20, p < 0.0001].
Figure 2. Sleep spindle parameter differences between PSY and HC & patient groups. Differences in sleep spindle parameters between the collective patient groups (PSY) and collective healthy control group (HC) are show in the left panel and between individual patient groups composing the PSY group [Childhood-Onset Schizophrenia (COS), Early-Onset Schizophrenia (EOS), anti-N-methyl-D-aspartate receptor encephalitis (NMDARE)] in the right panel. Note that significant differences between individual patient groups and associated HC are not presented here. All spindle values here are taken from the central electrode from the frequency band of 13.75–14.75 Hz. The * denotes a significant difference between the groups according to p-values (<0.05 implies significance). ** = p ≤ 0.01; **** = p ≤ 0.0001.
ANOVA conducted comparing the three groups and SSD, SSMA, and SP are presented in the right column of Figure 2. No differences in SSD were found between the three groups [F(2, 31) = 2.59, p = 0.091]. For SSMA and SP, the COS group was significantly reduced compared to the EOS [SSMA: p < 0.0001, 95% C.I. = (−10.58, −3.49); SP: p < 0.001, 95% C.I. = (−4.03, −1.06)] and the NMDARE groups [SSMA: p < 0.0001, 95% C.I. = (−13.49, −4.79); SP: p < 0.001, 95% C.I. = (−4.94, −1.29)]. No significant differences in SSMA [F(2, 31) = 19.19, p = 0.51] and SP [F(2, 31) = 13.4, p = 0.75] were found between the EOS and NMDARE groups.
3.2.2. Topographical distribution of sleep spindles between patient groups and associated HC
Topographical comparisons of sleep spindle parameters between patient groups and associated HC are presented in Figure 3. As previously reported (33), the COS group shows global significant decreases in all sleep spindle parameters compared to HC-C. As expected, the EOS group shows significant reductions in SSD in the central region and left, posterior temporal region compared to HC-E (32). Sleep spindles for the NMDARE group were statistically comparable to the HC-N group.
Figure 3. Topographical distribution of sleep spindle parameters in patient groups compared to associated HC. Two-tailed Student's unpaired t-tests with standard alpha-values at 0.05 were used; significant differences between the groups in individual electrodes are designated by white dots. The colors represent the percent change of the given sleep spindle parameter when comparing the patient group [Childhood-Onset Schizophrenia (COS), Early-Onset Schizophrenia (EOS), and anti-N-methyl-D-aspartate receptor encephalitis (NMDARE)] to the associated healthy control groups (HC-C, HC-E, and HC-N, associated with COS, EOS and NMDARE, respectively) group; e.g., a percent increase represents values that are higher in the patient group compared to HC and a percent decrease represents values that are lower in the patient group compared to HC.
3.3. Association between sleep spindles and duration of the disorder
For all patients with SZ, SP derived from the mean of central electrodes (C3, Cz, C4) was found to be significantly negatively correlated with the duration of illness (r = −0.41, p = 0.035). Duration of illness was not significantly associated with SSD (r = −0.11, p = 0.603). There was a trend for a negative correlation between duration of illness and SSMA (r = −0.37, p = 0.061). All sleep spindle parameters and duration of illness are presented in Figure 4.
Figure 4. The relationship between sleep spindle parameters in patients with SZ and duration of illness. Rho and p-values are presented for each sleep spindle parameter when compared to duration of illness in months. All values come from the mean of the central electrodes (C3, C4, Cz) and pertain to either the Childhood-Onset Schizophrenia (COS) group (light blue) or the Early-Onset Schizophrenia (EOS) group (dark blue).
4. Discussion
The main finding of this study is that all sleep spindle parameters (SSD, SSMA, SP) were reduced in the pooled sample of patients affected by SZ or NMDARE compared to matched HC. SSD was statistically comparable between the three patient groups, however patients with NMDARE showed no significant alterations in sleep spindle parameters when compared to HC-N. Young patients affected by SZ showed reductions in SSD compared to HC; the COS group also exhibited reductions in SSMA and SP. Further, in patients with SZ, reduced SP was associated with a longer duration of illness.
4.1. Sleep spindle parameters in patients with SZ
Consistent with previous literature, reduced SSD was found in patients with SZ compared to HC (9, 29, 32, 33). Topographical analysis of SSD indicated widespread reductions in patients with COS, whereas deficits in patients with EOS were restricted to primarily central electrodes. No significant differences were found in SSD measured from the central electrode between the COS and EOS groups. Patients with COS also showed global reductions in SSMA and SP, whereas no changes in these spindle parameters were found in the EOS group when comparing to HC. The COS and EOS groups differ by definition in the age of disease onset, and in this study the COS group also had a longer duration of illness and higher chlorpromazine equivalent than the EOS group. The longer duration of illness is explained by the fact that sleep data from patients with EOS were recorded promptly following disease onset, whereas data from patients with COS were recorded during follow-up assessments years after the initial emergence of the disorder. This is further illustrated by the fact that the mean age at the time of the sleep assessment did not differ between the groups. Thus, in terms of duration of illness, patients with EOS rather resemble adult samples with early course psychosis or first episode SZ whereas young patients with COS may rather be compared to patients with a chronic course of SZ. Considered in this light, the findings are in line with previous literature, namely the vast reduction of SSD in COS corresponds to global reduced SSD in adults with chronic SZ (10, 50). Likewise, the local reduction of SSD in patients with EOS compared to HC-E is consistent with findings in adult patients experiencing their first episode of psychosis (51). Additionally, SSMA and SP reductions in patients with COS are similar to findings from adult samples with chronic SZ (11). When comparing adult patients with chronic SZ to HC, the effect size of the reduction was highest for SSD, followed by SP, and then SSMA. In contrast, in adult first episode psychosis or early course SZ patients there was no significant reduction in SSMA, similar to the results regarding the EOS group (29, 51). Many factors influence spindle characteristics, including genes (52), illness duration (34, 43), and severity of symptomatology (42). Given the current sample size and heterogeneity in our population, we were unable to individually examine the impact of these parameters on spindle characteristics, however, this may be a fruitful future avenue.
The potential impact of the duration of illness in the COS and EOS groups was further explored, resulting in a significant relationship between reduced SP and longer duration of illness. Although individual studies have reported no association between sleep spindle activity and duration of illness or chlorpromazine equivalent in their samples (10, 53), a recent meta-analysis including eight studies demonstrated that across all individuals, the longer the illness lasted, the lower the SSD (34). Specifically for every year of illness, SSD decreased by 0.2 spindles per minute. In conclusion, it is assumed that in this sample, as in previous single studies with similar numbers of patients, the correlation analysis of SSD with duration of illness may have been underpowered. Additionally, while age-specific aspects of sleep spindles may be accounted for by age- and sex-matching when comparing patients and HC, directly comparing the patient sample may be confounded by developmental aspects of sleep spindles across the whole age-range (14, 54).
4.2. Sleep spindle parameters in young patients with NMDARE
Central SSD was reduced in the collective patient group compared to HC and no significant differences were found between the patient groups. This suggests that all patient groups contribute to the low, central SSD. However, no significant differences in any of the sleep spindle parameters were found when comparing NMDARE to HC-N. This finding may at least in part be due to the sample size and inter-individual variability making any definite conclusions challenging. Nonetheless, strong and global deficits in sleep spindles as found in patients with COS might have been detected irrespective of these limitations.
The exact molecular mechanism and pathways of neurotransmission underlying spindle deficits in patients with SZ is unclear. It has been suggested that a reduced binding or expression of the glutamatergic NMDAR within the thalamo-cortical system may be the key element contributing to early reductions in SSD (51). Further, postmortem studies found reductions of the NMDAR in the thalamus and in the cortex of patients with SZ (55). On the other hand, based on the premise that the spindle-pacemaker relies on GABAergic neurons, it has been hypothesized that GABA deficits may trigger the reduction in SSD whereas a hypofunction of the NMDAR within the thalamo-cortical feedback loop may rather impact the maintenance and modulation of sleep spindles and their coordination with other oscillations (53). Accordingly, the administration of eszopiclone, a GABA agonist targeting the thalamic reticular nucleus, in patients affected by SZ results in increased SSD (56). Additionally, an in vitro model showed that a blockade of the NMDAR shortened the spindle-like oscillation, but only the combined blockade of the NMDAR and GABA receptor repressed the spindle-like oscillation (13). Thus, the findings suggesting that patients with NMDARE do not have large spindle deficits rather support previous literature and the hypothesis that NMDAR-mediated hypofunction alone does not account for pronounced deficits in SSD. Specifically, the explicit antibody-mediated reduction of functional NMDAR in NMDARE may not entirely mimic the far less clear beginning and timing of abnormalities on the NMDAR during development and its impact on GABAergic activity and thalamocortical integrity in SZ (57, 58). Although speculative, it is possible that further research in NMDARE with increased sample size, age range, and density of electrodes may uncover more subtle reductions of SSD, perhaps making it similar to early stages of SZ such as the patient group with EOS.
4.3. Limitations and outlook
The results of the present study should be interpreted within the confines of several limitations. The rarity of NMDARE and Childhood- and Early-Onset SZ make the data of the included individuals very valuable, but the low sample size renders the statistical analyses less robust and age- and sex-matching between the patient groups challenging. Additionally, sleep disorders in the NMDARE group of patients were not assessed systematically. Thus, it remains difficult to discern and weigh the impact of the sleep spindle parameter findings, as well as the relationship with duration of illness. Further, analysis of the sleep structure showed group-specific differences that may be due to different medication and recording settings. However, sleep spindles predominantly appear in sleep stage N2, which was comparable between all groups.
Assessing sleep spindles and their coordination with other oscillations, specifically slow waves, in longitudinal studies may be essential to further discern the pathomechanism of SZ and to identify targets for novel interventions to treat this debilitating disorder. Likewise, sleep studies in larger samples of patients with NMDARE or with pharmacological models mimicking the NMDAR hypofunction may further elucidate similarities and dissimilarities of underlying signaling and thalamo-cortical impairments.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving human participants were reviewed and approved by Ethics Committee of the Canton of Zurich and Bern, Switzerland, BASEC 2021-00215. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.
Author contributions
AM, SG, MN, KR, SS, AS, SK, BB, and MG previously collected, provided, and/or analyzed patient data. BB and MG contributed to the conception of the project and the study design. MG additionally contributed to the writing process. MD contributed to the project design, data analysis, and writing process. MG, BB, LT, and RH closely supervised the project. All authors reviewed and approved the article for submission.
Funding
This work was supported by the Frutiger Foundation (to MG), Swiss National Science Foundation 320030_153387 (to RH), PCEFP1-181279 (to SK), and 32003B_184943 (to LT), radiz—Rare Disease Initiative Zurich (to SG), Interfaculty Research Cooperation: Decoding Sleep (to LT), the Anna Müller Grocholski foundation (to BB), and the EMDO foundation (to BB and MG). This research was also supported in part by the Intramural Research Program of the National Institute of Mental Health (Annual Report Number ZIAMH002581, ClinicalTrials.gov identifier NCT00001198, protocol ID 84-M-0050). Additionally, MN and SS were supported by the project DOR2008453/20, University of Padova (systematic literature review and evidence-based international consensus recommendations for the treatment of pediatric anti-NMDAR encephalitis, MN).
Acknowledgments
We would like to thank the young participants who volunteered their time and/or data for the project, as well as their families for supporting their participation. We would also like to thank the whole of RH's lab for their crucial advice and support.
Conflict of interest
In the past five years, SW has received royalties from Thieme Hogrefe, Kohlhammer, Springer, and Beltz, and her work was supported by the Swiss National Science Foundation (SNF), diff. EU FP7s, Bfarm Germany, ZInEP, Hartmann Müller, Olga Mayenfisch, Gertrud Thalmann, Vontobel, Unicentia, and Erika Schwarz Funds.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Kaskie RE, Graziano B, Ferrarelli F. Schizophrenia and sleep disorders: links, risks, and management challenges. Nat Sci Sleep. (2017) 9:227–33. doi: 10.2147/NSS.S121076
2. Winsky-Sommerer R, de Oliveira P, Loomis S, Wafford K, Dijk DJ, Gilmour G. Disturbances of sleep quality, timing and structure and their relationship with other neuropsychiatric symptoms in Alzheimer's disease and schizophrenia: Insights from studies in patient populations and animal models. Neurosci Biobehav Rev. (2018) 97:112–37. doi: 10.1016/j.neubiorev.2018.09.027
3. Zhang Y, Quiñones GM, Ferrarelli F. Sleep spindle and slow wave abnormalities in schizophrenia and other psychotic disorders: Recent findings and future directions. Schizophr Res. (2020) 221:29–36. doi: 10.1016/j.schres.2019.11.002
4. Vukadinovic Z. Sleep abnormalities in schizophrenia may suggest impaired trans-thalamic cortico-cortical communication: towards a dynamic model of the illness. Eur J Neurosci. (2011) 34:1031–9. doi: 10.1111/j.1460-9568.2011.07822.x
5. Benson KL. Sleep in schizophrenia: impairments, correlates, and treatment. Psychiatr Clin North Am. (2006) 29:1033–45. doi: 10.1016/j.psc.2006.08.002
6. Lunsford-Avery JR, LeBourgeois MK, Gupta T, Mittal VA. Actigraphic-measured sleep disturbance predicts increased positive symptoms in adolescents at ultra high-risk for psychosis: a longitudinal study. Schizophr Res. (2015) 164:15–20. doi: 10.1016/j.schres.2015.03.013
7. Ruhrmann S, Schultze-Lutter F, Salokangas RK, Heinimaa M, Linszen D, Dingemans P, et al. Prediction of psychosis in adolescents and young adults at high risk: results from the prospective European prediction of psychosis study. Arch Gen Psychiatry. (2010) 67:241–51. doi: 10.1001/archgenpsychiatry.2009.206
8. Manoach DS, Pan JQ, Purcell SM, Stickgold R. Reduced sleep spindles in schizophrenia: a treatable endophenotype that links risk genes to impaired cognition? Biol Psychiatry. (2015) 80:1–10. doi: 10.1016/j.biopsych.2015.10.003
9. D'Agostino A, Castelnovo A, Cavallotti S. Sleep endophenotypes of schizophrenia: slow waves and sleep spindles in unaffected first-degree relatives. npj Schizophr. (2018) 4:2. doi: 10.1038/s41537-018-0045-9
10. Ferrarelli F, Huber R, Peterson MJ, Massimini M, Murphy M, Riedner BA, et al. Reduced sleep spindle activity in schizophrenia patients. Am J Psychiatry. (2007) 164:483–92. doi: 10.1176/ajp.2007.164.3.483
11. Castelnovo A, Graziano B, Ferrarelli F. D'Agostino A. Sleep spindles and slow waves in schizophrenia and related disorders: main findings, challenges and future perspectives. Eur J Neurosci. (2018) 48:2738–58. doi: 10.1111/ejn.13815
12. Steriade M, Domich L, Oakson G, Deschenes M. The deafferented reticular thalamic nucleus generates spindle rhythmicity. J Neurophysiol. (1987) 57:260–73. doi: 10.1152/jn.1987.57.1.260
13. Jacobsen RB, Ulrich D, Huguenard JR. GABA(B) and NMDA receptors contribute to spindle-like oscillations in rat thalamus in vitro. J Neurophysiol. (2001) 86:1365–75. doi: 10.1152/jn.2001.86.3.1365
14. Fernandez L, Lüthi A. Sleep spindles: mechanisms and functions. Physiol Rev. (2020) 100:805–68. doi: 10.1152/physrev.00042.2018
15. Astori S, Wimmer RD, Lüthi A. Manipulating sleep spindles - expanding views on sleep, memory, and disease. Trends Neurosci. (2013) 6:738–48. doi: 10.1016/j.tins.2013.10.001
16. Manoach DS, Stickgold R. Abnormal sleep spindles, memory consolidation, and schizophrenia. Annu Rev Clin Psychol. (2019) 15:451–79. doi: 10.1146/annurev-clinpsy-050718-095754
17. Coyle JT. NMDA receptor and schizophrenia: a brief history. Schizophr Bull. (2012) 38:920–6. doi: 10.1093/schbul/sbs076
18. Kayser MS, Dalmau J. Anti-NMDA receptor encephalitis, autoimmunity, and psychosis. Schizophr Res. (2016) 176:36–40. doi: 10.1016/j.schres.2014.10.007
19. Manoach DS, Mylonas D, Baxter B. Targeting sleep oscillations to improve memory in schizophrenia. Schizophr Res. (2020) 221:63–70. doi: 10.1016/j.schres.2020.01.010
20. Dalmau J, Gleichman AJ, Hughes EG, Rossi JE, Peng X, Lai M, et al. Anti-NMDA-receptor encephalitis: case series and analysis of the effects of antibodies. Neurology. (2008) 7:1091–8. doi: 10.1016/S1474-4422(08)70224-2
21. Steullet P, Cabungcal JH, Monin A, Dwir D, O'Donnell P, Cuenod M, et al. Redox dysregulation, neuroinflammation, and NMDA receptor hypofunction: A “central hub” in schizophrenia pathophysiology? Schizophr Res. (2016) 176:41–51. doi: 10.1016/j.schres.2014.06.021
22. Gefferie SR, Maric A, Critelli H, Gueden S, Kurlemann G, Kurth S, et al. Altered EEG markers of synaptic plasticity in a human model of NMDA receptor deficiency: anti-NMDA receptor encephalitis. NeuroImage. (2021) 239:118281. doi: 10.1016/j.neuroimage.2021.118281
23. Lee G, Zhou Y. NMDAR hypofunction animal models of schizophrenia. Front Mol Neurosci. (2019) 12:185. doi: 10.3389/fnmol.2019.00185
24. Marsman A, Van Den Heuvel MP, Klomp DW, Kahn RS, Luijten PR, Hulshoff Pol HE. Glutamate in schizophrenia: a focused review and meta-analysis of 1H-MRS studies. Schizophr Bull. (2013) 39:120–9. doi: 10.1093/schbul/sbr069
25. Uhlhaas PJ, Singer W. The development of neural synchrony and large-scale cortical networks during adolescence: relevance for the pathophysiology of schizophrenia and neurodevelopmental hypothesis. Schizophr Bull. (2011) 37:514–23. doi: 10.1093/schbul/sbr034
26. Moghaddam B, Javitt D. From revolution to evolution: the glutamate hypothesis of schizophrenia and its implication for treatment. Neuropsychopharmacol. (2012) 37:4–15. doi: 10.1038/npp.2011.181
27. McCutcheon RA, Krystal JH, Howes OD. Dopamine and glutamate in schizophrenia: biology, symptoms and treatment. World Psychiatry. (2020) 19:15–33. doi: 10.1002/wps.20693
28. Florance NR, Davis RL, Lam C, Szperka C, Zhou L, Ahmad S, et al. Anti-N-methyl-D-aspartate receptor (NMDAR) encephalitis in children and adolescents. Ann Neurol. (2009) 66:11–8. doi: 10.1002/ana.21756
29. Manoach DS, Demanuele C, Wamsley EJ, Vangel M, Montrose DM, Miewald J, et al. Sleep spindle deficits in antipsychotic-naïve early course schizophrenia and in non-psychotic first-degree relatives. Front Hum Neurosci. (2014) 8:762. doi: 10.3389/fnhum.2014.00762
30. Paus T, Keshavan M, Giedd JN. Why do many psychiatric disorders emerge during adolescence? Nat Rev Neurosci. (2008) 9:947–57. doi: 10.1038/nrn2513
31. Shaw P, Kabani NJ, Lerch JP, Eckstrand K, Lenroot R, Gogtay N, et al. Neurodevelopmental trajectories of the human cerebral cortex. J Neurosci. (2008) 28:3586–94. doi: 10.1523/JNEUROSCI.5309-07.2008
32. Gerstenberg M, Furrer M, Tesler N, Franscini M, Walitza S, Huber R. Reduced sleep spindle density in adolescent patients with early-onset schizophrenia compared to major depressive disorder and healthy controls. Schizophr Res. (2020) 221:20–8. doi: 10.1016/j.schres.2019.11.060
33. Markovic A, Buckley A, Driver DI, Dillard-Broadnax D, Gochman PA, Hoedlmoser K, et al. Sleep spindle activity in childhood onset schizophrenia: diminished and associated with clinical symptoms. Schizophr Res. (2020) 223:327–36. doi: 10.1016/j.schres.2020.08.022
34. Lai M, Hegde R, Kelly S, Bannai D, Lizano P, Stickgold R, et al. Investigating sleep spindle density and schizophrenia: a meta-analysis. Psychiatry Res. (2022) 307:114265. doi: 10.1016/j.psychres.2021.114265
35. American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders (3rd ed., revised) (1987).
36. American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders (4th ed.) (1994).
37. Sheehan DV, Sheehan KH, Shytle RD, Janavs J, Bannon Y, Rogers JE, et al. Reliability and validity of the Mini International Neuropsychiatric Interview for Children and Adolescents (MINI-KID). J Clin Psychiatry. (2010) 71:313–26. doi: 10.4088/JCP.09m05305whi
38. Iber C, Ancoli-Israel SAC. The AASM manual for the scoring of sleep and associated events: rules, terminology and technical specifications. Am Acad Sleep Med. (2007).
39. Huber R, Deboer T, Tobler I. Sleep deprivation in prion protein deficient mice and control mice: genotype dependent regional rebound. Neuroreport. (2002) 13:1–4. doi: 10.1097/00001756-200201210-00005
40. Buckelmüller J, Landolt HP, Stassen HH, Achermann P. Trait-like individual differences in the human sleep electroencephalogram. Neuroscience. (2006) 138:351–6. doi: 10.1016/j.neuroscience.2005.11.005
41. Cox R, Schapiro AC, Manoach DS, Stickgold R. Individual differences in frequency and topography of slow and fast sleep spindles. Front Hum Neurosci. (2017) 11:241–51. doi: 10.3389/fnhum.2017.00433
42. Ferrarelli F, Peterson MJ, Sarasso S, Riedner BA, Murphy MJ, Benca RM, et al. Thalamic dysfunction in schizophrenia suggested by whole-night deficits in slow and fast spindles. Am J Psychiatry. (2010) 167:1339–48. doi: 10.1176/appi.ajp.2010.09121731
43. Markovic A, Buckley A, Driver DI, Dillard-Broadnax D, Gochman PA, Hoedlmoser K, et al. Sleep neurophysiology in childhood onset schizophrenia. J Sleep Res. (2020) 30:e13039. doi: 10.1111/jsr.13039
44. Yao D, Qin Y, Hu S. Which reference should we use for EEG and ERP practice? Brain Topogr. (2019) 32:530–49. doi: 10.1007/s10548-019-00707-x
45. Leuchs L. Choosing your reference - and why it matters. Brain Prod. Solut Neurohysiol Res. (2019). Available online at: https://pressrelease.brainproducts.com/referencing/
46. Vukadinovic Z. Sleep spindle reductions in schizophrenia and its implications for the development of cortical body map. Schizophr Res. (2015) 168:589–90. doi: 10.1016/j.schres.2015.08.012
47. Lustenberger C, Wehrle F, Tüshaus L, Achermann P, Huber R. The multidimensional aspects of sleep spindles and their relationship to word-pair memory consolidation. Sleep. (2015) 38:1093–103. doi: 10.5665/sleep.4820
48. Dijk DJ, EEG. slow waves and sleep spindles: windows on the sleeping brain. Behav Brain Res. (1995) 69:109–16. doi: 10.1016/0166-4328(95)00007-G
49. Warby SC, Wendt SL, Welinder P, Munk EG, Carrillo O, Sorensen HB, et al. Sleep-spindle detection: crowdsourcing and evaluating performance of experts, non-experts and automated methods. Nat Methods. (2014) 11:385–92. doi: 10.1038/nmeth.2855
50. Manoach DS, Thakkar KN, Stroynowski E, Ely A, McKinley SK, Wamsley E, et al. Reduced overnight consolidation of procedural learning in chronic medicated schizophrenia is related to specific sleep stages. J Psychiatr Res. (2010) 44:112–20. doi: 10.1016/j.jpsychires.2009.06.011
51. Kaskie RE, Gill KM, Ferrarelli F. Reduced frontal slow wave density during sleep in first-episode psychosis. Schizophr Res. (2019) 206:318–24. doi: 10.1016/j.schres.2018.10.024
52. Rusterholz T, Hamann C, Markovic A, Schmidt SJ, Achermann P, Tarokh L. Nature and nurture: brain region specific inheritance of sleep neurophysiology in adolescence. J Neurosci. (2018) 38:9275–85. doi: 10.1523/JNEUROSCI.0945-18.2018
53. Wamsley EJ, Tucker MA, Shinn AK, Ono KE, McKinley SK, Ely AV, et al. Reduced sleep spindles and spindle coherence in schizophrenia: mechanisms of impaired memory consolidation? Biol Psychiat. (2012) 71:154–61. doi: 10.1016/j.biopsych.2011.08.008
54. Purcell S, Manoach D, Demanuele C. Characterizing sleep spindles in 11,630 individuals from the National Sleep Research Resource. Nat Commun. (2017) 8:15930. doi: 10.1038/ncomms15930
55. Pakkenberg B. Scheel-Krüger J, Kristiansen LV. Schizophrenia; from structure to function with special focus on the mediodorsal thalamic prefrontal loop. Acta Psychiatr Scand. (2009) 120:345–54. doi: 10.1111/j.1600-0447.2009.01447.x
56. Mylonas D, Baran B, Demanuele C, Cox R, Vuper TC, Seicol BJ, et al. The effects of eszopiclone on sleep spindles and memory consolidation in schizophrenia: a randomized clinical trial. Neuropsychopharmacol. (2020) 45:2189–97. doi: 10.1038/s41386-020-00833-2
57. Cohen SM, Tsien RW, Goff DC, Halassa MM. The impact of NMDA receptor hypofunction on GABAergic neurons in the pathophysiology of schizophrenia. Schizophr Res. (2015) 167:98–107. doi: 10.1016/j.schres.2014.12.026
Keywords: psychosis, schizophrenia, anti-NMDAR encephalitis, sleep EEG, sleep spindles, thalamocortical network
Citation: Dimitriades ME, Markovic A, Gefferie SR, Buckley A, Driver DI, Rapoport JL, Nosadini M, Rostasy K, Sartori S, Suppiej A, Kurth S, Franscini M, Walitza S, Huber R, Tarokh L, Bölsterli BK and Gerstenberg M (2023) Sleep spindles across youth affected by schizophrenia or anti-N-methyl-D-aspartate-receptor encephalitis. Front. Psychiatry 14:1055459. doi: 10.3389/fpsyt.2023.1055459
Received: 27 September 2022; Accepted: 19 May 2023;
Published: 07 June 2023.
Edited by:
Izabela Guimaraes Barbosa, Federal University of Minas Gerais, BrazilReviewed by:
Siyuan Fan, Peking Union Medical College Hospital (CAMS), ChinaTakashi Kanbayashi, University of Tsukuba, Japan
Copyright © 2023 Dimitriades, Markovic, Gefferie, Buckley, Driver, Rapoport, Nosadini, Rostasy, Sartori, Suppiej, Kurth, Franscini, Walitza, Huber, Tarokh, Bölsterli and Gerstenberg. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Miriam Gerstenberg, miriam.gerstenberg@pukzh.ch