Mindfulness Meditation Is Related to Long-Lasting Changes in Hippocampal Functional Topology during Resting State: A Magnetoencephalography Study

Abstract

It has been suggested that the practice of meditation is associated to neuroplasticity phenomena, reducing age-related brain degeneration and improving cognitive functions. Neuroimaging studies have shown that the brain connectivity changes in meditators. In the present work, we aim to describe the possible long-term effects of meditation on the brain networks. To this aim, we used magnetoencephalography to study functional resting-state brain networks in Vipassana meditators. We observed topological modifications in the brain network in meditators compared to controls. More specifically, in the theta band, the meditators showed statistically significant ($p$ corrected = 0.009) higher degree (a centrality index that represents the number of connections incident upon a given node) in the right hippocampus as compared to controls. Taking into account the role of the hippocampus in memory processes, and in the pathophysiology of Alzheimer’s disease, meditation might have a potential role in a panel of preventive strategies.

This work is licensed under http://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.

1. Introduction

Several studies have shown that the constant practice of meditation induces neuroplasticity phenomena, including the reduction of age-related brain degeneration [1–3] and the improvement of cognitive functions [4]. More specifically, the effects of meditation are correlated to improvements in attention [5], working memory [6], spatial abilities [7], and long-term memory [2].

There are several meditation traditions that share some similarities but differ in practices and theoretical orientation [8]. Recently, meditation practices have been divided into two broad categories: focused attention (or concentrative style) meditation and open awareness (or mindfulness) meditation, depending on how the attentional processes are directed [9]. Concentrative meditation techniques consist in focusing on specific mental or sensory activity, such as a repeating sound, a mental image, or specific bodily sensations such as breathing. On the other hand, the fundamental feature of mindfulness meditation is a particular type of attention characterized by a nonjudgmental thinking which allows the meditator to act in the daily life in a “nonreactive” manner [2, 5].

Over the past ten years, scientific studies have begun to focus on one of the main open awareness or mindfulness meditations, named Vipassana meditation (VM) [10]. VM is a Buddhist practice that consists in focusing on the sensory awareness of the moment with a mental condition of calm and nonreactivity [11, 12]. In the ancient Indian language Pali, Vipassana means “introspection, penetrative vision, observation and understanding of reality as it is” [13]. Through the VM, one may reach meta-awareness [14], having cognition of the mental and emotional processes occurring during meditation [15].

Neuroimaging and electroencephalographic studies have shown that the brain connectivity of meditators changes as they meditate [9, 15, 16] as well as in the resting state [17–19]. Furthermore, the constant practice of meditation causes widespread long-term changes in structural connectivity [20], suggesting that meditation might induce neural plasticity.

Magnetoencephalography (MEG) has been used to study the brain activity since it directly captures the neuronal oscillations [21]. While retaining the high temporal resolution of the EEG, MEG signals are not distorted by the layers surrounding the brain, allowing for a temporally and spatially precise reconstruction of the neural activity within the brain [22]. Using MEG, Dor-Ziderman et al. have documented a correlation between the distinct type of self-awareness, reached through mindfulness meditation, and the power decrease in gamma and beta bands in several specific brain areas [23, 24]. Furthermore, MEG (as well as EEG) signals allow the use of the phase in order to quantify synchronization between brain areas [25].

In recent years, a graph theoretical approach has been applied to functional signals in order to extract relevant features of the interactions among brain areas [26]. Applying network theory to MEG signals, Marzetti et al. have shown that meditation induces a coupling of the posterior cingulate cortex to nodes of the Default Mode Network and of the executive control brain network in the alpha frequency band [27]. More recently, in an EEG study, it has been shown that meditation is associated with increased brain network integration [28].

Despite the well-documented role of the meditative practice on brain functions, little attention has been granted to the long-lasting effects of meditative practice on synchronization between neuronal ensembles and on brain topology. To the best of our knowledge, there are no source-level MEG studies that have investigated the brain network synchronization and its topological features in relation to meditative practice.

In this work, we hypothesized that the routine practice of mindfulness meditation may induce long-lasting effects in the topology of brain areas involved in attention and/or memory. To test our hypothesis, we have considered the sources within the brain as nodes and the functional interactions between them as links. We based the link estimation on the Phase Lag Index (PLI), which quantifies phase synchronization between cerebral areas [29]. Subsequently, we obtained the minimum spanning tree (MST), which allows the construction of a unique subgraph that connects all the anatomical nodes, identifying the backbone of the network [30]. We applied such analysis to a cohort of routinely Vipassana meditators and to controls (people with no experience in meditation) in the resting-state condition, comparing their brain network topology.

2. Methods

2.1. Participants

Twenty-nine people who had been practicing VM for more than one year and 31 people who had never meditated in their life were recruited for the study. The meditators were from a Vipassana community located in Caivano (Naples). All participants had been meditating in the last year for 1 h or more, for at least five days a week.

Criteria of inclusion were the following: (1) no major internal illnesses, (2) no neurological or psychiatric illnesses, (3) left dominance, and (4) participation in Vipassana retreats in the last year. Three meditators were rejected: one left-handed, one presenting a neurological ailment, and one with insufficient quality of the MEG signals. Two controls were discarded due to excessive signal noise.

The final Vipassana meditator group was composed of twenty-six participants (eight males and eighteen females), mean age 42.6 years (±SE 2.4), mean years of meditation experience 6.41 (±SE 1.4). The final control group was composed of twenty-nine (nine males and twenty females) age-gender-education-race-matched controls, mean age 43 (±SE 1.98), no experience in meditation. All participants were native Italian speakers (Table 1).

Table 1

Cohort characterization. Means and standard error are reported.

Features	Size	Gender(M/F)	Age	Education	Years of meditation	Hours of meditation/day
Groups
Meditators	26	8/18	42.30 ± 2.4	15.12 ± 0.66	6.4 ± 1.4	1.33 ± 0.72
Controls	29	9/20	43 ± 1.98	15.72 ± 0.70	—	—

