editor@360hinduism.com
| | |

Experienced Mindfulness Meditators Exhibit Higher Parietal-Occipital EEG Gamma Activity during NREM Sleep

Fabio Ferrarelli, Richard Smith, Daniela Dentico, Brady A. Riedner, Corinna Zennig, Ruth M. Benca, Giulio Tononi
Department of Psychiatry, University of Wisconsin-Madison, Madison, Wisconsin, United States of America

Antoine Lutz
Waisman Center for Brain Imaging and Behavior, University of Wisconsin-Madison, Madison, Wisconsin, United States of America, & Lyon Neuroscience Research Center, Lyon 1 University, Lyon, France

Richard J. Davidson
Waisman Center for Brain Imaging and Behavior, University of Wisconsin-Madison, Madison, Wisconsin, United States of America & Department of Psychology, University of Wisconsin-Madison, Madison, Wisconsin, United States of America,

Introduction

Meditation can be conceptualized as a set of regulatory and selfinquiry mental training regimes cultivated for various ends, including the training of well-being and psychological health [1,2],[3]. Accumulating evidence suggests that meditation training induces functional and anatomical neuronal changes [1],[2,3]. These neuroplastic changes are linked to changes in behavior during cognitive and affective tasks [4],[3]. The investigation of spontaneous brain activity, at rest or during practice, is a sensitive approach to identify these neuroplastic changes. Some studies found meditation-related increases in alpha and theta frequency bands, particularly in frontal midline and prefrontal areas, which were associated with increased relaxation (reviewed in [5]). Another study showed that long-term meditators (LTM) had high-amplitude gamma band (25–40 Hz) oscillations during mental practice localized to lateral frontal and posterior parietal electrodes bilaterally [6]. A gamma power increase in a parietaloccipital region was also reported in another group of LTM during Vipassana meditation [7], whereas in a recent EEG study meditation experts showed higher parietal-occipital gamma power compared to novices during resting wakefulness preceding meditation practice [8].

An important question raised by these studies is whether enhanced EEG activity is state dependent (e.g.– occurring only during meditation and immediately after it), or trait dependent (e.g. – occurring outside of formal meditation practice) [5]. The latter would suggests that long-term meditation training causes lasting neuro-plastic changes in cortico-thalamic circuits, which should be detected in the spontaneous brain activity [9]. This question is difficult to address because meditation training occurs at multiple time-scales. For example it can happen intensively, as during meditation retreats, or less intensively but across a longer period, as during daily practice. Training also occurs, in different contexts, either during formal practice, or when the practitioner intentionally cultivates meditative qualities within her/his daily activities. Another critical methodological issue is that resting state in meditation experts may be an elusive concept. Indeed, it is possible that during wakefulness a meditator may spontaneously enter a meditative state while “at rest” as a consequence of his/her level of training. However, during sleep brain activity is directed by neither conscious effort nor attention, but rather reflects the intrinsic function of cortico-cortical and cortico-thalamic circuits [10]. Thus, the examination of sleep EEG is an effective means to identify individual trait differences in the brain unconfounded by the effect of meditation on waking mentation. [11].

The recent availability of high-density electroencephalography (hdEEG) has greatly enhanced the spatial resolution of standard EEG recordings, thus allowing for better characterization of local plastic changes [12],[13],[14],[15]. HdEEG studies of sleep rhythms, with good spatial and temporal resolution, are extremely sensitive measures of changes in the activity in corticothalamic circuits. For example, we previously showed that slow wave parameters, such as amplitude and slope, display local changes in the brain in response to learning [16],[17],[12]. Moreover, we have demonstrated that sleep spindle activity is decreased in schizophrenic subjects compared to controls [14]. Finally, gamma power has been shown to increase during NREM sleep in human subjects following declarative learning [18].

Here we performed whole night sleep hdEEG (256-channels) in experienced meditators (LTM) and meditation naive individuals matched for age and sex. Even if this is a correlative and not a causal approach, it paves the way for studies on long term plastic effects of meditation on brain activity in healthy humans.

Materials and Methods

Participants

Twenty-nine right-handed long-term meditators (LTM, mean age = 50.7 6 10.4, 15 female) and a group of twenty-nine meditation naive subjects matched for age and sex were recruited. LTM had a history of daily meditation practice of at least 3 years and had participated in at least 3 one-week intensive retreats. Mean duration of meditation training was 15.6 years (6 7.8, SD). Naive subjects had no previous experience with meditation. All LTM participants were proficient in meditation practices, as taught within the framework of either Theravada or Tibetan Buddhisms. These practices included two attention-based meditations, which we referred to as open monitoring (OM) and focused attention (FA), as well as one compassion/loving kindness meditation referred to as metta meditation [1] (Table 1). Briefly, FA meditation involves directing and sustaining attention on a selected object (e.g., breathing), detecting mind wandering and distractors (e.g., thoughts), as well as disengagement of attention from distractors and shifting of the focus of attention back to the selected object. By contrast, OM meditation has no explicit focus of attention, but rather requires nonreactive meta cognitive monitoring of anything that is experienced, thus replacing the “effortful” selection of an object as primary focus with an “effortless” sustained awareness of the rich features of each experience [1]. The practice of compassion/loving kindness meditation is a form of concentration practice where the practitioner focuses his/her mind on the suffering of oneself or others and then on the wish that the individual(s) in question may be happy and free from suffering. After an initial phone screening to collect the medical and psychiatric history, each subject underwent a thorough in-person screening, which included several questionnaires (see below). Sleep-disordered breathing and sleeprelated movement disorders were also established/excluded with in-laboratory polysomnography (see below). All subjects provided written informed consent and were instructed to maintain regular sleep-wake schedules in the week preceding EEG recordings. This study was approved by the Institutional Review Board of the University of Wisconsin-Madison.

