Reliability of Static and Dynamic Network Metrics in the Resting-State: A MEG-Beamformed Connectivity Analysis

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.

Multiplexity of human brain oscillations as a personal brain signature

Human individuality is likely underpinned by the constitution of functional brain networks that ensure consistency of each person's cognitive and behavioral profile. These functional networks should, in principle, be detectable by noninvasive neurophysiology. We use a method that enables the detection of dominant frequencies of the interaction between every pair of brain areas at every temporal segment of the recording period, the dominant coupling modes (DoCM). We apply this method to brain oscillations, measured with magnetoencephalography (MEG) at rest in two independent datasets, and show that the spatiotemporal evolution of DoCMs constitutes an individualized brain fingerprint. Based on this successful fingerprinting we suggest that DoCMs are important targets for the investigation of neural correlates of individual psychological parameters and can provide mechanistic insight into the underlying neurophysiological processes, as well as their disturbance in brain diseases.

Non-invasive mapping of epileptogenic networks predicts surgical outcome

Epilepsy is increasingly considered a disorder of brain networks. Studying these networks with functional connectivity can help identify hubs that facilitate the spread of epileptiform activity. Surgical resection of these hubs may lead patients who suffer from drug-resistant epilepsy to seizure freedom. Here, we aim to map non-invasively epileptogenic networks, through the virtual implantation of sensors estimated with electric and magnetic source imaging, in patients with drug-resistant epilepsy. We hypothesize that highly connected hubs identified non-invasively with source imaging can predict the epileptogenic zone and the surgical outcome better than spikes localized with conventional source localization methods (dipoles). We retrospectively analysed simultaneous high-density electroencephalography (EEG) and magnetoencephalography data recorded from 37 children and young adults with drug-resistant epilepsy who underwent neurosurgery. Using source imaging, we estimated virtual sensors at locations where intracranial EEG contacts were placed. On data with and without spikes, we computed undirected functional connectivity between sensors/contacts using amplitude envelope correlation and phase locking value for physiologically relevant frequency bands. From each functional connectivity matrix, we generated an undirected network containing the strongest connections within sensors/contacts using the minimum spanning tree. For each sensor/contact, we computed graph centrality measures. We compared functional connectivity and their derived graph centrality of sensors/contacts inside resection for good (n = 22, ILAE I) and poor (n = 15, ILAE II-VI) outcome patients, tested their ability to predict the epileptogenic zone in good-outcome patients, examined the association between highly connected hubs removal and surgical outcome and performed leave-one-out cross-validation to support their prognostic value. We also compared the predictive values of functional connectivity with those of dipoles. Finally, we tested the reliability of virtual sensor measures via Spearman's correlation with intracranial EEG at population- and patient-level. We observed higher functional connectivity inside than outside resection (P < 0.05, Wilcoxon signed-rank test) for good-outcome patients, on data with and without spikes across different bands for intracranial EEG and electric/magnetic source imaging and few differences for poor-outcome patients. These functional connectivity measures were predictive of both the epileptogenic zone and outcome (positive and negative predictive values ≥55%, validated using leave-one-out cross-validation) outperforming dipoles on spikes. Significant correlations were found between source imaging and intracranial EEG measures (0.4 ≤ rho ≤ 0.9, P < 0.05). Our findings suggest that virtual implantation of sensors through source imaging can non-invasively identify highly connected hubs in patients with drug-resistant epilepsy, even in the absence of frank epileptiform activity. Surgical resection of these hubs predicts outcome better than dipoles.

Disrupted dynamic functional network connectivity in fetal alcohol spectrum disorders

BACKGROUND

Prenatal alcohol exposure (PAE) can result in harmful and long-lasting neurodevelopmental changes. Children with PAE or a fetal alcohol spectrum disorder (FASD) have decreased white matter volume and resting-state spectral power compared to typically developing controls (TDC) and impaired resting-state static functional connectivity. The impact of PAE on resting-state dynamic functional network connectivity (dFNC) is unknown.

METHODS

Using eyes-closed and eyes-open magnetoencephalography (MEG) resting-state data, global dFNC statistics and meta-states were examined in 89 children aged 6-16 years (51 TDC, 38 with FASD). Source analyzed MEG data were used as input to group spatial independent component analysis to derive functional networks from which the dFNC was calculated.

RESULTS

During eyes-closed, relative to TDC, participants with FASD spent a significantly longer time in state 2, typified by anticorrelation (i.e., decreased connectivity) within and between default mode network (DMN) and visual network (VN), and state 4, typified by stronger internetwork correlation. The FASD group exhibited greater dynamic fluidity and dynamic range (i.e., entered more states, changed from one meta-state to another more often, and traveled greater distances) than TDC. During eyes-open, TDC spent significantly more time in state 1, typified by positive intra- and interdomain connectivity with modest correlation within the frontal network (FN), while participants with FASD spent a larger fraction of time in state 2, typified by anticorrelation within and between DMN and VN and strong correlation within and between FN, attention network, and sensorimotor network.

