Cardiac sympathetic-vagal activity initiates a functional brain-body response to emotional arousal.

A century-long debate on bodily states and emotions persists. While the involvement of bodily activity in emotion physiology is widely recognized, the specificity and causal role of such activity related to brain dynamics has not yet been demonstrated. We hypothesize that the peripheral neural control on cardiovascular activity prompts and sustains brain dynamics during an emotional experience, so these afferent inputs are processed by the brain by triggering a concurrent efferent information transfer to the body. To this end, we investigated the functional brain­heart interplay under emotion elicitation in publicly available data from 62 healthy subjects using a computational model based on synthetic data generation of electroencephalography and electrocardiography signals. Our findings show that sympathovagal activity plays a leading and causal role in initiating the emotional response, in which ascending modulations from vagal activity precede neural dynamics and correlate to the reported level of arousal. The subsequent dynamic interplay observed between the central and autonomic nervous systems sustains the processing of emotional arousal. These findings should be particularly revealing for the psychophysiology and neuroscience of emotions.

Stochastic brain dynamics exhibits differential regional distribution and maturation-related changes.

Functional magnetic resonance imaging (fMRI) is a powerful non-invasive method for studying brain function by analyzing blood oxygenation level-dependent (BOLD) signals. These signals arise from intricate interplays of deterministic and stochastic biological elements. Quantifying the stochastic part is challenging due to its reliance on assumptions about the deterministic segment. We present a methodological framework to estimate intrinsic stochastic brain dynamics in fMRI data without assuming deterministic dynamics. Our approach utilizes Approximate Entropy and its behavior in noisy series to identify and characterize dynamical noise in unobservable fMRI dynamics. Applied to extensive fMRI datasets (645 Cam-CAN, 1086 Human Connectome Project subjects), we explore lifelong maturation of intrinsic brain noise. Findings indicate 10% to 60% of fMRI signal power is due to intrinsic stochastic brain elements, varying by age. These components demonstrate a physiological role of neural noise which shows a distinct distributions across brain regions and increase linearly during maturation.

Intracortical brain-heart interplay: An EEG model source study of sympathovagal changes.

The interplay between cerebral and cardiovascular activity, known as the functional brain-heart interplay (BHI), and its temporal dynamics, have been linked to a plethora of physiological and pathological processes. Various computational models of the brain-heart axis have been proposed to estimate BHI non-invasively by taking advantage of the time resolution offered by electroencephalograph (EEG) signals. However, investigations into the specific intracortical sources responsible for this interplay have been limited, which significantly hampers existing BHI studies. This study proposes an analytical modeling framework for estimating the BHI at the source-brain level. This analysis relies on the low-resolution electromagnetic tomography sources localization from scalp electrophysiological recordings. BHI is then quantified as the functional correlation between the intracortical sources and cardiovascular dynamics. Using this approach, we aimed to evaluate the reliability of BHI estimates derived from source-localized EEG signals as compared with prior findings from neuroimaging methods. The proposed approach is validated using an experimental dataset gathered from 32 healthy individuals who underwent standard sympathovagal elicitation using a cold pressor test. Additional resting state data from 34 healthy individuals has been analysed to assess robustness and reproducibility of the methodology. Experimental results not only confirmed previous findings on activation of brain structures affecting cardiac dynamics (e.g., insula, amygdala, hippocampus, and anterior and mid-cingulate cortices) but also provided insights into the anatomical bases of brain-heart axis. In particular, we show that the bidirectional activity of electrophysiological pathways of functional brain-heart communication increases during cold pressure with respect to resting state, mainly targeting neural oscillations in the δ $$ \delta $$ , ß $$ \beta $$ , and γ $$ \gamma $$ bands. The proposed approach offers new perspectives for the investigation of functional BHI that could also shed light on various pathophysiological conditions.

Microstates of the cortical brain-heart axis.

Electroencephalographic (EEG) microstates are brain states with quasi-stable scalp topography. Whether such states extend to the body level, that is, the peripheral autonomic nerves, remains unknown. We hypothesized that microstates extend at the brain-heart axis level as a functional state of the central autonomic network. Thus, we combined the EEG and heartbeat dynamics series to estimate the directional information transfer originating in the cortex targeting the sympathovagal and parasympathetic activity oscillations and vice versa for the afferent functional direction. Data were from two groups of participants: 36 healthy volunteers who were subjected to cognitive workload induced by mental arithmetic, and 26 participants who underwent physical stress induced by a cold pressure test. All participants were healthy at the time of the study. Based on statistical testing and goodness-of-fit evaluations, we demonstrated the existence of microstates of the functional brain-heart axis, with emphasis on the cerebral cortex, since the microstates are derived from EEG. Such nervous-system microstates are spatio-temporal quasi-stable states that exclusively refer to the efferent brain-to-heart direction. We demonstrated brain-heart microstates that could be associated with specific experimental conditions as well as brain-heart microstates that are non-specific to tasks.