Study design

All subjects underwent in-laboratory hdEEG polysomnography (PSG) that utilized 256 channel hdEEG (Electrical Geodesics Inc., Eugene, OR), as well as standard sleep monitoring leads, including electrooculogram, sub-mental electromyogram, electrocardiogram, bilateral tibial electromyogram, respiratory inductance plethysmography, pulse oximetry, and a position sensor. Participants arrived at the laboratory between 4:00 and 5:00 p.m. for set-up that took approximately two hours. Baseline EEG and an attention and a fear conditioning task were performed. After 9pm the participants were allowed to sleep undisturbed in the laboratory beginning within one hour of their usual bedtime. Additional measures were collected the next day. In this report we will focus only on the EEG data during sleep.

Self reported measures

The socioeconomic status (SES) measure (Socioeconomic status was measured with the Hollingshead Index of Social Position, [19]) was administered to assess the level of education. The Quick Inventory of Depressive Symptoms [20] and the Symptom Checklist-90-Revised [21] were collected to exclude current depression and other mental health issues. The cut-off for exclusion was a score ,2 for both questionnaires. Additionally, each participant completed validated sleep rating scales, including the Insomnia Severity Index (ISI) [22], the Fatigue Severity Scale (FSS) [23], the Epworth Sleepiness Scale [24], a sleep history questionnaire, and the Stanford Sleepiness Scale [25] to assess for symptoms of common sleep disorders, such as restless legs syndrome and obstructive sleep apnea. Thresholds for exclusion were an ISI .10 and/or an FSS .4 and/or an ESS .9.

Meditation Practice

LTM had an average of 8762 lifetime hours of meditation practice, ranging from 1,526 to 32,349 total hours. The computation of lifetime hours of practice was based on subjects’ reports of their average hours of formal (sitting and walking) meditation practice per week, including time spent on meditation retreats (for details see Table 1).

Sleep PSG assessment and hdEEG data analysis

Sleep staging, which was based on six mastoid-referenced channels (F3, F4, C3, C4, O1, and O2), a sub-mental electromyogram and an electrooculogram, was performed by a registered polysomnographic technologist in 30-second epochs according to standard criteria [26] using AliceH Sleepware (Philips Respironics, Murrysville, PA) The sleep technician was blind to group assignment. PSG recordings were reviewed by a board certified sleep medicine physician, who was able to confirm the absence of sleep disorders in 26 out of 29 LTM. Two LTM met the clinical criteria for sleep-related movement disorders (periodic limb movement arousal index .10/h), while one had sleep-disordered breathing (apnea–hypopnea index .10/h). As a result, data from these three LTM (as well as from the age- and sex- matched meditation naives) were not further analyzed. All-night sleep hdEEG recordings were collected with vertex-referencing, using a NetAmps 300 amplifier and NetStation software (Electrical Geodesics Inc., Eugene, OR). EEG data were sampled at 500 Hz, and a first-order high-pass filter (Kaiser type, 0.1 Hz) was applied to eliminate the DC shift. Data were then band-pass filtered (1–50 Hz), down-sampled to 128 Hz and average-referenced to the mean power in all channels. Spectral analysis was performed for each channel in six-second epochs (Welch’s averaged modified periodogram with a Hamming window). For NREM sleep data, a semi-automatic artifact rejection was conducted to remove six-second sleep epochs which exceeded a threshold based on the mean power at either a low (1–4 Hz) or high (20–30 Hz) frequency band [27],[28]. EEG channels in which artifacts affected most of the recording were excluded. REM sleep epochs were visually analyzed and divided in tonic and phasic (i.e., characterized by rapid eye movements). To eliminate artifacts distinctively observed during REM sleep (i.e., muscle twitches, eye movements, heartbeats), independent component analysis (ICA) was performed [29]. After removal of ICA components, power spectral of tonic and phasic REM epochs was computed. We computed six frequency ranges (delta: 1– 4.5 Hz, theta: 4.5–8 Hz, alpha: 8–12 Hz, sigma; 12–15 Hz, beta: 15–25 Hz, and gamma: 25–40 Hz) consistent with previous studies from our lab [30]. In order to determine whether group differences in NREM gamma activity were related to increased muscle tone or variation in saccades, we compared the power in the gamma band in neck EMG and EOG derivations [30] between LTM and meditation naive participants.

Sleep cycles were defined according to the modified criteria [31] of Feinberg and Floyd [32] also consistent with our previous studies [33].

Statistics

Differences in clinical as well as sleep architecture variables were examined using 2-tailed, unpaired t-tests. Topographical analysis was computed after spatially normalizing each subject’s topography (z-score across all channels) within the frequency bands of interest as a means to reduce between-subject variance. Group differences in topographical NREM and REM sleep hdEEG power were assessed with statistical non-parametric mapping (SnPM) using a suprathreshold cluster test to identify significant groups of electrodes after accounting for the multiple comparisons due to the numerous electrodes [34]. Briefly, an appropriate threshold t-value was chosen (t =2, corresponding to a = 0.05 for the given degrees of freedom) before topographic power maps were randomly shuffled between groups (LTM and meditation naive individuals). The size of the largest cluster above the threshold for each reshuffling was then used to create a cluster size distribution. Given the impracticality of computing all possible combinations (4.9661014), 50000 unique combinations were run for each comparison in order to approximate the actual cluster distribution. The suprathreshold cluster p-value was then determined by comparison of the true cluster size against the approximate maximal cluster size distribution. Bonferroni correction of the p value was performed to account for the 18 tests deriving from separately evaluating 6 different frequency bands over the course of 3 distinct phases of sleep.