CONCLUSIONS

There are important resting-state dFNC differences between children with FASD and TDC. Participants with FASD exhibited greater dynamic fluidity and dynamic range and spent more time in states typified by anticorrelation within and between DMN and VN, and more time in a state typified by high internetwork connectivity. Taken together, these network aberrations indicate that prenatal alcohol exposure has a global effect on resting-state connectivity.

Assessing the Repeatability of Multi-Frequency Multi-Layer Brain Network Topologies Across Alternative Researcher's Choice Paths

There is a growing interest in the neuroscience community on the advantages of multilayer functional brain networks. Researchers usually treated different frequencies separately at distinct functional brain networks. However, there is strong evidence that these networks share complementary information while their interdependencies could reveal novel findings. For this purpose, neuroscientists adopt multilayer networks, which can be described mathematically as an extension of trivial single-layer networks. Multilayer networks have become popular in neuroscience due to their advantage to integrate different sources of information. Here, Ι will focus on the multi-frequency multilayer functional connectivity analysis on resting-state fMRI (rs-fMRI) recordings. However, constructing a multilayer network depends on selecting multiple pre-processing steps that can affect the final network topology. Here, I analyzed the rs-fMRI dataset from a single human performing scanning over a period of 18 months (84 scans in total), and the rs-fMRI dataset containing 25 subjects with 3 repeat scans. I focused on assessing the reproducibility of multi-frequency multilayer topologies exploring the effect of two filtering methods for extracting frequencies from BOLD activity, three connectivity estimators, with or without a topological filtering scheme, and two spatial scales. Finally, I untangled specific combinations of researchers' choices that yield consistently brain networks with repeatable topologies, giving me the chance to recommend best practices over consistent topologies.

The thresholding problem and variability in the EEG graph network parameters

Graph thresholding is a frequently used practice of eliminating the weak connections in brain functional connectivity graphs. The main aim of the procedure is to delete the spurious connections in the data. However, the choice of the threshold is arbitrary, and the effect of the threshold choice is not fully understood. Here we present the description of the changes in the global measures of a functional connectivity graph depending on the different proportional thresholds based on the 146 resting-state EEG recordings. The dynamics is presented in five different synchronization measures (wPLI, ImCoh, Coherence, ciPLV, PPC) in sensors and source spaces. The analysis shows significant changes in the graph's global connectivity measures as a function of the chosen threshold which may influence the outcome of the study. The choice of the threshold could lead to different study conclusions; thus it is necessary to improve the reasoning behind the choice of the different analytic options and consider the adoption of different analytic approaches. We also proposed some ways of improving the procedure of thresholding in functional connectivity research.

Universal Lifespan Trajectories of Source-Space Information Flow Extracted from Resting-State MEG Data

Source activity was extracted from resting-state magnetoencephalography data of 103 subjects aged 18-60 years. The directionality of information flow was computed from the regional time courses using delay symbolic transfer entropy and phase entropy. The analysis yielded a dynamic source connectivity profile, disentangling the direction, strength, and time delay of the underlying causal interactions, producing independent time delays for cross-frequency amplitude-to-amplitude and phase-to-phase coupling. The computation of the dominant intrinsic coupling mode (DoCM) allowed me to estimate the probability distribution of the DoCM independently of phase and amplitude. The results support earlier observations of a posterior-to-anterior information flow for phase dynamics in {α1, α2, β, γ} and an opposite flow (anterior to posterior) in θ. Amplitude dynamics reveal posterior-to-anterior information flow in {α1, α2, γ}, a sensory-motor β-oriented pattern, and an anterior-to-posterior pattern in {δ, θ}. The DoCM between intra- and cross-frequency couplings (CFC) are reported here for the first time and independently for amplitude and phase; in both domains {δ, θ, α1}, frequencies are the main contributors to DoCM. Finally, a novel brain age index (BAI) is introduced, defined as the ratio of the probability distribution of inter- over intra-frequency couplings. This ratio shows a universal age trajectory: a rapid rise from the end of adolescence, reaching a peak in adulthood, and declining slowly thereafter. The universal pattern is seen in the BAI of each frequency studied and for both amplitude and phase domains. No such universal age dependence was previously reported.

Cortical network formation based on subthalamic beta bursts in Parkinson's disease

