Lagged and instantaneous dynamical influences related to brain structural connectivity

The structure of anticorrelated networks in the human brain

During the performance of a specific task--or at rest--, the activity of different brain regions shares statistical dependencies that reflect functional connections. While these relationships have been studied intensely for positively correlated networks, considerably less attention has been paid to negatively correlated networks, a. k.a. anticorrelated networks (ACNs). Although the most celebrated of all ACNs is the default mode network (DMN), and has even been extensively studied in health and disease, for systematically all ACNs other than DMN, there is no comprehensive study yet. Here, we have addressed this issue by making use of three neuroimaging data sets: one of N = 192 healthy young adults to fully describe ACN, another of N = 40 subjects to compare ACN between two groups of young and old participants, and another of N = 1,000 subjects from the Human Connectome Project to evaluate the association between ACN and cognitive scores. We first provide a comprehensive description of the anatomical composition of all ACNs, each of which participated in distinct resting-state networks (RSNs). In terms of participation ranking, from highest to the lowest, the major anticorrelated brain areas are the precuneus, the anterior supramarginal gyrus and the central opercular cortex. Next, by evaluating a more detailed structure of ACN, we show it is possible to find significant differences in ACN between specific conditions, in particular, by comparing groups of young and old participants. Our main finding is that of increased anticorrelation for cerebellar interactions in older subjects. Finally, in the voxel-level association study with cognitive scores, we show that ACN has multiple clusters of significance, clusters that are different from those obtained from positive correlated networks, indicating a functional cognitive meaning of ACN. Overall, our results give special relevance to ACN and suggest their use to disentangle unknown alterations in certain conditions, as could occur in early-onset neurodegenerative diseases or in some psychiatric conditions.

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.

Individual-based morphological brain network organization and its association with autistic symptoms in young children with autism spectrum disorder

Individual-based morphological brain networks built from T1-weighted magnetic resonance imaging (MRI) reflect synchronous maturation intensities between anatomical regions at the individual level. Autism spectrum disorder (ASD) is a socio-cognitive and neurodevelopmental disorder with high neuroanatomical heterogeneity, but the specific patterns of morphological networks in ASD remain largely unexplored at the individual level. In this study, individual-based morphological networks were constructed by using high-resolution structural MRI data from 40 young children with ASD (age range: 2-8 years) and 38 age-, gender-, and handedness-matched typically developing children (TDC). Measurements were recorded as threefold. Results showed that compared with TDC, young children with ASD exhibited lower values of small-worldness (i.e., σ) of individual-level morphological brain networks, increased morphological connectivity in cortico-striatum-thalamic-cortical (CSTC) circuitry, and decreased morphological connectivity in the cortico-cortical network. In addition, morphological connectivity abnormalities can predict the severity of social communication deficits in young children with ASD, thus confirming an associational impact at the behavioral level. These findings suggest that the morphological brain network in the autistic developmental brain is inefficient in segregating and distributing information. The results also highlight the crucial role of abnormal morphological connectivity patterns in the socio-cognitive deficits of ASD and support the possible use of the aberrant developmental patterns of morphological brain networks in revealing new clinically-relevant biomarkers for ASD.

Metastable Resting State Brain Dynamics

Metastability refers to the fact that the state of a dynamical system spends a large amount of time in a restricted region of its available phase space before a transition takes place, bringing the system into another state from where it might recur into the previous one. beim Graben and Hutt (2013) suggested to use the recurrence plot (RP) technique introduced by Eckmann et al. (1987) for the segmentation of system's trajectories into metastable states using recurrence grammars. Here, we apply this recurrence structure analysis (RSA) for the first time to resting-state brain dynamics obtained from functional magnetic resonance imaging (fMRI). Brain regions are defined according to the brain hierarchical atlas (BHA) developed by Diez et al. (2015), and as a consequence, regions present high-connectivity in both structure (obtained from diffusion tensor imaging) and function (from the blood-level dependent-oxygenation-BOLD-signal). Remarkably, regions observed by Diez et al. were completely time-invariant. Here, in order to compare this static picture with the metastable systems dynamics obtained from the RSA segmentation, we determine the number of metastable states as a measure of complexity for all subjects and for region numbers varying from 3 to 100. We find RSA convergence toward an optimal segmentation of 40 metastable states for normalized BOLD signals, averaged over BHA modules. Next, we build a bistable dynamics at population level by pooling 30 subjects after Hausdorff clustering. In link with this finding, we reflect on the different modeling frameworks that can allow for such scenarios: heteroclinic dynamics, dynamics with riddled basins of attraction, multiple-timescale dynamics. Finally, we characterize the metastable states both functionally and structurally, using templates for resting state networks (RSNs) and the automated anatomical labeling (AAL) atlas, respectively.

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.

