Identifying good practices for detecting inter-regional linear functional connectivity from EEG.

Aggregating voxel-level statistical dependencies between multivariate time series is an important intermediate step when characterising functional connectivity (FC) between larger brain regions. However, there are numerous ways in which voxel-level data can be aggregated into inter-regional FC, and the advantages of each of these approaches are currently unclear. In this study we generate ground-truth data and compare the performances of various pipelines that estimate directed and undirected linear phase-to-phase FC between regions. We test the ability of several existing and novel FC analysis pipelines to identify the true regions within which connectivity was simulated. We test various inverse modelling algorithms, strategies to aggregate time series within regions, and connectivity metrics. Furthermore, we investigate the influence of the number of interactions, the signal-to-noise ratio, the noise mix, the interaction time delay, and the number of active sources per region on the ability of detecting phase-to-phase FC. Throughout all simulated scenarios, lowest performance is obtained with pipelines involving the absolute value of coherency. Further, the combination of dynamic imaging of coherent sources (DICS) beamforming with directed FC metrics that aggregate information across multiple frequencies leads to unsatisfactory results. Pipelines that show promising results with our simulated pseudo-EEG data involve the following steps: (1) Source projection using the linearly-constrained minimum variance (LCMV) beamformer. (2) Principal component analysis (PCA) using the same fixed number of components within every region. (3) Calculation of the multivariate interaction measure (MIM) for every region pair to assess undirected phase-to-phase FC, or calculation of time-reversed Granger Causality (TRGC) to assess directed phase-to-phase FC. We formulate recommendations based on these results that may increase the validity of future experimental connectivity studies. We further introduce the free ROIconnect plugin for the EEGLAB toolbox that includes the recommended methods and pipelines that are presented here. We show an exemplary application of the best performing pipeline to the analysis of EEG data recorded during motor imagery.

Capturing the nature of events and event context using hierarchical event descriptors (HED).

Event-related data analysis plays a central role in EEG and MEG (MEEG) and other neuroimaging modalities including fMRI. Choices about which events to report and how to annotate their full natures significantly influence the value, reliability, and reproducibility of neuroimaging datasets for further analysis and meta- or mega-analysis. A powerful annotation strategy using the new third-generation formulation of the Hierarchical Event Descriptors (HED) framework and tools (hedtags.org) combines robust event description with details of experiment design and metadata in a human-readable as well as machine-actionable form, making event annotation relevant to the full range of neuroimaging and other time series data. This paper considers the event design and annotation process using as a case study the well-known multi-subject, multimodal dataset of Wakeman and Henson made available by its authors as a Brain Imaging Data Structure (BIDS) dataset (bids.neuroimaging.io). We propose a set of best practices and guidelines for event annotation integrated in a natural way into the BIDS metadata file architecture, examine the impact of event design decisions, and provide a working example of organizing events in MEEG and other neuroimaging data. We demonstrate how annotations using HED can document events occurring during neuroimaging experiments as well as their interrelationships, providing machine-actionable annotation enabling automated within- and across-experiment analysis and comparisons. We discuss the evolution of HED software tools and have made available an accompanying HED-annotated BIDS-formated edition of the MEEG data of the Wakeman and Henson dataset (openneuro.org, ds003645).

The open EEGLAB portal Interface: High-Performance computing with EEGLAB.

EEGLAB signal processing environment is currently the leading open-source software for processing electroencephalographic (EEG) data. The Neuroscience Gateway (NSG, nsgportal.org) is a web and API-based portal allowing users to easily run a variety of neuroscience-related software on high-performance computing (HPC) resources in the U.S. XSEDE network. We have reported recently (Delorme et al., 2019) on the Open EEGLAB Portal expansion of the free NSG services to allow the neuroscience community to build and run MATLAB pipelines using the EEGLAB tool environment. We are now releasing an EEGLAB plug-in, nsgportal, that interfaces EEGLAB with NSG directly from within EEGLAB running on MATLAB on any personal lab computer. The plug-in features a flexible MATLAB graphical user interface (GUI) that allows users to easily submit, interact with, and manage NSG jobs, and to retrieve and examine their results. Command line nsgportal tools supporting these GUI functionalities allow EEGLAB users and plug-in tool developers to build largely automated functions and workflows that include optional NSG job submission and processing. Here we present details on nsgportal implementation and documentation, provide user tutorials on example applications, and show sample test results comparing computation times using HPC versus laptop processing.

