Frequency-dependent functional connectivity within resting-state networks: an atlas-based MEG beamformer solution

FREQ-NESS Reveals the Dynamic Reconfiguration of Frequency-Resolved Brain Networks During Auditory Stimulation

The brain is a dynamic system whose network organization is often studied by focusing on specific frequency bands or anatomical regions, leading to fragmented insights, or by employing complex and elaborate methods that hinder straightforward interpretations. To address this issue, a new analytical pipeline named FREQuency-resolved Network Estimation via Source Separation (FREQ-NESS) is introduced. This pipeline is designed to estimate the activation and spatial configuration of simultaneous brain networks across frequencies by analyzing the frequency-resolved multivariate covariance between whole-brain voxel time series. In this study, FREQ-NESS is applied to source-reconstructed magnetoencephalography (MEG) data during resting state and isochronous auditory stimulation. Our results reveal simultaneous, frequency-specific brain networks during resting state, such as the default mode, alpha-band, and motor-beta networks. During auditory stimulation, FREQ-NESS detects: 1) emergence of networks attuned to the stimulation frequency, 2) spatial reorganization of existing networks, such as alpha-band networks shifting from occipital to sensorimotor areas, 3) stability of networks unaffected by auditory stimuli. Furthermore, auditory stimulation significantly enhances cross-frequency coupling, with the phase of auditory networks attuned to the stimulation modulating gamma band amplitude in medial temporal lobe networks. In conclusion, FREQ-NESS effectively maps the brain's spatiotemporal dynamics, providing a comprehensive view of brain function by revealing a landscape of simultaneous, frequency-resolved networks and their interaction.

How do the resting EEG preprocessing states affect the outcomes of postprocessing?

Plenty of artifact removal tools and pipelines have been developed to correct the resting EEG waves and discover scientific values behind. Without expertised visual inspection, it is susceptible to derive improper preprocessing, resulting in either insufficient preprocessed EEG (IPE) or excessive preprocessed EEG (EPE). However, little is known about the impacts of IPE or EPE on postprocessing in the temporal, frequency, and spatial domains, particularly as to the spectra and the functional connectivity analysis. Here, the clean EEG (CE) with linear and quasi-stationary assumption was synthesized as ground truth based on the New-York head model and the multivariate autoregressive model. Later, IPE and EPE were simulated by injecting Gaussian noise and losing brain components, respectively. Spectral homogeneities of all EEGs were evaluated by the proposed Parallel LOg Spectra index (PaLOSi). Then, the impacts on postprocessing were quantified by the IPE/EPE deviation from CE as to the temporal statistics, multichannel power, cross spectra, scalp EEG network properties, and source dispersion. Lastly, the association between PaLOSi and varying trends of postprocessing outcomes was analyzed with evolutionary preprocessing states. We found that compared with CE: 1) IPE (EPE) temporal statistics deviated more greatly with more noise injected (brain activities discarded); 2) IPE (EPE) power was higher (lower), and IPE power was almost parallel to that of CE across frequencies, while EPE power deviation decreased with higher frequencies; IPE cross spectra deviated more greatly than EPE, except for β band; 3) derived from 7 coupling measures, IPE (EPE) network had lower (higher) transmission efficiency and worse (better) integration ability; 4) IPE sources distributed more dispersedly with greater strength while EPE sources activated more focally with lower amplitudes; 5) PaLOSi was consistently correlated with varying trends of investigated postprocessing for both simulated and real data. This study shed light on how the postprocessing outcomes are affected by the preprocessing states and PaLOSi is a promising quality control metric for creating normative EEG databases.

Frequency-band specific directed connectivity networks reveal functional disruptions and pathogenic patterns in temporal lobe epilepsy: a MEG study

This study investigates the network mechanisms of temporal lobe epilepsy (TLE) using MEG data, focusing on directed connectivity networks across different frequency bands. Unlike previous studies that primarily localize epileptogenic zones, this research aims to explore whole-brain network differences between left TLE (lTLE), right TLE (rTLE), and healthy controls (HCs). MEG data from 13 lTLE patients, 21 rTLE patients, and 14 HCs were source-reconstructed to 116 brain regions (AAL116). Directed Transfer Function (DTF) was used to construct directed connectivity networks, followed by networks and graph-theoretical analyses. The results indicate that, compared to HCs, TLE subjects exhibited a significant increase in average connectivity strength in the Low Gamma band. The connectivity patterns across frequency bands in TLE patients were found to be unstable. Both HC and TLE subjects demonstrated left hemisphere lateralization. In the mid-to-low frequency bands, TLE subjects showed increases in global clustering coefficient (GCC), global characteristic path length (GCPL), and local efficiency (LE) compared to HCs, which is attributed to enhanced synchronization between local brain regions in TLE subjects.

Functional connectotomy of a whole-brain model reveals tumor-induced alterations to neuronal dynamics in glioma patients

Brain tumors can induce pathological changes in neuronal dynamics that are reflected in functional connectivity measures. Here, we use a whole-brain modeling approach to investigate pathological alterations to neuronal activity in glioma patients. By fitting a Hopf whole-brain model to empirical functional connectivity, we investigate glioma-induced changes in optimal model parameters. We observe considerable differences in neuronal dynamics between glioma patients and healthy controls, both on an individual and population-based level. In particular, model parameter estimation suggests that local tumor pathology causes changes in brain dynamics by increasing the influence of interregional interactions on global neuronal activity. Our approach demonstrates that whole-brain models provide valuable insights for understanding glioma-associated alterations in functional connectivity.

The impact of channel density, inverse solutions, connectivity metrics and calibration errors on OPM-MEG connectivity analysis: A simulation study

Magnetoencephalography (MEG) systems based on optically pumped magnetometers (OPMs) have rapidly developed in the fields of brain function, health, and disease. Functional connectivity analysis related to the resting-state has gained popularity as a field of research in recent years. Several studies have attempted to use OPM-based MEG (OPM-MEG) for brain network estimation research; however, the choice of source connectivity analysis pipeline may lead to outcome variability. Several methods and related parameters must be selected carefully at each step of the analysis. Therefore, this study assessed the effect of such analytical variability on the OPM-MEG connectivity analysis by conducting simulations. Synthetic MEG data corresponding to two default mode networks (DMN) with six or ten DMN regions were generated using the Gaussian Graphical Spectral (GGS) model. Six intersensor spacings were constructed, and six inverse algorithms and six functional connectivity measures were selected to assess their impact on the network reconstruction accuracy. Three potential sources of error - errors in the sensor gain, crosstalk, and angular errors of the sensitive axis of the OPM - were also assessed. Analytical variability with regard to the tested intersensor spacings, inverse solutions, and functional connectivity measures led to high result variability. Crosstalk exerted a significant impact on the accuracy, which may lead to network reconstruction failure. The accuracy improvement caused by an increase in the sensor density may be reduced by gain and angular errors. The minimum norm estimate (MNE) and weighted minimum norm estimate (wMNE) exhibited low robustness to sensor noise and calibration errors. Hence, a calibration workflow for accurate sensor parameters, such as the gain and direction of the sensitive axis, before commencing OPM-MEG measurement and a careful choice of different method combinations play crucial roles in ensuring that OPMs yield optimal results for functional connectivity analysis. A thorough framework for analyzing brain connectivity networks was provided herein.

Rhythm Facilitates Auditory Working Memory via Beta-Band Encoding and Theta-Band Maintenance

Rhythm, as a prominent characteristic of auditory experiences such as speech and music, is known to facilitate attention, yet its contribution to working memory (WM) remains unclear. Here, human participants temporarily retained a 12-tone sequence presented rhythmically or arrhythmically in WM and performed a pitch change-detection task. Behaviorally, while having comparable accuracy, rhythmic tone sequences showed a faster response time and lower response boundaries in decision-making. Electroencephalographic recordings revealed that rhythmic sequences elicited enhanced non-phase-locked beta-band (16 Hz-33 Hz) and theta-band (3 Hz-5 Hz) neural oscillations during sensory encoding and WM retention periods, respectively. Importantly, the two-stage neural signatures were correlated with each other and contributed to behavior. As beta-band and theta-band oscillations denote the engagement of motor systems and WM maintenance, respectively, our findings imply that rhythm facilitates auditory WM through intricate oscillation-based interactions between the motor and auditory systems that facilitate predictive attention to auditory sequences.

Functional network disruption in cognitively unimpaired autosomal dominant Alzheimer's disease: a magnetoencephalography study

