Neural Intrinsic Timescales in the Macaque Dorsal Premotor Cortex Predict the Strength of Spatial Response Coding

Response inhibition in premotor cortex corresponds to a complex reshuffle of the mesoscopic information network

Recent studies have explored functional and effective neural networks in animal models; however, the dynamics of information propagation among functional modules under cognitive control remain largely unknown. Here, we addressed the issue using transfer entropy and graph theory methods on mesoscopic neural activities recorded in the dorsal premotor cortex of rhesus monkeys. We focused our study on the decision time of a Stop-signal task, looking for patterns in the network configuration that could influence motor plan maturation when the Stop signal is provided. When comparing trials with successful inhibition to those with generated movement, the nodes of the network resulted organized into four clusters, hierarchically arranged, and distinctly involved in information transfer. Interestingly, the hierarchies and the strength of information transmission between clusters varied throughout the task, distinguishing between generated movements and canceled ones and corresponding to measurable levels of network complexity. Our results suggest a putative mechanism for motor inhibition in premotor cortex: a topological reshuffle of the information exchanged among ensembles of neurons.

The laminar organization of cell types in macaque cortex and its relationship to neuronal oscillations

The canonical microcircuit (CMC) has been hypothesized to be the fundamental unit of information processing in cortex. Each CMC unit is thought to be an interconnected column of neurons with specific connections between excitatory and inhibitory neurons across layers. Recently, we identified a conserved spectrolaminar motif of oscillatory activity across the primate cortex that may be the physiological consequence of the CMC. The spectrolaminar motif consists of local field potential (LFP) gamma-band power (40-150 Hz) peaking in superficial layers 2 and 3 and alpha/beta-band power (8-30 Hz) peaking in deep layers 5 and 6. Here, we investigate whether specific conserved cell types may produce the spectrolaminar motif. We collected laminar histological and electrophysiological data in 11 distinct cortical areas spanning the visual hierarchy: V1, V2, V3, V4, TEO, MT, MST, LIP, 8A/FEF, PMD, and LPFC (area 46), and anatomical data in DP and 7A. We stained representative slices for the three main inhibitory subtypes, Parvalbumin (PV), Calbindin (CB), and Calretinin (CR) positive neurons, as well as pyramidal cells marked with Neurogranin (NRGN). We found a conserved laminar structure of PV, CB, CR, and pyramidal cells. We also found a consistent relationship between the laminar distribution of inhibitory subtypes with power in the local field potential. PV interneuron density positively correlated with gamma (40-150 Hz) power. CR and CB density negatively correlated with alpha (8-12 Hz) and beta (13-30 Hz) oscillations. The conserved, layer-specific pattern of inhibition and excitation across layers is therefore likely the anatomical substrate of the spectrolaminar motif.

Significance Statement

Neuronal oscillations emerge as an interplay between excitatory and inhibitory neurons and underlie cognitive functions and conscious states. These oscillations have distinct expression patterns across cortical layers. Does cellular anatomy enable these oscillations to emerge in specific cortical layers? We present a comprehensive analysis of the laminar distribution of the three main inhibitory cell types in primate cortex (Parvalbumin, Calbindin, and Calretinin positive) and excitatory pyramidal cells. We found a canonical relationship between the laminar anatomy and electrophysiology in 11 distinct primate areas spanning from primary visual to prefrontal cortex. The laminar anatomy explained the expression patterns of neuronal oscillations in different frequencies. Our work provides insight into the cortex-wide cellular mechanisms that generate neuronal oscillations in primates.

Neural timescales reflect behavioral demands in freely moving rhesus macaques

Previous work demonstrated a highly reproducible cortical hierarchy of neural timescales at rest, with sensory areas displaying fast, and higher-order association areas displaying slower timescales. The question arises how such stable hierarchies give rise to adaptive behavior that requires flexible adjustment of temporal coding and integration demands. Potentially, this lack of variability in the hierarchical organization of neural timescales could reflect the structure of the laboratory contexts. We posit that unconstrained paradigms are ideal to test whether the dynamics of neural timescales reflect behavioral demands. Here we measured timescales of local field potential activity while male rhesus macaques foraged in an open space. We found a hierarchy of neural timescales that differs from previous work. Importantly, although the magnitude of neural timescales expanded with task engagement, the brain areas' relative position in the hierarchy was stable. Next, we demonstrated that the change in neural timescales is dynamic and contains functionally-relevant information, differentiating between similar events in terms of motor demands and associated reward. Finally, we demonstrated that brain areas are differentially affected by these behavioral demands. These results demonstrate that while the space of neural timescales is anatomically constrained, the observed hierarchical organization and magnitude is dependent on behavioral demands.