The role of midfrontal theta oscillations across the development of cognitive control in preschoolers and school-age children.

The development of cognitive control enables children to better resist acting based on distracting information that interferes with the current action. Cognitive control improvement serves different functions that differ in part by the type of interference to resolve. Indeed, resisting to interference at the task-set level or at the response-preparation level is, respectively, associated with cognitive flexibility and inhibition. It is, however, unknown whether the same neural mechanism underlies these two functions across development. Studies in adults have revealed the contribution of midfrontal theta (MFT) oscillations in interference resolution. This study investigated whether MFT is involved in the resolution of different types of interference in two age groups identified as corresponding to different latent structures of executive functions. Preschool (4-6 years) and school children (6-8 years) were tested with a task involving interference at the response level and/or the task-set level while (electroencephalogram) EEG was recorded. Behaviorally, response time and accuracy were affected by task-set. Both age groups were less accurate when the interference occurred at the task-set level and only the younger group showed decreased accuracy when interference was presented at the response-preparation level. Furthermore, MFT power was increased, relative to the baseline, during the resolution of both types of interference and in both age groups. These findings suggest that MFT is involved in immature cognitive control (i.e., preschool and school-ages), by orchestrating its different cognitive processes, irrespective of the interference to resolve and of the level of cognitive control development (i.e., the degree of differentiation of executive functions).

What Can Local Transfer Entropy Tell Us about Phase-Amplitude Coupling in Electrophysiological Signals?

Modulation of the amplitude of high-frequency cortical field activity locked to changes in the phase of a slower brain rhythm is known as phase-amplitude coupling (PAC). The study of this phenomenon has been gaining traction in neuroscience because of several reports on its appearance in normal and pathological brain processes in humans as well as across different mammalian species. This has led to the suggestion that PAC may be an intrinsic brain process that facilitates brain inter-area communication across different spatiotemporal scales. Several methods have been proposed to measure the PAC process, but few of these enable detailed study of its time course. It appears that no studies have reported details of PAC dynamics including its possible directional delay characteristic. Here, we study and characterize the use of a novel information theoretic measure that may address this limitation: local transfer entropy. We use both simulated and actual intracranial electroencephalographic data. In both cases, we observe initial indications that local transfer entropy can be used to detect the onset and offset of modulation process periods revealed by mutual information estimated phase-amplitude coupling (MIPAC). We review our results in the context of current theories about PAC in brain electrical activity, and discuss technical issues that must be addressed to see local transfer entropy more widely applied to PAC analysis. The current work sets the foundations for further use of local transfer entropy for estimating PAC process dynamics, and extends and complements our previous work on using local mutual information to compute PAC (MIPAC).

Human or not human? Performance monitoring ERPs during human agent and machine supervision.

Performance monitoring is a critical process which allows us to both learn from our own errors, and also interact with other human beings. However, our increasingly automated world requires us to interact more and more with automated systems, especially in risky environments. The present EEG study aimed at investigating and comparing the neuro-functional correlates associated with performance monitoring of an automated system and a human agent using a vertically-oriented arrowhead version of the flanker task. Given the influence of task difficulty on performance monitoring, two levels of difficulty were considered in order to assess their impact on supervision activity. A large N2-P3 complex in fronto-central regions was observed for both human agent error detection and system error detection during supervision. Using a cluster-based permutation analysis, a significantly decreased P3-like component was found for system compared to human agent error detection. This variation is in line with various psychosocial behavioral studies showing a difference between human-human and human-machine interactions, even though it was not clearly anticipated. Finally, the activity observed during error detection was significantly reduced in the difficult condition compared to the easy one, for both system and human agent supervision. Overall, this study is a first step towards the characterization of the neurophysiological correlates underlying system supervision, and a better understanding of their evolution in more complex environments. To go further, these results need to be replicated in other experiments with various paradigms to assess the robustness of the pattern and decrease during system supervision.

Measuring transient phase-amplitude coupling using local mutual information.

