Globally conditioned Granger causality in brain-brain and brain-heart interactions: a combined heart rate variability/ultra-high-field (7 T) functional magnetic resonance imaging study

Facial emotion processing in patients with borderline personality disorder as compared with healthy controls: an fMRI and ECG study

BACKGROUND

Maladaptive behaviors and interpersonal difficulties in patients with borderline personality disorder (BPD) seem connected to biased facial emotion processing. This bias is often accompanied by heightened amygdala activity in patients with BPD as compared to healthy controls. However, functional magnetic resonance imaging (fMRI) studies exploring differences between patients and healthy controls in facial emotion processing have produced divergent results. The current study explored fMRI and heart rate variability (HRV) correlates of negative facial emotion processing in patients with BPD and healthy controls.

METHODS

The study included 30 patients with BPD (29 females; age: M = 24.22, SD = 5.22) and 30 healthy controls (29 females; M = 24.66, SD = 5.28). All participants underwent the "faces" task, an emotional face perception task, in an fMRI session simultaneously with ECG. In this task, participants are presented with emotional expressions of disgust, sadness, and fear (as a negative condition) and with the same pictures in a scrambled version (as a neutral condition).

RESULTS

We found no differences in brain activity between patients with BPD and healthy controls when processing negative facial expressions as compared to neutral condition. We observed activation in large-scale brain areas in both groups when presented with negative facial expressions as compared to neutral condition. Patients with BPD displayed lower HRV than healthy controls in both conditions. However, there were no significant associations between HRV and amygdala activity and BPD symptoms.

CONCLUSION

The results of this study indicate no abnormal brain activity during emotional facial processing in patients with BPD. This result contrasts with previous studies and more studies are needed to clarify the relationship between facial emotion processing and brain activity in patients with BPD. Possible reasons for the absence of brain activity differences are discussed in the study. Consistent with previous findings, patients showed lower HRV than healthy controls. However, HRV was not associated with amygdala activity and BPD symptoms.

Usefulness and limitations of convergent cross sorting and continuity scaling methods for their application in simulated and real-world time series

Causality detection methods are valuable tools for detecting causal links in complex systems. The efficiency of continuity scaling (CS) and the convergent cross sorting (CSS) methods to detect causality was analysed. Usefulness and limitations of both methods in their application to simulated and real-world time series was explored under different scenarios. We find that CS is more robust and efficient than the CSS method for all simulated systems, even when increasing noise levels were considered. Both methods were not able to infer causality when time series with a marked difference in their main frequencies were analysed. Minimum time-series length required for the detection of a causal link depends on intrinsic system dynamics and on the method selected to detect it. Using simulated time series, only the CS method was capable to detect bidirectional causality. Causality detection, using the CS method, should at least include: (i) causality strength convergence analysis, (ii) statistical tests of significance, (iii) time-series standardization, and (iv) causality strength ratios as a strength indicator of relative causality between systems. Causality cannot be detected by either method in simulated time series that exhibit generalized synchronization.

Neuroimaging Studies of the Neural Correlates of Heart Rate Variability: A Systematic Review

Direct and indirect links between brain regions and cardiac function have been reported. We performed a systematic literature review to summarize current knowledge regarding the associations of heart rate variability (HRV) and brain region morphology, activity and connectivity involved in autonomic control at rest in healthy subjects. Both positive and negative correlations of cortical thickness and gray matter volumes of brain structures with HRV were observed. The strongest were found for a cluster located within the cingulate cortex. A decline in HRV, as well as cortical thickness with increasing age, especially in the orbitofrontal cortex were noted. When associations of region-specific brain activity with HRV were examined, HRV correlated most strongly with activity in the insula, cingulate cortex, frontal and prefrontal cortices, hippocampus, thalamus, striatum and amygdala. Furthermore, significant correlations, largely positive, between HRV and brain region connectivity (in the amygdala, cingulate cortex and prefrontal cortex) were observed. Notably, right-sided neural structures may be preferentially involved in heart rate and HRV control. However, the evidence for left hemispheric control of cardiac vagal function has also been reported. Our findings provide support for the premise that the brain and the heart are interconnected by both structural and functional networks and indicate complex multi-level interactions. Further studies of brain-heart associations promise to yield insights into their relationship to health and disease.

Pairwise and higher-order measures of brain-heart interactions in children with temporal lobe epilepsy

OBJECTIVE

While it is well-known that epilepsy has a clear impact on the activity of both the central nervous system (CNS) and the autonomic nervous system (ANS), its role on the complex interplay between CNS and ANS has not been fully elucidated yet. In this work, pairwise and higher-order predictability measures based on the concepts of Granger causality (GC) and Partial Information Decomposition (PID) were applied on time series of electroencephalographic (EEG) brain wave amplitude and heart rate variability (HRV) in order to investigate directed brain-heart interactions associated with the occurrence of focal epilepsy.

APPROACH

HRV and the envelopes of δ and α EEG activity recorded from ipsilateral (ipsi-EEG) and contralateral (contra-EEG) scalp regions were analyzed in 18 children suffering from temporal lobe epilepsy monitored during pre-ictal, ictal and post-ictal periods. After linear parametric model identification, we compared pairwise GC measures computed between HRV and a single EEG component with PID measures quantifying the unique, redundant and synergistic information transferred from ipsi-EEG and contra-EEG to HRV.

MAIN RESULTS

