Causal Interactions between Frontal(θ) - Parieto-Occipital(α2) Predict Performance on a Mental Arithmetic Task

Effects of Mental Workload Manipulation on Electroencephalography Spectrum Oscillation and Microstates in Multitasking Environments

INTRODUCTION

Multitasking during flights leads to a high mental workload, which is detrimental for maintaining task performance. Electroencephalography (EEG) power spectral analysis based on frequency-band oscillations and microstate analysis based on global brain network activation can be used to evaluate mental workload. This study explored the effects of a high mental workload during simulated flight multitasking on EEG frequency-band power and microstate parameters.

METHODS

Thirty-six participants performed multitasking with low and high mental workloads after 4 consecutive days of training. Two levels of mental workload were set by varying the number of subtasks. EEG signals were acquired during the task. Power spectral and microstate analyses were performed on the EEG. The indices of four frequency bands (delta, theta, alpha, and beta) and four microstate classes (A-D) were calculated, changes in the frequency-band power and microstate parameters under different mental workloads were compared, and the relationships between the two types of EEG indices were analyzed.

RESULTS

The theta-, alpha-, and beta-band powers were higher under the high than under the low mental workload condition. Compared with the low mental workload condition, the high mental workload condition had a lower global explained variance and time parameters of microstate B but higher time parameters of microstate D. Less frequent transitions between microstates A and B and more frequent transitions between microstates C and D were observed during high mental workload conditions. The time parameters of microstate B were positively correlated with the delta-, theta-, and beta-band powers, whereas the duration of microstate C was negatively correlated with the beta-band power.

CONCLUSION

EEG frequency-band power and microstate parameters can be used to detect a high mental workload. Power spectral analyses based on frequency-band oscillations and microstate analyses based on global brain network activation were not completely isolated during multitasking.

Unveiling Glutamate Dynamics: Cognitive Demands in Human Short-Term Memory Learning Across Frontal and Parieto-Occipital Cortex: A Functional MRS Study

Background

Acquiring new knowledge necessitates alterations at the synaptic level within the brain. Glutamate, a pivotal neurotransmitter, plays a critical role in these processes, particularly in learning and memory formation. Although previous research has explored glutamate's involvement in cognitive functions, a comprehensive understanding of its real-time dynamics remains elusive during memory tasks.

Objective

This study aimed to investigate glutamate modulation during memory tasks in the right Dorsolateral Prefrontal Cortex (DLPFC) and parieto-occipital regions using functional Magnetic Resonance Spectroscopy (fMRS).

Material and Methods

This experimental research applied fMRS acquisition concurrently with a modified Sternberg's verbal working memory task for fourteen healthy right-handed participants (5 females, mean age=30.64±4.49). The glutamate/total-creatine (Glu/tCr) ratio was quantified by LCModel in the DLPFC and parieto-occipital voxels while applying the tissue corrections.

Results

The significantly higher Glu/tCr modulation was observed during the task with a trend of increased modulation with memory load in both the DLPFC (19.9% higher, P-value=0.018) and parieto-occipital (33% higher, P-value=0.046) regions compared to the rest.

Conclusion

Our pioneering fMRS study has yielded groundbreaking insights into brain functions during S-term Memory (STM) and learning. This research provides valuable methodological advancements for investigating the metabolic functions of both healthy and disordered brains. Based on the findings, cognitive demands directly correlate with glutamate levels, highlighting the neurochemical underpinnings of cognitive processing. Additionally, the obtained results potentially challenge the traditional left-hemisphere-centric model of verbal working memory, leading to the deep vision of hemispheric contributions to cognitive functions.

Study on different brain activation rearrangement during cognitive workload from ERD/ERS and coherence analysis

