Invariant timescale hierarchy across the cortical somatosensory network

Multimodal mapping of macaque monkey somatosensory cortex

The somatosensory cortex is a brain region responsible for receiving and processing sensory information from across the body and is structurally and functionally heterogeneous. Since the chemoarchitectonic segregation of the cerebral cortex can be revealed by transmitter receptor distribution patterns, by using a quantitative multireceptor architectonical analysis, we determined the number and extent of distinct areas of the macaque somatosensory cortex. We identified three architectonically distinct cortical entities within the primary somatosensory cortex (i.e., 3bm, 3bli, 3ble), four within the anterior parietal cortex (i.e., 3am, 3al, 1 and 2) and six subdivisions (i.e., S2l, S2m, PVl, PVm, PRl and PRm) within the lateral fissure. We provide an ultra-high resolution 3D atlas of macaque somatosensory areas in stereotaxic space, which integrates cyto- and receptor architectonic features of identified areas. Multivariate analyses of the receptor fingerprints revealed four clusters of identified areas based on the degree of (dis)similarity of their receptor architecture. Each of these clusters can be associated with distinct levels of somatosensory processing, further demonstrating that the functional segregation of cortical areas is underpinned by differences in their molecular organization.

Thalamocortical interactions shape hierarchical neural variability during stimulus perception

The brain is organized hierarchically to process sensory signals. But, how do functional connections within and across areas contribute to this hierarchical order? We addressed this problem in the thalamocortical network, while monkeys detected vibrotactile stimulus. During this task, we quantified neural variability and directed functional connectivity in simultaneously recorded neurons sharing the cutaneous receptive field within and across VPL and areas 3b and 1. Before stimulus onset, VPL and area 3b exhibited similar fast dynamics while area 1 showed slower timescales. During the stimulus presence, inter-trial neural variability increased along the network VPL-3b-1 while VPL established two main feedforward pathways with areas 3b and 1 to process the stimulus. This lower variability of VPL and area 3b was found to regulate feedforward thalamocortical pathways. Instead, intra-cortical interactions were only anticipated by higher intrinsic timescales in area 1. Overall, our results provide evidence of hierarchical functional roles along the thalamocortical network.

Interacting ventral temporal gradients of timescales and functional connectivity and their relationships to visual behavior

Cortical gradients in endogenous and stimulus-evoked neurodynamic timescales, and long-range cortical interactions, provide organizational constraints to the brain and influence neural populations' roles in cognition. It is unclear how these functional gradients interrelate and which influence behavior. Here, intracranial recordings from 4,090 electrode contacts in 35 individuals map gradients of neural timescales and functional connectivity to assess their interactions along category-selective ventral temporal cortex. Endogenous and stimulus-evoked information processing timescales were not significantly correlated with one another suggesting that local neural timescales are context dependent and may arise through distinct neurophysiological mechanisms. Endogenous neural timescales correlated with functional connectivity even after removing the effects of shared anatomical gradients. Neural timescales and functional connectivity correlated with how strongly a population's activity predicted behavior in a simple visual task. These results suggest both interrelated and distinct neurophysiological processes give rise to different functional connectivity and neural timescale gradients, which together influence behavior.

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.

Evidence for Multiscale Multiplexed Representation of Visual Features in EEG

Distinct neural processes such as sensory and memory processes are often encoded over distinct timescales of neural activations. Animal studies have shown that this multiscale coding strategy is also implemented for individual components of a single process, such as individual features of a multifeature stimulus in sensory coding. However, the generalizability of this encoding strategy to the human brain has remained unclear. We asked if individual features of visual stimuli were encoded over distinct timescales. We applied a multiscale time-resolved decoding method to electroencephalography (EEG) collected from human subjects presented with grating visual stimuli to estimate the timescale of individual stimulus features. We observed that the orientation and color of the stimuli were encoded in shorter timescales, whereas spatial frequency and the contrast of the same stimuli were encoded in longer timescales. The stimulus features appeared in temporally overlapping windows along the trial supporting a multiplexed coding strategy. These results provide evidence for a multiplexed, multiscale coding strategy in the human visual system.

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.

Distributed Coding of Evidence Accumulation across the Mouse Brain Using Microcircuits with a Diversity of Timescales

