Impact of brain parcellation on prediction performance in models of cognition and demographics.

Brain connectivity analysis begins with the selection of a parcellation scheme that will define brain regions as nodes of a network whose connections will be studied. Brain connectivity has already been used in predictive modelling of cognition, but it remains unclear if the resolution of the parcellation used can systematically impact the predictive model performance. In this work, structural, functional and combined connectivity were each defined with five different parcellation schemes. The resolution and modality of the parcellation schemes were varied. Each connectivity defined with each parcellation was used to predict individual differences in age, education, sex, executive function, self-regulation, language, encoding and sequence processing. It was found that low-resolution functional parcellation consistently performed above chance at producing generalisable models of both demographics and cognition. However, no single parcellation scheme showed a superior predictive performance across all cognitive domains and demographics. In addition, although parcellation schemes impacted the graph theory measures of each connectivity type (structural, functional and combined), these differences did not account for the out-of-sample predictive performance of the models. Taken together, these findings demonstrate that while high-resolution parcellations may be beneficial for modelling specific individual differences, partial voluming of signals produced by the higher resolution of the parcellation likely disrupts model generalisability.

A mean-field model of glutamate and GABA synaptic dynamics for functional MRS.

Advances in functional magnetic resonance spectroscopy (fMRS) have enabled the quantification of activity-dependent changes in neurotransmitter concentrations in vivo. However, the physiological basis of the large changes in GABA and glutamate observed by fMRS (>10%) over short time scales of less than a minute remain unclear as such changes cannot be accounted for by known synthesis or degradation metabolic pathways. Instead, it has been hypothesized that fMRS detects shifts in neurotransmitter concentrations as they cycle from presynaptic vesicles, where they are largely invisible, to extracellular and cytosolic pools, where they are detectable. The present paper uses a computational modelling approach to demonstrate the viability of this hypothesis. A new mean-field model of the neural mechanisms generating the fMRS signal in a cortical voxel is derived. The proposed macroscopic mean-field model is based on a microscopic description of the neurotransmitter dynamics at the level of the synapse. Specifically, GABA and glutamate are assumed to cycle between three metabolic pools: packaged in the vesicles; active in the synaptic cleft; and undergoing recycling and repackaging in the astrocytic or neuronal cytosol. Computational simulations from the model are used to generate predicted changes in GABA and glutamate concentrations in response to different types of stimuli including pain, vision, and electric current stimulation. The predicted changes in the extracellular and cytosolic pools corresponded to those reported in empirical fMRS data. Furthermore, the model predicts a selective control mechanism of the GABA/glutamate relationship, whereby inhibitory stimulation reduces both neurotransmitters, whereas excitatory stimulation increases glutamate and decreases GABA. The proposed model bridges between neural dynamics and fMRS and provides a mechanistic account for the activity-dependent changes in the glutamate and GABA fMRS signals. Lastly, these results indicate that echo-time may be an important timing parameter that can be leveraged to maximise fMRS experimental outcomes.

Using graph theory as a common language to combine neural structure and function in models of healthy cognitive performance.

Graph theory has been used in cognitive neuroscience to understand how organisational properties of structural and functional brain networks relate to cognitive function. Graph theory may bridge the gap in integration of structural and functional connectivity by introducing common measures of network characteristics. However, the explanatory and predictive value of combined structural and functional graph theory have not been investigated in modelling of cognitive performance of healthy adults. In this work, a Principal Component Regression approach with embedded Step-Wise Regression was used to fit multiple regression models of Executive Function, Self-regulation, Language, Encoding and Sequence Processing with a collection of 20 different graph theoretic measures of structural and functional network organisation used as regressors. The predictive ability of graph theory-based models was compared to that of connectivity-based models. The present work shows that using combinations of graph theory metrics to predict cognition in healthy populations does not produce a consistent benefit relative to making predictions based on structural and functional connectivity values directly.

Combination of structural and functional connectivity explains unique variation in specific domains of cognitive function.