Understanding the nature and onset of neurophysiological changes, and the selective vulnerability of central hub regions in the functional network, may aid in managing the growing impact of Alzheimer's disease on society. However, the precise neurophysiological alterations occurring in the pre-clinical stage of human Alzheimer's disease remain controversial. This study aims to provide increased insights on quantitative neurophysiological alterations during a true early stage of Alzheimer's disease. Using high spatial resolution source-reconstructed magnetoencephalography, we investigated regional and whole-brain neurophysiological changes in a unique cohort of 11 cognitively unimpaired individuals with pathogenic mutations in the presenilin-1 or amyloid precursor protein gene and a 1:3 matched control group (n = 33) with a median age of 49 years. We examined several quantitative magnetoencephalography measures that have been shown robust in detecting differences in sporadic Alzheimer's disease patients and are sensitive to excitation-inhibition imbalance. This includes spectral power and functional connectivity in different frequency bands. We also investigated hub vulnerability using the hub disruption index. To understand how magnetoencephalography measures change as the disease progresses through its pre-clinical stage, correlations between magnetoencephalography outcomes and various clinical variables like age were analysed. A comparison of spectral power between mutation carriers and controls revealed oscillatory slowing, characterized by widespread higher theta (4-8 Hz) power, a lower posterior peak frequency and lower occipital alpha 2 (10-13 Hz) power. Functional connectivity analyses presented a lower whole-brain (amplitude-based) functional connectivity in the alpha (8-13 Hz) and beta (13-30 Hz) bands, predominantly located in parieto-temporal hub regions. Furthermore, we found a significant hub disruption index for (phase-based) functional connectivity in the theta band, attributed to both higher functional connectivity in 'non-hub' regions alongside a hub disruption. Neurophysiological changes did not correlate with indicators of pre-clinical disease progression in mutation carriers after multiple comparisons correction. Our findings provide evidence that oscillatory slowing and functional connectivity differences occur before cognitive impairment in individuals with autosomal dominant mutations leading to early onset Alzheimer's disease. The nature and direction of these alterations are comparable to those observed in the clinical stages of Alzheimer's disease, suggest an early excitation-inhibition imbalance, and fit with the activity-dependent functional degeneration hypothesis. These insights may prove useful for early diagnosis and intervention in the future.

Filtered Point Processes Tractably Capture Rhythmic And Broadband Power Spectral Structure in Neural Electrophysiological Recordings

Neural electrophysiological recordings arise from interacting rhythmic (oscillatory) and broadband (aperiodic) biological subprocesses. Both rhythmic and broadband processes contribute to the neural power spectrum, which decomposes the variance of a neural recording across frequencies. Although an extensive body of literature has successfully studied rhythms in various diseases and brain states, researchers only recently have systematically studied the characteristics of broadband effects in the power spectrum. Broadband effects can generally be categorized as 1) shifts in power across all frequencies, which correlate with changes in local firing rates and 2) changes in the overall shape of the power spectrum, such as the spectral slope or power law exponent. Shape changes are evident in various conditions and brain states, influenced by factors such as excitation to inhibition balance, age, and various diseases. It is increasingly recognized that broadband and rhythmic effects can interact on a sub-second timescale. For example, broadband power is time-locked to the phase of <1 Hz rhythms in propofol induced unconsciousness. Modeling tools that explicitly deal with both rhythmic and broadband contributors to the power spectrum and that capture their interactions are essential to help improve the interpretability of power spectral effects. Here, we introduce a tractable stochastic forward modeling framework designed to capture both narrowband and broadband spectral effects when prior knowledge or theory about the primary biophysical processes involved is available. Population-level neural recordings are modeled as the sum of filtered point processes (FPPs), each representing the contribution of a different biophysical process such as action potentials or postsynaptic potentials of different types. Our approach builds on prior neuroscience FPP work by allowing multiple interacting processes and time-varying firing rates and by deriving theoretical power spectra and cross-spectra. We demonstrate several properties of the models, including that they divide the power spectrum into frequency ranges dominated by rhythmic and broadband effects, and that they can capture spectral effects across multiple timescales, including sub-second cross-frequency coupling. The framework can be used to interpret empirically observed power spectra and cross-frequency coupling effects in biophysical terms, which bridges the gap between theoretical models and experimental results.

Corticocortical and corticomuscular connectivity dynamics in standing posture: electroencephalography study

Cortical mechanism is necessary for human standing control. Previous research has demonstrated that cortical oscillations and corticospinal excitability respond flexibly to postural demands. However, it is unclear how corticocortical and corticomuscular connectivity changes dynamically during standing with spontaneous postural sway and over time. This study investigated the dynamics of sway- and time-varying connectivity using electroencephalography and electromyography. Electroencephalography and electromyography were recorded in sitting position and 3 standing postures with varying base-of-support: normal standing, one-leg standing, and standing on a piece of wood. For sway-varying connectivity, corticomuscular connectivity was calculated based on the timing of peak velocity in anteroposterior sway. For time-varying connectivity, corticocortical connectivity was measured using the sliding-window approach. This study found that corticomuscular connectivity was strengthened at the peak velocity of postural sway in the γ- and β-frequency bands. For time-varying corticocortical connectivity, the θ-connectivity in all time-epoch was classified into 7 clusters including posture-relevant component. In one of the 7 clusters, strong connectivity pairs were concentrated in the mid-central region, and the proportion of epochs under narrow-base standing conditions was significantly higher, indicating a functional role for posture balance. These findings shed light on the connectivity dynamics and cortical oscillation that govern standing balance.

EEG resting state alpha dynamics predict an individual's vulnerability to auditory hallucinations

Task-free brain activity exhibits spontaneous fluctuations between functional states, characterized by synchronized activation patterns in distributed resting-state (RS) brain networks. The temporal dynamics of the networks' electrophysiological signatures reflect individual variations in brain activity and connectivity linked to mental states and cognitive functions and can predict or monitor vulnerability to develop psychiatric or neurological disorders. In particular, RS alpha fluctuations modulate perceptual sensitivity, attentional shifts, and cognitive control, and could therefore reflect a neural correlate of increased vulnerability to sensory distortions, including the proneness to hallucinatory experiences. We recorded 5 min of RS EEG from 33 non-clinical individuals varying in hallucination proneness (HP) to investigate links between task-free alpha dynamics and vulnerability to hallucinations. To this end, we used a dynamic brain state allocation method to identify five recurrent alpha states together with their spatiotemporal dynamics and most active brain areas through source reconstruction. The dynamical features of a state marked by activation in somatosensory, auditory, and posterior default-mode network areas predicted auditory and auditory-verbal HP, but not general HP, such that individuals with higher vulnerability to auditory hallucinations spent more time in this state. The temporal dynamics of spontaneous alpha activity might reflect individual differences in attention to internally generated sensory events and altered auditory perceptual sensitivity. Altered RS alpha dynamics could therefore instantiate a neural marker of increased vulnerability to auditory hallucinations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11571-024-10093-1.

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.

Functionally annotated electrophysiological neuromarkers of healthy ageing and memory function

The unprecedented increase in life expectancy presents a unique opportunity and the necessity to explore both healthy and pathological aspects of ageing. Electroencephalography (EEG) has been widely used to identify neuromarkers of cognitive ageing due to its affordability and richness in information. However, despite the growing volume of data and methodological advancements, the abundance of contradictory and non-reproducible findings has hindered clinical translation. To address these challenges, our study introduces a comprehensive workflow expanding on previous EEG studies and investigates various static and dynamic power and connectivity estimates as potential neuromarkers of cognitive ageing in a large dataset. We also assess the robustness of our findings by testing their susceptibility to band specification. Finally, we characterise our findings using functionally annotated brain networks to improve their interpretability and multi-modal integration. Our analysis demonstrates the effect of methodological choices on findings and that dynamic rather than static neuromarkers are not only more sensitive but also more robust. Consequently, they emerge as strong candidates for cognitive ageing neuromarkers. Moreover, we were able to replicate the most established EEG findings in cognitive ageing, such as alpha oscillation slowing, increased beta power, reduced reactivity across multiple bands, and decreased delta connectivity. Additionally, when considering individual variations in the alpha band, we clarified that alpha power is characteristic of memory performance rather than ageing, highlighting its potential as a neuromarker for cognitive ageing. Finally, our approach using functionally annotated source reconstruction allowed us to provide insights into domain-specific electrophysiological mechanisms underlying memory performance and ageing. HIGHLIGHTS: We provide an open and reproducible pipeline with a comprehensive workflow to investigate static and dynamic EEG neuromarkers. Neuromarkers related to neural dynamics are sensitive and robust. Individualised alpha power characterises cognitive performance rather than ageing. Functional annotation allows cross-modal interpretation of EEG findings.

Neurophysiological brain function predicts response to cognitive rehabilitation and mindfulness in multiple sclerosis: a randomized trial

BACKGROUND

Cognitive treatment response varies highly in people with multiple sclerosis (PwMS). Identification of mechanisms is essential for predicting response.

OBJECTIVES

This study aimed to investigate whether brain network function predicts response to cognitive rehabilitation therapy (CRT) and mindfulness-based cognitive therapy (MBCT).

METHODS

PwMS with cognitive complaints completed CRT, MBCT, or enhanced treatment as usual (ETAU) and performed three measurements (baseline, post-treatment, 6-month follow-up). Baseline magnetoencephalography (MEG) measures were used to predict treatment effects on cognitive complaints, personalized cognitive goals, and information processing speed (IPS) using mixed models (secondary analysis REMIND-MS study).