The functional activities of the brain during any task like imaginary, motor, or cognitive are different in pattern as well as their area of activation in the brain is also different. This variation in pattern is also found in the brain's electrical variations that can be measured from the scalp of the brain using an electroencephalogram (EEG). This work exclusively studied a group of subjects' EEG data (available at: https://archive.physionet.org/physiobank/database/eegmat/) to unravel the activation pattern of the human brain during a mental arithmetic task. Since any cognitive task creates variations in EEG signal pattern, the relative changes in the signal power also occur which is also known as event-related desynchronization/synchronization (ERD/ERS). In this work, ERD/ERS have calculated the band-wise power spectral density (PSD) using Welch's method from the EEG signals. Besides, the coherence analysis was also performed to verify the results of ERD/ERS analysis from several randomly chosen subjects' EEG data. Here, subjects performing mental arithmetic tasks were grouped based on their performances: good (subtraction solved > 10 on average) and bad (subtraction solved ≤ 10 on average) to conduct group-specific ERD/ERS analysis regarding their performance in cognitive tasks. It was found that when the brain is on count condition, the variations found in the band power of theta and beta. The amounts of ERS in the left hemisphere are increased. When the task complexity increases, it contributes to an increase in relative ERD/ERS amounts and durations. The left and right hemispheres were asymmetrically distributed by both the pre-stimulus stages of the corresponding band power and relative ERD/ERS.

Electroencephalogram-based objective assessment of cognitive function level associated with age-related hearing loss

Age-Related Hearing Loss (ARHL) is a common problem in aging. Numerous longitudinal cohort studies have revealed that ARHL is closely related to cognitive function, leading to a significant risk of cognitive decline and dementia. This risk gradually increases with the severity of hearing loss. We designed dual auditory Oddball and cognitive task paradigms for the ARHL subjects, then obtained the Montreal Cognitive Assessment (MoCA) scale evaluation results for all the subjects. Multi-dimensional EEG characteristics helped explore potential biomarkers to evaluate the cognitive level of the ARHL group, having a significantly lower P300 peak amplitude coupled with a prolonged latency. Moreover, visual memory, auditory memory, and logical calculation were investigated during the cognitive task paradigm. In the ARHL groups, the alpha-to-beta rhythm energy ratio in the visual and auditory memory retention period and the wavelet packet entropy value within the logical calculation period were significantly reduced. Correlation analysis between the above specificity indicators and the subjective scale results of the ARHL group revealed that the auditory P300 component characteristics could assess attention resources and information processing speed. The alpha and beta rhythm energy ratio and wavelet packet entropy can become potential indicators to determine working memory and logical cognitive computation-related cognitive ability.

Dynamic Reconfiguration of Dominant Intrinsic Coupling Modes in Elderly at Prodromal Alzheimer's Disease Risk

Large-scale human brain networks interact across both spatial and temporal scales. Especially for electro- and magnetoencephalography (EEG/MEG), there are many evidences that there is a synergy of different subnetworks that oscillate on a dominant frequency within a quasi-stable brain temporal frame. Intrinsic cortical-level integration reflects the reorganization of functional brain networks that support a compensation mechanism for cognitive decline. Here, a computerized intervention integrating different functions of the medial temporal lobes, namely, object-level and scene-level representations, was conducted. One hundred fifty-eight patients with mild cognitive impairment underwent 90 min of training per day over 10 weeks. An active control (AC) group of 50 subjects was exposed to documentaries, and a passive control group of 55 subjects did not engage in any activity. Following a dynamic functional source connectivity analysis, the dynamic reconfiguration of intra- and cross-frequency coupling mechanisms before and after the intervention was revealed. After the neuropsychological and resting state electroencephalography evaluation, the ratio of inter versus intra-frequency coupling modes and also the contribution of β1 frequency was higher for the target group compared to its pre-intervention period. These frequency-dependent contributions were linked to neuropsychological estimates that were improved due to intervention. Additionally, the time-delays of the cortical interactions were improved in {δ, θ, α2, β1} compared to the pre-intervention period. Finally, dynamic networks of the target group further improved their efficiency over the total cost of the network. This is the first study that revealed a dynamic reconfiguration of intrinsic coupling modes and an improvement of time-delays due to a target intervention protocol.

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.

Multilayer Network Framework Reveals Cross-Frequency Coupling Hubs in Cortical Olfactory Perception

Olfactory perception is shaped by dynamic in-teractions among networks of widely distributed brain regions involved in several neurocognitive processes. However, the neural mechanisms that enable effective coordination and integrative processing across these brain regions, which have different functions and operating characteristics, are not yet fully understood. In this study we use electroencephalography (EEG) signals and a multilayer network formalism to model cross-frequency coupling across the brain and identify brain regions that operate as connecting hubs, thus facilitating inte-grative function. To this goal, we investigate α-γ coupling and θ-γ coupling during exposure to olfactory stimuli of different pleasantness levels. We found that a wider distributed network of hubs emerges in the higher pleasantness condition and that significant differences in the hub connectivity are located in the middle frontal and central regions. Our results indicate the consistent functional role that γ band activity plays in information integration in olfactory perception.

Working Memory and Cross-Frequency Coupling of Neuronal Oscillations

Working memory (WM) is the active retention and processing of information over a few seconds and is considered an essential component of cognitive function. The reduced WM capacity is a common feature in many diseases, such as schizophrenia, attention deficit hyperactivity disorder (ADHD), mild cognitive impairment (MCI), and Alzheimer's disease (AD). The theta-gamma neural code is an essential component of memory representations in the multi-item WM. A large body of studies have examined the association between cross-frequency coupling (CFC) across the cerebral cortices and WM performance; electrophysiological data together with the behavioral results showed the associations between CFC and WM performance. The oscillatory entrainment (sensory, non-invasive electrical/magnetic, and invasive electrical) remains the key method to investigate the causal relationship between CFC and WM. The frequency-tuned non-invasive brain stimulation is a promising way to improve WM performance in healthy and non-healthy patients with cognitive impairment. The WM performance is sensitive to the phase and rhythm of externally applied stimulations. CFC-transcranial-alternating current stimulation (CFC-tACS) is a recent approach in neuroscience that could alter cognitive outcomes. The studies that investigated (1) the association between CFC and WM and (2) the brain stimulation protocols that enhanced WM through modulating CFC by the means of the non-invasive brain stimulation techniques have been included in this review. In principle, this review can guide the researchers to identify the most prominent form of CFC associated with WM processing (e.g., theta/gamma phase-amplitude coupling), and to define the previously published studies that manipulate endogenous CFC externally to improve WM. This in turn will pave the path for future studies aimed at investigating the CFC-tACS effect on WM. The CFC-tACS protocols need to be thoroughly studied before they can be considered as therapeutic tools in patients with WM deficits.

Mental Arithmetic Task Recognition Using Effective Connectivity and Hierarchical Feature Selection From EEG Signals

Introduction

Mental arithmetic analysis based on Electroencephalogram (EEG) signals can help understand disorders, such as attention-deficit hyperactivity, dyscalculia, or autism spectrum disorder where the difficulty in learning or understanding the arithmetic exists. Most mental arithmetic recognition systems rely on features of a single channel of EEG; however, the relationships between EEG channels in the form of effective brain connectivity analysis can contain valuable information. This study aims to find distinctive, effective brain connectivity features and create a hierarchical feature selection for effectively classifying mental arithmetic and baseline tasks.

Methods

We estimated effective connectivity using Directed Transfer Function (DTF), direct DTF (dDTF) and Generalized Partial Directed Coherence (GPDC) methods. These measures determine the causal relationship between different brain areas. A hierarchical feature subset selection method selects the most significant effective connectivity features. Initially, Kruskal- Wallis test was performed. Consequently, five feature selection algorithms, namely, Support Vector Machine (SVM) method based on Recursive Feature Elimination, Fisher score, mutual information, minimum Redundancy Maximum Relevance (RMR), and concave minimization and SVM are used to select the best discriminative features. Finally, the SVM method was used for classification.

Results

The obtained results indicated that the best EEG classification performance in 29 participants and 60 trials is obtained using GPDC and feature selection via concave minimization method in Beta2 (15-22Hz) frequency band with 89% accuracy.

Conclusion

This new hierarchical automated system could be helpful in the discrimination of mental arithmetic and baseline tasks from EEG signals effectively.

Highlights

Propose effective connectivity to describe EEG signals during mental arithmetic task.Most significant connectivity features from generalized partial directed coherence method.Hierarchical feature selection from Kruskal-Wallis test and concave minimization method.

Plain Language Summary

Brain analysis methods by Electroencephalogram (EEG) signals provide a suitable method to monitor human brain activity due to having high temporal resolution, being noninvasive, inexpensive, and portable method. Analysis of mental arithmetic based EEG signal is helpful for psychological disorders like dyscalculia where they have learning understanding arithmetic, attention deficit hyperactivity, and autism spectrum disorders with attention deficit problem. This study finds distinctive effective brain connectivity features and creates a hierarchical feature selection for classification of mental arithmetic and baseline tasks effectively. Best EEG classification performance in 29 participants and 60 trials is obtained using Generalized Partial Directed Coherence (GPDC) methods and feature selection via concave minimization method in Beta2 (15-22Hz) frequency band with 89% accuracy. Thus, this new hierarchical automated system is useful for discrimination of mental arithmetic and baseline tasks from EEG signal effectively.

Propofol Anesthesia Decreased the Efficiency of Long-Range Cortical Interaction in Humans

OBJECTIVE

Phase-amplitude coupling (PAC) has recently been used to illuminate cross-frequency coordination in neurophysiological activity of electroencephalogram. However, the PAC at a meso-scale (electrocorticogram, ECoG) and PAC between different areas have still not been fully clarified.

METHODS

In this study, we analyzed ECoG data recorded from surgical patients (n = 9) with pharmaco-resistant epilepsy during the surgical treatment. The modulogram and a genuine modulation index, based on a Kullback-Leibler distance and permutation test method, were developed and used to measure the slow oscillation (SO) (0.15-1 Hz)-α (8-13 Hz) PAC of within-lead and cross-lead during transitions from states of wakefulness to unconsciousness during propofol induced general anesthesia.

RESULTS

In within-lead SO-α PAC, the modulation index increased in the unconscious state (p < 0.05, Tukey's test), the percentages of genuine modulation indices also increased in the unconscious state (p < 0.001 in the frontal area and p < 0.01 in the parietal area), and distinct PAC patterns emerged more often. In cross-lead SO-α PAC, there are fewer PAC patterns compared to within-lead, and the percentages of genuine modulation indices decreased significantly (p < 0.001).

CONCLUSION

The increased modulation index of within-lead and cross-lead SO-α PAC is associated with a reduction of information integration and the efficiency of long distance synchronization. These findings demonstrate that the propofol causes the neuronal populations to enter a 'busy' state in a local scale, which prevents the information integration in long-range areas.

EEG correlation during the solving of simple and complex logical-mathematical problems

Solving logical-mathematical word problems is a complex task that requires numerous cognitive operations, including comprehension, reasoning, and calculation. These abilities have been associated with activation of the parietal, temporal, and prefrontal cortices. It has been suggested that the reasoning involved in solving logical-mathematical problems requires the coordinated functionality of all these cortical areas. In this study was evaluated the activation and electroencephalographic (EEG) correlation of the prefrontal, temporal, and parietal regions in young men while solving logical-mathematical word problems with two degrees of difficulty: simple and complex. During the solving of complex problems, higher absolute power and EEG correlation of the alpha and fast bands between the left frontal and parietal cortices were observed. A temporal deactivation and functional decoupling of the right parietal-temporal cortices also were obtained. Solving complex problems probably require activation of a left prefrontal-parietal circuit to maintain and manipulate multiple pieces of information. The temporal deactivation and decreased parietal-temporal correlation could be associated to text processing and suppression of the content-dependent reasoning to focus cognitive resources on the mathematical reasoning. Together, these findings support a pivotal role for the left prefrontal and parietal cortices in mathematical reasoning and of the temporal regions in text processing required to understand and solve written mathematical problems.

Neural Processing Mechanism of Mental Calculation Based on Cerebral Oscillatory Changes: A Comparison Between Abacus Experts and Novices

Background: Abacus experts could mentally calculate fast some mathematical operations using multi-digit numbers. The temporal dynamics of abacus mental calculation are still unknown although some behavioral and neuroimaging studies have suggested a visuospatial and visuomotor neural process during abacus mental calculation. Therefore, this contribution aims to clarify the significant similarities and the differences between experts and novices by investigating calculation-induced neuromagnetic responses based on cerebral oscillatory changes.

Methods: Twelve to 13 healthy abacus experts and 17 non-experts participated in two experimental paradigms using non-invasive neuromagnetic measurements. In experiments 1 and 2, the spatial distribution of oscillatory changes presented mental calculations and temporal frequency profiles during addition while examining multiplication tasks. The MEG data were analyzed using synthetic aperture magnetometry (SAM) with an adaptive beamformer to calculate the group average of the spatial distribution of oscillatory changes and their temporal frequency profiles in source-level analyses.

Results: Using a group average of the spatial distribution of oscillatory changes, we observed some common brain activities in both right-handed abacus experts and non-experts. In non-experts, we detected the right dorsolateral prefrontal cortex (DLPFC) and bilateral Intraparietal sulcus (IPS); whereas in experts, detected the bilateral parieto-occipital sulcus (POS), right inferior frontal gyrus (IFG), and left sensorimotor areas mainly. Based on the findings generated, we could propose calculation processing models for both abacus experts and non- experts conveniently.

Conclusion: The proposed model of calculation processing in abacus experts and novices revealed that the novices could calculate logically depending on numerical processing in the left IPS. In contrast, abacus experts are utilizing spatial processing using a memorized imaginary abacus, which distributed over the bilateral hemispheres in the IFG and sensorimotor areas.

From thoughtless awareness to effortful cognition: alpha - theta cross-frequency dynamics in experienced meditators during meditation, rest and arithmetic

Neural activity is known to oscillate within discrete frequency bands and the synchronization between these rhythms is hypothesized to underlie information integration in the brain. Since strict synchronization is only possible for harmonic frequencies, a recent theory proposes that the interaction between different brain rhythms is facilitated by transient harmonic frequency arrangements. In this line, it has been recently shown that the transient occurrence of 2:1 harmonic cross-frequency relationships between alpha and theta rhythms (i.e. f_alpha≈ 12 Hz; f_theta≈ 6 Hz) is enhanced during effortful cognition. In this study, we tested whether achieving a state of ‘mental emptiness’ during meditation is accompanied by a relative decrease in the occurrence of 2:1 harmonic cross-frequency relationships between alpha and theta rhythms. Continuous EEG recordings (19 electrodes) were obtained from 43 highly experienced meditators during meditation practice, rest and an arithmetic task. We show that the occurrence of transient alpha:theta 2:1 harmonic relationships increased linearly from a meditative to an active cognitive processing state (i.e. meditation < rest < arithmetic task). It is argued that transient EEG cross-frequency arrangements that prevent alpha:theta cross-frequency coupling could facilitate the experience of ‘mental emptiness’ by avoiding the interaction between the memory and executive components of cognition.

A brain connectivity characterization of children with different levels of mathematical achievement based on graph metrics

Recent studies aiming to facilitate mathematical skill development in primary school children have explored the electrophysiological characteristics associated with different levels of arithmetic achievement. The present work introduces an alternative EEG signal characterization using graph metrics and, based on such features, a classification analysis using a decision tree model. This proposal aims to identify group differences in brain connectivity networks with respect to mathematical skills in elementary school children. The methods of analysis utilized were signal-processing (EEG artifact removal, Laplacian filtering, and magnitude square coherence measurement) and the characterization (Graph metrics) and classification (Decision Tree) of EEG signals recorded during performance of a numerical comparison task. Our results suggest that the analysis of quantitative EEG frequency-band parameters can be used successfully to discriminate several levels of arithmetic achievement. Specifically, the most significant results showed an accuracy of 80.00% (α band), 78.33% (δ band), and 76.67% (θ band) in differentiating high-skilled participants from low-skilled ones, averaged-skilled subjects from all others, and averaged-skilled participants from low-skilled ones, respectively. The use of a decision tree tool during the classification stage allows the identification of several brain areas that seem to be more specialized in numerical processing.

Typical and Aberrant Functional Brain Flexibility: Lifespan Development and Aberrant Organization in Traumatic Brain Injury and Dyslexia

Intrinsic functional connectivity networks derived from different neuroimaging methods and connectivity estimators have revealed robust developmental trends linked to behavioural and cognitive maturation. The present study employed a dynamic functional connectivity approach to determine dominant intrinsic coupling modes in resting-state neuromagnetic data from 178 healthy participants aged 8-60 years. Results revealed significant developmental trends in three types of dominant intra- and inter-hemispheric neuronal population interactions (amplitude envelope, phase coupling, and phase-amplitude synchronization) involving frontal, temporal, and parieto-occipital regions. Multi-class support vector machines achieved 89% correct classification of participants according to their chronological age using dynamic functional connectivity indices. Moreover, systematic temporal variability in functional connectivity profiles, which was used to empirically derive a composite flexibility index, displayed an inverse U-shaped curve among healthy participants. Lower flexibility values were found among age-matched children with reading disability and adults who had suffered mild traumatic brain injury. The importance of these results for normal and abnormal brain development are discussed in light of the recently proposed role of cross-frequency interactions in the fine-grained coordination of neuronal population activity.

Detrended Fluctuation, Coherence, and Spectral Power Analysis of Activation Rearrangement in EEG Dynamics During Cognitive Workload

In the study of human cognitive activity using electroencephalogram (EEG), the brain dynamics parameters and characteristics play a crucial role. They allow to investigate the changes in functionality depending on the environment and task performance process, and also to access the intensity of the brain activity in various locations of the cortex and its dependencies. Usually, the dynamics of activation of different brain areas during the cognitive tasks are being studied by spectral analysis based on power spectral density (PSD) estimation, and coherence analysis, which are de facto standard tools in quantitative characterization of brain activity. PSD and coherence reflect the strength of oscillations and similarity of the emergence of these oscillations in the brain, respectively, while the concept of stability of brain activity over time is not well defined and less formalized. We propose to employ the detrended fluctuation analysis (DFA) as a measure of the EEG persistence over time, and use the DFA scaling exponent as its quantitative characteristics. We applied DFA to the study of the changes in activation in brain dynamics during mental calculations and united it with PSD and coherence estimation. In the experiment, EEGs during resting state and mental serial subtraction from 36 subjects were recorded and analyzed in four frequency ranges: θ1 (4.1-5.8 Hz), θ2 (5.9-7.4 Hz), β1 (13-19.9 Hz), and β2 (20-25 Hz). PSD maps to access the intensity of cortex activation and coherence to quantify the connections between different brain areas were calculated, the distribution of DFA scaling exponent over the head surface was exploited to measure the time characteristics of the dynamics of brain activity. Obtained arrangements of DFA scaling exponent suggest that normal functioning of the brain is characterized by long-term temporal correlations in the cortex. Topographical distribution of the DFA scaling exponent was comparable for θ and β frequency bands, demonstrating the largest values of DFA scaling exponent during cognitive activation. The study shows that the long-term temporal correlations evaluated by DFA can be of great interest for diagnosis of the variety of brain dysfunctions of different etiology in the future.

Tracking Transient Changes in the Neural Frequency Architecture: Harmonic Relationships between Theta and Alpha Peaks Facilitate Cognitive Performance

The synchronization between neural oscillations at different frequencies has been proposed as a core mechanism for the coordination and integration of neural systems at different spatiotemporal scales. Because neural oscillations of different frequencies can only fully synchronize when their "peak" frequencies form harmonic relationships (e.g., f ₂ = f ₁/2), the present study explored whether the transient occurrence of harmonic cross-frequency relationship between task-relevant rhythms underlies efficient cognitive processing. Continuous EEG recordings (51 human participants; 14 males) were obtained during an arithmetic task, rest and breath focus. In two separate experiments, we consistently show that the proportion of epochs displaying a 2:1 harmonic relationship between alpha (8-14 Hz) and theta (4-8 Hz) peak frequencies (i.e., alpha_peak ≈ 10.6 Hz; theta_peak ≈ 5.3 Hz), was significantly higher when cognitive demands increased. In addition, a higher incidence of 2:1 harmonic cross-frequency relationships was significantly associated with increased alpha-theta phase synchrony and improved arithmetic task performance, thereby underlining the functional relevance of this cross-frequency configuration. Notably, opposite dynamics were identified for a specific range of "nonharmonic" alpha-theta cross-frequency relationships (i.e., alpha_peak/theta_peak = 1.1-1.6), which showed a higher incidence during rest compared with the arithmetic task. The observation that alpha and theta rhythms shifted into harmonic versus nonharmonic cross-frequency relationships depending on (cognitive) task demands is in line with the notion that the neural frequency architecture entails optimal frequency arrangements to facilitate cross-frequency "coupling" and "decoupling".SIGNIFICANCE STATEMENT Neural activity is known to oscillate within discrete frequency bands and the interplay between these brain rhythms is hypothesized to underlie cognitive functions. A recent theory posits that shifts in the peak frequencies of oscillatory rhythms form the principal mechanism by which cross-frequency coupling and decoupling is implemented in the brain. In line with this notion, we show that the occurrence of a cross-frequency arrangement that mathematically enables coupling between alpha and theta rhythms is more prominent during active cognitive processing (compared with rest and non-cognitively demanding tasks) and is associated with improved cognitive performance. Together, our results open new vistas for future research on cross-frequency dynamics in the brain and their functional role in cognitive processing.

Data-Driven Topological Filtering Based on Orthogonal Minimal Spanning Trees: Application to Multigroup Magnetoencephalography Resting-State Connectivity

In the present study, a novel data-driven topological filtering technique is introduced to derive the backbone of functional brain networks relying on orthogonal minimal spanning trees (OMSTs). The method aims to identify the essential functional connections to ensure optimal information flow via the objective criterion of global efficiency minus the cost of surviving connections. The OMST technique was applied to multichannel, resting-state neuromagnetic recordings from four groups of participants: healthy adults (n = 50), adults who have suffered mild traumatic brain injury (n = 30), typically developing children (n = 27), and reading-disabled children (n = 25). Weighted interactions between network nodes (sensors) were computed using an integrated approach of dominant intrinsic coupling modes based on two alternative metrics (symbolic mutual information and phase lag index), resulting in excellent discrimination of individual cases according to their group membership. Classification results using OMST-derived functional networks were clearly superior to results using either relative power spectrum features or functional networks derived through the conventional minimal spanning tree algorithm.

How to Build a Functional Connectomic Biomarker for Mild Cognitive Impairment From Source Reconstructed MEG Resting-State Activity: The Combination of ROI Representation and Connectivity Estimator Matters

Our work aimed to demonstrate the combination of machine learning and graph theory for the designing of a connectomic biomarker for mild cognitive impairment (MCI) subjects using eyes-closed neuromagnetic recordings. The whole analysis based on source-reconstructed neuromagnetic activity. As ROI representation, we employed the principal component analysis (PCA) and centroid approaches. As representative bi-variate connectivity estimators for the estimation of intra and cross-frequency interactions, we adopted the phase locking value (PLV), the imaginary part (iPLV) and the correlation of the envelope (CorrEnv). Both intra and cross-frequency interactions (CFC) have been estimated with the three connectivity estimators within the seven frequency bands (intra-frequency) and in pairs (CFC), correspondingly. We demonstrated how different versions of functional connectivity graphs single-layer (SL-FCG) and multi-layer (ML-FCG) can give us a different view of the functional interactions across the brain areas. Finally, we applied machine learning techniques with main scope to build a reliable connectomic biomarker by analyzing both SL-FCG and ML-FCG in two different options: as a whole unit using a tensorial extraction algorithm and as single pair-wise coupling estimations. We concluded that edge-weighed feature selection strategy outperformed the tensorial treatment of SL-FCG and ML-FCG. The highest classification performance was obtained with the centroid ROI representation and edge-weighted analysis of the SL-FCG reaching the 98% for the CorrEnv in α1:α2 and 94% for the iPLV in α2. Classification performance based on the multi-layer participation coefficient, a multiplexity index reached 52% for iPLV and 52% for CorrEnv. Selected functional connections that build the multivariate connectomic biomarker in the edge-weighted scenario are located in default-mode, fronto-parietal, and cingulo-opercular network. Our analysis supports the notion of analyzing FCG simultaneously in intra and cross-frequency whole brain interactions with various connectivity estimators in beamformed recordings.

Improving the Reliability of Network Metrics in Structural Brain Networks by Integrating Different Network Weighting Strategies into a Single Graph

Structural brain networks estimated from diffusion MRI (dMRI) via tractography have been widely studied in healthy controls and patients with neurological and psychiatric diseases. However, few studies have addressed the reliability of derived network metrics both node-specific and network-wide. Different network weighting strategies (NWS) can be adopted to weight the strength of connection between two nodes yielding structural brain networks that are almost fully-weighted. Here, we scanned five healthy participants five times each, using a diffusion-weighted MRI protocol and computed edges between 90 regions of interest (ROI) from the Automated Anatomical Labeling (AAL) template. The edges were weighted according to nine different methods. We propose a linear combination of these nine NWS into a single graph using an appropriate diffusion distance metric. We refer to the resulting weighted graph as an Integrated Weighted Structural Brain Network (ISWBN). Additionally, we consider a topological filtering scheme that maximizes the information flow in the brain network under the constraint of the overall cost of the surviving connections. We compared each of the nine NWS and the ISWBN based on the improvement of: (a) intra-class correlation coefficient (ICC) of well-known network metrics, both node-wise and per network level; and (b) the recognition accuracy of each subject compared to the remainder of the cohort, as an attempt to access the uniqueness of the structural brain network for each subject, after first applying our proposed topological filtering scheme. Based on a threshold where the network level ICC should be >0.90, our findings revealed that six out of nine NWS lead to unreliable results at the network level, while all nine NWS were unreliable at the node level. In comparison, our proposed ISWBN performed as well as the best performing individual NWS at the network level, and the ICC was higher compared to all individual NWS at the node level. Importantly, both network and node-wise ICCs of network metrics derived from the topologically filtered ISBWN (ISWBN^TF), were further improved compared to the non-filtered ISWBN. Finally, in the recognition accuracy tests, we assigned each single ISWBN^TF to the correct subject. We also applied our methodology to a second dataset of diffusion-weighted MRI in healthy controls and individuals with psychotic experience. Following a binary classification scheme, the classification performance based on ISWBN^TF outperformed the nine different weighting strategies and the ISWBN. Overall, these findings suggest that the proposed methodology results in improved characterization of genuine between-subject differences in connectivity leading to the possibility of network-based structural phenotyping.

Mining Time-Resolved Functional Brain Graphs to an EEG-Based Chronnectomic Brain Aged Index (CBAI)

The brain at rest consists of spatially and temporal distributed but functionally connected regions that called intrinsic connectivity networks (ICNs). Resting state electroencephalography (rs-EEG) is a way to characterize brain networks without confounds associated with task EEG such as task difficulty and performance. A novel framework of how to study dynamic functional connectivity under the notion of functional connectivity microstates (FCμstates) and symbolic dynamics is further discussed. Furthermore, we introduced a way to construct a single integrated dynamic functional connectivity graph (IDFCG) that preserves both the strength of the connections between every pair of sensors but also the type of dominant intrinsic coupling modes (DICM). The whole methodology is demonstrated in a significant and unexplored task for EEG which is the definition of an objective Chronnectomic Brain Aged index (CBAI) extracted from resting-state data (N = 94 subjects) with both eyes-open and eyes-closed conditions. Novel features have been defined based on symbolic dynamics and the notion of DICM and FCμstates. The transition rate of FCμstates, the symbolic dynamics based on the evolution of FCμstates (the Markovian Entropy, the complexity index), the probability distribution of DICM, the novel Flexibility Index that captures the dynamic reconfiguration of DICM per pair of EEG sensors and the relative signal power constitute a valuable pool of features that can build the proposed CBAI. Here we applied a feature selection technique and Extreme Learning Machine (ELM) classifier to discriminate young adults from middle-aged and a Support Vector Regressor to build a linear model of the actual age based on EEG-based spatio-temporal features. The most significant type of features for both prediction of age and discrimination of young vs. adults age groups was the dynamic reconfiguration of dominant coupling modes derived from a subset of EEG sensor pairs. Specifically, our results revealed a very high prediction of age for eyes-open (R² = 0.60; y = 0.79x + 8.03) and lower for eyes-closed (R² = 0.48; y = 0.71x + 10.91) while we succeeded to correctly classify young vs. middle-age group with 97.8% accuracy in eyes-open and 87.2% for eyes-closed. Our results were reproduced also in a second dataset for further external validation of the whole analysis. The proposed methodology proved valuable for the characterization of the intrinsic properties of dynamic functional connectivity through the age untangling developmental differences using EEG resting-state recordings.