Recent evidence suggests that beta bursts in subthalamic nucleus (STN) play an important role in Parkinsonian pathophysiology. We studied the spatio-temporal relationship between STN beta bursts and cortical activity in 26 Parkinson's disease (PD) patients undergoing deep brain stimulation (DBS) surgery. Postoperatively, we simultaneously recorded STN local field potentials (LFP) from externalized DBS leads and cortical activity using whole-brain magnetoencephalography. Event-related magnetic fields (ERF) were averaged time-locked to STN beta bursts and subjected to source localization. Our results demonstrate that ERF exhibiting activity significantly different from baseline activity were localized within areas functionally related to associative, limbic, and motor systems as well as regions pertinent for visual and language processing. Our data suggest that STN beta bursts are involved in network formation between STN and cortex. This interaction is in line with the idea of parallel processing within the basal ganglia-cortex loop, specifically within the functional subsystems of the STN (i.e., associative, limbic, motor, and the related cortical areas). ERFs within visual and language-related cortical areas indicate involvement of beta bursts in STN-cortex networks beyond the associative, limbic, and motor loops. In sum, our results highlight the involvement of STN beta bursts in the formation of multiple STN - cortex loops in patients with PD.

Abnormal EEG Signal Energy in the Elderly: A Wavelet Analysis of Event-related Potentials During a Stroop Task

BACKGROUND

Previous work showed that elderly with excess in theta activity in their resting state electroencephalogram (EEG) are at higher risk of cognitive decline than those with a normal EEG. By using event-related potentials (ERP) during a counting Stroop task, our prior work showed that elderly with theta excess have a large P300 component compared with normal EEG group. This increased activity could be related to a higher EEG signal energy used during this task.

NEW METHOD

By wavelet analysis applied to ERP obtained during a counting Stroop task we quantified the energy in the different frequency bands of a group of elderly with altered EEG.

RESULTS

In theta and alpha bands, the total energy was higher in elderly subjects with theta excess, specifically in the stimulus categorization window (258-516 ms). Both groups solved the task with similar efficiency.

COMPARISON WITH EXISTING METHODS

The traditional ERP analysis in elderly compares voltage among conditions and groups for a given time windows, while the frequency composition is not usually examined. We complemented our previous ERP analysis using a wavelet methodology. Furthermore, we showed the advantages of wavelet analysis over Short Time Fourier Transform when exploring EEG signal during this task.

CONCLUSIONS

The higher EEG signal energy in ERP might reflect undergoing neurobiological mechanisms that allow the elderly with theta excess to cope with the cognitive task with similar behavioral results as the normal EEG group. This increased energy could promote a metabolic and cellular dysregulation causing a greater decline in cognitive function.

Reproducibility of graph measures derived from resting-state MEG functional connectivity metrics in sensor and source spaces

Prior studies have used graph analysis of resting‐state magnetoencephalography (MEG) to characterize abnormal brain networks in neurological disorders. However, a present challenge for researchers is the lack of guidance on which network construction strategies to employ. The reproducibility of graph measures is important for their use as clinical biomarkers. Furthermore, global graph measures should ideally not depend on whether the analysis was performed in the sensor or source space. Therefore, MEG data of the 89 healthy subjects of the Human Connectome Project were used to investigate test–retest reliability and sensor versus source association of global graph measures. Atlas‐based beamforming was used for source reconstruction, and functional connectivity (FC) was estimated for both sensor and source signals in six frequency bands using the debiased weighted phase lag index (dwPLI), amplitude envelope correlation (AEC), and leakage‐corrected AEC. Reliability was examined over multiple network density levels achieved with proportional weight and orthogonal minimum spanning tree thresholding. At a 100% density, graph measures for most FC metrics and frequency bands had fair to excellent reliability and significant sensor versus source association. The greatest reliability and sensor versus source association was obtained when using amplitude metrics. Reliability was similar between sensor and source spaces when using amplitude metrics but greater for the source than the sensor space in higher frequency bands when using the dwPLI. These results suggest that graph measures are useful biomarkers, particularly for investigating functional networks based on amplitude synchrony.

MEG-Derived Symptom-Sensitive Biomarkers with Long-Term Test-Retest Reliability

Neuroelectric measures derived from human magnetoencephalographic (MEG) recordings hold promise as aides to diagnosis and treatment monitoring and targeting for chronic sequelae of traumatic brain injury (TBI). This study tests novel MEG-derived regional brain measures of tonic neuroelectric activation for long-term test-retest reliability and sensitivity to symptoms. Resting state MEG recordings were obtained from a normative cohort (CamCAN, baseline: n = 613; mean 16-month follow-up: n = 245) and a chronic symptomatic TBI cohort (TEAM-TBI, baseline: n = 62; mean 6-month follow-up: n = 40). The MEG-derived neuroelectric measures were corrected for the empty-room contribution using a random forest classifier. The mean 16-month correlation between baseline and 16-month follow-up CamCAN measures was 0.67; test-retest reliability was markedly improved in this study compared with previous work. The TEAM-TBI cohort was screened for depression, somatization, and anxiety with the Brief Symptom Inventory and for insomnia with the Insomnia Severity Index and was assessed via adjudication for six clinical syndromes: chronic pain, psychological health, and oculomotor, vestibular, cognitive, and sleep dysfunction. Linear classifiers constructed from the 136 regional measures from each TEAM-TBI cohort member distinguished those with and without each symptom, p < 0.0003 for each, i.e., the tonic regional neuroelectric measures of activation are sensitive to the presence/absence of these symptoms and clinical syndromes. The novel regional MEG-derived neuroelectric measures obtained and tested in this study demonstrate the necessary and sufficient properties to be clinically useful, i.e., good test-retest reliability, sensitivity to symptoms in each individual, and obtainable using automatic processing without human judgement or intervention.

