Measuring functional connectivity using MEG: methodology and comparison with fcMRI

Using a linear dynamic system to measure functional connectivity from M/EEG

Objective.Measures of functional connectivity (FC) can elucidate which cortical regions work together in order to complete a variety of behavioral tasks. This study's primary objective was to expand a previously published model of measuring FC to include multiple subjects and several regions of interest. While FC has been more extensively investigated in vision and other sensorimotor tasks, it is not as well understood in audition. The secondary objective of this study was to investigate how auditory regions are functionally connected to other cortical regions when attention is directed to different distinct auditory stimuli.Approach.This study implements a linear dynamic system (LDS) to measure the structured time-lagged dependence across several cortical regions in order to estimate their FC during a dual-stream auditory attention task.Results.The model's output shows consistent functionally connected regions across different listening conditions, indicative of an auditory attention network that engages regardless of endogenous switching of attention or different auditory cues being attended.Significance.The LDS implemented in this study implements a multivariate autoregression to infer FC across cortical regions during an auditory attention task. This study shows how a first-order autoregressive function can reliably measure functional connectivity from M/EEG data. Additionally, the study shows how auditory regions engage with the supramodal attention network outlined in the visual attention literature.

Optimal gamma-band entrainment of visual cortex

Visual entrainment is a powerful and widely used research tool to study visual information processing in the brain. While many entrainment studies have focused on frequencies around 14-16 Hz, there is renewed interest in understanding visual entrainment at higher frequencies (e.g., gamma-band entrainment). Notably, recent groundbreaking studies have demonstrated that gamma-band visual entrainment at 40 Hz may have therapeutic effects in the context of Alzheimer's disease (AD) by stimulating specific neural ensembles, which utilize GABAergic signaling. Despite such promising findings, few studies have investigated the optimal parameters for gamma-band visual entrainment. Herein, we examined whether visual stimulation at 32, 40, or 48 Hz produces optimal visual entrainment responses using high-density magnetoencephalography (MEG). Our results indicated strong entrainment responses localizing to the primary visual cortex in each condition. Entrainment responses were stronger for 32 and 40 Hz relative to 48 Hz, indicating more robust synchronization of neural ensembles at these lower gamma-band frequencies. In addition, 32 and 40 Hz entrainment responses showed typical patterns of habituation across trials, but this effect was absent for 48 Hz. Finally, connectivity between visual cortex and parietal and prefrontal cortices tended to be strongest for 40 relative to 32 and 48 Hz entrainment. These results suggest that neural ensembles in the visual cortex may resonate at around 32 and 40 Hz and thus entrain more readily to photic stimulation at these frequencies. Emerging AD therapies, which have focused on 40 Hz entrainment to date, may be more effective at lower relative to higher gamma frequencies, although additional work in clinical populations is needed to confirm these findings. PRACTITIONER POINTS: Gamma-band visual entrainment has emerged as a therapeutic approach for eliminating amyloid in Alzheimer's disease, but its optimal parameters are unknown. We found stronger entrainment at 32 and 40 Hz compared to 48 Hz, suggesting neural ensembles prefer to resonate around these relatively lower gamma-band frequencies. These findings may inform the development and refinement of innovative AD therapies and the study of GABAergic visual cortical functions.

Spectral-switching analysis reveals real-time neuronal network representations of concurrent spontaneous naturalistic behaviors in human brain

Despite abundant evidence of functional networks in the human brain, their neuronal underpinnings, and relationships to real-time behavior have been challenging to resolve. Analyzing brain-wide intracranial-EEG recordings with video monitoring, acquired in awake subjects during clinical epilepsy evaluation, we discovered the tendency of each brain region to switch back and forth between 2 distinct power spectral densities (PSDs 2-55Hz). We further recognized that this 'spectral switching' occurs synchronously between distant sites, even between regions with differing baseline PSDs, revealing long-range functional networks that would be obscured in analysis of individual frequency bands. Moreover, the real-time PSD-switching dynamics of specific networks exhibited striking alignment with activities such as conversation and hand movements, revealing a multi-threaded functional network representation of concurrent naturalistic behaviors. Network structures and their relationships to behaviors were stable across days, but were altered during N3 sleep. Our results provide a new framework for understanding real-time, brain-wide neural-network dynamics.

Origin coordinate influence on performance of temporally extended signal space separation in magnetoencephalography

OBJECTIVE

Temporally extended signal space separation (tSSS) is a powerful method for artifact suppression in magnetoencephalography (MEG). Because tSSS first separates MEG signals coming from inside and outside a certain sphere, definition of the sphere origin is important. For this study, we explored the influence of origin choice on tSSS performance in spontaneous and evoked activity from epilepsy patients.

METHODS

Interictal epileptiform discharges (IEDs) and somatosensory evoked fields (SEFs) were processed with two tSSSs: one with the default origin of (0, 0, 40 mm) in the head coordinate, and the other with an individual origin estimated using each patient's anatomical magnetic resonance imaging (MRI). Equivalent current dipoles (ECDs) were calculated for the data. The ECD location and quality of estimation were compared across conditions.

RESULTS

MEG data from 21 patients revealed marginal differences in ECD location, but the estimation quality inferred from goodness of fit (GOF) and confidence volume (CV) was better for the tSSS with individual origins. This choice affected IEDs more than it affected SEFs.

CONCLUSIONS

Individual sphere model resulted in better GOF and CV.

SIGNIFICANCE

Application of tSSS using an individual origin would be more desirable when available. This parameter might influence spontaneous activity more strongly.

Beta: bursts of cognition

Beta oscillations are linked to the control of goal-directed processing of sensory information and the timing of motor output. Recent evidence demonstrates they are not sustained but organized into intermittent high-power bursts mediating timely functional inhibition. This implies there is a considerable moment-to-moment variation in the neural dynamics supporting cognition. Beta bursts thus offer new opportunities for studying how sensory inputs are selectively processed, reshaped by inhibitory cognitive operations and ultimately result in motor actions. Recent method advances reveal diversity in beta bursts that provide deeper insights into their function and the underlying neural circuit activity motifs. We propose that brain-wide, spatiotemporal patterns of beta bursting reflect various cognitive operations and that their dynamics reveal nonlinear aspects of cortical processing.

Individualized epidemic spreading models predict epilepsy surgery outcomes: A pseudo-prospective study

Epilepsy surgery is the treatment of choice for drug-resistant epilepsy patients, but up to 50% of patients continue to have seizures one year after the resection. In order to aid presurgical planning and predict postsurgical outcome on a patient-by-patient basis, we developed a framework of individualized computational models that combines epidemic spreading with patient-specific connectivity and epileptogeneity maps: the Epidemic Spreading Seizure and Epilepsy Surgery framework (ESSES). ESSES parameters were fitted in a retrospective study (N = 15) to reproduce invasive electroencephalography (iEEG)-recorded seizures. ESSES reproduced the iEEG-recorded seizures, and significantly better so for patients with good (seizure-free, SF) than bad (nonseizure-free, NSF) outcome. We illustrate here the clinical applicability of ESSES with a pseudo-prospective study (N = 34) with a blind setting (to the resection strategy and surgical outcome) that emulated presurgical conditions. By setting the model parameters in the retrospective study, ESSES could be applied also to patients without iEEG data. ESSES could predict the chances of good outcome after any resection by finding patient-specific model-based optimal resection strategies, which we found to be smaller for SF than NSF patients, suggesting an intrinsic difference in the network organization or presurgical evaluation results of NSF patients. The actual surgical plan overlapped more with the model-based optimal resection, and had a larger effect in decreasing modeled seizure propagation, for SF patients than for NSF patients. Overall, ESSES could correctly predict 75% of NSF and 80.8% of SF cases pseudo-prospectively. Our results show that individualised computational models may inform surgical planning by suggesting alternative resections and providing information on the likelihood of a good outcome after a proposed resection. This is the first time that such a model is validated with a fully independent cohort and without the need for iEEG recordings.

Neuronal Coupling Modes Show Differential Development in the Early Cortical Activity Networks of Human Newborns

The third trimester is a critical period for the development of functional networks that support the lifelong neurocognitive performance, yet the emergence of neuronal coupling in these networks is poorly understood. Here, we used longitudinal high-density electroencephalographic recordings from preterm infants during the period from 33 to 45 weeks of conceptional age (CA) to characterize early spatiotemporal patterns in the development of local cortical function and the intrinsic coupling modes [ICMs; phase-phase (PPCs), amplitude-amplitude (AACs), and phase-amplitude correlations (PACs)]. Absolute local power showed a robust increase with CA across the full frequency spectrum, while local PACs showed sleep state-specific, biphasic development that peaked a few weeks before normal birth. AACs and distant PACs decreased globally at nearly all frequencies. In contrast, the PPCs showed frequency- and region-selective development, with an increase of coupling strength with CA between frontal, central, and occipital regions at low-delta and alpha frequencies together with a wider-spread decrease at other frequencies. Our findings together present the spectrally and spatially differential development of the distinct ICMs during the neonatal period and provide their developmental templates for future basic and clinical research.