MRI Assessment of Intrinsic Neural Timescale and Gray Matter Volume in Parkinson's Disease

BACKGROUND

Numerous studies have indicated altered temporal features of the brain function in Parkinson's disease (PD), and the autocorrelation magnitude of intrinsic neural signals, called intrinsic neural timescales, were often applied to estimate how long neural information stored in local brain areas. However, it is unclear whether PD patients at different disease stages exhibit abnormal timescales accompanied with abnormal gray matter volume (GMV).

PURPOSE

To assess the intrinsic timescale and GMV in PD.

STUDY TYPE

Prospective.

POPULATION

74 idiopathic PD patients (44 early stage (PD-ES) and 30 late stage (PD-LS), as determined by the Hoehn and Yahr (HY) severity classification scale), and 73 healthy controls (HC).

FIELD STRENGTH/SEQUENCE

3.0 T MRI scanner; magnetization prepared rapid acquisition gradient echo and echo planar imaging sequences.

ASSESSMENT

The timescales were estimated by using the autocorrelation magnitude of neural signals. Voxel-based morphometry was performed to calculate GMV in the whole brain. Severity of motor symptoms and cognitive impairments were assessed using the unified PD rating scale, the HY scale, the Montreal cognitive assessment, and the mini-mental state examination.

STATISTICAL TEST

Analysis of variance; two-sample t-test; Spearman rank correlation analysis; Mann-Whitney U test; Kruskal-Wallis' H test. A P value <0.05 was considered statistically significant.

RESULTS

The PD group had significantly abnormal intrinsic timescales in the sensorimotor, visual, and cognitive-related areas, which correlated with the symptom severity (ρ = -0.265, P = 0.022) and GMV (ρ = 0.254, P = 0.029). Compared to the HC group, the PD-ES group had significantly longer timescales in anterior cortical regions, whereas the PD-LS group had significantly shorter timescales in posterior cortical regions.

CONCLUSION

This study suggested that PD patients have abnormal timescales in multisystem and distinct patterns of timescales and GMV in cerebral cortex at different disease stages. This may provide new insights for the neural substrate of PD.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 1.

Cortical timescales and the modular organization of structural and functional brain networks

Recent years have seen growing interest in characterizing the properties of regional brain dynamics and their relationship to other features of brain structure and function. In particular, multiple studies have observed regional differences in the "timescale" over which activity fluctuates during periods of quiet rest. In the cerebral cortex, these timescales have been associated with both local circuit properties as well as patterns of inter-regional connectivity, including the extent to which each region exhibits widespread connectivity to other brain areas. In the current study, we build on prior observations of an association between connectivity and dynamics in the cerebral cortex by investigating the relationship between BOLD fMRI timescales and the modular organization of structural and functional brain networks. We characterize network community structure across multiple scales and find that longer timescales are associated with greater within-community functional connectivity and diverse structural connectivity. We also replicate prior observations of a positive correlation between timescales and structural connectivity degree. Finally, we find evidence for preferential functional connectivity between cortical areas with similar timescales. We replicate these findings in an independent dataset. These results contribute to our understanding of functional brain organization and structure-function relationships in the human brain, and support the notion that regional differences in cortical dynamics may in part reflect the topological role of each region within macroscale brain networks.

Cortical timescales and the modular organization of structural and functional brain networks

