Next Article in Journal
A Clarion Call for Neuroinflammatory Assays in Brain Stimulation for Chronic Pain
Next Article in Special Issue
Mindfulness, Interoception, and the Body
Previous Article in Journal
Biomarkers of Brain Injury: A Window on Mechanisms of Injury and Recovery in the Brain
Previous Article in Special Issue
Mindfulness-Based Interventions and Body Awareness
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Neurofeedback-Augmented Mindfulness Training Elicits Distinct Responses in the Subregions of the Insular Cortex in Healthy Adolescents

1
Laureate Institute for Brain Research, Tulsa, OK 74136, USA
2
Department of Community Medicine, University of Tulsa, Tulsa, OK 74104, USA
3
Department of Psychiatry and Behavioral Sciences, Stanford University, Stanford, CA 94305, USA
*
Author to whom correspondence should be addressed.
Brain Sci. 2022, 12(3), 363; https://doi.org/10.3390/brainsci12030363
Submission received: 17 January 2022 / Revised: 25 February 2022 / Accepted: 3 March 2022 / Published: 9 March 2022
(This article belongs to the Special Issue Mindfulness, Interoception, and the Body)

Abstract

:
Mindfulness training (MT) reduces self-referential processing and promotes interoception, the perception of sensations from inside the body, by increasing one’s awareness of and regulating responses to them. The posterior cingulate cortex (PCC) and the insular cortex (INS) are considered hubs for self-referential processing and interoception, respectively. Although MT has been consistently found to decrease PCC, little is known about how MT relates to INS activity. Understanding links between mindfulness and interoception may be particularly important for informing mental health in adolescence, when neuroplasticity and emergence of psychopathology are heightened. We examined INS activity during real-time functional magnetic resonance imaging neurofeedback-augmented mindfulness training (NAMT) targeting the PCC. Healthy adolescents (N = 37; 16 female) completed the NAMT task, including Focus-on-Breath (MT), Describe (self-referential processing), and Rest conditions, across three neurofeedback runs and two non-neurofeedback runs (Observe, Transfer). Regression coefficients estimated from the generalized linear model were extracted from three INS subregions: anterior (aINS), mid (mINS), and posterior (pINS). Mixed model analyses revealed the main effect of run for Focus-on-Breath vs. Describe contrast in aINS [R2 = 0.39] and pINS [R2 = 0.33], but not mINS [R2 = 0.34]. Post hoc analyses revealed greater aINS activity and reduced pINS activity during neurofeedback runs, and such activities were related to lower self-reported life satisfaction and less pain behavior, respectively. These findings revealed the specific involvement of insula subregions in rtfMRI-nf MT.

1. Introduction

Mindfulness refers to the moment-to-moment, non-judgmental awareness that is cultivated by paying attention in a particular way: on purpose, in the present moment, and non-reactively [1]. Dispositional mindfulness appears to be positively related to psychological health [2,3], and mindfulness training (MT) has been shown to enhance psychological well-being in healthy populations [4,5], as well as benefit various health conditions including chronic pain, stress, anxiety, depression, addiction, and suicidality [6,7,8]. Adolescence represents a sensitive developmental period for emergence of psychopathology, including anxiety disorders [9], depression [10], and suicidal behavior [11]. Although early interventions such as school-based MT have demonstrated effectiveness in enhancing well-being in healthy adolescents [12,13], there is still a need to delineate the precise mechanisms by which MT exerts its mental health benefits to optimize prevention and early intervention efforts in youth.

1.1. The Role of Interoception in Mindfulness and Its Key Hub Insula

By focusing attention on the body and on the “here and now”, MT directly cultivates interoception [14,15,16,17], the individual’s ability to focus attention on the current body state to select and motivate most appropriate regulatory behaviors [18]. In fact, interoception has been proposed as one of the core mechanisms of action of MT [17,19,20,21]. MT ranging from 3 to 12 weeks has been found to improve interoceptive awareness as measured by the self-reported Multidimensional Assessment of Interoceptive Awareness (MAIA) in various clinical populations, including patients with chronic pain and comorbid depression, who reported greater increase in the self-regulation and non-distracting subscales after the 8 week Mindfulness-Based Cognitive Therapy (MBCT) [16]; veterans with post-traumatic stress disorder (PTSD) who reported increases in the self-regulation, body listening, and emotional awareness subscales after a 12 week integrative exercise program [22]; and depressed patients who showed improved regulatory and belief-related aspects of interoception as measured by the attention regulation, self-regulation, body listening, and trusting subscales after a 3 week MBCT [23].
The insular cortex (INS) has been identified as the key hub for interoception [24,25,26,27,28]. Characterizing increased complexity in INS functioning, INS can be subdivided into three subregions from posterior to anterior orientation (i.e., posterior-to-anterior progression hypothesis) [25,29]. Specifically, the posterior INS (pINS) is primarily involved in encoding interoceptive signals, the mid-INS (mINS) is important in integrating interoception with motivated behavior, and the anterior INS (aINS) encodes both introspective emotional awareness and bodily sensations [26,30]. Although several neuroimaging studies have examined the relation between MT and INS activity, the majority of them have reported INS activity as a whole, while extant findings are largely mixed. Earlier research suggested that MT reduced self-referential processing that was accompanied by decreased aINS activations during meditation among experienced meditators [31]. However, others have proposed that INS would show increased activation due to the stronger interoception that is cultivated through MT [19]. Indeed, a review of nine neuroimaging studies found increased INS activity following either the eight-week Mindfulness-Based Stress Reduction (MBSR) course [32] or mindfulness-based tasks among stressed, anxious, and healthy participants [33]. The majority of the literature evaluating insula subregions has reported effects on aINS, rather than mINS or pINS [34]. The mixed INS findings may relate to there being functionally distinct INS subregions involved in affective relative to somatosensory processing or may be explained by unaccounted individual, measurement, or treatment differences. Nevertheless, mixed findings and limited empirical data linking MT and activity in INS subregions require further clarification of the role of INS and interoception in MT.

1.2. Real-Time Functional Magnetic Resonance Imaging Neurofeedback-Augmented MT Targeting the Posterior Cingulate Cortex

The use of real-time functional magnetic resonance imaging neurofeedback (rtfMRI-nf) provides the means of experimentally testing and modulating brain networks central to MT, thereby allowing for a more direct assessment of the role INS plays in this process. MT engages distributed network of brain regions, particularly the regions of the default mode network (DMN) [35,36,37,38], whose central hub is the posterior cingulate cortex (PCC) that supports self-referential processing, the cognitive process of relating information to the self [39]. Indeed, changes in the PCC are among the most robust findings in mindfulness studies, as demonstrated by a recent review where near half of 49 neuroimaging studies reported significant effects of MT on the PCC [40]. PCC has been shown to be activated during mind wandering and self-referential processing [41,42], and deactivated during various mindfulness/meditation practices [43]. Importantly, both INS and PCC play interconnected roles in MT that PCC is also known to be involved in self-awareness [44,45], a construct is often used to measure self-reported interoception (i.e., the Self-Awareness Questionnaire [46] and the How do you feel questionnaire [47]). Given that PCC may be the key modulatory target of MT, neurofeedback targeting the PCC during MT has been successfully implemented in adults [48] and in our recent work with adolescents [49]. Specifically, we report the feasibility and tolerability of PCC-targeted rtfMRI-nf-augmented mindfulness training (MT) in the current sample as established by self-report and successful downregulation of PCC as a function of rtfMRI-nf-augmented mindfulness training and above and beyond mindfulness training alone. In addition, the whole-brain and functional connectivity analyses found co-modulation of a range of regions in the default mode and salience networks, including with the pINS.

1.3. The Present Study

Based on previous literature and our findings pointing to co-modulation of the pINS consequent to rtfMRI-nf targeting the PCC during mindfulness training [49], the present study aimed to clarify the role of INS in MT. We examined activity in distinct INS subregions during MT and its modulation with rtfMRI-nf-augmented MT (NAMT) targeting PCC in healthy adolescents. We selected this developmental period because adolescent brain neuroplasticity allows for effective attempts at improving learning and performance, and therefore makes it a critical period to study neural correlates of psychological and behavioral strategies and their optimization [50,51]. Because aINS is better studied than mINS and pINS in mindfulness research [34], the primary outcome was aINS activity during Focus-on-Breath (i.e., MT) relative to Describe (i.e., self-referential processing). The secondary outcomes included mINS and pINS activity under the same condition. Given that insula activation has been found to be the most consistent effect observed following MT in meditation-naïve participants [34], we hypothesized that (H1) relative to self-referential processing, aINS would show increased activity during MT in healthy adolescents, and (H2): this activation will be further modulated consequent to rtfMRI-nf. Our approach regarding the activity in mINS and pINS in relation to MT was exploratory given the limited research in this area [34]. For exploratory outcomes, we collected self-reported measures concerning task adherence (task ratings) measured immediately following each task run, the State Mindfulness Scale (SMS) measured pre- and post-NAMT [52], as well as affective and sensory awareness assessed by the Patient-Reported Outcomes Measurement Information System (PROMIS) Pediatric scales measured pre-NAMT, including life satisfaction [53], meaning and purpose [54], positive affect [55], pain behavior [56], pain interference [57], and fatigue [58]. Based on the posterior-to-anterior progression hypothesis [25,29] and awareness of emotion and sensory states being core to interoception and mindfulness [14,59], we hypothesized that aINS activity would be related to affective awareness (i.e., life satisfaction, meaning and purpose, and positive affect), while mINS and pINS activity would be related to sensory awareness such as pain behavior and pain interference (H3).

2. Materials and Methods

2.1. Participants

Adolescents in the present study took part in a larger ongoing longitudinal study, with recruitment taking place between September 2019 and July 2021. Adolescents were recruited from the community using flyers, radio and social media advertisements, billboards, and a school-based messaging platform (i.e., PeachJar). A phone screen determined initial eligibility. Remote and in-person visits with adolescents and primary caregivers provided demographic information, medical and psychiatry history, pubertal status, family history of psychiatric illness, and an MRI safety questionnaire. Eligible adolescents were between 13 and 17 years of age at the time of enrollment, had a parent or a legal guardian able to provide consent, were psychiatrically and physically healthy, and were able to validly and safely complete baseline assessments. All races and genders were included. Adolescents were excluded if diagnosed with a neurological or developmental disorder, were currently being managed for migraines (e.g., daily prophylactic medication), had history of traumatic brain injury, had a lifetime history of psychopathology, were currently using medications with major effects on brain function or blood flow (e.g., acne medication), and/or reported MRI contraindications. Forty adolescents were consented for the present study, with two adolescents withdrawn due to repeated missed appointments following consent procedures, and one participant not having usable data due to technical difficulties, for a total of 37 healthy adolescents included in current analyses. According to Desmond and Glover [60], a sample size of 24 is recommended for typical within-group fMRI experiments, in which inferences regarding the differences in activation between two or more conditions are intended to be made in a single population. Therefore, with N = 37, we were 98% powered to detect medium-size effects (f = 0.25) between conditions. Thirty-four subjects overlapped with the sample in [49].

2.2. Experimental Procedures

This study involves analysis of an ongoing study focused on the impact of self-regulation of PCC using rtfMRI-nf. The current study presents a tangential analysis to the primary outcomes as delineated on Clinical Trials (www.clinicaltrials.gov NCT04053582).

2.2.1. Neurofeedback-Augmented Mindfulness Training Task (NAMT)

The NAMT task (Figure 1) has been previously described [49], and further task details can be found in Appendix A. Briefly, adolescents were first given a brief psychoeducational introduction into mindfulness, followed by a guided traditional mindfulness practice focused on the breath [43,48]. Next, adolescents completed the same mindfulness practice and assessment with MRI noises in the background in the mock scanner. MT was manualized to ensure fidelity across participants. MT was delivered by a trained research assistant under the supervision of a licensed clinical psychologist. Training sessions were audio recorded and up to 20% sessions were randomly selected for fidelity ratings by research staff using an unpublished measure developed by NK in consultation with RLA for the purposes of this study. On 3-point Likert scale (0 = no adherence, 1 = adherence identified but weak or flawed, 2 = good adherence) of how closely the research assistant followed the manualized mindfulness training, the fidelity ratings indicate that the manualized training was delivered with satisfactory adherence (M = 95.81%, SD = 3.33%) (the percentage was obtained by summed score/highest total score possible).
The neuroimaging session included 8 runs (Figure 1a), including an anatomical scan, Resting State scan 1 (Rest-1), Observe (OBS), three neurofeedback runs (NF-1, NF-2, NF-3), Transfer run (TRS), and Resting State scan 2 (Rest-2). During Rest-1 and Rest-2 (6 min each), participants were instructed to clear their mind and not think about anything while fixating upon a fixation cross. OBS, NF-1, NF-2, NF-3, and TRS runs each lasted 6 min and 56 s. Runs started with a 66 s rest block, followed by alternating Describe (Active Control condition without neurofeedback; 20 s), Focus-on-Breath (MT condition with PCC neurofeedback; 70 s), and Rest (Baseline condition; 30 s) blocks. OBS and TRS runs did not involve neurofeedback (no bar displayed) during the Focus-on-Breath condition. During the Focus-on-Breath condition [48], adolescents were instructed to pay attention to the physical sensations of their breath, not trying to change their breathing in any way, and if their attention were to wander to something else, to gently bring their attention back to their breath [43]. In the Describe condition, adolescents were presented with various adjectives, which they had to mentally categorize as descriptive or not descriptive of them for the entire duration the word was displayed on the screen [61]. During neurofeedback runs, adolescents were told that they would see a bar displayed on the screen, representing the relative brain activity in a particular brain region in real time (Figure 1b). The instructions further indicated that the bar may change with the experience of focusing on the breath (i.e., the bar may go blue if they are fully concentrating on their breath, and red if their mind wanders elsewhere). The consensus on the reporting and experimental design of clinical and cognitive-behavioral neurofeedback studies (CRED-nf checklist) [62] is included in Appendix B Table A1.

