Spectral analysis for nonstationary and nonlinear systems: a discrete-time-model-based approach

Dementia classification using a graph neural network on imaging of effective brain connectivity

Alzheimer's disease (AD) and Parkinson's disease (PD) are two of the most common forms of neurodegenerative diseases. The literature suggests that effective brain connectivity (EBC) has the potential to track differences between AD, PD and healthy controls (HC). However, how to effectively use EBC estimations for the research of disease diagnosis remains an open problem. To deal with complex brain networks, graph neural network (GNN) has been increasingly popular in very recent years and the effectiveness of combining EBC and GNN techniques has been unexplored in the field of dementia diagnosis. In this study, a novel directed structure learning GNN (DSL-GNN) was developed and performed on the imaging of EBC estimations and power spectrum density (PSD) features. In comparison to the previous studies on GNN, our proposed approach enhanced the functionality for processing directional information, which builds the basis for more efficiently performing GNN on EBC. Another contribution of this study is the creation of a new framework for applying univariate and multivariate features simultaneously in a classification task. The proposed framework and DSL-GNN are validated in four discrimination tasks and our approach exhibited the best performance, against the existing methods, with the highest accuracy of 94.0% (AD vs. HC), 94.2% (PD vs. HC), 97.4% (AD vs. PD) and 93.0% (AD vs. PD vs. HC). In a word, this research provides a robust analytical framework to deal with complex brain networks containing causal directional information and implies promising potential in the diagnosis of two of the most common neurodegenerative conditions.

Brain functional and effective connectivity based on electroencephalography recordings: A review

Functional connectivity and effective connectivity of the human brain, representing statistical dependence and directed information flow between cortical regions, significantly contribute to the study of the intrinsic brain network and its functional mechanism. Many recent studies on electroencephalography (EEG) have been focusing on modeling and estimating brain connectivity due to increasing evidence that it can help better understand various brain neurological conditions. However, there is a lack of a comprehensive updated review on studies of EEG-based brain connectivity, particularly on visualization options and associated machine learning applications, aiming to translate those techniques into useful clinical tools. This article reviews EEG-based functional and effective connectivity studies undertaken over the last few years, in terms of estimation, visualization, and applications associated with machine learning classifiers. Methods are explored and discussed from various dimensions, such as either linear or nonlinear, parametric or nonparametric, time-based, and frequency-based or time-frequency-based. Then it is followed by a novel review of brain connectivity visualization methods, grouped by Heat Map, data statistics, and Head Map, aiming to explore the variation of connectivity across different brain regions. Finally, the current challenges of related research and a roadmap for future related research are presented.

Nonlinear System Identification of Neural Systems from Neurophysiological Signals

The human nervous system is one of the most complicated systems in nature. Complex nonlinear behaviours have been shown from the single neuron level to the system level. For decades, linear connectivity analysis methods, such as correlation, coherence and Granger causality, have been extensively used to assess the neural connectivities and input-output interconnections in neural systems. Recent studies indicate that these linear methods can only capture a certain amount of neural activities and functional relationships, and therefore cannot describe neural behaviours in a precise or complete way. In this review, we highlight recent advances in nonlinear system identification of neural systems, corresponding time and frequency domain analysis, and novel neural connectivity measures based on nonlinear system identification techniques. We argue that nonlinear modelling and analysis are necessary to study neuronal processing and signal transfer in neural systems quantitatively. These approaches can hopefully provide new insights to advance our understanding of neurophysiological mechanisms underlying neural functions. These nonlinear approaches also have the potential to produce sensitive biomarkers to facilitate the development of precision diagnostic tools for evaluating neurological disorders and the effects of targeted intervention.

Detection of time-, frequency- and direction-resolved communication within brain networks

Electroencephalography (EEG) records fast-changing neuronal signalling and communication and thus can offer a deep understanding of cognitive processes. However, traditional data analyses which employ the Fast-Fourier Transform (FFT) have been of limited use as they do not allow time- and frequency-resolved tracking of brain activity and detection of directional connectivity. Here, we applied advanced qEEG tools using autoregressive (AR) modelling, alongside traditional approaches, to murine data sets from common research scenarios: (a) the effect of age on resting EEG; (b) drug actions on non-rapid eye movement (NREM) sleep EEG (pharmaco-EEG); and (c) dynamic EEG profiles during correct vs incorrect spontaneous alternation responses in the Y-maze. AR analyses of short data strips reliably detected age- and drug-induced spectral EEG changes, while renormalized partial directed coherence (rPDC) reported direction- and time-resolved connectivity dynamics in mice. Our approach allows for the first time inference of behaviour- and stage-dependent data in a time- and frequency-resolved manner, and offers insights into brain networks that underlie working memory processing beyond what can be achieved with traditional methods.

Estimation of Time-Varying, Intrinsic and Reflex Dynamic Joint Stiffness during Movement. Application to the Ankle Joint