RESULTS

We included 105 PwMS (96 included in prediction analyses; 32 CRT, 31 MBCT, 33 ETAU), and 56 healthy controls with baseline MEG. MEG did not predict reductions in complaints. Higher connectivity predicted better goal achievement after MBCT (p = 0.010) and CRT (p = 0.018). Lower gamma power (p = 0.006) and higher connectivity (p = 0.020) predicted larger IPS benefits after MBCT. These MEG predictors indicated worse brain function compared to healthy controls (p < 0.05).

CONCLUSIONS

Brain network function predicted better cognitive goal achievement after MBCT and CRT, and IPS improvements after MBCT. PwMS with neuronal slowing and hyperconnectivity were most prone to show treatment response, making network function a promising tool for personalized treatment recommendations.

TRIAL REGISTRATION

The REMIND-MS study was prospectively registered in the Dutch Trial registry (NL6285; https://trialsearch.who.int/Trial2.aspx?TrialID=NTR6459 ).

Data-driven MEG analysis to extract fMRI resting-state networks

The electrophysiological basis of resting-state networks (RSN) is still under debate. In particular, no principled mechanism has been determined that is capable of explaining all RSN equally well. While magnetoencephalography (MEG) and electroencephalography are the methods of choice to determine the electrophysiological basis of RSN, no standard analysis pipeline of RSN yet exists. In this article, we compare the two main existing data-driven analysis strategies for extracting RSNs from MEG data and introduce a third approach. The first approach uses phase-amplitude coupling to determine the RSN. The second approach extracts RSN through an independent component analysis of the Hilbert envelope in different frequency bands, while the third new approach uses a singular value decomposition instead. To evaluate these approaches, we compare the MEG-RSN to the functional magnetic resonance imaging (fMRI)-RSN from the same subjects. Overall, it was possible to extract RSN with MEG using all three techniques, which matched the group-specific fMRI-RSN. Interestingly the new approach based on SVD yielded significantly higher correspondence to five out of seven fMRI-RSN than the two existing approaches. Importantly, with this approach, all networks-except for the visual network-had the highest correspondence to the fMRI networks within one frequency band. Thereby we provide further insights into the electrophysiological underpinnings of the fMRI-RSNs. This knowledge will be important for the analysis of the electrophysiological connectome.

Deep brain stimulation of the central thalamus restores arousal and motivation in a zolpidem-responsive patient with akinetic mutism after severe brain injury

After severe brain injury, zolpidem is known to cause spectacular, often short-lived, restorations of brain functions in a small subgroup of patients. Previously, we showed that these zolpidem-induced neurological recoveries can be paralleled by significant changes in functional connectivity throughout the brain. Deep brain stimulation (DBS) is a neurosurgical intervention known to modulate functional connectivity in a wide variety of neurological disorders. In this study, we used DBS to restore arousal and motivation in a zolpidem-responsive patient with severe brain injury and a concomitant disorder of diminished motivation, more than 10 years after surviving hypoxic ischemia. We found that DBS of the central thalamus, targeted at the centromedian-parafascicular complex, immediately restored arousal and was able to transition the patient from a state of deep sleep to full wakefulness. Moreover, DBS was associated with temporary restoration of communication and ability to walk and eat in an otherwise wheelchair-bound and mute patient. With the use of magnetoencephalography (MEG), we revealed that DBS was generally associated with a marked decrease in aberrantly high levels of functional connectivity throughout the brain, mimicking the effects of zolpidem. These results imply that 'pathological hyperconnectivity' after severe brain injury can be associated with reduced arousal and behavioral performance and that DBS is able to modulate connectivity towards a 'healthier baseline' with lower synchronization, and, can restore functional brain networks long after severe brain injury. The presence of hyperconnectivity after brain injury may be a possible future marker for a patient's responsiveness for restorative interventions, such as DBS, and suggests that lower degrees of overall brain synchronization may be conducive to cognition and behavioral responsiveness.

The relationship between pathological brain activity and functional network connectivity in glioma patients

PURPOSE

Glioma is associated with pathologically high (peri)tumoral brain activity, which relates to faster progression. Functional connectivity is disturbed locally and throughout the entire brain, associating with symptomatology. We, therefore, investigated how local activity and network measures relate to better understand how the intricate relationship between the tumor and the rest of the brain may impact disease and symptom progression.

METHODS

We obtained magnetoencephalography in 84 de novo glioma patients and 61 matched healthy controls. The offset of the power spectrum, a proxy of neuronal activity, was calculated for 210 cortical regions. We calculated patients' regional deviations in delta, theta and lower alpha network connectivity as compared to controls, using two network measures: clustering coefficient (local connectivity) and eigenvector centrality (integrative connectivity). We then tested group differences in activity and connectivity between (peri)tumoral, contralateral homologue regions, and the rest of the brain. We also correlated regional offset to connectivity.

RESULTS

As expected, patients' (peri)tumoral activity was pathologically high, and patients showed higher clustering and lower centrality than controls. At the group-level, regionally high activity related to high clustering in controls and patients alike. However, within-patient analyses revealed negative associations between regional deviations in brain activity and clustering, such that pathologically high activity coincided with low network clustering, while regions with 'normal' activity levels showed high network clustering.

CONCLUSION

Our results indicate that pathological activity and connectivity co-localize in a complex manner in glioma. This insight is relevant to our understanding of disease progression and cognitive symptomatology.

Functional brain alterations in symptomatic dermographism patients-An exploratory magnetoencephalography study

Symptomatic dermographism (SD) is a common form of urticaria, which is triggered by stroking the skin. Brain involvement in its aetiology was investigated by means of magnetoencephalography (MEG) after provocation with histamine and dermography. Wheals were induced by histamine skin prick test and dermography in twelve SD patients and fourteen controls. Itch severity was scored on a Visual Analogue Scale (VAS). Relative power and functional connectivity (FC) were measured using a 306-channel whole-head MEG system at baseline and 10 min after histamine and dermography, and contrasted between groups and conditions. Furthermore, wheal diameter and itch scores after these procedures were correlated with the MEG values. SD patients had higher itch scores after histamine and dermography. No significant group-differences were observed in relative power or FC for any condition. In both groups, power decreases were mostly observed in the beta band, and power increases in the alpha bands, after provocation, with more regions involved in patients compared to controls. Increased FC was seen after histamine in patients, and after dermography in controls. In patients only, dermography and histamine wheal size correlated with the alpha2 power in the regions of interest that showed significant condition effects after these procedures. Our findings may be cautiously interpreted as aberrant itch processing, and suggest involvement of the central nervous system in the aetiology of SD.

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.

Cortical Synchrony and Information Flow during Transition from Wakefulness to Light Non-Rapid Eye Movement Sleep

Sleep is a highly stereotyped phenomenon, requiring robust spatiotemporal coordination of neural activity. Understanding how the brain coordinates neural activity with sleep onset can provide insights into the physiological functions subserved by sleep and the pathologic phenomena associated with sleep onset. We quantified whole-brain network changes in synchrony and information flow during the transition from wakefulness to light non-rapid eye movement (NREM) sleep, using MEG imaging in a convenient sample of 14 healthy human participants (11 female; mean 63.4 years [SD 11.8 years]). We furthermore performed computational modeling to infer excitatory and inhibitory properties of local neural activity. The transition from wakefulness to light NREM was identified to be encoded in spatially and temporally specific patterns of long-range synchrony. Within the delta band, there was a global increase in connectivity from wakefulness to light NREM, which was highest in frontoparietal regions. Within the theta band, there was an increase in connectivity in fronto-parieto-occipital regions and a decrease in temporal regions from wakefulness to Stage 1 sleep. Patterns of information flow revealed that mesial frontal regions receive hierarchically organized inputs from broad cortical regions upon sleep onset, including direct inflow from occipital regions and indirect inflow via parieto-temporal regions within the delta frequency band. Finally, biophysical neural mass modeling demonstrated changes in the anterior-to-posterior distribution of cortical excitation-to-inhibition with increased excitation-to-inhibition model parameters in anterior regions in light NREM compared with wakefulness. Together, these findings uncover whole-brain corticocortical structure and the orchestration of local and long-range, frequency-specific cortical interactions in the sleep-wake transition.SIGNIFICANCE STATEMENT Our work uncovers spatiotemporal cortical structure of neural synchrony and information flow upon the transition from wakefulness to light non-rapid eye movement sleep. Mesial frontal regions were identified to receive hierarchically organized inputs from broad cortical regions, including both direct inputs from occipital regions and indirect inputs via the parieto-temporal regions within the delta frequency range. Biophysical neural mass modeling revealed a spatially heterogeneous, anterior-posterior distribution of cortical excitation-to-inhibition. Our findings shed light on the orchestration of local and long-range cortical neural structure that is fundamental to sleep onset, and support an emerging view of cortically driven regulation of sleep homeostasis.

The impact of ROI extraction method for MEG connectivity estimation: practical recommendations for the study of resting state data

Magnetoencephalography and electroencephalography (M/EEG) seed-based connectivity analysis requires the extraction of measures from regions of interest (ROI). M/EEG ROI-derived source activity can be treated in different ways. It is possible, for instance, to average each ROI's time series prior to calculating connectivity measures. Alternatively, one can compute connectivity maps for each element of the ROI prior to dimensionality reduction to obtain a single map. The impact of these different strategies on connectivity results is still unclear. Here, we address this question within a large MEG resting state cohort (N=113) and within simulated data. We consider 68 ROIs (Desikan-Kiliany atlas), two measures of connectivity (phase locking value-PLV, and its imaginary counterpart- ciPLV), and three frequency bands (theta 4-8 Hz, alpha 9-12 Hz, beta 15-30 Hz). We compare four extraction methods: (i) mean, or (ii) PCA of the activity within the seed or ROI before computing connectivity, map of the (iii) average, or (iv) maximum connectivity after computing connectivity for each element of the seed. Hierarchical clustering is then applied to compare connectivity outputs across multiple strategies, followed by direct contrasts across extraction methods. Finally, the results are validated by using a set of realistic simulations. We show that ROI-based connectivity maps vary remarkably across strategies in terms of connectivity magnitude and spatial distribution. Dimensionality reduction procedures conducted after computing connectivity are more similar to each-other, while PCA before approach is the most dissimilar to other approaches. Although differences across methods are consistent across frequency bands, they are influenced by the connectivity metric and ROI size. Greater differences were observed for ciPLV than PLV, and in larger ROIs. Realistic simulations confirmed that after aggregation procedures are generally more accurate but have lower specificity (higher rate of false positive connections). Though computationally demanding, after dimensionality reduction strategies should be preferred when higher sensitivity is desired. Given the remarkable differences across aggregation procedures, caution is warranted in comparing results across studies applying different methods.

Connectivity-based Meta-Bands: A new approach for automatic frequency band identification in connectivity analyses

The majority of electroencephalographic (EEG) and magnetoencephalographic (MEG) studies filter and analyse neural signals in specific frequency ranges, known as "canonical" frequency bands. However, this segmentation, is not exempt from limitations, mainly due to the lack of adaptation to the neural idiosyncrasies of each individual. In this study, we introduce a new data-driven method to automatically identify frequency ranges based on the topological similarity of the frequency-dependent functional neural network. The resting-state neural activity of 195 cognitively healthy subjects from three different databases (MEG: 123 subjects; EEG1: 27 subjects; EEG2: 45 subjects) was analysed. In a first step, MEG and EEG signals were filtered with a narrow-band filter bank (1 Hz bandwidth) from 1 to 70 Hz with a 0.5 Hz step. Next, the connectivity in each of these filtered signals was estimated using the orthogonalized version of the amplitude envelope correlation to obtain the frequency-dependent functional neural network. Finally, a community detection algorithm was used to identify communities in the frequency domain showing a similar network topology. We have called this approach the "Connectivity-based Meta-Bands" (CMB) algorithm. Additionally, two types of synthetic signals were used to configure the hyper-parameters of the CMB algorithm. We observed that the classical approaches to band segmentation are partially aligned with the underlying network topologies at group level for the MEG signals, but they are missing individual idiosyncrasies that may be biasing previous studies, as revealed by our methodology. On the other hand, the sensitivity of EEG signals to reflect this underlying frequency-dependent network structure is limited, revealing a simpler frequency parcellation, not aligned with that defined by the "canonical" frequency bands. To the best of our knowledge, this is the first study that proposes an unsupervised band segmentation method based on the topological similarity of functional neural network across frequencies. This methodology fully accounts for subject-specific patterns, providing more robust and personalized analyses, and paving the way for new studies focused on exploring the frequency-dependent structure of brain connectivity.

Tau protein spreads through functionally connected neurons in Alzheimer's disease: a combined MEG/PET study

Recent studies on Alzheimer's disease (AD) suggest that tau proteins spread through the brain following neuronal connections. Several mechanisms could be involved in this process: spreading between brain regions that interact strongly (functional connectivity); through the pattern of anatomical connections (structural connectivity); or simple diffusion. Using magnetoencephalography (MEG), we investigated which spreading pathways influence tau protein spreading by modelling the tau propagation process using an epidemic spreading model. We compared the modelled tau depositions with 18F-flortaucipir PET binding potentials at several stages of the AD continuum. In this cross-sectional study, we analysed source-reconstructed MEG data and dynamic 100-min 18F-flortaucipir PET from 57 subjects positive for amyloid-β pathology [preclinical AD (n = 16), mild cognitive impairment (MCI) due to AD (n = 16) and AD dementia (n = 25)]. Cognitively healthy subjects without amyloid-β pathology were included as controls (n = 25). Tau propagation was modelled as an epidemic process (susceptible-infected model) on MEG-based functional networks [in alpha (8-13 Hz) and beta (13-30 Hz) bands], a structural or diffusion network, starting from the middle and inferior temporal lobe. The group-level network of the control group was used as input for the model to predict tau deposition in three stages of the AD continuum. To assess performance, model output was compared to the group-specific tau deposition patterns as measured with 18F-flortaucipir PET. We repeated the analysis by using networks of the preceding disease stage and/or using regions with most observed tau deposition during the preceding stage as seeds. In the preclinical AD stage, the functional networks predicted most of the modelled tau-PET binding potential, with best correlations between model and tau-PET [corrected amplitude envelope correlation (AEC-c) alpha C = 0.584; AEC-c beta C = 0.569], followed by the structural network (C = 0.451) and simple diffusion (C = 0.451). Prediction accuracy declined for the MCI and AD dementia stages, although the correlation between modelled tau and tau-PET binding remained highest for the functional networks (C = 0.384; C = 0.376). Replacing the control-network with the network from the preceding disease stage and/or alternative seeds improved prediction accuracy in MCI but not in the dementia stage. These results suggest that in addition to structural connections, functional connections play an important role in tau spread, and highlight that neuronal dynamics play a key role in promoting this pathological process. Aberrant neuronal communication patterns should be taken into account when identifying targets for future therapy. Our results also suggest that this process is more important in earlier disease stages (preclinical AD/MCI); possibly, in later stages, other processes may be influential.

Comparing individual and group-level simulated neurophysiological brain connectivity using the Jansen and Rit neural mass model

Computational models are often used to assess how functional connectivity (FC) patterns emerge from neuronal population dynamics and anatomical brain connections. It remains unclear whether the commonly used group-averaged data can predict individual FC patterns. The Jansen and Rit neural mass model was employed, where masses were coupled using individual structural connectivity (SC). Simulated FC was correlated to individual magnetoencephalography-derived empirical FC. FC was estimated using phase-based (phase lag index (PLI), phase locking value (PLV)), and amplitude-based (amplitude envelope correlation (AEC)) metrics to analyze their goodness of fit for individual predictions. Individual FC predictions were compared against group-averaged FC predictions, and we tested whether SC of a different participant could equally well predict participants' FC patterns. The AEC provided a better match between individually simulated and empirical FC than phase-based metrics. Correlations between simulated and empirical FC were higher using individual SC compared to group-averaged SC. Using SC from other participants resulted in similar correlations between simulated and empirical FC compared to using participants' own SC. This work underlines the added value of FC simulations using individual instead of group-averaged SC for this particular computational model and could aid in a better understanding of mechanisms underlying individual functional network trajectories.

Neurophysiological alterations in mice and humans carrying mutations in APP and PSEN1 genes

BACKGROUND

Studies in animal models of Alzheimer's disease (AD) have provided valuable insights into the molecular and cellular processes underlying neuronal network dysfunction. Whether and how AD-related neurophysiological alterations translate between mice and humans remains however uncertain.

METHODS

We characterized neurophysiological alterations in mice and humans carrying AD mutations in the APP and/or PSEN1 genes, focusing on early pre-symptomatic changes. Longitudinal local field potential recordings were performed in APP/PS1 mice and cross-sectional magnetoencephalography recordings in human APP and/or PSEN1 mutation carriers. All recordings were acquired in the left frontal cortex, parietal cortex, and hippocampus. Spectral power and functional connectivity were analyzed and compared with wildtype control mice and healthy age-matched human subjects.

RESULTS

APP/PS1 mice showed increased absolute power, especially at higher frequencies (beta and gamma) and predominantly between 3 and 6 moa. Relative power showed an overall shift from lower to higher frequencies over almost the entire recording period and across all three brain regions. Human mutation carriers, on the other hand, did not show changes in power except for an increase in relative theta power in the hippocampus. Mouse parietal cortex and hippocampal power spectra showed a characteristic peak at around 8 Hz which was not significantly altered in transgenic mice. Human power spectra showed a characteristic peak at around 9 Hz, the frequency of which was significantly reduced in mutation carriers. Significant alterations in functional connectivity were detected in theta, alpha, beta, and gamma frequency bands, but the exact frequency range and direction of change differed for APP/PS1 mice and human mutation carriers.

CONCLUSIONS

Both mice and humans carrying APP and/or PSEN1 mutations show abnormal neurophysiological activity, but several measures do not translate one-to-one between species. Alterations in absolute and relative power in mice should be interpreted with care and may be due to overexpression of amyloid in combination with the absence of tau pathology and cholinergic degeneration. Future studies should explore whether changes in brain activity in other AD mouse models, for instance, those also including tau pathology, provide better translation to the human AD continuum.

Top-down reconfiguration of SMA cortical connectivity during action preparation

The Bereitschaftspotential (BP), a scalp potential recorded in humans during action preparation, is characterized by a slow amplitude increase over fronto-central regions as action execution approaches. We recorded TMS evoked-potentials (TEP) stimulating the supplementary motor area (SMA) at different time-points during a Go/No-Go task to assess whether and how cortical excitability and connectivity of this region change as the BP increases. When approaching BP peak, left SMA reactivity resulted greater. Concurrently, its effective connectivity increased with the left occipital areas, while it decreased with the right inferior frontal gyrus, indicating a fast reconfiguration of cortical networks during the preparation of the forthcoming action. Functional connectivity patterns supported these findings, suggesting a critical role of frequency-specific inter-areal interactions in implementing top-down mechanisms in the sensorimotor system prior to action. These findings reveal that BP time-course reflects quantitative and qualitative changes in SMA communication patterns that shape mechanisms involved in motor readiness.

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.

The longitudinal relation between executive functioning and multilayer network topology in glioma patients

Many patients with glioma, primary brain tumors, suffer from poorly understood executive functioning deficits before and/or after tumor resection. We aimed to test whether frontoparietal network centrality of multilayer networks, allowing for integration across multiple frequencies, relates to and predicts executive functioning in glioma. Patients with glioma (n = 37) underwent resting-state magnetoencephalography and neuropsychological tests assessing word fluency, inhibition, and set shifting before (T1) and one year after tumor resection (T2). We constructed binary multilayer networks comprising six layers, with each layer representing frequency-specific functional connectivity between source-localized time series of 78 cortical regions. Average frontoparietal network multilayer eigenvector centrality, a measure for network integration, was calculated at both time points. Regression analyses were used to investigate associations with executive functioning. At T1, lower multilayer integration (p = 0.017) and epilepsy (p = 0.006) associated with poorer set shifting (adj. R2 = 0.269). Decreasing multilayer integration (p = 0.022) and not undergoing chemotherapy at T2 (p = 0.004) related to deteriorating set shifting over time (adj. R2 = 0.283). No significant associations were found for word fluency or inhibition, nor did T1 multilayer integration predict changes in executive functioning. As expected, our results establish multilayer integration of the frontoparietal network as a cross-sectional and longitudinal correlate of executive functioning in glioma patients. However, multilayer integration did not predict postoperative changes in executive functioning, which together with the fact that this correlate is also found in health and other diseases, limits its specific clinical relevance in glioma.

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.

Uncorrelated bilateral cortical input becomes timed across hippocampal subfields for long waves whereas gamma waves are largely ipsilateral

The role of interhemispheric connections along successive segments of cortico-hippocampal circuits is poorly understood. We aimed to obtain a global picture of spontaneous transfer of activity during non-theta states across several nodes of the bilateral circuit in anesthetized rats. Spatial discrimination techniques applied to bilateral laminar field potentials (FP) across the CA1/Dentate Gyrus provided simultaneous left and right readouts in five FP generators that reflect activity in specific hippocampal afferents and associative pathways. We used a battery of correlation and coherence analyses to extract complementary aspects at different time scales and frequency bands. FP generators exhibited varying bilateral correlation that was high in CA1 and low in the Dentate Gyrus. The submillisecond delays indicate coordination but not support for synaptic dependence of one side on another. The time and frequency characteristics of bilateral coupling were specific to each generator. The Schaffer generator was strongly bilaterally coherent for both sharp waves and gamma waves, although the latter maintained poor amplitude co-variation. The lacunosum-moleculare generator was composed of up to three spatially overlapping activities, and globally maintained high bilateral coherence for long but not short (gamma) waves. These two CA1 generators showed no ipsilateral relationship in any frequency band. In the Dentate Gyrus, strong bilateral coherence was observed only for input from the medial entorhinal areas, while those from the lateral entorhinal areas were largely asymmetric, for both alpha and gamma waves. Granger causality testing showed strong bidirectional relationships between all homonymous bilateral generators except the lateral entorhinal input and a local generator in the Dentate Gyrus. It also revealed few significant relationships between ipsilateral generators, most notably the anticipation of lateral entorhinal cortex toward all others. Thus, with the notable exception of the lateral entorhinal areas, there is a marked interhemispheric coherence primarily for slow envelopes of activity, but not for pulse-like gamma waves, except in the Schafer segment. The results are consistent with essentially different streams of activity entering from and returning to the cortex on each side, with slow waves reflecting times of increased activity exchange between hemispheres and fast waves generally reflecting ipsilateral processing.

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.

Network Hyperexcitability in Early Alzheimer's Disease: Is Functional Connectivity a Potential Biomarker?

Network hyperexcitability (NH) is an important feature of the pathophysiology of Alzheimer's disease. Functional connectivity (FC) of brain networks has been proposed as a potential biomarker for NH. Here we use a whole brain computational model and resting-state MEG recordings to investigate the relation between hyperexcitability and FC. Oscillatory brain activity was simulated with a Stuart Landau model on a network of 78 interconnected brain regions. FC was quantified with amplitude envelope correlation (AEC) and phase coherence (PC). MEG was recorded in 18 subjects with subjective cognitive decline (SCD) and 18 subjects with mild cognitive impairment (MCI). Functional connectivity was determined with the corrected AECc and phase lag index (PLI), in the 4-8 Hz and the 8-13 Hz bands. The excitation/inhibition balance in the model had a strong effect on both AEC and PC. This effect was different for AEC and PC, and was influenced by structural coupling strength and frequency band. Empirical FC matrices of SCD and MCI showed a good correlation with model FC for AEC, but less so for PC. For AEC the fit was best in the hyperexcitable range. We conclude that FC is sensitive to changes in E/I balance. The AEC was more sensitive than the PLI, and results were better for the thetaband than the alpha band. This conclusion was supported by fitting the model to empirical data. Our study justifies the use of functional connectivity measures as surrogate markers for E/I balance.

Magnetoencephalography to measure the effect of contact point-specific deep brain stimulation in Parkinson's disease: A proof of concept study

BACKGROUND

Deep brain stimulation (DBS) of the subthalamic nucleus (STN) is an effective treatment for disabling fluctuations in motor symptoms in Parkinson's disease (PD) patients. However, iterative exploration of all individual contact points (four in each STN) by the clinician for optimal clinical effects may take months.

OBJECTIVE

In this proof of concept study we explored whether magnetoencephalography (MEG) has the potential to noninvasively measure the effects of changing the active contact point of STN-DBS on spectral power and functional connectivity in PD patients, with the ultimate aim to aid in the process of selecting the optimal contact point, and perhaps reduce the time to achieve optimal stimulation settings.

METHODS

The study included 30 PD patients who had undergone bilateral DBS of the STN. MEG was recorded during stimulation of each of the eight contact points separately (four on each side). Each stimulation position was projected on a vector running through the longitudinal axis of the STN, leading to one scalar value indicating a more dorsolateral or ventromedial contact point position. Using linear mixed models, the stimulation positions were correlated with band-specific absolute spectral power and functional connectivity of i) the motor cortex ipsilateral tot the stimulated side, ii) the whole brain.