Greek High Phenolic Early Harvest Extra Virgin Olive Oil Reduces the Over-Excitation of Information-Flow Based on Dominant Coupling Mode (DoCM) Model in Patients with Mild Cognitive Impairment: An EEG Resting-State Validation Approach

BACKGROUND

Extra virgin olive oil (EVOO) constitutes a natural compound with high protection over cognitive function that could positively alter brain dynamics and the mixture of within and between-frequency connectivity.

OBJECTIVE

The balance of cross-frequency coupling over within-frequency coupling can build a nonlinearity index (NI) that encapsulates the over-excitation of information flow between brain areas and across experimental time. The present study investigated for the very first time how the Greek High Phenolic Early Harvest Extra Virgin Olive Oil (HP-EH-EVOO) versus Moderate Phenolic (MP-EVOO) and Mediterranean Diet (MeDi) intervention in people with mild cognitive impairment (MCI) could affect their spontaneous EEG dynamic connectivity.

METHODS

Forty-three subjects (14 in MeDi, 16 in MP-EVOO, and 13 in HP-EH-EVOO) followed an EEG resting-state recording session (eyes-open and closed) before and after the treatment. Following our dominant coupling mode model, we built a dynamic integrated dynamic functional connectivity graph that tabulates the functional strength and the dominant coupling mode model of every pair of brain areas.

RESULTS

Signal spectrum within 1-13 Hz and theta/beta ratio have decreased in the HP-EH-EVOO group in the eyes-open condition. The intervention improved the FIDoCM across groups and conditions but was more prominent in the HP-EH-EVOO group (p < 0.001). Finally, we revealed a significant higher post-intervention reduction of NI (ΔNITotal and α) for the HP-EH-EVOO compared to the MP-EVOO and MeDi groups (p < 0.0001).

CONCLUSION

Long-term intervention with HP-EH-EVOO reduced the over-excitation of information flow in spontaneous brain activity and altered the signal spectrum of EEG rhythms.

Dyconnmap: Dynamic connectome mapping-A neuroimaging python module

Despite recent progress in the analysis of neuroimaging data sets, our comprehension of the main mechanisms and principles which govern human brain cognition and function remains incomplete. Network neuroscience makes substantial efforts to manipulate these challenges and provide real answers. For the last decade, researchers have been modelling brain structure and function via a graph or network that comprises brain regions that are either anatomically connected via tracts or functionally via a more extensive repertoire of functional associations. Network neuroscience is a relatively new multidisciplinary scientific avenue of the study of complex systems by pursuing novel ways to analyze, map, store and model the essential elements and their interactions in complex neurobiological systems, particularly the human brain, the most complex system in nature. Due to a rapid expansion of neuroimaging data sets' size and complexity, it is essential to propose and adopt new empirical tools to track dynamic patterns between neurons and brain areas and create comprehensive maps. In recent years, there is a rapid growth of scientific interest in moving functional neuroimaging analysis beyond simplified group or time‐averaged approaches and sophisticated algorithms that can capture the time‐varying properties of functional connectivity. We describe algorithms and network metrics that can capture the dynamic evolution of functional connectivity under this perspective. We adopt the word ‘chronnectome’ (integration of the Greek word ‘Chronos’, which means time, and connectome) to describe this specific branch of network neuroscience that explores how mutually informed brain activity correlates across time and brain space in a functional way. We also describe how good temporal mining of temporally evolved dynamic functional networks could give rise to the detection of specific brain states over which our brain evolved. This characteristic supports our complex human mind. The temporal evolution of these brain states and well‐known network metrics could give rise to new analytic trends. Functional brain networks could also increase the multi‐faced nature of the dynamic networks revealing complementary information. Finally, we describe a python module (https://github.com/makism/dyconnmap) which accompanies this article and contains a collection of dynamic complex network analytics and measures and demonstrates its great promise for the study of a healthy subject's repeated fMRI scans.

Presurgical Evaluation of Epilepsy Using Resting-State MEG Functional Connectivity