This study was carried out in accordance with the recommendations of Comitato Etico Campania Centro. The protocol was approved by the “Comitato Etico Campania Centro” (Prot.n.93C.E./Reg. n.14-17OSS). All subjects gave written informed consent in accordance with the Declaration of Helsinki.

2.2. MEG Acquisition

The MEG system was developed by the National Research Council, Pozzuoli, Naples, at the Institute of Applied Sciences and Intelligent Systems “E. Caianiello.” The system is equipped with 163 magnetometers, including 9 reference sensors that are located further away from the helmet [31], and it was placed in a magnetically shielded room (ATB, Ulm, Germany) to reduce background noise. Acquisition, preprocessing, and source reconstruction were carried out as previously described [32]. More in details, before the acquisition, four positions coils were placed on the subject’s head and were digitalized using Fastrak (Polhemus®). The coils were activated, and localized, at the beginning of each segment of registration. The subject was seated on a comfortable armchair placed in the shielded room. Electrocardiographic and electrooculographic signals were corecorded to aid artefact removal [33]. The brain activity (Figure 1(a)) was recorded during resting state for two distinct segments of 2.5 minutes with the eyes closed. MEG data, after an antialiasing filter, were acquired with a sampling frequency of 1024 Hz. The signal was then filtered using a fourth order Butterworth IIR band-pass filter in the 0.5–48 Hz band.

[figure omitted; refer to PDF]

2.3. Preprocessing

Firstly, Principal Component Analysis (PCA) was used to reduce environmental noise [34, 35]. Specifically, the filter is obtained by orthogonalizing the reference signals to obtain a base, projecting the brain sensors on the base of the noise and removing the projections to obtain clean data [35]. We adopted the PCA filtering implementation available within the Fieldtrip toolbox [36]. Subsequently, noisy channels were removed manually through visual inspection by an experienced rater [33]. Finally, for each subject, supervised Independent Component Analysis (ICA) [37] was used to remove physiological, cardiac (generally one component), and blinking (if present) artefacts from the MEG signals.

The first ten epochs of 8 seconds for each subject that did not contain artefacts (either system related or physiological) or strong environmental noise were selected. The length of 8 seconds is a trade-off between the need to have enough cleaned epochs, to avoid drowsiness [33] and to obtain a reliable estimate of the connectivity measure [38].

2.4. Source Reconstruction

All the processing related to the beamforming procedure has been done using the Fieldtrip toolbox [36]. Based on a MRI template, the volume conduction model proposed by Nolte [39] was considered and the Linearly Constrained Minimum Variance (LCMV) beamformer [40] was implemented to reconstruct time series related to the centroids of 116 regions-of-interest (ROIs), derived from the Automated Anatomical Labeling (AAL) atlas [41–43] (Figure 1(b)). Although cerebellar areas have been reconstructed previously in some EEG studies [44, 45], we decided not to include them (hence taking into account 90 ROIs) given their lower reliability [46, 47]. For each source, we projected the time series along the dipole direction that explains the most variance by means of singular value decomposition (SVD). Source time series were resampled at 512 Hz (Figure 1(c)).

2.5. Connectivity and Network Analysis