RESULTS

At group level, more dorsolateral stimulation was associated with lower low-beta absolute band power in the ipsilateral motor cortex (p = .019). More ventromedial stimulation was associated with higher whole-brain absolute delta (p = .001) and theta (p = .005) power, as well as higher whole-brain theta band functional connectivity (p = .040). At the level of the individual patient, switching the active contact point caused significant changes in spectral power, but the results were highly variable.

CONCLUSIONS

We demonstrate for the first time that stimulation of the dorsolateral (motor) STN in PD patients is associated with lower low-beta power values in the motor cortex. Furthermore, our group-level data show that the location of the active contact point correlates with whole-brain brain activity and connectivity. As results in individual patients were quite variable, it remains unclear if MEG is useful in the selection of the optimal DBS contact point.

Cortical and Subcortical Changes in MEG Activity Reflect Parkinson's Progression over a Period of 7 Years

In this study of early functional changes in Parkinson's disease (PD), we aimed to provide a comprehensive assessment of the development of changes in both cortical and subcortical neurophysiological brain activity, including their association with clinical measures of disease severity. Repeated resting-state MEG recordings and clinical assessments were obtained in the context of a unique longitudinal cohort study over a seven-year period using a multiple longitudinal design. We used linear mixed-models to analyze the relationship between neurophysiological (spectral power and functional connectivity) and clinical data. At baseline, early-stage (drug-naïve) PD patients demonstrated spectral slowing compared to healthy controls in both subcortical and cortical brain regions, most outspoken in the latter. Over time, spectral slowing progressed in strong association with clinical measures of disease progression (cognitive and motor). Global functional connectivity was not different between groups at baseline and hardly changed over time. Therefore, investigation of associations with clinical measures of disease progression were not deemed useful. An analysis of individual connections demonstrated differences between groups at baseline (higher frontal theta, lower parieto-occipital alpha2 band functional connectivity) and over time in PD patients (increase in frontal delta and theta band functional connectivity). Our results suggest that spectral measures are promising candidates in the search for non-invasive markers of both early-stage PD and of the ongoing disease process.