The gradual accumulation of noisy evidence for or against options is the main step in the perceptual decision-making process. Using brain-wide electrophysiological recording in mice (Steinmetz et al., 2019), we examined neural correlates of evidence accumulation across brain areas. We demonstrated that the neurons with drift-diffusion model (DDM)-like firing rate activity (i.e., evidence-sensitive ramping firing rate) were distributed across the brain. Exploring the underlying neural mechanism of evidence accumulation for the DDM-like neurons revealed different accumulation mechanisms (i.e., single and race) both within and across the brain areas. Our findings support the hypothesis that evidence accumulation is happening through multiple integration mechanisms in the brain. We further explored the timescale of the integration process in the single and race accumulator models. The results demonstrated that the accumulator microcircuits within each brain area had distinct properties in terms of their integration timescale, which were organized hierarchically across the brain. These findings support the existence of evidence accumulation over multiple timescales. Besides the variability of integration timescale across the brain, a heterogeneity of timescales was observed within each brain area as well. We demonstrated that this variability reflected the diversity of microcircuit parameters, such that accumulators with longer integration timescales had higher recurrent excitation strength.

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.

Hierarchical unimodal processing within the primary somatosensory cortex during a bimodal detection task

Do sensory cortices process more than one sensory modality? To answer these questions, scientists have generated a wide variety of studies at distinct space-time scales in different animal models, and often shown contradictory conclusions. Some conclude that this process occurs in early sensory cortices, but others that this occurs in areas central to sensory cortices. Here, we sought to determine whether sensory neurons process and encode physical stimulus properties of different modalities (tactile and acoustic). For this, we designed a bimodal detection task where the senses of touch and hearing compete from trial to trial. Two Rhesus monkeys performed this novel task, while neural activity was recorded in areas 3b and 1 of the primary somatosensory cortex (S1). We analyzed neurons' coding properties and variability, organizing them by their receptive field's position relative to the stimulation zone. Our results indicate that neurons of areas 3b and 1 are unimodal, encoding only the tactile modality in both the firing rate and variability. Moreover, we found that neurons in area 3b carried more information about the periodic stimulus structure than those in area 1, possessed lower response and coding latencies, and had a lower intrinsic time scale. In sum, these differences reveal a hidden processing-based hierarchy. Finally, using a powerful nonlinear dimensionality reduction algorithm, we show that the activity from areas 3b and 1 can be separated, establishing a clear division in the functionality of these two subareas of S1.

Dimensionality reduction and recurrence analysis reveal hidden structures of striatal pathological states

A pipeline is proposed here to describe different features to study brain microcircuits on a histological scale using multi-scale analyses, including the uniform manifold approximation and projection (UMAP) dimensional reduction technique and modularity algorithm to identify neuronal ensembles, Runs tests to show significant ensembles activation, graph theory to show trajectories between ensembles, and recurrence analyses to describe how regular or chaotic ensembles dynamics are. The data set includes ex-vivo NMDA-activated striatal tissue in control conditions as well as experimental models of disease states: decorticated, dopamine depleted, and L-DOPA-induced dyskinetic rodent samples. The goal was to separate neuronal ensembles that have correlated activity patterns. The pipeline allows for the demonstration of differences between disease states in a brain slice. First, the ensembles were projected in distinctive locations in the UMAP space. Second, graphs revealed functional connectivity between neurons comprising neuronal ensembles. Third, the Runs test detected significant peaks of coactivity within neuronal ensembles. Fourth, significant peaks of coactivity were used to show activity transitions between ensembles, revealing recurrent temporal sequences between them. Fifth, recurrence analysis shows how deterministic, chaotic, or recurrent these circuits are. We found that all revealed circuits had recurrent activity except for the decorticated circuits, which tended to be divergent and chaotic. The Parkinsonian circuits exhibit fewer transitions, becoming rigid and deterministic, exhibiting a predominant temporal sequence that disrupts transitions found in the controls, thus resembling the clinical signs of rigidity and paucity of movements. Dyskinetic circuits display a higher recurrence rate between neuronal ensembles transitions, paralleling clinical findings: enhancement in involuntary movements. These findings confirm that looking at neuronal circuits at the histological scale, recording dozens of neurons simultaneously, can show clear differences between control and diseased striatal states: "fingerprints" of the disease states. Therefore, the present analysis is coherent with previous ones of striatal disease states, showing that data obtained from the tissue are robust. At the same time, it adds heuristic ways to interpret circuitry activity in different states.

