Phasic Dopamine Changes and Hebbian Mechanisms during Probabilistic Reversal Learning in Striatal Circuits: A Computational Study.

Cognitive flexibility is essential to modify our behavior in a non-stationary environment and is often explored by reversal learning tasks. The basal ganglia (BG) dopaminergic system, under a top-down control of the pre-frontal cortex, is known to be involved in flexible action selection through reinforcement learning. However, how adaptive dopamine changes regulate this process and learning mechanisms for training the striatal synapses remain open questions. The current study uses a neurocomputational model of the BG, based on dopamine-dependent direct (Go) and indirect (NoGo) pathways, to investigate reinforcement learning in a probabilistic environment through a task that associates different stimuli to different actions. Here, we investigated: the efficacy of several versions of the Hebb rule, based on covariance between pre- and post-synaptic neurons, as well as the required control in phasic dopamine changes crucial to achieving a proper reversal learning. Furthermore, an original mechanism for modulating the phasic dopamine changes is proposed, assuming that the expected reward probability is coded by the activity of the winner Go neuron before a reward/punishment takes place. Simulations show that this original formulation for an automatic phasic dopamine control allows the achievement of a good flexible reversal even in difficult conditions. The current outcomes may contribute to understanding the mechanisms for active control of dopamine changes during flexible behavior. In perspective, it may be applied in neuropsychiatric or neurological disorders, such as Parkinson's or schizophrenia, in which reinforcement learning is impaired.

An integrative model of Parkinson's disease treatment including levodopa pharmacokinetics, dopamine kinetics, basal ganglia neurotransmission and motor action throughout disease progression.

Levodopa is considered the gold standard treatment of Parkinson's disease. Although very effective in alleviating symptoms at their onset, its chronic use with the progressive neuronal denervation in the basal ganglia leads to a decrease in levodopa's effect duration and to the appearance of motor complications. This evolution challenges the establishment of optimal regimens to manage the symptoms as the disease progresses. Based on up-to-date pathophysiological and pharmacological knowledge, we developed an integrative model for Parkinson's disease to evaluate motor function in response to levodopa treatment as the disease progresses. We combined a pharmacokinetic model of levodopa to a model of dopamine's kinetics and a neurocomputational model of basal ganglia. The parameter values were either measured directly or estimated from human and animal data. The concentrations and behaviors predicted by our model were compared to available information and data. Using this model, we were able to predict levodopa plasma concentration, its related dopamine concentration in the brain and the response performance of a motor task for different stages of disease.

Alpha and theta mechanisms operating in internal-external attention competition.

Attention is the ability to prioritize a set of information at expense of others and can be internally- or externally-oriented. Alpha and theta oscillations have been extensively implicated in attention. However, it is unclear how these oscillations operate when sensory distractors are presented continuously during task-relevant internal processes, in close-to-real-life conditions. Here, EEG signals from healthy participants were obtained at rest and in three attentional conditions, characterized by the execution of a mental math task (internal attention), presentation of pictures on a monitor (external attention), and task execution under the distracting action of picture presentation (internal-external competition). Alpha and theta power were investigated at scalp level and at some cortical regions of interest (ROIs); moreover, functional directed connectivity was estimated via spectral Granger Causality. Results show that frontal midline theta was distinctive of mental task execution and was more prominent during competition compared to internal attention alone, possibly reflecting higher executive control; anterior cingulate cortex appeared as mainly involved and causally connected to distant (temporal/occipital) regions. Alpha power in visual ROIs strongly decreased in external attention alone, while it assumed values close to rest during competition, reflecting reduced visual engagement against distractors; connectivity results suggested that bidirectional alpha influences between frontal and visual regions could contribute to reduce visual interference in internal attention. This study can help to understand how our brain copes with internal-external attention competition, a condition intrinsic in the human sensory-cognitive interplay, and to elucidate the relationships between brain oscillations and attentional functions/dysfunctions in daily tasks.

A non-linear deterministic model of action selection in the basal ganglia to simulate motor fluctuations in Parkinson's disease.

