Causal interactions in brain networks predict pain levels in trigeminal neuralgia

Trigeminal neuralgia (TN) is a highly debilitating facial pain condition. Magnetic resonance imaging (MRI) is the main method for generating insights into the central mechanisms of TN pain in humans. Studies have found both structural and functional abnormalities in various brain structures in TN patients as compared with healthy controls. Whereas studies have also examined aberrations in brain networks in TN, no studies have to date investigated causal interactions in these brain networks and related these causal interactions to the levels of TN pain. We recorded fMRI data from 39 TN patients who either rested comfortably in the scanner during the resting state session or tracked their pain levels during the pain tracking session. Applying Granger causality to analyze the data and requiring consistent findings across the two scanning sessions, we found 5 causal interactions, including: (1) Thalamus → dACC, (2) Caudate → Inferior temporal gyrus, (3) Precentral gyrus → Inferior temporal gyrus, (4) Supramarginal gyrus → Inferior temporal gyrus, and (5) Bankssts → Inferior temporal gyrus, that were consistently associated with the levels of pain experienced by the patients. Utilizing these 5 causal interactions as predictor variables and the pain score as the predicted variable in a linear multiple regression model, we found that in both pain tracking and resting state sessions, the model was able to explain ∼36 % of the variance in pain levels, and importantly, the model trained on the 5 causal interaction values from one session was able to predict pain levels using the 5 causal interaction values from the other session, thereby cross-validating the models. These results, obtained by applying novel analytical methods to neuroimaging data, provide important insights into the pathophysiology of TN and could inform future studies aimed at developing innovative therapies for treating TN.

An information-theoretic quantification of the content of communication between brain regions

Quantifying the amount, content and direction of communication between brain regions is key to understanding brain function. Traditional methods to analyze brain activity based on the Wiener-Granger causality principle quantify the overall information propagated by neural activity between simultaneously recorded brain regions, but do not reveal the information flow about specific features of interest (such as sensory stimuli). Here, we develop a new information theoretic measure termed Feature-specific Information Transfer (FIT), quantifying how much information about a specific feature flows between two regions. FIT merges the Wiener-Granger causality principle with information-content specificity. We first derive FIT and prove analytically its key properties. We then illustrate and test them with simulations of neural activity, demonstrating that FIT identifies, within the total information flowing between regions, the information that is transmitted about specific features. We then analyze three neural datasets obtained with different recording methods, magneto- and electro-encephalography, and spiking activity, to demonstrate the ability of FIT to uncover the content and direction of information flow between brain regions beyond what can be discerned with traditional anaytical methods. FIT can improve our understanding of how brain regions communicate by uncovering previously hidden feature-specific information flow.

Pathway-resolved decomposition demonstrates correlation and noise dependencies of redundant information processing in recurrent feed-forward topologies

In a biochemical assay that converts fan-in networks into feed-forward loops (FFLs), we show that the inter-regulator redundant information about the output gene product can be decomposed into finer components, mediated by the constituent pathways. Variance-based information within the linear noise regime facilitates quantifying these submodular redundancies. Contrary to the conventional wisdom on information decomposition, we report that information redundancy depends nontrivially on inter-regulator correlation. For the type-1 coherent (C1) and incoherent (I1) FFLs, the direct regulatory path-mediated redundancy is certainly correlation independent. However, components induced by the indirect regulatory path and interpathway interference are correlation dependent in (non)linear fashion. The trade-off between information redundancy and similarly decomposable extrinsic noise from input to output node has been demonstrated for the pathways and full motifs. Our analyses suggest that the interpathway cross redundancy positively and negatively influences the superposition of elementary redundancies in the C1- and I1-FFLs, respectively. Their corresponding total extrinsic noise is produced by the weighted sum and difference of the pathway-specific components. We find that the I1-FFL is able to manufacture more varied redundancy and extrinsic noise responses compared to the C1-FFL. Underlying the differing characteristics of the composite metrics across FFL variants, there exist uniformly behaving pathway-dependent elements. The decomposition framework has been meticulously explored in biologically rational parametric realizations through analytical estimates and stochastic simulations.

Categorizing the Role of Respiration in Cardiovascular and Cerebrovascular Variability Interactions

OBJECTIVE