Tracking the neurodevelopmental trajectory of beta band oscillations with optically pumped magnetometer-based magnetoencephalography

Neural oscillations mediate the coordination of activity within and between brain networks, supporting cognition and behaviour. How these processes develop throughout childhood is not only an important neuroscientific question but could also shed light on the mechanisms underlying neurological and psychiatric disorders. However, measuring the neurodevelopmental trajectory of oscillations has been hampered by confounds from instrumentation. In this paper, we investigate the suitability of a disruptive new imaging platform - optically pumped magnetometer-based magnetoencephalography (OPM-MEG) - to study oscillations during brain development. We show how a unique 192-channel OPM-MEG device, which is adaptable to head size and robust to participant movement, can be used to collect high-fidelity electrophysiological data in individuals aged between 2 and 34 years. Data were collected during a somatosensory task, and we measured both stimulus-induced modulation of beta oscillations in sensory cortex, and whole-brain connectivity, showing that both modulate significantly with age. Moreover, we show that pan-spectral bursts of electrophysiological activity drive task-induced beta modulation, and that their probability of occurrence and spectral content change with age. Our results offer new insights into the developmental trajectory of beta oscillations and provide clear evidence that OPM-MEG is an ideal platform for studying electrophysiology in neurodevelopment.

The effects of spatial leakage correction on the reliability of EEG-based functional connectivity networks

Electroencephalography (EEG) functional connectivity (FC) estimates are confounded by the volume conduction problem. This effect can be greatly reduced by applying FC measures insensitive to instantaneous, zero-lag dependencies (corrected measures). However, numerous studies showed that FC measures sensitive to volume conduction (uncorrected measures) exhibit higher reliability and higher subject-level identifiability. We tested how source reconstruction contributed to the reliability difference of EEG FC measures on a large (n = 201) resting-state data set testing eight FC measures (including corrected and uncorrected measures). We showed that the high reliability of uncorrected FC measures in resting state partly stems from source reconstruction: idiosyncratic noise patterns define a baseline resting-state functional network that explains a significant portion of the reliability of uncorrected FC measures. This effect remained valid for template head model-based, as well as individual head model-based source reconstruction. Based on our findings we made suggestions how to best use spatial leakage corrected and uncorrected FC measures depending on the main goals of the study.

The Association Between Task Complexity and Cortical Language Mapping Accuracy

BACKGROUND AND OBJECTIVES

Direct cortical stimulation (DCS) mapping enables the identification of functional language regions within and around gliomas before tumor resection. Intraoperative mapping is required because glioma-infiltrated cortex engages in synchronous activity during task performance in a manner similar to normal-appearing cortex but has decreased ability to encode information for complex tasks. It is unknown whether task complexity influenced DCS mapping results. We aim to understand correlations between audiovisual picture naming (PN) task complexity and DCS error rate. We also asked what functional and oncological factors might be associated with higher rates of erroneous responses.

METHODS

We retrospectively reviewed intraoperative PN and word reading (WR) task performance during awake DCS language mapping for resection of dominant hemisphere World Health Organization grade 2 to 4 gliomas. The complexity of word tested in PN/WR tasks, patient characteristics, and tumor characteristics were compared between correct and incorrect trials.

RESULTS

Between 2017 and 2021, 74 patients met inclusion criteria. At median 18.6 months of follow-up, 73.0% were alive and 52.7% remained recurrence-free. A total of 2643 PN and 978 WR trials were analyzed. A greater number of syllables in PN was associated with a higher DCS error rate (P = .001). Multivariate logistic regression found that each additional syllable in PN tasks independently increased odds of error by 2.40 (P < .001). Older age was also an independent correlate of higher error rate (P < .043). World Health Organization grade did not correlate with error rate (P = .866). More severe language impairment before surgery correlated with worse performance on more complex intraoperative tasks (P < .001). A higher error rate on PN testing did not correlate with lower extent of glioma resection (P = .949).

CONCLUSION

Word complexity, quantified by the number of syllables, is associated with higher error rates for intraoperative PN tasks but does not affect extent of resection.

Exploring Motor Network Connectivity in State-Dependent Transcranial Magnetic Stimulation: A Proof-of-Concept Study

State-dependent non-invasive brain stimulation (NIBS) informed by electroencephalography (EEG) has contributed to the understanding of NIBS inter-subject and inter-session variability. While these approaches focus on local EEG characteristics, it is acknowledged that the brain exhibits an intrinsic long-range dynamic organization in networks. This proof-of-concept study explores whether EEG connectivity of the primary motor cortex (M1) in the pre-stimulation period aligns with the Motor Network (MN) and how the MN state affects responses to the transcranial magnetic stimulation (TMS) of M1. One thousand suprathreshold TMS pulses were delivered to the left M1 in eight subjects at rest, with simultaneous EEG. Motor-evoked potentials (MEPs) were measured from the right hand. The source space functional connectivity of the left M1 to the whole brain was assessed using the imaginary part of the phase locking value at the frequency of the sensorimotor μ-rhythm in a 1 s window before the pulse. Group-level connectivity revealed functional links between the left M1, left supplementary motor area, and right M1. Also, pulses delivered at high MN connectivity states result in a greater MEP amplitude compared to low connectivity states. At the single-subject level, this relation is more highly expressed in subjects that feature an overall high cortico-spinal excitability. In conclusion, this study paves the way for MN connectivity-based NIBS.

CiftiStorm pipeline: facilitating reproducible EEG/MEG source connectomics

We present CiftiStorm, an electrophysiological source imaging (ESI) pipeline incorporating recently developed methods to improve forward and inverse solutions. The CiftiStorm pipeline produces Human Connectome Project (HCP) and megconnectome-compliant outputs from dataset inputs with varying degrees of spatial resolution. The input data can range from low-sensor-density electroencephalogram (EEG) or magnetoencephalogram (MEG) recordings without structural magnetic resonance imaging (sMRI) to high-density EEG/MEG recordings with an HCP multimodal sMRI compliant protocol. CiftiStorm introduces a numerical quality control of the lead field and geometrical corrections to the head and source models for forward modeling. For the inverse modeling, we present a Bayesian estimation of the cross-spectrum of sources based on multiple priors. We facilitate ESI in the T1w/FSAverage32k high-resolution space obtained from individual sMRI. We validate this feature by comparing CiftiStorm outputs for EEG and MRI data from the Cuban Human Brain Mapping Project (CHBMP) acquired with technologies a decade before the HCP MEG and MRI standardized dataset.

Multispectral brain connectivity during visual attention distinguishes controlled from uncontrolled hypertension

Hypertension-related changes in brain function place individuals at higher risk for cognitive impairment and Alzheimer's disease. The existing functional neuroimaging literature has identified important neural and behavioural differences between normotensive and hypertensive individuals. However, previously-used methods (i.e. magnetic resonance imaging, functional near-infrared spectroscopy) rely on neurovascular coupling, which is a useful but indirect measure of neuronal activity. Furthermore, most studies fail to distinguish between controlled and uncontrolled hypertensive individuals, who exhibit significant behavioural and clinical differences. To partially remedy this gap in the literature, we used magnetoencephalography (MEG) to directly examine neuronal activity that is invariant to neurovascular coupling changes induced by hypertension. Our study included 52 participants (19 healthy controls, 15 controlled hypertensives, 18 uncontrolled hypertensives) who completed a modified flanker attention task during MEG. We identified significant oscillatory neural responses in two frequencies (alpha: 8-14 Hz, gamma: 48-60 Hz) for imaging and used grand-averaged images to determine seeds for whole-brain connectivity analysis. We then conducted Fisher-z tests for each pair of groups, using the relationship between the neural connectivity and behavioural attention effects. This highlighted a distributed network of regions associated with cognitive control and selective attention, including frontal-occipital and interhemispheric occipital connections. Importantly, the inferior frontal cortex exhibited a unique neurobehavioural relationship that distinguished the uncontrolled hypertensive group from the controlled hypertensive and normotensive groups. This is the first investigation of hypertension using MEG and identifies critical whole-brain connectivity differences based on hypertension profiles. KEY POINTS: Structural and functional changes in brain circuitry scale with hypertension severity and increase the risk of cognitive impairment and Alzheimer's disease. We harness the excellent spatiotemporal precision of magnetoencephalography (MEG) to directly quantify dynamic functional connectivity in healthy control, controlled hypertensive and uncontrolled hypertensive groups during a flanker task. In the first MEG study of hypertension, we show that there are neurobehavioural relationships that distinguish the uncontrolled hypertensive group from healthy and controlled hypertensive group in the prefrontal cortex. These results provide novel insights into the differential impact of hypertension on brain dynamics underlying selective attention.

