Sensitive and reproducible MEG resting-state metrics of functional connectivity in Alzheimer's disease

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.

A Multiscale Closed-Loop Neurotoxicity Model of Alzheimer's Disease Progression Explains Functional Connectivity Alterations

The accumulation of amyloid-β (Aβ) and hyperphosphorylated-tau (hp-tau) are two classical histopathological biomarkers in Alzheimer's disease (AD). However, their detailed interactions with the electrophysiological changes at the meso- and macroscale are not yet fully understood. We developed a mechanistic multiscale model of AD progression, linking proteinopathy to its effects on neural activity and vice-versa. We integrated a heterodimer model of prion-like protein propagation and a brain network model of Jansen-Rit neural masses derived from human neuroimaging data whose parameters varied due to neurotoxicity. Results showed that changes in inhibition guided the electrophysiological alterations found in AD, and these changes were mainly attributed to Aβ effects. Additionally, we found a causal disconnection between cellular hyperactivity and interregional hypersynchrony contrary to previous beliefs. Finally, we demonstrated that early Aβ and hp-tau depositions' location determine the spatiotemporal profile of the proteinopathy. The presented model combines the molecular effects of both Aβ and hp-tau together with a mechanistic protein propagation model and network effects within a closed-loop model. This holds the potential to enlighten the interplay between AD mechanisms on various scales, aiming to develop and test novel hypotheses on the contribution of different AD-related variables to the disease evolution.

Local signal variability and functional connectivity: Sensitive measures of the excitation-inhibition ratio?

A novel network version of permutation entropy, the inverted joint permutation entropy (JPEinv), holds potential as non-invasive biomarker of abnormal excitation-inhibition (E-I) ratio in Alzheimer's disease (AD). In this computational modelling study, we test the hypotheses that this metric, and related measures of signal variability and functional connectivity, are sensitive to altered E-I ratios. The E-I ratio in each neural mass of a whole-brain computational network model was systematically varied. We evaluated whether JPEinv, local signal variability (by permutation entropy) and functional connectivity (by weighted symbolic mutual information (wsMI)) were related to E-I ratio, on whole-brain and regional level. The hub disruption index can identify regions primarily affected in terms of functional connectivity strength (or: degree) by the altered E-I ratios. Analyses were performed for a range of coupling strengths, filter and time-delay settings. On whole-brain level, higher E-I ratios were associated with higher functional connectivity (by JPEinv and wsMI) and lower local signal variability. These relationships were nonlinear and depended on the coupling strength, filter and time-delay settings. On regional level, hub-like regions showed a selective decrease in functional degree (by JPEinv and wsMI) upon a lower E-I ratio, and non-hub-like regions showed a selective increase in degree upon a higher E-I ratio. These results suggest that abnormal functional connectivity and signal variability, as previously reported in patients across the AD continuum, can inform us about altered E-I ratios.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11571-023-10003-x.

Neurophysiological trajectories in Alzheimer's disease progression

Alzheimer's disease (AD) is characterized by the accumulation of amyloid-β and misfolded tau proteins causing synaptic dysfunction, and progressive neurodegeneration and cognitive decline. Altered neural oscillations have been consistently demonstrated in AD. However, the trajectories of abnormal neural oscillations in AD progression and their relationship to neurodegeneration and cognitive decline are unknown. Here, we deployed robust event-based sequencing models (EBMs) to investigate the trajectories of long-range and local neural synchrony across AD stages, estimated from resting-state magnetoencephalography. The increases in neural synchrony in the delta-theta band and the decreases in the alpha and beta bands showed progressive changes throughout the stages of the EBM. Decreases in alpha and beta band synchrony preceded both neurodegeneration and cognitive decline, indicating that frequency-specific neuronal synchrony abnormalities are early manifestations of AD pathophysiology. The long-range synchrony effects were greater than the local synchrony, indicating a greater sensitivity of connectivity metrics involving multiple regions of the brain. These results demonstrate the evolution of functional neuronal deficits along the sequence of AD progression.

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.

Network Hyperexcitability in Early-Stage Alzheimer's Disease: Evaluation of Functional Connectivity Biomarkers in a Computational Disease Model

Background

There is increasing evidence from animal and clinical studies that network hyperexcitability (NH) may be an important pathophysiological process and potential target for treatment in early Alzheimer's disease (AD). Measures of functional connectivity (FC) have been proposed as promising biomarkers for NH, but it is unknown which measure has the highest sensitivity for early-stage changes in the excitation/inhibition balance.

Objective

We aim to test the performance of different FC measures in detecting NH at the earliest stage using a computational approach.

Methods

We use a whole brain computational model of activity dependent degeneration to simulate progressive AD pathology and NH. We investigate if and at what stage four measures of FC (amplitude envelope correlation corrected [AECc], phase lag index [PLI], joint permutation entropy [JPE] and a new measure: phase lag time [PLT]) can detect early-stage AD pathophysiology.

Results

The activity dependent degeneration model replicates spectral changes in line with clinical data and demonstrates increasing NH. Compared to relative theta power as a gold standard the AECc and PLI are shown to be less sensitive in detecting early-stage NH and AD-related neurophysiological abnormalities, while the JPE and the PLT show more sensitivity with excellent test characteristics.

Conclusions

Novel FC measures, which are better in detecting rapid fluctuations in neural activity and connectivity, may be superior to well-known measures such as the AECc and PLI in detecting early phase neurophysiological abnormalities and in particular NH in AD. These markers could improve early diagnosis and treatment target identification.