Dynamic joint stiffness determines the relation between joint position and torque, and plays a vital role in the control of posture and movement. Dynamic joint stiffness can be quantified during quasi-stationary conditions using disturbance experiments, where small position perturbations are applied to the joint and the torque response is recorded. Dynamic joint stiffness is composed of intrinsic and reflex mechanisms that act and change together, so that nonlinear, mathematical models and specialized system identification techniques are necessary to estimate their relative contributions to overall joint stiffness. Quasi-stationary experiments have demonstrated that dynamic joint stiffness is heavily modulated by joint position and voluntary torque. Consequently, during movement, when joint position and torque change rapidly, dynamic joint stiffness will be Time-Varying (TV). This paper introduces a new method to quantify the TV intrinsic and reflex components of dynamic joint stiffness during movement. The algorithm combines ensemble and deterministic approaches for estimation of TV systems; and uses a TV, parallel-cascade, nonlinear system identification technique to separate overall dynamic joint stiffness into intrinsic and reflex components from position and torque records. Simulation studies of a stiffness model, whose parameters varied with time as is expected during walking, demonstrated that the new algorithm accurately tracked the changes in dynamic joint stiffness using as little as 40 gait cycles. The method was also used to estimate the intrinsic and reflex dynamic ankle stiffness from an experiment with a healthy subject during which ankle movements were imposed while the subject maintained a constant muscle contraction. The method identified TV stiffness model parameters that predicted the measured torque very well, accounting for more than 95% of its variance. Moreover, both intrinsic and reflex dynamic stiffness were heavily modulated through the movement in a manner that could not be predicted from quasi-stationary experiments. The new method provides the tool needed to explore the role of dynamic stiffness in the control of movement.

Nonlinear interactions in the thalamocortical loop in essential tremor: A model-based frequency domain analysis

There is increasing evidence to suggest that essential tremor has a central origin. Different structures appear to be part of the central tremorogenic network, including the motor cortex, the thalamus and the cerebellum. Some studies using electroencephalogram (EEG) and magnetoencephalography (MEG) show linear association in the tremor frequency between the motor cortex and the contralateral tremor electromyography (EMG). Additionally, high thalamomuscular coherence is found with the use of thalamic local field potential (LFP) recordings and tremulous EMG in patients undergoing surgery for deep brain stimulation (DBS). Despite a well-established reciprocal anatomical connection between the thalamus and cortex, the functional association between the two structures during "tremor-on" periods remains elusive. Thalamic (Vim) LFPs, ipsilateral scalp EEG from the sensorimotor cortex and contralateral tremor arm EMG recordings were obtained from two patients with essential tremor who had undergone successful surgery for DBS. Coherence analysis shows a strong linear association between thalamic LFPs and contralateral tremor EMG, but the relationship between the EEG and the thalamus is much less clear. These measurements were then analyzed by constructing a novel parametric nonlinear autoregressive with exogenous input (NARX) model. This new approach uncovered two distinct and not overlapping frequency "channels" of communication between Vim thalamus and the ipsilateral motor cortex, defining robustly "tremor-on" versus "tremor-off" states. The associated estimated nonlinear time lags also showed non-overlapping values between the two states, with longer corticothalamic lags (exceeding 50ms) in the tremor active state, suggesting involvement of an indirect multisynaptic loop. The results reveal the importance of the nonlinear interactions between cortical and subcortical areas in the central motor network of essential tremor. This work is important because it demonstrates for the first time that in essential tremor the functional interrelationships between the cortex and thalamus should not be sought exclusively within individual frequencies but more importantly between cross-frequency nonlinear interactions. Should our results be successfully reproduced on a bigger cohort of patients with essential tremor, our approach could be used to create an on-demand closed-loop DBS device, able to automatically activate when the tremor is on.

A nonlinear generalization of spectral Granger causality

Spectral measures of linear Granger causality have been widely applied to study the causal connectivity between time series data in neuroscience, biology, and economics. Traditional Granger causality measures are based on linear autoregressive with exogenous (ARX) inputs models of time series data, which cannot truly reveal nonlinear effects in the data especially in the frequency domain. In this study, it is shown that the classical Geweke's spectral causality measure can be explicitly linked with the output spectra of corresponding restricted and unrestricted time-domain models. The latter representation is then generalized to nonlinear bivariate signals and for the first time nonlinear causality analysis in the frequency domain. This is achieved by using the nonlinear ARX (NARX) modeling of signals, and decomposition of the recently defined output frequency response function which is related to the NARX model.

Assessing dynamic spectral causality by lagged adaptive directed transfer function and instantaneous effect factor

It is of significance to assess the dynamic spectral causality among physiological signals. Several practical estimators adapted from spectral Granger causality have been exploited to track dynamic causality based on the framework of time-varying multivariate autoregressive (tvMVAR) models. The nonzero covariance of the model's residuals has been used to describe the instantaneous effect phenomenon in some causality estimators. However, for the situations with Gaussian residuals in some autoregressive models, it is challenging to distinguish the directed instantaneous causality if the sufficient prior information about the "causal ordering" is missing. Here, we propose a new algorithm to assess the time-varying causal ordering of tvMVAR model under the assumption that the signals follow the same acyclic causal ordering for all time lags and to estimate the instantaneous effect factor (IEF) value in order to track the dynamic directed instantaneous connectivity. The time-lagged adaptive directed transfer function (ADTF) is also estimated to assess the lagged causality after removing the instantaneous effect. In this study, we first investigated the performance of the causal-ordering estimation algorithm and the accuracy of IEF value. Then, we presented the results of IEF and time-lagged ADTF method by comparing with the conventional ADTF method through simulations of various propagation models. Statistical analysis results suggest that the new algorithm could accurately estimate the causal ordering and give a good estimation of the IEF values in the Gaussian residual conditions. Meanwhile, the time-lagged ADTF approach is also more accurate in estimating the time-lagged dynamic interactions in a complex nervous system after extracting the instantaneous effect. In addition to the simulation studies, we applied the proposed method to estimate the dynamic spectral causality on real visual evoked potential (VEP) data in a human subject. Its usefulness in time-variant spectral causality assessment was demonstrated through the mutual causality investigation of brain activity during the VEP experiments.