Developmentally sensitive multispectral cortical connectivity profiles serving visual selective attention

Throughout childhood and adolescence, the brain undergoes significant structural and functional changes that contribute to the maturation of multiple cognitive domains, including selective attention. Selective attention is crucial for healthy executive functioning and while key brain regions serving selective attention have been identified, their age-related changes in neural oscillatory dynamics and connectivity remain largely unknown. We examined the developmental sensitivity of selective attention circuitry in 91 typically developing youth aged 6 - 13 years old. Participants completed a number-based Simon task while undergoing magnetoencephalography (MEG) and the resulting data were preprocessed and transformed into the time-frequency domain. Significant oscillatory brain responses were imaged using a beamforming approach, and task-related peak voxels in the occipital, parietal, and cerebellar cortices were used as seeds for subsequent whole-brain connectivity analyses in the alpha and gamma range. Our key findings revealed developmentally sensitive connectivity profiles in multiple regions crucial for selective attention, including the temporoparietal junction (alpha) and prefrontal cortex (gamma). Overall, these findings suggest that brain regions serving selective attention are highly sensitive to developmental changes during the pubertal transition period.

Motor learning- and consolidation-related resting state fast and slow brain dynamics across wake and sleep

Motor skills dynamically evolve during practice and after training. Using magnetoencephalography, we investigated the neural dynamics underpinning motor learning and its consolidation in relation to sleep during resting-state periods after the end of learning (boost window, within 30 min) and at delayed time scales (silent 4 h and next day 24 h windows) with intermediate daytime sleep or wakefulness. Resting-state neural dynamics were investigated at fast (sub-second) and slower (supra-second) timescales using Hidden Markov modelling (HMM) and functional connectivity (FC), respectively, and their relationship to motor performance. HMM results show that fast dynamic activities in a Temporal/Sensorimotor state network predict individual motor performance, suggesting a trait-like association between rapidly recurrent neural patterns and motor behaviour. Short, post-training task re-exposure modulated neural network characteristics during the boost but not the silent window. Re-exposure-related induction effects were observed on the next day, to a lesser extent than during the boost window. Daytime naps did not modulate memory consolidation at the behavioural and neural levels. These results emphasise the critical role of the transient boost window in motor learning and memory consolidation and provide further insights into the relationship between the multiscale neural dynamics of brain networks, motor learning, and consolidation.

The neurodevelopmental trajectory of beta band oscillations: an OPM-MEG study

Neural oscillations mediate the coordination of activity within and between brain networks, supporting cognition and behaviour. How these processes develop throughout childhood is not only an important neuroscientific question but could also shed light on the mechanisms underlying neurological and psychiatric disorders. However, measuring the neurodevelopmental trajectory of oscillations has been hampered by confounds from instrumentation. In this paper, we investigate the suitability of a disruptive new imaging platform - Optically Pumped Magnetometer-based magnetoencephalography (OPM-MEG) - to study oscillations during brain development. We show how a unique 192-channel OPM-MEG device, which is adaptable to head size and robust to participant movement, can be used to collect high-fidelity electrophysiological data in individuals aged between 2 and 34 years. Data were collected during a somatosensory task, and we measured both stimulus-induced modulation of beta oscillations in sensory cortex, and whole-brain connectivity, showing that both modulate significantly with age. Moreover, we show that pan-spectral bursts of electrophysiological activity drive task-induced beta modulation, and that their probability of occurrence and spectral content change with age. Our results offer new insights into the developmental trajectory of beta oscillations and provide clear evidence that OPM-MEG is an ideal platform for studying electrophysiology in neurodevelopment.

Healthy aging alters the oscillatory dynamics and fronto-parietal connectivity serving fluid intelligence

Fluid intelligence (Gf) involves logical reasoning and novel problem-solving abilities. Often, abstract reasoning tasks like Raven's progressive matrices are used to assess Gf. Prior work has shown an age-related decline in fluid intelligence capabilities, and although many studies have sought to identify the underlying mechanisms, our understanding of the critical brain regions and dynamics remains largely incomplete. In this study, we utilized magnetoencephalography (MEG) to investigate 78 individuals, ages 20-65 years, as they completed an abstract reasoning task. MEG data was co-registered with structural MRI data, transformed into the time-frequency domain, and the resulting neural oscillations were imaged using a beamformer. We found worsening behavioral performance with age, including prolonged reaction times and reduced accuracy. MEG analyses indicated robust oscillations in the theta, alpha/beta, and gamma range during the task. Whole brain correlation analyses with age revealed relationships in the theta and alpha/beta frequency bands, such that theta oscillations became stronger with increasing age in a right prefrontal region and alpha/beta oscillations became stronger with increasing age in parietal and right motor cortices. Follow-up connectivity analyses revealed increasing parieto-frontal connectivity with increasing age in the alpha/beta frequency range. Importantly, our findings are consistent with the parieto-frontal integration theory of intelligence (P-FIT). These results further suggest that as people age, there may be alterations in neural responses that are spectrally specific, such that older people exhibit stronger alpha/beta oscillations across the parieto-frontal network during abstract reasoning tasks.

Mindfulness meditation styles differently modulate source-level MEG microstate dynamics and complexity

Background

The investigation of mindfulness meditation practice, classically divided into focused attention meditation (FAM), and open monitoring meditation (OMM) styles, has seen a long tradition of theoretical, affective, neurophysiological and clinical studies. In particular, the high temporal resolution of magnetoencephalography (MEG) or electroencephalography (EEG) has been exploited to fill the gap between the personal experience of meditation practice and its neural correlates. Mounting evidence, in fact, shows that human brain activity is highly dynamic, transiting between different brain states (microstates). In this study, we aimed at exploring MEG microstates at source-level during FAM, OMM and in the resting state, as well as the complexity and criticality of dynamic transitions between microstates.

Methods

Ten right-handed Theravada Buddhist monks with a meditative expertise of minimum 2,265 h participated in the experiment. MEG data were acquired during a randomized block design task (6 min FAM, 6 min OMM, with each meditative block preceded and followed by 3 min resting state). Source reconstruction was performed using eLORETA on individual cortical space, and then parcellated according to the Human Connect Project atlas. Microstate analysis was then applied to parcel level signals in order to derive microstate topographies and indices. In addition, from microstate sequences, the Hurst exponent and the Lempel-Ziv complexity (LZC) were computed.

Results

Our results show that the coverage and occurrence of specific microstates are modulated either by being in a meditative state or by performing a specific meditation style. Hurst exponent values in both meditation conditions are reduced with respect to the value observed during rest, LZC shows significant differences between OMM, FAM, and REST, with a progressive increase from REST to FAM to OMM.

Discussion

Importantly, we report changes in brain criticality indices during meditation and between meditation styles, in line with a state-like effect of meditation on cognitive performance. In line with previous reports, we suggest that the change in cognitive state experienced in meditation is paralleled by a shift with respect to critical points in brain dynamics.

Penalized canonical correlation analysis reveals a relationship between temperament traits and brain oscillations during mind wandering

INTRODUCTION

There has been a growing interest in studying brain activity under naturalistic conditions. However, the relationship between individual differences in ongoing brain activity and psychological characteristics is not well understood. We investigated this connection, focusing on the association between oscillatory activity in the brain and individually characteristic dispositional traits. Given the variability of unconstrained resting states among individuals, we devised a paradigm that could harmonize the state of mind across all participants.

METHODS

We constructed task contrasts that included focused attention (FA), self-centered future planning, and rumination on anxious thoughts triggered by visual imagery. Magnetoencephalography was recorded from 28 participants under these 3 conditions for a duration of 16 min. The oscillatory power in the alpha and beta bands was converted into spatial contrast maps, representing the difference in brain oscillation power between the two conditions. We performed permutation cluster tests on these spatial contrast maps. Additionally, we applied penalized canonical correlation analysis (CCA) to study the relationship between brain oscillation patterns and behavioral traits.

RESULTS

The data revealed that the FA condition, as compared to the other conditions, was associated with higher alpha and beta power in the temporal areas of the left hemisphere and lower alpha and beta power in the parietal areas of the right hemisphere. Interestingly, the penalized CCA indicated that behavioral inhibition was positively correlated, whereas anxiety was negatively correlated, with a pattern of high oscillatory power in the bilateral precuneus and low power in the bilateral temporal regions. This unique association was found in the anxious-thoughts condition when contrasted with the focused-attention condition.

CONCLUSION

Our findings suggest individual temperament traits significantly affect brain engagement in naturalistic conditions. This research underscores the importance of considering individual traits in neuroscience and offers an effective method for analyzing brain activity and psychological differences.

Functional connectivity of sensorimotor network is enhanced in spastic diplegic cerebral palsy: A multimodal study using fMRI and MEG

OBJECTIVE

To assess the effects to functional connectivity (FC) caused by lesions related to spastic diplegic cerebral palsy (CP) in children and adolescents using multiple imaging modalities.

METHODS

