Concentration Meditation Block Data.
In the Med. block paradigm, participants performed concentration meditation, focusing on a simple visual stimulus, alternating with a specific form of a neutral, resting state while brain function was recorded with functional MRI. The patterns of significant activation for the Med. blocks vs. the Rest blocks are shown for EMs (see
Fig. 1A) and NMs (
Fig. 1B) on cortical surface models (
21). EMs and NMs activated a large overlapping network of attention-related brain regions, including frontal parietal regions, lateral occipital (LO), insula (Ins), multiple thalamic nuclei, basal ganglia, and cerebellar regions (
Tables 1 and
2). Only NMs showed negative activation (Rest > Med.) in anterior temporal lobe bilaterally (blue hues in
Fig. 1B).
As predicted in our hypothesis, in Med. vs. Rest, EMs showed greater activation than NMs in multiple attentional and other regions including frontoparietal regions, cerebellar, temporal, parahippocampal, and posterior (P.) occipital cortex, likely including the foveal visual cortex of the attended dot (red in
Fig. 1C and
Tables 1 and
2). NMs showed more activation than EMs in medial frontal gyrus (MeFG)/anterior cingulate (Acc) and in the right mid-Ins to P. Ins (
Fig. 1C and
Table 3), regions that have been shown to negatively correlate with performance in a sustained attention task (
16,
22).
We were concerned that these differences may have resulted in part from structural differences between participant-group brains, because seven of 12 EMs were Asian (five Caucasian), and all NMs were Caucasian. Therefore, we performed a separate analysis in which structural differences were taken to account by using probability of gray matter maps as voxel-wise covariates in a
t test comparison between groups (
23). All significant regions remained significant in this analysis, and several regions just below threshold became larger and thus survived multiple correction [
supporting information (SI) Fig. 3A]. In addition, we were concerned with possible motivation differences between groups. Therefore, to better match motivational arousal, we collected data from a set of 10 INMs who were told they would receive a monetary award ($50) if they were in the top one-third of the INMs in activating attention-related regions.
We examined all participant groups, including the INMs, using
a priori regions of interest (ROIs) from a metaanalysis of 31 studies involving attention-shifting paradigms (
24). The EMs showed significantly more activation (two-tailed
t test) than NMs in all attention ROIs except the thalamus (red vs. dark blue in
Fig. 1D). However, the INMs (light blue) showed more activation than the NMs and were not significantly different from the EMs in these ROIs. In the
t test of EMs vs. INMs, EMs had more activation in the left superior frontal gyrus (SFG)/middle frontal gyrus (MFG), and INMs had more activation in left P. Ins, left inferior frontal gyrus (IFG), and LO (
SI Fig. 3 B and C and
Table 2).
Next, because we predicted that these results would correlate with hours of practice, we split the EM group into those with the most hours of practice (top four MHEMs, mean hours = 44,000, range 37,000–52,000, mean age 52.3 years) and those with the least hours of practice (lower four LHEMs, mean hours = 19,000, range 10,000–24,000, mean age 48.8 years, youngest participant not included to ensure age-matching). Two Asians and two Caucasians were in each group. Consistent with an inverted u-shaped function, we found that the LHEMs (brown) had the strongest activation, significantly higher than both sets of NMs in all attention ROIs except left intraparietal sulcus (IPS) and LO (
SI Fig. 3D) and significantly higher than MHEMs (orange) in all ROIs except LO. Results were not significantly different when the top five MHEMs were used (rather than the top one-third; data not shown) nor when the youngest LHEM was used (making mean age 42.3 years), except in thalamus ROI, in which LHEMs were not significantly different (same trend) from INMs or MHEMs (the thalamus ROI was more posterior than the thalamus cluster activated in our study).
In addition, we performed correlations with hours of practice within the EM group. Because age was a potential confound, we calculated the correlation between a participant's age and hours of practice. This was not significant (
r = 0.22,
P < 0.44), but it had a positive trend of older participants having more hours. Thus, we list partial
r values for activation vs. hours of practice, accounting for age. Many regions, including those in the attention network, showed significant negative correlation with hours, whereas no regions showed positive correlation with hours (see last columns of
Tables 1 and
2,
SI Table 4, and
SI Fig. 4A), consistent with the view that expertise may lead to decreased activation, possibly because of increased processing efficiency. The notion of increased processing efficiency in long-term practitioners is consistent with recent evidence from our laboratory using another task, the attentional blink task, where we found that a 3-month period of intensive meditation led to decreased amplitude of the late component of the event-related potential to an initial target, a marker of increased processing efficiency that predicted improved behavioral performance on a subsequent target (
5).
We reasoned that if these results could be explained by differences in the amount of effort required to maintain attentional focus with expertise, one should see differences in the time courses of the hemodynamic response. In the left dorsal lateral prefrontal cortex (DLPFC) ROI, MHEMs had only a short activation period at the beginning of the Med. block (
P < 0.02) that returned to baseline within the first 10–20 sec (significantly less than the LHEMs;
P < 0.001). In contrast, LHEMs had a larger, sustained response over the duration of the block (
Fig. 1E). This short vs. sustained response contributed in part to the decreased activation for MHEMs vs. LHEMs in the attention-related ROIs (
Fig. 1D) because the hemodynamic response function we used in our analysis modeled a continuous response over the entire block. “Meditation startup” increases occurred in most attention ROIs except for the thalamus and left anterior Ins and were also seen in right fusiform gyrus and bilateral caudate. Several other types of responses were seen in MHEMs, including suppression in regions like MeFG/Acc and P. Ins and more sustained responses in IPS, LO, inferior occipital, SFG, and MFG (regions with activation in last 80 sec of Med.) [
SI Fig. 3E (
P < 0.05 uncorrected); also see representative time course plots in
SI Fig. 4B]. The left SFG/MFG region overlapped with the only region that was significantly greater in the 12 EMs vs. INMs (compare
SI Fig. 3 C and E; see also
Table 2).
If the hemodynamic time course is influenced by effort, one should also see a more sustained response in the highly motivated INMs compared with the regular NMs. Indeed, INMs had a greater sustained response than the NMs in which activation at times fell within baseline levels. However, both NM groups had reduced sustained activation over time compared with LHEMs and also showed a delay in the amount of time it took to reach maximum activation in these regions, typically 10–20 sec longer. These results are presented for the DLPFC ROI in
Fig. 1F. All groups had significant (NMs and LHEMs) or near significant (INM and MHEMs;
P < 0.06) activation in the first 10 sec of the meditation block (LHEMs significantly greater than all other groups). However, for the last 80 sec of the block, there was an inverted u-shaped curve in which activation for NM < INM < LHEM > MHEM (all groups significantly different from each other;
P < 0.001). However, whether these activation differences are due to skill learning or strategy and task performance differences cannot be definitely resolved in this study.
Because MHEMs may have been able to reach a less effortful tranquil meditation state within these short blocks, it is possible that regions that remained active in the latter part (last 80 sec) of the meditation block for the MHEMs may be the minimal brain regions necessary to sustain attention on a visual object.
Distracting Sound Data.
In addition to looking at the brain regions involved in generating and sustaining the meditation state, we examined event-related neural responses during presentation of distracting sounds, presented at 2-sec intervals during the last two-thirds of the Med. and Rest blocks. These sounds could be neutral (restaurant ambiance), positive (baby cooing), or negative (woman screaming) and were contrasted with randomly presented silent, null events with the same timing. In this paradigm, 13 EMs, 13 NMs, and 10 INMs were included (see
Methods). General auditory processing pathways (temporal cortex and Ins) were commonly activated for all participant groups in response to distracting sounds during both states (data not shown). A state ANOVA (sounds in Med. vs. Rest) revealed that participants showed an overall “active response” (no suppressed regions) in response to the sounds in Med., involving regions such as right intraparietal lobule/temporal parietal junction, bilateral pre- and post-central sulci, DLPFC, Ins, and anterior SFG (see
SI Table 5 for state effects for all three groups; also see
SI Fig. 5 A–C).
Next we looked for differences between the groups. Our hypothesis predicted that NMs would be more distracted by the sounds and thus would show more activation in default-mode regions related to task-irrelevant thoughts and in emotion regions. First, NMs did not have any regions that were more active than either the EMs or INMs [
SI Fig. 5C vs.
A and
B; see also state-by-group (EM vs. NM) ANOVA in
SI Table 6]. These reduced differences for NMs may have been due to the greater similarity between Med. and Rest states for these participants, as we saw in the Med. block data. Therefore, we viewed the better motivated INMs as the more appropriate control group who would more accurately demonstrate the full potential of novices. As predicted, EMs had less involvement than INMs in medial “default-mode network” regions such as P. cingulate (P. Cing)/precuneus and MeFG/Acc [
Fig. 2A,
SI Fig. 4C (state by group, EM vs. INM), and
Table 3]. EMs also had less activation in left DLPFC, caudate, and pulvinar (
Table 3). In contrast, EMs showed more activation than INMs in bilateral dorsal IPS extending into post-central sulcus, visual cortex, and left, IFG (area 47) (
Table 3).
According to our hypothesis, areas that showed differential effects for EMs vs. NMs should show similar trends when comparing MHEMs vs. LHEMs. A voxel-wise analysis identified multiple regions in which activation in response to sounds correlated with hours of practice (see
Fig. 2 B and
C,
SI Table 7, and
SI Fig. 4 D and E). When all sounds were included together (positive, negative, and neutral), the voxel-wise regression identified negative correlation with hours of practice in multiple regions including right amygdala (Amyg), MeFG/Acc, and P. Cing (
19,
25) (see
Fig. 2 B and
C and
SI Table 7). This P. Cing cluster partially overlapped the P. Cing region more active in INM vs. EMs (compare
A and
B in
Fig. 2). In addition, there was negative correlation with hours of practice in intraparietal lobule, fusiform, and P. temporal regions. There were also several regions with positive correlation with hours of practice, including Ins, subthalamic, left IFG, supplementary motor area, and others; however, slopes of these correlations were usually less steep than areas showing negative correlation (see
Fig. 2B,
SI Table 7, and examples in
SI Fig. 4 D and E). Partial correlations are reported here because the participants included in these analyses showed a substantial but nonsignificant positive association between age and hours of practice (
r = 0.53,
P < 0.08).
Voxel-wise regression of brain responses of each sound valence separately vs. hours of practice identified similar regions (compared with all sounds together) for positive and neutral sounds (data not shown). In response to negative sounds in the EMs, there was a significant inverse correlation between MR signal change in the Amyg and MeFG/Acc and hours; a greater number of hours was associated with less activation to negative sounds in these brain regions (
SI Table 8). These regions overlapped with results from a state by group (EMs vs. INMs) ANOVA for negative sounds, in which INMs also showed more activation than EMs in default network regions (compare
F and
G in
SI Fig. 4) and in right Amyg (compare
D and
E in
Fig. 2). The correlation with hours for negative sounds within the EMs was significantly greater than the correlation for positive (happy baby) sounds in the Amyg (negative sounds, partial
r = −0.64; positive sounds, partial
r = −0.13; difference, Steiger's
Z = 2.6 and
P < 0.04) and in MeFG/Acc (left MeFG, negative sounds, partial
r = −0.86, positive sounds, partial
r = 0.33, Steiger's
Z = 3.3,
P < 0.01; right MeFG, negative sounds, partial
r = −0.81, positive sounds, partial
r = 0.41, Steiger's
Z = 2.4,
P < 0.05). Differences between zero order
r values (without age statistically removed) are also significant (data not shown).
The only positive correlations between response to the negative sounds in Med. and hours of practice were seen in left cerebellar tonsil and subthalamic regions (
SI Fig. 4 G and H and
SI Table 8).