Longitudinal resting-state EEG in amyloid-positive patients along the Alzheimer's disease continuum: considerations for clinical trials

BACKGROUND

To enable successful inclusion of electroencephalography (EEG) outcome measures in Alzheimer's disease (AD) clinical trials, we retrospectively mapped the progression of resting-state EEG measures over time in amyloid-positive patients with mild cognitive impairment (MCI) or dementia due to AD.

METHODS

Resting-state 21-channel EEG was recorded in 148 amyloid-positive AD patients (MCI, n = 88; dementia due to AD, n = 60). Two or more EEG recordings were available for all subjects. We computed whole-brain and regional relative power (i.e., theta (4-8 Hz), alpha1 (8-10 Hz), alpha2 (10-13 Hz), beta (13-30 Hz)), peak frequency, signal variability (i.e., theta permutation entropy), and functional connectivity values (i.e., alpha and beta corrected amplitude envelope correlation, theta phase lag index, weighted symbolic mutual information, inverted joint permutation entropy). Whole-group linear mixed effects models were used to model the development of EEG measures over time. Group-wise analysis was performed to investigate potential differences in change trajectories between the MCI and dementia subgroups. Finally, we estimated the minimum sample size required to detect different treatment effects (i.e., 50% less deterioration, stabilization, or 50% improvement) on the development of EEG measures over time, in hypothetical clinical trials of 1- or 2-year duration.

RESULTS

Whole-group analysis revealed significant regional and global oscillatory slowing over time (i.e., increased relative theta power, decreased beta power), with strongest effects for temporal and parieto-occipital regions. Disease severity at baseline influenced the EEG measures' rates of change, with fastest deterioration reported in MCI patients. Only AD dementia patients displayed a significant decrease of the parieto-occipital peak frequency and theta signal variability over time. We estimate that 2-year trials, focusing on amyloid-positive MCI patients, require 36 subjects per arm (2 arms, 1:1 randomization, 80% power) to detect a stabilizing treatment effect on temporal relative theta power.

CONCLUSIONS

Resting-state EEG measures could facilitate early detection of treatment effects on neuronal function in AD patients. Their sensitivity depends on the region-of-interest and disease severity of the study population. Conventional spectral measures, particularly recorded from temporal regions, present sensitive AD treatment monitoring markers.

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.

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.

Functional connectivity discriminates epileptogenic states and predicts surgical outcome in children with drug resistant epilepsy

Normal brain functioning emerges from a complex interplay among regions forming networks. In epilepsy, these networks are disrupted causing seizures. Highly connected nodes in these networks are epilepsy surgery targets. Here, we assess whether functional connectivity (FC) using intracranial electroencephalography can quantify brain regions epileptogenicity and predict surgical outcome in children with drug resistant epilepsy (DRE). We computed FC between electrodes on different states (i.e. interictal without spikes, interictal with spikes, pre-ictal, ictal, and post-ictal) and frequency bands. We then estimated the electrodes' nodal strength. We compared nodal strength between states, inside and outside resection for good- (n = 22, Engel I) and poor-outcome (n = 9, Engel II-IV) patients, respectively, and tested their utility to predict the epileptogenic zone and outcome. We observed a hierarchical epileptogenic organization among states for nodal strength: lower FC during interictal and pre-ictal states followed by higher FC during ictal and post-ictal states (p < 0.05). We further observed higher FC inside resection (p < 0.05) for good-outcome patients on different states and bands, and no differences for poor-outcome patients. Resection of nodes with high FC was predictive of outcome (positive and negative predictive values: 47-100%). Our findings suggest that FC can discriminate epileptogenic states and predict outcome in patients with DRE.

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.

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.

Source space connectomics of neurodegeneration: One-metric approach does not fit all

Brain functional connectivity in dementia has been assessed with dissimilar EEG connectivity metrics and estimation procedures, thereby increasing results' heterogeneity. In this scenario, joint analyses integrating information from different metrics may allow for a more comprehensive characterization of brain functional interactions in different dementia subtypes. To test this hypothesis, resting-state electroencephalogram (rsEEG) was recorded in individuals with Alzheimer's Disease (AD), behavioral variant frontotemporal dementia (bvFTD), and healthy controls (HCs). Whole-brain functional connectivity was estimated in the EEG source space using 101 different types of functional connectivity, capturing linear and nonlinear interactions in both time and frequency-domains. Multivariate machine learning and progressive feature elimination was run to discriminate AD from HCs, and bvFTD from HCs, based on joint analyses of i) EEG frequency bands, ii) complementary frequency-domain metrics (e.g., instantaneous, lagged, and total connectivity), and iii) time-domain metrics with different linearity assumption (e.g., Pearson correlation coefficient and mutual information). <10% of all possible connections were responsible for the differences between patients and controls, and atypical connectivity was never captured by >1/4 of all possible connectivity measures. Joint analyses revealed patterns of hypoconnectivity (patientsHCs) in both groups was mainly identified in frontotemporal regions. These atypicalities were differently captured by frequency- and time-domain connectivity metrics, in a bandwidth-specific fashion. The multi-metric representation of source space whole-brain functional connectivity evidenced the inadequacy of single-metric approaches, and resulted in a valid alternative for the selection problem in EEG connectivity. These joint analyses reveal patterns of brain functional interdependence that are overlooked with single metrics approaches, contributing to a more reliable and interpretable description of atypical functional connectivity in neurodegeneration.