Search

Spatial source phase, the phase information of spatial maps extracted from functional magnetic resonance imaging (fMRI) data by data-driven methods such as independent component analysis (ICA), has rarely been studied. While the observed phase has been shown to convey unique brain information, the role of spatial source phase in representing the intrinsic activity of the brain is yet not clear. This study explores the spatial source phase for identifying spatial differences between patients with schizophrenia (SZs) and healthy controls (HCs) using complex-valued resting-state fMRI data from 82 individuals. ICA is first applied to preprocess fMRI data, and post-ICA phase de-ambiguity and denoising are then performed. The ability of spatial source phase to characterize spatial differences is examined by the homogeneity of variance test (voxel-wise F-test) with false discovery rate correction. Resampling techniques are performed to ensure that the observations are significant and reliable. We focus on two components of interest widely used in analyzing SZs, including the default mode network (DMN) and auditory cortex. Results show that the spatial source phase exhibits more significant variance changes and higher sensitivity to the spatial differences between SZs and HCs in the anterior areas of DMN and the left auditory cortex, compared to the magnitude of spatial activations. Our findings show that the spatial source phase can potentially serve as a new brain imaging biomarker and provide a novel perspective on differences in SZs compared to HCs, consistent with but extending previous work showing increased variability in patient data.

Previous researches demonstrate that major depression disorder (MDD) is associated with widespread network dysconnectivity, and the dynamics of functional connectivity networks are important to delineate the neural mechanisms of MDD. Cortical electroencephalography (EEG) oscillations act as coordinators to connect different brain regions, and various assemblies of oscillations can form different networks to support different cognitive tasks. Studies have demonstrated that the dysconnectivity of EEG oscillatory networks is related with MDD. In this study, we investigated the oscillatory hyperconnectivity and hypoconnectivity networks in MDD under a naturalistic and continuous stimuli condition of music listening. With the assumption that the healthy group and the MDD group share similar brain topology from the same stimuli and also retain individual brain topology for group differences, we applied the coupled nonnegative tensor decomposition algorithm on two adjacency tensors with the dimension of time × frequency × connectivity × subject, and imposed double-coupled constraints on spatial and spectral modes. The music-induced oscillatory networks were identified by a correlation analysis approach based on the permutation test between extracted temporal factors and musical features. We obtained three hyperconnectivity networks from the individual features of MDD and three hypoconnectivity networks from common features. The results demonstrated that the dysfunction of oscillation-modulated networks could affect the involvement in music perception for MDD patients. Those oscillatory dysconnectivity networks may provide promising references to reveal the pathoconnectomics of MDD and potential biomarkers for the diagnosis of MDD.

Canonical polyadic decomposition (CPD) of multi-subject complex-valued fMRI data can be used to provide spatially and temporally shared components among groups with both magnitude and phase information. However, the CPD model is not well formulated due to the large subject variability in the spatial and temporal modalities, as well as the high noise level in complex-valued fMRI data. Considering that the shift-invariant CPD can model temporal variability across subjects, we propose to further impose a phase sparsity constraint on the shared spatial maps to denoise the complex-valued components and to model the inter-subject spatial variability as well. More precisely, subject-specific time delays are first estimated for the complex-valued shared time courses in the framework of real-valued shift-invariant CPD. Source phase sparsity is then imposed on the complex-valued shared spatial maps. A smoothed $\ell _{\mathbf {{0}}}$ norm is specifically used to reduce voxels with large phase values after phase de-ambiguity based on the small phase characteristic of BOLD-related voxels. The results from both the simulated and experimental fMRI data demonstrate improvements of the proposed method over three complex-valued algorithms, namely, tensor-based spatial ICA, shift-invariant CPD and CPD without spatiotemporal constraints. When comparing with a real-valued algorithm combining shift-invariant CPD and ICA, the proposed method detects 178.7% more contiguous task-related activations. peerReviewed

Background Ongoing EEG data are recorded as mixtures of stimulus-elicited EEG, spontaneous EEG and noises, which require advanced signal processing techniques for separation and analysis. Existing methods cannot simultaneously consider common and individual characteristics among/within subjects when extracting stimulus-elicited brain activities from ongoing EEG elicited by 512-s long modern tango music. New method Aiming to discover the commonly music-elicited brain activities among subjects, we provide a comprehensive framework based on fast double-coupled nonnegative tensor decomposition (FDC-NTD) algorithm. The proposed algorithm with a generalized model is capable of simultaneously decomposing EEG tensors into common and individual components. Results With the proposed framework, the brain activities can be effectively extracted and sorted into the clusters of interest. The proposed algorithm based on the generalized model achieved higher fittings and stronger robustness. In addition to the distribution of centro-parietal and occipito-parietal regions with theta and alpha oscillations, the music-elicited brain activities were also located in the frontal region and distributed in the 4∼11 Hz band. Comparison with existing method(s) The present study, by providing a solution of how to separate common stimulus-elicited brain activities using coupled tensor decomposition, has shed new light on the processing and analysis of ongoing EEG data in multi-subject level. It can also reveal more links between brain responses and the continuous musical stimulus. Conclusions The proposed framework based on coupled tensor decomposition can be successfully applied to group analysis of ongoing EEG data, as it can be reliably inferred that those brain activities we obtained are associated with musical stimulus.

Sustained attention encompasses a cascade of fundamental functions. The human ability to implement a sustained attention task is supported by brain networks that dynamically formed and dissolved through oscillatory synchronization. The decrement of vigilance induced by prolonged task engagement affects sustained attention. However, little is known about which stage or combinations are affected by vigilance decrement. Here, we applied an analysis framework composed of weighted phase lag index (wPLI) and tensor component analysis (TCA) to an EEG dataset collected during 80 min sustained attention task to examine the electrophysiological basis of such effect. We aimed to characterize the phase-coupling networks to untangle different phases involved in sustained attention and study how they are modulated by vigilance decrement. We computed the time–frequency domain wPLI from each block and subject and constructed a fourth-order tensor, containing the time, frequency, functional connectivity (FC), and blocks × subjects. This tensor was subjected to the TCA to identify the interacted and low-dimensional components representing the frequency-specific dynamic FC (fdFC). We extracted four types of neuromakers during a sustained attention task, namely the pre-stimulus alpha right-lateralized parieto-occipital FC, the post-stimulus theta fronto-parieto-occipital FC, delta fronto-parieto-occipital FC, and beta right/left sensorimotor FCs. All these fdFCs were impaired by vigilance decrement. These fdFCs, except for the beta left sensorimotor network, were restored by rewards, although the restoration by reward in the beta right sensorimotor network was transient. These findings provide implications for dissociable effects of vigilance decrement on sustained attention by utilizing the tensor-based framework. Electronic supplementary material The online version of this article (10.1007/s10548-020-00795-0) contains supplementary material, which is available to authorized users.

Deep learning for the automated detection of epileptic seizures has received much attention during recent years. In this work, one dimensional convolutional neural network (1D-CNN) and two dimensional convolutional neural network (2D-CNN) are simultaneously used on electroencephalogram (EEG) data for seizure detection. Firstly, using sliding windows without overlap on raw EEG to obtain the definite one-dimension time EEG segments (1D-T), and continuous wavelet transform (CWT) for 1D-T signals to obtain the two-dimension time-frequency representations (2D-TF). Then, 1D-CNN and 2D-CNN model architectures are used on 1D-T and 2D-TF signals for automatic classification, respectively. Finally, the classification results from 1D-CNN and 2D-CNN are showed. In the two-classification and three-classification problems of seizure detection, the highest accuracy can reach 99.92% and 99.55%, respectively. It shows that the proposed method for a benchmark clinical dataset can achieve good performance in terms of seizure detection. peerReviewed

Sleep scoring is a fundamental but time-consuming process in any sleep laboratory. To speed up the process of sleep scoring without compromising accuracy, this paper develops an automatic sleep scoring toolbox with the capability of multi-signal processing. It allows the user to choose signal types and the number of target classes. Then, an automatic process containing signal pre-processing, feature extraction, classifier training (or prediction) and result correction will be performed. Finally, the application interface displays predicted sleep structure, related sleep parameters and the sleep quality index for reference. To improve the identification accuracy of minority stages, a layer-wise classification strategy is proposed according to the signal characteristics of sleep stages. The context of the current stage is taken into consideration in the correction phase by employing a Hidden Markov Model to study the transition rules of sleep stages in the training dataset. These transition rules will be used for logic classification results. The performance of proposed toolbox has been tested on 100 subjects with an average accuracy of 85.76%. The proposed automatic scoring toolbox would alleviate the burden of the physicians, speed up sleep scoring, and expedite sleep research. peerReviewed

Background: Sleep scoring is an essential but time-consuming process, and therefore automatic sleep scoring is crucial and urgent to help address the growing unmet needs for sleep research. This paper aims to develop a versatile deep-learning architecture to automate sleep scoring using raw polysomnography recordings. Method: The model adopts a linear function to address different numbers of inputs, thereby extending model applications. Two-dimensional convolution neural networks are used to learn features from multi-modality polysomnographic signals, a “squeeze and excitation” block to recalibrate channel-wise features, together with a long short-term memory module to exploit long-range contextual relation. The learnt features are finally fed to the decision layer to generate predictions for sleep stages. Result: Model performance is evaluated on three public datasets. For all tasks with different available channels, our model achieves outstanding performance not only on healthy subjects but even on patients with sleep disorders (SHHS: Acc-0.87, K-0.81; ISRUC: Acc-0.86, K-0.82; Sleep-EDF: Acc-0.86, K-0.81). The highest classification accuracy is achieved by a fusion of multiple polysomnographic signals. Comparison: Compared to state-of-the-art methods that use the same dataset, the proposed model achieves a comparable or better performance, and exhibits low computational cost. Conclusions: The model demonstrates its transferability among different datasets, without changing model architecture or hyper-parameters across tasks. Good model transferability promotes the application of transfer learning on small group studies with mismatched channels. Due to demonstrated availability and versatility, the proposed method can be integrated with diverse polysomnography systems, thereby facilitating sleep monitoring in clinical or routine care. peerReviewed

To examine the electrophysiological underpinnings of the functional networks involved in music listening, previous approaches based on spatial independent component analysis (ICA) have recently been used to ongoing electroencephalography (EEG) and magnetoencephalography (MEG). However, those studies focused on healthy subjects, and failed to examine the group-level comparisons during music listening. Here, we combined group-level spatial Fourier ICA with acoustic feature extraction, to enable group comparisons in frequency-specific brain networks of musical feature processing. It was then applied to healthy subjects and subjects with major depressive disorder (MDD). The music-induced oscillatory brain patterns were determined by permutation correlation analysis between individual time courses of Fourier-ICA components and musical features. We found that (1) three components, including a beta sensorimotor network, a beta auditory network and an alpha medial visual network, were involved in music processing among most healthy subjects; and that (2) one alpha lateral component located in the left angular gyrus was engaged in music perception in most individuals with MDD. The proposed method allowed the statistical group comparison, and we found that: (1) the alpha lateral component was activated more strongly in healthy subjects than in the MDD individuals, and that (2) the derived frequency-dependent networks of musical feature processing seemed to be altered in MDD participants compared to healthy subjects. The proposed pipeline appears to be valuable for studying disrupted brain oscillations in psychiatric disorders during naturalistic paradigms.

Background Component splitting at higher model orders is a widely accepted finding for independent component analysis (ICA) of functional magnetic resonance imaging (fMRI) data. However, our recent study found that intact components occurred with subcomponents at higher model orders. New method This study investigated model order effects on ICA of resting-state complex-valued fMRI data from 82 subjects, which included 40 healthy controls (HCs) and 42 schizophrenia patients. In addition, we explored underlying causes for distinct component splitting between complex-valued data and magnitude-only data by examining model order effects on ICA of phase fMRI data. A best run selection method was proposed to combine subject averaging and a one-sample t-test. We selected the default mode network (DMN)-, visual-, and sensorimotor-related components from the best run of ICA at varying model orders from 10 to 140. Results Results show that component integration occurred in complex-valued and phase analyses, whereas component splitting emerged in magnitude-only analysis with increasing model order. Incorporation of phase data appears to play a complementary role in preserving integrity of brain networks. Comparison with existing method(s) When compared with magnitude-only analysis, the intact DMN component obtained in complex-valued analysis at higher model orders exhibited highly significant subject-level differences between HCs and patients with schizophrenia. We detected significantly higher activity and variation in anterior areas for HCs and in posterior areas for patients with schizophrenia. Conclusions These results demonstrate the potential of complex-valued fMRI data to contribute generally and specifically to brain network analysis in identification of schizophrenia-related changes.