The relationship between structural and functional brain networks has been characterised as complex: the two networks mirror each other and show mutual influence but they also diverge in their organisation. This work explored whether a combination of structural and functional connectivity can improve the fit of regression models of cognitive performance. Principal Component Analysis (PCA) was first applied to cognitive data from the Human Connectome Project to identify latent cognitive components: Executive Function, Self-regulation, Language, Encoding and Sequence Processing. A Principal Component Regression approach with embedded Step-Wise Regression (SWR-PCR) was then used to fit regression models of each cognitive domain based on structural (SC), functional (FC) or combined structural-functional (CC) connectivity. Executive Function was best explained by the CC model. Self-regulation was equally well explained by SC and FC. Language was equally well explained by CC and FC models. Encoding and Sequence Processing were best explained by SC. Evaluation of out-of-sample models' skill via cross-validation showed that SC, FC and CC produced generalisable models of Language performance. SC models performed most effectively at predicting Language performance in unseen sample. Executive Function was most effectively predicted by SC models, followed only by CC models. Self-regulation was only effectively predicted by CC models and Sequence Processing was only effectively predicted by FC models. The present study demonstrates that integrating structural and functional connectivity can help explaining cognitive performance, but that the added explanatory value (in-sample) may be domain-specific and can come at the expense of reduced generalisation performance (out-of-sample).

A tutorial and tool for exploring feature similarity gradients with MRI data.

There has been an increasing interest in examining organisational principles of the cerebral cortex (and subcortical regions) using different MRI features such as structural or functional connectivity. Despite the widespread interest, introductory tutorials on the underlying technique targeted for the novice neuroimager are sparse in the literature. Articles that investigate various "neural gradients" (for example based on region studied "cortical gradients," "cerebellar gradients," "hippocampal gradients" etc … or feature of interest "functional gradients," "cytoarchitectural gradients," "myeloarchitectural gradients" etc …) have increased in popularity. Thus, we believe that it is opportune to discuss what is generally meant by "gradient analysis". We introduce basics concepts in graph theory, such as graphs themselves, the degree matrix, and the adjacency matrix. We discuss how one can think about gradients of feature similarity (the similarity between timeseries in fMRI, or streamline in tractography) using graph theory and we extend this to explore such gradients across the whole MRI scale; from the voxel level to the whole brain level. We proceed to introduce a measure for quantifying the level of similarity in regions of interest. We propose the term "the Vogt-Bailey index" for such quantification to pay homage to our history as a brain mapping community. We run through the techniques on sample datasets including a brain MRI as an example of the application of the techniques on real data and we provide several appendices that expand upon details. To maximise intuition, the appendices contain a didactic example describing how one could use these techniques to solve a particularly pernicious problem that one may encounter at a wedding. Accompanying the article is a tool, available in both MATLAB and Python, that enables readers to perform the analysis described in this article on their own data. We refer readers to the graphical abstract as an overview of the analysis pipeline presented in this work.

Evidence for frequency-dependent cortical plasticity in the human brain.

Frequency-dependent plasticity (FDP) describes adaptation at the synapse in response to stimulation at different frequencies. Its consequence on the structure and function of cortical networks is unknown. We tested whether cortical "resonance," favorable stimulation frequencies at which the sensory cortices respond maximally, influenced the impact of FDP on perception, functional topography, and connectivity of the primary somatosensory cortex using psychophysics and functional imaging (fMRI). We costimulated two digits on the hand synchronously at, above, or below the resonance frequency of the somatosensory cortex, and tested subjects' accuracy and speed on tactile localization before and after costimulation. More errors and slower response times followed costimulation at above- or below-resonance, respectively. Response times were faster after at-resonance costimulation. In the fMRI, the cortical representations of the two digits costimulated above-resonance shifted closer, potentially accounting for the poorer performance. Costimulation at-resonance did not shift the digit regions, but increased the functional coupling between them, potentially accounting for the improved response time. To relate these results to synaptic plasticity, we simulated a network of oscillators incorporating Hebbian learning. Two neighboring patches embedded in a cortical sheet, mimicking the two digit regions, were costimulated at different frequencies. Network activation outside the stimulated patches was greatest at above-resonance frequencies, reproducing the spread of digit representations seen with fMRI. Connection strengths within the patches increased following at-resonance costimulation, reproducing the increased fMRI connectivity. We show that FDP extends to the cortical level and is influenced by cortical resonance.

Cortical Resonance Frequencies Emerge from Network Size and Connectivity.