2.2.2. Psychological Measurements

Following each fMRI run, participants completed task ratings, including the following items: “How much did your mind wander while you were asked to focus on your breath?”, “How easy did you find it to focus on your breath?”, “How do you feel right now?”, each rated from 1 to 10. Participants completed the SMS at T1 (pre-) and T2 (post-NAMT) to quantify adolescents’ perceived level of attention to and awareness of their present experience (i.e., mind, body, the pleasant/unpleasant/neutral hedonic tones of these objects of awareness, and the qualities thought to characterize mindful awareness) [52]. Additionally, adolescents completed the following PROMIS Pediatric scales pre-NAMT: life satisfaction (example items: Thinking about the past 4 weeks, “I was satisfied with my life”, “I had a good life”, 1—Not at all, 5—Very much) [53], meaning and purpose (example items: Thinking about my life “I can reach my goals in life”, “I have a reason for living”, 1—Not at all, 5—Very much) [54], positive affect (example items: In the past 7 days, “I felt happy”, “I felt peaceful”,1—Never, 5—Always) [55], pain behavior (example items: In the last 7 days, when I was in pain, “It showed on my face”, “I talked about my pain.” 1—Never, 5—Always) [56], pain interference (example items: In the past 7 days, “I had trouble sleeping when I had pain”, “I felt angry when I had pain”, 1—Not at all, 5—Very much) [57], and fatigue (example items: In the past 7 days, “How often did you have to push yourself to get things done because of your fatigue?”, “I have trouble starting things because I am tired”, 1—Not at all, 5—Very much) [58].

2.2.3. Data Acquisition

Neuroimaging was performed using a GE MR750 3T MRI scanner with the 8-channel receive-only head coil. To acquire T1-weighted anatomical images, a 3D magnetization-prepared rapid gradient echo (MPRAGE) pulse sequence accelerated with sensitivity encoding (SENSE) [63] was used. The MPRAGE parameters were as follows: FOV/slice thickness = 240/1.2 mm, axial slices per slab = 128, image matrix size = 256 × 256, TR/TE = 5.0/1.9 ms, SENSE acceleration factor R = 2, flip angle = 8°, delay/inversion times TD/TI = 1400/725 ms, sampling band- width = 31.2 kHz, and scan time = 5 min 33 s.
For the whole-brain fMRI recording, an accelerated single-shot gradient EPI with SENSE was used. EPI sequence parameters were optimized to maximize sensitivity to BOLD contrast and minimalize image distortion and susceptibility dropouts [64,65]. EPI parameters were as follows: FOV/slice = 240/2.9 mm, TR/TE = 2000/25 ms, SENSE acceleration R = 2, acquisition matrix: 96 × 96, flip angle = 90°, image matrix: 128 × 128, 46 axial slices, and voxel volume: 1.9 × 1.9 × 2.9 mm3. To allow the fMRI signal to reach a steady state, three EPI volumes (6 s) were added at the beginning of each run and were excluded from data analysis. Physiological pulse oximetry and respiration waveforms were recorded simultaneously with fMRI (with 25 ms sampling interval, i.e., the sampling rate of the pulse oximeter and respiration measurements) using a photoplethysmograph placed on the subject’s finger and a pneumatic respiration belt, respectively. rtfMRI-nf procedures are described elsewhere [49] and in Appendix A.

2.2.4. Data Processing and Analysis

AFNI [66] was used for data image analysis. The first 5 fMRI volumes were discarded to wait for a steady state. fMRI data preprocessing included despiking, RETROICOR [67], respiration volume per time correction [68], slice-timing and motion corrections, non-linear warping to the Montreal Neurological Institute (MNI) template brain with resampling to 2 mm3 voxels using the ANTs [69], spatial smoothing with a 6 mm FWHM Gaussian kernel, and scaling signal to percent change relative to the mean in each voxel. The general linear model (GLM) analysis was used for independently evaluating the brain response in the OBS, NF-1, NF-2, NF-3, and TRS runs. One participant included in analysis had no NF-2 due to technical difficulties. The design matrix included a modeled response to the Focus-on-Breath block (boxcar function convolved with hemodynamic response function), 12 motion parameters (3 shift and 3 rotation parameters with their temporal derivatives), three principal components of the ventricle signal, local white matter average signal (ANATICOR) [70], and low-frequency fluctuation (fourth-order Legendre polynomial model).
Our study is one of the first to include real-time physiological noise correction, and then offline RETROICOR correction in our preprocessing pipeline was followed by further physiological noise correction during the subsequent GLM analysis. Such a sequential approach of artifact removal can potentially reintroduce previously removed noise artifacts in later steps [71], thus it is important to evaluate the efficacy of physiological noise correction. We calculated the signal variance ratio (R2 value) explained by the physiological noise regressors (RETROICOR) for the real-time processed and offline processed signals in the PCC region. This measure indicates an amount of residual physiological noise effect in the processed signal [72]. The R2 was small for both real-time and offline-processed signals on average (Appendix B Figure A1). While some subjects’ runs showed a relatively high residual noise variance ratio, there was no significant correlation between the mean neurofeedback signal amplitude and the residual noise variance ratio for the real-time processed data (Spearman’s rho = −0.091, p = 0.361). Although a relatively high residual physiological noise variance ratio was seen for the offline-processed data, the correlation between the PCC parameter estimates (beta value) and the physiological noise variance ratio was not significant for the offline-processed data (Spearman’s rho = 0.031, p = 0.744). These indicate that the physiological noise effect was not significant on the neurofeedback signal in the real-time analysis and the PCC parameter estimation in the offline analysis.
Regression coefficients estimated from the GLM were extracted from both hemispheres of all three probabilistic cytoarchitectonic segmentations of INS regions defined by the Brainnetome atlas [73]: aINS (ventral and dorsal agranular), mINS (ventral dysgranular and granular, and dorsal dysgranular), and pINS (hypergranular and dorsal granular) averaged across left and right hemispheres (Figure 2a). See Appendix B Table A2 for peak coordinates of each subregion. The average parameter estimate (beta coefficient) of the Focus-on-Breath vs. Describe block regressor was extracted to evaluate brain activation during each run (OBS, NF-1, NF-2, NF-3, and TRS; H1).
All remaining statistical analyses were performed using the R statistical package [74]. Descriptive statistics regarding participant characteristics and PROMIS pediatric measures were obtained using the R package ‘psych’ [75]. To test H2, first, Pearson’s correlation was conducted to examine the relation between parameter estimate (Focus-on Breath vs. Describe) in PCC and INS subregions; then separate linear mixed-effects models (LMEs) were conducted to examine task ratings and INS subregions activity across runs (OBS, NF-1, NF-2, NF-3, and TRS) using the ‘lmer’ function in the R package ‘lme4′ [76], with Run entered as a fixed effect and subjects as a random effect. Regarding exploratory outcomes (H3), LME was conducted for SMS, where Time (T1: pre-, T2: post-NAMT) was entered as a fixed effect and subjects as a random effect. Follow-up pairwise comparisons for LMEs were conducted using the ‘glht’ function in R package [77] and corrected for multiple comparisons with Tukey’s Honestly Significant Difference test. Spearman’s correlation analysis was conducted to examine the relation between parameter estimate (Focus-on Breath vs. Describe) in each INS subregion during neurofeedback runs (averaged across NF-1, NF-2, NF-3) and PROMIS pediatric measures. Bonferroni correction was used to correct for multiple correlation comparisons.

2.2.5. Data and Code Availability Statement

The data and data analysis scripts that support the findings of this study are available on request from the corresponding author after a formal data sharing agreement has been signed. The data are not publicly available due to privacy or ethical restrictions.

3. Results

3.1. Demographic, Task, and Clinical Characteristics

Participants were 37 adolescents (age 14.61 ± 1.25 years, 16 female), and the majority were White (71%). Table 1 provides additional demographic information. LME results of self-report data are in Table 2. In general, adolescents reported moderate to high ability to focus on their breath during the Focus-on-Breath condition, moderate mind wandering during the Focus-on-Breath condition, and moderate calmness during the task. The scores on these measures did not differ across fMRI runs (OBS, NF-1, NF- 2, NF-3, and TRS; Focus-on-Breath: [F(4, 135) = 0.93, p = 0.45; R2 = 0.31]; mind wandering: [F(4, 135) = 2.03, p = 0.09; R2 = 0.45]; and current feeling: [F(4, 135) = 0.65, p = 0.63; R2 = 0.56]). State mindfulness (SMS) increased from pre- to post-NAMT session [F(1, 36) = 5.82, p = 0.02; R2 = 0.79]. The PROMIS scales are anchored with a mean of 50 for the United States general population, thus the current sample displayed average levels of life satisfaction (M ± SD = 49.05 ± 8.8), slightly lower levels of meaning and purpose (M ± SD = 47.41 ± 7.59), and slightly higher levels of positive affect (M ± SD = 55.42 ± 6.15), as well as lower levels of fatigue (M ± SD = 43.8 ± 10.62), pain behavior (M ± SD = 35.75 ± 10.52), and pain interference (M ± SD = 40.65 ± 7.48).

3.2. Insula Region of Interest (ROI) Results

Table 3 shows the uncorrected correlations between parameter estimate (Focus-on-Breath vs. Describe) in PCC and insula subregions across runs. Overall, INS subregions were positively correlated with PCC. pINS showed moderate to strong correlations with PCC across runs, while aINS and mINS showed moderate correlation with PCC in all runs except for NF-1, where no correlations were found. Notably, only the correlations between INS subregions and PCC in OBS, as well as pINS and PCC in NF-2 and NF-3 survived the multiple comparison threshold (p = 0.05/15 = 0.003).
LME analyses revealed the main effect of run for parameter estimate (Focus-on-Breath vs. Describe) in aINS [F(4, 143) = 4.98, p < 0.001; R2 = 0.39] and pINS [F(4, 143) = 11.80, p < 0.001; R2 = 0.33], but not mINS [F(4, 143) = 1.56, p = 0.18; R2 = 0.34] (Figure 2; Table 4). Post hoc analyses performed for aINS and pINS subregions revealed significant differences between neurofeedback and non-neurofeedback runs. Specifically, parameter estimate (Focus-on-Breath vs. Describe) was greater in neurofeedback runs than OBS in aINS (p < 0.01 for NF-1 and NF-2, p < 0.05 for NF-3), but was lower in neurofeedback runs than OBS (all ps < 0.001) and TRS (p < 0.01 for NF-1 and NF-3, p < 0.05 for NF-2) in pINS (Table 5). Additional analyses in parameter estimate (Focus-on-Breath vs. Rest) confirmed activation in aINS and deactivation in pINS during neurofeedback runs. For parameter estimate (Focus-on-Breath vs. Rest) across each run for insula subregions, the LME results and post hoc tests are summarized in Appendix B Table A3 and Table A4, respectively. Line graphs are displayed in Appendix B Figure A2.
We also tested parameter estimate (Focus-on-Breath vs. Describe) in aINS, mINS, and pINS across experimental runs, separately (Figure 2). In aINS, parameter estimate for Focus-on-Breath was significantly lower than that in Describe during OBS and TRS [OBS: t(36) = −5.21, p < 0.001, Cohen’s d = −0.86; TRS: t(36) = −2.98, p < 0.01, Cohen’s d = −0.49], but not NF-1, NF-2, and NF-3 runs (all ps > 0.10), in mINS during NF-1 and NF-3 [NF-1: t(36) = −2.13, p < 0.05, Cohen’s d = −0.35; NF-3: t(36) = −2.50, p < 0.05, Cohen’s d = −0.41], but not OBS, NF-2, and TRS (all ps > 0.05), and in pINS during NF-1, NF-2, and NF-3 [NF-1: t(36) = −7.23, p < 0.001, Cohen’s d = −1.19; NF- 2: t(35) = −6.14, p < 0.001, Cohen’s d = −1.01; NF-3: t(36) = −5.83, p < 0.001, Cohen’s d = −0.96], but not OBS and TRS runs (all ps > 0.05).

3.3. Self-Reported Questionnaire Responses and Insula ROI Results

Relationship between insula subregions and self-reported task ratings and psychological function are reported in Appendix B Table A5 and Figure A3. None of the correlations met Bonferroni correction; p = 0.05/9 = 0.005. pINS activity was negatively correlated with mind wandering during OBS, r = −0.34, p < 0.05, while current feeling (1 = perfectly calm, 10 = very anxious) was negatively correlated with aINS, r = −0.43, p < 0.01. and mINS, r = −0.37, p < 0.05 in TRS. During neurofeedback runs, aINS activity was negatively correlated with PROMIS life satisfaction (r = −0.37, p < 0.05), whereas pINS activity was positively correlated with PROMIS pain behavior (r = 0.33, p < 0.05).

4. Discussion