Recent years have seen growing interest in characterizing the properties of regional brain dynamics and their relationship to other features of brain structure and function. In particular, multiple studies have observed regional differences in the "timescale" over which activity fluctuates during periods of quiet rest. In the cerebral cortex, these timescales have been associated with both local circuit properties as well as patterns of inter-regional connectivity, including the extent to which each region exhibits widespread connectivity to other brain areas. In the current study, we build on prior observations of an association between connectivity and dynamics in the cerebral cortex by investigating the relationship between BOLD fMRI timescales and the modular organization of structural and functional brain networks. We characterize network community structure across multiple scales and find that longer timescales are associated with greater within-community functional connectivity and diverse structural connectivity. We also replicate prior observations of a positive correlation between timescales and structural connectivity degree. Finally, we find evidence for preferential functional connectivity between cortical areas with similar timescales. We replicate these findings in an independent dataset. These results contribute to our understanding of functional brain organization and structure-function relationships in the human brain, and support the notion that regional differences in cortical dynamics may in part reflect the topological role of each region within macroscale brain networks.

Neural timescales reflect behavioral demands in freely moving rhesus macaques

Previous work has demonstrated remarkably reproducible and consistent hierarchies of neural timescales across cortical areas at rest. The question arises how such stable hierarchies give rise to adaptive behavior that requires flexible adjustment of temporal coding and integration demands. Potentially, this previously found lack of variability in the hierarchical organization of neural timescales could be a reflection of the structure of the laboratory contexts in which they were measured. Indeed, computational work demonstrates the existence of multiple temporal hierarchies within the same anatomical network when the input structure is altered. We posit that unconstrained behavioral environments where relatively little temporal demands are imposed from the experimenter are an ideal test bed to address the question of whether the hierarchical organization and the magnitude of neural timescales reflect ongoing behavioral demands. To tackle this question, we measured timescales of local field potential activity while rhesus macaques were foraging freely in a large open space. We find a hierarchy of neural timescales that is unique to this foraging environment. Importantly, although the magnitude of neural timescales generally expanded with task engagement, the brain areas' relative position in the hierarchy was stable across the recording sessions. Notably, the magnitude of neural timescales monotonically expanded with task engagement across a relatively long temporal scale spanning the duration of the recording session. Over shorter temporal scales, the magnitude of neural timescales changed dynamically around foraging events. Moreover, the change in the magnitude of neural timescales contained functionally relevant information, differentiating between seemingly similar events in terms of motor demands and associated reward. That is, the patterns of change were associated with the cognitive and behavioral meaning of these events. Finally, we demonstrated that brain areas were differentially affected by these behavioral demands - i.e., the expansion of neural timescales was not the same across all areas. Together, these results demonstrate that the observed hierarchy of neural timescales is context-dependent and that changes in the magnitude of neural timescales are closely related to overall task engagement and behavioral demands.

Ongoing Brain Activity and Its Role in Cognition: Dual versus Baseline Models

What is the role of the brain's ongoing activity for cognition? The predominant perspectives associate ongoing brain activity with resting state, the default-mode network (DMN), and internally oriented mentation. This triad is often contrasted with task states, non-DMN brain networks, and externally oriented mentation, together comprising a "dual model" of brain and cognition. In opposition to this duality, however, we propose that ongoing brain activity serves as a neuronal baseline; this builds upon Raichle's original search for the default mode of brain function that extended beyond the canonical default-mode brain regions. That entails what we refer to as the "baseline model." Akin to an internal biological clock for the rest of the organism, the ongoing brain activity may serve as an internal point of reference or standard by providing a shared neural code for the brain's rest as well as task states, including their associated cognition. Such shared neural code is manifest in the spatiotemporal organization of the brain's ongoing activity, including its global signal topography and dynamics like intrinsic neural timescales. We conclude that recent empirical evidence supports a baseline model over the dual model; the ongoing activity provides a global shared neural code that allows integrating the brain's rest and task states, its DMN and non-DMN, and internally and externally oriented cognition.

Intrinsic timescales as an organizational principle of neural processing across the whole rhesus macaque brain

Hierarchical temporal dynamics are a fundamental computational property of the brain; however, there are no whole-brain, noninvasive investigations into timescales of neural processing in animal models. To that end, we used the spatial resolution and sensitivity of ultrahigh field fMRI performed at 10.5 Tesla to probe timescales across the whole macaque brain. We uncovered within-species consistency between timescales estimated from fMRI and electrophysiology. Crucially, we extended existing electrophysiological hierarchies to whole brain topographies. Our results validate the complementary use of hemodynamic and electrophysiological intrinsic timescales, establishing a basis for future translational work. Further, with these results in hand, we were able to show that one facet of the high-dimensional functional connectivity topography of any region in the brain is closely related to hierarchical temporal dynamics. We demonstrated that intrinsic timescales are organized along spatial gradients that closely match functional connectivity gradient topographies across the whole brain. We conclude that intrinsic timescales are a unifying organizational principle of neural processing across the whole brain.