The magnitude of NREM sleep EEG gamma power differences between LTM and meditation naive individuals was assessed with the Cohen’s d, a measure of the effect size [35]. We also performed correlation analysis between duration of meditation practice in LTM and NREM sleep EEG gamma power, as well as TST and WASO. Statistical analyses were performed using MATLAB (The MathWorks Inc., Natick, MA) and STATISTICA (StatSoft Inc., Tulsa, OK). Non-parametric statistics were used to assess group differences in EMG, EOG, and global EEG gamma power (Wilcoxon rank sum test), as well as to correlate gamma power with variables of interest (Spearman rank correlation). An outlier subject was excluded from the correlation between EEG gamma power and hours of daily practice in LTM based on a suspected over inflation of daily practice (value above 1.5 interquartile ranges from the 75th percentile). Specifically, this LTM had the least amount of retreat time compared with all the LTMs, but the highest amount of non-retreat hours.

Results

Demographics and sleep variables

The LTM and the meditation naive groups were matched for age and sex (see Table 2). Additionally, the groups did not differ in education level, as assessed by the socioeconomic status questionnaire (Table 2). LTM had a significantly reduced total sleep time (TST) as well as an increased wake after sleep onset (WASO) compared to meditation naive subjects (Table 2). These sleep parameters did not show a correlation with daily meditation practice. By contrast, the two groups did not differ in sleep onset latency or in the relative time (% of TST) spent in each sleep stage and did not differ in the number of arousals (Table 2).

LTM had higher NREM sleep gamma power compared to meditation naives in a parietal-occipital brain region We first checked for global, frequency non-specific differences, by comparing the absolute average power across the spectrum between the groups. This analysis revealed no difference between LTM and meditation naive individuals (p = 0.094). We then focused on local, frequency specific effects. Whole night NREM sleep EEG power was topographically compared between LTM and meditation naives by spatially normalizing each subject’s map to the average power in six frequency ranges (delta: 1–4.5 Hz, theta: 4–8 Hz, alpha: 8–12 Hz, sigma; 12–15 Hz, beta: 15– 25 Hz, and gamma: 25–40 Hz). A frequency specific increase in EEG gamma power was found in LTM. Whereas NREM sleep gamma power in both LTM and meditation naives was strongest in frontal/prefrontal areas, and weakest in the temporal regions bilaterally (Figure 1, top row), LTM showed higher relative EEG gamma activity in a parietal-occipital region compared to meditation naives (Figure 1, bottom left). Suprathreshold cluster analysis, a statistical non parametric mapping test (SnPM) which corrects for the multiple comparison problem resulting from the numerous individual electrode tests inherent in hdEEG analysis, confirmed that a parieto-occipital cluster (N= 39 electrodes) showed significantly higher relative gamma power in LTM compared to meditation naive subjects (SnPM, p =0.002, Figure 1, bottom right). This result survived Bonferroni correction for the 18 tests derived from separately evaluating 6 different frequency bands over the course of 3 distinct phases of sleep (p = 0.05/18= 0.003). Therefore, we wanted to check whether this topographical increase in gamma power was paralleled by a difference in the non-normalized gamma power between groups. We found a parieto-occipital cluster (N = 46 electrodes) largely overlapping the normalized cluster, that showed significantly higher absolute gamma power in LTM compared to meditation naive subjects (SnPM, p =0.036). We also examined whether this result was consistent across the night and found similar results in the first three NREM cycles, although the magnitude of the result in the first NREM cycle did not survive multiple comparison corrections after normalization (SnPM, normalized data, cycle 2: N= 20, p =0.012, cycle 3: N= 30, p =0.002; absolute data, cycle 1: N= 37, p = 0.0250, cycle 2: N= 44, p= 0.021, cycle 3: N= 54, p = 0.020). Similar results were obtained when breaking down NREM by stages into N2 and N3, with only N2 surviving multiple comparisons corrections after normalization (SnPM, normalized data, N2: N= 23, p= 0.010; absolute data, N2: N= 49, p = 0.019, N3, N= 51, p =0.014). We suspect that the low values of gamma at the beginning of the night as well as during deep NREM sleep make it difficult to appreciate a difference between groups and therefore the effect is evident but less robust during these times.
We found no gamma power difference in either the electromyographic (EMG) or electrooculographic (EOG) derivations between the two groups (p= 0.589 and p= 0.493 respectively) (Table 2), suggesting that the local gamma increase was not artifactual. SnPM topographic analysis found no significant differences between the two groups for the other five frequency ranges.

LTM NREM sleep gamma power correlated with meditation training

To further characterize the gamma increase at the individual subject level, we calculated the average NREM gamma power of the parietal-occipital cluster for each subject. The average cluster difference in gamma power between groups was 35% (Figure 2, left panel) and the between group Cohen’s d was 0.8, corresponding to a large effect size (more than 50% separation between LTM and meditation naives). We next investigated whether the cluster gamma power was significantly correlated with the overall duration of meditation practice. We found a significant correlation between parietal-occipital NREM gamma power and daily meditation practice, but not retreat time (rho= 0.475, p = 0.017, and rho = 0.029, p =0.887, for daily practice and retreat time respectively). As gamma power in LTM was correlated with age (rho =0.487, p= 0.012), we further tested whether the age of meditation naive individuals predicted also EEG gamma power, but no correlation was found (rho =20.035, p = 0.866). For each group there was also no correlation between gamma power in the cluster and the significantly different sleep architecture variables (TST, rho =20.195, p = 0.339, rho =0.287, p = 0.156; WASO, rho = 0.170, p= 0.408, rho =20.319, p = 0.112, for LTM and meditation naive individuals, respectively). Thus, daily practice was the most sensitive predictor of the correlation with gamma power activity.