Based on our previous findings that established the successful downregulation of PCC as a function of rtfMRI-nf-augmented mindfulness training, as well as the co-modulation of the posterior INS [49], the present study sought to directly examine the role of INS in mindfulness in healthy adolescents during the PCC-targeted rtfMRI-nf. Given that different subregions of INS specialize in integrating across the body, cognitive, affective, and awareness domains, we focused our examination on activity in the anterior (aINS), mid- (mINS), and posterior (pINS) subregions and their associations with self-reported affective and sensory measures. The observed correlations between PCC activity and insula subregions confirmed previously reported whole-brain and functional connectivity findings. Contrary to our hypothesis, relative to self-referential processing (i.e., Describe), MT (i.e., Focus-on-Breath) significantly reduced aINS activity during non-neurofeedback runs (OBS and TRS). However, this deactivation in aINS was not further modulated consequent to rtfMRI-nf. Instead, as hypothesized, aINS showed more activation during MT when PCC neurofeedback was given (NF-1 to NF-3). Second, we observed an opposite effect in pINS, such that no differences between conditions were present during non-neurofeedback runs, while significant deactivations occurred for MT during neurofeedback runs. Finally, mINS activity did not differ between conditions during non-neurofeedback runs but reduced for MT during the first and the last neurofeedback runs. Correlational analyses showed that pINS activity was negatively correlated with mind wandering during OBS, while aINS and mINS were negatively correlated with current feeling (calm vs. anxious) during TRS. During neurofeedback runs, aINS activity was negatively correlated with life satisfaction and pINS activity was positively correlated with pain behavior, while no associations were found for mINS activity. These findings provide direct evidence for the involvement of INS in MT and demonstrate that MT elicited distinct responses across INS subregions that may partially explain the previously reported mixed findings regarding INS activity during MT.
INS subregions have distinct functions. An activation-likelihood-estimation meta-analysis of 1768 neuroimaging studies revealed that sensorimotor tasks consistently activate mid-posterior INS, whereas social-emotional and cognitive tasks activate aINS [78]. Given that aINS is an important hub for interoceptive awareness, accumulated studies have shown aINS to be sensitive to MT relative to mINS and pINS [79,80,81]. For example, following 7- to 8-week MT, aINS activation was observed in healthy participants while focusing attention on breathing, body, and thoughts [36], in patients with general anxiety disorder performing affect labeling of emotional facial expressions [82], and in elite athletes experiencing an interoceptive challenge [83]. Consistent with prior research, we found that, relative to self-referential processing, mindful attention to breathing led to an increase in aINS activity, however, only during neurofeedback runs.
Recent fMRI studies have suggested aINS to be part of the “salience network”, important in attentional control and detecting behaviorally salient stimuli [84,85,86]. Various factors may contribute to activation in the salience network, including increases in cognitive effort, reward anticipation [87,88], heightened sympathetic arousal [89], and the pressure to perform self-regulation within limited time period [90]. Moreover, aINS is critical in detecting discrepancies between actual and desired states [91,92], as well as integrating the feedback display (external) with brain activity (internal) [93]. Neurofeedback training may enhance these processes because participants must integrate information from the feedback and regulate their breathing in order to obtain the desired state. For example, aINS engagement is observed when participants actively tried to move the neurofeedback signal relative to a passively watching condition [94]. Furthermore, neurofeedback training may enhance self-referential processing due to self-evaluation of one’s performance while attempting to control the neurofeedback bar. Therefore, aINS activity during neurofeedback relative to non-neurofeedback MT runs in the present study may reflect the complex process of integrating state awareness of physical sensations of breath with increased attentional control and self-evaluation. A modest trend toward greater deactivation across neurofeedback runs potentially suggests habituation to the neurofeedback signal bar. Together with the minimal difference in aINS activity between OBS and TRS runs, it is possible that aINS activity during neurofeedback training in part reflects the introduced cognitive load of the feedback signal.
aINS, particularly its dorsal region, has long been linked to subjective well-being, a multidimensional construct that involves both cognitive and affective evaluation of life satisfaction [95,96,97,98,99]. Associations have been found between well-being and INS gray matter volume [100], as well as functional connectivity of the dorsal aINS [98]. Remarkably, aINS activation has been reported to be positively correlated with momentary happiness ratings [95]. This is consistent with our findings which showed a relationship between higher activity in aINS and mINS feeling more calm during the task, as well as between higher activity in pINS and less mind wandering.
It is well known that MT increases subjective well-being by promoting greater awareness of the present moment [101]. Therefore, greater interference of external stimuli (i.e., attentional and evaluative engagement with the neurofeedback signal bar) during MT may indicate overall propensity toward less present moment awareness in the presence of salient and self-relevant stimuli. This may not only further explain the difference in aINS activity between neurofeedback and non-neurofeedback runs, but also the association between mindful awareness and psychological well-being.
pINS is considered the primary interoceptive cortex and plays a crucial role in integrating somatic processing and pain perception [85,102,103,104]. Further, pINS has been shown to play a role in momentary self-referential processing [105]. Previous literature is scant on pINS activity in response to MT. In contrast to one study that reported greater pINS activation during an interoceptive breath-focused task in participants who completed MBSR [102], we found that augmenting MT with neurofeedback decreased pINS activation. Deactivation in pINS during neurofeedback runs might be a result of decreased momentary self-referential processing given that participants were engaged in monitoring the feedback signal and with the physical sensations of their breath. This assumption was further supported by the fact that pINS activity was increased during the TRS where the neurofeedback signal (and thus, attentional engagement with the signal bar) was removed, albeit to a lesser extent than during the OBS run. Consistent with past studies showing increased INS activity during processing of painful stimuli [106,107,108], we found that adolescents who reported more pain behaviors also showed increased pINS activity during MT relative to self-referential processing in neurofeedback runs. pINS activity during NAMT may be related to the individual’s levels of sensory processing of pain (e.g., individuals who have greater sensory awareness of pain may have poor inhibited sensory processing of pain), which may lead to the relatively increased pINS activity during MT compared to self-referential processing; however, the mechanism underlying this is unclear. These findings are further linking INS activity to attentional deployment, engagement with physical sensations, and self-referential processing.
Although mINS activity reduced during NF-1 and NF-3 of MT relative to self-referential processing, it did not differ significantly across experimental runs, nor was it related to any PROMIS pediatric measures. The lack of consistent change in mINS activity during neurofeedback runs could potentially be explained by its unique anatomical position, a connecting structure between aINS and pINS. Thus, mINS activity might be affected by both aINS activation and pINS deactivation during neurofeedback runs [85]. Furthermore, the absence of change in mINS activity across all experimental runs during MT could be due to its main role being integrating interoception with motivated behavior, or selecting action-outcome behaviors [109]. As a result, the current study might fail to engage mINS because it did not provide motivated signals or require action selection. Further research that aims to determine the role of mINS in MT might consider these factors in the experimental design.
INS is involved in a wide variety of functions ranging from sensory and affective processing to high-level cognition [27,85], which may be implicated in the rich cortical connections of INS. For example, aINS has connections with the cingulate cortex, frontal, orbitofrontal, anterior temporal and limbic areas [85,110,111,112,113], affording its role underlying various cognitive and affective functions; as well as the pregenual anterior cingulate and the anterior mid-cingulate that are related with emotional and pain processing [114,115]. mINS mainly projects to the mid-cingulate cortex [116,117], a region important for integrating bodily sensations and exteroceptive sensory afferents [118], supporting its main role as discussed above; while pINS is primarily connected to regions for sensorimotor processing such as posterior temporal, parietal, and sensorimotor areas. Recent rodent model suggests that pINS is also involved in top-down modulation of behavior upon the detection of internal aversive state [119].
Given that attention control and acceptance are the two core components of mindfulness [120,121,122], INS function is particularly relevant for MT due to its involvement in interoception [30], self-awareness [123], attention control, and emotion regulation. Attention control enables attending to the moment-to-moment experience (i.e., breathing), which, in turn, cultivating interoceptive awareness and provides an integrated representation of the present moment [93]. In addition, achieving a non-judgmental attitude toward the present experience requires both cognitive reappraisal and emotion regulation [124]. Thus, INS is integral to mindfulness in that it senses, interprets, integrates, and regulates internal and external inputs [125]. Finally, INS-dependent interoceptive regulation is particularly important during adolescence when self-regulatory abilities rapidly develop [126]. Atypical interoception may contribute to onset of psychopathology and decreased socio-emotional competence in late adulthood [127]. Heightened interoceptive reactivity to pleasant stimuli accompanied by increased pINS activity in adolescents relative to young and mature adults may explain developmentally appropriate increased levels of risky decision making during this period of life [128,129,130]. Taken together, findings from our studies and others suggest that MT may influence INS function and interoceptive processing in clinically relevant ways.

Limitations

Several limitations should be mentioned for this study. First, although this study was well powered to detect differences in activation across experimental conditions in a single population, the sample size of the present study remains relatively small (n = 37). Second, the pain behavior rating is of low variability in nearly half of this healthy sample (n = 15 rated 23.9). Future studies would benefit from more diverse samples across the gender, age, race/ethnicity, and clinical domains. Importantly, clinical samples will give evidence to whether INS-level changes translate into modifications in symptom measures. Third, the present study employed rtfMRI-nf targeting PCC and not INS directly. Future work employing rtfMRI-nf targeting INS might provide a more accurate investigation of MT on INS activity, thus to directly evaluate the success of modulating INS activity with neurofeedback and the downstream effects of NAMT targeting INS. Fourth, we did not employ a sham condition against which the effects of PCC-targeted rtfMRI-nf on INS activity could be evaluated. Furthermore, although we provided the associations between PCC and INS subregions, future work should examine changes in functional connectivity between the INS subregions and PCC as a function of rtfMRI-nf using resting state data to provide more details regarding this relationship. Lastly, it is worth noting that the self-referential task might elicit emotional responses that further induce bodily change because participants were instructed to determine whether the adjectives described themselves or not [131]. Future studies could improve the control task to provide a more stable reference by avoiding instructions that elicit potential bodily changes.

5. Conclusions

This is the first study to examine activity in INS subregions during PCC-targeted rtfMRI-nf MT in healthy adolescents. These findings add to the existing literature for the integral role of INS in MT. The data also showed a relation between INS activity during NAMT and self-reported cognitive/affective/sensory processing. The divergent effect of PCC rtfMRI-nf on anterior vs. posterior INS subregions during mindfulness practice relative to self-referential processing may support previous findings whereby aINS is involved in the experience of cognitive-affective states, while pINS plays a more prominent role in somatosensory processes. Future studies are needed to directly examine how distinct mindfulness practices modulate INS along the proposed subregion specializations given that we focused only on focused-attention MT. Finally, studies with larger and clinical samples will determine whether MT impacts INS activity and interoceptive processes to improve clinical outcomes.

Author Contributions