Dynamic fluctuations in ascending heart-to-brain communication under mental stress.

Dynamical information exchange between central and autonomic nervous systems, as referred to functional brain-heart interplay, occurs during emotional and physical arousal. It is well documented that physical and mental stress lead to sympathetic activation. Nevertheless, the role of autonomic inputs in nervous system-wise communication under mental stress is yet unknown. In this study, we estimated the causal and bidirectional neural modulations between electroencephalogram (EEG) oscillations and peripheral sympathetic and parasympathetic activities using a recently proposed computational framework for a functional brain-heart interplay assessment, namely the sympathovagal synthetic data generation model. Mental stress was elicited in 37 healthy volunteers by increasing their cognitive demands throughout three tasks associated with increased stress levels. Stress elicitation induced an increased variability in sympathovagal markers, as well as increased variability in the directional brain-heart interplay. The observed heart-to-brain interplay was primarily from sympathetic activity targeting a wide range of EEG oscillations, whereas variability in the efferent direction seemed mainly related to EEG oscillations in the γ band. These findings extend current knowledge on stress physiology, which mainly referred to top-down neural dynamics. Our results suggest that mental stress may not cause an increase in sympathetic activity exclusively as it initiates a dynamic fluctuation within brain-body networks including bidirectional interactions at a brain-heart level. We conclude that directional brain-heart interplay measurements may provide suitable biomarkers for a quantitative stress assessment and bodily feedback may modulate the perceived stress caused by increased cognitive demand.

Functional assessment of bidirectional cortical and peripheral neural control on heartbeat dynamics: A brain-heart study on thermal stress.

The study of functional Brain-Heart Interplay (BHI) from non-invasive recordings has gained much interest in recent years. Previous endeavors aimed at understanding how the two dynamical systems exchange information, providing novel holistic biomarkers and important insights on essential cognitive aspects and neural system functioning. However, the interplay between cardiac sympathovagal and cortical oscillations still has much room for further investigation. In this study, we introduce a new computational framework for a functional BHI assessment, namely the Sympatho-Vagal Synthetic Data Generation Model, combining cortical (electroencephalography, EEG) and peripheral (cardiac sympathovagal) neural dynamics. The causal, bidirectional neural control on heartbeat dynamics was quantified on data gathered from 26 human volunteers undergoing a cold-pressor test. Results show that thermal stress induces heart-to-brain functional interplay sustained by EEG oscillations in the delta and gamma bands, primarily originating from sympathetic activity, whereas brain-to-heart interplay originates over central brain regions through sympathovagal control. The proposed methodology provides a viable computational tool for the functional assessment of the causal interplay between cortical and cardiac neural control.

Activation of brain-heart axis during REM sleep: a trigger for dreaming.

Dreams may be recalled after awakening from sleep following a defined electroencephalographic pattern that involves local decreases in low-frequency activity in the posterior cortical regions. Although a dreaming experience implies bodily changes at many organ, system, and timescale levels, the entity and causal role of such peripheral changes in a conscious dream experience are unknown. We performed a comprehensive, causal, multivariate analysis of physiological signals acquired during rapid eye movement (REM) sleep at night, including high-density electroencephalography and peripheral dynamics including electrocardiography and blood pressure. In this preliminary study, we investigated multiple recalls and nonrecalls of dream experiences using data from nine healthy volunteers. The aim was not only to investigate the changes in central and autonomic dynamics associated with dream recalls and nonrecalls, but also to characterize the central-peripheral dynamical and (causal) directional interactions, and the temporal relations of the related arousals upon awakening. We uncovered a brain-body network that drives a conscious dreaming experience that acts with specific interaction and time delays. Such a network is sustained by the blood pressure dynamics and the increasing functional information transfer from the neural heartbeat regulation to the brain. We conclude that bodily changes play a crucial and causative role in a conscious dream experience during REM sleep.

Heart rate variability and the dawn of complex physiological signal analysis: methodological and clinical perspectives.