Patient-specific computational modeling of cortical spreading depression via diffusion tensor imaging

Cortical spreading depression, a depolarization wave originating in the visual cortex and traveling towards the frontal lobe, is commonly accepted as a correlate of migraine visual aura. As of today, little is known about the mechanisms that can trigger or stop such phenomenon. However, the complex and highly individual characteristics of the brain cortex suggest that the geometry might have a significant impact in supporting or contrasting the propagation of cortical spreading depression. Accurate patient-specific computational models are fundamental to cope with the high variability in cortical geometries among individuals, but also with the conduction anisotropy induced in a given cortex by the complex neuronal organisation in the grey matter. In this paper, we integrate a distributed model for extracellular potassium concentration with patient-specific diffusivity tensors derived locally from diffusion tensor imaging data.

Group-Level Progressive Alterations in Brain Connectivity Patterns Revealed by Diffusion-Tensor Brain Networks across Severity Stages in Alzheimer's Disease

Alzheimer's disease (AD) is a chronically progressive neurodegenerative disease highly correlated to aging. Whether AD originates by targeting a localized brain area and propagates to the rest of the brain across disease-severity progression is a question with an unknown answer. Here, we aim to provide an answer to this question at the group-level by looking at differences in diffusion-tensor brain networks. In particular, making use of data from Alzheimer's Disease Neuroimaging Initiative (ADNI), four different groups were defined (all of them matched by age, sex and education level): G1 (N1 = 36, healthy control subjects, Control), G2 (N2 = 36, early mild cognitive impairment, EMCI), G3 (N3 = 36, late mild cognitive impairment, LMCI) and G4 (N4 = 36, AD). Diffusion-tensor brain networks were compared across three disease stages: stage I (Control vs. EMCI), stage II (Control vs. LMCI) and stage III (Control vs. AD). The group comparison was performed using the multivariate distance matrix regression analysis, a technique that was born in genomics and was recently proposed to handle brain functional networks, but here applied to diffusion-tensor data. The results were threefold: First, no significant differences were found in stage I. Second, significant differences were found in stage II in the connectivity pattern of a subnetwork strongly associated to memory function (including part of the hippocampus, amygdala, entorhinal cortex, fusiform gyrus, inferior and middle temporal gyrus, parahippocampal gyrus and temporal pole). Third, a widespread disconnection across the entire AD brain was found in stage III, affecting more strongly the same memory subnetwork appearing in stage II, plus the other new subnetworks, including the default mode network, medial visual network, frontoparietal regions and striatum. Our results are consistent with a scenario where progressive alterations of connectivity arise as the disease severity increases and provide the brain areas possibly involved in such a degenerative process. Further studies applying the same strategy to longitudinal data are needed to fully confirm this scenario.

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.

Correcting for Blood Arrival Time in Global Mean Regression Enhances Functional Connectivity Analysis of Resting State fMRI-BOLD Signals