Motor fluctuations and dyskinesias are severe complications of Parkinson's disease (PD), especially evident at its advanced stage, under long-term levodopa therapy. Despite their strong clinical prevalence, the neural origin of these motor symptoms is still a subject of intense debate. In this work, a non-linear deterministic neurocomputational model of the basal ganglia (BG), inspired by biology, is used to provide more insights into possible neural mechanisms at the basis of motor complications in PD. In particular, the model is used to simulate the finger tapping task. The model describes the main neural pathways involved in the BG to select actions [the direct or Go, the indirect or NoGo, and the hyperdirect pathways via the action of the sub-thalamic nucleus (STN)]. A sensitivity analysis is performed on some crucial model parameters (the dopamine level, the strength of the STN mechanism, and the strength of competition among different actions in the motor cortex) at different levels of synapses, reflecting major or minor motor training. Depending on model parameters, results show that the model can reproduce a variety of clinically relevant motor patterns, including normokinesia, bradykinesia, several attempts before movement, freezing, repetition, and also irregular fluctuations. Motor symptoms are, especially, evident at low or high dopamine levels, with excessive strength of the STN and with weak competition among alternative actions. Moreover, these symptoms worsen if the synapses are subject to insufficient learning. The model may help improve the comprehension of motor complications in PD and, ultimately, may contribute to the treatment design.

Nonlinear pharmacodynamics of levodopa through Parkinson's disease progression.

The effect of levodopa in alleviating the symptoms of Parkinson's disease is altered in a highly nonlinear manner as the disease progresses. This can be attributed to different compensation mechanisms taking place in the basal ganglia where the dopaminergic neurons are progressively lost. This alteration in the effect of levodopa complicates the optimization of a drug regimen. The present work aims at investigating the nonlinear dynamics of Parkinson's disease and its therapy through mechanistic mathematical modeling. Using a holistic approach, a pharmacokinetic model of levodopa was combined to a dopamine dynamics and a neurocomputational model of basal ganglia. The influence of neuronal death on these different mechanisms was also integrated. Using this model, we were able to investigate the nonlinear relationships between the levodopa plasma concentration, the dopamine brain concentration, and a response to a motor task. Variations in dopamine concentrations in the brain for different levodopa doses were also studied. Finally, we investigated the narrowing of a levodopa therapeutic index with the progression of the disease as a result of these nonlinearities. In conclusion, various consequences of nonlinear dynamics in Parkinson's disease treatment were studied by developing an integrative model. This model paves the way toward individualization of a dosing regimen. Using sensor based information, the parameters of the model could be fitted to individual data to propose optimal individual regimens.

Aberrant learning in Parkinson's disease: A neurocomputational study on bradykinesia.

Parkinson's disease (PD) is a neurodegenerative disorder characterized by a progressive decline in motor functions, such as bradykinesia, caused by the pathological denervation of nigrostriatal dopaminergic neurons within the basal ganglia (BG). It is acknowledged that dopamine (DA) directly affects the modulatory role of BG towards the cortex. However, a growing body of literature is suggesting that DA-induced aberrant synaptic plasticity could play a role in the core symptoms of PD, thus recalling for a "reconceptualization" of the pathophysiology. The aim of this work was to investigate DA-driven aberrant learning as a concurrent cause of bradykinesia, using a comprehensive, biologically inspired neurocomputational model of action selection in the BG. The model includes the three main pathways operating in the BG circuitry, that is the direct, indirect and hyperdirect pathways, and use a two-term Hebb rule to train synapses in the striatum, based on previous history of rewards and punishments. Levodopa pharmacodynamics is also incorporated. Through model simulations of the Alternate Finger Tapping motor task, we assessed the role of aberrant learning on bradykinesia. The results show that training under drug medication (levodopa) provides not only immediate but also delayed benefit lasting in time. Conversely, if performed in conditions of vanishing levodopa efficacy, training may result in dysfunctional corticostriatal synaptic plasticity, further worsening motor performances in PD subjects. This suggests that bradykinesia may result from the concurrent effects of low DA levels and dysfunctional plasticity and that training can be exploited in medicated subjects to improve levodopa treatment.

Sodium Prescription in the Prevention of Intradialytic Hypotension: New Insights into an Old Concept.