Neural oscillations occur within a wide frequency range with different brain regions exhibiting resonance-like characteristics at specific points in the spectrum. At the microscopic scale, single neurons possess intrinsic oscillatory properties, such that is not yet known whether cortical resonance is consequential to neural oscillations or an emergent property of the networks that interconnect them. Using a network model of loosely-coupled Wilson-Cowan oscillators to simulate a patch of cortical sheet, we demonstrate that the size of the activated network is inversely related to its resonance frequency. Further analysis of the parameter space indicated that the number of excitatory and inhibitory connections, as well as the average transmission delay between units, determined the resonance frequency. The model predicted that if an activated network within the visual cortex increased in size, the resonance frequency of the network would decrease. We tested this prediction experimentally using the steady-state visual evoked potential where we stimulated the visual cortex with different size stimuli at a range of driving frequencies. We demonstrate that the frequency corresponding to peak steady-state response inversely correlated with the size of the network. We conclude that although individual neurons possess resonance properties, oscillatory activity at the macroscopic level is strongly influenced by network interactions, and that the steady-state response can be used to investigate functional networks.

Hearing silences: human auditory processing relies on preactivation of sound-specific brain activity patterns.

The remarkable capabilities displayed by humans in making sense of an overwhelming amount of sensory information cannot be explained easily if perception is viewed as a passive process. Current theoretical and computational models assume that to achieve meaningful and coherent perception, the human brain must anticipate upcoming stimulation. But how are upcoming stimuli predicted in the brain? We unmasked the neural representation of a prediction by omitting the predicted sensory input. Electrophysiological brain signals showed that when a clear prediction can be formulated, the brain activates a template of its response to the predicted stimulus before it arrives to our senses.

Temporal regularity facilitates higher-order sensory predictions in fast auditory sequences.

Does temporal regularity facilitate prediction in audition? To test this, we recorded human event-related potentials to frequent standard tones and infrequent pitch deviant tones, pre-attentively delivered within isochronous and anisochronous (20% onset jitter) rapid sequences. Deviant tones were repeated, either with high or low probability. Standard tone repetition sets a first-order prediction, which is violated by deviant tone onset, leading to a first-order prediction error response (Mismatch Negativity). The response to highly probable deviant repetitions is, however, attenuated relative to less probable repetitions, reflecting the formation of higher-order sensory predictions. Results show that temporal regularity is required for higher-order predictions, but does not modulate first-order prediction error responses. Inverse solution analyses (Variable Resolution Electrical Tomography; VARETA) localized the error response attenuation to posterior regions of the left superior temporal gyrus. In a control experiment with a slower stimulus rate, we found no evidence for higher-order predictions, and again no effect of temporal information on first-order prediction error. We conclude that: (i) temporal regularity facilitates the establishing of higher-order sensory predictions, i.e. 'knowing what next', in fast auditory sequences; (ii) first-order prediction error relies predominantly on stimulus feature mismatch, reflecting the adaptive fit of fast deviance detection processes.

Test-retest reliability of Bayesian estimations of the effects of stimulation, prior information and individual traits on pain perception.

BACKGROUND: There is inter-individual variability in the influence of different components (e.g. nociception and expectations) on pain perception. Identifying the individual effect of these components could serve for patient stratification, but only if these influences are stable in time. METHODS: In this study, 30 healthy participants underwent a cognitive pain paradigm in which they rated pain after viewing a probabilistic cue informing of forthcoming pain intensity and then receiving electrical stimulation. The trial information was then used in a Bayesian probability model to compute the relative weight each participant put on stimulation, cue, cue uncertainty and trait-like bias. The same procedure was repeated 2 weeks later. Relative and absolute test-retest reliability of all measures was assessed. RESULTS: Intraclass correlation results showed good reliability for the effect of the stimulation (0.83), the effect of the cue (0.75) and the trait-like bias (0.75 and 0.75), and a moderate reliability for the effect of the cue uncertainty (0.55). Absolute reliability measures also supported the temporal stability of the results and indicated that a change in parameters corresponding to a difference in pain ratings ranging between 0.47 and 1.45 (depending on the parameters) would be needed to consider differences in outcomes significant. The comparison of these measures with the closest clinical data we possess supports the reliability of our results. CONCLUSIONS: These findings support the hypothesis that inter-individual differences in the weight placed on different pain factors are stable in time and could therefore be a possible target for patient stratification. SIGNIFICANCE: Our results demonstrate the temporal stability of the weight healthy individuals place on the different factors leading to the pain response. These findings give validity to the idea of using Bayesian estimations of the influence of different factors on pain as a way to stratify patients for treatment personalization.