The analysis of GC revealed a dominance of the information transfer from EEG to HRV and negligible transfer from HRV to EEG, suggesting that CNS activities drive the ANS modulation of the heart rhythm, but did not evidence clear differences between δ and α rhythms, ipsi-EEG and contra-EEG, or pre- and post-ictal periods. On the contrary, PID revealed that epileptic seizures induce a reorganization of the interactions from brain to heart, as the unique predictability of HRV originated from the ipsi-EEG for the δ waves and from the contra-EEG for the α waves in the pre-ictal phase, while these patterns were reversed after the seizure.

SIGNIFICANCE

These results highlight the importance of considering higher-order interactions elicited by PID for the study of the neuro-autonomic effects of focal epilepsy, and may have neurophysiological and clinical implications.

Assessing rheoencephalography dynamics through analysis of the interactions among brain and cardiac networks during general anesthesia

Cerebral blood flow (CBF) reflects the rate of delivery of arterial blood to the brain. Since no nutrients, oxygen or water can be stored in the cranial cavity due to space and pressure restrictions, a continuous perfusion of the brain is critical for survival. Anesthetic procedures are known to affect cerebral hemodynamics, but CBF is only monitored in critical patients due, among others, to the lack of a continuous and affordable bedside monitor for this purpose. A potential solution through bioelectrical impedance technology, also known as rheoencephalography (REG), is proposed, that could fill the existing gap for a low-cost and effective CBF monitoring tool. The underlying hypothesis is that REG signals carry information on CBF that might be recovered by means of the application of advanced signal processing techniques, allowing to track CBF alterations during anesthetic procedures. The analysis of REG signals was based on geometric features extracted from the time domain in the first place, since this is the standard processing strategy for this type of physiological data. Geometric features were tested to distinguish between different anesthetic depths, and they proved to be capable of tracking cerebral hemodynamic changes during anesthesia. Furthermore, an approach based on Poincaré plot features was proposed, where the reconstructed attractors form REG signals showed significant differences between different anesthetic states. This was a key finding, providing an alternative to standard processing of REG signals and supporting the hypothesis that REG signals do carry CBF information. Furthermore, the analysis of cerebral hemodynamics during anesthetic procedures was performed by means of studying causal relationships between global hemodynamics, cerebral hemodynamics and electroencephalogram (EEG) based-parameters. Interactions were detected during anesthetic drug infusion and patient positioning (Trendelenburg positioning and passive leg raise), providing evidence of the causal coupling between hemodynamics and brain activity. The provided alternative of REG signal processing confirmed the hypothesis that REG signals carry information on CBF. The simplicity of the technology, together with its low cost and easily interpretable outcomes, should provide a new opportunity for REG to reach standard clinical practice. Moreover, causal relationships among the hemodynamic physiological signals and brain activity were assessed, suggesting that the inclusion of REG information in depth of anesthesia monitors could be of valuable use to prevent unwanted CBF alterations during anesthetic procedures.

Disruption of the Cortical-Vagal Communication Network in Parkinson's Disease

Parkinson's disease (PD) is a neuropathy characterized by motor disorders, but it has also been associated with the presence of autonomic alterations as a result of degradation of the dopaminergic system. Studying the relation between Band Power time series (BPts) and Heart Rate Variability (HRV), has been proposed as a tool to explore the bidirectional communication pathways between cortex and autonomic control. This work presents a primer analysis on study brain ↔ heart interaction on a databse of PD patients under two conditions: without and after levadopa (L-dopa) intake. Additionally a healthy control population was also analyzed, and used as comparison level between both conditions. Results show PD affects pathways by reducing the number of connections, specially association of beta and power and the second faster component of HRV seems to be more sensitive to L-dopa administration.

Multivariate Correlation Measures Reveal Structure and Strength of Brain-Body Physiological Networks at Rest and During Mental Stress

In this work, we extend to the multivariate case the classical correlation analysis used in the field of network physiology to probe dynamic interactions between organ systems in the human body. To this end, we define different correlation-based measures of the multivariate interaction (MI) within and between the brain and body subnetworks of the human physiological network, represented, respectively, by the time series of δ, θ, α, and β electroencephalographic (EEG) wave amplitudes, and of heart rate, respiration amplitude, and pulse arrival time (PAT) variability (η, ρ, π). MI is computed: (i) considering all variables in the two subnetworks to evaluate overall brain-body interactions; (ii) focusing on a single target variable and dissecting its global interaction with all other variables into contributions arising from the same subnetwork and from the other subnetwork; and (iii) considering two variables conditioned to all the others to infer the network topology. The framework is applied to the time series measured from the EEG, electrocardiographic (ECG), respiration, and blood volume pulse (BVP) signals recorded synchronously via wearable sensors in a group of healthy subjects monitored at rest and during mental arithmetic and sustained attention tasks. We find that the human physiological network is highly connected, with predominance of the links internal of each subnetwork (mainly η-ρ and δ-θ, θ-α, α-β), but also statistically significant interactions between the two subnetworks (mainly η-β and η-δ). MI values are often spatially heterogeneous across the scalp and are modulated by the physiological state, as indicated by the decrease of cardiorespiratory interactions during sustained attention and by the increase of brain-heart interactions and of brain-brain interactions at the frontal scalp regions during mental arithmetic. These findings illustrate the complex and multi-faceted structure of interactions manifested within and between different physiological systems and subsystems across different levels of mental stress.

Autonomic modulation networks in schizophrenia: The relationship between heart rate variability and functional and structural connectivity in the brain