BACKGROUND: Sodium prescription in patients with intradialytic hypotension remains a challenge for the attending nephrologist, as it increases dialysate conductivity in hypotension-prone patients, thereby adding to dietary sodium levels. METHODS: New sodium prescription strategies are now available, including the use of a mathematical model to compute the sodium mass to be removed during dialysis as a physiological controller. RESULTS: This review describes the sodium load of patients with end-stage renal disease on chronic hemodialysis (HD) and discusses 2 strategies to remove excess sodium in patients prone to intradialytic hypotension, namely, Profiled HD and the hemodiafiltration Aequilibrium System. CONCLUSION: The Profiled HD and Aequilibrium System trial both proved effective in counteracting intradialytic hypotension.

A biologically inspired neurocomputational model for audiovisual integration and causal inference.

Recently, experimental and theoretical research has focused on the brain's abilities to extract information from a noisy sensory environment and how cross-modal inputs are processed to solve the causal inference problem to provide the best estimate of external events. Despite the empirical evidence suggesting that the nervous system uses a statistically optimal and probabilistic approach in addressing these problems, little is known about the brain's architecture needed to implement these computations. The aim of this work was to realize a mathematical model, based on physiologically plausible hypotheses, to analyze the neural mechanisms underlying multisensory perception and causal inference. The model consists of three layers topologically organized: two encode auditory and visual stimuli, separately, and are reciprocally connected via excitatory synapses and send excitatory connections to the third downstream layer. This synaptic organization realizes two mechanisms of cross-modal interactions: the first is responsible for the sensory representation of the external stimuli, while the second solves the causal inference problem. We tested the network by comparing its results to behavioral data reported in the literature. Among others, the network can account for the ventriloquism illusion, the pattern of sensory bias and the percept of unity as a function of the spatial auditory-visual distance, and the dependence of the auditory error on the causal inference. Finally, simulations results are consistent with probability matching as the perceptual strategy used in auditory-visual spatial localization tasks, agreeing with the behavioral data. The model makes untested predictions that can be investigated in future behavioral experiments.

Multisensory Bayesian Inference Depends on Synapse Maturation during Training: Theoretical Analysis and Neural Modeling Implementation.

Recent theoretical and experimental studies suggest that in multisensory conditions, the brain performs a near-optimal Bayesian estimate of external events, giving more weight to the more reliable stimuli. However, the neural mechanisms responsible for this behavior, and its progressive maturation in a multisensory environment, are still insufficiently understood. The aim of this letter is to analyze this problem with a neural network model of audiovisual integration, based on probabilistic population coding-the idea that a population of neurons can encode probability functions to perform Bayesian inference. The model consists of two chains of unisensory neurons (auditory and visual) topologically organized. They receive the corresponding input through a plastic receptive field and reciprocally exchange plastic cross-modal synapses, which encode the spatial co-occurrence of visual-auditory inputs. A third chain of multisensory neurons performs a simple sum of auditory and visual excitations. The work includes a theoretical part and a computer simulation study. We show how a simple rule for synapse learning (consisting of Hebbian reinforcement and a decay term) can be used during training to shrink the receptive fields and encode the unisensory likelihood functions. Hence, after training, each unisensory area realizes a maximum likelihood estimate of stimulus position (auditory or visual). In cross-modal conditions, the same learning rule can encode information on prior probability into the cross-modal synapses. Computer simulations confirm the theoretical results and show that the proposed network can realize a maximum likelihood estimate of auditory (or visual) positions in unimodal conditions and a Bayesian estimate, with moderate deviations from optimality, in cross-modal conditions. Furthermore, the model explains the ventriloquism illusion and, looking at the activity in the multimodal neurons, explains the automatic reweighting of auditory and visual inputs on a trial-by-trial basis, according to the reliability of the individual cues.

An integrated mathematical model of the human cardiopulmonary system: model development.

Several cardiovascular and pulmonary models have been proposed in the last few decades. However, very few have addressed the interactions between these two systems. Our group has developed an integrated cardiopulmonary model (CP Model) that mathematically describes the interactions between the cardiovascular and respiratory systems, along with their main short-term control mechanisms. The model has been compared with human and animal data taken from published literature. Due to the volume of the work, the paper is divided in two parts. The present paper is on model development and normophysiology, whereas the second is on the model's validation on hypoxic and hypercapnic conditions. The CP Model incorporates cardiovascular circulation, respiratory mechanics, tissue and alveolar gas exchange, as well as short-term neural control mechanisms acting on both the cardiovascular and the respiratory functions. The model is able to simulate physiological variables typically observed in adult humans under normal and pathological conditions and to explain the underlying mechanisms and dynamics.