Modular Organization of Signal Transmission in Primate Somatosensory Cortex

Axonal patches are known as the major sites of synaptic connections in the cerebral cortex of higher order mammals. However, the functional role of these patches is highly debated. Patches are formed by populations of nearby neurons in a topographic manner and are recognized as the termination fields of long-distance lateral connections within and between cortical areas. In addition, axons form numerous boutons that lie outside the patches, whose function is also unknown. To better understand the functional roles of these two distinct populations of boutons, we compared individual and collective morphological features of axons within and outside the patches of intra-areal, feedforward, and feedback pathways by way of tract tracing in the somatosensory cortex of New World monkeys. We found that, with the exception of tortuosity, which is an invariant property, bouton spacing and axonal convergence properties differ significantly between axons within patch and no-patch domains. Principal component analyses corroborated the clustering of axons according to patch formation without any additional effect by the type of pathway or laminar distribution. Stepwise logistic regression identified convergence and bouton density as the best predictors of patch formation. These findings support that patches are specific sites of axonal convergence that promote the synchronous activity of neuronal populations. On the other hand, no-patch domains could form a neuroanatomical substrate to diversify the responses of cortical neurons.

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.

The diversity and specificity of functional connectivity across spatial and temporal scales

Macroscopic neuroimaging modalities in humans have revealed the organization of brain-wide activity into distributed functional networks that re-organize according to behavioral demands. However, the inherent coarse-graining of macroscopic measurements conceals the diversity and specificity in responses and connectivity of many individual neurons contained in each local region. New invasive approaches in animals enable recording and manipulating neural activity at meso- and microscale resolution, with cell-type specificity and temporal precision down to milliseconds. Determining how brain-wide activity patterns emerge from interactions across spatial and temporal scales will allow us to identify the key circuit mechanisms contributing to global brain states and how the dynamic activity of these states enables adaptive behavior.

Neural Basis of Somatosensory Spatial and Temporal Discrimination in Humans: The Role of Sensory Detection

While detecting somatic stimuli from the external environment, an accurate determination of their spatial and temporal properties is essential for human behavior. Whether and how detection relates to human capacity for somatosensory spatial discrimination (SD) and temporal discrimination (TD) remains unclear. Here, participants underwent functional magnetic resonance imaging scanning when simply detecting vibrotactile stimuli of the leg, judging their location (SD), or deciding their number in time (TD). By conceptualizing tactile discrimination as consisting of detection and determination processes, we found that tactile detection elicited activation specifically involved in SD within the right inferior and superior parietal lobules, 2 regions previously implicated in the control of spatial attention. These 2 regions remained activated in the determination process, during which functional connectivity between these 2 regions predicted individual SD ability. In contrast, tactile detection produced little activation specifically related to TD. Participants' TD ability was implemented in brain regions implicated in coding temporal structures of somatic stimuli (primary somatosensory cortex) and time estimation (anterior cingulate, pre-supplementary motor area, and putamen). Together, our findings indicate a close link between somatosensory detection and SD (but not TD) at the neural level, which aids in explaining why we can promptly respond toward detected somatic stimuli.

A continuum of invariant sensory and behavioral-context perceptual coding in secondary somatosensory cortex

A crucial role of cortical networks is the conversion of sensory inputs into perception. In the cortical somatosensory network, neurons of the primary somatosensory cortex (S1) show invariant sensory responses, while frontal lobe neuronal activity correlates with the animal's perceptual behavior. Here, we report that in the secondary somatosensory cortex (S2), neurons with invariant sensory responses coexist with neurons whose responses correlate with perceptual behavior. Importantly, the vast majority of the neurons fall along a continuum of combined sensory and categorical dynamics. Furthermore, during a non-demanding control task, the sensory responses remain unaltered while the sensory information exhibits an increase. However, perceptual responses and the associated categorical information decrease, implicating a task context-dependent processing mechanism. Conclusively, S2 neurons exhibit intriguing dynamics that are intermediate between those of S1 and frontal lobe. Our results contribute relevant evidence about the role that S2 plays in the conversion of touch into perception.