We used resting state magnetoencephalography (MEG) envelope signals in alpha, beta and gamma ranges and resting state functional magnetic resonance imaging (fMRI) signals to quantify FC between selected sensorimotor regions of interest (ROIs) in 11 adolescents with spastic diplegic cerebral palsy and 24 typically developing controls. Motor performance of the hands was quantified with gross motor, fine motor and kinesthesia tests.

RESULTS

In fMRI, participants with CP showed enhanced FC within posterior parietal regions; in MEG, they showed enhanced interhemispheric FC between sensorimotor regions and posterior parietal regions both in alpha and lower beta bands. There was a correlation between the kinesthesia score and fronto-parietal connectivity in the control population.

CONCLUSIONS

CP is associated with enhanced FC in sensorimotor network. This difference is not correlated with hand coordination performance. The effect of the lesion is likely not fully captured by temporal correlation of ROI signals.

SIGNIFICANCE

Brain lesions can show as increased temporal correlation of activity between remote brain areas. We suggest this effect is likely separate from typical physiological correlates of functional connectivity.

Solving large-scale MEG/EEG source localisation and functional connectivity problems simultaneously using state-space models

State-space models are widely employed across various research disciplines to study unobserved dynamics. Conventional estimation techniques, such as Kalman filtering and expectation maximisation, offer valuable insights but incur high computational costs in large-scale analyses. Sparse inverse covariance estimators can mitigate these costs, but at the expense of a trade-off between enforced sparsity and increased estimation bias, necessitating careful assessment in low signal-to-noise ratio (SNR) situations. To address these challenges, we propose a three-fold solution: (1) Introducing multiple penalised state-space (MPSS) models that leverage data-driven regularisation; (2) Developing novel algorithms derived from backpropagation, gradient descent, and alternating least squares to solve MPSS models; (3) Presenting a K-fold cross-validation extension for evaluating regularisation parameters. We validate this MPSS regularisation framework through lower and more complex simulations under varying SNR conditions, including a large-scale synthetic magneto- and electro-encephalography (MEG/EEG) data analysis. In addition, we apply MPSS models to concurrently solve brain source localisation and functional connectivity problems for real event-related MEG/EEG data, encompassing thousands of sources on the cortical surface. The proposed methodology overcomes the limitations of existing approaches, such as constraints to small-scale and region-of-interest analyses. Thus, it may enable a more accurate and detailed exploration of cognitive brain functions.

Characterizing systemic physiological effects on the blood oxygen level dependent signal of resting-state fMRI in time-frequency space using wavelets

Systemic physiological dynamics, such as heart rate variability (HRV) and respiration volume per time (RVT), are known to account for significant variance in the blood oxygen level dependent (BOLD) signal of resting-state functional magnetic resonance imaging (rsfMRI). However, synchrony between these cardiorespiratory changes and the BOLD signal could be due to neuronal (i.e., autonomic activity inducing changes in heart rate and respiration) or vascular (i.e., cardiorespiratory activity facilitating hemodynamic changes and thus the BOLD signal) effects and the contributions of these effects may differ spatially, temporally, and spectrally. In this study, we characterize these brain-body dynamics using a wavelet analysis in rapidly sampled rsfMRI data with simultaneous pulse oximetry and respiratory monitoring of the Human Connectome Project. Our time-frequency analysis across resting-state networks (RSNs) revealed differences in the coherence of the BOLD signal and heartbeat interval (HBI)/RVT dynamics across frequencies, with unique profiles per network. Somatomotor (SMN), visual (VN), and salience (VAN) networks demonstrated the greatest synchrony with both systemic physiological signals when compared to other networks; however, significant coherence was observed in all RSNs regardless of direct autonomic involvement. Our phase analysis revealed distinct frequency profiles of percentage of time with significant coherence between BOLD and systemic physiological signals for different phase offsets across RSNs, suggesting that the phase offset and temporal order of signals varies by frequency. Lastly, our analysis of temporal variability of coherence provides insight on potential influence of autonomic state on brain-body communication. Overall, the novel wavelet analysis enables an efficient characterization of the dynamic relationship between cardiorespiratory activity and the BOLD signal in spatial, temporal, and spectral dimensions to inform our understanding of autonomic states and improve our interpretation of the BOLD signal.

Towards a biologically annotated brain connectome

The brain is a network of interleaved neural circuits. In modern connectomics, brain connectivity is typically encoded as a network of nodes and edges, abstracting away the rich biological detail of local neuronal populations. Yet biological annotations for network nodes - such as gene expression, cytoarchitecture, neurotransmitter receptors or intrinsic dynamics - can be readily measured and overlaid on network models. Here we review how connectomes can be represented and analysed as annotated networks. Annotated connectomes allow us to reconceptualize architectural features of networks and to relate the connection patterns of brain regions to their underlying biology. Emerging work demonstrates that annotated connectomes help to make more veridical models of brain network formation, neural dynamics and disease propagation. Finally, annotations can be used to infer entirely new inter-regional relationships and to construct new types of network that complement existing connectome representations. In summary, biologically annotated connectomes offer a compelling way to study neural wiring in concert with local biological features.

Classification of bipolar disorders using the multilayer modularity in dynamic minimum spanning tree from resting state fMRI

The diagnosis of bipolar disorders (BD) mainly depends on the clinical history and behavior observation, while only using clinical tools often limits the diagnosis accuracy. The study aimed to create a novel BD diagnosis framework using multilayer modularity in the dynamic minimum spanning tree (MST). We collected 45 un-medicated BD patients and 47 healthy controls (HC). The sliding window approach was utilized to construct dynamic MST via resting-state functional magnetic resonance imaging (fMRI) data. Firstly, we used three null models to explore the effectiveness of multilayer modularity in dynamic MST. Furthermore, the module allegiance exacted from dynamic MST was applied to train a classifier to discriminate BD patients. Finally, we explored the influence of the FC estimator and MST scale on the performance of the model. The findings indicated that multilayer modularity in the dynamic MST was not a random process in the human brain. And the model achieved an accuracy of 83.70% for identifying BD patients. In addition, we found the default mode network, subcortical network (SubC), and attention network played a key role in the classification. These findings suggested that the multilayer modularity in dynamic MST could highlight the difference between HC and BD patients, which opened up a new diagnostic tool for BD patients.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11571-022-09907-x.

Altered age-related alpha and gamma prefrontal-occipital connectivity serving distinct cognitive interference variants

The presence of conflicting stimuli adversely affects behavioral outcomes, which could either be at the level of stimulus (Flanker), response (Simon), or both (Multisource). Briefly, flanker interference involves conflicting stimuli requiring selective attention, Simon interference is caused by an incongruity between the spatial location of the task-relevant stimulus and prepotent motor mapping, and multisource is combination of both. Irrespective of the variant, interference resolution necessitates cognitive control to filter irrelevant information and allocate neural resources to task-related goals. Though previously studied in healthy young adults, the direct quantification of changes in oscillatory activity serving such cognitive control and associated inter-regional interactions in healthy aging are poorly understood. Herein, we used an adapted version of the multisource interference task and magnetoencephalography to investigate age-related alterations in the neural dynamics governing both divergent and convergent cognitive interference in 78 healthy participants (age range: 20-66 years). We identified weaker alpha connectivity between bilateral visual and right dorsolateral prefrontal cortices (DLPFC) and left dorsomedial prefrontal cortices (dmPFC), as well as weaker gamma connectivity between bilateral occipital regions and the right dmPFC during flanker interference with advancing age. Further, an age-related decrease in gamma power was observed in the left cerebellum and parietal region for Simon and differential interference effects (i.e., flanker-Simon), respectively. Moreover, the superadditivity model showed decreased gamma power in the right temporoparietal junction (TPJ) with increasing age. Overall, our findings suggest age-related declines in the engagement of top-down attentional control secondary to reduced alpha and gamma coupling between prefrontal and occipital cortices.

Variation in functional networks between clinical and subclinical discharges in childhood absence epilepsy: A multi-frequency MEG study

OBJECTIVE

Two types of spike-and-wave discharges (SWDs) exist in childhood absence epilepsy (CAE): clinical discharges are prolonged and manifest primarily as impaired consciousness, whereas subclinical discharges are brief with few objectively visible symptoms. This study aimed to compare neural functional network and default mode network (DMN) activity between clinical and subclinical discharges to better understand the underlying mechanism of CAE.

METHODS

Using magnetoencephalography (MEG) data from 21 patients, we obtained 25 segments each of clinical discharges and subclinical discharges. Amplitude envelope correlation analysis was used to construct functional networks and graph theory was used to calculate network topological data. We then compared differences in functional connectivity within the DMN between clinical and subclinical discharges. All statistical comparisons were performed using paired-sample tests.

RESULTS