Resting-state oscillations reveal disturbed excitation-inhibition ratio in Alzheimer's disease patients

An early disruption of neuronal excitation-inhibition (E-I) balance in preclinical animal models of Alzheimer's disease (AD) has been frequently reported, but is difficult to measure directly and non-invasively in humans. Here, we examined known and novel neurophysiological measures sensitive to E-I in patients across the AD continuum. Resting-state magnetoencephalography (MEG) data of 86 amyloid-biomarker-confirmed subjects across the AD continuum (17 patients diagnosed with subjective cognitive decline, 18 with mild cognitive impairment (MCI) and 51 with dementia due to probable AD (AD dementia)), 46 healthy elderly and 20 young control subjects were reconstructed to source-space. E-I balance was investigated by detrended fluctuation analysis (DFA), a functional E/I (fE/I) algorithm, and the aperiodic exponent of the power spectrum. We found a disrupted E-I ratio in AD dementia patients specifically, by a lower DFA, and a shift towards higher excitation, by a higher fE/I and a lower aperiodic exponent. Healthy subjects showed lower fE/I ratios (< 1.0) than reported in previous literature, not explained by age or choice of an arbitrary threshold parameter, which warrants caution in interpretation of fE/I results. Correlation analyses showed that a lower DFA (E-I imbalance) and a lower aperiodic exponent (more excitation) was associated with a worse cognitive score in AD dementia patients. In contrast, a higher DFA in the hippocampi of MCI patients was associated with a worse cognitive score. This MEG-study showed E-I imbalance, likely due to increased excitation, in AD dementia, but not in early stage AD patients. To accurately determine the direction of shift in E-I balance, validations of the currently used markers and additional in vivo markers of E-I are required.