Intrinsic neural timescales: temporal integration and segregation

We are continuously bombarded by external inputs of various timescales from the environment. How does the brain process this multitude of timescales? Recent resting state studies show a hierarchy of intrinsic neural timescales (INT) with a shorter duration in unimodal regions (e.g., visual cortex and auditory cortex) and a longer duration in transmodal regions (e.g., default mode network). This unimodal-transmodal hierarchy is present across acquisition modalities [electroencephalogram (EEG)/magnetoencephalogram (MEG) and fMRI] and can be found in different species and during a variety of different task states. Together, this suggests that the hierarchy of INT is central to the temporal integration (combining successive stimuli) and segregation (separating successive stimuli) of external inputs from the environment, leading to temporal segmentation and prediction in perception and cognition.

Intrinsic timescales across the basal ganglia

Recent studies have shown that temporal stability of the neuronal activity over time can be estimated by the structure of the spike-count autocorrelation of neuronal populations. This estimation, called the intrinsic timescale, has been computed for several cortical areas and can be used to propose a cortical hierarchy reflecting a scale of temporal receptive windows between areas. In this study, we performed an autocorrelation analysis on neuronal populations of three basal ganglia (BG) nuclei, including the striatum and the subthalamic nucleus (STN), the input structures of the BG, and the external globus pallidus (GPe). The analysis was performed during the baseline period of a motivational visuomotor task in which monkeys had to apply different amounts of force to receive different amounts of reward. We found that the striatum and the STN have longer intrinsic timescales than the GPe. Moreover, our results allow for the placement of these subcortical structures within the already-defined scale of cortical temporal receptive windows. Estimates of intrinsic timescales are important in adding further constraints in the development of computational models of the complex dynamics among these nuclei and throughout cortico-BG-thalamo-cortical loops.

The brain and its time: intrinsic neural timescales are key for input processing

We process and integrate multiple timescales into one meaningful whole. Recent evidence suggests that the brain displays a complex multiscale temporal organization. Different regions exhibit different timescales as described by the concept of intrinsic neural timescales (INT); however, their function and neural mechanisms remains unclear. We review recent literature on INT and propose that they are key for input processing. Specifically, they are shared across different species, i.e., input sharing. This suggests a role of INT in encoding inputs through matching the inputs' stochastics with the ongoing temporal statistics of the brain's neural activity, i.e., input encoding. Following simulation and empirical data, we point out input integration versus segregation and input sampling as key temporal mechanisms of input processing. This deeply grounds the brain within its environmental and evolutionary context. It carries major implications in understanding mental features and psychiatric disorders, as well as going beyond the brain in integrating timescales into artificial intelligence.

Extracting single-trial neural interaction using latent dynamical systems model

In systems neuroscience, advances in simultaneous recording technology have helped reveal the population dynamics that underlie the complex neural correlates of animal behavior and cognitive processes. To investigate these correlates, neural interactions are typically abstracted from spike trains of pairs of neurons accumulated over the course of many trials. However, the resultant averaged values do not lead to understanding of neural computation in which the responses of populations are highly variable even under identical external conditions. Accordingly, neural interactions within the population also show strong fluctuations. In the present study, we introduce an analysis method reflecting the temporal variation of neural interactions, in which cross-correlograms on rate estimates are applied via a latent dynamical systems model. Using this method, we were able to predict time-varying neural interactions within a single trial. In addition, the pairwise connections estimated in our analysis increased along behavioral epochs among neurons categorized within similar functional groups. Thus, our analysis method revealed that neurons in the same groups communicate more as the population gets involved in the assigned task. We also showed that the characteristics of neural interaction from our model differ from the results of a typical model employing cross-correlation coefficients. This suggests that our model can extract nonoverlapping information about network topology, unlike the typical model.

Macaque monkeys learn and perform a non-match-to-goal task using an automated home cage training procedure