Compared to subclinical discharges, the functional network of clinical discharges exhibited higher synchronization - particularly in the parahippocampal gyrus - as early as 10 s before the seizure. Additionally, the functional network of clinical SWDs presented an anterior shift of key nodes in the alpha frequency band. Regarding clinical discharge progression, there were gradual increases in the parameter node strengths (S), clustering coefficients (C), and global efficiency (E) of the functional networks, while the path lengths (L) decreased. These changes were most prominent at the onset of discharges and followed by some recovery in the high-frequency bands, but no significant change in the low-frequency bands. Furthermore, connections within the DMN during the discharge period were significantly stronger for clinical discharge compared to subclinical discharges.

CONCLUSIONS

These findings suggest that a more regular network before abnormal discharges in clinical discharges contributes to SWD explosion and that the parahippocampal gyrus plays an important role in maintaining oscillations. An absence seizure is a gradual process and the emergence of SWDs may be accompanied by initiation of inhibitory mechanisms. Enhanced functional connectivity among DMN brain regions may indicate that patients have entered a state of introspection, and functional abnormalities in the parahippocampal gyrus may be associated with patients' transient memory loss.

Texture-based brain networks for characterization of healthy subjects from MRI

Brain networks have been widely used to study the relationships between brain regions based on their dynamics using, e.g. fMRI or EEG, and to characterize their real physical connections using DTI. However, few studies have investigated brain networks derived from structural properties; and those have been based on cortical thickness or gray matter volume. The main objective of this work was to investigate the feasibility of obtaining useful information from brain networks derived from structural MRI, using texture features. We also wanted to verify if texture brain networks had any relation with established functional networks. T1-MR images were segmented using AAL and texture parameters from the gray-level co-occurrence matrix were computed for each region, for 760 subjects. Individual texture networks were used to evaluate the structural connections between regions of well-established functional networks; assess possible gender differences; investigate the dependence of texture network measures with age; and single out brain regions with different texture-network characteristics. Although around 70% of texture connections between regions belonging to the default mode, attention, and visual network were greater than the mean connection value, this effect was small (only between 7 and 15% of these connections were larger than one standard deviation), implying that texture-based morphology does not seem to subside function. This differs from cortical thickness-based morphology, which has been shown to relate to functional networks. Seventy-five out of 86 evaluated regions showed significant (ANCOVA, p < 0.05) differences between genders. Forty-four out of 86 regions showed significant (ANCOVA, p < 0.05) dependence with age; however, the R2 indicates that this is not a linear relation. Thalamus and putamen showed a very unique texture-wise structure compared to other analyzed regions. Texture networks were able to provide useful information regarding gender and age-related differences, as well as for singling out specific brain regions. We did not find a morphological texture-based subsidy for the evaluated functional brain networks. In the future, this approach will be extended to neurological patients to investigate the possibility of extracting biomarkers to help monitor disease evolution or treatment effectiveness.

Integrating multimodal and multiscale connectivity blueprints of the human cerebral cortex in health and disease

The brain is composed of disparate neural populations that communicate and interact with one another. Although fiber bundles, similarities in molecular architecture, and synchronized neural activity all reflect how brain regions potentially interact with one another, a comprehensive study of how all these interregional relationships jointly reflect brain structure and function remains missing. Here, we systematically integrate 7 multimodal, multiscale types of interregional similarity ("connectivity modes") derived from gene expression, neurotransmitter receptor density, cellular morphology, glucose metabolism, haemodynamic activity, and electrophysiology in humans. We first show that for all connectivity modes, feature similarity decreases with distance and increases when regions are structurally connected. Next, we show that connectivity modes exhibit unique and diverse connection patterns, hub profiles, spatial gradients, and modular organization. Throughout, we observe a consistent primacy of molecular connectivity modes-namely correlated gene expression and receptor similarity-that map onto multiple phenomena, including the rich club and patterns of abnormal cortical thickness across 13 neurological, psychiatric, and neurodevelopmental disorders. Finally, to construct a single multimodal wiring map of the human cortex, we fuse all 7 connectivity modes and show that the fused network maps onto major organizational features of the cortex including structural connectivity, intrinsic functional networks, and cytoarchitectonic classes. Altogether, this work contributes to the integrative study of interregional relationships in the human cerebral cortex.

Comparison of Resting-State EEG Network Analyses With and Without Parallel MRI in Genetic Generalized Epilepsy

Genetic generalized epilepsy (GGE) is conceptualized as a brain disorder involving distributed bilateral networks. To study these networks, simultaneous EEG-fMRI measurements can be used. However, inside-MRI EEG suffers from strong MR-related artifacts; it is not established whether EEG-based metrics in EEG-fMRI resting-state measurements are suitable for the analysis of group differences at source-level. We evaluated the impact of the inside-MR measurement condition on statistical group comparisons of EEG on source-level power and functional connectivity in patients with GGE versus healthy controls. We studied the cross-modal spatial relation of statistical group differences in seed-based FC derived from EEG and parallel fMRI. We found a significant increase in power and a frequency-specific change in functional connectivity for the inside MR-scanner compared to the outside MR-scanner condition. For power, we found reduced group difference between GGE and controls both in terms of statistical significance as well as effect size. Group differences for ImCoh remained similar both in terms of statistical significance as well as effect size. We found increased seed-based FC for GGE patients from the thalamus to the precuneus cortex region in fMRI, and in the theta band of simultaneous EEG. Our findings suggest that the analysis of EEG functional connectivity based on ImCoh is suitable for MR-EEG, and that relative group difference in a comparison of patients with GGE against controls are preserved. Spatial correspondence of seed-based FC group differences between the two modalities was found for the thalamus.

Spatially heterogeneous structure-function coupling in haemodynamic and electromagnetic brain networks

The relationship between structural and functional connectivity in the brain is a key question in connectomics. Here we quantify patterns of structure-function coupling across the neocortex, by comparing structural connectivity estimated using diffusion MRI with functional connectivity estimated using both neurophysiological (MEG-based) and haemodynamic (fMRI-based) recordings. We find that structure-function coupling is heterogeneous across brain regions and frequency bands. The link between structural and functional connectivity is generally stronger in multiple MEG frequency bands compared to resting state fMRI. Structure-function coupling is greater in slower and intermediate frequency bands compared to faster frequency bands. We also find that structure-function coupling systematically follows the archetypal sensorimotor-association hierarchy, as well as patterns of laminar differentiation, peaking in granular layer IV. Finally, structure-function coupling is better explained using structure-informed inter-regional communication metrics than using structural connectivity alone. Collectively, these results place neurophysiological and haemodynamic structure-function relationships in a common frame of reference and provide a starting point for a multi-modal understanding of structure-function coupling in the brain.

Multi-modal and multi-model interrogation of large-scale functional brain networks

Existing whole-brain models are generally tailored to the modelling of a particular data modality (e.g., fMRI or MEG/EEG). We propose that despite the differing aspects of neural activity each modality captures, they originate from shared network dynamics. Building on the universal principles of self-organising delay-coupled nonlinear systems, we aim to link distinct features of brain activity - captured across modalities - to the dynamics unfolding on a macroscopic structural connectome. To jointly predict connectivity, spatiotemporal and transient features of distinct signal modalities, we consider two large-scale models - the Stuart Landau and Wilson and Cowan models - which generate short-lived 40 Hz oscillations with varying levels of realism. To this end, we measure features of functional connectivity and metastable oscillatory modes (MOMs) in fMRI and MEG signals - and compare them against simulated data. We show that both models can represent MEG functional connectivity (FC), functional connectivity dynamics (FCD) and generate MOMs to a comparable degree. This is achieved by adjusting the global coupling and mean conduction time delay and, in the WC model, through the inclusion of balance between excitation and inhibition. For both models, the omission of delays dramatically decreased the performance. For fMRI, the SL model performed worse for FCD and MOMs, highlighting the importance of balanced dynamics for the emergence of spatiotemporal and transient patterns of ultra-slow dynamics. Notably, optimal working points varied across modalities and no model was able to achieve a correlation with empirical FC higher than 0.4 across modalities for the same set of parameters. Nonetheless, both displayed the emergence of FC patterns that extended beyond the constraints of the anatomical structure. Finally, we show that both models can generate MOMs with empirical-like properties such as size (number of brain regions engaging in a mode) and duration (continuous time interval during which a mode appears). Our results demonstrate the emergence of static and dynamic properties of neural activity at different timescales from networks of delay-coupled oscillators at 40 Hz. Given the higher dependence of simulated FC on the underlying structural connectivity, we suggest that mesoscale heterogeneities in neural circuitry may be critical for the emergence of parallel cross-modal functional networks and should be accounted for in future modelling endeavours.

Identification of Parkinson's Disease Subtypes from Resting-State Electroencephalography

BACKGROUND

Parkinson's disease (PD) patients present with a heterogeneous clinical phenotype, including motor, cognitive, sleep, and affective disruptions. However, this heterogeneity is often either ignored or assessed using only clinical assessments.

OBJECTIVES