An integrated mathematical model of the human cardiopulmonary system: model validation under hypercapnia and hypoxia.

A novel integrated physiological model of the interactions between the cardiovascular and respiratory systems has been in development for the past few years. The model has hundreds of parameters and variables representing the physical and physiological properties of the human cardiopulmonary system. It can simulate many dynamic states and scenarios. The description of the model and the results in normal resting conditions were presented in a companion paper (Albanese A, Cheng L, Ursino M, Chbat NW.Am J Physiol Heart Circ Physiol 310: 2016; doi:10.1152/ajpheart.00230.2014), where model predictions were compared against average population data from literature. However, it is also essential to test the model in abnormal or pathological conditions to prove its consistency. Hence, in this paper, we concentrate on testing the cardiopulmonary model under hypercapnic and hypoxic conditions, by comparing model's outputs to population-averaged cardiorespiratory data reported in the literature. The utility of this comprehensive model is demonstrated by testing the internal consistency of the simulated responses of a significant number of cardiovascular variables (heart rate, arterial pressure, and cardiac output) and respiratory variables (tidal volume, respiratory rate, minute ventilation, alveolar O2 and CO2 partial pressures) over a wide range of perturbations and conditions; namely, hypercapnia at 3-7% CO2 levels and hypoxia at 7-9% O2 levels with controlled CO2(isocapnic hypoxia) and without controlled CO2(hypocapnic hypoxia). Finally, a sensitivity analysis is performed to analyze the role of the main cardiorespiratory control mechanisms triggered by hypercapnia and hypoxia.

A neural mass model of place cell activity: theta phase precession, replay and imagination of never experienced paths.

Recent results on hippocampal place cells show that the replay of behavioral sequences does not simply reflect previously experienced trajectories, but may also occur in the reverse direction, or may even include never experienced paths. In order to elucidate the possible mechanisms at the basis of this phenomenon, we have developed a model of sequence learning. The present model consists of two layers of place cell units. Long-range connections among units implement heteroassociation between the two layers, trained with a temporal Hebb rule. The network was trained assuming that a virtual rat moves within a virtual maze. This training leads to the formation of bidirectional synapses between the two layers, i.e. synapses connecting a neuron both with its previous and subsequent element in the path. Subsequently, two distinct conditions were simulated with the trained network. During an exploratory phase, characterized by a similar consideration to the external environment and to the internal representation, the model simulates the occurrence of theta precession in the forward path and the temporal compression. During an imagination phase, when there is no consideration to the external location, the model produces trains of gamma oscillations, without the presence of a theta rhythm, and simulates the occurrence of both direct and reverse replay, and the imagination of never experienced paths. The new paths are built by combining bunches of previous trajectories. The main mechanisms at the basis of this behavior are explained in detail, and lines for future improvements (e.g., to simulate preplay) are discussed.

A neurocomputational analysis of the sound-induced flash illusion.

Perception of the external world is based on the integration of inputs from different sensory modalities. Recent experimental findings suggest that this phenomenon is present in lower-level cortical areas at early processing stages. The mechanisms underlying these early processes and the organization of the underlying circuitries are still a matter of debate. Here, we investigate audiovisual interactions by means of a simple neural network consisting of two layers of visual and auditory neurons. We suggest that the spatial and temporal aspects of audio-visual illusions can be explained within this simple framework, based on two main assumptions: auditory and visual neurons communicate via excitatory synapses; and spatio-temporal receptive fields are different in the two modalities, auditory processing exhibiting a higher temporal resolution, while visual processing a higher spatial acuity. With these assumptions, the model is able: i) to simulate the sound-induced flash fission illusion; ii) to reproduce psychometric curves assuming a random variability in some parameters; iii) to account for other audio-visual illusions, such as the sound-induced flash fusion and the ventriloquism illusions; and iv) to predict that visual and auditory stimuli are combined optimally in multisensory integration. In sum, the proposed model provides a unifying summary of spatio-temporal audio-visual interactions, being able to both account for a wide set of empirical findings, and be a framework for future experiments. In perspective, it may be used to understand the neural basis of Bayesian audio-visual inference.

Modeling the contribution of theta-gamma coupling to sequential memory, imagination, and dreaming.