Heart rate variability (HRV), a measurement of autonomic nervous system (ANS) activity, has been found reduced in schizophrenia. The anterior cingulate cortex (ACC), which is important in regulating the ANS, is structurally and functionally affected in schizophrenia. We investigate the relationship between HRV and functional and structural connectivity of the ACC in patients with schizophrenia and healthy controls. Ten patients with a diagnosis of schizophrenia and ten healthy controls were recruited. Heart rate was monitored in a naturalistic out-of-clinic setting. Magnetic resonance imaging (MRI) was performed, including resting-state functional MRI and diffusion tensor imaging. Patients with schizophrenia had significantly lower HRV compared to controls. A positive correlation between ACC connectivity with the bilateral cerebellum and HRV was found in the patients. HRV was also positively correlated with amplitude of low frequency fluctuations (ALFF) in the cerebellum, and with axial diffusivity in the middle cerebellar peduncle, in the patients. There was a significant negative relationship between antipsychotic medication dosage, HRV and all neuroimaging measures related to HRV. We conclude that ACC connectivity seems to be affected in schizophrenia, both structurally and functionally, and that the ACC-cerebellum connectivity, as well as cerebellar function, is associated with ANS regulation in patients with schizophrenia.

Signal Information Prediction of Mortality Identifies Unique Patient Subsets after Severe Traumatic Brain Injury: A Decision-Tree Analysis Approach

Nonlinear physiological signal features that reveal information content and causal flow have recently been shown to be predictors of mortality after severe traumatic brain injury (TBI). The extent to which these features interact together, and with traditional measures to describe patients in a clinically meaningful way remains unclear. In this study, we incorporated basic demographics (age and initial Glasgow Coma Scale [GCS]) with linear and non-linear signal information based features (approximate entropy [ApEn], and multivariate conditional Granger causality [GC]) to evaluate their relative contributions to mortality using cardio-cerebral monitoring data from 171 severe TBI patients admitted to a single neurocritical care center over a 10 year period. Beyond linear modelling, we employed a decision tree analysis approach to define a predictive hierarchy of features. We found ApEn (p = 0.009) and GC (p = 0.004) based features to be independent predictors of mortality at a time when mean intracranial pressure (ICP) was not. Our combined model with both signal information-based features performed the strongest (area under curve = 0.86 vs. 0.77 for linear features only). Although low "intracranial" complexity (ApEn-ICP) outranked both age and GCS as crucial drivers of mortality (fivefold increase in mortality where ApEn-ICP <1.56, 36.2% vs. 7.8%), decision tree analysis revealed clear subsets of patient populations using all three predictors. Patients with lower ApEn-ICP who were >60 years of age died, whereas those with higher ApEn-ICP and GCS ≥5 all survived. Yet, even with low initial intracranial complexity, as long as patients maintained robust GC and "extracranial" complexity (ApEn of mean arterial pressure), they all survived. Incorporating traditional linear and novel, non-linear signal information features, particularly in a framework such as decision trees, may provide better insight into "health" status. However, caution is required when interpreting these results in a clinical setting prior to external validation.

Uncovering complex central autonomic networks at rest: a functional magnetic resonance imaging study on complex cardiovascular oscillations

This study aims to uncover brain areas that are functionally linked to complex cardiovascular oscillations in resting-state conditions. Multi-session functional magnetic resonance imaging (fMRI) and cardiovascular data were gathered from 34 healthy volunteers recruited within the human connectome project (the '100-unrelated subjects' release). Group-wise multi-level fMRI analyses in conjunction with complex instantaneous heartbeat correlates (entropy and Lyapunov exponent) revealed the existence of a specialized brain network, i.e. a complex central autonomic network (CCAN), reflecting what we refer to as complex autonomic control of the heart. Our results reveal CCAN areas comprised the paracingulate and cingulate gyri, temporal gyrus, frontal orbital cortex, planum temporale, temporal fusiform, superior and middle frontal gyri, lateral occipital cortex, angular gyrus, precuneous cortex, frontal pole, intracalcarine and supracalcarine cortices, parahippocampal gyrus and left hippocampus. The CCAN visible at rest does not include the insular cortex, thalamus, putamen, amygdala and right caudate, which are classical CAN regions peculiar to sympatho-vagal control. Our results also suggest that the CCAN is mainly involved in complex vagal control mechanisms, with possible links with emotional processing networks.

Resting State Neural Correlates of Cardiac Sympathetic Dynamics in Healthy Subjects

Recent advances in functional Magnetic Resonance Imaging (fMRI) research have uncovered the existence of the central autonomic network (CAN), which comprises brain regions whose activity correlates with autonomic nervous system dynamics. By exploiting the spectral paradigm of heartbeat dynamics, cortical and sub-cortical areas functionally linked to vagal activity have been identified. However, due to methodological limitations, functional neural correlates of cardiac sympathetic dynamics remain uncharacterized. To this extent, we exploit the high spatiotemporal resolution of fMRI data from the Human Connectome Project to study the CAN activity by correlating a recently proposed instantaneous characterization of sympathetic activity (the sympathetic activity index - SAI) from heartbeat dynamics. SAI estimates are embedded into the probabilistic modeling of inhomogeneous point-processes, and are derived from a combination of disentangling coefficients linked to the orthonormal Laguerre functions. By analyzing resting state recordings from 34 young healthy people, we obtain positive correlations between instantaneous SAI estimates and a number of brain regions including frontal pole, insular cortex, frontal and temporal gyri, lateral occipital cortex, paracingulate and cingulate gyri, precuneus and temporal fusiform cortices, as well as thalamus, caudate nucleus, putamen, brain-stem, hippocampus, amygdala, and nucleus accumbens. Our findings significantly extend current knowledge on the CAN, opening new avenues in the characterization of healthy and pathological states in humans.