Here we demonstrate the suitability of a local mutual information measure for estimating the temporal dynamics of cross-frequency coupling (CFC) in brain electrophysiological signals. In CFC, concurrent activity streams in different frequency ranges interact and transiently couple. A particular form of CFC, phase-amplitude coupling (PAC), has raised interest given the growing amount of evidence of its possible role in healthy and pathological brain information processing. Although several methods have been proposed for PAC estimation, only a few have addressed the estimation of the temporal evolution of PAC, and these typically require a large number of experimental trials to return a reliable estimate. Here we explore the use of mutual information to estimate a PAC measure (MIPAC) in both continuous and event-related multi-trial data. To validate these two applications of the proposed method, we first apply it to a set of simulated phase-amplitude modulated signals and show that MIPAC can successfully recover the temporal dynamics of the simulated coupling in either continuous or multi-trial data. Finally, to explore the use of MIPAC to analyze data from human event-related paradigms, we apply it to an actual event-related human electrocorticographic (ECoG) data set that exhibits strong PAC, demonstrating that the MIPAC estimator can be used to successfully characterize amplitude-modulation dynamics in electrophysiological data.

Applying dimension reduction to EEG data by Principal Component Analysis reduces the quality of its subsequent Independent Component decomposition.

Independent Component Analysis (ICA) has proven to be an effective data driven method for analyzing EEG data, separating signals from temporally and functionally independent brain and non-brain source processes and thereby increasing their definition. Dimension reduction by Principal Component Analysis (PCA) has often been recommended before ICA decomposition of EEG data, both to minimize the amount of required data and computation time. Here we compared ICA decompositions of fourteen 72-channel single subject EEG data sets obtained (i) after applying preliminary dimension reduction by PCA, (ii) after applying no such dimension reduction, or else (iii) applying PCA only. Reducing the data rank by PCA (even to remove only 1% of data variance) adversely affected both the numbers of dipolar independent components (ICs) and their stability under repeated decomposition. For example, decomposing a principal subspace retaining 95% of original data variance reduced the mean number of recovered 'dipolar' ICs from 30 to 10 per data set and reduced median IC stability from 90% to 76%. PCA rank reduction also decreased the numbers of near-equivalent ICs across subjects. For instance, decomposing a principal subspace retaining 95% of data variance reduced the number of subjects represented in an IC cluster accounting for frontal midline theta activity from 11 to 5. PCA rank reduction also increased uncertainty in the equivalent dipole positions and spectra of the IC brain effective sources. These results suggest that when applying ICA decomposition to EEG data, PCA rank reduction should best be avoided.

Reduced mind wandering in experienced meditators and associated EEG correlates.

One outstanding question in the contemplative science literature relates to the direct impact of meditation experience on the monitoring of internal states and its respective correspondence with neural activity. In particular, to what extent does meditation influence the awareness, duration and frequency of the tendency of the mind to wander. To assess the relation between mind wandering and meditation, we tested 2 groups of meditators, one with a moderate level of experience (non-expert) and those who are well advanced in their practice (expert). We designed a novel paradigm using self-reports of internal mental states based on an experiential sampling probe paradigm presented during ~1 h of seated concentration meditation to gain insight into the dynamic measures of electroencephalography (EEG) during absorption in meditation as compared to reported mind wandering episodes. Our results show that expert meditation practitioners report a greater depth and frequency of sustained meditation, whereas non-expert practitioners report a greater depth and frequency of mind wandering episodes. This is one of the first direct behavioral indices of meditation expertise and its associated impact on the reduced frequency of mind wandering, with corresponding EEG activations showing increased frontal midline theta and somatosensory alpha rhythms during meditation as compared to mind wandering in expert practitioners. Frontal midline theta and somatosensory alpha rhythms are often observed during executive functioning, cognitive control and the active monitoring of sensory information. Our study thus provides additional new evidence to support the hypothesis that the maintenance of both internal and external orientations of attention may be maintained by similar neural mechanisms and that these mechanisms may be modulated by meditation training.

Influence of automation on mind wandering frequency in sustained attention.

There is accumulating evidence which shows that mind wandering may be increased within automated environments. This is particularly concerning when considering the negative effect of mind wandering on short-term performance. Seventeen participants performed an obstacle avoidance task under two conditions, manual and automated, each lasting 40 min. Subjects perceived the manual condition as more demanding than the automated one. We noted a significant increase of mind wandering frequency after only approximately 20 min under the automated condition. While learning and workload alone cannot explain these results, more automation-related phenomena, such as complacency or loss of agency, could play a role. Pupil diameter decreased during mind wandering compared to focus periods, revealing a decoupling from the task. The decrease remained stable in amplitude across different times and conditions. Research on mind wandering could be used to characterize an operator's state of mind regarding issues related to system interactions.