Altered directional functional connectivity underlies post-stroke cognitive recovery

Cortical ischaemic strokes result in cognitive deficits depending on the area of the affected brain. However, we have demonstrated that difficulties with attention and processing speed can occur even with small subcortical infarcts. Symptoms appear independent of lesion location, suggesting they arise from generalized disruption of cognitive networks. Longitudinal studies evaluating directional measures of functional connectivity in this population are lacking. We evaluated six patients with minor stroke exhibiting cognitive impairment 6-8 weeks post-infarct and four age-similar controls. Resting-state magnetoencephalography data were collected. Clinical and imaging evaluations of both groups were repeated 6- and 12 months later. Network Localized Granger Causality was used to determine differences in directional connectivity between groups and across visits, which were correlated with clinical performance. Directional connectivity patterns remained stable across visits for controls. After the stroke, inter-hemispheric connectivity between the frontoparietal cortex and the non-frontoparietal cortex significantly increased between visits 1 and 2, corresponding to uniform improvement in reaction times and cognitive scores. Initially, the majority of functional links originated from non-frontal areas contralateral to the lesion, connecting to ipsilesional brain regions. By visit 2, inter-hemispheric connections, directed from the ipsilesional to the contralesional cortex significantly increased. At visit 3, patients demonstrating continued favourable cognitive recovery showed less reliance on these inter-hemispheric connections. These changes were not observed in those without continued improvement. Our findings provide supporting evidence that the neural basis of early post-stroke cognitive dysfunction occurs at the network level, and continued recovery correlates with the evolution of inter-hemispheric connectivity.

Validating MEG source imaging of resting state oscillatory patterns with an intracranial EEG atlas

BACKGROUND

Magnetoencephalography (MEG) is a widely used non-invasive tool to estimate brain activity with high temporal resolution. However, due to the ill-posed nature of the MEG source imaging (MSI) problem, the ability of MSI to identify accurately underlying brain sources along the cortical surface is still uncertain and requires validation.

METHOD

We validated the ability of MSI to estimate the background resting state activity of 45 healthy participants by comparing it to the intracranial EEG (iEEG) atlas (https://mni-open-ieegatlas.

RESEARCH

mcgill.ca/). First, we applied wavelet-based Maximum Entropy on the Mean (wMEM) as an MSI technique. Next, we converted MEG source maps into intracranial space by applying a forward model to the MEG-reconstructed source maps, and estimated virtual iEEG (ViEEG) potentials on each iEEG channel location; we finally quantitatively compared those with actual iEEG signals from the atlas for 38 regions of interest in the canonical frequency bands.

RESULTS

The MEG spectra were more accurately estimated in the lateral regions compared to the medial regions. The regions with higher amplitude in the ViEEG than in the iEEG were more accurately recovered. In the deep regions, MEG-estimated amplitudes were largely underestimated and the spectra were poorly recovered. Overall, our wMEM results were similar to those obtained with minimum norm or beamformer source localization. Moreover, the MEG largely overestimated oscillatory peaks in the alpha band, especially in the anterior and deep regions. This is possibly due to higher phase synchronization of alpha oscillations over extended regions, exceeding the spatial sensitivity of iEEG but detected by MEG. Importantly, we found that MEG-estimated spectra were more comparable to spectra from the iEEG atlas after the aperiodic components were removed.

CONCLUSION

This study identifies brain regions and frequencies for which MEG source analysis is likely to be reliable, a promising step towards resolving the uncertainty in recovering intracerebral activity from non-invasive MEG studies.

Effects of gamma frequency binaural beats on attention and anxiety

Binaural beats (BB) are an auditory phenomenon produced from a combination of two sine waves with slightly different frequencies presented to each ear. Previous research has implicated the role of BBs through brainwave entrainment in potentially giving rise to benefits ranging from enhanced memory and attention to reduced anxiety and stress. Here, we investigated the effect of gamma (40-Hz) BBs on attention using the attention network test (ANT), a previously unused task that assesses three subtypes of attention: Alerting, Orienting, and Executive Control. Fifty-eight healthy adults performed the ANT remotely under the exposure of 340-Hz BBs and a 380-Hz control tone. All completed a rating scale for levels of anxiety before and after each exposure. Performance on the ANT task (reaction time and error rates) between BB and control groups was evaluated using Wilcoxon signed-rank tests. We found no significant differences in Reaction Time (RT), Error Rate (ER), or the efficacy of the Attention Networks (AN) between the experimental and control conditions (p > 0.05). We found no effect of BB on self-rated measures of anxiety. Our findings do not provide evidence for improvement in attention with gamma BB.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12144-023-04681-3.

Biomagnetism: The First Sixty Years

Biomagnetism is the measurement of the weak magnetic fields produced by nerves and muscle. The magnetic field of the heart-the magnetocardiogram (MCG)-is the largest biomagnetic signal generated by the body and was the first measured. Magnetic fields have been detected from isolated tissue, such as a peripheral nerve or cardiac muscle, and these studies have provided insights into the fundamental properties of biomagnetism. The magnetic field of the brain-the magnetoencephalogram (MEG)-has generated much interest and has potential clinical applications to epilepsy, migraine, and psychiatric disorders. The biomagnetic inverse problem, calculating the electrical sources inside the brain from magnetic field recordings made outside the head, is difficult, but several techniques have been introduced to solve it. Traditionally, biomagnetic fields are recorded using superconducting quantum interference device (SQUID) magnetometers, but recently, new sensors have been developed that allow magnetic measurements without the cryogenic technology required for SQUIDs.

Brain-wide network analysis of resting-state neuromagnetic data

The present study performed a brain-wide network analysis of resting-state magnetoencephalograms recorded from 53 healthy participants to visualize elaborate brain maps of phase- and amplitude-derived graph-theory metrics at different frequencies. To achieve this, we conducted a vertex-wise computation of threshold-independent graph metrics by combining proportional thresholding and a conjunction analysis and applied them to a correlation analysis of age and brain networks. Source power showed a frequency-dependent cortical distribution. Threshold-independent graph metrics derived from phase- and amplitude-based connectivity showed similar or different distributions depending on frequency. Vertex-wise age-brain correlation maps revealed that source power at the beta band and the amplitude-based degree at the alpha band changed with age in local regions. The present results indicate that a brain-wide analysis of neuromagnetic data has the potential to reveal neurophysiological network features in the human brain in a resting state.

Non-invasive measurements of ictal and interictal epileptiform activity using optically pumped magnetometers