The central autonomic network at rest: Uncovering functional MRI correlates of time-varying autonomic outflow

Peripheral measures of autonomic nervous system (ANS) activity at rest have been extensively employed as putative biomarkers of autonomic cardiac control. However, a comprehensive characterization of the brain-based central autonomic network (CAN) sustaining cardiovascular oscillations at rest is missing, limiting the interpretability of these ANS measures as biomarkers of cardiac control. We evaluated combined cardiac and fMRI data from 34 healthy subjects from the Human Connectome Project to detect brain areas functionally linked to cardiovagal modulation at rest. Specifically, we combined voxel-wise fMRI analysis with instantaneous heartbeat and spectral estimates obtained from inhomogeneous linear point-process models. We found exclusively negative associations between cardiac parasympathetic activity at rest and a widespread network including bilateral anterior insulae, right dorsal middle and left posterior insula, right parietal operculum, bilateral medial dorsal and ventrolateral posterior thalamic nuclei, anterior and posterior mid-cingulate cortex, medial frontal gyrus/pre-supplementary motor area. Conversely, we found only positive associations between instantaneous heart rate and brain activity in areas including frontopolar cortex, dorsomedial prefrontal cortex, anterior, middle and posterior cingulate cortices, superior frontal gyrus, and precuneus. Taken together, our data suggests a much wider involvement of diverse brain areas in the CAN at rest than previously thought, which could reflect a differential (both spatially and directionally) CAN activation according to the underlying task. Our insight into CAN activity at rest also allows the investigation of its impairment in clinical populations in which task-based fMRI is difficult to obtain (e.g., comatose patients or infants).

Heart rate variability as a biomarker in health and affective disorders: A perspective on neuroimaging studies

The dynamic embodiment of psychological processes is evident in the association of health outcomes, behavioural traits and psychological functioning with Heart Rate Variability (HRV). The dominant high-frequency component of HRV is an index of the central neural control of heart rhythm, mediated via the parasympathetic vagus nerve. HRV provides a potential objective measure of action policies for the adaptive and predictive allostatic regulation of homeostasis within the cardiovascular system. In its support, a network of brain regions (referred to as the 'central autonomic network') maps internal state, and controls autonomic responses. This network includes regions of prefrontal cortex, anterior cingulate cortex, insula, amygdala, periaqueductal grey, pons and medulla. Human neuroimaging studies of neural activation and functional connectivity broadly endorse this architecture, and its link with cardiac regulation at rest and dysregulation in clinical states that include affective disorders. In this review, we appraise neuroimaging research and related evidence for HRV as an informative marker of autonomic integration with affect and cognition, taking a perspective on function and organisation. We consider evidence for the utility of HRV as a metric to inform targeted interventions to improve autonomic and affective dysregulation, and suggest research questions for further investigation.

Variability and Reproducibility of Directed and Undirected Functional MRI Connectomes in the Human Brain

A growing number of studies are focusing on methods to estimate and analyze the functional connectome of the human brain. Graph theoretical measures are commonly employed to interpret and synthesize complex network-related information. While resting state functional MRI (rsfMRI) is often employed in this context, it is known to exhibit poor reproducibility, a key factor which is commonly neglected in typical cohort studies using connectomics-related measures as biomarkers. We aimed to fill this gap by analyzing and comparing the inter- and intra-subject variability of connectivity matrices, as well as graph-theoretical measures, in a large (n = 1003) database of young healthy subjects which underwent four consecutive rsfMRI sessions. We analyzed both directed (Granger Causality and Transfer Entropy) and undirected (Pearson Correlation and Partial Correlation) time-series association measures and related global and local graph-theoretical measures. While matrix weights exhibit a higher reproducibility in undirected, as opposed to directed, methods, this difference disappears when looking at global graph metrics and, in turn, exhibits strong regional dependence in local graphs metrics. Our results warrant caution in the interpretation of connectivity studies, and serve as a benchmark for future investigations by providing quantitative estimates for the inter- and intra-subject variabilities in both directed and undirected connectomic measures.

Disruption of brain network organization in primary open angle glaucoma

Resting-state functional magnetic resonance imaging (rs-fMRI) is commonly employed to study changes in functional brain connectivity. Recently, the hypothesis of a brain involvement in primary open angle glaucoma has sprung interest for neuroimaging studies in this pathology. The purpose of this study is to evaluate a putative reorganization of brain networks in glaucomatous patients through graph-theoretical measures of integration, segregation and centrality by exploiting a multivariate networks association measure and a recently introduced global and local brain network disruption index. Nineteen glaucoma patients and sixteen healthy control subjects (age: 50 - 76, mean 61 years) underwent rs-fMRI examination at 3T. After preprocessing, rs-fMRI time series were parcellated into 116 regions (AAL atlas), adjacency matrices were computed based on partial correlations and graph-theoretical measures of integration, segregation and centrality as well as group-wise and subject-wise disruption index estimates were generated for all subjects. We found that the group-wise disruption index was negative and statistically different from 0 in for all graph theoretical metrics. Additionally, statistically significant group-wise differences in subject-wise disruption indexes were found in all local metrics. The differences in local network measures highlight cerebral reorganization of brain networks in glaucoma patients, supporting the interpretation of glaucoma as central nervous system disease, likely part of the heterogeneous group of recently described disconnection syndromes.

Network physiology of 'fight or flight' response in facial superficial blood vessels

OBJECTIVE