Resting state functional connectivity analysis is a widely used method for mapping intrinsic functional organization of the brain. Global signal regression (GSR) is commonly employed for removing systemic global variance from resting state BOLD-fMRI data; however, recent studies have demonstrated that GSR may introduce spurious negative correlations within and between functional networks, calling into question the meaning of anticorrelations reported between some networks. In the present study, we propose that global signal from resting state fMRI is composed primarily of systemic low frequency oscillations (sLFOs) that propagate with cerebral blood circulation throughout the brain. We introduce a novel systemic noise removal strategy for resting state fMRI data, “dynamic global signal regression” (dGSR), which applies a voxel-specific optimal time delay to the global signal prior to regression from voxel-wise time series. We test our hypothesis on two functional systems that are suggested to be intrinsically organized into anticorrelated networks: the default mode network (DMN) and task positive network (TPN). We evaluate the efficacy of dGSR and compare its performance with the conventional “static” global regression (sGSR) method in terms of (i) explaining systemic variance in the data and (ii) enhancing specificity and sensitivity of functional connectivity measures. dGSR increases the amount of BOLD signal variance being modeled and removed relative to sGSR while reducing spurious negative correlations introduced in reference regions by sGSR, and attenuating inflated positive connectivity measures. We conclude that incorporating time delay information for sLFOs into global noise removal strategies is of crucial importance for optimal noise removal from resting state functional connectivity maps.

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

The causal, directed interactions between brain regions at rest (brain-brain networks) and between resting-state brain activity and autonomic nervous system (ANS) outflow (brain-heart links) have not been completely elucidated. We collected 7 T resting-state functional magnetic resonance imaging (fMRI) data with simultaneous respiration and heartbeat recordings in nine healthy volunteers to investigate (i) the causal interactions between cortical and subcortical brain regions at rest and (ii) the causal interactions between resting-state brain activity and the ANS as quantified through a probabilistic, point-process-based heartbeat model which generates dynamical estimates for sympathetic and parasympathetic activity as well as sympathovagal balance. Given the high amount of information shared between brain-derived signals, we compared the results of traditional bivariate Granger causality (GC) with a globally conditioned approach which evaluated the additional influence of each brain region on the causal target while factoring out effects concomitantly mediated by other brain regions. The bivariate approach resulted in a large number of possibly spurious causal brain-brain links, while, using the globally conditioned approach, we demonstrated the existence of significant selective causal links between cortical/subcortical brain regions and sympathetic and parasympathetic modulation as well as sympathovagal balance. In particular, we demonstrated a causal role of the amygdala, hypothalamus, brainstem and, among others, medial, middle and superior frontal gyri, superior temporal pole, paracentral lobule and cerebellar regions in modulating the so-called central autonomic network (CAN). In summary, we show that, provided proper conditioning is employed to eliminate spurious causalities, ultra-high-field functional imaging coupled with physiological signal acquisition and GC analysis is able to quantify directed brain-brain and brain-heart interactions reflecting central modulation of ANS outflow.

A posteriori model validation for the temporal order of directed functional connectivity maps

A posteriori model validation for the temporal order of neural directed functional connectivity maps is rare. This is striking because models that require sequential independence among residuals are regularly implemented. The aim of the current study was (a) to apply to directed functional connectivity maps of functional magnetic resonance imaging data an a posteriori model validation procedure (i.e., white noise tests of one-step-ahead prediction errors combined with decision criteria for revising the maps based upon Lagrange Multiplier tests), and (b) to demonstrate how the procedure applies to single-subject simulated, single-subject task-related, and multi-subject resting state data. Directed functional connectivity was determined by the unified structural equation model family of approaches in order to map contemporaneous and first order lagged connections among brain regions at the group- and individual-levels while incorporating external input, then white noise tests were run. Findings revealed that the validation procedure successfully detected unmodeled sequential dependencies among residuals and recovered higher order (greater than one) simulated connections, and that the procedure can accommodate task-related input. Findings also revealed that lags greater than one were present in resting state data: With a group-level network that contained only contemporaneous and first order connections, 44% of subjects required second order, individual-level connections in order to obtain maps with white noise residuals. Results have broad methodological relevance (e.g., temporal validation is necessary after directed functional connectivity analyses because the presence of unmodeled higher order sequential dependencies may bias parameter estimates) and substantive implications (e.g., higher order lags may be common in resting state data).