Epilepsy is caused by abnormal electrical discharges (clinically identified by electrophysiological recording) in a specific part of the brain [originating in only one part of the brain, namely, the epileptogenic zone (EZ)]. Epilepsy is now defined as an archetypical hyperexcited neural network disorder. It can be investigated through the network analysis of interictal discharges, ictal discharges, and resting-state functional connectivity. Currently, there is an increasing interest in embedding resting-state connectivity analysis into the preoperative evaluation of epilepsy. Among the various neuroimaging technologies employed to achieve brain functional networks, magnetoencephalography (MEG) with the excellent temporal resolution is an ideal tool for estimating the resting-state connectivity between brain regions, which can reveal network abnormalities in epilepsy. What value does MEG resting-state functional connectivity offer for epileptic presurgical evaluation? Regarding this topic, this paper introduced the origin of MEG and the workflow of constructing source-space functional connectivity based on MEG signals. Resting-state functional connectivity abnormalities correlate with epileptogenic networks, which are defined by the brain regions involved in the production and propagation of epileptic activities. This paper reviewed the evidence of altered epileptic connectivity based on low- or high-frequency oscillations (HFOs) and the evidence of the advantage of using simultaneous MEG and intracranial electroencephalography (iEEG) recordings. More importantly, this review highlighted that MEG-based resting-state functional connectivity has the potential to predict postsurgical outcomes. In conclusion, resting-state MEG functional connectivity has made a substantial progress toward serving as a candidate biomarker included in epileptic presurgical evaluations.

Reconfiguration of αmplitude driven dominant coupling modes (DoCM) mediated by α-band in adolescents with schizophrenia spectrum disorders

Electroencephalography (EEG) based biomarkers have been shown to correlate with the presence of psychotic disorders. Increased delta and decreased alpha power in psychosis indicate an abnormal arousal state. We investigated brain activity across the basic EEG frequencies and also dynamic functional connectivity of both intra and cross-frequency coupling that could reveal a neurophysiological biomarker linked to an aberrant modulating role of alpha frequency in adolescents with schizophrenia spectrum disorders (SSDs). A dynamic functional connectivity graph (DFCG) has been estimated using the imaginary part of phase lag value (iPLV) and correlation of the envelope (corrEnv). We analyzed DFCG profiles of electroencephalographic resting state (eyes closed) recordings of healthy controls (HC) (n = 39) and SSDs subjects (n = 45) in basic frequency bands {δ,θ,α₁,α₂,β₁,β₂,γ}. In our analysis, we incorporated both intra and cross-frequency coupling modes. Adopting our recent Dominant Coupling Mode (DοCM) model leads to the construction of an integrated DFCG (iDFCG) that encapsulates the functional strength and the DοCM of every pair of brain areas. We revealed significantly higher ratios of delta/alpha₁,₂ power spectrum in SSDs subjects versus HC. The probability distribution (PD) of amplitude driven DoCM mediated by alpha frequency differentiated SSDs from HC with absolute accuracy (100%). The network Flexibility Index (FI) was significantly lower for subjects with SSDs compared to the HC group. Our analysis supports the central role of alpha frequency alterations in the neurophysiological mechanisms of SSDs. Currents findings open up new diagnostic pathways to clinical detection of SSDs and support the design of rational neurofeedback training.

Symptom-Dependent Changes in MEG-Derived Neuroelectric Brain Activity in Traumatic Brain Injury Patients with Chronic Symptoms

Neuroelectric measures derived from human magnetoencephalographic (MEG) recordings hold promise as aides to diagnosis and treatment monitoring and targeting for chronic sequelae of traumatic brain injury (TBI). This study tests novel MEG-derived regional brain measures of tonic neuroelectric activation for long-term test-retest reliability and sensitivity to symptoms. Resting state MEG recordings were obtained from a normative cohort, Cambridge Centre for Ageing and Neuroscience (CamCAN), baseline: n = 619; mean 16-month follow-up: n = 253) and a chronic symptomatic TBI cohort, Targeted Evaluation, Action and Monitoring of Traumatic Brain Injury (TEAM-TBI), baseline: n = 64; mean 6-month follow-up: n = 39). For the CamCAN cohort, MEG-derived neuroelectric measures showed good long-term test-retest reliability for most of the 103 automatically identified stereotypic regions. The TEAM-TBI cohort was screened for depression, somatization, and anxiety with the Brief Symptom Inventory and for insomnia with the Insomnia Severity Index. Linear classifiers constructed from the 103 regional measures from each TEAM-TBI cohort member distinguished those with and without each symptom, with p < 0.01 for each-i.e., the tonic regional neuroelectric measures of activation are sensitive to the presence/absence of these symptoms. The novel regional MEG-derived neuroelectric measures obtained and tested in this study demonstrate the necessary and sufficient properties to be clinically useful-i.e., good test-retest reliability, sensitivity to symptoms in each individual, and obtainable using automatic processing without human judgement or intervention.

Mild traumatic brain injury is associated with dysregulated neural network functioning in children and adolescents