Electroencephalography correlates of spatial working memory deficits in attention-deficit/hyperactivity disorder: vigilance, encoding, and maintenance.

In the current study we sought to dissociate the component processes of working memory (WM) (vigilance, encoding and maintenance) that may be differentially impaired in attention-deficit/ hyperactivity disorder (ADHD). We collected electroencephalographic (EEG) data from 52 children with ADHD and 47 typically developing (TD) children, ages 7-14 years, while they performed a spatial Sternberg working memory task. We used independent component analysis and time-frequency analysis to identify midoccipital alpha (8-12 Hz) to evaluate encoding processes and frontal midline theta (4-7 Hz) to evaluate maintenance processes. We tested for effects of task difficulty and cue processing to evaluate vigilance. Children with ADHD showed attenuated alpha band event-related desynchronization (ERD) during encoding. This effect was more pronounced when task difficulty was low (consistent with impaired vigilance) and was predictive of memory task performance and symptom severity. Correlated with alpha ERD during encoding were alpha power increases during the maintenance period (relative to baseline), suggesting a compensatory effort. Consistent with this interpretation, midfrontal theta power increases during maintenance were stronger in ADHD and in high-load memory conditions. Furthermore, children with ADHD exhibited a maturational lag in development of posterior alpha power whereas age-related changes in frontal theta power deviated from the TD pattern. Last, subjects with ADHD showed age-independent attenuation of evoked responses to warning cues, suggesting low vigilance. Combined, these three EEG measures predicted diagnosis with 70% accuracy. We conclude that the interplay of impaired vigilance and encoding in ADHD may compromise maintenance and lead to impaired WM performance in this group.

RELICA: a method for estimating the reliability of independent components.

Independent Component Analysis (ICA) is a widely applied data-driven method for parsing brain and non-brain EEG source signals, mixed by volume conduction to the scalp electrodes, into a set of maximally temporally and often functionally independent components (ICs). Many ICs may be identified with a precise physiological or non-physiological origin. However, this process is hindered by partial instability in ICA results that can arise from noise in the data. Here we propose RELICA (RELiable ICA), a novel method to characterize IC reliability within subjects. RELICA first computes IC "dipolarity" a measure of physiological plausibility, plus a measure of IC consistency across multiple decompositions of bootstrap versions of the input data. RELICA then uses these two measures to visualize and cluster the separated ICs, providing a within-subject measure of IC reliability that does not involve checking for its occurrence across subjects. We demonstrate the use of RELICA on EEG data recorded from 14 subjects performing a working memory experiment and show that many brain and ocular artifact ICs are correctly classified as "stable" (highly repeatable across decompositions of bootstrapped versions of the input data). Many stable ICs appear to originate in the brain, while other stable ICs account for identifiable non-brain processes such as line noise. RELICA might be used with any linear blind source separation algorithm to reduce the risk of basing conclusions on unstable or physiologically un-interpretable component processes.

HRV and EEG correlates of well-being using ultra-short, portable, and low-cost measurements.

Wearable electroencephalography (EEG) and electrocardiography (ECG) devices may offer a non-invasive, user-friendly, and cost-effective approach for assessing well-being (WB) in real-world settings. However, challenges remain in dealing with signal artifacts (such as environmental noise and movements) and identifying robust biomarkers. We evaluated the feasibility of using portable hardware to identify potential EEG and heart-rate variability (HRV) correlates of WB. We collected simultaneous ultrashort (2-min) EEG and ECG data from 60 individuals in real-world settings using a wrist ECG electrode connected to a 4-channel wearable EEG headset. These data were processed, assessed for signal quality, and analyzed using the open-source EEGLAB BrainBeats plugin to extract several theory-driven metrics as potential correlates of WB. Namely, the individual alpha frequency (IAF), frontal and posterior alpha asymmetry, and signal entropy for EEG. SDNN, the low/high frequency (LF/HF) ratio, the Poincaré SD1/SD2 ratio, and signal entropy for HRV. We assessed potential associations between these features and the main WB dimensions (hedonic, eudaimonic, global, physical, and social) implementing a pairwise correlation approach, robust Spearman's correlations, and corrections for multiple comparisons. Only eight files showed poor signal quality and were excluded from the analysis. Eudaimonic (psychological) WB was positively correlated with SDNN and the LF/HF ratio. EEG posterior alpha asymmetry was positively correlated with Physical WB (i.e., sleep and pain levels). No relationships were found with the other metrics, or between EEG and HRV metrics. These physiological metrics enable a quick, objective assessment of well-being in real-world settings using scalable, user-friendly tools.