Spontaneous beat-to-beat variations of heart rate (HR) have intrigued scientists and casual observers for centuries; however, it was not until the 1970s that investigators began to apply engineering tools to the analysis of these variations, fostering the field we now know as heart rate variability or HRV. Since then, the field has exploded to not only include a wide variety of traditional linear time and frequency domain applications for the HR signal, but also more complex linear models that include additional physiological parameters such as respiration, arterial blood pressure, central venous pressure and autonomic nerve signals. Most recently, the field has branched out to address the nonlinear components of many physiological processes, the complexity of the systems being studied and the important issue of specificity for when these tools are applied to individuals. When the impact of all these developments are combined, it seems likely that the field of HRV will soon begin to realize its potential as an important component of the toolbox used for diagnosis and therapy of patients in the clinic. This article is part of the theme issue 'Advanced computation in cardiovascular physiology: new challenges and opportunities'.

Advanced computation in cardiovascular physiology: new challenges and opportunities.

Recent developments in computational physiology have successfully exploited advanced signal processing and artificial intelligence tools for predicting or uncovering characteristic features of physiological and pathological states in humans. While these advanced tools have demonstrated excellent diagnostic capabilities, the high complexity of these computational 'black boxes' may severely limit scientific inference, especially in terms of biological insight about both physiology and pathological aberrations. This theme issue highlights current challenges and opportunities of advanced computational tools for processing dynamical data reflecting autonomic nervous system dynamics, with a specific focus on cardiovascular control physiology and pathology. This includes the development and adaptation of complex signal processing methods, multivariate cardiovascular models, multiscale and nonlinear models for central-peripheral dynamics, as well as deep and transfer learning algorithms applied to large datasets. The width of this perspective highlights the issues of specificity in heartbeat-related features and supports the need for an imminent transition from the black-box paradigm to explainable and personalized clinical models in cardiovascular research. This article is part of the theme issue 'Advanced computation in cardiovascular physiology: new challenges and opportunities'.

Functional brain-heart interplay extends to the multifractal domain.

The study of functional brain-heart interplay has provided meaningful insights in cardiology and neuroscience. Regarding biosignal processing, this interplay involves predominantly neural and heartbeat linear dynamics expressed via time and frequency domain-related features. However, the dynamics of central and autonomous nervous systems show nonlinear and multifractal behaviours, and the extent to which this behaviour influences brain-heart interactions is currently unknown. Here, we report a novel signal processing framework aimed at quantifying nonlinear functional brain-heart interplay in the non-Gaussian and multifractal domains that combines electroencephalography (EEG) and heart rate variability series. This framework relies on a maximal information coefficient analysis between nonlinear multiscale features derived from EEG spectra and from an inhomogeneous point-process model for heartbeat dynamics. Experimental results were gathered from 24 healthy volunteers during a resting state and a cold pressor test, revealing that synchronous changes between brain and heartbeat multifractal spectra occur at higher EEG frequency bands and through nonlinear/complex cardiovascular control. We conclude that significant bodily, sympathovagal changes such as those elicited by cold-pressure stimuli affect the functional brain-heart interplay beyond second-order statistics, thus extending it to multifractal dynamics. These results provide a platform to define novel nervous-system-targeted biomarkers. This article is part of the theme issue 'Advanced computation in cardiovascular physiology: new challenges and opportunities'.

A protocol for a multicentre, parallel-group, pragmatic randomised controlled trial to evaluate the NEVERMIND system in preventing and treating depression in patients with severe somatic conditions.