Mild traumatic brain injury is highly prevalent in paediatric populations, and can result in chronic physical, cognitive and emotional impairment, known as persistent post-concussive symptoms. Magnetoencephalography has been used to investigate neurophysiological dysregulation in mild traumatic brain injury in adults; however, whether neural dysrhythmia persists in chronic mild traumatic brain injury in children and adolescents is largely unknown. We predicted that children and adolescents would show similar dysfunction as adults, including pathological slow-wave oscillations and maladaptive, frequency-specific, alterations to neural connectivity. Using magnetoencephalography, we investigated regional oscillatory power and distributed brain-wide networks in a cross-sectional sample of children and adolescents in the chronic stages of mild traumatic brain injury. Additionally, we used a machine learning pipeline to identify the most relevant magnetoencephalography features for classifying mild traumatic brain injury and to test the relative classification performance of regional power versus functional coupling. Results revealed that the majority of participants with chronic mild traumatic brain injury reported persistent post-concussive symptoms. For neurophysiological imaging, we found increased regional power in the delta band in chronic mild traumatic brain injury, predominantly in bilateral occipital cortices and in the right inferior temporal gyrus. Those with chronic mild traumatic brain injury also showed dysregulated neuronal coupling, including decreased connectivity in the delta range, as well as hyper-connectivity in the theta, low gamma and high gamma bands, primarily involving frontal, temporal and occipital brain areas. Furthermore, our multivariate classification approach combined with functional connectivity data outperformed regional power in terms of between-group classification accuracy. For the first time, we establish that local and large-scale neural activity are altered in youth in the chronic phase of mild traumatic brain injury, with the majority presenting persistent post-concussive symptoms, and that dysregulated interregional neural communication is a reliable marker of lingering paediatric ‘mild’ traumatic brain injury.

Magnetoencephalography: physics, techniques, and applications in the basic and clinical neurosciences

Magnetoencephalography (MEG) is a technique used to measure the magnetic fields generated from neuronal activity in the brain. MEG has a high temporal resolution on the order of milliseconds and provides a more direct measure of brain activity when compared with hemodynamic-based neuroimaging methods such as magnetic resonance imaging and positron emission tomography. The current review focuses on basic features of MEG such as the instrumentation and the physics that are integral to the signals that can be measured, and the principles of source localization techniques, particularly the physics of beamforming and the techniques that are used to localize the signal of interest. In addition, we review several metrics that can be used to assess functional coupling in MEG and describe the advantages and disadvantages of each approach. Lastly, we discuss the current and future applications of MEG.

Quantitative Identification of Functional Connectivity Disturbances in Neuropsychiatric Lupus Based on Resting-State fMRI: A Robust Machine Learning Approach

Neuropsychiatric systemic lupus erythematosus (NPSLE) is an autoimmune entity comprised of heterogenous syndromes affecting both the peripheral and central nervous system. Research on the pathophysiological substrate of NPSLE manifestations, including functional neuroimaging studies, is extremely limited. The present study examined person-specific patterns of whole-brain functional connectivity in NPSLE patients (n = 44) and age-matched healthy control participants (n = 39). Static functional connectivity graphs were calculated comprised of connection strengths between 90 brain regions. These connections were subsequently filtered through rigorous surrogate analysis, a technique borrowed from physics, novel to neuroimaging. Next, global as well as nodal network metrics were estimated for each individual functional brain network and were input to a robust machine learning algorithm consisting of a random forest feature selection and nested cross-validation strategy. The proposed pipeline is data-driven in its entirety, and several tests were performed in order to ensure model robustness. The best-fitting model utilizing nodal graph metrics for 11 brain regions was associated with 73.5% accuracy (74.5% sensitivity and 73% specificity) in discriminating NPSLE from healthy individuals with adequate statistical power. Closer inspection of graph metric values suggested an increased role within the functional brain network in NSPLE (indicated by higher nodal degree, local efficiency, betweenness centrality, or eigenvalue efficiency) as compared to healthy controls for seven brain regions and a reduced role for four areas. These findings corroborate earlier work regarding hemodynamic disturbances in these brain regions in NPSLE. The validity of the results is further supported by significant associations of certain selected graph metrics with accumulated organ damage incurred by lupus, with visuomotor performance and mental flexibility scores obtained independently from NPSLE patients.

Aberrant Whole-Brain Transitions and Dynamics of Spontaneous Network Microstates in Mild Traumatic Brain Injury