Magneto- and electroencephalography (MEG/EEG) are important techniques for the diagnosis and pre-surgical evaluation of epilepsy. Yet, in current cryogen-based MEG systems the sensors are offset from the scalp, which limits the signal-to-noise ratio (SNR) and thereby the sensitivity to activity from deep structures such as the hippocampus. This effect is amplified in children, for whom adult-sized fixed-helmet systems are typically too big. Moreover, ictal recordings with fixed-helmet systems are problematic because of limited movement tolerance and/or logistical considerations. Optically Pumped Magnetometers (OPMs) can be placed directly on the scalp, thereby improving SNR and enabling recordings during seizures. We aimed to demonstrate the performance of OPMs in a clinical population. Seven patients with challenging cases of epilepsy underwent MEG recordings using a 12-channel OPM-system and a 306-channel cryogen-based whole-head system: three adults with known deep or weak (low SNR) sources of interictal epileptiform discharges (IEDs), along with three children with focal epilepsy and one adult with frequent seizures. The consistency of the recorded IEDs across the two systems was assessed. In one patient the OPMs detected IEDs that were not found with the SQUID-system, and in two patients no IEDs were found with either system. For the other patients the OPM data were remarkably consistent with the data from the cryogenic system, noting that these were recorded in different sessions, with comparable SNRs and IED-yields overall. Importantly, the wearability of OPMs enabled the recording of seizure activity in a patient with hyperkinetic movements during the seizure. The observed ictal onset and semiology were in agreement with previous video- and stereo-EEG recordings. The relatively affordable technology, in combination with reduced running and maintenance costs, means that OPM-based MEG could be used more widely than current MEG systems, and may become an affordable alternative to scalp EEG, with the potential benefits of increased spatial accuracy, reduced sensitivity to volume conduction/field spread, and increased sensitivity to deep sources. Wearable MEG thus provides an unprecedented opportunity for epilepsy, and given its patient-friendliness, we envisage that it will not only be used for presurgical evaluation of epilepsy patients, but also for diagnosis after a first seizure.

Alpha and gamma EEG coherence during on-task and mind wandering states in schizophrenia

OBJECTIVE

Electroencephalographic (EEG) coherence is one of the most relevant physiological measures used to detect abnormalities in patients with schizophrenia. The present study applies a task-related EEG coherence approach to understand cognitive processing in patients with schizophrenia and healthy controls.

METHODS

EEG coherence for alpha and gamma frequency bands was analyzed in a group of patients with schizophrenia and a group of healthy controls during the performance of an ecological task of sustained attention. We compared EEG coherence when participants presented externally directed cognitive states (On-Task) and when they presented cognitive distraction episodes (Mind-Wandering).

RESULTS

Results reflect cortical differences between groups (higher coherence for schizophrenia in the frontocentral and fronto-temporal regions, and higher coherence for healthy-controls in the postero-central regions), especially in the On-Task condition for the alpha band, compared to Mind-Wandering episodes. Few individual differences in gamma coherence were found.

CONCLUSIONS

The current study provides evidence of neurophysiological differences underlying different cognitive states in schizophrenia and healthy controls.

SIGNIFICANCE

Differences between groups may reflect inhibitory processes necessary for the successful processing of information, especially in the alpha band, given its role in cortical inhibition processes. Patients may activate compensatory inhibitory mechanisms when performing the task, reflected in increased coherence in fronto-temporal regions.

Association between cognitive impairments and aberrant dynamism of overlapping brain sub-networks in unmedicated major depressive disorder: A resting-state MEG study

OBJECTIVE

Little is known about the pathogenesis underlying cognitive impairment in major depressive disorder (MDD). We aimed to explore the mechanisms of cognitive impairments among patients with MDD by investigating the dynamics of overlapping brain sub-networks.

METHODS

Forty unmedicated patients with MDD and 28 healthy controls (HC) were enrolled in this study. Cognitive function was measured using the Chinese versions of MATRICS Consensus Cognitive Battery (MCCB). All participants were scanned using a whole-head resting-state magnetoencephalography (MEG) machine. The dynamism of neural sub-networks was analyzed based on the detection of overlapping communities in five frequency bands of oscillatory brain signals.

RESULTS

MDD demonstrated poorer cognitive performance in six domains compared to HC. The difference in community detection (functional integration mode) in MDD was frequency-dependent. MDD showed significantly decreased community dynamics in all frequency bands compared to HC. Specifically, differences in the visual network (VN) and default mode network (DMN) were detected in all frequency bands, differences in the cognitive control network (CCN) were detected in the alpha2 and beta frequency bands, and differences in the bilateral limbic network (BLN) were only detected in the beta frequency band. Moreover, community dynamics in the alpha2 frequency band were positively correlated with verbal learning and reasoning problem solving abilities in MDD.

CONCLUSIONS

Our study found that decreasing in the dynamics of overlapping sub-networks may differ by frequency bands. The aberrant dynamics of overlapping neural sub-networks revealed by frequency-specific MEG signals may provide new information on the mechanism of cognitive impairments that result from MDD.

Multimodal multilayer network centrality relates to executive functioning

Executive functioning (EF) is a higher order cognitive process that is thought to depend on a network organization facilitating integration across subnetworks, in the context of which the central role of the fronto-parietal network (FPN) has been described across imaging and neurophysiological modalities. However, the potentially complementary unimodal information on the relevance of the FPN for EF has not yet been integrated. We employ a multilayer framework to allow for integration of different modalities into one 'network of networks.' We used diffusion MRI, resting-state functional MRI, MEG, and neuropsychological data obtained from 33 healthy adults to construct modality-specific single-layer networks as well as a single multilayer network per participant. We computed single-layer and multilayer eigenvector centrality of the FPN as a measure of integration in this network and examined their associations with EF. We found that higher multilayer FPN centrality, but not single-layer FPN centrality, was related to better EF. We did not find a statistically significant change in explained variance in EF when using the multilayer approach as compared to the single-layer measures. Overall, our results show the importance of FPN integration for EF and underline the promise of the multilayer framework toward better understanding cognitive functioning.

Preliminary Observations of Resting-State Magnetoencephalography in Nonmedicated Children with Obsessive-Compulsive Disorder

Background: Cortico-striato-thalamo-cortical (CSTC) network alterations are hypothesized to contribute to symptoms of obsessive-compulsive disorder (OCD). To date, very few studies have examined whether CSTC network alterations are present in children with OCD, who are medication naive. Medication-naive pediatric imaging samples may be optimal to study neural correlates of illness and identify brain-based markers, given the proximity to illness onset. Methods: Magnetoencephalography (MEG) data were analyzed at rest, in 18 medication-naive children with OCD (M = 12.1 years ±2.0 standard deviation [SD]; 10 M/8 F) and 13 typically developing children (M = 12.3 years ±2.2 SD; 6 M/7 F). Whole-brain MEG-derived resting-state functional connectivity (rs-fc), for alpha- and gamma-band frequencies were compared between OCD and typically developing (control) groups. Results: Increased MEG-derived rs-fc across alpha- and gamma-band frequencies was found in the OCD group compared to the control group. Increased MEG-derived rs-fc at alpha-band frequencies was evident across a number of regions within the CSTC circuitry and beyond, including the cerebellum and limbic regions. Increased MEG-derived rs-fc at gamma-band frequencies was restricted to the frontal and temporal cortices. Conclusions: This MEG study provides preliminary evidence of altered alpha and gamma networks, at rest, in medication-naive children with OCD. These results support prior findings pointing to the relevance of CSTC circuitry in pediatric OCD and further support accumulating evidence of altered connectivity between regions that extend beyond this network, including the cerebellum and limbic regions. Given the substantial portion of children and youth whose OCD symptoms do not respond to conventional treatments, our findings have implications for future treatment innovation research aiming to target and track whether brain patterns associated with having OCD may change with treatment and/or predict treatment response.

A multiscale brain network model links Alzheimer’s disease-mediated neuronal hyperactivity to large-scale oscillatory slowing

Background

Neuronal hyperexcitability and inhibitory interneuron dysfunction are frequently observed in preclinical animal models of Alzheimer’s disease (AD). This study investigates whether these microscale abnormalities explain characteristic large-scale magnetoencephalography (MEG) activity in human early-stage AD patients.

Methods

To simulate spontaneous electrophysiological activity, we used a whole-brain computational network model comprised of 78 neural masses coupled according to human structural brain topology. We modified relevant model parameters to simulate six literature-based cellular scenarios of AD and compare them to one healthy and six contrast (non-AD-like) scenarios. The parameters include excitability, postsynaptic potentials, and coupling strength of excitatory and inhibitory neuronal populations. Whole-brain spike density and spectral power analyses of the simulated data reveal mechanisms of neuronal hyperactivity that lead to oscillatory changes similar to those observed in MEG data of 18 human prodromal AD patients compared to 18 age-matched subjects with subjective cognitive decline.

Results

All but one of the AD-like scenarios showed higher spike density levels, and all but one of these scenarios had a lower peak frequency, higher spectral power in slower (theta, 4–8Hz) frequencies, and greater total power. Non-AD-like scenarios showed opposite patterns mainly, including reduced spike density and faster oscillatory activity. Human AD patients showed oscillatory slowing (i.e., higher relative power in the theta band mainly), a trend for lower peak frequency and higher total power compared to controls. Combining model and human data, the findings indicate that neuronal hyperactivity can lead to oscillatory slowing, likely due to hyperexcitation (by hyperexcitability of pyramidal neurons or greater long-range excitatory coupling) and/or disinhibition (by reduced excitability of inhibitory interneurons or weaker local inhibitory coupling strength) in early AD.