BACKGROUND: Depressive symptoms are common in individuals suffering from severe somatic conditions. There is a lack of interventions and evidence-based interventions aiming to reduce depressive symptoms in patients with severe somatic conditions. The aim of the NEVERMIND project is to address these issues and provide evidence by testing the NEVERMIND system, designed to reduce and prevent depressive symptoms in comparison to treatment as usual. METHODS: The NEVERMIND study is a parallel-groups, pragmatic randomised controlled trial to assess the effectiveness of the NEVERMIND system in reducing depressive symptoms among individuals with severe somatic conditions. The NEVERMIND system comprises a smart shirt and a user interface, in the form of a mobile application. The system is a real-time decision support system, aiming to predict the severity and onset of depressive symptoms by modelling the well-being condition of patients based on physiological data, body movement, and the recurrence of social interactions. The study includes 330 patients who have a diagnosis of myocardial infarction, breast cancer, prostate cancer, kidney failure, or lower limb amputation. Participants are randomised in blocks of ten to either the NEVERMIND intervention or treatment as usual as the control group. Clinical interviews and structured questionnaires are administered at baseline, at 12 weeks, and 24 weeks to assess whether the NEVERMIND system is superior to treatment as usual. The endpoint of primary interest is Beck Depression Inventory II (BDI-II) at 12 weeks defined as (i) the severity of depressive symptoms as measured by the BDI-II. Secondary outcomes include prevention of the onset of depressive symptoms, changes in quality of life, perceived stigma, and self-efficacy. DISCUSSION: There is a lack of evidence-based interventions aiming to reduce and prevent depressive symptoms in patients with severe somatic conditions. If the NEVERMIND system is effective, it will provide healthcare systems with a novel and innovative method to attend to depressive symptoms in patients with severe somatic conditions. TRIAL REGISTRATION: DRKS00013391. Registered 23 November 2017.

Exploiting upper-limb functional principal components for human-like motion generation of anthropomorphic robots.

BACKGROUND: Human-likeliness of robot movements is a key component to enable a safe and effective human-robot interaction, since it contributes to increase acceptance and motion predictability of robots that have to closely interact with people, e.g. for assistance and rehabilitation purposes. Several parameters have been used to quantify how much a robot behaves like a human, which encompass aspects related to both the robot appearance and motion. The latter point is fundamental to allow the operator to interpret robotic actions, and plan a meaningful reactions. While different approaches have been presented in literature, which aim at devising bio-aware control guidelines, a direct implementation of human actions for robot planning is not straightforward, still representing an open issue in robotics. METHODS: We propose to embed a synergistic representation of human movements for robot motion generation. To do this, we recorded human upper-limb motions during daily living activities. We used functional Principal Component Analysis (fPCA) to extract principal motion patterns. We then formulated the planning problem by optimizing the weights of a reduced set of these components. For free-motions, our planning method results into a closed form solution which uses only one principal component. In case of obstacles, a numerical routine is proposed, incrementally enrolling principal components until the problem is solved with a suitable precision. RESULTS: Results of fPCA show that more than 80% of the observed variance can be explained by only three functional components. The application of our method to different meaningful movements, with and without obstacles, show that our approach is able to generate complex motions with a very reduced number of functional components. We show that the first synergy alone accounts for the 96% of cost reduction and that three components are able to achieve a satisfactory motion reconstruction in all the considered cases. CONCLUSIONS: In this work we moved from the analysis of human movements via fPCA characterization to the design of a novel human-like motion generation algorithm able to generate, efficiently and with a reduced set of basis elements, several complex movements in free space, both in free motion and in case of obstacle avoidance tasks.

fNIRS Complexity Analysis for the Assessment of Motor Imagery and Mental Arithmetic Tasks.

Conventional methods for analyzing functional near-infrared spectroscopy (fNIRS) signals primarily focus on characterizing linear dynamics of the underlying metabolic processes. Nevertheless, linear analysis may underrepresent the true physiological processes that fully characterizes the complex and nonlinear metabolic activity sustaining brain function. Although there have been recent attempts to characterize nonlinearities in fNIRS signals in various experimental protocols, to our knowledge there has yet to be a study that evaluates the utility of complex characterizations of fNIRS in comparison to standard methods, such as the mean value of hemoglobin. Thus, the aim of this study was to investigate the entropy of hemoglobin concentration time series obtained from fNIRS signals and perform a comparitive analysis with standard mean hemoglobin analysis of functional activation. Publicly available data from 29 subjects performing motor imagery and mental arithmetics tasks were exploited for the purpose of this study. The experimental results show that entropy analysis on fNIRS signals may potentially uncover meaningful activation areas that enrich and complement the set identified through a traditional linear analysis.

Assessing Complexity in Physiological Systems through Biomedical Signals Analysis.

The idea that most physiological systems are complex has become increasingly popular in recent decades [...].

Lateralization of directional brain-heart information transfer during visual emotional elicitation.