Gamma oscillations nested in a theta rhythm are observed in the hippocampus, where are assumed to play a role in sequential episodic memory, i.e., memorization and retrieval of events that unfold in time. In this work, we present an original neurocomputational model based on neural masses, which simulates the encoding of sequences of events in the hippocampus and subsequent retrieval by exploiting the theta-gamma code. The model is based on a three-layer structure in which individual Units oscillate with a gamma rhythm and code for individual features of an episode. The first layer (working memory in the prefrontal cortex) maintains a cue in memory until a new signal is presented. The second layer (CA3 cells) implements an auto-associative memory, exploiting excitatory and inhibitory plastic synapses to recover an entire episode from a single feature. Units in this layer are disinhibited by a theta rhythm from an external source (septum or Papez circuit). The third layer (CA1 cells) implements a hetero-associative net with the previous layer, able to recover a sequence of episodes from the first one. During an encoding phase, simulating high-acetylcholine levels, the network is trained with Hebbian (synchronizing) and anti-Hebbian (desynchronizing) rules. During retrieval (low-acetylcholine), the network can correctly recover sequences from an initial cue using gamma oscillations nested inside the theta rhythm. Moreover, in high noise, the network isolated from the environment simulates a mind-wandering condition, randomly replicating previous sequences. Interestingly, in a state simulating sleep, with increased noise and reduced synapses, the network can "dream" by creatively combining sequences, exploiting features shared by different episodes. Finally, an irrational behavior (erroneous superimposition of features in various episodes, like "delusion") occurs after pathological-like reduction in fast inhibitory synapses. The model can represent a straightforward and innovative tool to help mechanistically understand the theta-gamma code in different mental states.

Convolutional neural networks reveal properties of reach-to-grasp encoding in posterior parietal cortex.

Deep neural networks (DNNs) are widely adopted to decode motor states from both non-invasively and invasively recorded neural signals, e.g., for realizing brain-computer interfaces. However, the neurophysiological interpretation of how DNNs make the decision based on the input neural activity is limitedly addressed, especially when applied to invasively recorded data. This reduces decoder reliability and transparency, and prevents the exploitation of decoders to better comprehend motor neural encoding. Here, we adopted an explainable artificial intelligence approach - based on a convolutional neural network and an explanation technique - to reveal spatial and temporal neural properties of reach-to-grasping from single-neuron recordings of the posterior parietal area V6A. The network was able to accurately decode 5 different grip types, and the explanation technique automatically identified the cells and temporal samples that most influenced the network prediction. Grip encoding in V6A neurons already started at movement preparation, peaking during movement execution. A difference was found within V6A: dorsal V6A neurons progressively encoded more for increasingly advanced grips, while ventral V6A neurons for increasingly rudimentary grips, with both subareas following a linear trend between the amount of grip encoding and the level of grip skills. By revealing the elements of the neural activity most relevant for each grip with no a priori assumptions, our approach supports and advances current knowledge about reach-to-grasp encoding in V6A, and it may represent a general tool able to investigate neural correlates of motor or cognitive tasks (e.g., attention and memory tasks) from single-neuron recordings.

Graphene oxide electrodes enable electrical stimulation of distinct calcium signalling in brain astrocytes.

Astrocytes are responsible for maintaining homoeostasis and cognitive functions through calcium signalling, a process that is altered in brain diseases. Current bioelectronic tools are designed to study neurons and are not suitable for controlling calcium signals in astrocytes. Here, we show that electrical stimulation of astrocytes using electrodes coated with graphene oxide and reduced graphene oxide induces respectively a slow response to calcium, mediated by external calcium influx, and a sharp one, exclusively due to calcium release from intracellular stores. Our results suggest that the different conductivities of the substrate influence the electric field at the cell-electrolyte or cell-material interfaces, favouring different signalling events in vitro and ex vivo. Patch-clamp, voltage-sensitive dye and calcium imaging data support the proposed model. In summary, we provide evidence of a simple tool to selectively control distinct calcium signals in brain astrocytes for straightforward investigations in neuroscience and bioelectronic medicine.

The formation of categories and the representation of feature saliency: analysis with a computational model trained with an Hebbian paradigm.