Respiration disturbs cardiovascular and cerebrovascular controls but its role is not fully elucidated.

METHODS

Respiration can be classified as a confounder if its observation reduces the strength of the causal relationship from source to target. Respiration is a suppressor if the opposite situation holds. We prove that a confounding/suppression (C/S) test can be accomplished by evaluating the sign of net redundancy/synergy balance in the predictability framework based on multivariate autoregressive modelling. In addition, we suggest that, under the hypothesis of Gaussian processes, the C/S test can be given in the transfer entropy decomposition framework as well. Experimental protocols: We applied the C/S test to variability series of respiratory movements, heart period, systolic arterial pressure, mean arterial pressure, and mean cerebral blood flow recorded in 17 pathological individuals (age: 648 yrs; 17 males) before and after induction of propofol-based general anesthesia prior to coronary artery bypass grafting, and in 13 healthy subjects (age: 278 yrs; 5 males) at rest in supine position and during head-up tilt with a table inclination of 60.

RESULTS

Respiration behaved systematically as a confounder for cardiovascular and cerebrovascular controls. In addition, its role was affected by propofol-based general anesthesia but not by a postural stimulus of limited intensity.

CONCLUSION

The C/S test can be fruitfully exploited to categorize the role of respiration over causal variability interactions.

SIGNIFICANCE

The application of the C/S test could favor the comprehension of the role of respiration in cardiovascular and cerebrovascular regulations.

Stroke-related alterations in inter-areal communication

•

We used covariance-based Granger Causality on resting-state fMRI of stroke patients.

•

Stroke determines an overall decrease of homotopic Granger causality (GC)

•

Stroke determines a decrease of GC from and within the lesioned hemisphere.

•

Stroke causes imbalances in GC between the lesioned and the healthy hemisphere.

•

GC anomalies correlate with impaired performance in several behavioral domains.

Beyond causing local ischemia and cell damage at the site of injury, stroke strongly affects long-range anatomical connections, perturbing the functional organization of brain networks. Several studies reported functional connectivity abnormalities parallelling both behavioral deficits and functional recovery across different cognitive domains. FC alterations suggest that long-range communication in the brain is altered after stroke. However, standard FC analyses cannot reveal the directionality and time scale of inter-areal information transfer. We used resting-state fMRI and covariance-based Granger causality analysis to quantify network-level information transfer and its alteration in stroke. Two main large-scale anomalies were observed in stroke patients. First, inter-hemispheric information transfer was significantly decreased with respect to healthy controls. Second, stroke caused inter-hemispheric asymmetries, as information transfer within the affected hemisphere and from the affected to the intact hemisphere was significantly reduced. Both anomalies were more prominent in resting-state networks related to attention and language, and they correlated with impaired performance in several behavioral domains. Overall, our findings support the hypothesis that stroke provokes asymmetries between the affected and spared hemisphere, with different functional consequences depending on which hemisphere is lesioned.

A Controlled Thermoalgesic Stimulation Device for Exploring Novel Pain Perception Biomarkers

OBJECTIVE

To develop a new device for identifying physiological markers of pain perception by reading the brain's electrical activity and hemodynamic interactions while applying thermoalgesic stimulation.

METHODS

We designed a compact prototype that generates well-controlled thermal stimuli using a computer-driven Peltier cell while simultaneously capturing electroencephalography (EEG) and photoplethysmography (PPG) signals. The study was performed on 35 healthy subjects (mean age 30.46 years, SD 4.93 years; 20 males, 15 females). We first determined the heat pain threshold (HPT) for each subject, defined as the maximum temperature that the subject can withstand when the Peltier cell gradually increased the temperature. Next, we defined the painful condition as the one occurring at temperature equal to 90% of the HPT, comparing this to the no-pain state (control) in the absence of thermoalgesic stimulation.

RESULTS

Both the one-dimensional and the two-dimensional spectral entropy (SE) obtained from both the EEG and PPG signals differentiated the condition of pain. In particular, the SE for PPG was significantly reduced in association with pain, while the SE for EEG increased slightly. Moreover, significant discrimination occurred within a specific range of frequencies, 26-30 Hz for EEG and about 5-10 Hz for PPG.

CONCLUSION

Hemodynamics, brain dynamics and their interactions can discriminate thermal pain perception.

SIGNIFICANCE