Previous studies have characterized the physiological interactions between central nervous system (brain) and peripheral cardiovascular system (heart) during affective elicitation in healthy subjects; however, questions related to the directionality of this functional interplay have been gaining less attention from the scientific community. Here, we explore brain-heart interactions during visual emotional elicitation in healthy subjects using measures of Granger causality (GC), a widely used descriptor of causal influences between two dynamical systems. The proposed approach inferences causality between instantaneous cardiovagal dynamics estimated from inhomogeneous point-process models of the heartbeat and high-density electroencephalogram (EEG) dynamics in 22 healthy subjects who underwent pleasant/unpleasant affective elicitation by watching pictures from the International Affective Picture System database. Particularly, we calculated the GC indexes between the EEG spectrogram in the canonical θ-, α-, ß-, and γ-bands and both the instantaneous mean heart rate and its continuous parasympathetic modulations (i.e., the instantaneous HF power). Thus we looked for significant statistical differences among GC values estimated during the resting state, neutral elicitation, and pleasant/unpleasant arousing elicitation. As compared with resting state, coupling strength increases significantly in the left hemisphere during positive stimuli and in the right hemisphere during negative stimuli. Our results further reveal a correlation between emotional valence and lateralization of the dynamical information transfer going from brain-to-heart, mainly localized in the prefrontal, somatosensory, and posterior cortexes, and of the information transfer from heart-to-brain, mainly reflected into the fronto-parietal cortex oscillations in the γ-band (30-45 Hz).

A Parsimonious Granger Causality Formulation for Capturing Arbitrarily Long Multivariate Associations.

High-frequency neuroelectric signals like electroencephalography (EEG) or magnetoencephalography (MEG) provide a unique opportunity to infer causal relationships between local activity of brain areas. While causal inference is commonly performed through classical Granger causality (GC) based on multivariate autoregressive models, this method may encounter important limitations (e.g., data paucity) in the case of high dimensional data from densely connected systems like the brain. Additionally, physiological signals often present long-range dependencies which commonly require high autoregressive model orders/number of parameters. We present a generalization of autoregressive models for GC estimation based on Wiener-Volterra decompositions with Laguerre polynomials as basis functions. In this basis, the introduction of only one additional global parameter allows to capture arbitrary long dependencies without increasing model order, hence retaining model simplicity, linearity and ease of parameters estimation. We validate our method in synthetic data generated from families of complex, densely connected networks and demonstrate superior performance as compared to classical GC. Additionally, we apply our framework to studying the directed human brain connectome through MEG data from 89 subjects drawn from the Human Connectome Project (HCP) database, showing that it is able to reproduce current knowledge as well as to uncover previously unknown directed influences between cortical and limbic brain regions.

Neuron-glia cross talk revealed in reverberating networks by simultaneous extracellular recording of spikes and astrocytes' glutamate transporter and K+ currents.

Astrocytes uptake synaptically released glutamate with electrogenic transporters (GluT) and buffer the spike-dependent extracellular K+ excess with background K+ channels. We studied neuronal spikes and the slower astrocytic signals on reverberating neocortical cultures and organotypic slices from mouse brains. Spike trains and glial responses were simultaneously captured from individual sites of multielectrode arrays (MEA) by splitting the recorded traces into appropriate filters and reconstructing the original signal by deconvolution. GluT currents were identified by using dl-threo-ß-benzyloxyaspartate (TBOA). K+ currents were blocked by 30 µM Ba2+, suggesting a major contribution of inwardly rectifying K+ currents. Both types of current were tightly correlated with the spike rate, and their astrocytic origin was tested in primary cultures by blocking glial proliferation with cytosine ß-d-arabinofuranoside (AraC). The spike-related, time-locked inward and outward K+ currents in different regions of the astrocyte syncytium were consistent with the assumptions of the spatial K+ buffering model. In organotypic slices from ventral tegmental area and prefrontal cortex, the GluT current amplitudes exceeded those observed in primary cultures by several orders of magnitude, which allowed to directly measure transporter currents with a single electrode. Simultaneously measuring cell signals displaying widely different amplitudes and kinetics will help clarify the neuron-glia interplay and make it possible to follow the cross talk between different cell types in excitable as well as nonexcitable tissue.

Physiological Noise: Definition, Estimation, and Characterization in Complex Biomedical Signals.

BACKGROUND: Nonlinear physiological systems exhibit complex dynamics driven by intrinsic dynamical noise. In cases where there is no specific knowledge or assumption about system dynamics, such as in physiological systems, it is not possible to formally estimate noise. AIM: We introduce a formal method to estimate the power of dynamical noise, referred to as physiological noise, in a closed form, without specific knowledge of the system dynamics. METHODOLOGY: Assuming that noise can be modeled as a sequence of independent, identically distributed (IID) random variables on a probability space, we demonstrate that physiological noise can be estimated through a nonlinear entropy profile. We estimated noise from synthetic maps that included autoregressive, logistic, and Pomeau-Manneville systems under various conditions. Noise estimation is performed on 70 heart rate variability series from healthy and pathological subjects, and 32 electroencephalographic (EEG) healthy series. RESULTS: Our results showed that the proposed model-free method can discern different noise levels without any prior knowledge of the system dynamics. Physiological noise accounts for around 11% of the overall power observed in EEG signals and approximately 32% to 65% of the power related to heartbeat dynamics. Cardiovascular noise increases in pathological conditions compared to healthy dynamics, and cortical brain noise increases during mental arithmetic computations over the prefrontal and occipital regions. Brain noise is differently distributed across cortical regions. CONCLUSION: Physiological noise is very part neurobiological dynamics and can be measured using the proposed framework in any biomedical series.

Uncovering Hemispheric Asymmetry and Directed Oscillatory Brain-Heart Interplay in Anxiety Processing: An fMRI Study.

Brain-heart interactions (BHI) are critical for generating and processing emotions, including anxiety. Understanding specific neural correlates would be instrumental for greater comprehension and potential therapeutic interventions of anxiety disorders. While prior work has implicated the pontine structure as a central processor in cardiac regulation in anxiety, the distributed nature of anxiety processing across the cortex remains elusive. To address this, we performed a whole-brain-heart analysis using the full frequency directed transfer function to study resting-state spectral differences in BHI between high and low anxiety groups undergoing fMRI scans. Our findings revealed a hemispheric asymmetry in low-frequency interplay (0.05 Hz - 0.15 Hz) characterized by ascending BHI to the left insula and descending BHI from the right insula. Furthermore, we provide evidence supporting the "pacemaker hypothesis", highlighting the pons' function in regulating cardiac activity. Higher frequency interplay (0.2 Hz - 0.4Hz) demonstrate a preference for ascending interactions, particularly towards ventral prefrontal cortical activity in high anxiety groups, suggesting the heart's role in triggering a cognitive response to regulate anxiety. These findings highlight the impact of anxiety on BHI, contributing to a better understanding of its effect on the resting-state fMRI signal, with further implications for potential therapeutic interventions in treating anxiety disorders.

Impaired brain-heart axis in focal epilepsy: Alterations in information flow and implications for seizure dynamics.

This study delves into functional brain-heart interplay (BHI) dynamics during interictal periods before and after seizure events in focal epilepsy. Our analysis focuses on elucidating the causal interaction between cortical and autonomic nervous system (ANS) oscillations, employing electroencephalography and heart rate variability series. The dataset for this investigation comprises 47 seizure events from 14 independent subjects, obtained from the publicly available Siena Dataset. Our findings reveal an impaired brain-heart axis especially in the heart-to-brain functional direction. This is particularly evident in bottom-up oscillations originating from sympathovagal activity during the transition between preictal and postictal periods. These results indicate a pivotal role of the ANS in epilepsy dynamics. Notably, the brain-to-heart information flow targeting cardiac oscillations in the low-frequency band does not display significant changes. However, there are noteworthy changes in cortical oscillations, primarily originating in central regions, influencing heartbeat oscillations in the high-frequency band. Our study conceptualizes seizures as a state of hyperexcitability and a network disease affecting both cortical and peripheral neural dynamics. Our results pave the way for a deeper understanding of BHI in epilepsy, which holds promise for the development of advanced diagnostic and therapeutic approaches also based on bodily neural activity for individuals living with epilepsy.

This study focuses on brain-heart interplay (BHI) during pre- and postictal periods surrounding seizures. Employing multichannel EEG and heart rate variability data from subjects with focal epilepsy, our analysis reveals a disrupted brain-heart axis dynamic, particularly in the heart-to-brain direction. Notably, sympathovagal activity alterations during preictal to postictal transitions underscore the autonomic nervous system's pivotal role in epilepsy dynamics. While brain-to-heart information flow targeting low-frequency band cardiac oscillations remains stable, significant changes occur in cortical oscillations, predominantly in central regions, influencing high-frequeny-band heartbeat oscillations, that is, vagal activity. Viewing seizures as states of hyperexcitability and confirming focal epilepsy as a network disease affecting both central and peripheral neural dynamics, our study enhances understanding of BHI in epilepsy. These findings offer potential for advanced diagnostic and therapeutic approaches grounded in bodily neural activity for individuals with epilepsy.