An important issue in semantic memory models is the formation of categories and taxonomies, and the different role played by shared vs. distinctive and salient vs. marginal features. Aim of this work is to extend our previous model to critically discuss the mechanisms leading to the formation of categories, and to investigate how feature saliency can be learned from past experience. The model assumes that an object is represented as a collection of features, which belong to different cortical areas and are topologically organized. Excitatory synapses among features are created on the basis of past experience of object presentation, with a Hebbian paradigm, including the use of potentiation and depression of synapses, and thresholding for the presynaptic and postsynaptic. The model was trained using simple schematic objects as input (i.e., vector of features) having some shared features (so as to realize a simple category) and some distinctive features with different frequency. Three different taxonomies of objects were separately trained and tested, which differ as to the number of correlated features and the structure of categories. Results show that categories can be formed from past experience, using Hebbian rules with a different threshold for postsynaptic and presynaptic activity. Furthermore, features have a different saliency, as a consequence of their different frequency during training. The trained network is able to solve simple object recognition tasks, by maintaining a distinction between categories and individual members in the category, and providing a different role for salient features vs. not-salient features. In particular, not-salient features are not evoked in memory when thinking about the object, but they facilitate the reconstruction of objects when provided as input to the model. The results can provide indications on which neural mechanisms can be exploited to form robust categories among objects and on which mechanisms could be implemented in artificial connectionist systems to extract concepts and categories from a continuous stream of input objects (each represented as a vector of features).

A model of working memory for encoding multiple items and ordered sequences exploiting the theta-gamma code.

Recent experimental evidence suggests that oscillatory activity plays a pivotal role in the maintenance of information in working memory, both in rodents and humans. In particular, cross-frequency coupling between theta and gamma oscillations has been suggested as a core mechanism for multi-item memory. The aim of this work is to present an original neural network model, based on oscillating neural masses, to investigate mechanisms at the basis of working memory in different conditions. We show that this model, with different synapse values, can be used to address different problems, such as the reconstruction of an item from partial information, the maintenance of multiple items simultaneously in memory, without any sequential order, and the reconstruction of an ordered sequence starting from an initial cue. The model consists of four interconnected layers; synapses are trained using Hebbian and anti-Hebbian mechanisms, in order to synchronize features in the same items, and desynchronize features in different items. Simulations show that the trained network is able to desynchronize up to nine items without a fixed order using the gamma rhythm. Moreover, the network can replicate a sequence of items using a gamma rhythm nested inside a theta rhythm. The reduction in some parameters, mainly concerning the strength of GABAergic synapses, induce memory alterations which mimic neurological deficits. Finally, the network, isolated from the external environment ("imagination phase") and stimulated with high uniform noise, can randomly recover sequences previously learned, and link them together by exploiting the similarity among items.

Motor decoding from the posterior parietal cortex using deep neural networks.

Objective.Motor decoding is crucial to translate the neural activity for brain-computer interfaces (BCIs) and provides information on how motor states are encoded in the brain. Deep neural networks (DNNs) are emerging as promising neural decoders. Nevertheless, it is still unclear how different DNNs perform in different motor decoding problems and scenarios, and which network could be a good candidate for invasive BCIs.Approach.Fully-connected, convolutional, and recurrent neural networks (FCNNs, CNNs, RNNs) were designed and applied to decode motor states from neurons recorded from V6A area in the posterior parietal cortex (PPC) of macaques. Three motor tasks were considered, involving reaching and reach-to-grasping (the latter under two illumination conditions). DNNs decoded nine reaching endpoints in 3D space or five grip types using a sliding window approach within the trial course. To evaluate decoders simulating a broad variety of scenarios, the performance was also analyzed while artificially reducing the number of recorded neurons and trials, and while performing transfer learning from one task to another. Finally, the accuracy time course was used to analyze V6A motor encoding.Main results.DNNs outperformed a classic Naïve Bayes classifier, and CNNs additionally outperformed XGBoost and Support Vector Machine classifiers across the motor decoding problems. CNNs resulted the top-performing DNNs when using less neurons and trials, and task-to-task transfer learning improved performance especially in the low data regime. Lastly, V6A neurons encoded reaching and reach-to-grasping properties even from action planning, with the encoding of grip properties occurring later, closer to movement execution, and appearing weaker in darkness.Significance.Results suggest that CNNs are effective candidates to realize neural decoders for invasive BCIs in humans from PPC recordings also reducing BCI calibration times (transfer learning), and that a CNN-based data-driven analysis may provide insights about the encoding properties and the functional roles of brain regions.