We introduced a novel framework to identify the dynamic pattern of blood flow changes in the cutaneous superficial blood vessels of the face for 'fight or flight' responses through facial thermal imaging.

APPROACH

For this purpose, a thermal dataset was collected from 41 subjects in a mock crime scenario. Five facial areas including periorbital, forehead, perinasal, cheek and chin were selected on the face. Due to the cause and effect movement of blood in the facial cutaneous vasculature, the effective connectivity approach and graph analysis were used to extract causality features. The effective connectivity was quantified using a modified version of the multivariate Granger causality (GC) method among each pair of facial region of interests.

MAIN RESULTS

Validation was performed using statistical analysis, and the results demonstrated that the proposed method was statistically significant in detecting the physiological pattern of deceptive anxiety on the face. Moreover, the obtained graph is visualized by different schemes to show these interactions more effectively. We used machine learning techniques to classify our data based on the GC values, which result in a greater than 87% accuracy rate in discriminating between deceptive and truthful subjects.

Functional neuroimaging of the central autonomic network: recent developments and clinical implications

Purpose

The central autonomic network (CAN) is an intricate system of brainstem, subcortical, and cortical structures that play key roles in the function of the autonomic nervous system. Prior to the advent of functional neuroimaging, in vivo studies of the human CAN were limited. The purpose of this review is to highlight the contribution of functional neuroimaging, specifically functional magnetic resonance imaging (fMRI), to the study of the CAN, and to discuss recent advances in this area. Additionally, we aim to emphasize exciting areas for future research.

Methods

We reviewed the existing literature in functional neuroimaging of the CAN. Here, we focus on fMRI research conducted in healthy human subjects, as well as research that has been done in disease states, to understand CAN function. To minimize confounding, papers examining CAN function in the context of cognition, emotion, pain, and affective disorders were excluded.

Results

fMRI has led to significant advances in the understanding of human CAN function. The CAN is composed of widespread brainstem and forebrain structures that are intricately connected and play key roles in reflexive and modulatory control of autonomic function.

Conclusions

fMRI technology has contributed extensively to current knowledge of CAN function. It holds promise to serve as a biomarker in disease states. With ongoing advancements in fMRI technology, there is great opportunity and need for future research involving the CAN.

A realistic neuronal network and neurovascular coupling model for the study of multivariate directed connectivity in fMRI data

The use of Multivariate Granger Causality (MVGC) in estimating directed Blood-Oxygen-Level- Dependant (BOLD) connectivity is still controversial. This is mostly due to the short data Ienghts typically available in func- tional MRI (fMRI) acquisitions, to the very nature of the BOLD acquisition strategy (which yields extremely low signal- to-noise-ratio) and importantly to the fact that neuronal activi- ty is convolved with a slow-varying haemodynamic response function (HRF) which therefore generates a temporal confound which is arduous to account for when basing MVGC estimates on vector autoregressive models (VAR). In this paper, we em- ploy realistic complex network models based on Izhikevich neuronal populations, interlinked by realistic neuronal fiber bundles which exert compounded directed influences and cas- cade into Baloon-model-like neurovascular coupling, to explore and validate the MVGC approach to directed connectivity es- timation in realistic fMRI conditions and in a complex directed network setting. In particular, we show in silico that the top 1 percentile of a BOLD connectivity matrix estimated with MVGC from BOLD data similar to the one provided by the Human Connectome Project (HCP) has a Positive Predictive Value very close to 1, hence corroborating the evidence that the "strongest" connections can be safely studied with this method in fMRI.

A Heartbeat Away From Consciousness: Heart Rate Variability Entropy Can Discriminate Disorders of Consciousness and Is Correlated With Resting-State fMRI Brain Connectivity of the Central Autonomic Network

Background: Disorders of consciousness are challenging to diagnose, with inconsistent behavioral responses, motor and cognitive disabilities, leading to approximately 40% misdiagnoses. Heart rate variability (HRV) reflects the complexity of the heart-brain two-way dynamic interactions. HRV entropy analysis quantifies the unpredictability and complexity of the heart rate beats intervals. We here investigate the complexity index (CI), a score of HRV complexity by aggregating the non-linear multi-scale entropies over a range of time scales, and its discriminative power in chronic patients with unresponsive wakefulness syndrome (UWS) and minimally conscious state (MCS), and its relation to brain functional connectivity. Methods: We investigated the CI in short (CIs) and long (CIl) time scales in 14 UWS and 16 MCS sedated. CI for MCS and UWS groups were compared using a Mann-Whitney exact test. Spearman's correlation tests were conducted between the Coma Recovery Scale-revised (CRS-R) and both CI. Discriminative power of both CI was assessed with One-R machine learning model. Correlation between CI and brain connectivity (detected with functional magnetic resonance imagery using seed-based and hypothesis-free intrinsic connectivity) was investigated using a linear regression in a subgroup of 10 UWS and 11 MCS patients with sufficient image quality. Results: Higher CIs and CIl values were observed in MCS compared to UWS. Positive correlations were found between CRS-R and both CI. The One-R classifier selected CIl as the best discriminator between UWS and MCS with 90% accuracy, 7% false positive and 13% false negative rates after a 10-fold cross-validation test. Positive correlations were observed between both CI and the recovery of functional connectivity of brain areas belonging to the central autonomic networks (CAN). Conclusion: CI of MCS compared to UWS patients has high discriminative power and low false negative rate at one third of the estimated human assessors' misdiagnosis, providing an easy, inexpensive and non-invasive diagnostic tool. CI reflects functional connectivity changes in the CAN, suggesting that CI can provide an indirect way to screen and monitor connectivity changes in this neural system. Future studies should assess the extent of CI's predictive power in a larger cohort of patients and prognostic power in acute patients.

Dynamic inter-network connectivity in the human brain

Recently, the field of functional brain connectivity has shifted its attention on studying how functional connectivity (FC) between remote regions changes over time. It is becoming increasingly evident that the human "connectome" is a dynamical entity whose variations are effected over very short timescales and reflect crucial mechanisms which underline the physiological functioning of the brain. In this study, we employ ad-hoc statistical and surrogate data generation methods to quantify whether and which brain networks displayed dynamic behaviors in a very large sample of healthy subjects provided by the Human Connectome Project (HCP). Our findings provided evidences that there are specific pairs of networks and specific networks within the healthy brain that are more likely to display dynamic behaviors. This new set of findings supports the notion that studying the time-variant connectivity in the brain could reveal useful and important properties about brain functioning in health and disease.

Resting-state brain correlates of instantaneous autonomic outflow

A prominent pathway of brain-heart interaction is represented by autonomic nervous system (ANS) heartbeat modulation. While within-brain resting state networks have been the object of intense functional Magnetic Resonance Imaging (fMRI) research, technological and methodological limitations have hampered research on the central correlates of cardiovascular control dynamics. Here we combine the high temporal and spatial resolution as well as data volume afforded by the Human Connectome Project with a probabilistic model of heartbeat dynamics to characterize central correlates of sympathetic and parasympathetic ANS activity at rest. We demonstrate an involvement of a number of brain regions such as the Insular cortex, Frontal Gyrus, Lateral Occipital Cortex, Paracingulate and Cingulate Gyrus and Precuneous Cortex, as well as subcortical structures (Thalamus, Putamen, Pallidum, Brain-Stem, Hippocampus, Amygdala, and Right Caudate) in the modulation of ANS-mediated cardiovascular control, possibly indicating a broader definition of the central autonomic network (CAN). Our findings provide a basis for an informed neurobiological interpretation of the numerous studies which employ HRV-related measures as standalone biomarkers in health and disease.

Resting-state brain correlates of cardiovascular complexity

While estimates of complex heartbeat dynamics have provided effective prognostic and diagnostic markers for a wide range of pathologies, brain correlates of complex cardiac measures in general and of complex sympatho-vagal dynamics in particular are still unknown. In this study we combine resting state functional Magnetic Resonance Imaging (fMRI) and physiological signal acquisition from 34 healthy subjects selected from the Human Connectome Project (HCP) repository with inhomogeneous point-process approximate and sample heartbeat entropy measures (ipApEn and ipSampEn) to investigate brain areas involved in complex cardiovascular control. Our results show that activity in the Temporal Gyrus, Frontal Orbital Cortex, Temporal Fusiform and Opercular cortices, Planum Temporale, and Paracingulate cortex are negatively correlated with ipApEn dynamics. Activity in the same cortical areas as well as in the Temporal Fusiform cortex are negatively correlated with ipSampEn dynamics. No significant positive correlations were found. These pioneering results suggest that cardiovascular complexity at rest is linked to a few specific cortical brain structures, including crucial areas connected with parasympathetic outflow. This corroborates the hypothesis of a multidimensional central network which controls nonlinear cardiac dynamics under a predominantly vagally-driven tone.

Dynamical brain connectivity estimation using GARCH models: An application to personality neuroscience

It has recently become evident that the functional connectome of the human brain is a dynamical entity whose time evolution carries important information underpinning physiological brain function as well as its disease-related aberrations. While simple sliding window approaches have had some success in estimating dynamical brain connectivity in a functional MRI (fMRI) context, these methods suffer from limitations related to the arbitrary choice of window length and limited time resolution. Recently, Generalized autoregressive conditional heteroscedastic (GARCH) models have been employed to generate dynamical covariance models which can be applied to fMRI. Here, we employ a GARCH-based method (dynamic conditional correlation - DCC) to estimate dynamical brain connectivity in the Human Connectome Project (HCP) dataset and study how the dynamic functional connectivity behaviors related to personality as described by the five-factor model. Openness, a trait related to curiosity and creativity, is the only trait associated with significant differences in the amount of time-variability (but not in absolute median connectivity) of several inter-network functional connections in the human brain. The DCC method offers a novel window to extract dynamical information which can aid in elucidating the neurophysiological underpinning of phenomena to which conventional static brain connectivity estimates are insensitive.

Sparse Estimation of Resting-State Effective Connectivity From fMRI Cross-Spectra

In functional magnetic resonance imaging (fMRI), functional connectivity is conventionally characterized by correlations between fMRI time series, which are intrinsically undirected measures of connectivity. Yet, some information about the directionality of network connections can nevertheless be extracted from the matrix of pairwise temporal correlations between all considered time series, when expressed in the frequency-domain as a cross-spectral density matrix. Using a sparsity prior, it then becomes possible to determine a unique directed network topology that best explains the observed undirected correlations, without having to rely on temporal precedence relationships that may not be valid in fMRI. Applying this method on simulated data with 100 nodes yielded excellent retrieval of the underlying directed networks under a wide variety of conditions. Importantly, the method did not depend on temporal precedence to establish directionality, thus reducing susceptibility to hemodynamic variability. The computational efficiency of the algorithm was sufficient to enable whole-brain estimations, thus circumventing the problem of missing nodes that otherwise occurs in partial-brain analyses. Applying the method to real resting-state fMRI data acquired with a high temporal resolution, the inferred networks showed good consistency with structural connectivity obtained from diffusion tractography in the same subjects. Interestingly, this agreement could also be seen when considering high-frequency rather than low-frequency connectivity (average correlation: r = 0.26 for f < 0.3 Hz, r = 0.43 for 0.3 < f < 5 Hz). Moreover, this concordance was significantly better (p < 0.05) than for networks obtained with conventional functional connectivity based on correlations (average correlation r = 0.18). The presented methodology thus appears to be well-suited for fMRI, particularly given its lack of explicit dependence on temporal lag structure, and is readily applicable to whole-brain effective connectivity estimation.