Conceptualization, N.K. and X.Y.; methodology, N.K., J.B., M.M. and A.T.; software, M.M.; validation, N.K., M.M. and X.Y.; formal analysis, N.K., X.Y. and M.M.; investigation: G.C., Z.P.C. and N.K.; resources, J.B. and M.P.P.; data curation, M.M.; writing—original draft preparation, X.Y. and N.K.; writing—review and editing, X.Y., Z.P.C., A.T., G.C., R.L.A., J.L.S., M.K.S., M.M., M.P.P. and N.K. (J.B. deceased during the preparation of the manuscript); visualization, X.Y. and N.K.; supervision, N.K.; project administration, N.K.; funding acquisition, N.K. and M.P.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the P20 GM121312 award from the National Institute of General Medical Sciences and the Laureate Institute for Brain Research and the William K. Warren Foundation. Authors further report funding support from the U.S. Department of Defense (W81XWH-12-1-0697 award [JB]), National Institute on Drug Abuse (01 DA050677 [JLS]) and National Institutes of Mental Health (R01 MH123691 and K23 MH108707 awards [RA]). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki, and approved by Western Institutional Review Board (protocol code: LIBR # 2019-003-04 and date of approval: 14 July 2020). This study is registered at the US National Institutes of Health (ClinicalTrials.gov identifier #NCT04053582; 12 August 2019).

Informed Consent Statement

Informed consent was obtained from all subjects involved in this study. Parents provided written informed consent, while adolescents provided written assent for study participation.

Data Availability Statement

Data will be available upon request from the corresponding author.

Conflicts of Interest

M.K.S. has received research support from Stanford’s Maternal Child Health Research Institute and Stanford’s Department of Psychiatry and Behavioral Sciences, National Institute of Mental Health, National Institute of Aging, Patient Centered Outcomes Research Institute, Johnson and Johnson, and the Brain and Behavior Research Foundation. She is on the advisory board for Sunovion and Skyland Trail, is a consultant for Johnson and Johnson, has previously consulted for X, moonshot factory, Alphabet Inc. and Limbix Health, and has received honoraria from the American Academy of Child and Adolescent Psychiatry, and royalties from American Psychiatric Association Publishing and Thrive Global. All other authors declare no conflict of interest. The funders had no role in the design of this study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Appendix A.1. Neurofeedback-Augmented Mindfulness Training Task (NAMT)

Prior to the rtfMRI-nf session, adolescents underwent brief mindfulness training (MT). Participants were first given a brief psychoeducational introduction into mindfulness, including that (1) mindfulness refers to paying attention to thoughts, feelings, and physical sensations in the present moment without any judgment, and (2) mindfulness can reduce stress and increase attention. Next, participants were guided through a traditional mindfulness practice focused on the breath [43,48], that is: “Please pay attention to the physical sensations of your breath where you most strongly feel it. Follow the natural and spontaneous movement of the breath, not trying to change it any way. Just pay attention to it. If you find that your attention wanders to something else, gently but firmly bring it back to the physical sensations of the breath.” Difficulty of performing the task and how mindful they currently feel of their body and mind were assessed. Following practice, adolescents were provided with an opportunity to ask clarification questions. Next, adolescents went into the mock scanner and completed the same mindfulness practice and assessment with MRI noises in the background. Adolescents were also given instructions and feedback around minimizing motion while in the scanner. Finally, adolescents were given instructions for the neuroimaging session. Training was manualized to ensure fidelity across participants. MT was delivered by a trained research assistant under the supervision of a licensed clinical psychologist. Training sessions were audio recorded and up to 20% sessions were randomly selected for fidelity ratings by research staff.
The neuroimaging session included 8 runs (Figure 1a), including an anatomical scan, Resting State scan 1 (Rest-1), Observe (OBS), three neurofeedback runs (NF-1, NF-2, NF-3), Transfer run (TRS), and Resting State scan 2 (Rest-2). During Rest-1 and Rest-2 (6 min each), participants were instructed to clear their mind and not think about anything while fixating upon a fixation cross. OBS, NF-1, NF-2, NF-3, and TRS runs each lasted 6 min and 56 s. Runs started with a 66 s rest block, followed by alternating Describe (Active Control condition without neurofeedback; 20 s), Focus-on-Breath (MT condition with PCC neurofeedback; 70 s), and Rest (Baseline condition; 30 s) blocks. OBS and TRS runs did not involve neurofeedback (no bar displayed) during the Focus-on-Breath condition. The initial long Rest block was required for obtaining enough samples for real-time noise regression analysis [132].
During the Focus-on-Breath condition [48], adolescents were instructed to pay attention to the physical sensations of their breath, not trying to change their breathing in any way, and if their attention were to wander to something else, to gently bring their attention back to their breath [43]. To aid in MT, numerous useful strategies were provided prior to scanning, including “Notice the feeling of your belly rising when you breath in, and gently falling when you breath out”; “Notice if it enters and leaves through your nose or your mouth.” In the Describe condition, adolescents were presented with various adjectives, which they had to mentally categorize as descriptive or not descriptive of them for the entire duration the word was displayed on the screen [61]. The Describe condition is designed to elicit self-referential thinking, and therefore is thought to be a better comparator to Focus-on-Breath than Rest [43]. During the Rest condition, adolescents were presented with the cue “Rest” and asked to relax while looking at the display screen.
During neurofeedback runs, adolescents were told that they would see a bar displayed on the screen, representing the relative brain activity in a particular brain region in real time (Figure 1b). The instructions further indicated that the bar may change with the experience of focusing on the breath (i.e., the bar may go blue if they are fully concentrating on their breath, and red if their mind wanders elsewhere). The green bar represented the target to attain, and adolescents’ goal was to try and see how much they could make the bar change to blue to match the green bar. The target levels were −0.5%, −0.75%, and −1.0% (percent signal change is relative to the preceding rest block) for the NF-1, NF-2, and NF-3 runs, respectively. Adolescents were told that there might be a 5–6 s delay between their experience and the change in the blue bar. To assess for aspects of feasibility and tolerability, adolescents answered the following questions via a response box and visual analog scale after reach run: (1) How well were you able to follow instructions on the screen? (2) How easy did you find it to focus on your breath? (3) How much did your mind wander while you were asked to focus on your breath? (4) How easy did you find it to mentally decide whether or not the words described you? (5) How easy did you find it to clear your mind while resting? (6) How do you feel right now (1 = perfectly calm, 10 = very anxious)? Two additional questions followed the neurofeedback runs only: (1) How well did the blue bar correspond with your experience of focusing on your breath? (2) How well did the red bar correspond with the experience of your mind wandering elsewhere?

Appendix A.2. Real-Time fMRI Processing

The region of interest (ROI) and rtfMRI-nf target location (spherical ROI, 7mm radius, [MNI coordinates: x = −5, y = −55, z = 23]; Figure 1a) were selected based on a meta-analysis investigating functional neuroimaging studies of the DMN [133], mindfulness meditation studies, including neurofeedback [43,48], and conducted pilot testing.
We employed an advanced fMRI real-time processing (RTP) protocol which included slice-timing correction, motion correction, spatial smoothing with 6 mm-FWHM Gaussian kernel within the brain mask, scaling to a percent change relative to the average for the first 19 TRs (in the initial rest period), and regressing out noise components [72,132]. The noise regressors were six motion parameters, eight RETROICOR [67] regressors (four cardiac and four respiration), white matter mean signal, ventricle mean signal, and Legendre polynomial models of slow signal fluctuation. This comprehensive noise reduction was performed in real-time (less than 400 ms [132]). This fMRI RTP system operates real-time motion tracking, alignment, and motion parameter regression, thus allowing for suppression of head motion effects, and importantly, providing physiological noise correction (RETROICOR) in real time before the PCC-based neurofeedback signal computation and visual presentation to the adolescent [132,134,135]. This ensured that the PCC neurofeedback signal reflected the largest possible extent the underlying neuronal activity and does not reflect head motion, heart rate, and/or respiratory motions, all of which overlap with DMN neural activity [136,137]. After real-time noise regression, the fMRI RTP system exports the mean value of the noise-reduced signals for the PCC ROI for each acquired data volume.
The PCC ROI (i.e., neurofeedback target signal) in the MNI space was warped into the individual brain space using the Advanced Normalization Tools (ANTs) software [69] (http://stnava.github.io/ANTs/). The neurofeedback stimulus was delivered via custom-developed software using PsychoPy [138]. The neurofeedback value was a signal change relative to the baseline obtained by averaging the preceding 30 s long Rest condition. The two initial volumes in the Rest condition were excluded from the baseline calculation to avoid the delayed hemodynamic response effect of the preceding Describe block. The neurofeedback started from the third volume in the Focus block to wait for a hemodynamic response delay. The bar height was updated at every TR as a moving average of the current and up to the two available preceding values to reduce the bar fluctuation [139]. The consensus on the reporting and experimental design of clinical and cognitive-behavioral neurofeedback studies (CRED-nf checklist) is shown in Appendix B Table A1 [62].

Appendix B

Table A1. Consensus on the reporting and experimental design of clinical and cognitive-behavioral neurofeedback studies (CRED-nf) best practices checklist 2020.
Table A1. Consensus on the reporting and experimental design of clinical and cognitive-behavioral neurofeedback studies (CRED-nf) best practices checklist 2020.
DomainItem #Checklist ItemReported on Page #
Pre-experiment
1aPre-register experimental protocol and planned analysesn/a
1bJustify sample size4
Control groups
2aEmploy control group(s) or control condition(s)5
2bWhen leveraging experimental designs where a double-blind is possible, use a double-blindn/a
2cBlind those who rate the outcomes, and when possible, the statisticians involvedn/a
2dExamine to what extent participants and experimenters remain blindedn/a
2eIn clinical efficacy studies, employ a standard-of-care intervention group as a benchmark for improvementn/a
Control measures
3aCollect data on psychosocial factors5–6
3bReport whether participants were provided with a strategy5
3cReport the strategies participants usedn/a
3dReport methods used for online-data processing and artifact correction6–7
3eReport condition and group effects for artifactsn/a
Feedback specifications
4aReport how the online-feature extraction was defined6
4bReport and justify the reinforcement schedulen/a
4cReport the feedback modality and content4–5
4dCollect and report all brain activity variable(s) and/or contrasts used for feedback, as displayed to experimental participants4–5
4eReport the hardware and software used6–9
Outcome measures
Brain5aReport neurofeedback regulation success based on the feedback signaln/a
5bPlot within-session and between-session regulation blocks of feedback variable(s), as well as pre-to-post resting baselines or contrastsn/a
5cStatistically compare the experimental condition/group to the control condition(s)/group(s) (not only each group to baseline measures)n/a
Behavior6aInclude measures of clinical or behavioral significance, defined a priori, and describe whether they were reachedn/a
6bRun correlational analyses between regulation success and behavioral outcomesn/a
Data storage
7aUpload all materials, analysis scripts, code, and raw data used for analyses, as well as final values, to an open access data repository, when feasiblen/a
Note: Darker shaded boxes represent Essential checklist items; lightly shaded boxes represent Encouraged checklist items. We recommend using this checklist in conjunction with the standardized CRED-nf online tool (rtfin.org/CREDnf) and the CRED-nf article, which explains the motivation behind this checklist and provides details regarding many of the checklist items.
Table A2. Peak coordinates of the clusters from group fMRI analysis for the mean Focus-on-Breath vs. Describe contrast. Insular regions defined by the Brainnetome atlas [73].
Table A2. Peak coordinates of the clusters from group fMRI analysis for the mean Focus-on-Breath vs. Describe contrast. Insular regions defined by the Brainnetome atlas [73].
GyrusROILabel ID.LLabel ID.RAnatomical and Modified Cyto-Architectonic DescriptionsLeft Hemisphere MNI CoordinatesRight Hemisphere MNI Coordinates
Insular GyrusAnterior Insula165166vIa, ventral agranular insula−32, 14, −13 33, 14, −13
167168dIa, dorsal agranular insula−34, 18, 1 36, 18, 1
Mid-Insula169170vId/vIg, ventral dysgranular and granular insula−38, −4, −9 39, −2, −9
173174dId, dorsal dysgranular insula−38, 5, 5 38, 5, 5
Posterior Insula163164G, hypergranular insula−36, −20, 10 37, −18, 8
171172dIg, dorsal granular insula−38, −8, 8 39, −7, 8
Note. Anterior insula was obtained by combining left and right hemisphere, dorsal and ventral agranular insula. Mid-insula was obtained by combining left and right hemisphere, ventral dysgranular and granular insula, and dorsal dysgranular insula. Posterior insula was obtained by combining left and right hemisphere, hypergranular insula, and dorsal granular insula. See Figure 2a.
Table A3. Unadjusted means, standard deviations, effect sizes, and main analyses of parameter estimate (Focus-on-Breath vs. Rest) in insula subregions across runs.
Table A3. Unadjusted means, standard deviations, effect sizes, and main analyses of parameter estimate (Focus-on-Breath vs. Rest) in insula subregions across runs.
RunMeanSDEstimateSEtpCohen’s d
Anterior insula
 OBS−0.120.17
 NF-10.030.180.160.044.40<0.0010.74
 NF-2−0.020.170.100.042.74<0.010.46
 NF-3−0.030.230.090.042.59<0.050.43
 TRS−0.120.180.010.040.230.810.04
Mid-insula
 OBS−0.070.14
 NF-1−0.090.15−0.020.03−0.750.45−0.13
 NF-2−0.100.14−0.030.03−1.220.23−0.20
 NF-3−0.090.16−0.020.03−0.710.48−0.12
 TRS−0.110.14−0.040.03−1.530.13−0.26
Posterior insula
 OBS−0.080.13
 NF-1−0.200.160.030.031.110.270.19
 NF-2−0.170.150.070.032.25<0.050.38
 NF-3−0.130.130.120.034.16<0.0010.70
 TRS−0.090.150.110.033.56<0.0010.60
Abbreviations: OBS, Observe; NF, neurofeedback run; TRS, Transfer.
Table A4. Post hoc comparisons for parameter estimate (Focus-on-Breath vs. Rest) in insula subregions across runs.
Table A4. Post hoc comparisons for parameter estimate (Focus-on-Breath vs. Rest) in insula subregions across runs.
RunEstimateStd. Errorz Statisticp Value
Anterior insular cortex (aINS)
NF-1:OBS0.160.044.40<0.001
NF-2:OBS0.100.042.740.05
NF-3:OBS0.090.042.590.07
TR:OBS0.010.040.231.00
NF-2:NF-1−0.060.04−1.630.48
NF-3:NF-1−0.060.04−1.810.37
TR:NF-1−0.150.04−4.17<0.001
NF-3:NF-2−0.010.04−0.171.00
TR:NF-2−0.090.04−2.500.09
TR:NF-3−0.080.04−2.360.13
Mid-insular cortex (mINS)
NF-1:OBS−0.020.03−0.750.94
NF-2:OBS−0.030.03−1.220.74
NF-3:OBS−0.020.03−0.710.96
TR:OBS−0.040.03−1.530.54
NF-2:NF-1−0.010.03−0.470.99
NF-3:NF-10.000.030.041.00
TR:NF-1−0.020.03−0.780.94
NF-3:NF-20.010.030.520.99
TR:NF-2−0.010.03−0.301.00
TR:NF-3−0.020.03−0.820.92
Posterior insular cortex (pINS)
NF-1:OBS−0.120.03−4.16<0.001
NF-2:OBS−0.090.03−3.02<0.05
NF-3:OBS−0.060.03−1.920.31
TR:OBS−0.020.03−0.600.98
NF-2:NF-10.030.031.110.80
NF-3:NF-10.070.032.250.16
TR:NF-10.110.033.56<0.01
NF-3:NF-20.030.031.120.80
TR:NF-20.070.032.430.11
TR:NF-30.040.031.320.68
Abbreviations: NF, Neurofeedback; OBS, Observe; TR, Transfer.
Table A5. Correlations between parameter estimate (Focus-on-Breath vs. Describe) in insula subregions during neu-rofeedback runs and variables of interest.
Table A5. Correlations between parameter estimate (Focus-on-Breath vs. Describe) in insula subregions during neu-rofeedback runs and variables of interest.
Neurofeedback RunsObserve RunTransfer Run
aINSmINSpINSaINSmINSpINSaINSmINSpINS
PROMIS Positive Affect−0.2−0.15−0.1−0.24−0.16−0.190.080.150.16
PROMIS Meaning and Purpose−0.23−0.24−0.06−0.23−0.09−0.190.070.090.22
PROMIS Life Satisfaction−0.37 *−0.32−0.19−0.21−0.2−0.170.060.000.11
PROMIS Pain Interference0.20.180.290.120.030.150.150.090.12
PROMIS Pain Behavior0.180.320.33 *−0.09−0.02−0.040.130.200.18
PROMIS Fatigue0.280.30.230.270.210.290.060.050.10
Task Ratings Current Feeling−0.13−0.22−0.250.01−0.25−0.18−0.43 **−0.37 *−0.19
Task Ratings Mind Wander−0.31−0.36−0.24−0.11−0.24−0.34 *−0.20−0.28−0.22
Task Ratings Focus-on-Breath−0.03−0.14−0.18−0.14−0.17−0.040.110.260.16
Note. aINS: anterior insular cortex, mINS: mid-insular cortex, pINS: posterior insular cortex. PROMIS = Patient-Reported Outcomes Measurement Information System. Task rating scores were averaged across three neurofeedback runs.* p < 0.05. ** p < 0.01.
Figure A1. Barplot of the signal variance ratio (R^2 value) explained by the physiological noise regressors (RETROICOR) for the real-time processed and offline processed signals in the PCC region.
Figure A1. Barplot of the signal variance ratio (R^2 value) explained by the physiological noise regressors (RETROICOR) for the real-time processed and offline processed signals in the PCC region.
Brainsci 12 00363 g0a1
Figure A2. (ac), parameter estimate (Focus-on-Breath vs. Rest) across runs for insula subregions. Abbreviations: OBS, Observe; NF, neurofeedback run; TRS, Transfer * p < 0.05. ** p < 0.01. *** p < 0.001.
Figure A2. (ac), parameter estimate (Focus-on-Breath vs. Rest) across runs for insula subregions. Abbreviations: OBS, Observe; NF, neurofeedback run; TRS, Transfer * p < 0.05. ** p < 0.01. *** p < 0.001.
Brainsci 12 00363 g0a2
Figure A3. Scatterplot of significant correlations between parameter estimate (Focus-on-Breath vs. Describe) in insula subregions during neurofeedback runs and PROMIS measures.
Figure A3. Scatterplot of significant correlations between parameter estimate (Focus-on-Breath vs. Describe) in insula subregions during neurofeedback runs and PROMIS measures.
Brainsci 12 00363 g0a3

References

  1. Kabat-Zinn, J. Mindfulness. Mindfulness 2015, 6, 1481–1483. [Google Scholar] [CrossRef]
  2. Tomlinson, E.R.; Yousaf, O.; Vittersø, A.D.; Jones, L. Dispositional Mindfulness and Psychological Health: A Systematic Review. Mindfulness 2018, 9, 23–43. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Yu, X.; Langberg, J.M.; Becker, S.P. Dispositional mindfulness moderates the relation between brooding rumination and sleep problems in adolescents. Sleep Med. 2022, 90, 131–134. [Google Scholar] [CrossRef] [PubMed]
  4. Khoury, B.; Sharma, M.; Rush, S.E.; Fournier, C. Mindfulness-Based stress reduction for healthy individuals: A meta-analysis. J. Psychosom. Res. 2015, 78, 519–528. [Google Scholar] [CrossRef]
  5. Guo, D.; Sun, L.; Yu, X.; Liu, T.; Wu, L.; Sun, Z.; Zhang, F.; Zhou, Y.; Shen, M.; Liu, W. Mindfulness-based stress reduction improves the general health and stress of Chinese military recruits: A pilot study. Psychiatry Res. 2019, 281, 112571. [Google Scholar] [CrossRef]
  6. Blanck, P.; Perleth, S.; Heidenreich, T.; Kröger, P.; Ditzen, B.; Bents, H.; Mander, J. Effects of mindfulness exercises as stand-alone intervention on symptoms of anxiety and depression: Systematic review and meta-analysis. Behav. Res. Ther. 2018, 102, 25–35. [Google Scholar] [CrossRef]
  7. Goldberg, S.B.; Tucker, R.P.; Greene, P.A.; Davidson, R.J.; Wampold, B.E.; Kearney, D.J.; Simpson, T.L. Mindfulness-Based interventions for psychiatric disorders: A systematic review and meta-analysis. Clin. Psychol. Rev. 2018, 59, 52–60. [Google Scholar] [CrossRef]
  8. Raj, S.; Ghosh, D.; Verma, S.K.; Singh, T. The mindfulness trajectories of addressing suicidal behaviour: A systematic review. Int. J. Soc. Psychiatry 2021, 67, 507–519. [Google Scholar] [CrossRef]
  9. Essau, C.A.; Lewinsohn, P.M.; Lim, J.X.; Moon-ho, R.H.; Rohde, P. Incidence, recurrence and comorbidity of anxiety disorders in four major developmental stages. J. Affect. Disord. 2018, 228, 248–253. [Google Scholar] [CrossRef] [Green Version]
  10. Keyes, K.M.; Gary, D.; O’Malley, P.M.; Hamilton, A.; Schulenberg, J. Recent increases in depressive symptoms among US adolescents: Trends from 1991 to 2018. Soc. Psychiatry Psychiatr. Epidemiol. 2019, 54, 987–996. [Google Scholar] [CrossRef]
  11. Burstein, B.; Agostino, H.; Greenfield, B. Suicidal attempts and ideation among children and adolescents in US emergency departments, 2007–2015. JAMA Pediatrics 2019, 173, 598–600. [Google Scholar] [CrossRef] [PubMed]
  12. Sanger, K.L.; Thierry, G.; Dorjee, D. Effects of school-based mindfulness training on emotion processing and well-being in adolescents: Evidence from event-related potentials. Dev. Sci. 2018, 21, e12646. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Tan, L.B.; Martin, G. Mind full or mindful: A report on mindfulness and psychological health in healthy adolescents. Int. J. Adolesc. Youth 2016, 21, 64–74. [Google Scholar] [CrossRef]
  14. Gibson, J. Mindfulness, Interoception, and the Body: A Contemporary Perspective. Front. Psychol. 2019, 10, 2012. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Haase, L.; Thom, N.J.; Shukla, A.; Davenport, P.W.; Simmons, A.N.; Stanley, E.A.; Paulus, M.P.; Johnson, D.C. Mindfulness-based training attenuates insula response to an aversive interoceptive challenge. Soc. Cogn. Affect. Neurosci. 2015, 11, 182–190. [Google Scholar] [CrossRef] [Green Version]
  16. De Jong, M.; Lazar, S.W.; Hug, K.; Mehling, W.E.; Hölzel, B.K.; Sack, A.T.; Peeters, F.; Ashih, H.; Mischoulon, D.; Gard, T. Effects of mindfulness-based cognitive therapy on body awareness in patients with chronic pain and comorbid depression. Front. Psychol. 2016, 7, 967. [Google Scholar] [CrossRef] [Green Version]
  17. Bornemann, B.; Singer, T. Taking time to feel our body: Steady increases in heartbeat perception accuracy and decreases in alexithymia over 9 months of contemplative mental training. Psychophysiology 2017, 54, 469–482. [Google Scholar] [CrossRef]
  18. Khalsa, S.S.; Adolphs, R.; Cameron, O.G.; Critchley, H.D.; Davenport, P.W.; Feinstein, J.S.; Feusner, J.D.; Garfinkel, S.N.; Lane, R.D.; Mehling, W.E.; et al. Interoception and Mental Health: A Roadmap. Biol. Psychiatry Cogn. Neurosci. Neuroimaging 2018, 3, 501–513. [Google Scholar] [CrossRef]
  19. Hölzel, B.K.; Lazar, S.W.; Gard, T.; Schuman-Olivier, Z.; Vago, D.R.; Ott, U. How does mindfulness meditation work? Proposing mechanisms of action from a conceptual and neural perspective. Perspect. Psychol. Sci. 2011, 6, 537–559. [Google Scholar] [CrossRef]
  20. Farb, N.; Daubenmier, J.; Price, C.J.; Gard, T.; Kerr, C.; Dunn, B.D.; Klein, A.C.; Paulus, M.P.; Mehling, W.E. Interoception, contemplative practice, and health. Front. Psychol. 2015, 6, 763. [Google Scholar] [CrossRef] [Green Version]
  21. Mehling, W.E.; Price, C.; Daubenmier, J.J.; Acree, M.; Bartmess, E.; Stewart, A. The multidimensional assessment of interoceptive awareness (MAIA). PLoS ONE 2012, 7, e48230. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Mehling, W.E.; Chesney, M.A.; Metzler, T.J.; Goldstein, L.A.; Maguen, S.; Geronimo, C.; Agcaoili, G.; Barnes, D.E.; Hlavin, J.A.; Neylan, T.C. A 12-week integrative exercise program improves self-reported mindfulness and interoceptive awareness in war veterans with posttraumatic stress symptoms. J. Clin. Psychol. 2018, 74, 554–565. [Google Scholar] [CrossRef] [PubMed]
  23. Fissler, M.; Winnebeck, E.; Schroeter, T.; Gummersbach, M.; Huntenburg, J.M.; Gaertner, M.; Barnhofer, T. An investigation of the effects of brief mindfulness training on self-reported interoceptive awareness, the ability to decenter, and their role in the reduction of depressive symptoms. Mindfulness 2016, 7, 1170–1181. [Google Scholar] [CrossRef]
  24. Fogel, A. The Psychophysiology of Self-Awareness: Rediscovering the Lost Art of Body Sense; WW Norton & Co.: New York, NY, USA, 2009. [Google Scholar]
  25. Craig, A.D.; Craig, A.D. How do you feel? Interoception: The sense of the physiological condition of the body. Nat. Rev. Neurosci. 2002, 3, 655–666. [Google Scholar] [CrossRef]
  26. Critchley, H.D.; Wiens, S.; Ohman, P.R.A.; Dolan, R.J. Neural systems supporting interoceptive awareness. Psychol. Sect. Karolinska Hosp. 2004, 7, 189–195. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Gogolla, N. The insular cortex. Curr. Biol. 2017, 27, R580–R586. [Google Scholar] [CrossRef] [Green Version]
  28. Khalsa, S.S.; Rudrauf, D.; Feinstein, J.S.; Tranel, D. The pathways of interoceptive awareness. Nat. Neurosci. 2009, 12, 1494–1496. [Google Scholar] [CrossRef] [Green Version]
  29. Khalsa, S.S.; Rudrauf, D.; Damasio, A.R.; Davidson, R.J.; Lutz, A.; Tranel, D. Interoceptive awareness in experienced meditators. Psychophysiology 2008, 45, 671–677. [Google Scholar] [CrossRef] [Green Version]
  30. Paulus, M.P.; Stein, M.B. An insular view of anxiety. Biol. Psychiatry 2006, 60, 383–387. [Google Scholar] [CrossRef]
  31. Ives-Deliperi, V.L.; Solms, M.; Meintjes, E.M. The neural substrates of mindfulness: An fMRI investigation. Soc. Neurosci. 2011, 6, 231–242. [Google Scholar] [CrossRef]
  32. Kabat-Zinn, J. Full Catastrophe Living, Revised Edition: How to Cope With Stress, Pain and Illness Using Mindfulness Meditation; Hachette: Paris, France, 2013. [Google Scholar]
  33. Gotink, R.A.; Meijboom, R.; Vernooij, M.W.; Smits, M.; Hunink, M.G. 8-Week Mindfulness Based Stress Reduction induces brain changes similar to traditional long-term meditation practice—A systematic review. Brain Cogn. 2016, 108, 32–41. [Google Scholar] [CrossRef] [PubMed]
  34. Young, K.S.; van der Velden, A.M.; Craske, M.G.; Pallesen, K.J.; Fjorback, L.; Roepstorff, A.; Parsons, C.E. The impact of Mindfulness-Based Interventions on brain activity: A systematic review of functional magnetic resonance imaging studies. Neurosci. Biobehav. Rev. 2017, 84, 424–433. [Google Scholar] [CrossRef] [PubMed]
  35. Hölzel, B.K.; Ott, U.; Hempel, H.; Hackl, A.; Wolf, K.; Stark, R.; Vaitl, D. Differential engagement of anterior cingulate and adjacent medial frontal cortex in adept meditators and non-meditators. Neurosci. Lett. 2007, 421, 16–21. [Google Scholar] [CrossRef] [PubMed]
  36. Tomasino, B.; Fabbro, F. Increases in the right dorsolateral prefrontal cortex and decreases the rostral prefrontal cortex activation after-8 weeks of focused attention based mindfulness meditation. Brain Cogn. 2016, 102, 46–54. [Google Scholar] [CrossRef]
  37. Zeidan, F.; Emerson, N.M.; Farris, S.R.; Ray, J.N.; Jung, Y.; McHaffie, J.G.; Coghill, R.C. Mindfulness meditation-based pain relief employs different neural mechanisms than placebo and sham mindfulness meditation-induced analgesia. J. Neurosci. 2015, 35, 15307–15325. [Google Scholar] [CrossRef]
  38. Tang, Y.-Y.; Hölzel, B.K.; Posner, M.I. The neuroscience of mindfulness meditation. Nat. Rev. Neurosci. 2015, 16, 213. [Google Scholar] [CrossRef]
  39. Kabat-Zinn, J. Wherever You Go, There You Are: Mindfulness Meditation in Everyday Life; Hachette Books: New York, NY, USA, 2009. [Google Scholar]
  40. Zsadanyi, S.E.; Kurth, F.; Luders, E. The Effects of Mindfulness and Meditation on the Cingulate Cortex in the Healthy Human Brain: A Review. Mindfulness 2021, 12, 2371–2387. [Google Scholar] [CrossRef]
  41. Buckner, R.L.; Andrews-Hanna, J.R.; Schacter, D.L. The brain’s default network: Anatomy, function, and relevance to disease. Ann. N. Y. Acad. Sci. 2008, 1124, 1–38. [Google Scholar] [CrossRef] [Green Version]
  42. Whitfield-Gabrieli, S.; Moran, J.M.; Nieto-Castañón, A.; Triantafyllou, C.; Saxe, R.; Gabrieli, J.D. Associations and dissociations between default and self-reference networks in the human brain. Neuroimage 2011, 55, 225–232. [Google Scholar] [CrossRef]
  43. Brewer, J.A.; Worhunsky, P.D.; Gray, J.R.; Tang, Y.-Y.; Weber, J.; Kober, H. Meditation experience is associated with differences in default mode network activity and connectivity. Proc. Natl. Acad. Sci. USA 2011, 108, 20254–20259. [Google Scholar] [CrossRef] [Green Version]
  44. Gusnard, D.A.; Akbudak, E.; Shulman, G.L.; Raichle, M.E. Medial prefrontal cortex and self-referential mental activity: Relation to a default mode of brain function. Proc. Natl. Acad. Sci. USA 2001, 98, 4259–4264. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Raichle, M.E.; MacLeod, A.M.; Snyder, A.Z.; Powers, W.J.; Gusnard, D.A.; Shulman, G.L. A default mode of brain function. Proc. Natl. Acad. Sci. USA 2001, 98, 676–682. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Longarzo, M.; D’Olimpio, F.; Chiavazzo, A.; Santangelo, G.; Trojano, L.; Grossi, D. The relationships between interoception and alexithymic trait. The Self-Awareness Questionnaire in healthy subjects. Front. Psychol. 2015, 6, 1149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Grossi, D.; Di Vita, A.; Palermo, L.; Sabatini, U.; Trojano, L.; Guariglia, C. The brain network for self-feeling: A symptom-lesion mapping study. Neuropsychologia 2014, 63, 92–98. [Google Scholar] [CrossRef]
  48. Garrison, K.A.; Scheinost, D.; Worhunsky, P.D.; Elwafi, H.M.; Thornhill, T.A.; Thompson, E.; Saron, C.; Desbordes, G.; Kober, H.; Hampson, M. Real-Time fMRI links subjective experience with brain activity during focused attention. Neuroimage 2013, 81, 110–118. [Google Scholar] [CrossRef] [Green Version]
  49. Kirlic, N.; Cohen, Z.P.; Tsuchiyagaito, A.; Misaki, M.; McDermott, T.J.; Aupperle, R.L.; Stewart, J.L.; Singh, M.K.; Paulus, M.P.; Bodurka, J. Self-Regulation of the Posterior Cingulate Cortex with Real-Time fMRI Neurofeedback in Healthy Adolescents. Cogn. Affect. Behav. Neurosci. 2021; in press. [Google Scholar]
  50. Takesian, A.E.; Hensch, T.K. Balancing plasticity/stability across brain development. Prog. Brain Res. 2013, 207, 3–34. [Google Scholar]
  51. Fuhrmann, D.; Knoll, L.J.; Blakemore, S.-J. Adolescence as a sensitive period of brain development. Trends Cogn. Sci. 2015, 19, 558–566. [Google Scholar] [CrossRef] [Green Version]
  52. Tanay, G.; Bernstein, A. State Mindfulness Scale (SMS): Development and initial validation. Psychol. Assess. 2013, 25, 1286. [Google Scholar] [CrossRef] [Green Version]
  53. Cella, D.; Riley, W.; Stone, A.; Rothrock, N.; Reeve, B.; Yount, S.; Amtmann, D.; Bode, R.; Buysse, D.; Choi, S. The Patient-Reported Outcomes Measurement Information System (PROMIS) developed and tested its first wave of adult self-reported health outcome item banks: 2005–2008. J. Clin. Epidemiol. 2010, 63, 1179–1194. [Google Scholar] [CrossRef] [Green Version]
  54. Forrest, C.B.; Bevans, K.B.; Filus, A.; Devine, J.; Becker, B.D.; Carle, A.C.; Teneralli, R.E.; Moon, J.; Ravens-Sieberer, U. Assessing children’s eudaimonic well-being: The PROMIS pediatric meaning and purpose item banks. J. Pediatr. Psychol. 2019, 44, 1074–1082. [Google Scholar] [CrossRef] [PubMed]
  55. Forrest, C.B.; Ravens-Sieberer, U.; Devine, J.; Becker, B.D.; Teneralli, R.E.; Moon, J.; Carle, A.C.; Tucker, C.A.; Bevans, K.B. Development and evaluation of the PROMIS® pediatric positive affect item bank, child-report and parent-proxy editions. J. Happiness Stud. 2018, 19, 699–718. [Google Scholar] [CrossRef] [PubMed]
  56. Cunningham, N.R.; Kashikar-Zuck, S.; Mara, C.; Goldschneider, K.R.; Revicki, D.A.; Dampier, C.; Sherry, D.D.; Crosby, L.; Carle, A.; Cook, K.F. Development and validation of the self-reported PROMIS pediatric pain behavior item bank and short form scale. Pain 2017, 158, 1323. [Google Scholar] [CrossRef]
  57. Amtmann, D.; Cook, K.F.; Jensen, M.P.; Chen, W.-H.; Choi, S.; Revicki, D.; Cella, D.; Rothrock, N.; Keefe, F.; Callahan, L. Development of a PROMIS item bank to measure pain interference. Pain 2010, 150, 173–182. [Google Scholar] [CrossRef] [Green Version]
  58. Lai, J.-S.; Stucky, B.D.; Thissen, D.; Varni, J.W.; DeWitt, E.M.; Irwin, D.E.; Yeatts, K.B.; DeWalt, D.A. Development and psychometric properties of the PROMIS® pediatric fatigue item banks. Qual. Life Res. 2013, 22, 2417–2427. [Google Scholar] [CrossRef] [Green Version]
  59. Rudkin, E.; Medvedev, O.N.; Siegert, R.J. The five-facet mindfulness questionnaire: Why the observing subscale does not predict psychological symptoms. Mindfulness 2018, 9, 230–242. [Google Scholar] [CrossRef]
  60. Desmond, J.E.; Glover, G.H. Estimating sample size in functional MRI (fMRI) neuroimaging studies: Statistical power analyses. J. Neurosci. Methods 2002, 118, 115–128. [Google Scholar] [CrossRef] [Green Version]
  61. Kelley, W.M.; Macrae, C.N.; Wyland, C.L.; Caglar, S.; Inati, S.; Heatherton, T.F. Finding the self? An event-related fMRI study. J. Cogn. Neurosci. 2002, 14, 785–794. [Google Scholar] [CrossRef] [Green Version]
  62. Ros, T.; Enriquez-Geppert, S.; Zotev, V.; Young, K.D.; Wood, G.; Whitfield-Gabrieli, S.; Wan, F.; Vuilleumier, P.; Vialatte, F.; Van De Ville, D. Consensus on the Reporting and Experimental Design of Clinical and Cognitive-Behavioural Neurofeedback Studies (CRED-nf Checklist); Oxford University Press: Oxford, UK, 2020. [Google Scholar]
  63. Pruessmann, K.P.; Weiger, M.; Scheidegger, M.B.; Boesiger, P. SENSE: Sensitivity encoding for fast MRI. Magn. Reson. Med. Off. J. Int. Soc. Magn. Reson. Med. 1999, 42, 952–962. [Google Scholar] [CrossRef]
  64. Bellgowan, P.S.; Bandettini, P.A.; van Gelderen, P.; Martin, A.; Bodurka, J. Improved BOLD detection in the medial temporal region using parallel imaging and voxel volume reduction. Neuroimage 2006, 29, 1244–1251. [Google Scholar] [CrossRef]
  65. Bodurka, J.; Ye, F.; Petridou, N.; Murphy, K.; Bandettini, P.A. Mapping the MRI voxel volume in which thermal noise matches physiological noise—Implications for fMRI. Neuroimage 2007, 34, 542–549. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Cox, R.W. AFNI: Software for analysis and visualization of functional magnetic resonance neuroimages. Comput. Biomed. Res. Int. J. 1996, 29, 162–173. [Google Scholar] [CrossRef] [PubMed]
  67. Glover, G.H.; Li, T.Q.; Ress, D. Image-Based method for retrospective correction of physiological motion effects in fMRI: RETROICOR. Magn. Reson. Med. Off. J. Int. Soc. Magn. Reson. Med. 2000, 44, 162–167. [Google Scholar] [CrossRef]
  68. Birn, R.M.; Smith, M.A.; Jones, T.B.; Bandettini, P.A. The respiration response function: The temporal dynamics of fMRI signal fluctuations related to changes in respiration. Neuroimage 2008, 40, 644–654. [Google Scholar] [CrossRef] [Green Version]
  69. Avants, B.B.; Epstein, C.L.; Grossman, M.; Gee, J.C. Symmetric diffeomorphic image registration with cross-correlation: Evaluating automated labeling of elderly and neurodegenerative brain. Med. Image Anal. 2008, 12, 26–41. [Google Scholar] [CrossRef] [Green Version]
  70. Jo, H.J.; Saad, Z.S.; Simmons, W.K.; Milbury, L.A.; Cox, R.W. Mapping sources of correlation in resting state FMRI, with artifact detection and removal. Neuroimage 2010, 52, 571–582. [Google Scholar] [CrossRef] [Green Version]
  71. Lindquist, M.A.; Geuter, S.; Wager, T.D.; Caffo, B.S. Modular preprocessing pipelines can reintroduce artifacts into fMRI data. Hum. Brain Mapp. 2019, 40, 2358–2376. [Google Scholar] [CrossRef] [Green Version]
  72. Misaki, M.; Bodurka, J. The impact of real-time fMRI denoising on online evaluation of brain activity and functional connectivity. J. Neural Eng. 2021, 18, 046092. [Google Scholar] [CrossRef]
  73. Fan, L.; Li, H.; Zhuo, J.; Zhang, Y.; Wang, J.; Chen, L.; Yang, Z.; Chu, C.; Xie, S.; Laird, A.R. The human brainnetome atlas: A new brain atlas based on connectional architecture. Cereb. Cortex 2016, 26, 3508–3526. [Google Scholar] [CrossRef]
  74. RCoreTeam. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria, 2013. [Google Scholar]
  75. Revelle, W.R. Psych: Procedures for Personality and Psychological Research; Northwestern University Press: Evanston, IL, USA, 2017. [Google Scholar]
  76. Bates, D.; Sarkar, D.; Bates, M.D.; Matrix, L. The lme4 Package. R Package Version. 2007, Volume 2, p. 74. Available online: https://cran.r-project.org/web/packages/lme4/lme4.pdf (accessed on 2 March 2022).
  77. Hothorn, T.; Bretz, F.; Westfall, P. Simultaneous inference in general parametric models. Biom. J. J. Math. Methods Biosci. 2008, 50, 346–363. [Google Scholar] [CrossRef] [Green Version]
  78. Kurth, F.; Zilles, K.; Fox, P.T.; Laird, A.R.; Eickhoff, S.B. A link between the systems: Functional differentiation and integration within the human insula revealed by meta-analysis. Brain Struct. Funct. 2010, 214, 519–534. [Google Scholar] [CrossRef] [PubMed]
  79. Allen, M.; Dietz, M.; Blair, K.S.; van Beek, M.; Rees, G.; Vestergaard-Poulsen, P.; Lutz, A.; Roepstorff, A. Cognitive-Affective neural plasticity following active-controlled mindfulness intervention. J. Neurosci. Off. J. Soc. Neurosci. 2012, 32, 15601–15610. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Kirk, U.; Gu, X.; Sharp, C.; Hula, A.; Fonagy, P.; Montague, P.R. Mindfulness training increases cooperative decision making in economic exchanges: Evidence from fMRI. NeuroImage 2016, 138, 274–283. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Sevinc, G.; Hölzel, B.K.; Hashmi, J.; Greenberg, J.; McCallister, A.; Treadway, M.; Schneider, M.L.; Dusek, J.A.; Carmody, J.; Lazar, S.W. Common and dissociable neural activity after mindfulness-based stress reduction and relaxation response programs. Psychosom. Med. 2018, 80, 439. [Google Scholar] [CrossRef] [Green Version]
  82. Holzel, B.K.; Hoge, E.A.; Greve, D.N.; Gard, T.; Creswell, J.D.; Brown, K.W.; Barrett, L.F.; Schwartz, C.; Vaitl, D.; Lazar, S.W. Neural mechanisms of symptom improvements in generalized anxiety disorder following mindfulness training. Neuroimage Clin. 2013, 2, 448–458. [Google Scholar] [CrossRef] [Green Version]
  83. Haase, L.; May, A.C.; Falahpour, M.; Isakovic, S.; Simmons, A.N.; Hickman, S.D.; Liu, T.T.; Paulus, M.P. A pilot study investigating changes in neural processing after mindfulness training in elite athletes. Front. Behav. Neurosci. 2015, 9, 229. [Google Scholar] [CrossRef] [Green Version]
  84. Azzalini, D.; Rebollo, I.; Tallon-Baudry, C. Visceral signals shape brain dynamics and cognition. Trends Cogn. Sci. 2019, 23, 488–509. [Google Scholar] [CrossRef] [Green Version]
  85. Uddin, L.Q.; Nomi, J.S.; Hébert-Seropian, B.; Ghaziri, J.; Boucher, O. Structure and function of the human insula. J. Clin. Neurophysiol. Off. Publ. Am. Electroencephalogr. Soc. 2017, 34, 300. [Google Scholar] [CrossRef]
  86. Kleckner, I.R.; Zhang, J.; Touroutoglou, A.; Chanes, L.; Xia, C.; Simmons, W.K.; Quigley, K.S.; Dickerson, B.C.; Barrett, L.F. Evidence for a large-scale brain system supporting allostasis and interoception in humans. Nat. Hum. Behav. 2017, 1, 0069. [Google Scholar] [CrossRef] [Green Version]
  87. Sitaram, R.; Ros, T.; Stoeckel, L.; Haller, S.; Scharnowski, F.; Lewis-Peacock, J.; Weiskopf, N.; Blefari, M.L.; Rana, M.; Oblak, E. Closed-Loop brain training: The science of neurofeedback. Nat. Rev. Neurosci. 2017, 18, 86–100. [Google Scholar] [CrossRef] [Green Version]
  88. Vassena, E.; Silvetti, M.; Boehler, C.N.; Achten, E.; Fias, W.; Verguts, T. Overlapping neural systems represent cognitive effort and reward anticipation. PLoS ONE 2014, 9, e91008. [Google Scholar] [CrossRef] [PubMed]
  89. Westbrook, A.; Braver, T.S. Cognitive effort: A neuroeconomic approach. Cogn. Affect. Behav. Neurosci. 2015, 15, 395–415. [Google Scholar] [CrossRef] [PubMed]
  90. Krause, F.; Kogias, N.; Krentz, M.; Luhrs, M.; Goebel, R.; Hermans, E.J. Self-Regulation of stress-related large-scale brain network balance using real-time fMRI neurofeedback. Neuroimage 2021, 243, 118527. [Google Scholar] [CrossRef] [PubMed]
  91. Carter, C.S.; Macdonald, A.M.; Botvinick, M.; Ross, L.L.; Stenger, V.A.; Noll, D.; Cohen, J.D. Parsing executive processes: Strategic vs. evaluative functions of the anterior cingulate cortex. Proc. Natl. Acad. Sci. USA 2000, 97, 1944–1948. [Google Scholar] [CrossRef] [Green Version]
  92. Kerns, J.G.; Cohen, J.D.; MacDonald, A.W.; Cho, R.Y.; Stenger, V.A.; Carter, C.S. Anterior cingulate conflict monitoring and adjustments in control. Science 2004, 303, 1023–1026. [Google Scholar] [CrossRef] [Green Version]
  93. Craig, A.D. How do you feel—Now? The anterior insula and human awareness. Nat. Rev. Neurosci. 2009, 10, 59–70. [Google Scholar] [CrossRef]
  94. Ninaus, M.; Kober, S.; Witte, M.; Koschutnig, K.; Stangl, M.; Neuper, C.; Wood, G. Neural substrates of cognitive control under the belief of getting neurofeedback training. Front. Hum. Neurosci. 2013, 7, 914. [Google Scholar] [CrossRef] [Green Version]
  95. Rutledge, R.B.; Skandali, N.; Dayan, P.; Dolan, R.J. A computational and neural model of momentary subjective well-being. Proc. Natl. Acad. Sci. USA 2014, 111, 12252–12257. [Google Scholar] [CrossRef] [Green Version]
  96. Nardo, D.; Högberg, G.; Flumeri, F.; Jacobsson, H.; Larsson, S.; Hällström, T.; Pagani, M. Self-Rating scales assessing subjective well-being and distress correlate with rCBF in PTSD-sensitive regions. Psychol. Med. 2011, 41, 2549–2561. [Google Scholar] [CrossRef]
  97. Cabeen, R.P.; Toga, A.W.; Allman, J.M. Frontoinsular cortical microstructure is linked to life satisfaction in young adulthood. Brain Imaging Behav. 2021, 15, 2775–2789. [Google Scholar] [CrossRef]
  98. Li, R.; Zhu, X.; Zheng, Z.; Wang, P.; Li, J. Subjective well-being is associated with the functional connectivity network of the dorsal anterior insula. Neuropsychologia 2020, 141, 107393. [Google Scholar] [CrossRef] [PubMed]
  99. Jones, T.G.; Rapport, L.J.; Hanks, R.A.; Lichtenberg, P.A.; Telmet, K. Cognitive and psychosocial predictors of subjective well-being in urban older adults. Clin. Neuropsychol. 2003, 17, 3–18. [Google Scholar] [CrossRef] [PubMed]
  100. Lewis, G.J.; Kanai, R.; Rees, G.; Bates, T.C. Neural correlates of the ‘good life’: Eudaimonic well-being is associated with insular cortex volume. Soc. Cogn. Affect. Neurosci. 2014, 9, 615–618. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  101. Cash, M.; Whittingham, K. What facets of mindfulness contribute to psychological well-being and depressive, anxious, and stress-related symptomatology? Mindfulness 2010, 1, 177–182. [Google Scholar] [CrossRef]
  102. Farb, N.; Segal, Z.; Anderson, A. Mindfulness meditation training alters cortical representations of interoceptive attention. Soc. Cogn. Affect. Neurosci. 2013, 8, 15–26. [Google Scholar] [CrossRef] [Green Version]
  103. Segerdahl, A.R.; Mezue, M.; Okell, T.W.; Farrar, J.T.; Tracey, I. The dorsal posterior insula is not an island in pain but subserves a fundamental role-Response to:“Evidence against pain specificity in the dorsal posterior insula” by Davis et al. F1000Research 2015, 4, 1207. [Google Scholar] [CrossRef] [Green Version]
  104. Bergeron, D.; Obaid, S.; Fournier-Gosselin, M.P.; Bouthillier, A.; Nguyen, D.K. Deep Brain Stimulation of the Posterior Insula in Chronic Pain: A Theoretical Framework. Brain Sci. 2021, 11, 639. [Google Scholar] [CrossRef]
  105. Farb, N.; Segal, Z.V.; Mayberg, H.; Bean, J.; McKeon, D.; Fatima, Z.; Anderson, A.K. Attending to the present: Mindfulness meditation reveals distinct neural modes of self-reference. Soc. Cogn. Affect. Neurosci. 2007, 2, 313–322. [Google Scholar] [CrossRef]
  106. Grant, J.A.; Courtemanche, J.; Rainville, P. A non-elaborative mental stance and decoupling of executive and pain-related cortices predicts low pain sensitivity in Zen meditators. Pain 2011, 152, 150–156. [Google Scholar] [CrossRef]
  107. Gard, T.; Hölzel, B.K.; Sack, A.T.; Hempel, H.; Lazar, S.W.; Vaitl, D.; Ott, U. Pain attenuation through mindfulness is associated with decreased cognitive control and increased sensory processing in the brain. Cereb. Cortex 2012, 22, 2692–2702. [Google Scholar] [CrossRef] [Green Version]
  108. Kirlic, N.; Aupperle, R.L.; Rhudy, J.L.; Misaki, M.; Kuplicki, R.; Sutton, A.; Alvarez, R.P. Latent variable analysis of negative affect and its contributions to neural responses during shock anticipation. Neuropsychopharmacology 2019, 44, 695–702. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  109. Paulus, M.P.; Khalsa, S. When you don’t feel right inside: Homeostatic dysregulation and the mid-insular cortex in psychiatric disorders. Am. Psychiatric. Assoc. 2021, 178, 683–685. [Google Scholar] [CrossRef] [PubMed]
  110. Jakab, A.; Molnár, P.P.; Bogner, P.; Béres, M.; Berényi, E.L. Connectivity-Based parcellation reveals interhemispheric differences in the insula. Brain Topogr. 2012, 25, 264–271. [Google Scholar] [CrossRef] [PubMed]
  111. Ghaziri, J.; Tucholka, A.; Girard, G.; Houde, J.-C.; Boucher, O.; Gilbert, G.; Descoteaux, M.; Lippé, S.; Rainville, P.; Nguyen, D.K. The corticocortical structural connectivity of the human insula. Cereb. Cortex 2017, 27, 1216–1228. [Google Scholar] [CrossRef] [Green Version]
  112. Cerliani, L.; Thomas, R.M.; Jbabdi, S.; Siero, J.C.; Nanetti, L.; Crippa, A.; Gazzola, V.; D’Arceuil, H.; Keysers, C. Probabilistic tractography recovers a rostrocaudal trajectory of connectivity variability in the human insular cortex. Hum. Brain Mapp. 2012, 33, 2005–2034. [Google Scholar] [CrossRef]
  113. Uddin, L.Q.; Supekar, K.S.; Ryali, S.; Menon, V. Dynamic reconfiguration of structural and functional connectivity across core neurocognitive brain networks with development. J. Neurosci. 2011, 31, 18578–18589. [Google Scholar] [CrossRef]
  114. Apkarian, A.V.; Bushnell, M.C.; Treede, R.-D.; Zubieta, J.-K. Human brain mechanisms of pain perception and regulation in health and disease. Eur. J. Pain 2005, 9, 463–484. [Google Scholar] [CrossRef]
  115. Taylor, K.S.; Seminowicz, D.A.; Davis, K.D. Two systems of resting state connectivity between the insula and cingulate cortex. Hum. Brain Mapp. 2009, 30, 2731–2745. [Google Scholar] [CrossRef]
  116. Vogt, B.A. Cingulate cortex in the three limbic subsystems. Handb. Clin. Neurol. 2019, 166, 39–51. [Google Scholar]
  117. Vogt, B.A.; Vogt, L.; Farber, N.B.; Bush, G. Architecture and neurocytology of monkey cingulate gyrus. J. Comp. Neurol. 2005, 485, 218–239. [Google Scholar] [CrossRef] [Green Version]
  118. Maldonado, I.L.; de Matos, V.P.; Cuesta, T.A.C.; Herbet, G.; Destrieux, C. The human cingulum: From the limbic tract to the connectionist paradigm. Neuropsychologia 2020, 144, 107487. [Google Scholar] [CrossRef] [PubMed]
  119. Gehrlach, D.A.; Dolensek, N.; Klein, A.S.; Chowdhury, R.R.; Matthys, A.; Junghänel, M.; Gaitanos, T.N.; Podgornik, A.; Black, T.D.; Vaka, N.R. Aversive state processing in the posterior insular cortex. Nat. Neurosci. 2019, 22, 1424–1437. [Google Scholar] [CrossRef] [PubMed]
  120. Bishop, S.R.; Lau, M.; Shapiro, S.; Carlson, L.; Anderson, N.D.; Carmody, J.; Segal, Z.V.; Abbey, S.; Speca, M.; Velting, D.; et al. Mindfulness: A Proposed Operational Definition. Clin. Psychol. Sci. Pract. 2004, 11, 230–241. [Google Scholar] [CrossRef]
  121. Kabat-Zinn, J. An outpatient program in behavioral medicine for chronic pain patients based on the practice of mindfulness meditation: Theoretical considerations and preliminary results. Gen. Hosp. Psychiatry 1982, 4, 33–47. [Google Scholar] [CrossRef]
  122. Lindsay, E.K.; Creswell, J.D. Mechanisms of mindfulness training: Monitor and Acceptance Theory (MAT). Clin. Psychol. Rev. 2017, 51, 48–59. [Google Scholar] [CrossRef] [Green Version]
  123. Landtblom, A.-M.; Lindehammar, H.; Karlsson, H. Insular cortex activation in a patient with “sensed presence”/ecstatic seizures. Epilepsy Behav. 2011, 20, 714–718. [Google Scholar] [CrossRef]
  124. Chambers, R.; Gullone, E.; Allen, N.B. Mindful emotion regulation: An integrative review. Clin. Psychol. Rev. 2009, 29, 560–572. [Google Scholar] [CrossRef]
  125. Chen, W.G.; Schloesser, D.; Arensdorf, A.M.; Simmons, J.M.; Cui, C.; Valentino, R.; Gnadt, J.W.; Nielsen, L.; Hillaire-Clarke, C.S.; Spruance, V. The emerging science of interoception: Sensing, integrating, interpreting, and regulating signals within the self. Trends Neurosci. 2021, 44, 3–16. [Google Scholar] [CrossRef]
  126. Li, D.; Zucker, N.L.; Kragel, P.A.; Covington, V.E.; LaBar, K.S. Adolescent development of insula-dependent interoceptive regulation. Dev. Sci. 2017, 20, e12438. [Google Scholar] [CrossRef] [Green Version]
  127. Murphy, J.; Brewer, R.; Catmur, C.; Bird, G. Interoception and psychopathology: A developmental neuroscience perspective. Dev. Cogn. Neurosci. 2017, 23, 45–56. [Google Scholar] [CrossRef]
  128. May, A.C.; Stewart, J.L.; Paulus, M.P.; Tapert, S.F. The effect of age on neural processing of pleasant soft touch stimuli. Front. Behav. Neurosci. 2014, 8, 52. [Google Scholar] [CrossRef] [Green Version]
  129. Van Hoorn, J.; Crone, E.A.; Van Leijenhorst, L. Hanging out with the right crowd: Peer influence on risk-taking behavior in adolescence. J. Res. Adolesc. 2017, 27, 189–200. [Google Scholar] [CrossRef] [Green Version]
  130. Smith, D.G.; Xiao, L.; Bechara, A. Decision making in children and adolescents: Impaired Iowa Gambling Task performance in early adolescence. Dev. Psychol. 2012, 48, 1180. [Google Scholar] [CrossRef]
  131. Northoff, G.; Heinzel, A.; De Greck, M.; Bermpohl, F.; Dobrowolny, H.; Panksepp, J. Self-Referential processing in our brain—A meta-analysis of imaging studies on the self. Neuroimage 2006, 31, 440–457. [Google Scholar] [CrossRef]
  132. Misaki, M.; Barzigar, N.; Zotev, V.; Phillips, R.; Cheng, S.; Bodurka, J. Real-Time fMRI processing with physiological noise correction—Comparison with off-line analysis. J. Neurosci. Methods 2015, 256, 117–121. [Google Scholar] [CrossRef] [Green Version]
  133. Wang, S.; Taren, A.A.; Smith, D.V. Functional Parcellation of the Default Mode Network: A Large-Scale Meta-Analysis. Sci. Rep. 2020, 10, 16096. [Google Scholar] [CrossRef]
  134. Zotev, V.; Phillips, R.; Yuan, H.; Misaki, M.; Bodurka, J. Self-Regulation of human brain activity using simultaneous real-time fMRI and EEG neurofeedback. NeuroImage 2014, 85, 985–995. [Google Scholar] [CrossRef] [Green Version]
  135. Wong, C.-K.; Zotev, V.; Misaki, M.; Phillips, R.; Luo, Q.; Bodurka, J. Automatic EEG-assisted retrospective motion correction for fMRI (aE-REMCOR). Neuroimage 2016, 129, 133. [Google Scholar] [CrossRef] [Green Version]
  136. Birn, R.M.; Diamond, J.B.; Smith, M.A.; Bandettini, P.A. Separating respiratory-variation-related fluctuations from neuronal-activity-related fluctuations in fMRI. Neuroimage 2006, 31, 1536–1548. [Google Scholar] [CrossRef]
  137. Chang, C.; Cunningham, J.P.; Glover, G.H. Influence of heart rate on the BOLD signal: The cardiac response function. Neuroimage 2009, 44, 857–869. [Google Scholar] [CrossRef] [Green Version]
  138. Peirce, J.; Gray, J.R.; Simpson, S.; MacAskill, M.; Höchenberger, R.; Sogo, H.; Kastman, E.; Lindeløv, J.K. PsychoPy2: Experiments in behavior made easy. Behav. Res. Methods 2019, 51, 195–203. [Google Scholar] [PubMed] [Green Version]
  139. Zotev, V.; Krueger, F.; Phillips, R.; Alvarez, R.P.; Simmons, W.K.; Bellgowan, P.; Drevets, W.C.; Bodurka, J. Self-Regulation of amygdala activation using real-time fMRI neurofeedback. PLoS ONE 2011, 6, e24522. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Real-time fMRI neurofeedback-augmented mindfulness training task [49]. (a) Posterior cingulate cortex (PCC, MNI coordinates: x = −5, y = −55, z = 23) was selected as the target (ROI, spheres of 7 mm radius) for the real-time fMRI neurofeedback (rtfMRI-nf) training. (b) The experimental protocol consisted of eight fMRI runs, including an anatomical scan, Resting State scan 1 (Rest-1), Observe (OBS), three neurofeedback runs (NF-1, NF-2, NF-3), a Transfer run (TRS), and Resting State scan 2 (Rest-2). During Rest runs (lasting 6 min), the participants were instructed to clear their minds and not to think about anything in particular while fixating at the display screen. OBS, NF-1, NF-2, NF-3, and TRS runs each lasted 6 min and 56 s. They started with a 66 s rest block, followed by alternating Focus-on-Breath (mindfulness training condition; 70 s), Describe (Active baseline condition; 20 s), and Rest (Baseline condition; 30 s) blocks. During the Focus condition, participants were instructed to pay attention to the physical sensations of their breath, not trying to change it in any way, and if their attention were to wander to something else, to gently bring it back to their breath. In the Describe condition, participants were presented with various adjectives, which they had to mentally categorize as descriptive or not descriptive of them. During the Rest condition, the participants were presented with the cue “Rest” and asked to relax while looking at the display screen. No neurofeedback was provided (no bars displayed) during the Rest and Describe conditions or during the entire OBS and TR runs. (c) During the Focus condition, participants viewed a graphical user interface (GUI) screen with neurofeedback bars (blue) and target bars (green). The participants were told that the blue bar may change with their experience of focusing on the breath, and that their goal was to make the blue bar match the green bar as often as possible. The target levels were −0.5%, −0.75%, and −1.0% (% signal change relative to the preceding rest block) for NF-1, NF-2, and NF-3, respectively.
Figure 1. Real-time fMRI neurofeedback-augmented mindfulness training task [49]. (a) Posterior cingulate cortex (PCC, MNI coordinates: x = −5, y = −55, z = 23) was selected as the target (ROI, spheres of 7 mm radius) for the real-time fMRI neurofeedback (rtfMRI-nf) training. (b) The experimental protocol consisted of eight fMRI runs, including an anatomical scan, Resting State scan 1 (Rest-1), Observe (OBS), three neurofeedback runs (NF-1, NF-2, NF-3), a Transfer run (TRS), and Resting State scan 2 (Rest-2). During Rest runs (lasting 6 min), the participants were instructed to clear their minds and not to think about anything in particular while fixating at the display screen. OBS, NF-1, NF-2, NF-3, and TRS runs each lasted 6 min and 56 s. They started with a 66 s rest block, followed by alternating Focus-on-Breath (mindfulness training condition; 70 s), Describe (Active baseline condition; 20 s), and Rest (Baseline condition; 30 s) blocks. During the Focus condition, participants were instructed to pay attention to the physical sensations of their breath, not trying to change it in any way, and if their attention were to wander to something else, to gently bring it back to their breath. In the Describe condition, participants were presented with various adjectives, which they had to mentally categorize as descriptive or not descriptive of them. During the Rest condition, the participants were presented with the cue “Rest” and asked to relax while looking at the display screen. No neurofeedback was provided (no bars displayed) during the Rest and Describe conditions or during the entire OBS and TR runs. (c) During the Focus condition, participants viewed a graphical user interface (GUI) screen with neurofeedback bars (blue) and target bars (green). The participants were told that the blue bar may change with their experience of focusing on the breath, and that their goal was to make the blue bar match the green bar as often as possible. The target levels were −0.5%, −0.75%, and −1.0% (% signal change relative to the preceding rest block) for NF-1, NF-2, and NF-3, respectively.
Brainsci 12 00363 g001
Figure 2. (a) Three insular cortex regions of interest (ROIs) extracted from the Brainnetome Atlas. (bd) Parameter estimate (Focus-on-Breath vs. Describe) across each run in each insula subregion. See Appendix B Table A2 for peak coordinates for each insula subregion. Abbreviations: OBS, Observe; NF, Neurofeedback; TRS, Transfer.
Figure 2. (a) Three insular cortex regions of interest (ROIs) extracted from the Brainnetome Atlas. (bd) Parameter estimate (Focus-on-Breath vs. Describe) across each run in each insula subregion. See Appendix B Table A2 for peak coordinates for each insula subregion. Abbreviations: OBS, Observe; NF, Neurofeedback; TRS, Transfer.
Brainsci 12 00363 g002
Table 1. Sample demographic information.
Table 1. Sample demographic information.
Demographic%
Race
 White71
 Black2
 Asian5
 American Indian/Alaska Native11
 Biracial/Multiracial11
Education
 7th grade18
 8th grade38
 9th grade18
 10th grade16
 11th grade7
 13th grade 3
Family Income
$0–$49,9997
$50,000–$99,99938
$100,000–$149,99921
$150,000–$199,99917
 >$200,000 17
Table 2. Unadjusted means, standard deviations, effect sizes, and main analyses of task ratings and symptom measures.
Table 2. Unadjusted means, standard deviations, effect sizes, and main analyses of task ratings and symptom measures.
Task RatingsMeanSDEstimateSEtpCohen’s d
Focus-on-Breath
 OBS7.091.54
 NF-16.561.67−0.510.36−1.420.16−0.25
 NF-26.451.70−0.630.36−1.720.09−0.30
 NF-36.702.04−0.370.35−1.050.30−0.18
 TRS6.531.61−0.500.35−1.420.16−0.25
Mind Wander
 OBS5.232.12
 NF-14.851.86−0.470.36−1.310.19−0.23
 NF-25.212.09−0.110.36−0.310.76−0.05
 NF-35.831.830.100.350.290.770.05
 TRS5.751.840.450.361.260.210.22
Current Feeling
 OBS2.891.53
 NF-13.411.760.380.291.300.200.22
 NF-23.061.690.150.290.510.610.09
 NF-33.031.780.050.290.160.870.03
 TRS3.062.040.100.290.340.740.06
MeasureMeanSDEstimateSEtpCohen’s d
State Mindfulness Scale (SMS)
 T171.2214.09
 T274.6812.693.461.432.41<0.050.53
Note. Task ratings were answered following the completion of each run. Questions included “ How easy did you find it to focus on your breath?” (1 = not easy at all; 10 = very easy); “How much did your mind wander while you were asked to focus on your breath?” (1 = not at all; 10 = all the time); “How do you feel right now?” (1 = perfectly calm; 10 = very anxious). Abbreviations: OBS, Observe; NF, Neurofeedback; TR, Transfer; T1, Pre-NAMT; T2, Post-NAMT.
Table 3. Uncorrected correlations between parameter estimate (Focus-on-Breath vs. Describe) in PCC and insula subregions across runs.
Table 3. Uncorrected correlations between parameter estimate (Focus-on-Breath vs. Describe) in PCC and insula subregions across runs.
INS SubregionsRunPCC_OBSPCC_NF-1PCC_NF-2PCC_NF-3PCC_TRS
aINSOBS0.65 ***----
NF-1-0.12---
NF-2--0.35 *--
NF-3---0.37 *-
TRS----0.34 *
mINSOBS0.69 ***----
NF-1-0.09---
NF-2--0.36 *--
NF-3---0.42 **-
TRS----0.35 *
pINSOBS0.71 ***----
NF-1-0.45 **---
NF-2--0.56 ***--
NF-3---0.53 ***-
TRS----0.42 **
Abbreviations: NF, Neurofeedback; OBS, Observe; TR, Transfer. * p < 0.05. ** p < 0.01. *** p < 0.001.
Table 4. Unadjusted means, standard deviations, effect sizes, and main analyses of fMRI coefficient for the Focus-on-Breath vs. Describe contrast in subregions of insula across runs.
Table 4. Unadjusted means, standard deviations, effect sizes, and main analyses of fMRI coefficient for the Focus-on-Breath vs. Describe contrast in subregions of insula across runs.
RunMeanSDEstimateSEtpCohen’s d
Anterior insular cortex (aINS)
 OBS−0.160.19
 NF-1−0.010.180.160.043.81<0.0010.63
 NF-2−0.020.210.140.043.26<0.010.54
 NF-3−0.050.280.110.042.75<0.010.46
 TRS−0.120.240.050.041.16<0.0010.19
Mid-insular cortex (mINS)
 OBS00.18
 NF-1−0.060.17−0.060.03−1.700.09−0.28
 NF-2−0.050.15−0.050.03−1.470.15−0.24
 NF-3−0.080.2−0.080.03−2.420.05−0.40
 TRS−0.060.19−0.060.03−1.630.11−0.27
Posterior insular cortex (pINS)
 OBS0.010.17
 NF-1−0.170.14−0.180.03−5.25<0.001−0.88
 NF-2−0.160.15−0.170.03−4.91<0.001−0.82
 NF-3−0.170.18−0.180.03−5.27<0.001−0.88
 TRS−0.050.17−0.060.03−1.720.09−0.29
Abbreviations: OBS, Observe; NF, Neurofeedback; TRS, Transfer.
Table 5. Post hoc comparisons for parameter estimates (Focus-on-Breath vs. Describe) in INS subregions across runs.
Table 5. Post hoc comparisons for parameter estimates (Focus-on-Breath vs. Describe) in INS subregions across runs.
RunEstimateStd. Errorz Statisticp Value
Anterior insular cortex (aINS)
NF-1:OBS0.160.043.81<0.01
NF-2:OBS0.140.043.26<0.05
NF-3:OBS0.110.042.75<0.05
TR:OBS0.050.041.160.78
NF-2:NF-1−0.020.04−0.520.98
NF-3:NF-1−0.040.04−1.060.83
TR:NF-1−0.110.04−2.660.06
NF-3:NF-2−0.020.04−0.530.98
TR:NF-2−0.090.04−2.110.22
TR:NF-3−0.070.04−1.600.50
Mid-insular cortex (mINS)
NF-1:OBS−0.060.03−1.700.43
NF-2:OBS−0.050.03−1.470.59
NF-3:OBS−0.080.03−2.420.11
TR:OBS−0.060.03−1.630.48
NF-2:NF-10.010.030.221.00
NF-3:NF-1−0.020.03−0.720.95
TR:NF-10.000.030.071.00
NF-3:NF-2−0.030.03−0.930.89
TR:NF-2−0.010.03−0.151.00
TR:NF-30.030.030.790.93
Posterior insular cortex (pINS)
NF-1:OBS−0.180.03−5.25<0.001
NF-2:OBS−0.170.03−4.91<0.001
NF-3:OBS−0.180.03−5.27<0.001
TR:OBS−0.060.03−1.720.42
NF-2:NF-10.010.030.301.00
NF-3:NF-10.000.03−0.021.00
TR:NF-10.120.033.52<0.01
NF-3:NF-2−0.010.03−0.321.00
TR:NF-20.110.033.20<0.05
TR:NF-30.120.033.55<0.01
Abbreviations: NF, Neurofeedback; OBS, Observe; TR, Transfer.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yu, X.; Cohen, Z.P.; Tsuchiyagaito, A.; Cochran, G.; Aupperle, R.L.; Stewart, J.L.; Singh, M.K.; Misaki, M.; Bodurka, J.; Paulus, M.P.; et al. Neurofeedback-Augmented Mindfulness Training Elicits Distinct Responses in the Subregions of the Insular Cortex in Healthy Adolescents. Brain Sci. 2022, 12, 363. https://doi.org/10.3390/brainsci12030363

AMA Style

Yu X, Cohen ZP, Tsuchiyagaito A, Cochran G, Aupperle RL, Stewart JL, Singh MK, Misaki M, Bodurka J, Paulus MP, et al. Neurofeedback-Augmented Mindfulness Training Elicits Distinct Responses in the Subregions of the Insular Cortex in Healthy Adolescents. Brain Sciences. 2022; 12(3):363. https://doi.org/10.3390/brainsci12030363

Chicago/Turabian Style

Yu, Xiaoqian, Zsofia P. Cohen, Aki Tsuchiyagaito, Gabriella Cochran, Robin L. Aupperle, Jennifer L. Stewart, Manpreet K. Singh, Masaya Misaki, Jerzy Bodurka, Martin P. Paulus, and et al. 2022. "Neurofeedback-Augmented Mindfulness Training Elicits Distinct Responses in the Subregions of the Insular Cortex in Healthy Adolescents" Brain Sciences 12, no. 3: 363. https://doi.org/10.3390/brainsci12030363

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop