Previous studies have shown that emotions elicit responses in various physiological systems that underlie emotional experience. However, these studies have focused almost solely on how emotions alter the activity of individual physiological systems. Moreover, recent research has highlighted that the human organism is an integrated entity, wherein physiological systems constantly interact to form a dynamic network. Therefore, in this study, we attempted to evaluate whether the physiological network changes with different emotions. We employed multivariate transfer entropy to assess the interaction between six types of emotion-related physiological signals and subsequently constructed physiological networks across forty emotional trials. Additionally, topological parameters were applied to each physiological network to quantify its overall structure. The results revealed that emotion has a significant impact on the topology of physiological networks, indicating the integrated behavior of physiological systems in response to emotion at the macro-scale level. Furthermore, different physiological systems may play diverse roles within the physiological network. In summary, we propose a system-level integrative approach, referred to as the physiological network, to study the physiological responses to emotion. This not only contributes to a deeper understanding of emotional mechanisms but also demonstrates potential for disease research and engineering applications.

The power spectra estimated from the brain recordings are the mixed representation of aperiodic transient activity and periodic oscillations, i.e., aperiodic component (AC) and periodic component (PC). Quantitative neurophysiology requires precise decomposition preceding parameterizing each component. However, the shape, statistical distribution, scale, and mixing mechanism of AC and PCs are unclear, challenging the effectiveness of current popular parametric models such as FOOOF, IRASA, BOSC, etc. Here, ξ- π was proposed to decompose the neural spectra by embedding the nonparametric spectra estimation with penalized Whittle likelihood and the shape language modeling into the expectation maximization framework. ξ- π was validated on the synthesized spectra with loss statistics and on the sleep EEG and the large sample iEEG with evaluation metrics and neurophysiological evidence. Compared to FOOOF, both the simulation presenting shape irregularities and the batch simulation with multiple isolated peaks indicated that ξ- π improved the fit of AC and PCs with less loss and higher F1-score in recognizing the centering frequencies and the number of peaks; the sleep EEG revealed that ξ- π produced more distinguishable AC exponents and improved the sleep state classification accuracy; the iEEG showed that ξ- π approached the clinical findings in peak discovery. Overall, ξ- π offered good performance in the spectra decomposition, which allows flexible parameterization using descriptive statistics or kernel functions. ξ- π is a seminal tool for brain signal decoding in fields such as cognitive neuroscience, brain-computer interface, neurofeedback, and brain diseases.

Current high-throughput data acquisition technologies probe dynamical systems with different imaging modalities, generating massive data sets at different spatial and temporal resolutions-posing challenging problems in multimodal data fusion. A case in point is the attempt to parse out the brain structures and networks that underpin human cognitive processes by analysis of different neuroimaging modalities (functional MRI, EEG, NIRS, etc.). We emphasize that the multimodal, multiscale nature of neuroimaging data is well reflected by a multiway (tensor) structure where the underlying processes can be summarized by a relatively small number of components or “atoms.” We introduce Markov-Penrose diagrams-an integration of Bayesian DAG and tensor network notation in order to analyze these models. These diagrams not only clarify matrix and tensor EEG and fMRI time/frequency analysis and inverse problems, but also help understand multimodal fusion via multiway partial least squares and coupled matrix-tensor factorization. We show here, for the first time, that Granger causal analysis of brain networks is a tensor regression problem, thus allowing the atomic decomposition of brain networks. Analysis of EEG and fMRI recordings shows the potential of the methods and suggests their use in other scientific domains.

Absence epilepsy is believed to be associated with the abnormal interactions between the cerebral cortex and thalamus. Besides the direct coupling, anatomical evidence indicates that the cerebral cortex and thalamus also communicate indirectly through an important intermediate bridge-basal ganglia. It has been thus postulated that the basal ganglia might play key roles in the modulation of absence seizures, but the relevant biophysical mechanisms are still not completely established. Using a biophysically based model, we demonstrate here that the typical absence seizure activities can be controlled and modulated by the direct GABAergic projections from the substantia nigra pars reticulata (SNr) to either the thalamic reticular nucleus (TRN) or the specific relay nuclei (SRN) of thalamus, through different biophysical mechanisms. Under certain conditions, these two types of seizure control are observed to coexist in the same network. More importantly, due to the competition between the inhibitory SNr-TRN and SNr-SRN pathways, we find that both decreasing and increasing the activation of SNr neurons from the normal level may considerably suppress the generation of spike-and-slow wave discharges in the coexistence region. Overall, these results highlight the bidirectional functional roles of basal ganglia in controlling and modulating absence seizures, and might provide novel insights into the therapeutic treatments of this brain disorder.