Task-induced deactivation in diverse brain systems correlates with interindividual differences in distinct autonomic indices

Neuroimaging research has shown that different cognitive tasks induce relatively specific activation patterns, as well as less task-specific deactivation patterns. Here we examined whether individual differences in Autonomic Nervous System (ANS) activity during task performance correlate with the magnitude of task-induced deactivation. In an fMRI study, participants performed a continuous mental arithmetic task in a task/rest block design, while undergoing combined fMRI and heart/respiration rate acquisitions using photoplethysmograph and respiration belt. As expected, task performance increased heart-rate and reduced the RMSSD, a cardiac index related to vagal tone. Across participants, higher heart rate during task was linked to increased activation in fronto-parietal regions, as well as to stronger deactivation in ventromedial prefrontal regions. Respiration frequency during task was associated with similar patterns, but in different regions than those identified for heart-rate. Finally, in a large set of regions, almost exclusively limited to the Default Mode Network, lower RMSSD was associated with greater deactivation, and furthermore, the vast majority of these regions were task-deactivated at the group level. Together, our findings show that inter-individual differences in ANS activity are strongly linked to task-induced deactivation. Importantly, our findings suggest that deactivation is a multifaceted construct potentially linked to ANS control, because distinct ANS measures correlate with deactivation in different regions. We discuss the implications for current theories of cortical control of the ANS and for accounts of deactivation, with particular reference to studies documenting a "failure to deactivate" in multiple clinical states.

Revolution of Alzheimer Precision Neurology. Passageway of Systems Biology and Neurophysiology

The Precision Neurology development process implements systems theory with system biology and neurophysiology in a parallel, bidirectional research path: a combined hypothesis-driven investigation of systems dysfunction within distinct molecular, cellular, and large-scale neural network systems in both animal models as well as through tests for the usefulness of these candidate dynamic systems biomarkers in different diseases and subgroups at different stages of pathophysiological progression. This translational research path is paralleled by an "omics"-based, hypothesis-free, exploratory research pathway, which will collect multimodal data from progressing asymptomatic, preclinical, and clinical neurodegenerative disease (ND) populations, within the wide continuous biological and clinical spectrum of ND, applying high-throughput and high-content technologies combined with powerful computational and statistical modeling tools, aimed at identifying novel dysfunctional systems and predictive marker signatures associated with ND. The goals are to identify common biological denominators or differentiating classifiers across the continuum of ND during detectable stages of pathophysiological progression, characterize systems-based intermediate endophenotypes, validate multi-modal novel diagnostic systems biomarkers, and advance clinical intervention trial designs by utilizing systems-based intermediate endophenotypes and candidate surrogate markers. Achieving these goals is key to the ultimate development of early and effective individualized treatment of ND, such as Alzheimer's disease. The Alzheimer Precision Medicine Initiative (APMI) and cohort program (APMI-CP), as well as the Paris based core of the Sorbonne University Clinical Research Group "Alzheimer Precision Medicine" (GRC-APM) were recently launched to facilitate the passageway from conventional clinical diagnostic and drug development toward breakthrough innovation based on the investigation of the comprehensive biological nature of aging individuals. The APMI movement is gaining momentum to systematically apply both systems neurophysiology and systems biology in exploratory translational neuroscience research on ND.

Simultaneous estimation of the in-mean and in-variance causal connectomes of the human brain

In recent years, the study of the human connectome (i.e. of statistical relationships between non spatially contiguous neurophysiological events in the human brain) has been enormously fuelled by technological advances in high-field functional magnetic resonance imaging (fMRI) as well as by coordinated world wide data-collection efforts like the Human Connectome Project (HCP). In this context, Granger Causality (GC) approaches have recently been employed to incorporate information about the directionality of the influence exerted by a brain region on another. However, while fluctuations in the Blood Oxygenation Level Dependent (BOLD) signal at rest also contain important information about the physiological processes that underlie neurovascular coupling and associations between disjoint brain regions, so far all connectivity estimation frameworks have focused on central tendencies, hence completely disregarding so-called in-variance causality (i.e. the directed influence of the volatility of one signal on the volatility of another). In this paper, we develop a framework for simultaneous estimation of both in-mean and in-variance causality in complex networks. We validate our approach using synthetic data from complex ensembles of coupled nonlinear oscillators, and successively employ HCP data to provide the very first estimate of the in-variance connectome of the human brain.

Estimating directed brain-brain and brain-heart connectivity through globally conditioned Granger causality approaches

While a large body of research has focused on the study of within-brain physiological networks (i.e. brain connectivity) as well as their disease-related aberration, few investigators have focused on estimating the directionality of these brain-brain interaction which, given the complexity of brain networks, should be properly conditioned in order to avoid the high number of false positives commonly encountered when using bivariate approaches to brain connectivity estimation. Additionally, the constituents of a number of brain subnetworks, and in particular of the central autonomic network (CAN), are still not completely determined. In this study we present and validate a global conditioning approach to reconstructing directed networks using complex synthetic networks of nonlinear oscillators. We then employ our framework, along with a probabilistic model for heartbeat generation, to characterize the directed functional connectome of the human brain and to establish which parts of this connectome effect the directed central modulation of peripheral autonomic cardiovascular control. We demonstrate the effectiveness of our conditioning approach and unveil a top-down directed influence of the default mode network on the salience network, which in turn is seen to be the strongest modulator of directed autonomic cardiovascular control.