We aimed to identify different PD sub-phenotypes in a longitudinal follow-up analysis and their electrophysiological profile based on resting-state electroencephalography (RS-EEG) and to assess their clinical significance over the course of the disease.

METHODS

Using electrophysiological features obtained from RS-EEG recordings and data-driven methods (similarity network fusion and source-space spectral analysis), we have performed a clustering analysis to identify disease sub-phenotypes and we examined whether their different patterns of disruption are predictive of disease outcome.

RESULTS

We showed that PD patients (n = 44) can be sub-grouped into three phenotypes with distinct electrophysiological profiles. These clusters are characterized by different levels of disruptions in the somatomotor network (Δ and β band), the frontotemporal network (α2 band) and the default mode network (α1 band), which consistently correlate with clinical profiles and disease courses. These clusters are classified into either moderate (only-motor) or mild-to-severe (diffuse) disease. We showed that EEG features can predict cognitive evolution of PD patients from baseline, when the cognitive clinical scores were overlapped.

CONCLUSIONS

The identification of novel PD subtypes based on electrical brain activity signatures may provide a more accurate prognosis in individual patients in clinical practice and help to stratify subgroups in clinical trials. Innovative profiling in PD can also support new therapeutic strategies that are brain-based and designed to modulate brain activity disruption. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Non-reversibility outperforms functional connectivity in characterisation of brain states in MEG data

Characterising brain states during tasks is common practice for many neuroscientific experiments using electrophysiological modalities such as electroencephalography (EEG) and magnetoencephalography (MEG). Brain states are often described in terms of oscillatory power and correlated brain activity, i.e. functional connectivity. It is, however, not unusual to observe weak task induced functional connectivity alterations in the presence of strong task induced power modulations using classical time-frequency representation of the data. Here, we propose that non-reversibility, or the temporal asymmetry in functional interactions, may be more sensitive to characterise task induced brain states than functional connectivity. As a second step, we explore causal mechanisms of non-reversibility in MEG data using whole brain computational models. We include working memory, motor, language tasks and resting-state data from participants of the Human Connectome Project (HCP). Non-reversibility is derived from the lagged amplitude envelope correlation (LAEC), and is based on asymmetry of the forward and reversed cross-correlations of the amplitude envelopes. Using random forests, we find that non-reversibility outperforms functional connectivity in the identification of task induced brain states. Non-reversibility shows especially better sensitivity to capture bottom-up gamma induced brain states across all tasks, but also alpha band associated brain states. Using whole brain computational models we find that asymmetry in the effective connectivity and axonal conduction delays play a major role in shaping non-reversibility across the brain. Our work paves the way for better sensitivity in characterising brain states during both bottom-up as well as top-down modulation in future neuroscientific experiments.

Identifying oscillatory brain networks with hidden Gaussian graphical spectral models of MEEG

Identifying the functional networks underpinning indirectly observed processes poses an inverse problem for neurosciences or other fields. A solution of such inverse problems estimates as a first step the activity emerging within functional networks from EEG or MEG data. These EEG or MEG estimates are a direct reflection of functional brain network activity with a temporal resolution that no other in vivo neuroimage may provide. A second step estimating functional connectivity from such activity pseudodata unveil the oscillatory brain networks that strongly correlate with all cognition and behavior. Simulations of such MEG or EEG inverse problem also reveal estimation errors of the functional connectivity determined by any of the state-of-the-art inverse solutions. We disclose a significant cause of estimation errors originating from misspecification of the functional network model incorporated into either inverse solution steps. We introduce the Bayesian identification of a Hidden Gaussian Graphical Spectral (HIGGS) model specifying such oscillatory brain networks model. In human EEG alpha rhythm simulations, the estimation errors measured as ROC performance do not surpass 2% in our HIGGS inverse solution and reach 20% in state-of-the-art methods. Macaque simultaneous EEG/ECoG recordings provide experimental confirmation for our results with 1/3 times larger congruence according to Riemannian distances than state-of-the-art methods.

Developmental differences in functional organization of multispectral networks

Assessing brain connectivity during rest has become a widely used approach to identify changes in functional brain organization during development. Generally, previous works have demonstrated that brain activity shifts from more local to more distributed processing from childhood into adolescence. However, the majority of those works have been based on functional magnetic resonance imaging measures, whereas multispectral functional connectivity, as measured using magnetoencephalography (MEG), has been far less characterized. In our study, we examined spontaneous cortical activity during eyes-closed rest using MEG in 101 typically developing youth (9-15 years old; 51 females, 50 males). Multispectral MEG images were computed, and connectivity was estimated in the canonical delta, theta, alpha, beta, and gamma bands using the imaginary part of the phase coherence, which was computed between 200 brain regions defined by the Schaefer cortical atlas. Delta and alpha connectivity matrices formed more communities as a function of increasing age. Connectivity weights predominantly decreased with age in both frequency bands; delta-band differences largely implicated limbic cortical regions and alpha band differences in attention and cognitive networks. These results are consistent with previous work, indicating the functional organization of the brain becomes more segregated across development, and highlight spectral specificity across different canonical networks.

Effects of acute psychosocial stress on source level EEG power and functional connectivity measures

The usage of EEG to uncover the influence of psychosocial stressors (PSSs) on neural activity has gained significant attention throughout recent years, but the results are often troubled by confounding stressor types. To investigate the effect of PSSs alone on neural activity, we employed a paradigm where participants are exposed to negative peer comparison as PSS, while other possible stressors are kept constant, and compared this with a condition where participants received neutral feedback. We analyzed commonly used sensor level EEG indices (frontal theta, alpha, and beta power) and further investigated whether source level power and functional connectivity (i.e., the temporal dependence between spatially seperated brain regions) measures, which have to our knowledge not yet been used, are more sensitive to PSSs than sensor level-derived EEG measures. Our results show that on sensor level, no significant frontal power changes are present (all p's > 0.16), indicating that sensor level frontal power measures are not sensitive enough to be affected by only PSSs. On source level, we find increased alpha power (indicative of decreased cortical activity) in the left- and right precuneus and right posterior cingulate cortex (all p's < 0.03) and increased functional connectivity between the left- and right precuneus (p < 0.001), indicating that acute, trial based PSSs lead to decreased precuneus/PCC activity, and possibly indicates a temporary disruption in the self-referential neural processes of an individual.

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.

Glioblastoma remodelling of human neural circuits decreases survival

Gliomas synaptically integrate into neural circuits1,2. Previous research has demonstrated bidirectional interactions between neurons and glioma cells, with neuronal activity driving glioma growth1-4 and gliomas increasing neuronal excitability2,5-8. Here we sought to determine how glioma-induced neuronal changes influence neural circuits underlying cognition and whether these interactions influence patient survival. Using intracranial brain recordings during lexical retrieval language tasks in awake humans together with site-specific tumour tissue biopsies and cell biology experiments, we find that gliomas remodel functional neural circuitry such that task-relevant neural responses activate tumour-infiltrated cortex well beyond the cortical regions that are normally recruited in the healthy brain. Site-directed biopsies from regions within the tumour that exhibit high functional connectivity between the tumour and the rest of the brain are enriched for a glioblastoma subpopulation that exhibits a distinct synaptogenic and neuronotrophic phenotype. Tumour cells from functionally connected regions secrete the synaptogenic factor thrombospondin-1, which contributes to the differential neuron-glioma interactions observed in functionally connected tumour regions compared with tumour regions with less functional connectivity. Pharmacological inhibition of thrombospondin-1 using the FDA-approved drug gabapentin decreases glioblastoma proliferation. The degree of functional connectivity between glioblastoma and the normal brain negatively affects both patient survival and performance in language tasks. These data demonstrate that high-grade gliomas functionally remodel neural circuits in the human brain, which both promotes tumour progression and impairs cognition.

Altered functional connectivity and oscillatory dynamics in polysubstance and cannabis only users during visuospatial processing

RATIONALE AND OBJECTIVES

Cannabis use is often associated with the use of other psychoactive substances, which is subsequently linked to an increased risk for addiction. While there is a growing body of neuroimaging literature investigating the cognitive effect of long-term cannabis use, very little is known about the potential additive effects of cannabis polysubstance use.

METHODS

Fifty-six adults composed of 18 polysubstance users (i.e., cannabis plus at least one other illicit substance), 19 cannabis-only users, and 19 nonusers completed a visuospatial attention task while undergoing magnetoencephalography. A data-driven approach was used to identify oscillatory neural responses, which were imaged using a beamforming approach. The resulting cortical regions were probed for group differences and used as seeds for whole-brain connectivity analysis.

RESULTS

Participants exhibited robust theta, alpha, beta, and gamma responses during visuospatial processing. Statistical analyses indicated that the cannabis-only group had weaker occipital theta relative to the nonusers, and that both polysubstance and cannabis-only users had reduced spontaneous gamma in the occipital cortices during the pre-stimulus baseline period relative to nonusers. Finally, functional connectivity analyses revealed that polysubstance users had sharply reduced beta connectivity between occipital and prefrontal, as well as occipital and left temporal cortices.