The possibility of monitoring on-line variations in thermal pain perception using a similar device and algorithms may be of interest to study different pathologies that affect the peripheral nervous system, such as small fiber neuropathies, fibromyalgia or painful diabetic neuropathy.

High-Order Interdependencies in the Aging Brain

Background: Brain interdependencies can be studied from either a structural/anatomical perspective ("structural connectivity") or by considering statistical interdependencies ("functional connectivity" [FC]). Interestingly, while structural connectivity is by definition pairwise (white-matter fibers project from one region to another), FC is not. However, most FC analyses only focus on pairwise statistics and they neglect higher order interactions. A promising tool to study high-order interdependencies is the recently proposed O-Information, which can quantify the intrinsic statistical synergy and the redundancy in groups of three or more interacting variables. Methods: We analyzed functional magnetic resonance imaging (fMRI) data obtained at rest from 164 healthy subjects with ages ranging in 10 to 80 years and used O-Information to investigate how high-order statistical interdependencies are affected by age. Results: Older participants (from 60 to 80 years old) exhibited a higher predominance of redundant dependencies compared with younger participants, an effect that seems to be pervasive as it is evident for all orders of interaction. In addition, while there is strong heterogeneity across brain regions, we found a "redundancy core" constituted by the prefrontal and motor cortices in which redundancy was evident at all the interaction orders studied. Discussion: High-order interdependencies in fMRI data reveal a dominant redundancy in functions such as working memory, executive, and motor functions. Our methodology can be used for a broad range of applications, and the corresponding code is freely available.

Human brain networks in physiological and pathological aging: reproducibility of EEG graph theoretical analysis in cortical connectivity

BACKGROUND

Physiological and pathological brain aging plays a central role in brain networks modulation. The aim of the present paper was to assess the stability of a proposed method for the evaluation of Small World (SW) characteristics for the study of Human Connectome.

METHODS

80 subjects were recruited: 36 young healthy controls, 32 elderly healthy controls, and 12 patients affected by Alzheimer's disease. Electroencephalograms (EEG) were recorded during six separate sessions (480 recordings) at an average inter-session interval of 3.8±0.2 days. Graph theory functions were applied to the undirected and weighted networks obtained by the lagged linear coherence evaluated by exact Low Resolution Electromagnetic Tomography (eLORETA). Were explored the following frequency bands: delta (2-4Hz), theta (4-8Hz), alpha1 (8-10.5Hz), alpha2 (10.5-13Hz), beta1 (13-20Hz), beta2 (20-30Hz) and gamma (30-40Hz).

RESULTS

The proposed method for the evaluation of Small World (SW) characteristics showed good reproducibility and stability. Furthermore, the results showed the pattern Young>Elderly>AD in low frequency delta and theta bands and vice versa in the higher alpha band. Finally, the correlation with age was confirmed in healthy subjects showing that older the age higher the SW values for alpha2.

DISCUSSION

Evidences from the present study confirm the stability of the Small World index and suggest that graph theory can support the analysis of connectivity patterns estimated from EEG. The proposed method for the evaluation of the characteristics of the Small World (SW) has shown good reproducibility and stability and applied to patient data, this technique could provide more information on the pathophysiological processes underlying the age-related brain disconnection, as well as on the administration of rehabilitation treatments at the right time that could allow to avoid unnecessary interventions.

MAXENT3D_PID: An Estimator for the Maximum-Entropy Trivariate Partial Information Decomposition

Partial information decomposition (PID) separates the contributions of sources about a target into unique, redundant, and synergistic components of information. In essence, PID answers the question of “who knows what” of a system of random variables and hence has applications to a wide spectrum of fields ranging from social to biological sciences. The paper presents MaxEnt3D_Pid, an algorithm that computes the PID of three sources, based on a recently-proposed maximum entropy measure, using convex optimization (cone programming). We describe the algorithm and its associated software utilization and report the results of various experiments assessing its accuracy. Moreover, the paper shows that a hierarchy of bivariate and trivariate PID allows obtaining the finer quantities of the trivariate partial information measure.

Multivariate information processing characterizes fitness of a cascaded gene-transcription machinery