BrainBeats as an Open-Source EEGLAB Plugin to Jointly Analyze EEG and Cardiovascular Signals.

The interplay between the brain and the cardiovascular systems is garnering increased attention for its potential to advance our understanding of human physiology and improve health outcomes. However, the multimodal analysis of these signals is challenging due to the lack of guidelines, standardized signal processing and statistical tools, graphical user interfaces (GUIs), and automation for processing large datasets or increasing reproducibility. A further void exists in standardized EEG and heart-rate variability (HRV) feature extraction methods, undermining clinical diagnostics or the robustness of machine learning (ML) models. In response to these limitations, we introduce the BrainBeats toolbox. Implemented as an open-source EEGLAB plugin, BrainBeats integrates three main protocols: 1) Heartbeat-evoked potentials (HEP) and oscillations (HEO) for assessing time-locked brain-heart interplay at the millisecond accuracy; 2) EEG and HRV feature extraction for examining associations/differences between various brain and heart metrics or for building robust feature-based ML models; 3) Automated extraction of heart artifacts from EEG signals to remove any potential cardiovascular contamination while conducting EEG analysis. We provide a step-by-step tutorial for applying these three methods to an open-source dataset containing simultaneous 64-channel EEG, ECG, and PPG signals. Users can easily fine-tune parameters to tailor their unique research needs using the graphical user interface (GUI) or the command line. BrainBeats should make brain-heart interplay research more accessible and reproducible.

Who's calling? Evaluating the accuracy of guessing who is on the phone.

Some people claim to occasionally know who is calling them without using traditional means. Controlled experiments testing these claims report mixed results. We conducted a cross-sectional study of triads examining the accuracy of knowing who was calling using two randomly selected designs: 1) a web server randomly chose the caller before the callee's guess (telepathic/pre-selected trials), and 2) a web server randomly chose the caller after the callee's guess (precognitive/post-selected trials). We also performed exploratory multilevel mixed-effects logistic regressions on the relationship of genetic relationships, emotional closeness, communication frequency, and physical distance data with accuracy. A total of 177 participants completed at least one trial (105 "completers" completed all 12 trials). Accuracy was significantly above chance for the 210 completers telepathic/pre-selected trials (50.0% where the chance expectation was 33.3%, p<.001) but not the 630 completers precognitive/post-selected trials (31.9% where the chance expectation was 33.3%, p = .61). We discuss how these results favor the psi hypothesis, although conventional explanations cannot be completely excluded. Genetic relatedness significantly predicted accuracy in the regression model (Wald χ2 = 53.0, P < .001) for all trials. Compared to 0% genetic relatedness, the odds of accurately identifying the caller was 2.88 times (188%) higher for 25% genetic relatedness (Grandparent/Grandchild or Aunt/Uncle or Niece/Nephew or Half Sibling; ß = 1.06, z = 2.10, P = .04), but the other genetic relatedness levels were not significant. In addition, communication frequency was significant (ß = 0.006, z = 2.19, P = .03) but physical distance (ß = 0.0002, z = 1.56, P = .12) and emotional closeness (ß = 0.005, z = 1.87, P = .06) were not for all trials. To facilitate study recruitment and completion, unavoidable changes to the protocol were made during the study due to persistent recruitment difficulties, including changing inclusion/exclusion criteria, increasing total call attempts to participants, adjusting trial type randomization schema to ensure trial type balance, and participant compensation. Thus, future research will be needed to continue to improve the methodology and examine the mechanism by which people claim to know who is calling, as well as factors that may moderate the effects.

The Lab Streaming Layer for Synchronized Multimodal Recording.