CONCLUSIONS

Cannabis use should be considered in a polysubstance context, as our correlational design suggests differences in functional connectivity among those who reported cannabis-only versus polysubstance use in occipital to prefrontal pathways critical to visuospatial processing and attention function. Future work should distinguish the effect of different polysubstance combinations and use more causal designs.

Dissociation between phase and power correlation networks in the human brain is driven by co-occurrent bursts

Well-known haemodynamic resting-state networks are better mirrored in power correlation networks than phase coupling networks in electrophysiological data. However, what do these power correlation networks reflect? We address this long-outstanding question in neuroscience using rigorous mathematical analysis, biophysical simulations with ground truth and application of these mathematical concepts to empirical magnetoencephalography (MEG) data. Our mathematical derivations show that for two non-Gaussian electrophysiological signals, their power correlation depends on their coherence, cokurtosis and conjugate-coherence. Only coherence and cokurtosis contribute to power correlation networks in MEG data, but cokurtosis is less affected by artefactual signal leakage and better mirrors haemodynamic resting-state networks. Simulations and MEG data show that cokurtosis may reflect co-occurrent bursting events. Our findings shed light on the origin of the complementary nature of power correlation networks to phase coupling networks and suggests that the origin of resting-state networks is partly reflected in co-occurent bursts in neuronal activity.

Altered neuromagnetic activity in default mode network in childhood absence epilepsy

Purpose

The electrophysiological characterization of resting state oscillatory functional connectivity within the default mode network (DMN) during interictal periods in childhood absence epilepsy (CAE) remains unclear. Using magnetoencephalographic (MEG) recordings, this study investigated how the connectivity within the DMN was altered in CAE.

Methods

Using a cross-sectional design, we analyzed MEG data from 33 children newly diagnosed with CAE and 26 controls matched for age and sex. The spectral power and functional connectivity of the DMN were estimated using minimum norm estimation combined with the Welch technique and corrected amplitude envelope correlation.

Results

Default mode network showed stronger activation in the delta band during the ictal period, however, the relative spectral power in other bands was significantly lower than that in the interictal period (pcorrected &lt; 0.05 for DMN regions, except bilateral medial frontal cortex, left medial temporal lobe, left posterior cingulate cortex in the theta band, and the bilateral precuneus in the alpha band). It should be noted that the significant power peak in the alpha band was lost compared with the interictal data. Compared with controls, the interictal relative spectral power of DMN regions (except bilateral precuneus) in CAE patients was significantly increased in the delta band (pcorrected &lt; 0.01), whereas the values of all DMN regions in the beta-gamma 2 band were significantly decreased (pcorrected &lt; 0.01). In the higher frequency band (alpha-gamma1), especially in the beta and gamma1 band, the ictal node strength of DMN regions except the left precuneus was significantly higher than that in the interictal periods (pcorrected &lt; 0.01), and the node strength of the right inferior parietal lobe increased most significantly in the beta band (Ictal: 3.8712 vs. Interictal: 0.7503, pcorrected &lt; 0.01). Compared with the controls, the interictal node strength of DMN increased in all frequency bands, especially the right medial frontal cortex in the beta band (Controls: 0.1510 vs. Interictal: 3.527, pcorrected &lt; 0.01). Comparing relative node strength between groups, the right precuneus in CAE children decreased significantly (β: Controls: 0.1009 vs. Interictal: 0.0475; γ 1: Controls:0.1149 vs. Interictal:0.0587, pcorrected &lt; 0.01) such that it was no longer the central hub.

Conclusion

These findings indicated DMN abnormalities in CAE patients, even in interictal periods without interictal epileptic discharges. Abnormal functional connectivity in CAE may reflect abnormal anatomo-functional architectural integration in DMN, as a result of cognitive mental impairment and unconsciousness during absence seizure. Future studies are needed to examine if the altered functional connectivity can be used as a biomarker for treatment responses, cognitive dysfunction, and prognosis in CAE patients.

Minimizing the distortions in electrophysiological source imaging of cortical oscillatory activity via Spectral Structured Sparse Bayesian Learning

Oscillatory processes at all spatial scales and on all frequencies underpin brain function. Electrophysiological Source Imaging (ESI) is the data-driven brain imaging modality that provides the inverse solutions to the source processes of the EEG, MEG, or ECoG data. This study aimed to carry out an ESI of the source cross-spectrum while controlling common distortions of the estimates. As with all ESI-related problems under realistic settings, the main obstacle we faced is a severely ill-conditioned and high-dimensional inverse problem. Therefore, we opted for Bayesian inverse solutions that posited a priori probabilities on the source process. Indeed, rigorously specifying both the likelihoods and a priori probabilities of the problem leads to the proper Bayesian inverse problem of cross-spectral matrices. These inverse solutions are our formal definition for cross-spectral ESI (cESI), which requires a priori of the source cross-spectrum to counter the severe ill-condition and high-dimensionality of matrices. However, inverse solutions for this problem were NP-hard to tackle or approximated within iterations with bad-conditioned matrices in the standard ESI setup. We introduce cESI with a joint a priori probability upon the source cross-spectrum to avoid these problems. cESI inverse solutions are low-dimensional ones for the set of random vector instances and not random matrices. We achieved cESI inverse solutions through the variational approximations via our Spectral Structured Sparse Bayesian Learning (ssSBL) algorithm https://github.com/CCC-members/Spectral-Structured-Sparse-Bayesian-Learning. We compared low-density EEG (10-20 system) ssSBL inverse solutions with reference cESIs for two experiments: (a) high-density MEG that were used to simulate EEG and (b) high-density macaque ECoG that were recorded simultaneously with EEG. The ssSBL resulted in two orders of magnitude with less distortion than the state-of-the-art ESI methods. Our cESI toolbox, including the ssSBL method, is available at https://github.com/CCC-members/BC-VARETA_Toolbox.

Effect of channel density, inverse solutions and connectivity measures on EEG resting-state networks reconstruction: A simulation study

Along with the study of brain activity evoked by external stimuli, the past two decades witnessed an increased interest in characterizing the spontaneous brain activity occurring during resting conditions. The identification of connectivity patterns in this so-called "resting-state" has been the subject of a great number of electrophysiology-based studies, using the Electro/Magneto-Encephalography (EEG/MEG) source connectivity method. However, no consensus has been reached yet regarding a unified (if possible) analysis pipeline, and several involved parameters and methods require cautious tuning. This is particularly challenging when different analytical choices induce significant discrepancies in results and drawn conclusions, thereby hindering the reproducibility of neuroimaging research. Hence, our objective in this study was to shed light on the effect of analytical variability on outcome consistency by evaluating the implications of parameters involved in the EEG source connectivity analysis on the accuracy of resting-state networks (RSNs) reconstruction. We simulated, using neural mass models, EEG data corresponding to two RSNs, namely the default mode network (DMN) and dorsal attentional network (DAN). We investigated the impact of five channel densities (19, 32, 64, 128, 256), three inverse solutions (weighted minimum norm estimate (wMNE), exact low-resolution brain electromagnetic tomography (eLORETA), and linearly constrained minimum variance (LCMV) beamforming) and four functional connectivity measures (phase-locking value (PLV), phase-lag index (PLI), and amplitude envelope correlation (AEC) with and without source leakage correction), on the correspondence between reconstructed and reference networks. We showed that, with different analytical choices related to the number of electrodes, source reconstruction algorithm, and functional connectivity measure, high variability is present in the results. More specifically, our results show that a higher number of EEG channels significantly increased the accuracy of the reconstructed networks. Additionally, our results showed significant variability in the performance of the tested inverse solutions and connectivity measures. Such methodological variability and absence of analysis standardization represent a critical issue for neuroimaging studies that should be prioritized. We believe that this work could be useful for the field of electrophysiology connectomics, by increasing awareness regarding the challenge of variability in methodological approaches and its implications on reported results.

Functional brain network properties of healthy full-term newborns quantified by scalp and source-reconstructed EEG

OBJECTIVE

Identifying the functional brain network properties of term low-risk newborns using high-density EEG (HD-EEG) and comparing these properties with those of established functional magnetic resonance image (fMRI) - based networks.

METHODS

HD-EEG was collected from 113 low-risk term newborns before delivery hospital discharge and within 72 hours of birth. Functional brain networks were reconstructed using coherence at the scalp and source levels in delta, theta, alpha, beta, and gamma frequency bands. These networks were characterized for the global and local network architecture.

RESULTS

Source-level networks in all the frequency bands identified the presence of the efficient small world (small-world propensity (SWP) > 0.6) architecture with four distinct modules linked by hub regions and rich-club (coefficient > 1) topology. The modular regions included primary, association, limbic, paralimbic, and subcortical regions, which have been demonstrated in fMRI studies. In contrast, scalp-level networks did not display consistent small world architecture (SWP < 0.6), and also identified only 2-3 modules in each frequency band.The modular regions of the scalp-network primarily included frontal and occipital regions.