Dynamic Functional Connectivity (DFC) analysis is a promising approach for the characterization of brain electrophysiological activity. In this study, we investigated abnormal alterations due to mild Traumatic Brain Injury (mTBI) using DFC of the source reconstructed magnetoencephalographic (MEG) resting-state recordings. Brain activity in several well-known frequency bands was first reconstructed using beamforming of the MEG data to determine ninety anatomical brain regions of interest. A DFC graph was formulated using the imaginary part of phase-locking values, which were obtained from 30 mTBI patients and 50 healthy controls (HC). Subsequently, we estimated normalized Laplacian transformations of individual, statistically and topologically filtered quasi-static graphs. The corresponding eigenvalues of each node synchronization were then computed and through the neural-gas algorithm, we quantized the evolution of the eigenvalues resulting in distinct network microstates (NMstates). The discrimination level between the two groups was assessed using an iterative cross-validation classification scheme with features either the NMstates in each frequency band, or the combination of the so-called chronnectomics (flexibility index, occupancy time of NMstate, and Dwell time) with the complexity index over the evolution of the NMstates across all frequency bands. Classification performance based on chronnectomics showed 80% accuracy, 99% sensitivity, and 49% specificity. However, performance was much higher (accuracy: 91-97%, sensitivity: 100%, and specificity: 77-93%) when focusing on the microstates. Exploring the mean node degree within and between brain anatomical networks (default mode network, frontoparietal, occipital, cingulo-opercular, and sensorimotor), a reduced pattern occurred from lower to higher frequency bands, with statistically significant stronger degrees for the HC than the mTBI group. A higher entropic profile on the temporal evolution of the modularity index was observed for both NMstates for the mTBI group across frequencies. A significant difference in the flexibility index was observed between the two groups for the β frequency band. The latter finding may support a central role of the thalamus impairment in mTBI. The current study considers a complete set of frequency-dependent connectomic markers of mTBI-caused alterations in brain connectivity that potentially could serve as markers to assess the return of an injured subject back to normality.

Aberrant MEG multi-frequency phase temporal synchronization predicts conversion from mild cognitive impairment-to-Alzheimer's disease

Many neuroimaging studies focus on a frequency-specific or a multi-frequency network analysis showing that functional brain networks are disrupted in patients with Alzheimer's disease (AD). Although those studies enriched our knowledge of the impact of AD in brain's functionality, our goal is to test the effectiveness of combining neuroimaging with network neuroscience to predict with high accuracy subjects with mild cognitive impairment (MCI) that will convert to AD.

In this study, eyes-closed resting-state magnetoencephalography (MEG) recordings from 27 stable MCI (sMCI) and 27 progressive MCI (pMCI) from two scan sessions (baseline and follow-up after approximately 3 years) were projected via beamforming onto an atlas-based set of regions of interest (ROIs). Dynamic functional connectivity networks were constructed independently for the five classical frequency bands while a multivariate phase-based coupling metric was adopted. Thus, computing the distance between the fluctuation of functional strength of every pair of ROIs between the two conditions with dynamic time wrapping (DTW), a large set of features was extracted. A machine learning algorithm revealed 30 DTW-based features in the five frequency bands that can distinguish the sMCI from pMCI with absolute accuracy (100%). Further analysis of the selected links revealed that most of the connected ROIs were part of the default mode network (DMN), the cingulo-opercular (CO), the fronto-parietal and the sensorimotor network.

Overall, our dynamic network multi-frequency analysis approach provides an effective framework of constructing a sensitive MEG-based connectome biomarker for the prediction of conversion from MCI to Alzheimer's disease.

•

49 features in the five frequency bands discriminated the two groups with 100% accuracy.

•

The analysis revealed most of the links were part of important brain networks as the DMN.

•

The alteration of these links, usually affected in AD, is related to synaptic loss.

•

Our multi-frequency approach provided a MEG connectome biomarker to predict the conversion to AD.

How Sensitive Are Conventional MEG Functional Connectivity Metrics With Sliding Windows to Detect Genuine Fluctuations in Dynamic Functional Connectivity?

Despite advances in the field of dynamic connectivity, fixed sliding window approaches for the detection of fluctuations in functional connectivity are still widely used. The use of conventional connectivity metrics in conjunction with a fixed sliding window comes with the arbitrariness of the chosen window lengths. In this paper we use multivariate autoregressive and neural mass models with a priori defined ground truths to systematically analyze the sensitivity of conventional metrics in combination with different window lengths to detect genuine fluctuations in connectivity for various underlying state durations. Metrics of interest are the coherence, imaginary coherence, phase lag index, phase locking value and the amplitude envelope correlation. We performed analysis for two nodes and at the network level. We demonstrate that these metrics show indeed higher variability for genuine temporal fluctuations in connectivity compared to a static connectivity state superimposed by noise. Overall, the error of the connectivity estimates themselves decreases for longer state durations (order of seconds), while correlations of the connectivity fluctuations with the ground truth was higher for longer state durations. In general, metrics, in combination with a sliding window, perform poorly for very short state durations. Increasing the SNR of the system only leads to a moderate improvement. In addition, at the network level, only longer window widths were sufficient to detect plausible resting state networks that matched the underlying ground truth, especially for the phase locking value, amplitude envelope correlation and coherence. The length of these longer window widths did not necessarily correspond to the underlying state durations. For short window widths resting state network connectivity patterns could not be retrieved. We conclude that fixed sliding window approaches for connectivity can detect modulations of connectivity, but mostly if the underlying dynamics operate on moderate to slow timescales. In practice, this can be a drawback, as state durations can vary significantly in empirical data.