Accurately recording the interactions of humans or other organisms with their environment or other agents requires synchronized data access via multiple instruments, often running independently using different clocks. Active, hardware-mediated solutions are often infeasible or prohibitively costly to build and run across arbitrary collections of input systems. The Lab Streaming Layer (LSL) offers a software-based approach to synchronizing data streams based on per-sample time stamps and time synchronization across a common LAN. Built from the ground up for neurophysiological applications and designed for reliability, LSL offers zero-configuration functionality and accounts for network delays and jitters, making connection recovery, offset correction, and jitter compensation possible. These features ensure precise, continuous data recording, even in the face of interruptions. The LSL ecosystem has grown to support over 150 data acquisition device classes as of Feb 2024, and establishes interoperability with and among client software written in several programming languages, including C/C++, Python, MATLAB, Java, C#, JavaScript, Rust, and Julia. The resilience and versatility of LSL have made it a major data synchronization platform for multimodal human neurobehavioral recording and it is now supported by a wide range of software packages, including major stimulus presentation tools, real-time analysis packages, and brain-computer interfaces. Outside of basic science, research, and development, LSL has been used as a resilient and transparent backend in scenarios ranging from art installations to stage performances, interactive experiences, and commercial deployments. In neurobehavioral studies and other neuroscience applications, LSL facilitates the complex task of capturing organismal dynamics and environmental changes using multiple data streams at a common timebase while capturing time details for every data frame.

The Past, Present, and Future of the Brain Imaging Data Structure (BIDS).

The Brain Imaging Data Structure (BIDS) is a community-driven standard for the organization of data and metadata from a growing range of neuroscience modalities. This paper is meant as a history of how the standard has developed and grown over time. We outline the principles behind the project, the mechanisms by which it has been extended, and some of the challenges being addressed as it evolves. We also discuss the lessons learned through the project, with the aim of enabling researchers in other domains to learn from the success of BIDS.

The past, present, and future of the brain imaging data structure (BIDS).

The Brain Imaging Data Structure (BIDS) is a community-driven standard for the organization of data and metadata from a growing range of neuroscience modalities. This paper is meant as a history of how the standard has developed and grown over time. We outline the principles behind the project, the mechanisms by which it has been extended, and some of the challenges being addressed as it evolves. We also discuss the lessons learned through the project, with the aim of enabling researchers in other domains to learn from the success of BIDS.

EEG is better left alone.

Automated preprocessing methods are critically needed to process the large publicly-available EEG databases, but the optimal approach remains unknown because we lack data quality metrics to compare them. Here, we designed a simple yet robust EEG data quality metric assessing the percentage of significant channels between two experimental conditions within a 100 ms post-stimulus time range. Because of volume conduction in EEG, given no noise, most brain-evoked related potentials (ERP) should be visible on every single channel. Using three publicly available collections of EEG data, we showed that, with the exceptions of high-pass filtering and bad channel interpolation, automated data corrections had no effect on or significantly decreased the percentage of significant channels. Referencing and advanced baseline removal methods were significantly detrimental to performance. Rejecting bad data segments or trials could not compensate for the loss in statistical power. Automated Independent Component Analysis rejection of eyes and muscles failed to increase performance reliably. We compared optimized pipelines for preprocessing EEG data maximizing ERP significance using the leading open-source EEG software: EEGLAB, FieldTrip, MNE, and Brainstorm. Only one pipeline performed significantly better than high-pass filtering the data.

Evaluating brain spectral and connectivity differences between silent mind-wandering and trance states.

Trance is an altered state of consciousness characterized by alterations in cognition. In general, trance states induce mental silence (i.e., cognitive thought reduction), and mental silence can induce trance states. Conversely, mind-wandering is the mind's propensity to stray its attention away from the task at hand and toward content irrelevant to the current moment, and its main component is inner speech. Building on the previous literature on mental silence and trance states and incorporating inverse source reconstruction advances, the study's objectives were to evaluate differences between trance and mind-wandering states using: (1) electroencephalography (EEG) power spectra at the electrode level, (2) power spectra at the area level (source reconstructed signal), and (3) EEG functional connectivity between these areas (i.e., how they interact). The relationship between subjective trance depths ratings and whole-brain connectivity during trance was also evaluated. Spectral analyses revealed increased delta and theta power in the frontal region and increased gamma in the centro-parietal region during mind-wandering, whereas trance showed increased beta and gamma power in the frontal region. Power spectra at the area level and pairwise comparisons of the connectivity between these areas demonstrated no significant difference between the two states. However, subjective trance depth ratings were inversely correlated with whole-brain connectivity in all frequency bands (i.e., deeper trance is associated with less large-scale connectivity). Trance allows one to enter mentally silent states and explore their neurophenomenological processes. Limitations and future directions are discussed.