Early Asymmetric Cardio-Cerebral Causality and Outcome after Severe Traumatic Brain Injury

The brain and heart are two vital systems in health and disease, increasingly recognized as a complex, interdependent network with constant information flow in both directions. After severe traumatic brain injury (TBI), the causal, directed interactions between the brain, heart, and autonomic nervous system have not been well established. Novel methods are needed to probe unmeasured, potentially prognostic information in complex biological networks that are not revealed by traditional means. In this study, we examined potential bidirectional causality between intracranial pressure (ICP), mean arterial pressure (MAP), and heart rate (HR) and its relationship to mortality in a 24-h period early post-TBI. We applied Granger causality (GC) analysis to cardio-cerebral monitoring data from 171 severe TBI patients admitted to a single neurocritical care center over a 10-year period. There was significant bidirectional causality between ICP and MAP, MAP and HR, and ICP and HR in the majority of patients (p < 0.01). MAP influenced both ICP and HR to a greater extent (higher GC, p < 0. 00001), but there was no dominant unidirectional causality between ICP and HR (p = 0.85). Those who died had significantly lower GC for ICP causing MAP and HR causing ICP (p = 0.006 and p = 0.004, respectively) and were predictors of mortality independent of age, sex, and traditional intracranial variables (ICP, cerebral perfusion pressure, GCS, and pressure reactivity index). Examining the brain and heart with GC-based features for the first time in severe TBI patients has confirmed strong interdependence and reveals a significant relationship between select causality pairs and mortality. These results support the notion that impaired causal information flow between the cerebrovascular, autonomic, and cardiovascular systems are of central importance in severe TBI.

Functional connectivity in amygdalar-sensory/(pre)motor networks at rest: new evidence from the Human Connectome Project

The word 'e-motion' derives from the Latin word 'ex-moveo' which literally means 'moving away from something/somebody'. Emotions are thus fundamental to prime action and goal-directed behavior with obvious implications for individual's survival. However, the brain mechanisms underlying the interactions between emotional and motor cortical systems remain poorly understood. A recent diffusion tensor imaging study in humans has reported the existence of direct anatomical connections between the amygdala and sensory/(pre)motor cortices, corroborating an initial observation in animal research. Nevertheless, the functional significance of these amygdala-sensory/(pre)motor pathways remain uncertain. More specifically, it is currently unclear whether a distinct amygdala-sensory/(pre)motor circuit can be identified with resting-state functional magnetic resonance imaging (rs-fMRI). This is a key issue, as rs-fMRI offers an opportunity to simultaneously examine distinct neural circuits that underpin different cognitive, emotional and motor functions, while minimizing task-related performance confounds. We therefore tested the hypothesis that the amygdala and sensory/(pre)motor cortices could be identified as part of the same resting-state functional connectivity network. To this end, we examined independent component analysis results in a very large rs-fMRI data-set drawn from the Human Connectome Project (n = 820 participants, mean age: 28.5 years). To our knowledge, we report for the first time the existence of a distinct amygdala-sensory/(pre)motor functional network at rest. rs-fMRI studies are now warranted to examine potential abnormalities in this circuit in psychiatric and neurological diseases that may be associated with alterations in the amygdala-sensory/(pre)motor pathways (e.g. conversion disorders, impulse control disorders, amyotrophic lateral sclerosis and multiple sclerosis).

Reconstructing multivariate causal structure between functional brain networks through a Laguerre-Volterra based Granger causality approach

Classical multivariate approaches based on Granger causality (GC) which estimate functional connectivity in the brain are almost exclusively based on autoregressive models. Nevertheless, information available from past samples is limited due to both signal autocorrelation and necessarily low model orders. Consequently, multiple time-scales interactions are usually unaccounted for. To overcome these limitations, in this study we propose the use of discrete-time orthogonal Laguerre basis functions within a Wiener-Volterra decomposition of the BOLD signals to perform effective GC assessments of brain functional connectivity. We validate our method in synthetic noisy oscillator networks, and analyze experimental fMRI data from 30 healthy subjects publicly available within the Human Connectome Project (HCP). Synthetic results demonstrate that our Laguerre-Volterra based GC estimates outperform classical approaches in terms of accuracy in detecting true causal links while rejecting false causal links in complex nonlinear networks. Human data analysis shows for the first time that the default mode network modulates both the salience network as well as fronto-temporal circuits in a causal fashion.

Uncovering brain-heart information through advanced signal and image processing

Through their dynamical interplay, the brain and the heart ensure fundamental homeostasis and mediate a number of physiological functions as well as their disease-related aberrations. Although a vast number of ad hoc analytical and computational tools have been recently applied to the non-invasive characterization of brain and heart dynamic functioning, little attention has been devoted to combining information to unveil the interactions between these two physiological systems. This theme issue collects contributions from leading experts dealing with the development of advanced analytical and computational tools in the field of biomedical signal and image processing. It includes perspectives on recent advances in 7 T magnetic resonance imaging as well as electroencephalogram, electrocardiogram and cerebrovascular flow processing, with the specific aim of elucidating methods to uncover novel biological and physiological correlates of brain-heart physiology and physiopathology.