Magnetic Source Imaging and Infant MEG: Current Trends and Technical Advances

Magnetoencephalography (MEG) is known for its temporal precision and good spatial resolution in cognitive brain research. Nonetheless, it is still rarely used in developmental research, and its role in developmental cognitive neuroscience is not adequately addressed. The current review focuses on the source analysis of MEG measurement and its potential to answer critical questions on neural activation origins and patterns underlying infants’ early cognitive experience. The advantages of MEG source localization are discussed in comparison with functional magnetic resonance imaging (fMRI) and functional near-infrared spectroscopy (fNIRS), two leading imaging tools for studying cognition across age. Challenges of the current MEG experimental protocols are highlighted, including measurement and data processing, which could potentially be resolved by developing and improving both software and hardware. A selection of infant MEG research in auditory, speech, vision, motor, sleep, cross-modality, and clinical application is then summarized and discussed with a focus on the source localization analyses. Based on the literature review and the advancements of the infant MEG systems and source analysis software, typical practices of infant MEG data collection and analysis are summarized as the basis for future developmental cognitive research.

Is brain connectome research the future frontier for subjective cognitive decline? A systematic review

OBJECTIVE

We performed a systematic literature review on Subjective Cognitive Decline (SCD) in order to examine whether the resemblance of brain connectome and functional connectivity (FC) alterations in SCD with respect to MCI, AD and HC can help us draw conclusions on the progression of SCD to more advanced stages of dementia.

METHODS

We searched for studies that used any neuroimaging tool to investigate potential differences/similarities of brain connectome in SCD with respect to HC, MCI, and AD.

RESULTS

Sixteen studies were finally included in the review. Apparent FC connections and disruptions were observed in the white matter, default mode and gray matter networks in SCD with regards to HC, MCI, and AD. Interestingly, more apparent connections in SCD were located over the posterior regions, while an increase of FC over anterior regions was observed as the disease progressed.

CONCLUSIONS

Elders with SCD display a significant disruption of the brain network, which in most of the cases is worse than HC across multiple network parameters.

SIGNIFICANCE

The present review provides comprehensive and balanced coverage of a timely target research activity around SCD with the intention to identify similarities/differences across patient groups on the basis of brain connectome properties.

Modeling the Switching Behavior of Functional Connectivity Microstates (FCμstates) as a Novel Biomarker for Mild Cognitive Impairment

The need for designing and validating novel biomarkers for the detection of mild cognitive impairment (MCI) is evident. MCI patients have a high risk of developing Alzheimer's disease (AD), and for that reason the introduction of novel and reliable biomarkers is of significant clinical importance. Motivated by recent findings on the rich information of dynamic functional connectivity graphs (DFCGs) about brain (dys) function, we introduced a novel approach of identifying MCI based on magnetoencephalographic (MEG) resting state recordings. The activity of different brain rhythms {δ, 𝜃, α1, α2, β1, β2, γ1, γ2} was first beamformed with linear constrained minimum norm variance in the MEG data to determine 90 anatomical regions of interest (ROIs). A DFCG was then estimated using the imaginary part of phase lag value (iPLV) for both intra-frequency coupling (8) and cross-frequency coupling pairs (28). We analyzed DFCG profiles of neuromagnetic resting state recordings of 18 MCI patients and 22 healthy controls. We followed our model of identifying the dominant intrinsic coupling mode (DICM) across MEG sources and temporal segments, which further leads to the construction of an integrated DFCG (iDFCG). We then filtered statistically and topologically every snapshot of the iDFCG with data-driven approaches. An estimation of the normalized Laplacian transformation for every temporal segment of the iDFCG and the related eigenvalues created a 2D map based on the network metric time series of the eigenvalues (NMTSeigs). The NMTSeigs preserves the non-stationarity of the fluctuated synchronizability of iDCFG for each subject. Employing the initial set of 20 healthy elders and 20 MCI patients, as training set, we built an overcomplete dictionary set of network microstates (n μstates). Afterward, we tested the whole procedure in an extra blind set of 20 subjects for external validation. We succeeded in gaining a high classification accuracy on the blind dataset (85%), which further supports the proposed Markovian modeling of the evolution of brain states. The adaptation of appropriate neuroinformatic tools that combine advanced signal processing and network neuroscience tools could properly manipulate the non-stationarity of time-resolved FC patterns revealing a robust biomarker for MCI.