Previous Article |
Table of Contents
| Next Article
NEUROSCIENCE
Long-term meditators self-induce high-amplitude gamma synchrony during mental practice



*W. M. Keck Laboratory for Functional Brain Imaging and Behavior, Waisman Center, and Laboratory for Affective Neuroscience, Department of Psychology, University of Wisconsin, 1500 Highland Avenue, Madison, WI 53705; and
Shechen Monastery, P.O. Box 136, Kathmandu, Nepal
Communicated by Burton H. Singer, Princeton University, Princeton, NJ, October 6, 2004 (received for review August 26, 2004)
| Abstract |
|---|
|
|
|---|
electroencephalogram synchrony | gamma activity | meditation
| Methods |
|---|
|
|
|---|
We first collected an initial electroencephalogram (EEG) baseline consisting of four 60-s blocks of ongoing activity with a balanced random ordering of eyes open or closed for each block. Then, subjects generated three meditative states, only one of which will be described in this report. During each meditative session, a 30-s block of resting activity and a 60-s block of meditation were collected four times sequentially. The subjects were verbally instructed to begin the meditation and meditated at least 20 s before the start of the meditation block. We focus here on the last objectless meditative practice during which both the controls and Buddhist practitioners generated a state of "unconditional loving-kindness and compassion."
Meditative Instruction. The state of unconditional loving-kindness and compassion is described as an "unrestricted readiness and availability to help living beings." This practice does not require concentration on particular objects, memories, or images, although in other meditations that are also part of their long-term training, practitioners focus on particular persons or groups of beings. Because "benevolence and compassion pervades the mind as a way of being," this state is called "pure compassion" or "nonreferential compassion" (dmigs med snying rje in Tibetan). A week before the collection of the data, meditative instructions were given to the control subjects, who were asked to practice daily for 1 h. The quality of their training was verbally assessed before EEG collection. During the training session, the control subjects were asked to think of someone they care about, such as their parents or beloved, and to let their mind be invaded by a feeling of love or compassion (by imagining a sad situation and wishing freedom from suffering and well being for those involved) toward these persons. After some training, the subjects were asked to generate such feeling toward all sentient beings without thinking specifically about anyone in particular. During the EEG data collection period, both controls and long-term practitioners tried to generate this nonreferential state of loving-kindness and compassion. During the neutral states, all of the subjects were asked to be in a nonmeditative, relaxed state.
EEG Recordings and Protocol. EEG data were recorded at standard extended 10/20 positions with a 128-channel Geodesic Sensor Net (Electrical Geodesics, Eugene, OR), sampled at 500 Hz, and referenced to the vertex (Cz) with analog band-pass filtering between 0.1 and 200 Hz. EEG signals showing eye movements or muscular artifacts were manually excluded from the study. A digital notch filter was applied to the data at 60 Hz to remove any artifacts caused by alternating current line noise.
Bad channels were replaced by using spherical spline interpolation (12). Two-second epochs without artifact were extracted after the digital rereferencing to the average reference.
Spectral Analysis. For each electrode and for each 2-s epoch, the power spectral distribution was computed by using Welch's method (13), which averages power values across sliding and overlapping 512-ms time windows. To compute the relative gamma activity, the power spectral distribution was computed on the z-transformed EEG by using the mean and SD of the signal in each 2-s window. This distribution was averaged through all electrodes, and the ratio between gamma and slow rhythms was computed. Intraindividual analyses were run on this measure and a group analysis was run on the average ratio across 2-s windows. The group analysis of the topography was performed by averaging the power spectral distribution for each electrode in each block and then computing the ratio of gamma to slow rhythms before averaging across blocks.
Despite careful visual examination, the electroencephalographic spectral analysis was hampered by the possible contamination of brain signals by muscle activity. Here we assume that the spectral emission between 80 and 120 Hz provided an adequate measure of the muscle activity (14, 15). The muscle EEG signature is characterized by a broad-band spectrum profile (8-150 Hz) peaking at 70-80 Hz (16). Thus, the variation through time of the average spectral power in the 80-120 Hz frequency band provided a way to quantify the variations of the muscle contribution to the EEG gamma activity through time. To estimate the gamma activity, adjusted for the very high frequencies, we performed a covariance analysis for each region of interest (ROI) for each subject. The dependent variable was the average gamma activity (25-42 Hz) in each ROI. The continuous predictor was the electromyogram activity (80-120 Hz power). The categorical predictors were the blocks (initial baseline with eyes open and neutral blocks from 2 to 4) and the mental states (ongoing neutral versus meditation).
For the group analysis, separate repeated ANOVAs were then performed on the relative gamma and adjusted gamma variation between states, with the blocks as the within factor and the group (practitioners versus controls) as the categorical predictor. For the intrasubject analysis, we compared separately the relative gamma and the raw gamma activity averaged within the ROIs in the initial baseline state versus the meditative state.
Phase-Synchrony Detection. Electrodes of interest were referenced to a local average potential defined as the average potential of its six surrounding neighbors. This referencing montage restricted the electrical measurement to local sources only and prevented spurious long-range synchrony from being detected if the muscle activity over one electrode propagated to another distant electrode. The methods used to measure long-range synchronization are described in detail in Supporting Methods, which is published as supporting information on the PNAS web site. In summary, for each epoch and electrode, the instantaneous phase of the signal was extracted at each frequency band between 25 and 42 Hz in 2-Hz steps by using a convolution with Morlet wavelets. The stability through time of their phase difference was quantified in comparison with white-noise signals as independent surrogates. A measure of synchronous activity was defined as the number of electrode pairs among the 294 studied combinations that had higher synchrony density on average across frequencies than would be expected to occur between independent signals. The electrode pairs were taken between the ROIs when we measured the scalp distribution of gamma activity (see Fig. 3a). A repeated-measures ANOVA was performed on the average size of the synchrony pattern across all frequency bands and epochs in each block with the original resting state and the meditative state as the within factors and the group (practitioners versus controls) as the between-groups factor.
|
| Results |
|---|
|
|
|---|
|
|
These data suggest that the two groups had different electrophysiological spectral profiles in the baseline, which are characterized by a higher ratio of gamma-band oscillatory rhythm to slow oscillatory rhythms for the long-term practitioners than for the controls. This group difference is enhanced during the meditative practice and continues into the postmeditative resting blocks.
Absolute Gamma Power. We then studied the variation through time of the ongoing gamma-band activity itself. The gamma-band activity (25-42 Hz) was first z-transformed in each block and compared over each electrode with the mean and SD of their respective neutral block (ongoing baseline). The normalized gamma activity was then averaged across the blocks. Fig. 3a shows the percentage of subjects presenting an increase of at least 1 SD during meditation compared with neutral state. A common topographical pattern of gamma activity emerged across the long-term practitioners but not across the control subjects. This pattern was located bilaterally over the parieto-temporal and midfrontal electrodes. Fig. 3a shows four ROIs containing seven electrodes each and located around F3-8, Fc3-6, T7-8, Tp7-10, and P7-10. Hereafter, we focus on the electrodes activated in these ROIs.
Intraindividual analyses similar to those for relative gamma activity were run on the average gamma power across these ROIs and exhibited the same pattern as that found for relative gamma. It is possible that these high-amplitude oscillations are partially contaminated by muscle activity (18). Because we found increases in gamma activity during the postmeditative resting baseline compared with the initial resting baseline, it is unlikely that the changes we reported could be solely caused by muscle activity, because there was little evidence of any muscle activity during these baseline periods. (Fig. 2e). Secondly, we showed that the meditative state and nonmeditative state that mimicked and exaggerated the possible muscle activity during meditation exhibit significantly different spectral profiles (Fig. 4, which is published as supporting information on the PNAS web site). Furthermore, for the two subjects showing the highest gamma activity, we showed that amplitude of the gamma-band activity before external stimulation predicts the amplitude of high fast-frequency oscillations (20-45 Hz) evoked by auditory stimuli (Fig. 5, which is published as supporting information on the PNAS web site). Because the evoked activity is relatively independent of muscle activity, the relationship between the pre-stimulation fast-frequency oscillation and the evoked activity suggests that these high-amplitude gamma rhythms are not muscle artifacts (Fig. 5 and Fig. 6, which is published as supporting information on the PNAS web site). This claim is further supported by the localization within the brain of the dipole sources of these fast-frequency-evoked oscillations (Figs. 7-9, which are published as supporting information on the PNAS web site).
Yet we still chose to cautiously interpret the raw values of these gamma oscillations because of the concomitant increase of spectral power >80 Hz during meditation. This increase could also reflect a change in muscle activity rather than high-frequency, gamma-band oscillations [70-105 Hz (19)], which are mostly low-pass filtered by the skull at >80 Hz. Thus, we chose to conservatively interpret the activity at >80 Hz as indicating muscle activity.
To remove the contribution of putative muscle activity, we quantified the increase in the average amplitude of gamma oscillation (25-42 Hz) adjusted for the effect of the very high-frequency variation (80-120 Hz) (see Methods and ref. 20). The adjusted average variation in gamma activity was >30-fold greater among practitioners compared with controls (Fig. 3b). Group analysis was run on the average adjusted gamma activity over these ROIs. Gamma activity increased for both the long-term practitioners and controls from neutral to meditation states [F(1, 16) = 5.2, P < 0.05; ANOVA], yet this increase was higher for the long-time practitioners than for the controls [F(1, 16) = 4.6, P < 0.05; interaction between the state and group factors ANOVA] (Fig. 3b). In summary, the generation of this meditative state was associated with gamma oscillations that were significantly higher in amplitude for the group of practitioners than for the group of control subjects.
Long-Distance Gamma Synchrony. Finally, a long-distance synchrony analysis was conducted between electrodes from the ROIs found in Fig. 3a. Long-distance synchrony is thought to reflect large-scale neural coordination (9) and can occur when two neural populations recorded by two distant electrodes oscillate with a precise phase relationship that remains constant during a certain number of oscillation cycles. This approach is illustrated in Fig. 1c for selected electrodes (F3/4, Fc5/6, and Cp5/6). For subject S4, the density of cross-hemisphere, long-distance synchrony increases by
30% on average during meditation and follows a pattern similar to the oscillatory gamma activity.
For all subjects, locally referenced, long-distance synchronies were computed for each 2-s epoch during the neutral and meditative states between all electrode pairs and across eight frequencies ranging from 25 to 42 Hz. In each meditative or neutral block, we counted the number of electrode pairs (294 electrode pairs maximum) that had an average density of synchrony higher than those derived from noise surrogates (see Methods). We ran a group analysis on the size of the synchronous pattern and found that its size was greater for long-time practitioners than for controls [F(1, 16) = 10.3, P < 0.01; ANOVA] and increased from neutral to meditation states [F(1, 16) = 8.2, P < 0.02; ANOVA]. Fig. 3c shows that the group and state factors interacted on long-distance synchrony [F(1, 16) = 6.5, P < 0.05; ANOVA]: The size of synchrony patterns increased more for the long-time practitioners than for the controls from neutral to meditation states. These data suggest that large-scale brain coordination increases during mental practice.
Finally, we investigated whether there was a correlation between the hours of formal sitting meditation (for subjects S1-S8, 9,855-52,925 h) and these electrophysiological measures for the long-term practitioners, in either the initial or meditative states (same values as in Figs. 2 and 3). The correlation coefficients for the relative, absolute, and phase-synchrony gamma measures were positive: r = 0.79, 0.63, and 0.64, respectively, in the initial state, and r = 0.66, 0.62, and 0.43, respectively, in the meditative state. A significant positive correlation was found only in the initial baseline for the relative gamma (r = 0.79, P < 0.02) (Fig. 3d). These data suggest that the degree of training can influence the spectral distribution of the ongoing baseline EEG. The age of the subject was not a confounding factor in this effect as suggested by the low correlation between the practitioner age and the relative gamma (r = 0.23).
| Discussion |
|---|
|
|
|---|
The high-amplitude gamma activity found in some of these practitioners are, to our knowledge, the highest reported in the literature in a nonpathological context (23). Assuming that the amplitude of the gamma oscillation is related to the size of the oscillating neural population and the degree of precision with which cells oscillate, these data suggest that massive distributed neural assemblies are synchronized with a high temporal precision in the fast frequencies during this state. The gradual increase of gamma activity during meditation is in agreement with the view that neural synchronization, as a network phenomenon, requires time to develop (24), proportional to the size of the synchronized neural assembly (25). But this increase could also reflect an increase in the temporal precision of the thalamo-cortical and corticocortical interactions rather than a change in the size of the assemblies (8). This gradual increase also corroborates the Buddhist subjects' verbal report of the chronometry of their practice. Typically, the transition from the neutral state to this meditative state is not immediate and requires 5-15 s, depending on the subject. The endogenous gamma-band synchrony found here could reflect a change in the quality of moment-to-moment awareness, as claimed by the Buddhist practitioners and as postulated by many models of consciousness (26, 27).
In addition to the meditation-induced effects, we found a difference in the normative EEG spectral profile between the two populations during the resting state before meditation. It is not unexpected that such differences would be detected during a resting baseline, because the goal of meditation practice is to transform the baseline state and to diminish the distinction between formal meditation practice and everyday life. Moreover, Gusnard and Raichle (28) have highlighted the importance of characteristic patterns of brain activity during the resting state and argue that such patterns affect the nature of task-induced changes. The differences in baseline activity reported here suggest that the resting state of the brain may be altered by long-term meditative practice and imply that such alterations may affect task-related changes. Our practitioners and control subjects differed in many respects, including age, culture of origin, and first language, and they likely differed in many more respects, including diet and sleep. We examined whether age was an important factor in producing the baseline differences we observed by comparing the three youngest practitioners with the controls and found that the mean age difference between groups is unlikely the sole factor responsible for this baseline difference. Moreover, hours of practice but not age significantly predicted relative gamma activity during the initial baseline period. Whether other demographic factors are important in producing these effects will necessarily require further research, particularly longitudinal research that follows individuals over time in response to mental training.
Our study is consistent with the idea that attention and affective processes, which gamma-band EEG synchronization may reflect, are flexible skills that can be trained (29). It remains for future studies to show that these EEG signatures are caused by long-term training itself and not by individual differences before the training, although the positive correlation that we found with hours of training and other randomized controlled trials suggest that these are training-related effects (2). The functional consequences of sustained gamma-activity during mental practice are not currently known but need to be studied in the future. The study of experts in mental training may offer a promising research strategy to investigate high-order cognitive and affective processes (30).
| Acknowledgements |
|---|
| Footnotes |
|---|
Freely available online through the PNAS open access option.
Abbreviations: ROI, region of interest; EEG, electroencephalogram.
To whom correspondence may be addressed. E-mail: alutz{at}wisc.edu or rjdavids{at}wisc.edu.
© 2004 by The National Academy of Sciences of the USA
| References |
|---|
|
|
|---|
This article has been cited by other articles in HighWire Press-hosted journals:
![]() |
L. Zylowska, D. L. Ackerman, M. H. Yang, J. L. Futrell, N. L. Horton, T. S. Hale, C. Pataki, and S. L. Smalley Mindfulness Meditation Training in Adults and Adolescents With ADHD: A Feasibility Study J Atten Disord, May 1, 2008; 11(6): 737 - 746. [Abstract] [PDF] |
||||
![]() |
W. J. Rejeski Mindfulness: Reconnecting the Body and Mind in Geriatric Medicine and Gerontology Gerontologist, April 1, 2008; 48(2): 135 - 141. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. D. Kelly Buddhist Psychology, Psychotherapy and the Brain: A Critical Introduction Transcultural Psychiatry, March 1, 2008; 45(1): 5 - 30. [Abstract] [PDF] |
||||
![]() |
R. C. Muresan, O. F. Jurjut, V. V. Moca, W. Singer, and D. Nikolic The Oscillation Score: An Efficient Method for Estimating Oscillation Strength in Neuronal Activity J Neurophysiol, March 1, 2008; 99(3): 1333 - 1353. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. N. Singh, G. E. Lancioni, A. S. W. Winton, J. Singh, W. J. Curtis, R. G. Wahler, and K. M. McAleavey Mindful Parenting Decreases Aggression and Increases Social Behavior in Children With Developmental Disabilities Behav Modif, November 1, 2007; 31(6): 749 - 771. [Abstract] [PDF] |
||||
![]() |
J. A. Brefczynski-Lewis, A. Lutz, H. S. Schaefer, D. B. Levinson, and R. J. Davidson Neural correlates of attentional expertise in long-term meditation practitioners PNAS, July 3, 2007; 104(27): 11483 - 11488. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. D. Koss-Chioino Spiritual Transformation, Relation and Radical Empathy: Core Components of the Ritual Healing Process Transcultural Psychiatry, December 1, 2006; 43(4): 652 - 670. [Abstract] [PDF] |
||||
![]() |
A. Hankey Studies of Advanced Stages of Meditation in the Tibetan Buddhist and Vedic Traditions. I: A Comparison of General Changes Evid. Based Complement. Altern. Med., December 1, 2006; 3(4): 513 - 521. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. D. Wright Meditation: A New Role for an Old Friend American Journal of Hospice and Palliative Medicine, August 1, 2006; 23(4): 323 - 327. [Abstract] [PDF] |
||||
![]() |
M. Galanter Innovations: Alcohol & Drug Abuse: Spirituality in Alcoholics Anonymous: A Valuable Adjunct to Psychiatric Services Psychiatr Serv, March 1, 2006; 57(3): 307 - 309. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. B. Strand Enhanced Communication by Developing a Non-anxious Presence: A Key Attribute for the Successful Veterinarian J Vet Med Educ, January 1, 2006; 33(1): 65 - 70. [Abstract] [Full Text] [PDF] |
||||
![]() |
Practice Makes Perfect? Journal Watch Psychiatry, December 8, 2004; 2004(1208): 6 - 6. [Full Text] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||