Conclusions

Using a computational brain network model, we link findings from different scales and models and support the hypothesis of early-stage neuronal hyperactivity underlying E/I imbalance and whole-brain network dysfunction in prodromal AD.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13195-022-01041-4.

Network alterations in eating epilepsy during resting state and during eating using Magnetoencephalography

OBJECTIVE

Eating epilepsy presents various imaging and electrophysiological features along with various seizure triggers. As such, network changes in eating epilepsy have not been comprehensively explored. This study was conducted to illustrate resting state network changes in eating epilepsy and to study the changes in network configurations during eating.

METHODS

Magnetoencephalography recordings of nineteen patients with drug-resistant eating epilepsy were compared with healthy controls during resting state. A subgroup of nine patients and 12 controls had MEG recordings during eating. Network changes were analyzed using phase lag index across 5 frequency bands [delta, theta, alpha, beta, and gamma] using clustering coefficient (CC), betweenness centrality (BC), path length (PL), modularity (Q), and small worldness (SW).

RESULTS

During the resting state, PL was decreased in patients with epilepsy in the delta, theta, and gamma band. Q was lower in patients with epilepsy in the beta and gamma bands. During eating, in patients with epilepsy, PL and SW were increased in all frequency bands, and Q was decreased in the beta band and increased in the rest of the frequency bands. Patients with mixed types of seizures showed higher PL in all bands except alpha, higher Q in all bands, and higher SW in the alpha and beta bands. Node-wise changes in CC and BC implicated changes in DMN and 'eating' networks.

CONCLUSION

Reflex Eating epilepsy presents with a hyperconnected network that exacerbates during eating. The cause of seizure onset and loss of consciousness in eating epilepsy might be due to aberrant network interaction between the regions of the brain involved with eating, such as the sensorimotor cortex, lateral parietal cortex, and insula with the limbic cortex and default mode network across multiple frequency bands.

Modulation of cortical beta oscillations influences motor vigor: A rhythmic TMS-EEG study

Previous electro- or magnetoencephalography (Electro/Magneto EncephaloGraphic; E/MEG) studies using a correlative approach have shown that β (13-30 Hz) oscillations emerging in the primary motor cortex (M1) are implicated in regulating motor response vigor and associated with an anti-kinetic role, that is, slowness of movement. However, the functional role of M1 β oscillations in regulation of motor responses remains unclear. To address this gap, we combined EEG with rhythmic TMS (rhTMS) delivered to M1 at the β (20 Hz) frequency shortly before subjects performed an isometric ramp-and-hold finger force production task at three force levels. rhTMS is a novel approach that can modulate rhythmic patterns of neural activity. β-rhTMS over M1 induced a modulation of neural oscillations to β frequency in the sensorimotor area and reduced peak force rate during the ramp-up period relative to sham and catch trials. Interestingly, this rhTMS effect occurred only in the large force production condition. To distinguish whether the effects of rhTMS on EEG and behavior stemmed from phase-resetting by each magnetic pulse or neural entrainment by the periodicity of rhTMS, we performed a control experiment using arrhythmic TMS (arTMS). arTMS did not induce changes in EEG oscillations nor peak force rate during the rump-up period. Our results provide novel evidence that β neural oscillations emerging the sensorimotor area influence the regulation of motor response vigor. Furthermore, our findings further demonstrate that rhTMS is a promising tool for tuning neural oscillations to the target frequency.

Modulations of local synchrony over time lead to resting-state functional connectivity in a parsimonious large-scale brain model

Biophysical models of large-scale brain activity are a fundamental tool for understanding the mechanisms underlying the patterns observed with neuroimaging. These models combine a macroscopic description of the within- and between-ensemble dynamics of neurons within a single architecture. A challenge for these models is accounting for modulations of within-ensemble synchrony over time. Such modulations in local synchrony are fundamental for modeling behavioral tasks and resting-state activity. Another challenge comes from the difficulty in parametrizing large scale brain models which hinders researching principles related with between-ensembles differences. Here we derive a parsimonious large scale brain model that can describe fluctuations of local synchrony. Crucially, we do not reduce within-ensemble dynamics to macroscopic variables first, instead we consider within and between-ensemble interactions similarly while preserving their physiological differences. The dynamics of within-ensemble synchrony can be tuned with a parameter which manipulates local connectivity strength. We simulated resting-state static and time-resolved functional connectivity of alpha band envelopes in models with identical and dissimilar local connectivities. We show that functional connectivity emerges when there are high fluctuations of local and global synchrony simultaneously (i.e. metastable dynamics). We also show that for most ensembles, leaning towards local asynchrony or synchrony correlates with the functional connectivity with other ensembles, with the exception of some regions belonging to the default-mode network.

A systematic data-driven approach to analyze sensor-level EEG connectivity: Identifying robust phase-synchronized network components using surface Laplacian with spectral-spatial PCA

Although conventional averaging across predefined frequency bands reduces the complexity of EEG functional connectivity (FC), it obscures the identification of resting-state brain networks (RSN) and impedes accurate estimation of FC reliability. Extending prior work, we combined scalp current source density (CSD; spherical spline surface Laplacian) and spectral-spatial PCA to identify FC components. Phase-based FC was estimated via debiased-weighted phase-locking index from CSD-transformed resting EEGs (71 sensors, 8 min, eyes open/closed, 35 healthy adults, 1-week retest). Spectral PCA extracted six robust alpha and theta components (86.6% variance). Subsequent spatial PCA for each spectral component revealed seven robust regionally focused (posterior, central, and frontal) and long-range (posterior-anterior) alpha components (peaks at 8, 10, and 13 Hz) and a midfrontal theta (6 Hz) component, accounting for 37.0% of FC variance. These spatial FC components were consistent with well-known networks (e.g., default mode, visual, and sensorimotor), and four were sensitive to eyes open/closed conditions. Most FC components had good-to-excellent internal consistency (odd/even epochs, eyes open/closed) and test-retest reliability (ICCs ≥ .8). Moreover, the FC component structure was generally present in subsamples (session × odd/even epoch, or smaller subgroups [n = 7-10]), as indicated by high similarity of component loadings across PCA solutions. Apart from systematically reducing FC dimensionality, our approach avoids arbitrary thresholds and allows quantification of meaningful and reliable network components that may prove to be of high relevance for basic and clinical research applications.

Regional healthy brain activity, glioma occurrence and symptomatology

It is unclear why exactly gliomas show preferential occurrence in certain brain areas. Increased spiking activity around gliomas leads to faster tumour growth in animal models, while higher non-invasively measured brain activity is related to shorter survival in patients. However, it is unknown how regional intrinsic brain activity, as measured in healthy controls, relates to glioma occurrence. We first investigated whether gliomas occur more frequently in regions with intrinsically higher brain activity. Second, we explored whether intrinsic cortical activity at individual patients’ tumour locations relates to tumour and patient characteristics.

Across three cross-sectional cohorts, 413 patients were included. Individual tumour masks were created. Intrinsic regional brain activity was assessed through resting-state magnetoencephalography acquired in healthy controls and source-localized to 210 cortical brain regions. Brain activity was operationalized as: (i) broadband power; and (ii) offset of the aperiodic component of the power spectrum, which both reflect neuronal spiking of the underlying neuronal population. We additionally assessed (iii) the slope of the aperiodic component of the power spectrum, which is thought to reflect the neuronal excitation/inhibition ratio. First, correlation coefficients were calculated between group-level regional glioma occurrence, as obtained by concatenating tumour masks across patients, and group-averaged regional intrinsic brain activity. Second, intrinsic brain activity at specific tumour locations was calculated by overlaying patients’ individual tumour masks with regional intrinsic brain activity of the controls and was associated with tumour and patient characteristics.

As proposed, glioma preferentially occurred in brain regions characterized by higher intrinsic brain activity in controls as reflected by higher offset. Second, intrinsic brain activity at patients’ individual tumour locations differed according to glioma subtype and performance status: the most malignant isocitrate dehydrogenase-wild-type glioblastoma patients had the lowest excitation/inhibition ratio at their individual tumour locations as compared to isocitrate dehydrogenase-mutant, 1p/19q-codeleted glioma patients, while a lower excitation/inhibition ratio related to poorer Karnofsky Performance Status, particularly in codeleted glioma patients.

In conclusion, gliomas more frequently occur in cortical brain regions with intrinsically higher activity levels, suggesting that more active regions are more vulnerable to glioma development. Moreover, indices of healthy, intrinsic excitation/inhibition ratio at patients’ individual tumour locations may capture both tumour biology and patients’ performance status. These findings contribute to our understanding of the complex and bidirectional relationship between normal brain functioning and glioma growth, which is at the core of the relatively new field of ‘cancer neuroscience’.