We report that a genetic two-step activation cascade processes diverse flavors of information, e.g., synergy, redundancy, and unique information. Our computations measuring reduction in Shannon entropies and reduction in variances produce differently behaving absolute magnitudes of these informational flavors. We find that similarity can be brought in if these terms are evaluated in fractions with respect to corresponding total information. Each of the input signal and final gene-product is found to generate common or redundant information fractions (mostly) to predict each other, whereas they also complement one another to harness synergistic information fraction, predicting the intermediate biochemical species. For an optimally growing signal to maintain fixed steady-state abundance of activated downstream gene-products, the interaction information fractions for this cascade module shift from net-redundancy to information-independence.

Electrophysiological Brain Connectivity: Theory and Implementation

We review the theory and algorithms of electrophysiological brain connectivity analysis. This tutorial is aimed at providing an introduction to brain functional connectivity from electrophysiological signals, including electroencephalography (EEG), magnetoencephalography (MEG), electrocorticography (ECoG), stereoelectroencephalography (SEEG). Various connectivity estimators are discussed, and algorithms introduced. Important issues for estimating and mapping brain functional connectivity with electrophysiology are discussed.

Paced Breathing Increases the Redundancy of Cardiorespiratory Control in Healthy Individuals and Chronic Heart Failure Patients

Synergy and redundancy are concepts that suggest, respectively, adaptability and fault tolerance of systems with complex behavior. This study computes redundancy/synergy in bivariate systems formed by a target X and a driver Y according to the predictive information decomposition approach and partial information decomposition framework based on the minimal mutual information principle. The two approaches assess the redundancy/synergy of past of X and Y in reducing the uncertainty of the current state of X. The methods were applied to evaluate the interactions between heart and respiration in healthy young subjects (n = 19) during controlled breathing at 10, 15 and 20 breaths/minute and in two groups of chronic heart failure patients during paced respiration at 6 (n = 9) and 15 (n = 20) breaths/minutes from spontaneous beat-to-beat fluctuations of heart period and respiratory signal. Both methods suggested that slowing respiratory rate below the spontaneous frequency increases redundancy of cardiorespiratory control in both healthy and pathological groups, thus possibly improving fault tolerance of the cardiorespiratory control. The two methods provide markers complementary to respiratory sinus arrhythmia and the strength of the linear coupling between heart period variability and respiration in describing the physiology of the cardiorespiratory reflex suitable to be exploited in various pathophysiological settings.

Interaction Information Along Lifespan of the Resting Brain Dynamics Reveals a Major Redundant Role of the Default Mode Network

Interaction Information (II) generalizes the univariate Shannon entropy to triplets of variables, allowing the detection of redundant (R) or synergetic (S) interactions in dynamical networks. Here, we calculated II from functional magnetic resonance imaging data and asked whether R or S vary across brain regions and along lifespan. Preserved along lifespan, we found high overlapping between the pattern of high R and the default mode network, whereas high values of S were overlapping with different cognitive domains, such as spatial and temporal memory, emotion processing and motor skills. Moreover, we have found a robust balance between R and S among different age intervals, indicating informational compensatory mechanisms in brain networks.

Structure-function multi-scale connectomics reveals a major role of the fronto-striato-thalamic circuit in brain aging

Physiological aging affects brain structure and function impacting morphology, connectivity, and performance. However, whether some brain connectivity metrics might reflect the age of an individual is still unclear. Here, we collected brain images from healthy participants (N = 155) ranging from 10 to 80 years to build functional (resting state) and structural (tractography) connectivity matrices, both data sets combined to obtain different connectivity features. We then calculated the brain connectome age-an age estimator resulting from a multi-scale methodology applied to the structure-function connectome, and compared it to the chronological age (ChA). Our results were twofold. First, we found that aging widely affects the connectivity of multiple structures, such as anterior cingulate and medial prefrontal cortices, basal ganglia, thalamus, insula, cingulum, hippocampus, parahippocampus, occipital cortex, fusiform, precuneus, and temporal pole. Second, we found that the connectivity between basal ganglia and thalamus to frontal areas, also known as the fronto-striato-thalamic (FST) circuit, makes the major contribution to age estimation. In conclusion, our results highlight the key role played by the FST circuit in the process of healthy aging. Notably, the same methodology can be generally applied to identify the structural-functional connectivity patterns correlating to other biomarkers than ChA.