CONCLUSIONS

Our findings show that EEG sources in low-risk newborns corroborate fMRI-based connectivity results.

SIGNIFICANCE

EEG source analysis characterizes functional connectivity at the bedside of low-risk newborn infants soon after birth.

Reliable evaluation of functional connectivity and graph theory measures in source-level EEG: How many electrodes are enough?

OBJECTIVE

Using EEG to characterise functional brain networks through graph theory has gained significant interest in clinical and basic research. However, the minimal requirements for reliable measures remain largely unaddressed. Here, we examined functional connectivity estimates and graph theory metrics obtained from EEG with varying electrode densities.

METHODS

EEG was recorded with 128 electrodes in 33 participants. The high-density EEG data were subsequently subsampled into three sparser montages (64, 32, and 19 electrodes). Four inverse solutions, four measures of functional connectivity, and five graph theory metrics were tested.

RESULTS

The correlation between the results obtained with 128-electrode and the subsampled montages decreased as a function of the number of electrodes. As a result of decreased electrode density, the network metrics became skewed: mean network strength and clustering coefficient were overestimated, while characteristic path length was underestimated.

CONCLUSIONS

Several graph theory metrics were altered when electrode density was reduced. Our results suggest that, for optimal balance between resource demand and result precision, a minimum of 64 electrodes should be utilised when graph theory metrics are used to characterise functional brain networks in source-reconstructed EEG data.

SIGNIFICANCE

Characterisation of functional brain networks derived from low-density EEG warrants careful consideration.

Functional connectivity fingerprints of the frontal eye field and inferior frontal junction suggest spatial versus nonspatial processing in the prefrontal cortex

Neuroimaging evidence suggests that the frontal eye field (FEF) and inferior frontal junction (IFJ) govern the encoding of spatial and nonspatial (such as feature- or object-based) representations, respectively, both during visual attention and working memory tasks. However, it is still unclear whether such contrasting functional segregation is also reflected in their underlying functional connectivity patterns. Here, we hypothesized that FEF has predominant functional coupling with spatiotopically organized regions in the dorsal ('where') visual stream whereas IFJ has predominant functional connectivity with the ventral ('what') visual stream. We applied seed-based functional connectivity analyses to temporally high-resolving resting-state magnetoencephalography (MEG) recordings. We parcellated the brain according to the multimodal Glasser atlas and tested, for various frequency bands, whether the spontaneous activity of each parcel in the ventral and dorsal visual pathway has predominant functional connectivity with FEF or IFJ. The results show that FEF has a robust power correlation with the dorsal visual pathway in beta and gamma bands. In contrast, anterior IFJ (IFJa) has a strong power coupling with the ventral visual stream in delta, beta and gamma oscillations. Moreover, while FEF is phase-coupled with the superior parietal lobe in the beta band, IFJa is phase-coupled with the middle and inferior temporal cortex in delta and gamma oscillations. We argue that these intrinsic connectivity fingerprints are congruent with each brain region's function. Therefore, we conclude that FEF and IFJ have dissociable connectivity patterns that fit their respective functional roles in spatial versus nonspatial top-down attention and working memory control.

The neural correlates of amplitude of low-frequency fluctuation: a multimodal resting-state MEG and fMRI-EEG study

The amplitude of low-frequency fluctuation (ALFF) describes the regional intensity of spontaneous blood-oxygen-level-dependent signal in resting-state functional magnetic resonance imaging (fMRI). How the fMRI-ALFF relates to the amplitude in electrophysiological signals remains unclear. We here aimed to investigate the neural correlates of fMRI-ALFF by comparing the spatial difference of amplitude between the eyes-closed (EC) and eyes-open (EO) states from fMRI and magnetoencephalography (MEG), respectively. By synthesizing MEG signal into amplitude-based envelope time course, we first investigated 2 types of amplitude in MEG, meaning the amplitude of neural activities from delta to gamma (i.e. MEG-amplitude) and the amplitude of their low-frequency modulation at the fMRI range (i.e. MEG-ALFF). We observed that the MEG-ALFF in EC was increased at parietal sensors, ranging from alpha to beta; whereas the MEG-amplitude in EC was increased at the occipital sensors in alpha. Source-level analysis revealed that the increased MEG-ALFF in the sensorimotor cortex overlapped with the most reliable EC-EO differences observed in fMRI at slow-3 (0.073-0.198 Hz), and these differences were more significant after global mean standardization. Taken together, our results support that (i) the amplitude at 2 timescales in MEG reflect distinct physiological information and that (ii) the fMRI-ALFF may relate to the ALFF in neural activity.

Brain network mapping and glioma pathophysiology

Adult diffuse gliomas are among the most difficult brain disorders to treat in part due to a lack of clarity regarding the anatomical origins and mechanisms of migration of the tumours. While the importance of studying networks of glioma spread has been recognized for at least 80 years, the ability to carry out such investigations in humans has emerged only recently. Here, we comprehensively review the fields of brain network mapping and glioma biology to provide a primer for investigators interested in merging these areas of inquiry for the purposes of translational research. Specifically, we trace the historical development of ideas in both brain network mapping and glioma biology, highlighting studies that explore clinical applications of network neuroscience, cells-of-origin of diffuse glioma and glioma-neuronal interactions. We discuss recent research that has merged neuro-oncology and network neuroscience, finding that the spatial distribution patterns of gliomas follow intrinsic functional and structural brain networks. Ultimately, we call for more contributions from network neuroimaging to realize the translational potential of cancer neuroscience.

Virtual MEG sensors based on beamformer and independent component analysis can reconstruct epileptic activity as measured on simultaneous intracerebral recordings

The prevailing gold standard for presurgical determination of epileptogenic brain networks is intracerebral EEG, a potent yet invasive approach. Magnetoencephalography (MEG) is a state-of-the art non-invasive method for investigating epileptiform discharges. However, it is not clear at what level the precision offered by MEG can reach that of SEEG. Here, we present a strategy for non-invasively retrieving the constituents of the interictal network, with high spatial and temporal precision. Our method is based on MEG and a combination of spatial filtering and independent component analysis (ICA). We validated this approach in twelve patients with drug-resistant focal epilepsy, thanks to the unprecedented ground truth provided by simultaneous recordings of MEG and SEEG. A minimum variance adaptive beamformer estimated the source time series and ICA was used to further decompose these time series into network constituents (MEG-ICs), each having a time series (virtual electrode) and a topography (spatial distribution of amplitudes in the brain). We show that MEG has a considerable sensitivity of 0.80 and 0.84 and a specificity of 0.93 and 0.91 for reconstructing deep and superficial sources, respectively, when compared to the ground truth (SEEG). For each epileptic MEG-IC (n = 131), we found at least one significantly correlating SEEG contact close to zero lag after correcting for multiple comparisons. All the patients except one had at least one epileptic component that was highly correlated (Spearman rho>0.3) with that of SEEG traces. MEG-ICs correlated well with SEEG traces. The strength of correlation coefficients did not depend on the depth of the SEEG contacts or the clinical outcome of the patient. A significant proportion of the MEG-ICs (n = 83/131) were localized in proximity with their maximally correlating SEEG, within a mean distance of 20±12.18mm. Our research is the first to validate the MEG-retrieved beamformer IC sources against SEEG-derived ground truth in a simultaneous MEG-SEEG framework. Observations from the present study suggest that non-invasive MEG source components may potentially provide additional information, comparable to SEEG in a number of instances.

An orderly sequence of autonomic and neural events at transient arousal changes

Resting-state functional magnetic resonance imaging (rsfMRI) allows the study of functional brain connectivity based on spatially structured variations in neuronal activity. Proper evaluation of connectivity requires removal of non-neural contributions to the fMRI signal, in particular hemodynamic changes associated with autonomic variability. Regression analysis based on autonomic indicator signals has been used for this purpose, but may be inadequate if neuronal and autonomic activities covary. To investigate this potential co-variation, we performed rsfMRI experiments while concurrently acquiring electroencephalography (EEG) and autonomic indicator signals, including heart rate, respiratory depth, and peripheral vascular tone. We identified a recurrent and systematic spatiotemporal pattern of fMRI (named as fMRI cascade), which features brief signal reductions in salience and default-mode networks and the thalamus, followed by a biphasic global change with a sensory-motor dominance. This fMRI cascade, which was mostly observed during eyes-closed condition, was accompanied by large EEG and autonomic changes indicative of arousal modulations. Importantly, the removal of the fMRI cascade dynamics from rsfMRI diminished its correlations with various signals. These results suggest that the rsfMRI correlations with various physiological and neural signals are not independent but arise, at least partly, from the fMRI cascades and associated neural and physiological changes at arousal modulations.