In neurophysiology, nonhuman primates represent an important model for studying the brain. Typically, monkeys are moved from their home cage to an experimental room daily, where they sit in a primate chair and interact with electronic devices. Refining this procedure would make the researchers' work easier and improve the animals' welfare. To address this issue, we used home-cage training to train two macaque monkeys in a non-match-to-goal task, where each trial required a switch from the choice made in the previous trial to obtain a reward. The monkeys were tested in two versions of the task, one in which they acted as the agent in every trial and one in which some trials were completed by a "ghost agent". We evaluated their involvement in terms of their performance and their interaction with the apparatus. Both monkeys were able to maintain a constant involvement in the task with good, stable performance within sessions in both versions of the task. Our study confirms the feasibility of home-cage training and demonstrates that even with challenging tasks, monkeys can complete a large number of trials at a high performance level, which is a prerequisite for electrophysiological studies of monkey behavior.

A Diversity of Intrinsic Timescales Underlie Neural Computations

Neural processing occurs across a range of temporal scales. To facilitate this, the brain uses fast-changing representations reflecting momentary sensory input alongside more temporally extended representations, which integrate across both short and long temporal windows. The temporal flexibility of these representations allows animals to behave adaptively. Short temporal windows facilitate adaptive responding in dynamic environments, while longer temporal windows promote the gradual integration of information across time. In the cognitive and motor domains, the brain sets overarching goals to be achieved within a long temporal window, which must be broken down into sequences of actions and precise movement control processed across much shorter temporal windows. Previous human neuroimaging studies and large-scale artificial network models have ascribed different processing timescales to different cortical regions, linking this to each region's position in an anatomical hierarchy determined by patterns of inter-regional connectivity. However, even within cortical regions, there is variability in responses when studied with single-neuron electrophysiology. Here, we review a series of recent electrophysiology experiments that demonstrate the heterogeneity of temporal receptive fields at the level of single neurons within a cortical region. This heterogeneity appears functionally relevant for the computations that neurons perform during decision-making and working memory. We consider anatomical and biophysical mechanisms that may give rise to a heterogeneity of timescales, including recurrent connectivity, cortical layer distribution, and neurotransmitter receptor expression. Finally, we reflect on the computational relevance of each brain region possessing a heterogeneity of neuronal timescales. We argue that this architecture is of particular importance for sensory, motor, and cognitive computations.

Brief localised monocular deprivation in adults alters binocular rivalry predominance retinotopically and reduces spatial inhibition

Short-term deprivation (2.5 h) of an eye has been shown to boost its relative ocular dominance in young adults. Here, we show that a much shorter deprivation period (3–6 min) produces a similar paradoxical boost that is retinotopic and reduces spatial inhibition on neighbouring, non-deprived areas. Partial deprivation was conducted in the left hemifield, central vision or in an annular region, later assessed with a binocular rivalry tracking procedure. Post-deprivation, dominance of the deprived eye increased when rivalling images were within the deprived retinotopic region, but not within neighbouring, non-deprived areas where dominance was dependent on the correspondence between the orientation content of the stimuli presented in the deprived and that of the stimuli presented in non-deprived areas. Together, these results accord with other deprivation studies showing V1 activity changes and reduced GABAergic inhibition.

Multiple timescales of neural dynamics and integration of task-relevant signals across cortex

A long-lasting challenge in neuroscience has been to find a set of principles that could be used to organize the brain into distinct areas with specific functions. Recent studies have proposed the orderly progression in the time constants of neural dynamics as an organizational principle of cortical computations. However, relationships between these timescales and their dependence on response properties of individual neurons are unknown, making it impossible to determine how mechanisms underlying such a computational principle are related to other aspects of neural processing. Here, we developed a comprehensive method to simultaneously estimate multiple timescales in neuronal dynamics and integration of task-relevant signals along with selectivity to those signals. By applying our method to neural and behavioral data during a dynamic decision-making task, we found that most neurons exhibited multiple timescales in their response, which consistently increased from parietal to prefrontal and cingulate cortex. While predicting rates of behavioral adjustments, these timescales were not correlated across individual neurons in any cortical area, resulting in independent parallel hierarchies of timescales. Additionally, none of these timescales depended on selectivity to task-relevant signals. Our results not only suggest the existence of multiple canonical mechanisms for increasing timescales of neural dynamics across cortex but also point to additional mechanisms that allow decorrelation of these timescales to enable more flexibility.