Pointwise Partial Information Decomposition Using the Specificity and Ambiguity Lattices

What are the distinct ways in which a set of predictor variables can provide information about a target variable? When does a variable provide unique information, when do variables share redundant information, and when do variables combine synergistically to provide complementary information? The redundancy lattice from the partial information decomposition of Williams and Beer provided a promising glimpse at the answer to these questions. However, this structure was constructed using a much criticised measure of redundant information, and despite sustained research, no completely satisfactory replacement measure has been proposed. In this paper, we take a different approach, applying the axiomatic derivation of the redundancy lattice to a single realisation from a set of discrete variables. To overcome the difficulty associated with signed pointwise mutual information, we apply this decomposition separately to the unsigned entropic components of pointwise mutual information which we refer to as the specificity and ambiguity. This yields a separate redundancy lattice for each component. Then based upon an operational interpretation of redundancy, we define measures of redundant specificity and ambiguity enabling us to evaluate the partial information atoms in each lattice. These atoms can be recombined to yield the sought-after multivariate information decomposition. We apply this framework to canonical examples from the literature and discuss the results and the various properties of the decomposition. In particular, the pointwise decomposition using specificity and ambiguity satisfies a chain rule over target variables, which provides new insights into the so-called two-bit-copy example.

The Identity of Information: How Deterministic Dependencies Constrain Information Synergy and Redundancy

Understanding how different information sources together transmit information is crucial in many domains. For example, understanding the neural code requires characterizing how different neurons contribute unique, redundant, or synergistic pieces of information about sensory or behavioral variables. Williams and Beer (2010) proposed a partial information decomposition (PID) that separates the mutual information that a set of sources contains about a set of targets into nonnegative terms interpretable as these pieces. Quantifying redundancy requires assigning an identity to different information pieces, to assess when information is common across sources. Harder et al. (2013) proposed an identity axiom that imposes necessary conditions to quantify qualitatively common information. However, Bertschinger et al. (2012) showed that, in a counterexample with deterministic target-source dependencies, the identity axiom is incompatible with ensuring PID nonnegativity. Here, we study systematically the consequences of information identity criteria that assign identity based on associations between target and source variables resulting from deterministic dependencies. We show how these criteria are related to the identity axiom and to previously proposed redundancy measures, and we characterize how they lead to negative PID terms. This constitutes a further step to more explicitly address the role of information identity in the quantification of redundancy. The implications for studying neural coding are discussed.

Enhanced prefrontal functional-structural networks to support postural control deficits after traumatic brain injury in a pediatric population

Traumatic brain injury (TBI) affects structural connectivity, triggering the reorganization of structural-functional circuits in a manner that remains poorly understood. We focus here on brain network reorganization in relation to postural control deficits after TBI. We enrolled young participants who had suffered moderate to severe TBI, comparing them to young, typically developing control participants. TBI patients (but not controls) recruited prefrontal regions to interact with two separated networks: (1) a subcortical network, including parts of the motor network, basal ganglia, cerebellum, hippocampus, amygdala, posterior cingulate gyrus, and precuneus; and (2) a task-positive network, involving regions of the dorsal attention system, together with dorsolateral and ventrolateral prefrontal regions. We also found that the increased prefrontal connectivity in TBI patients was correlated with some postural control indices, such as the amount of body sway, whereby patients with worse balance increased their connectivity in frontal regions more strongly. The increased prefrontal connectivity found in TBI patients may provide the structural scaffolding for stronger cognitive control of certain behavioral functions, consistent with the observations that various motor tasks are performed less automatically following TBI and that more cognitive control is associated with such actions.

The Power of Connecting Dots: Advanced Techniques to Evaluate Brain Functional Connectivity in Humans

Brain functional connectivity estimation allows us to depict patterns of cerebral activity not understandable otherwise with the standard brain imaging techniques such as functional magnetic resonance imaging (fMRI) as well as electro or magnetoencephalography (hr-EEG, MEG). This special issue of the IEEE Transactions on Biomedical Engineering reports a range of methodological innovations toward the estimation of functional connectivity from brain activity data, with emphasis on neuroelectric and hemodynamic imaging modalities. Functional connectivity methodologies enable "connecting of the dots" derived from brain activity observations over multiple distributed sites, as depicted by such fMRI and hr-EEG/MEG devices.