REM sleep gamma power did not differ between LTM and naives

To assess whether the gamma increase of LTM was specific for NREM sleep, topographic EEG analysis of REM sleep was also performed. As previous literature has suggested functional differences between tonic and phasic REM, we separated these two REM sleep patterns in our analysis [36],[37],[38]. Both tonic and phasic REM gamma power peaked in a frontal/prefrontal as well as in a parieto-occipital region in LTM and meditation naives (Figure S1). SNPM analyses revealed no topographical differences in any frequency band in REM sleep between the two groups.

Discussion

By performing sleep hdEEG (256 channels) recordings in this study we found that LTM had increased parietal-occipital EEG gamma power during NREM sleep compared to meditation naives. This increase was specific for the gamma frequency range, was not found during REM sleep, and was positively correlated with the length of daily meditation practice.

The LTMs recorded in this study were experienced in Focused Attention (FA), Open Monitoring (OM), as well as loving kindness / compassion meditations. An increase in gamma activity was recently demonstrated in expert Buddhist practitioners (.10.000 hours of practice) during a style of meditation which contains features of both loving kindness/compassion and OM meditations [39]. Compared to a group of novices, the expert practitioners showed self-induced, higher-amplitude, sustained EEG gamma-band oscillations, especially over lateral frontalparietal electrodes, while meditating as well as in the resting state immediately preceding and following meditation [6]. Notably, a link between higher gamma-band activity and stronger cognitive control has been reported by a variety of human electrophysiological techniques, including magnetoencephalography (MEG) [40], scalp EEG [41], and direct cortical recordings [42]. Furthermore, a recent EEG study has found that LTM had increased gamma power in a parietal-occipital area compared to meditation naive individuals during resting state, as well as during meditation practice as compared to baseline [8]. Based on these findings, the authors concluded that enhanced posterior EEG gamma power was a state (meditation-related), and to some same extent trait (resting state associated) feature of meditation practice.

Here we established that LTM had higher EEG gamma power over an extended period of spontaneous brain activity (i.e., whole night NREM sleep) compared to meditation naive individuals. This parietal-occipital gamma increase corresponded to a large effect size (ES =0.8, allowing more than 50% separation between meditation experts and naives). These findings strongly suggest that changes in EEG gamma activity related to meditation practice are trait (in addition to state) related. Comparing the activity of meditation experts and novices during an unambiguous resting condition is a challenging task [5]. This is because experienced long-term meditation experts are usually able to blend formal meditation session with daily life, and that a meditator might spontaneously generate a meditative state while at rest in the lab out of demand characteristic. In this regard, sleep provides an exquisite window to explore the spontaneous brain activity as well as the function of neural circuits at rest. Brain activity during sleep does not require conscious effort or attention, and a condition of physical immobility is obtained for several hours. Moreover, the recently available combination of standard polysomnography with hdEEG, which provides enhanced spatial and temporal resolution, offers the opportunity to analyze in greater details NREM/REM sleep activity as well as to observe local changes in brain function due to neuroplasticity [43],[16].

Daily practice and meditation retreat could contribute differently to the neuroplastic changes induced by meditation. For instance, meditation frequency (days per week with meditation practice) has recently been shown to reliably predict both higher mindfulness and psychological well-being [44]. In this study we showed that the daily practice, but not the retreat time, predicted the parietal-occipital gamma activity during NREM sleep in LTM. The differential effect of the amount of daily practice and of retreat time on localized NREM gamma activity is to our knowledge the first indication of a specific effect of constant meditation daily practice, but not of intensive retreat practice, on brain neuroplasticity. On one side, this finding raises some methodological issues about how to quantify meditation practice, suggesting the potential usefulness of differentially investigating the contribution of retreat time and daily practice on behavioral and physiological measures. On the other side, it enhances the effectiveness of our approach in revealing stable (trait-like) effects on brain functioning induced by prolonged training during waking.

In the present study LTM had a significantly reduced total sleep time (TST) and increased wake after sleep onset (WASO) compared to meditation naive individuals. A reduction in TST has been recently reported by another sleep study in LTM [45], and it indicates that meditation practice may decrease sleep needs. However other studies [46,47] investigating sleep architecture didn’t find a reduction of TST, suggesting that the sleep architecture is not the most reliable parameter to study the effect of meditation on neuronal plasticity during sleep. Consistent with this idea, we did not find a correlation between the changes in any traditional polysomnographic sleep parameters and meditation practice.

Only a handful of studies so far have investigated the sleep EEG activity of meditation experts beyond sleep architecture. One of these studies explored EEG differences in thirteen individuals with at least 2 years of meditation experience during Transcendental meditation (TM; a form of meditation different from that explored in the current study), resting wakefulness, drowsiness, and sleep and found a progressive slowing of the main EEG frequency from wakefulness to sleep, with no appreciable change in power between wakefulness and meditation EEG [48]. The authors also analyzed the EEG activity of meditation practitioners and of meditation naive control subjects during resting wakefulness, and found no difference in power but a slight slowing in the mean EEG frequency of the practitioners; however, they did not compare the sleep EEG of these two groups [48]. Another study found that eleven long term TM practitioners had increased theta-alpha power during slow wave sleep compared to nine short term practitioners as well as eleven experienced practitioners [49]. Here we found no difference in theta-alpha EEG power between LTM and meditation naive individuals during NREM sleep. Differences in style of meditation practices may account for the discrepancy of these findings. Furthermore, whereas we screened participants for sleep disorders and performed PSG recordings during the first hdEEG night, the LTM recorded by Mason et al. did not undergo such screening. An EEG pattern of alpha wave in delta wave sleep (alpha-delta sleep) is commonly reported in individuals with sleep disturbances, including restless leg syndrome [50] and sleep apnea [51], and is associated with increased arousability and lighter sleep [52]. Notably, their LTM spent a significantly higher amount of time in light (N1) NREM sleep, while we found no difference in N1 between LTM and meditation naive individuals. [45]And finally, the type of meditation practice examined in the current study differs from the practice studied in these other studies.