The PLI was used to estimate functional connectivity [29], using BrainWave software (CJS, version 09.152.1.23, available from http://home.kpn.nl/stam7883/brainwave.html). The PLI is based on the distribution of the differences of the instantaneous phases $\left(∆\mathrm{\Phi }\left(t\right)\right)$ (derived from the Hilbert transform of the time series) and for two-time series and is computed as $$\begin{array}{}\text{(1)}& \mathrm{P}\mathrm{L}\mathrm{I}=\left|⟨\mathrm{sign}\left[\sin \left(\Delta \mathrm{\Phi }\left(t_k\right)\right)\right]⟩\right|,\end{array}$$where “< >” indicates the mean value, “sign” stands for the signum function, “|.|” denotes the absolute value and “$t_k$” are the samples. The phase difference is defined in the [$-\pi$, $\pi$] range. This measure is insensitive to volume conduction (at the cost of discarding true zero–lag interactions). PLI values range between 0 and 1, where 1 indicates perfect synchronization and 0 indicates nonsynchronous activity. We obtained a 90 × 90 adjacency matrix for each epoch for each subject, in all the frequency bands [29]. For each epoch, the PLI matrix was computed, and after this step, they were merged by arithmetic average. Hence, by computing the PLI for each couple of brain regions, we obtained a 90 × 90 weighted adjacency matrix (Figure 1(d)) for each epoch and for each subject, in all of the frequency bands: delta (0.5–4.0 Hz), theta (4.0–8.0 Hz), alpha (8.0–13.0 Hz), beta (13.0–30.0 Hz), and gamma (30.0–48.0 Hz).

The weighted adjacency matrix was used to reconstruct a network, where the 90 areas of the AAL atlas are represented as nodes and the PLI values form the weighted edges. A frequency-specific minimum spanning tree was calculated for each epoch (Figure 1(e)). Since we were interested in the strongest connections, for the construction of the MST, the edge weight was defined as 1/PLI. In fact, Kruskal’s algorithm [48] first ranks the links in ascending order and then constructs the network by adding one link at a time, discarding links that would form a loop. The algorithm proceeds until all nodes are connected resulting in a loopless graph with $N$ nodes and $M=N-1$ links.

To avoid some of the biases in traditional network analyses [49], we used the minimum spanning tree (MST) that allows for an unbiased topological interpretation of the results [50, 51]. Based on the MST, we calculated both global (leaf fraction, degree divergence, and tree hierarchy) and nodal parameters (degree, betweenness centrality, and eccentricity) [51, 52]. The leaf fraction is defined as the fraction of nodes with a degree of 1 [52], providing an indication of the integration of the network. A higher leaf fraction implies that the network tends toward a starlike topology, where the nodes are on average closer to each other as compared to a more line-like topology. The degree divergence is a measure of the broadness of the degree distribution, related to the resilience against attacks [51]. Finally, the tree hierarchy is defined as the number of leaf over the maximum betweenness centrality. The idea behind this measure is that an optimal network should achieve efficient communication while avoiding hub overload. The tree hierarchy has been designed to quantify the balance between both features.

Furthermore, we calculated the degree, the betweenness centrality (BC), and the eccentricity for each node [53], in order to determine if specific regions differed according to the number of years of meditation. The degree is the number of connection incident on a given node. The BC is defined as the number of the shortest paths passing through a given node over the total of the shortest paths of the network [53]. The eccentricity is defined as the longest path between a node and any other node of the network. The lower the eccentricity, the more central the node is [51]. These metrics were calculated for each epoch and subsequently averaged across epochs for each subject separately.

2.6. Statistical Analysis

All statistical analyses were performed in Matlab (Mathworks®, version R2013a). To compare the two groups, for each frequency band and for each parameter, we used permutation testing [54] where the null distribution for between-group differences is derived from the data. In details, assuming no group differences, the labels of the subjects were permuted 10,000 times. Each time, the difference between the averages of the two groups was computed, obtaining the null distribution for between-group differences. Such distribution was used to define the statistical significance of the observed difference between the two groups. We used the False Discovery Rate (FDR) [55] to correct for multiple comparisons. A significance level of $p\le 0.05$ was used.

3. Results

After building a PLI-based adjacency matrix, we used an MST-based network approach in order to compare nodal and global topological parameters in meditators and controls.

Global parameters (leaf fraction, degree divergence, and tree hierarchy) were not significantly different between groups before correction for multiple comparisons between metrics.

With regard to nodal metrics, before correction for multiple comparisons, we found differences between meditators and controls in regions that are typically affected by the meditative practice. In particular, in the amygdala, in the gamma band, both the degree and BC were higher in meditators as compared to controls ($p=0.049$, $p=0.034$, respectively). The caudatus, another area known to be affected by VM, was shown to be more connected in the theta band in meditators ($p=0.0049$). In the prefrontal cortex, we found fairly widespread differences in areas such as (among others) the gyrus rectus in the alpha band ($p=0.0201$) and the cingulum in the theta band ($p=0181$). The calcarine cortex appeared less connected in the theta band ($p=0.011$) and more connected in the gamma band ($p=0.038$). For the reader’s interpretation, we report in Table 2 the uncorrected $p$ values. However, it is important to notice that these results are not FDR corrected and should be considered with extreme caution.

Table 2

Uncorrected $p$ values comparing nodal metrics between meditators and controls. “+” indicates that meditators display higher values as compared to controls, “-” signifies the opposite.

Degree

Betweenness centrality

Eccentricity

Delta

Theta

Alpha

Beta

Gamma

Delta

Theta

Alpha

Beta

Gamma

Delta

Theta

Alpha

Beta

Gamma

Frontal_Inf_Orb_L

0.0055-

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

0.0097-

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Frontal_Sup_Orb_L

n.s.

n.s.

n.s.

n.s.

n.s.

0.0415+

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Frontal_Mid_Orb_L

n.s.

n.s.

n.s.

n.s.

0.0263-

n.s.

n.s.

n.s.

n.s.

0.0319-

n.s.

n.s.

n.s.

n.s.

n.s.

Cingulum_Ant_L

n.s.

0.0337-

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

0.0118+

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Cingulum_Mid_L

0.0181-

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Olfactory_L

n.s.

n.s.

n.s.

n.s.

0.0292-

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Frontal_Sup_L

n.s.

n.s.

n.s.

0.0239-

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Frontal_Sup_Medial_L

n.s.

n.s.

0.0044-

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Frontal_Med_Orb_R

n.s.

n.s.

n.s.

n.s.

0.0338+

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Rectus_R

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

0.0201+

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Temporal_Pole_Mid_L

n.s.

n.s.

0.0456+

n.s.

n.s.

n.s.

n.s.

0.0085+

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Heschl_L

n.s.

n.s.

n.s.

0.0102+

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Hippocampus_L

n.s.

n.s.

n.s.

n.s.

0.0161-

n.s.

n.s.

n.s.

n.s.

0.0472-

n.s.

n.s.

n.s.

n.s.

n.s.

Amygdala_L

n.s.

n.s.

n.s.

n.s.

0.0495+

n.s.

n.s.

n.s.

n.s.

0.0347+

n.s.

n.s.

n.s.

n.s.

n.s.

Temporal_Pole_Mid_R

n.s.

n.s.

n.s.

n.s.

n.s.

0.0411+

n.s.

n.s.

n.s.

0.0183-

n.s.

n.s.

n.s.

n.s.

n.s.

Temporal_Sup_R

0.0064+

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Temporal_Inf_R

n.s.

n.s.

n.s.

n.s.

0.0115+

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

ParaHippocampal_R

n.s.

n.s.

n.s.

n.s.

0.0318+

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Hippocampus_R

n.s.

0.001+

n.s.

n.s.

n.s.

n.s.

0.0081+

n.s.

0.0029+

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Insula_R

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

0.0483+

n.s.

n.s.

Occipital_Sup_L

n.s.

n.s.

n.s.

0.0357+

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Occipital_Mid_L

n.s.

n.s.

n.s.

n.s.

n.s.

0.0155-

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Cuneus_L

n.s.

n.s.

n.s.

0.0062+

n.s.

n.s.

n.s.

0.007-

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Calcarine_L

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

0.0367-

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Occipital_Mid_R

n.s.

n.s.

n.s.

0.0289+

n.s.

n.s.

0.0418+

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Occipital_Sup_R

n.s.

n.s.

n.s.

n.s.

n.s.

0.0138-

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Cuneus_R

n.s.

n.s.

n.s.

n.s.

n.s.

0.03800-

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Calcarine_R

0.0115-

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

0.0381+

n.s.

n.s.

n.s.

n.s.

n.s.

SupraMarginal_R

n.s.

n.s.

n.s.

n.s.

n.s.

0.0034+

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Parietal_Inf_R

n.s.

n.s.

n.s.

n.s.

0.004+

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Pallidum_L

0.0434-

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Caudate_L

n.s.

0.0049+

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

Putamen_R

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

0.0221+

n.s.

n.s.

n.s.

n.s.

n.s.

After correction for multiple comparisons between areas, higher degree ($p$ corrected = 0.009) was found in meditators, in the theta band, in the right hippocampus (Figure 2). No other nodal significances held FDR in any other frequency bands.

[figure omitted; refer to PDF]

4. Discussion

The present study was addressed to verify whether mindfulness meditation may induce long-lasting effects in topological features in brain networks. To this aim, we applied network analysis to compare source-level MEG data obtained from a cohort of routinely Vipassana meditators to controls, during resting-state condition. No difference in the global metrics was observed. Moreover, before FDR correction, we observed altered topology in meditators in a number of brain areas, including the frontal lobes, the temporal lobes, the occipital lobes, the amygdala, and the caudatus. Importantly, after FDR correction, our results showed that the right hippocampus displays a higher degree in the theta band with respect to controls.

As expected, we did not find difference in the global metrics comparing meditators to controls. Global metrics capture widespread reorganization in the brain activity, and those are found to be altered typically in processes that cause diffuse damage to the brain, such as neurological diseases [32, 56–58]. Our results showed that the connectedness or the brain network (Leaf fraction), its scale-freeness (Degree divergence), and its resiliency versus targeted attacks (Tree hierarchy) are not affected by VM. Hence, we conclude that meditation does not induce a reorganization in the whole network.

When considering nodal results not corrected for multiple comparisons, according to the previous literature, we observed areas with altered topology in the frontal lobes [18, 59], in the temporal lobes [16, 60], in the amygdala [61], in the caudatus [19], and in the occipital lobes [62] in multiple frequency bands. The results we found are in multiple frequency bands, which is in line with the previous literature taken altogether [9, 63–65]. However, provided the lack of correction for multiple comparisons, such results are only explorative.

After correction for multiple comparisons between areas, the right hippocampus of meditators showed a higher degree in the theta band. This result is in line with previous evidence showing the central role of the hippocampus in meditative practices [10, 66, 67]. The higher degree implies that meditators have more link incident upon the right hippocampus as compared to controls. This can be interpreted as the right hippocampus being more connected to the rest of the network. More precisely, this is to say that the topology of the area (i.e., the properties of its logical relationships with the rest of the brain) is modified so that the role of that area is more prominent within the network. Provided we are exploring the long-term effects of meditation, one would expect that the areas that are involved more (and more often) in the meditative practice would be the ones undergoing long-lasting functional modifications. According to this evidence, it has been observed that meditators show a larger volume [67] and greater grey matter concentration [68] in the right hippocampus as compared to nonmeditators. Interestingly, the connectivity of the same area [60] is increased during the practice of meditation as well. Such results might be in line with our finding, showing that long-term changes are induced by VM in the right hippocampus.

Our results were specific to the right hippocampus. This might be related to the different functional specialization between the right and left hippocampus. It was observed that while the left hippocampus was activated both in the construction of past events and in the future ones, the right hippocampus was activated only during the creation of future events (prospective memory) [69]. It has been reported that, when new information is being processed, the synchronization between the hippocampus and the prefrontal cortex happens mainly in the theta band [70, 71]. More in detail, during the imagination of future events, novel memory and new elements are learned and prospective memory is built. Besides the involvement in prospective memory, the right hippocampus is also implicated in visuo-spatial memory processes [72–74], with a critical role in the consolidation processes.

Interestingly, we observed that the difference in the centrality of the hippocampus was specific to the theta band, which, as mentioned above, is involved in memory processes [75]. According to Buzsáki, theta oscillations are considered as an “essential temporal organizer,” helping the orientation in time during episodic memory as well as in the physical world [76]. Increases in the theta power during meditative practice have been widely reported [77, 78] also in the temporal areas [79].

In the light of the functional role of the right hippocampus and its activation in the theta band, our results indicate that meditation might have functional effects on prospective and spatial memories. These findings suggest that meditation, if constantly practiced, could be a candidate as a nonpharmacological intervention in pathologies characterized by alteration in the hippocampal areas.

It should finally be noted that in the literature there are often conflicting data. This may be explained by the different methodological approaches, such as the analysis of the meditative phase vs. resting state, the type of meditation explored (focused attention vs. mindfulness meditation), and finally, the modality of data acquisition and/or data analysis. Of note, our study is the only study to date which has analysed the difference of the brain network between Vipassana meditators and controls during the resting state and using source-level magnetoencephalography.

Our results seem to bear a promise of what the effects of meditation might be on brain networks. Nevertheless, it is important to underline that they are preliminary and that further analyses will have to be carried out. For example, it is essential to investigate perspective and spatial memory in meditators and nonmeditators during tasks in order to relate the functional connectivity with cognitive functioning. Furthermore, it will be helpful to analyse the brain networks in people that practice other forms of meditation, in order to better characterize the benefits of specific meditative practices.

Despite the aforementioned limits, we showed that routine meditative practice is associated with a long-lasting change in the topology of definite brain areas, suggesting that meditation might be able to induce brain plasticity.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

The present paper was supported by the University of Naples Parthenope Ricerca Competitiva 2017 (D.R. 289/2017) to LM and GF and by grants from the Regione Campania (Sportello Innovazione) to GS. We thank the “Prana Vidya Yoga” Association.

References

[1] E. Luders, "Exploring age-related brain degeneration in meditation practitioners," Annals of the New York Academy of Sciences, vol. 1307 no. 1, pp. 82-88, DOI: 10.1111/nyas.12217, 2014.

[2] A. B. Newberg, M. Serruya, N. Wintering, A. S. Moss, D. Reibel, D. A. Monti, "Meditation and neurodegenerative diseases," Annals of the New York Academy of Sciences, vol. 1307 no. 1, pp. 112-123, DOI: 10.1111/nyas.12187, 2014.

[3] D. S. Khalsa, "Stress, meditation, and Alzheimer’s disease prevention: where the evidence stands," Journal of Alzheimer's Disease, vol. 48 no. 1,DOI: 10.3233/JAD-142766, 2015.

[4] A. Chiesa, R. Calati, A. Serretti, "Does mindfulness training improve cognitive abilities? A systematic review of neuropsychological findings," Clinical Psychology Review, vol. 31 no. 3, pp. 449-464, DOI: 10.1016/j.cpr.2010.11.003, 2011.

[5] A. Lutz, H. A. Slagter, J. D. Dunne, R. J. Davidson, "Attention regulation and monitoring in meditation," Trends in Cognitive Sciences, vol. 12 no. 4, pp. 163-169, DOI: 10.1016/j.tics.2008.01.005, 2008.

[6] A. P. Jha, E. A. Stanley, A. Kiyonaga, L. Wong, L. Gelfand, "Examining the protective effects of mindfulness training on working memory capacity and affective experience," Emotion, vol. 10 no. 1, pp. 54-64, DOI: 10.1037/a0018438, 2010.

[7] L. Geng, L. Zhang, D. Zhang, "Improving spatial abilities through mindfulness: effects on the mental rotation task," Consciousness and Cognition, vol. 20 no. 3, pp. 801-806, DOI: 10.1016/j.concog.2011.02.004, 2011.

[8] D. Goleman, J. R. Horne, "The varieties of the meditative experience," Religious Studies, vol. 16 no. 3, pp. 365-366, 1980.

[9] B. R. Cahn, J. Polich, "Meditation states and traits: EEG, ERP, and neuroimaging studies," Psychological Bulletin, vol. 132 no. 2, pp. 180-211, DOI: 10.1037/0033-2909.132.2.180, 2006.

[10] A. Chiesa, "Vipassana meditation: systematic review of current evidence," The Journal of Alternative and Complementary Medicine, vol. 16 no. 1, pp. 37-46, DOI: 10.1089/acm.2009.0362, 2010.

[11] V. H. M. Gunaratana, Mindfulness in Plain English. 1990, 2002.

[12] W. Hart, The Art of Living: Vipassana Meditation as Taught by S. N, 1987.

[13] W. Hart, The Art of Living: Vipassana Meditation as Taught by SN Goenka, 2011.

[14] J. H. Flavell, "Metacognition and cognitive monitoring: a new area of cognitive-developmental inquiry," American Psychologist, vol. 34 no. 10, pp. 906-911, DOI: 10.1037/0003-066X.34.10.906, 1979.

[15] H. Wahbeh, A. Sagher, W. Back, P. Pundhir, F. Travis, "A systematic review of transcendent states across meditation and contemplative traditions," EXPLORE, vol. 14 no. 1, pp. 19-35, DOI: 10.1016/j.explore.2017.07.007, 2018.

[16] K. B. Baerentsen, N. V. Hartvig, H. Stødkilde-Jørgensen, J. Mammen, "Onset of meditation explored with fMRI," NeuroImage, vol. 13 no. 6,DOI: 10.1016/S1053-8119(01)91640-4, 2001.

[17] J. H. Jang, W. H. Jung, D.-H. Kang, M. S. Byun, S. J. Kwon, C. H. Choi, J. S. Kwon, "Increased default mode network connectivity associated with meditation," Neuroscience Letters, vol. 487 no. 3, pp. 358-362, DOI: 10.1016/j.neulet.2010.10.056, 2011.

[18] V. A. Taylor, V. Daneault, J. Grant, G. Scavone, E. Breton, S. Roffe-Vidal, J. Courtemanche, A. S. Lavarenne, G. Marrelec, H. Benali, M. Beauregard, "Impact of meditation training on the default mode network during a restful state," Social Cognitive and Affective Neuroscience, vol. 8 no. 1,DOI: 10.1093/scan/nsr087, 2013.

[19] T. Gard, M. Taquet, R. Dixit, B. K. Hölzel, B. C. Dickerson, S. W. Lazar, "Greater widespread functional connectivity of the caudate in older adults who practice kripalu yoga and vipassana meditation than in controls," Frontiers in Human Neuroscience, vol. 9,DOI: 10.3389/fnhum.2015.00137, 2015.

[20] E. Luders, K. Clark, K. L. Narr, A. W. Toga, "Enhanced brain connectivity in long-term meditation practitioners," NeuroImage, vol. 57 no. 4, pp. 1308-1316, DOI: 10.1016/j.neuroimage.2011.05.075, 2011.

[21] F. Lopes da Silva, "EEG and MEG: relevance to neuroscience," Neuron, vol. 80 no. 5, pp. 1112-1128, DOI: 10.1016/j.neuron.2013.10.017, 2013.

[22] S. Baillet, "Magnetoencephalography for brain electrophysiology and imaging," Nature Neuroscience, vol. 20 no. 3, pp. 327-339, DOI: 10.1038/nn.4504, 2017.

[23] Y. Dor-Ziderman, A. Berkovich-Ohana, J. Glicksohn, A. Goldstein, "Mindfulness-induced selflessness: a MEG neurophenomenological study," Frontiers in Human Neuroscience, vol. 7,DOI: 10.3389/fnhum.2013.00582, 2013.

[24] P. Sorrentino, D. Nieboer, J. W. R. Twisk, C. J. Stam, L. Douw, A. Hillebrand, "The hierarchy of brain networks is related to insulin growth factor-1 in a large, middle-aged, healthy cohort: an exploratory magnetoencephalography study," Brain Connectivity, vol. 7 no. 5, pp. 321-330, DOI: 10.1089/brain.2016.0469, 2017.

[25] P. Tass, M. G. Rosenblum, J. Weule, J. Kurths, A. Pikovsky, J. Volkmann, A. Schnitzler, H.-J. Freund, "Detection of n : m Phase Locking from Noisy Data: Application to Magnetoencephalography," Physical Review Letters, vol. 81, pp. 3291-3294, DOI: 10.1103/physrevlett.81.3291, 1998.

[26] E. Bullmore, O. Sporns, "Complex brain networks: graph theoretical analysis of structural and functional systems," Nature Reviews Neuroscience, vol. 10 no. 3, pp. 186-198, DOI: 10.1038/nrn2575, 2009.

[27] L. Marzetti, C. Di Lanzo, F. Zappasodi, F. Chella, A. Raffone, V. Pizzella, "Magnetoencephalographic alpha band connectivity reveals differential default mode network interactions during focused attention and open monitoring meditation," Frontiers in Human Neuroscience, vol. 8,DOI: 10.3389/fnhum.2014.00832, 2014.

[28] R. van Lutterveld, E. van Dellen, P. Pal, H. Yang, C. J. Stam, J. Brewer, "Meditation is associated with increased brain network integration," NeuroImage, vol. 158, pp. 18-25, DOI: 10.1016/j.neuroimage.2017.06.071, 2017.

[29] C. J. Stam, G. Nolte, A. Daffertshofer, "Phase lag index: assessment of functional connectivity from multi channel EEG and MEG with diminished bias from common sources," Human Brain Mapping, vol. 28 no. 11, pp. 1178-1193, DOI: 10.1002/hbm.20346, 2007.

[30] C. J. Stam, "Modern network science of neurological disorders," Nature Reviews Neuroscience, vol. 15 no. 10, pp. 683-695, DOI: 10.1038/nrn3801, 2014.

[31] S. Rombetto, C. Granata, A. Vettoliere, M. Russo, "Multichannel system based on a high sensitivity superconductive sensor for magnetoencephalography," Sensors, vol. 14 no. 7, pp. 12114-12126, DOI: 10.3390/s140712114, 2014.

[32] P. Sorrentino, R. Rucco, F. Jacini, F. Trojsi, A. Lardone, F. Baselice, C. Femiano, G. Santangelo, C. Granata, A. Vettoliere, M. R. Monsurrò, G. Tedeschi, G. Sorrentino, "Brain functional networks become more connected as amyotrophic lateral sclerosis progresses: a source level magnetoencephalographic study," NeuroImage: Clinical, vol. 20, pp. 564-571, DOI: 10.1016/j.nicl.2018.08.001, 2018.

[33] J. Gross, S. Baillet, G. R. Barnes, R. N. Henson, A. Hillebrand, O. Jensen, K. Jerbi, V. Litvak, B. Maess, R. Oostenveld, L. Parkkonen, J. R. Taylor, V. van Wassenhove, M. Wibral, J. M. Schoffelen, "Good practice for conducting and reporting MEG research," NeuroImage, vol. 65, pp. 349-363, DOI: 10.1016/j.neuroimage.2012.10.001, 2013.

[34] P. K. Sadasivan, D. N. Dutt, "SVD based technique for noise reduction in electroencephalographic signals," Signal Processing, vol. 55 no. 2, pp. 179-189, DOI: 10.1016/S0165-1684(96)00129-6, 1996.

[35] A. de Cheveigné, J. Z. Simon, "Denoising based on time-shift PCA," Journal of Neuroscience Methods, vol. 165 no. 2, pp. 297-305, DOI: 10.1016/j.jneumeth.2007.06.003, 2007.

[36] R. Oostenveld, P. Fries, E. Maris, J.-M. Schoffelen, "FieldTrip: open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data," Computational Intelligence and Neuroscience, vol. 2011,DOI: 10.1155/2011/156869, 2011.

[37] G. Barbati, C. Porcaro, F. Zappasodi, P. M. Rossini, F. Tecchio, "Optimization of an independent component analysis approach for artifact identification and removal in magnetoencephalographic signals," Clinical Neurophysiology, vol. 115 no. 5, pp. 1220-1232, DOI: 10.1016/j.clinph.2003.12.015, 2004.

[38] M. Fraschini, M. Demuru, A. Crobe, F. Marrosu, C. J. Stam, A. Hillebrand, "The effect of epoch length on estimated EEG functional connectivity and brain network organisation," Journal of Neural Engineering, vol. 13 no. 3, article 036015,DOI: 10.1088/1741-2560/13/3/036015, 2016.

[39] G. Nolte, "The magnetic lead field theorem in the quasi-static approximation and its use for magnetoencephalography forward calculation in realistic volume conductors," Physics in Medicine and Biology, vol. 48 no. 22, pp. 3637-3652, DOI: 10.1088/0031-9155/48/22/002, 2003.

[40] B. D. Van Veen, W. Van Drongelen, M. Yuchtman, A. Suzuki, "Localization of brain electrical activity via linearly constrained minimum variance spatial filtering," IEEE Transactions on Biomedical Engineering, vol. 44 no. 9, pp. 867-880, DOI: 10.1109/10.623056, 1997.

[41] G. Gong, Y. He, L. Concha, C. Lebel, D. W. Gross, A. C. Evans, C. Beaulieu, "Mapping anatomical connectivity patterns of human cerebral cortex using in vivo diffusion tensor imaging tractography," Cerebral Cortex, vol. 19 no. 3, pp. 524-536, DOI: 10.1093/cercor/bhn102, 2009.

[42] A. Hillebrand, P. Tewarie, E. van Dellen, M. Yu, E. W. S. Carbo, L. Douw, A. A. Gouw, E. C. W. van Straaten, C. J. Stam, "Direction of information flow in large-scale resting-state networks is frequency-dependent," Proceedings of the National Academy of Sciences of the United States of America, vol. 113 no. 14, pp. 3867-3872, DOI: 10.1073/pnas.1515657113, 2016.

[43] N. Tzourio-Mazoyer, B. Landeau, D. Papathanassiou, F. Crivello, O. Etard, N. Delcroix, B. Mazoyer, M. Joliot, "Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain," NeuroImage, vol. 15 no. 1, pp. 273-289, DOI: 10.1006/nimg.2001.0978, 2002.

[44] A.-M. Cebolla, E. Palmero-Soler, A. Leroy, G. Cheron, "EEG spectral generators involved in motor imagery: a swLORETA study," Frontiers in Psychology, vol. 8,DOI: 10.3389/fpsyg.2017.02133, 2017.

[45] A. M. Cebolla, M. Petieau, B. Dan, L. Balazs, J. McIntyre, G. Cheron, "Cerebellar contribution to visuo-attentional alpha rhythm: insights from weightlessness," Scientific Reports, vol. 6 no. 1, article 37824,DOI: 10.1038/srep37824, 2016.

[46] K. Jerbi, J. P. Lachaux, K. N'Diaye, D. Pantazis, R. M. Leahy, L. Garnero, S. Baillet, "Coherent neural representation of hand speed in humans revealed by MEG imaging," Proceedings of the National Academy of Sciences of the United States of America, vol. 104 no. 18, pp. 7676-7681, DOI: 10.1073/pnas.0609632104, 2007.

[47] M. J. Brookes, M. Woolrich, H. Luckhoo, D. Price, J. R. Hale, M. C. Stephenson, G. R. Barnes, S. M. Smith, P. G. Morris, "Investigating the electrophysiological basis of resting state networks using magnetoencephalography," Proceedings of the National Academy of Sciences of the United States of America, vol. 108 no. 40, pp. 16783-16788, DOI: 10.1073/pnas.1112685108, 2011.

[48] J. B. Kruskal, "On the shortest spanning subtree of a graph and the traveling salesman problem," Proceedings of the American Mathematical Society, vol. 7 no. 1,DOI: 10.1090/S0002-9939-1956-0078686-7, 1956.

[49] B. C. M. van Wijk, C. J. Stam, A. Daffertshofer, "Comparing brain networks of different size and connectivity density using graph theory," PLoS One, vol. 5 no. 10, article e13701,DOI: 10.1371/journal.pone.0013701, 2010.

[50] C. J. Stam, P. Tewarie, E. Van Dellen, E. C. W. van Straaten, A. Hillebrand, P. Van Mieghem, "The trees and the forest: characterization of complex brain networks with minimum spanning trees," International Journal of Psychophysiology, vol. 92 no. 3, pp. 129-138, DOI: 10.1016/j.ijpsycho.2014.04.001, 2014.

[51] P. Tewarie, E. van Dellen, A. Hillebrand, C. J. Stam, "The minimum spanning tree: an unbiased method for brain network analysis," NeuroImage, vol. 104, pp. 177-188, DOI: 10.1016/j.neuroimage.2014.10.015, 2015.

[52] M. Boersma, D. J. A. Smit, D. I. Boomsma, E. J. C. De Geus, H. A. Delemarre-van de Waal, C. J. Stam, "Growing trees in child brains: graph theoretical analysis of electroencephalography-derived minimum spanning tree in 5- and 7-year-old children reflects brain maturation," Brain Connectivity, vol. 3 no. 1, pp. 50-60, DOI: 10.1089/brain.2012.0106, 2013.

[53] L. C. Freeman, "A set of measures of centrality based on betweenness," Sociometry, vol. 40 no. 1,DOI: 10.2307/3033543, 1977.

[54] T. E. Nichols, A. P. Holmes, "Nonparametric permutation tests for functional neuroimaging: a primer with examples," Human Brain Mapping, vol. 15 no. 1,DOI: 10.1002/hbm.1058, 2002.

[55] Y. Benjamini, Y. Hochberg, "Controlling the false discovery rate: a practical and powerful approach to multiple testing," Journal of the Royal Statistical Society. Series B (Methodological), vol. 57 no. 1, pp. 289-300, 1995.

[56] M. Yu, A. A. Gouw, A. Hillebrand, B. M. Tijms, C. J. Stam, E. C. W. van Straaten, Y. A. L. Pijnenburg, "Different functional connectivity and network topology in behavioral variant of frontotemporal dementia and Alzheimer’s disease: an EEG study," Neurobiology of Aging, vol. 42, pp. 150-162, DOI: 10.1016/j.neurobiolaging.2016.03.018, 2016.

[57] P. Tewarie, A. Hillebrand, M. M. Schoonheim, B. W. van Dijk, J. J. G. Geurts, F. Barkhof, C. H. Polman, C. J. Stam, "Functional brain network analysis using minimum spanning trees in multiple sclerosis: an MEG source-space study," NeuroImage, vol. 88, pp. 308-318, DOI: 10.1016/j.neuroimage.2013.10.022, 2014.

[58] F. Trojsi, P. Sorrentino, G. Sorrentino, G. Tedeschi, "Neurodegeneration of brain networks in the amyotrophic lateral sclerosis–frontotemporal lobar degeneration (ALS–FTLD) continuum: evidence from MRI and MEG studies," CNS Spectrums,DOI: 10.1017/S109285291700075X, 2017.

[59] N. A. S. Farb, Z. V. Segal, H. Mayberg, J. Bean, D. McKeon, Z. Fatima, A. K. Anderson, "Attending to the present: mindfulness meditation reveals distinct neural modes of self-reference," Social Cognitive and Affective Neuroscience, vol. 2 no. 4, pp. 313-322, DOI: 10.1093/scan/nsm030, 2007.

[60] S. W. Lazar, I. S. Rosman, M. Vangel, V. Rao, H. Dusek, H. Benson, R. L. Gollub, Functional Brain Imaging of Mindfulness and Mantra-Based Meditation, 2003.

[61] B. K. Hölzel, S. W. Lazar, T. Gard, Z. Schuman-Olivier, D. R. Vago, U. Ott, "How does mindfulness meditation work? Proposing mechanisms of action from a conceptual and neural perspective," Perspectives on Psychological Science, vol. 6 no. 6, pp. 537-559, DOI: 10.1177/1745691611419671, 2011.

[62] B. R. Cahn, A. Delorme, J. Polich, "Occipital gamma activation during Vipassana meditation," Cognitive Processing, vol. 11 no. 1, pp. 39-56, DOI: 10.1007/s10339-009-0352-1, 2010.

[63] T. Lomas, I. Ivtzan, C. H. Y. Fu, "A systematic review of the neurophysiology of mindfulness on EEG oscillations," Neuroscience & Biobehavioral Reviews, vol. 57, pp. 401-410, DOI: 10.1016/j.neubiorev.2015.09.018, 2015.

[64] J. Fell, N. Axmacher, S. Haupt, "From alpha to gamma: electrophysiological correlates of meditation-related states of consciousness," Medical Hypotheses, vol. 75 no. 2, pp. 218-224, DOI: 10.1016/j.mehy.2010.02.025, 2010.

[65] B. R. Cahn, A. Delorme, J. Polich, "Event-related delta, theta, alpha and gamma correlates to auditory oddball processing during Vipassana meditation," Social Cognitive and Affective Neuroscience, vol. 8 no. 1, pp. 100-111, DOI: 10.1093/scan/nss060, 2013.

[66] R. E. Wells, G. Y. Yeh, C. E. Kerr, J. Wolkin, R. B. Davis, Y. Tan, R. Spaeth, R. B. Wall, J. Walsh, T. J. Kaptchuk, D. Press, R. S. Phillips, J. Kong, "Meditation’s impact on default mode network and hippocampus in mild cognitive impairment: a pilot study," Neuroscience Letters, vol. 556, pp. 15-19, DOI: 10.1016/j.neulet.2013.10.001, 2013.

[67] E. Luders, A. W. Toga, N. Lepore, C. Gaser, "The underlying anatomical correlates of long-term meditation: larger hippocampal and frontal volumes of gray matter," NeuroImage, vol. 45 no. 3, pp. 672-678, DOI: 10.1016/j.neuroimage.2008.12.061, 2009.

[68] B. K. Hölzel, U. Ott, T. Gard, H. Hempel, M. Weygandt, K. Morgen, D. Vaitl, "Investigation of mindfulness meditation practitioners with voxel-based morphometry," Social Cognitive and Affective Neuroscience, vol. 3 no. 1, pp. 55-61, DOI: 10.1093/scan/nsm038, 2008.

[69] D. R. Addis, A. T. Wong, D. L. Schacter, "Remembering the past and imagining the future: common and distinct neural substrates during event construction and elaboration," Neuropsychologia, vol. 45 no. 7, pp. 1363-1377, DOI: 10.1016/j.neuropsychologia.2006.10.016, 2007.

[70] M. W. Jones, M. A. Wilson, "Theta rhythms coordinate hippocampal–prefrontal interactions in a spatial memory task," PLoS Biology, vol. 3 no. 12, article e402,DOI: 10.1371/journal.pbio.0030402, 2005.

[71] C. K. Young, "Behavioral significance of hippocampal theta oscillations: looking elsewhere to find the right answers," Journal of Neurophysiology, vol. 106 no. 2, pp. 497-499, DOI: 10.1152/jn.00358.2011, 2011.

[72] V. D. Bohbot, M. Kalina, K. Stepankova, N. Spackova, M. Petrides, L. Nadel, "Spatial memory deficits in patients with lesions to the right hippocampus and to the right parahippocampal cortex," Neuropsychologia, vol. 36 no. 11, pp. 1217-1238, DOI: 10.1016/S0028-3932(97)00161-9, 1998.

[73] J. O'Keefe, N. Burgess, J. G. Donnett, K. J. Jeffery, E. A. Maguire, "Place cells, navigational accuracy, and the human hippocampus," Philosophical Transactions of the Royal Society B: Biological Sciences, vol. 353 no. 1373, pp. 1333-1340, DOI: 10.1098/rstb.1998.0287, 1998.

[74] H. J. Spiers, N. Burgess, E. A. Maguire, S. A. Baxendale, T. Hartley, P. J. Thompson, J. O'Keefe, "Unilateral temporal lobectomy patients show lateralized topographical and episodic memory deficits in a virtual town," Brain, vol. 124 no. 12, pp. 2476-2489, DOI: 10.1093/brain/124.12.2476, 2001.

[75] C. D. Tesche, J. Karhu, "Theta oscillations index human hippocampal activation during a working memory task," Proceedings of the National Academy of Sciences of the United States of America, vol. 97 no. 2, pp. 919-924, DOI: 10.1073/pnas.97.2.919, 2000.

[76] G. Buzsáki, "Theta rhythm of navigation: link between path integration and landmark navigation, episodic and semantic memory," Hippocampus, vol. 15 no. 7, pp. 827-840, DOI: 10.1002/hipo.20113, 2005.

[77] L. I. Aftanas, S. A. Golocheikine, "Human anterior and frontal midline theta and lower alpha reflect emotionally positive state and internalized attention: high-resolution EEG investigation of meditation," Neuroscience Letters, vol. 310 no. 1, pp. 57-60, DOI: 10.1016/S0304-3940(01)02094-8, 2001.

[78] D. L. Schacter, "EEG theta waves and psychological phenomena: a review and analysis," Biological Psychology, vol. 5 no. 1, pp. 47-82, DOI: 10.1016/0301-0511(77)90028-X, 1977.

[79] A. Ahani, H. Wahbeh, H. Nezamfar, M. Miller, D. Erdogmus, B. Oken, "Quantitative change of EEG and respiration signals during mindfulness meditation," Journal of NeuroEngineering and Rehabilitation, vol. 11 no. 1,DOI: 10.1186/1743-0003-11-87, 2014.