What is the functional meaning of the gamma increase in LTM that was found here during NREM sleep? A large body of evidence from animal and human recordings have suggested that gamma-frequency activity is implicated during wakefulness in plasticity-related processes, including attention, learning, as well as both working and long-term memory [53]. For instance, an increase in gamma activity occurs when sensory stimuli are attended [54], as well as during the active maintenance of representations during working memory tasks [55],[56]. Several studies employing EEG, MEG and intracranial EEG recordings have also shown that gamma-frequency activity during encoding predicts successful formation of long-term memory [57],[58].
Little is still known about the functional significance of gamma activity during sleep. In this study we found a gamma power increase in meditation experts during NREM sleep in a scalp region overlying posterior parietal and occipital cortical areas. Our finding is in line with the recent report of Valderrama et al., which suggested a functional significance for gamma oscillations during NREM sleep in humans [59]. A possible confound in the study of gamma activity is the presence of muscular or ocular artifacts in the scalp EEG. However, we found no difference in neck EMG and EOG gamma power between LTM and meditation naïve individuals. This result therefore strongly indicates that such artifacts do not contribute significantly to our finding.
Previous studies have reported an increase in fast-frequency activity during sleep in pathological conditions such as schizophrenia, depression and insomnia [60],[61],[62]. This increase has been associated with abnormalities in the arousal mechanisms [62]. However, this functional interpretation of gamma activity is not warranted for a number of reasons. First, we did not find a group difference in arousal. Our finding was specific to spatially normalized gamma activity. Furthermore, the topography of the current finding over parietal-occipital electrodes differed from those earlier findings. Finally inter-individual variability in gamma activity among meditators predicted daily meditation experience in life, which has been found to predict positive mental health outcomes [44],[1]. Instead, we speculate here the specific gamma increase is not pathological, but reflects the lasting, plastic effect on specific neuronal circuits of long-term meditation practice. For instance, the parietal cortex has been implicated in directing the focus of attention on a specific object [63], a cognitive skill perfected by attention-based meditation [1]. The finding of higher whole night NREM gamma power in LTM compared to meditation naives could therefore reflect an increase in activity and connectivity within these neuronal circuits due to extensive meditation training. Consistent with this assumption, here we established that NREM sleep gamma band power correlated with the duration of daily meditation practice in LTM. This finding not only confirms a previous report from our group that the amount of EEG gamma power generated during meditation by LTM correlated with the length of their meditation training [6], but also suggests that EEG gamma activity during sleep may represent a trait-like sensitive measure for the neuronal plastic changes determined over time by meditative training.

The increase in NREM sleep parieto-occipital gamma power reported here could also reflect an enhanced activity of the underlying cortical areas in LTM compared to meditation naives. Higher EEG gamma power reflects higher firing rates of the underlying neuronal populations [64], and local changes in gamma oscillations closely mirror underlying activity in both visual [65],[66] and parietal default network-associated cortical regions [67],[68]. During NREM sleep both neuronal firing and gamma power tend to decrease, as does the ability to process sensory information [69],[70]. Thus, a higher gamma activity in LTM could reflect a partially maintained capacity of parietooccipital sensory and default network-associated areas to process information and maintain some level of awareness, even during a state when usually these cognitive functions are greatly impaired. Consistent with this idea, a higher incidence of dream reports has been found in meditation experts compared to meditation naives even during the deepest stages of NREM sleep [71]. If experienced meditators retain a higher capacity for internal information processing and awareness during NREM sleep compared to meditation naives, such advantage should be reduced during REM sleep, when these functions are partially restored and spontaneous neuronal firing/gamma activity is enhanced compared to NREM sleep [72]. In this study we found that LTM had only a slight, non significant increase in REM sleep EEG gamma power compared to meditation naives in the same parietal-occipital regions (Figure S1).

Limitations of the study include a lack of an adaptation night, which could account for the truncated sleep time (, 7 hrs) in all participants, but it is unlikely to explain the observed group difference in the sleep EEG. Future work will also need to address some of the questions left unanswered in the present study. For example, the relationship found here between higher EEG gamma activity and longer meditation daily practice suggests that gamma power is a good correlate of meditation training. This correlation should be confirmed in longitudinal studies performing EEG recordings in meditation naive individuals before and after meditation training, ideally using only one style of meditation practice. It will also be important to investigate whether the observed gamma increase may be affected by pre-existing “baseline” gamma activity differences between groups (i.e., meditation experts and naives), as previously suggested [6]. Gamma activity has been shown to be influenced by several factors, including age [73], sex [74], and cognitive abilities [75]. However, these factors are unlikely to have contributed to the present findings, given that LTM and meditation naive subjects were matched for age, sex, and did not differ in education level. Future studies should investigate whether the group difference in NREM gamma activity in meditators is associated to a specific meditation practice (e.g. mindfulness meditation vs. compassion meditation) or style of meditation training (e.g. Tibetan Buddhism vs. Theravada Buddhism). Specifically, investigating the acute effect of an intense meditation session in LTM on sleep EEG patterns could help in establishing a causal relationship between meditation training and specific changes in EEG activity. Finally, future experiments combining fMRI with simultaneous hdEEG will be critical to fully characterize the cortical (and possibly subcortical) networks underlying the enhanced NREM sleep EEG gamma activity found in this study in meditation experts, whereas studies investigating the healing effects of meditation interventions could explore the ability of EEG gamma power to predict a beneficial effect of such interventions. This work would contribute to identify the neural circuits underlying the EEG correlates of meditation training. It will also help to establish whether EEG gamma activity represents a sensitive and objective measure of the effects of meditative practice on brain function in both healthy subjects and brain disordered patients.

Supporting Information
Figure S1 REM tonic as well as phasic gamma power did not differ between LTM and meditation naives.

Topographic color plots showed maximal REM tonic as well as phasic gamma power in frontal/prefrontal regions in both groups. Compared to meditation naives, LTM had a slightly higher power in the same parieto-occipital region significantly more active during NREM sleep, which however failed to reach significance in both tonic and phasic REM (white topographic plots, p =0.975, and p = 0.810, SnPM). Only the tonic REM topographies are shown. (TIF)

Acknowledgments
We would like to thank Me´lanie Boly for providing helpful comments on this manuscript and for helpful suggestion on EEG analysis and Donal MacCoon and Brianna Schuyler for helping analyzing the self-report questionnaires in this study.

Author Contributions
Conceived and designed the experiments: FF AL RJD GT. Performed the experiments: FF RS DD BR CZ. Analyzed the data: FF RS DD CZ RB. Contributed reagents/materials/analysis tools: FF BR RS DD. Wrote the paper: FF BR AL RJD DD GT.

References

1. Lutz A, Slagter HA, Dunne JD, Davidson RJ (2008) Attention regulation and monitoring in meditation. Trends in Cognitive Sciences 12: 163–169.
2. Chambers R, Gullone E, Allen NB (2009) Mindful emotion regulation: An integrative review. Clin Psychol Rev 29: 560–572.
3. Holzel BK, Carmody J, Vangel M, Congleton C, Yerramsetti SM, et al. (2011) Mindfulness practice leads to increases in regional brain gray matter density. Psychiatry Res 191: 36–43.
4. Lutz A, Slagter HA, Rawlings NB, Francis AD, Greischar LL, et al. (2009) Mental training enhances attentional stability: neural and behavioral evidence.
The Journal of Neuroscience: The Official Journal of the Society for Neuroscience 29: 13418–13427.
5. Cahn BR, Polich J (2006) Meditation states and traits: EEG, ERP, and neuroimaging studies. Psychological Bulletin 132: 180–211.
6. Lutz A, Greischar LL, Rawlings NB, Ricard M, Davidson RJ (2004) Long-term meditators self-induce high-amplitude gamma synchrony during mental
practice. Proceedings of the National Academy of Sciences of the United States of America 101: 16369–16373.
7. Cahn BR, Delorme A, Polich J (2010) Occipital gamma activation during Vipassana meditation. Cognitive Processing 11: 39–56.
8. Berkovich-Ohana A, Glicksohn J, Goldstein A (2012) Mindfulness-induced changes in gamma band activity – Implications for the default mode network, self-reference and attention. Clinical Neurophysiology: Official Journal of the International Federation of Clinical Neurophysiology 123: 700–710.
9. Davidson RJ, Lutz A (2008) Buddha’s Brain: Neuroplasticity and Meditation. IEEE Signal Process Mag 25: 176–174.
10. Steriade M (2003) The corticothalamic system in sleep. Front Biosci 8: d878– 899.
11. Finelli LA, Achermann P, Borbely AA (2001) Individual ’fingerprints’ in human sleep EEG topography. Neuropsychopharmacology 25: S57–62.
12. Huber R, Ghilardi MF, Massimini M, Tononi G (2004) Local sleep and learning. Nature 430: 78–81.
13. Huber R, Ghilardi MF, Massimini M, Ferrarelli F, Riedner BA, et al. (2006) Arm immobilization causes cortical plastic changes and locally decreases sleep slow wave activity. Nat Neurosci 9: 1169–1176.
14. Ferrarelli F, Huber R, Peterson MJ, Massimini M, Murphy M, et al. (2007) Reduced sleep spindle activity in schizophrenia patients. Am J Psychiatry 164: 483–492.
15. Ferrarelli F, Peterson MJ, Sarasso S, Riedner BA, Murphy MJ, et al. (2010) Thalamic dysfunction in schizophrenia suggested by whole-night deficits in slow and fast spindles. Am J Psychiatry 167: 1339–1348.
16. Sarasso S, Santhanam P, Maatta S, Poryazova R, Ferrarelli F, et al. (2010) Nonfluent aphasia and neural reorganization after speech therapy: insights from human sleep electrophysiology and functional magnetic resonance imaging. Arch Ital Biol 148: 271–278.
17. Huber R, Esser SK, Ferrarelli F, Massimini M, Peterson MJ, et al. (2007) TMSinduced cortical potentiation during wakefulness locally increases slow wave activity during sleep. PLoS ONE 2: e276.
18. Molle M, Marshall L, Gais S, Born J (2004) Learning increases human electroencephalographic coherence during subsequent slow sleep oscillations. Proc Natl Acad Sci U S A 101: 13963–13968.
19. Hollingshead AB (1957) Two Factor Index of Social Position. Mimeo. New Haven, Connecticut: Yale University.
20. Rush AJ, Trivedi MH, Ibrahim HM, Carmody TJ, Arnow B, et al. (2003) The 16-Item Quick Inventory of Depressive Symptomatology (QIDS), clinician rating (QIDS-C), and self-report (QIDS-SR): a psychometric evaluation in patients with chronic major depression. Biol Psychiatry 54: 573–583.
21. Derogatis LR, Cleary PA (1977) Factorial invariance across gender for the primary symptom dimensions of the SCL-90. Br J Soc Clin Psychol 16: 347–
356.
22. Bastien CH, Vallieres A, Morin CM (2001) Validation of the Insomnia Severity Index as an outcome measure for insomnia research. Sleep Med 2: 297–307.
23. Valko PO, Bassetti CL, Bloch KE, Held U, Baumann CR (2008) Validation of the fatigue severity scale in a Swiss cohort. Sleep 31: 1601–1607.
24. Johns MW (1991) A new method for measuring daytime sleepiness: the Epworth sleepiness scale. Sleep 14: 540–545.
25. MacLean AW, Fekken GC, Saskin P, Knowles JB (1992) Psychometric evaluation of the Stanford Sleepiness Scale. J Sleep Res 1: 35–39.
26. Silber MH, Ancoli-Israel S, Bonnet MH, Chokroverty S, Grigg-Damberger MM, et al. (2007) The visual scoring of sleep in adults. J Clin Sleep Med 3: 121–
131.
27. Landsness E, Bruno MA, Noirhomme Q, Riedner B, Gosseries O, et al. (2011) Electrophysiological correlates of behavioural changes in vigilance in vegetative state and minimally conscious state. Brain 134: 2222–2232.
28. Murphy M, Riedner BA, Huber R, Massimini M, Ferrarelli F, et al. (2009) Source modeling sleep slow waves. Proc Natl Acad Sci U S A 106: 1608–1613.
29. Jung TP, Makeig S, McKeown MJ, Bell AJ, Lee TW, et al. (2001) Imaging Brain Dynamics Using Independent Component Analysis. Proc IEEE Inst Electr
Electron Eng 89: 1107–1122.
30. Murphy M, Bruno MA, Riedner BA, Boveroux P, Noirhomme Q, et al. (2011) Propofol anesthesia and sleep: a high-density EEG study. Sleep 34: 283–291A.
31. Aeschbach D, Borbely AA (1993) All-night dynamics of the human sleep EEG. J Sleep Res 2: 70–81.
32. Feinberg I, Floyd TC (1979) Systematic trends across the night in human sleep cycles. Psychophysiology 16: 283–291.
33. Riedner BA, Vyazovskiy VV, Huber R, Massimini M, Esser S, et al. (2007) Sleep homeostasis and cortical synchronization: III. A high-density EEG study of sleep slow waves in humans. Sleep 30: 1643–1657.
34. Nichols TE, Holmes AP (2002) Nonparametric permutation tests for functional neuroimaging: a primer with examples. Hum Brain Mapp 15: 1–25.
35. Cohen J (1988) Statistical Power Ananlysis for the Behavioral Sciences. Hillsdale, NJ: Erlbaum.
36. Corsi-Cabrera M, Guevara MA, del Rio-Portilla Y (2008) Brain activity and temporal coupling related to eye movements during REM sleep: EEG and MEG results. Brain Res 1235: 82–91.
37. Ermis U, Krakow K, Voss U (2010) Arousal thresholds during human tonic and phasic REM sleep. J Sleep Res 19: 400–406.
38. Takahara M, Nittono H, Hori T (2002) Comparison of the event-related potentials between tonic and phasic periods of rapid eye movement sleep. Psychiatry Clin Neurosci 56: 257–258.
39. Lutz A, Greischar LL, Rawlings NB, Ricard M, Davidson RJ (2004) Long-term meditators self-induce high-amplitude gamma synchrony during mental practice. Proc Natl Acad Sci U S A 101: 16369–16373.
40. Kinsey K, Anderson SJ, Hadjipapas A, Holliday IE (2011) The role of oscillatory brain activity in object processing and figure-ground segmentation in human vision. Int J Psychophysiol 79: 392–400.
41. Schadow J, Lenz D, Dettler N, Frund I, Herrmann CS (2009) Early gammaband responses reflect anticipatory top-down modulation in the auditory cortex.
NeuroImage 47: 651–658.
42. Engell AD, McCarthy G (2010) Selective attention modulates face-specific induced gamma oscillations recorded from ventral occipitotemporal cortex.
J Neurosci 30: 8780–8786.
43. Slagter HA, Lutz A, Greischar LL, Francis AD, Nieuwenhuis S, et al. (2007) Mental training affects distribution of limited brain resources. PLoS Biology 5: e138.
44. Schoormans D, Nyklicek I (2011) Mindfulness and psychologic well-being: are they related to type of meditation technique practiced? J Altern Complement Med 17: 629–634.
45. Kaul P, Passafiume J, Sargent CR, O9Hara BF (2010) Meditation acutely improves psychomotor vigilance, and may decrease sleep need. Behavioral and Brain Functions: BBF 6: 47.
46. Pattanashetty R, Sathiamma S, Talakkad S, Nityananda P, Trichur R, et al. (2010) Practitioners of vipassana meditation exhibit enhanced slow wave sleep and REM sleep states across different age groups. Sleep and Biological Rhythms 8: 34–41.
47. Nagendra RP, Maruthai N, Kutty BM (2012) Meditation and its regulatory role on sleep. Frontiers in Neurology 3: 54.
48. Stigsby B, Rodenberg JC, Moth HB (1981) Electroencephalographic findings during mantra meditation (transcendental meditation). A controlled, quantitative study of experienced meditators. Electroencephalogr Clin Neurophysiol 51: 434–442.
49. Mason LI, Alexander CN, Travis FT, Marsh G, Orme-Johnson DW, et al. (1997) Electrophysiological correlates of higher states of consciousness during sleep in long-term practitioners of the Transcendental Meditation program. Sleep 20: 102–110.
50. Saskin P, Moldofsky H, Lue FA (1985) Periodic movements in sleep and sleepwake complaint. Sleep 8: 319–324.
51. Gold AR, Dipalo F, Gold MS, O9Hearn D (2003) The symptoms and signs of upper airway resistance syndrome: a link to the functional somatic syndromes. Chest 123: 87–95.
52. Martinez D, Breitenbach TC, Lenz Mdo C (2010) Light sleep and sleep time misperception – relationship to alpha-delta sleep. Clin Neurophysiol 121: 704– 711.
53. Jensen O, Kaiser J, Lachaux JP (2007) Human gamma-frequency oscillations associated with attention and memory. Trends Neurosci 30: 317–324.
54. Siegel M, Konig P (2003) A functional gamma-band defined by stimulus dependent synchronization in area 18 of awake behaving cats. J Neurosci 23: 4251–4260.
55. Jokisch D, Jensen O (2007) Modulation of gamma and alpha activity during a working memory task engaging the dorsal or ventral stream. J Neurosci 27: 3244–3251.
56. Mainy N, Kahane P, Minotti L, Hoffmann D, Bertrand O, et al. (2007) Neural correlates of consolidation in working memory. Hum Brain Mapp 28: 183–193.
57. Osipova D, Takashima A, Oostenveld R, Fernandez G, Maris E, et al. (2006) Theta and gamma oscillations predict encoding and retrieval of declarative memory. J Neurosci 26: 7523–7531.
58. Sederberg PB, Miller JF, Howard MW, Kahana MJ (2010) The temporal contiguity effect predicts episodic memory performance. Mem Cognit 38: 689– 699.
59. Valderrama M, Cre´pon B, Botella-Soler V, Martinerie J, Hasboun D, et al. (2012) Human Gamma Oscillations during Slow Wave Sleep. PloS One 7: e33477.
60. Tekell JL, Hoffmann R, Hendrickse W, Greene RW, Rush AJ, et al. (2005) High frequency EEG activity during sleep: characteristics in schizophrenia and depression. Clin EEG Neurosci 36: 25–35.
61. Perlis ML, Merica H, Smith MT, Giles DE (2001) Beta EEG activity and insomnia. Sleep Med Rev 5: 363–374.
62. Feige B, Scaal S, Hornyak M, Gann H, Riemann D (2007) Sleep electroencephalographic spectral power after withdrawal from alcohol in alcohol-dependent patients. Alcohol Clin Exp Res 31: 19–27.
63. Colby CL, Goldberg ME (1999) Space and attention in parietal cortex. Annu Rev Neurosci 22: 319–349.
64. Nir Y, Fisch L, Mukamel R, Gelbard-Sagiv H, Arieli A, et al. (2007) Coupling between neuronal firing rate, gamma LFP, and BOLD fMRI is related to interneuronal correlations. Curr Biol 17: 1275–1285.
65. Logothetis NK, Pauls J, Augath M, Trinath T, Oeltermann A (2001) Neurophysiological investigation of the basis of the fMRI signal. Nature 412: 150–157.
66. Shmuel A, Leopold DA (2008) Neuronal correlates of spontaneous fluctuations in fMRI signals in monkey visual cortex: Implications for functional connectivity at rest. Hum Brain Mapp 29: 751–761.
67. Miller KJ, Weaver KE, Ojemann JG (2009) Direct electrophysiological measurement of human default network areas. Proc Natl Acad Sci U S A 106: 12174 12177.
68. Ossandon T, Jerbi K, Vidal JR, Bayle DJ, Henaff MA, et al. (2011) Transient suppression of broadband gamma power in the default-mode network is correlated with task complexity and subject performance. J Neurosci 31: 14521–14530.
69. Born AP, Law I, Lund TE, Rostrup E, Hanson LG, et al. (2002) Cortical deactivation induced by visual stimulation in human slow-wave sleep. Neuro- Image 17: 1325–1335.
70. Czisch M, Wehrle R, Stiegler A, Peters H, Andrade K, et al. (2009) Acoustic oddball during NREM sleep: a combined EEG/fMRI study. PLoS ONE 4: e6749.
71. Faber PA, Saayman GS, Touyz SW (1978) Meditation and archetypal content of nocturnal dreams. J Anal Psychol 23: 1–22.
72. Massimini M, Ferrarelli F, Murphy M, Huber R, Riedner B, et al. (2010)
Cortical reactivity and effective connectivity during REM sleep in humans. Cogn Neurosci 1: 176–183.
73. Lu CB, Hamilton JB, Powell AD, Toescu EC, Vreugdenhil M (2011) Effect of ageing on CA3 interneuron sAHP and gamma oscillations is activity-dependent. Neurobiol Aging 32: 956–965.
74. Guntekin B, Basar E (2007) Brain oscillations are highly influenced by gender differences. Int J Psychophysiol 65: 294–299.
75. Jausovec N, Jausovec K (2007) Personality, gender and brain oscillations. Int J Psychophysiol 66: 215–224.

Copyright: _ 2013 Ferrarelli et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.