A Large-Scale Circuit Mechanism for Hierarchical Dynamical Processing in the Primate Cortex

Distribution of calbindin-positive neurons across areas and layers of the marmoset cerebral cortex

The diversity of the mammalian cerebral cortex demands technical approaches to map the spatial distribution of neurons with different biochemical identities. This issue is magnified in the case of the primate cortex, characterized by a large number of areas with distinctive cytoarchitectures. To date, no full map of the distribution of cells expressing a specific protein has been reported for the cortex of any primate. Here we have charted the 3-dimensional distribution of neurons expressing the calcium-binding protein calbindin (CB+ neurons) across the entire marmoset cortex, using a combination of immunohistochemistry, automated cell identification, computerized reconstruction, and cytoarchitecture-aware registration. CB+ neurons formed a heterogeneous population, which together corresponded to 10-20% of the cortical neurons. They occurred in higher proportions in areas corresponding to low hierarchical levels of processing, such as sensory cortices. Although CB+ neurons were concentrated in the supragranular and granular layers, there were clear global trends in their laminar distribution. For example, their relative density in infragranular layers increased with hierarchical level along sensorimotor processing streams, and their density in layer 4 was lower in areas involved in sensorimotor integration, action planning and motor control. These results reveal new quantitative aspects of the cytoarchitectural organization of the primate cortex, and demonstrate an approach to mapping the full distribution of neurochemically distinct cells throughout the brain which is readily applicable to most other mammalian species.

Model-agnostic neural mean field with a data-driven transfer function

As one of the most complex systems known to science, modeling brain behavior and function is both fascinating and extremely difficult. Empirical data is increasingly available from ex vivo human brain organoids and surgical samples, as well as in vivo animal models, so the problem of modeling the behavior of large-scale neuronal systems is more relevant than ever. The statistical physics concept of a mean-field model offers a tractable way to bridge the gap between single-neuron and population-level descriptions of neuronal activity, by modeling the behavior of a single representative neuron and extending this to the population. However, existing neural mean-field methods typically either take the limit of small interaction sizes, or are applicable only to the specific neuron models for which they were derived. This paper derives a mean-field model by fitting a transfer function called Refractory SoftPlus, which is simple yet applicable to a broad variety of neuron types. The transfer function is fitted numerically to simulated spike time data, and is entirely agnostic to the underlying neuronal dynamics. The resulting mean-field model predicts the response of a network of randomly connected neurons to a time-varying external stimulus with a high degree of accuracy. Furthermore, it enables an accurate approximate bifurcation analysis as a function of the level of recurrent input. This model does not assume large presynaptic rates or small postsynaptic potential size, allowing mean-field models to be developed even for populations with large interaction terms.

Hierarchical gradients of multiple timescales in the mammalian forebrain

Many anatomical and physiological features of cortical circuits, ranging from the biophysical properties of synapses to the connectivity patterns among different neuron types, exhibit consistent variation along the hierarchical axis from sensory to association areas. Notably, the scale of temporal correlation of neural activity at rest, known as the intrinsic timescale, increases systematically along this hierarchy in both primates and rodents, analogous to the growing scale and complexity of spatial receptive fields. However, how the timescales for task-related activity vary across brain regions and whether their hierarchical organization appears consistently across different mammalian species remain unexplored. Here, we show that both the intrinsic timescale and the timescales of task-related activity follow a similar hierarchical gradient in the cortices of monkeys, rats, and mice. We also found that these timescales covary similarly in both the cortex and basal ganglia, whereas the timescales of thalamic activity are shorter than cortical timescales and do not conform to the hierarchical order predicted by their cortical projections. These results suggest that the hierarchical gradient of cortical timescales might be a universal feature of intra-cortical circuits in the mammalian brain.

Medial to lateral frontal functional connectivity mapping reveals the organization of cingulate cortex

The functional organization of the frontal lobe is a source of debate, focusing on broad functional subdivisions, large-scale networks, or local refined specificities. Multiple neurocognitive models have tried to explain how functional interactions between cingulate and lateral frontal regions contribute to decision making and cognitive control, but their neuroanatomical bases remain unclear. We provide a detailed description of the functional connectivity between cingulate and lateral frontal regions using resting-state functional MRI in rhesus macaques. The analysis focuses on the functional connectivity of the rostral part of the cingulate sulcus with the lateral frontal cortex. Data-driven and seed-based analysis revealed three clusters within the cingulate sulcus organized along the rostro-caudal axis: the anterior, mid, and posterior clusters display increased functional connectivity with, respectively, the anterior lateral prefrontal regions, face-eye lateral frontal motor cortical areas, and hand lateral frontal motor cortex. The location of these clusters can be predicted in individual subjects based on morphological landmarks. These results suggest that the anterior cluster corresponds to the anterior cingulate cortex, whereas the posterior clusters correspond to the face-eye and hand cingulate motor areas within the anterior midcingulate cortex. These data provide a comprehensive framework to identify cingulate subregions based on functional connectivity and local organization.

Signatures of hierarchical temporal processing in the mouse visual system

A core challenge for the brain is to process information across various timescales. This could be achieved by a hierarchical organization of temporal processing through intrinsic mechanisms (e.g., recurrent coupling or adaptation), but recent evidence from spike recordings of the rodent visual system seems to conflict with this hypothesis. Here, we used an optimized information-theoretic and classical autocorrelation analysis to show that information- and correlation timescales of spiking activity increase along the anatomical hierarchy of the mouse visual system under visual stimulation, while information-theoretic predictability decreases. Moreover, intrinsic timescales for spontaneous activity displayed a similar hierarchy, whereas the hierarchy of predictability was stimulus-dependent. We could reproduce these observations in a basic recurrent network model with correlated sensory input. Our findings suggest that the rodent visual system employs intrinsic mechanisms to achieve longer integration for higher cortical areas, while simultaneously reducing predictability for an efficient neural code.

Dynamic Formation of a Posterior-to-Anterior Peak-Alpha-Frequency Gradient Driven by Two Distinct Processes

Peak-alpha frequency varies across individuals and mental states, but it also forms a negative gradient from posterior to anterior regions in association with increases in cortical thickness and connectivity, reflecting a cortical hierarchy in temporal integration. Tracking the spatial standard deviation of peak-alpha frequency in scalp EEG, we observed that a posterior-to-anterior gradient dynamically formed and dissolved. Periods of high spatial standard deviation yielded robustly negative posterior-to-anterior gradients-the "gradient state"-while periods of low spatial standard deviation yielded globally converged peak-alpha frequency-the "uniform state." The state variations were characterized by a combination of slow (0.3-0.5 Hz) oscillations and random-walk-like fluctuations. They were relatively independently correlated with peak-alpha frequency variations in anterior regions and peak-alpha power variations in central regions driven by posterior regions (together accounting for ∼50% of the state variations), suggesting that two distinct mechanisms modulate the state variations: an anterior mechanism that directly adjusts peak-alpha frequencies and a posterior-central mechanism that indirectly adjusts them by influencing synchronization. The state variations likely reflect general operations as their spatiotemporal characteristics remained unchanged while participants engaged in a variety of tasks (breath focus, vigilance, working memory, mental arithmetic, and generative thinking) with their eyes closed or watched a silent nature video. The ongoing state variations may dynamically balance two global processing modes, one that facilitates greater temporal integration (and potentially also information influx) toward anterior regions in the gradient state and the other that facilitates flexible global communication (via phase locking) in the uniform state.

Neural dynamics underlying minute-timescale persistent behavior in the human brain

The ability to pursue long-term goals relies on a representations of task context that can both be maintained over long periods of time and switched flexibly when goals change. Little is known about the neural substrate for such minute-scale maintenance of task sets. Utilizing recordings in neurosurgical patients, we examined how groups of neurons in the human medial frontal cortex and hippocampus represent task contexts. When cued explicitly, task context was encoded in both brain areas and changed rapidly at task boundaries. Hippocampus exhibited a temporally dynamic code with fast decorrelation over time, preventing cross-temporal generalization. Medial frontal cortex exhibited a static code that decorrelated slowly, allowing generalization across minutes of time. When task context needed to be inferred as a latent variable, hippocampus encoded task context with a static code. These findings reveal two possible regimes for encoding minute-scale task-context representations that were engaged differently based on task demands.

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.

Distinct feedforward and feedback pathways for cell-type specific attention effects

Selective attention is thought to depend on enhanced firing activity in extrastriate areas. Theories suggest that this enhancement depends on selective inter-areal communication via gamma (30-80 Hz) phase-locking. To test this, we simultaneously recorded from different cell types and cortical layers of macaque V1 and V4. We find that while V1-V4 gamma phase-locking between local field potentials increases with attention, the V1 gamma rhythm does not engage V4 excitatory-neurons, but only fast-spiking interneurons in L4 of V4. By contrast, attention enhances V4 spike-rates in both excitatory and inhibitory cells, most strongly in L2/3. The rate increase in L2/3 of V4 precedes V1 in time. These findings suggest enhanced signal transmission with attention does not depend on inter-areal gamma phase-locking and show that the endogenous gamma rhythm has cell-type- and layer-specific effects on downstream target areas. Similar findings were made in the mouse visual system, based on opto-tagging of identified interneurons.

Atypical intrinsic neural timescale in the left angular gyrus in Alzheimer's disease

Alzheimer's disease is characterized by cognitive impairment and progressive brain atrophy. Recent human neuroimaging studies reported atypical anatomical and functional changes in some regions in the default mode network in patients with Alzheimer's disease, but which brain area of the default mode network is the key region whose atrophy disturbs the entire network activity and consequently contributes to the symptoms of the disease remains unidentified. Here, in this case-control study, we aimed to identify crucial neural regions that mediated the phenotype of Alzheimer's disease, and as such, we examined the intrinsic neural timescales-a functional metric to evaluate the capacity to integrate diverse neural information-and grey matter volume of the regions in the default mode network using resting-state functional MRI images and structural MRI data obtained from individuals with Alzheimer's disease and cognitively typical people. After confirming the atypically short neural timescale of the entire default mode network in Alzheimer's disease and its link with the symptoms of the disease, we found that the shortened neural timescale of the default mode network was associated with the aberrantly short neural timescale of the left angular gyrus. Moreover, we revealed that the shortened neural timescale of the angular gyrus was correlated with the atypically reduced grey matter volume of this parietal region. Furthermore, we identified an association between the neural structure, brain function and symptoms and proposed a model in which the reduced grey matter volume of the left angular gyrus shortened the intrinsic neural time of the region, which then destabilized the entire neural timescale of the default mode network and resultantly contributed to cognitive decline in Alzheimer's disease. These findings highlight the key role of the left angular gyrus in the anatomical and functional aetiology of Alzheimer's disease.

The primate cortical LFP exhibits multiple spectral and temporal gradients and widespread task dependence during visual short-term memory

Although cognitive functions are hypothesized to be mediated by synchronous neuronal interactions in multiple frequency bands among widely distributed cortical areas, we still lack a basic understanding of the distribution and task dependence of oscillatory activity across the cortical map. Here, we ask how the spectral and temporal properties of the local field potential (LFP) vary across the primate cerebral cortex, and how they are modulated during visual short-term memory. We measured the LFP from 55 cortical areas in two macaque monkeys while they performed a visual delayed match to sample task. Analysis of peak frequencies in the LFP power spectra reveals multiple discrete frequency bands between 3 and 80 Hz that differ between the two monkeys. The LFP power in each band, as well as the sample entropy, a measure of signal complexity, display distinct spatial gradients across the cortex, some of which correlate with reported spine counts in cortical pyramidal neurons. Cortical areas can be robustly decoded using a small number of spectral and temporal parameters, and significant task-dependent increases and decreases in spectral power occur in all cortical areas. These findings reveal pronounced, widespread, and spatially organized gradients in the spectral and temporal activity of cortical areas. Task-dependent changes in cortical activity are globally distributed, even for a simple cognitive task.NEW & NOTEWORTHY We recorded extracellular electrophysiological signals from roughly the breadth and depth of a cortical hemisphere in nonhuman primates (NHPs) performing a visual memory task. Analyses of the band-limited local field potential (LFP) power displayed widespread, frequency-dependent cortical gradients in spectral power. Using a machine learning classifier, these features allowed robust cortical area decoding. Further task dependence in LFP power were found to be widespread, indicating large-scale gradients of LFP activity, and task-related activity.

Altered intrinsic neural timescales and neurotransmitter activity in males with tobacco use disorder

Previous researches of tobacco use disorder (TUD) has overlooked the hierarchy of cortical functions and single modality design separated the relationship between macroscopic neuroimaging aberrance and microscopic molecular basis. At present, intrinsic timescale gradient of TUD and its molecular features are not fully understood. Our study recruited 146 male subjects, including 44 heavy smokers, 50 light smokers and 52 non-smokers, then obtained their rs-fMRI data and clinical scales related to smoking. Intrinsic neural timescale (INT) method was performed to describe how long neural information was stored in a brain region by calculating the autocorrelation function (ACF) of each voxel to examine the difference in the ability of information integration among the three groups. Then, correlation analyses were conducted to explore the relationship between INT abnormalities and clinical scales of smokers. Finally, cross-modal JuSpace toolbox was used to investigate the association between INT aberrance and the expression of specific receptor/transporters. Compared to healthy controls, TUD subjects displayed decreased INT in control network (CN), default mode network (DMN), sensorimotor areas and visual cortex, and such trend of decreasing INT was more pronounced in heavy smokers. Moreover, various neurotransmitters (including dopaminergic, acetylcholine and μ-opioid receptors) were involved in the molecular mechanism of timescale decreasing and differed in heavy and light smokers. These findings supplied novel insights into the brain functional aberrance in TUD from an intrinsic neural dynamic perspective and confirm INT was a potential neurobiological marker. And also established the connection between macroscopic imaging aberrance and microscopic molecular changes in TUD.

A sensory-motor theory of the neocortex

Recent neurophysiological and neuroanatomical studies suggest a close interaction between sensory and motor processes across the neocortex. Here, I propose that the neocortex implements active predictive coding (APC): each cortical area estimates both latent sensory states and actions (including potentially abstract actions internal to the cortex), and the cortex as a whole predicts the consequences of actions at multiple hierarchical levels. Feedback from higher areas modulates the dynamics of state and action networks in lower areas. I show how the same APC architecture can explain (1) how we recognize an object and its parts using eye movements, (2) why perception seems stable despite eye movements, (3) how we learn compositional representations, for example, part-whole hierarchies, (4) how complex actions can be planned using simpler actions, and (5) how we form episodic memories of sensory-motor experiences and learn abstract concepts such as a family tree. I postulate a mapping of the APC model to the laminar architecture of the cortex and suggest possible roles for cortico-cortical and cortico-subcortical pathways.

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.

A measure centrality index for systematic empirical comparison of consciousness theories

Consciousness science is marred by disparate constructs and methodologies, making it challenging to systematically compare theories. This foundational crisis casts doubts on the scientific character of the field itself. Addressing it, we propose a framework for systematically comparing consciousness theories by introducing a novel inter-theory classification interface, the Measure Centrality Index (MCI). Recognizing its gradient distribution, the MCI assesses the degree of importance a specific empirical measure has for a given consciousness theory. We apply the MCI to probe how the empirical measures of the Global Neuronal Workspace Theory (GNW), Integrated Information Theory (IIT), and Temporospatial Theory of Consciousness (TTC) would fare within the context of the other two. We demonstrate that direct comparison of IIT, GNW, and TTC is meaningful and valid for some measures like Lempel-Ziv Complexity (LZC), Autocorrelation Window (ACW), and possibly Mutual Information (MI). In contrast, it is problematic for others like the anatomical and physiological neural correlates of consciousness (NCC) due to their MCI-based differential weightings within the structure of the theories. In sum, we introduce and provide proof-of-principle of a novel systematic method for direct inter-theory empirical comparisons, thereby addressing isolated evolution of theories and confirmatory bias issues in the state-of-the-art neuroscience of consciousness.

The meso-connectomes of mouse, marmoset, and macaque: network organization and the emergence of higher cognition

The recent publications of the inter-areal connectomes for mouse, marmoset, and macaque cortex have allowed deeper comparisons across rodent vs. primate cortical organization. In general, these show that the mouse has very widespread, "all-to-all" inter-areal connectivity (i.e. a "highly dense" connectome in a graph theoretical framework), while primates have a more modular organization. In this review, we highlight the relevance of these differences to function, including the example of primary visual cortex (V1) which, in the mouse, is interconnected with all other areas, therefore including other primary sensory and frontal areas. We argue that this dense inter-areal connectivity benefits multimodal associations, at the cost of reduced functional segregation. Conversely, primates have expanded cortices with a modular connectivity structure, where V1 is almost exclusively interconnected with other visual cortices, themselves organized in relatively segregated streams, and hierarchically higher cortical areas such as prefrontal cortex provide top-down regulation for specifying precise information for working memory storage and manipulation. Increased complexity in cytoarchitecture, connectivity, dendritic spine density, and receptor expression additionally reveal a sharper hierarchical organization in primate cortex. Together, we argue that these primate specializations permit separable deconstruction and selective reconstruction of representations, which is essential to higher cognition.

Temporal dynamics of short-term neural adaptation across human visual cortex

Neural responses in visual cortex adapt to prolonged and repeated stimuli. While adaptation occurs across the visual cortex, it is unclear how adaptation patterns and computational mechanisms differ across the visual hierarchy. Here we characterize two signatures of short-term neural adaptation in time-varying intracranial electroencephalography (iEEG) data collected while participants viewed naturalistic image categories varying in duration and repetition interval. Ventral- and lateral-occipitotemporal cortex exhibit slower and prolonged adaptation to single stimuli and slower recovery from adaptation to repeated stimuli compared to V1-V3. For category-selective electrodes, recovery from adaptation is slower for preferred than non-preferred stimuli. To model neural adaptation we augment our delayed divisive normalization (DN) model by scaling the input strength as a function of stimulus category, enabling the model to accurately predict neural responses across multiple image categories. The model fits suggest that differences in adaptation patterns arise from slower normalization dynamics in higher visual areas interacting with differences in input strength resulting from category selectivity. Our results reveal systematic differences in temporal adaptation of neural population responses between lower and higher visual brain areas and show that a single computational model of history-dependent normalization dynamics, fit with area-specific parameters, accounts for these differences.

Taking time to compose thoughts with prefrontal schemata

Under what conditions can prefrontal cortex direct the composition of brain states, to generate coherent streams of thoughts? Using a simplified Potts model of cortical dynamics, crudely differentiated into two halves, we show that once activity levels are regulated, so as to disambiguate a single temporal sequence, whether the contents of the sequence are mainly determined by the frontal or by the posterior half, or by neither, depends on statistical parameters that describe its microcircuits. The frontal cortex tends to lead if it has more local attractors, longer lasting and stronger ones, in order of increasing importance. Its guidance is particularly effective to the extent that posterior cortices do not tend to transition from state to state on their own. The result may be related to prefrontal cortex enforcing its temporally-oriented schemata driving coherent sequences of brain states, unlike the atemporal "context" contributed by the hippocampus. Modelling a mild prefrontal (vs. posterior) lesion offers an account of mind-wandering and event construction deficits observed in prefrontal patients.

Potential of focal cortical dysplasia in migraine pathogenesis

Focal cortical dysplasias are abnormalities of the cerebral cortex associated with an elevated risk of neurological disturbances. Cortical spreading depolarization/depression is a correlate of migraine aura/headache and a trigger of migraine pain mechanisms. However, cortical spreading depolarization/depression is associated with cortical structural changes, which can be classified as transient focal cortical dysplasias. Migraine is reported to be associated with changes in various brain structures, including malformations and lesions in the cortex. Such malformations may be related to focal cortical dysplasias, which may play a role in migraine pathogenesis. Results obtained so far suggest that focal cortical dysplasias may belong to the causes and consequences of migraine. Certain focal cortical dysplasias may lower the threshold of cortical excitability and facilitate the action of migraine triggers. Migraine prevalence in epileptic patients is higher than in the general population, and focal cortical dysplasias are an established element of epilepsy pathogenesis. In this narrative/hypothesis review, we present mainly information on cortical structural changes in migraine, but studies on structural alterations in deep white matter and other brain regions are also presented. We develop the hypothesis that focal cortical dysplasias may be causally associated with migraine and link pathogeneses of migraine and epilepsy.

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.

Decreased intrinsic neural timescale in treatment-naïve adolescent depression

BACKGROUND

Major depressive disorder (MDD) is mainly characterized by its core dysfunction in higher-order brain cortices involved in emotional and cognitive processes, whose neurobiological basis remains unclear. In this study, we applied a relatively new developed resting-state functional magnetic resonance imaging (rs-fMRI) method of intrinsic neural timescale (INT), which reflects how long neural information is stored in a local brain area and reflects an ability of information integration, to investigate the local intrinsic neural dynamics using univariate and multivariate analyses in adolescent depression.

METHOD

Based on the rs-fMRI data of sixty-six treatment-naïve adolescents with MDD and fifty-two well-matched healthy controls (HCs), we calculated an INT by assessing the magnitude of autocorrelation of the resting-state brain activity, and then compared the difference of INT between the two groups. Correlation between abnormal INT and clinical features was performed. We also utilized multivariate pattern analysis to determine whether INT could differentiate MDD patients from HCs at the individual level.

RESULT

Compared with HCs, patients with MDD showed shorter INT widely distributed in cortical and partial subcortical regions. Interestingly, the decreased INT in the left hippocampus was related to disease severity of MDD. Furthermore, INT can distinguish MDD patients from HCs with the most discriminative regions located in the dorsolateral prefrontal cortex, angular, middle occipital gyrus, and cerebellar posterior lobe.

CONCLUSION

Our research aids in advancing understanding the brain abnormalities of treatment-naïve adolescents with MDD from the perspective of the local neural dynamics, highlighting the significant role of INT in understanding neurophysiological mechanisms. This study shows that the altered intrinsic timescales of local neural signals widely distributed in higher-order brain cortices regions may be the neurodynamic basis of cognitive and emotional disturbances in MDD patients, and provides preliminary support for the suggestion that these could be used to aid the identification of MDD patients in clinical practice.

Temporal dynamics of short-term neural adaptation across human visual cortex

Neural responses in visual cortex adapt to prolonged and repeated stimuli. While adaptation occurs across the visual cortex, it is unclear how adaptation patterns and computational mechanisms differ across the visual hierarchy. Here we characterize two signatures of short-term neural adaptation in time-varying intracranial electroencephalography (iEEG) data collected while participants viewed naturalistic image categories varying in duration and repetition interval. Ventral- and lateral-occipitotemporal cortex exhibit slower and prolonged adaptation to single stimuli and slower recovery from adaptation to repeated stimuli compared to V1-V3. For category-selective electrodes, recovery from adaptation is slower for preferred than non-preferred stimuli. To model neural adaptation we augment our delayed divisive normalization (DN) model by scaling the input strength as a function of stimulus category, enabling the model to accurately predict neural responses across multiple image categories. The model fits suggest that differences in adaptation patterns arise from slower normalization dynamics in higher visual areas interacting with differences in input strength resulting from category selectivity. Our results reveal systematic differences in temporal adaptation of neural population responses across the human visual hierarchy and show that a single computational model of history-dependent normalization dynamics, fit with area-specific parameters, accounts for these differences.

Bounded contribution of human early visual cortex to the topographic anisotropy in spatial extent perception

To interact successfully with objects, it is crucial to accurately perceive their spatial extent, an enclosed region they occupy in space. Although the topographic representation of space in the early visual cortex (EVC) has been favored as a neural correlate of spatial extent perception, its exact nature and contribution to perception remain unclear. Here, we inspect the topographic representations of human individuals' EVC and perception in terms of how much their anisotropy is influenced by the orientation (co-axiality) and radial position (radiality) of stimuli. We report that while the anisotropy is influenced by both factors, its direction is primarily determined by radiality in EVC but by co-axiality in perception. Despite this mismatch, the individual differences in both radial and co-axial anisotropy are substantially shared between EVC and perception. Our findings suggest that spatial extent perception builds on EVC's spatial representation but requires an additional mechanism to transform its topographic bias.

Model-Agnostic Neural Mean Field With The Refractory SoftPlus Transfer Function

Due to the complexity of neuronal networks and the nonlinear dynamics of individual neurons, it is challenging to develop a systems-level model which is accurate enough to be useful yet tractable enough to apply. Mean-field models which extrapolate from single-neuron descriptions to large-scale models can be derived from the neuron's transfer function, which gives its firing rate as a function of its synaptic input. However, analytically derived transfer functions are applicable only to the neurons and noise models from which they were originally derived. In recent work, approximate transfer functions have been empirically derived by fitting a sigmoidal curve, which imposes a maximum firing rate and applies only in the diffusion limit, restricting applications. In this paper, we propose an approximate transfer function called Refractory SoftPlus, which is simple yet applicable to a broad variety of neuron types. Refractory SoftPlus activation functions allow the derivation of simple empirically approximated mean-field models using simulation results, which enables prediction of the response of a network of randomly connected neurons to a time-varying external stimulus with a high degree of accuracy. These models also support an accurate approximate bifurcation analysis as a function of the level of recurrent input. Finally, the model works without assuming large presynaptic rates or small postsynaptic potential size, allowing mean-field models to be developed even for populations with large interaction terms.

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.

Dynamic predictive coding: A model of hierarchical sequence learning and prediction in the neocortex

We introduce dynamic predictive coding, a hierarchical model of spatiotemporal prediction and sequence learning in the neocortex. The model assumes that higher cortical levels modulate the temporal dynamics of lower levels, correcting their predictions of dynamics using prediction errors. As a result, lower levels form representations that encode sequences at shorter timescales (e.g., a single step) while higher levels form representations that encode sequences at longer timescales (e.g., an entire sequence). We tested this model using a two-level neural network, where the top-down modulation creates low-dimensional combinations of a set of learned temporal dynamics to explain input sequences. When trained on natural videos, the lower-level model neurons developed space-time receptive fields similar to those of simple cells in the primary visual cortex while the higher-level responses spanned longer timescales, mimicking temporal response hierarchies in the cortex. Additionally, the network's hierarchical sequence representation exhibited both predictive and postdictive effects resembling those observed in visual motion processing in humans (e.g., in the flash-lag illusion). When coupled with an associative memory emulating the role of the hippocampus, the model allowed episodic memories to be stored and retrieved, supporting cue-triggered recall of an input sequence similar to activity recall in the visual cortex. When extended to three hierarchical levels, the model learned progressively more abstract temporal representations along the hierarchy. Taken together, our results suggest that cortical processing and learning of sequences can be interpreted as dynamic predictive coding based on a hierarchical spatiotemporal generative model of the visual world.

Hierarchical rhythmic propagation of corticothalamic interactions for consciousness: A computational study

Clarifying the mechanisms of loss and recovery of consciousness in the brain is a major challenge in neuroscience, and research on the spatiotemporal organization of rhythms at the brain region scale at different levels of consciousness remains scarce. By applying computational neuroscience, an extended corticothalamic network model was developed in this study to simulate the altered states of consciousness induced by different concentration levels of propofol. The cortex area containing oscillation spread from posterior to anterior in four successive time stages, defining four groups of brain regions. A quantitative analysis showed that hierarchical rhythm propagation was mainly due to heterogeneity in the inter-brain region connections. These results indicate that the proposed model is an anatomically data-driven testbed and a simulation platform with millisecond resolution. It facilitates understanding of activity coordination across multiple areas of the conscious brain and the mechanisms of action of anesthetics in terms of brain regions.

The Hopf whole-brain model and its linear approximation

Whole-brain models have proven to be useful to understand the emergence of collective activity among neural populations or brain regions. These models combine connectivity matrices, or connectomes, with local node dynamics, noise, and, eventually, transmission delays. Multiple choices for the local dynamics have been proposed. Among them, nonlinear oscillators corresponding to a supercritical Hopf bifurcation have been used to link brain connectivity and collective phase and amplitude dynamics in different brain states. Here, we studied the linear fluctuations of this model to estimate its stationary statistics, i.e., the instantaneous and lagged covariances and the power spectral densities. This linear approximation-that holds in the case of heterogeneous parameters and time-delays-allows analytical estimation of the statistics and it can be used for fast parameter explorations to study changes in brain state, changes in brain activity due to alterations in structural connectivity, and modulations of parameter due to non-equilibrium dynamics.

The Primate Cortical LFP Exhibits Multiple Spectral and Temporal Gradients and Widespread Task-Dependence During Visual Short-Term Memory

Although cognitive functions are hypothesized to be mediated by synchronous neuronal interactions in multiple frequency bands among widely distributed cortical areas, we still lack a basic understanding of the distribution and task dependence of oscillatory activity across the cortical map. Here, we ask how the spectral and temporal properties of the local field potential (LFP) vary across the primate cerebral cortex, and how they are modulated during visual short-term memory. We measured the LFP from 55 cortical areas in two macaque monkeys while they performed a visual delayed match to sample task. Analysis of peak frequencies in the LFP power spectra reveals multiple discrete frequency bands between 3-80 Hz that differ between the two monkeys. The LFP power in each band, as well as the Sample Entropy, a measure of signal complexity, display distinct spatial gradients across the cortex, some of which correlate with reported spine counts in layer 3 pyramidal neurons. Cortical areas can be robustly decoded using a small number of spectral and temporal parameters, and significant task dependent increases and decreases in spectral power occur in all cortical areas. These findings reveal pronounced, widespread and spatially organized gradients in the spectral and temporal activity of cortical areas. Task-dependent changes in cortical activity are globally distributed, even for a simple cognitive task.

Training-Dependent Gradients of Timescales of Neural Dynamics in the Primate Prefrontal Cortex and Their Contributions to Working Memory

Cortical neurons exhibit multiple timescales related to dynamics of spontaneous fluctuations (intrinsic timescales) and response to task events (seasonal timescales) in addition to selectivity to task-relevant signals. These timescales increase systematically across the cortical hierarchy, for example, from parietal to prefrontal and cingulate cortex, pointing to their role in cortical computations. It is currently unknown whether these timescales are inherent properties of neurons and/or depend on training in a specific task and if the latter, how their modulations contribute to task performance. To address these questions, we analyzed single-cell recordings within five subregions of the prefrontal cortex (PFC) of male macaques before and after training on a working-memory task. We found fine-grained but opposite gradients of intrinsic and seasonal timescales that mainly appeared after training. Intrinsic timescales decreased whereas seasonal timescales increased from posterior to anterior subregions within both dorsal and ventral PFC. Moreover, training was accompanied by increases in proportions of neurons that exhibited intrinsic and seasonal timescales. These effects were comparable to the emergence of response selectivity due to training. Finally, task selectivity accompanied opposite neural dynamics such that neurons with task-relevant selectivity exhibited longer intrinsic and shorter seasonal timescales. Notably, neurons with longer intrinsic and shorter seasonal timescales exhibited superior population-level coding, but these advantages extended to the delay period mainly after training. Together, our results provide evidence for plastic, fine-grained gradients of timescales within PFC that can influence both single-cell and population coding, pointing to the importance of these timescales in understanding cognition.

Characterizing Spatiotemporal Population Receptive Fields in Human Visual Cortex with fMRI

The use of fMRI and computational modeling has advanced understanding of spatial characteristics of population receptive fields (pRFs) in human visual cortex. However, we know relatively little about the spatiotemporal characteristics of pRFs because neurons' temporal properties are one to two orders of magnitude faster than fMRI BOLD responses. Here, we developed an image-computable framework to estimate spatiotemporal pRFs from fMRI data. First, we developed a simulation software that predicts fMRI responses to a time-varying visual input given a spatiotemporal pRF model and solves the model parameters. The simulator revealed that ground-truth spatiotemporal parameters can be accurately recovered at the millisecond resolution from synthesized fMRI responses. Then, using fMRI and a novel stimulus paradigm, we mapped spatiotemporal pRFs in individual voxels across human visual cortex in 10 participants (both females and males). We find that a compressive spatiotemporal (CST) pRF model better explains fMRI responses than a conventional spatial pRF model across visual areas spanning the dorsal, lateral, and ventral streams. Further, we find three organizational principles of spatiotemporal pRFs: (1) from early to later areas within a visual stream, spatial and temporal windows of pRFs progressively increase in size and show greater compressive nonlinearities, (2) later visual areas show diverging spatial and temporal windows across streams, and (3) within early visual areas (V1-V3), both spatial and temporal windows systematically increase with eccentricity. Together, this computational framework and empirical results open exciting new possibilities for modeling and measuring fine-grained spatiotemporal dynamics of neural responses using fMRI.

Computational components of visual predictive coding circuitry

If a full visual percept can be said to be a 'hypothesis', so too can a neural 'prediction' - although the latter addresses one particular component of image content (such as 3-dimensional organisation, the interplay between lighting and surface colour, the future trajectory of moving objects, and so on). And, because processing is hierarchical, predictions generated at one level are conveyed in a backward direction to a lower level, seeking to predict, in fact, the neural activity at that prior stage of processing, and learning from errors signalled in the opposite direction. This is the essence of 'predictive coding', at once an algorithm for information processing and a theoretical basis for the nature of operations performed by the cerebral cortex. Neural models for the implementation of predictive coding invoke specific functional classes of neuron for generating, transmitting and receiving predictions, and for producing reciprocal error signals. Also a third general class, 'precision' neurons, tasked with regulating the magnitude of error signals contingent upon the confidence placed upon the prediction, i.e., the reliability and behavioural utility of the sensory data that it predicts. So, what is the ultimate source of a 'prediction'? The answer is multifactorial: knowledge of the current environmental context and the immediate past, allied to memory and lifetime experience of the way of the world, doubtless fine-tuned by evolutionary history too. There are, in consequence, numerous potential avenues for experimenters seeking to manipulate subjects' expectation, and examine the neural signals elicited by surprising, and less surprising visual stimuli. This review focuses upon the predictive physiology of mouse and monkey visual cortex, summarising and commenting on evidence to date, and placing it in the context of the broader field. It is concluded that predictive coding has a firm grounding in basic neuroscience and that, unsurprisingly, there remains much to learn.

Cell type-specific connectome predicts distributed working memory activity in the mouse brain

Recent advances in connectomics and neurophysiology make it possible to probe whole-brain mechanisms of cognition and behavior. We developed a large-scale model of the multiregional mouse brain for a cardinal cognitive function called working memory, the brain's ability to internally hold and process information without sensory input. The model is built on mesoscopic connectome data for interareal cortical connections and endowed with a macroscopic gradient of measured parvalbumin-expressing interneuron density. We found that working memory coding is distributed yet exhibits modularity; the spatial pattern of mnemonic representation is determined by long-range cell type-specific targeting and density of cell classes. Cell type-specific graph measures predict the activity patterns and a core subnetwork for memory maintenance. The model shows numerous attractor states, which are self-sustained internal states (each engaging a distinct subset of areas). This work provides a framework to interpret large-scale recordings of brain activity during cognition, while highlighting the need for cell type-specific connectomics.

Stability and synchronization in neural network with delayed synaptic connections

In this paper, we investigate the stability of the synchronous state in a complex network using the master stability function technique. We use the extended Hindmarsh-Rose neuronal model including time delayed electrical, chemical, and hybrid couplings. We find the corresponding master stability equation that describes the whole dynamics for each coupling mode. From the maximum Lyapunov exponent, we deduce the stability state for each coupling mode. We observe that for electrical coupling, there exists a mixing between stable and unstable states. For a good setting of some system parameters, the position and the size of unstable areas can be modified. For chemical coupling, we observe difficulties in having a stable area in the complex plane. For hybrid coupling, we observe a stable behavior in the whole system compared to the case where these couplings are considered separately. The obtained results for each coupling mode help to analyze the stability state of some network topologies by using the corresponding eigenvalues. We observe that using electrical coupling can involve a full or partial stability of the system. In the case of chemical coupling, unstable states are observed whereas in the case of hybrid interactions a full stability of the network is obtained. Temporal analysis of the global synchronization is also done for each coupling mode, and the results show that when the network is stable, the synchronization is globally observed, while in the case when it is unstable, its nodes are not globally synchronized.

A reservoir of timescales emerges in recurrent circuits with heterogeneous neural assemblies

The temporal activity of many physical and biological systems, from complex networks to neural circuits, exhibits fluctuations simultaneously varying over a large range of timescales. Long-tailed distributions of intrinsic timescales have been observed across neurons simultaneously recorded within the same cortical circuit. The mechanisms leading to this striking temporal heterogeneity are yet unknown. Here, we show that neural circuits, endowed with heterogeneous neural assemblies of different sizes, naturally generate multiple timescales of activity spanning several orders of magnitude. We develop an analytical theory using rate networks, supported by simulations of spiking networks with cell-type specific connectivity, to explain how neural timescales depend on assembly size and show that our model can naturally explain the long-tailed timescale distribution observed in the awake primate cortex. When driving recurrent networks of heterogeneous neural assemblies by a time-dependent broadband input, we found that large and small assemblies preferentially entrain slow and fast spectral components of the input, respectively. Our results suggest that heterogeneous assemblies can provide a biologically plausible mechanism for neural circuits to demix complex temporal input signals by transforming temporal into spatial neural codes via frequency-selective neural assemblies.

Whole-brain modeling of the differential influences of amyloid-beta and tau in Alzheimer's disease

BACKGROUND

Alzheimer's disease is a neurodegenerative condition associated with the accumulation of two misfolded proteins, amyloid-beta (A[Formula: see text]) and tau. We study their effect on neuronal activity, with the aim of assessing their individual and combined impact.

METHODS

We use a whole-brain dynamic model to find the optimal parameters that best describe the effects of A[Formula: see text] and tau on the excitation-inhibition balance of the local nodes.

RESULTS

We found a clear dominance of A[Formula: see text] over tau in the early disease stages (MCI), while tau dominates over A[Formula: see text] in the latest stages (AD). We identify crucial roles for A[Formula: see text] and tau in complex neuronal dynamics and demonstrate the viability of using regional distributions to define models of large-scale brain function in AD.

CONCLUSIONS

Our study provides further insight into the dynamics and complex interplay between these two proteins, opening the path for further investigations on biomarkers and candidate therapeutic targets in-silico.

The virtual aging brain: Causal inference supports interhemispheric dedifferentiation in healthy aging

The mechanisms of cognitive decline and its variability during healthy aging are not fully understood, but have been associated with reorganization of white matter tracts and functional brain networks. Here, we built a brain network modeling framework to infer the causal link between structural connectivity and functional architecture and the consequent cognitive decline in aging. By applying in-silico interhemispheric degradation of structural connectivity, we reproduced the process of functional dedifferentiation during aging. Thereby, we found the global modulation of brain dynamics by structural connectivity to increase with age, which was steeper in older adults with poor cognitive performance. We validated our causal hypothesis via a deep-learning Bayesian approach. Our results might be the first mechanistic demonstration of dedifferentiation during aging leading to cognitive decline.

Auditory inputs modulate intrinsic neuronal timescales during sleep

Functional magnetic resonance imaging (fMRI) studies have demonstrated that intrinsic neuronal timescales (INT) undergo modulation by external stimulation during consciousness. It remains unclear if INT keep the ability for significant stimulus-induced modulation during primary unconscious states, such as sleep. This fMRI analysis addresses this question via a dataset that comprises an awake resting-state plus rest and stimulus states during sleep. We analyzed INT measured via temporal autocorrelation supported by median frequency (MF) in the frequency-domain. Our results were replicated using a biophysical model. There were two main findings: (1) INT prolonged while MF decreased from the awake resting-state to the N2 resting-state, and (2) INT shortened while MF increased during the auditory stimulus in sleep. The biophysical model supported these results by demonstrating prolonged INT in slowed neuronal populations that simulate the sleep resting-state compared to an awake state. Conversely, under sine wave input simulating the stimulus state during sleep, the model's regions yielded shortened INT that returned to the awake resting-state level. Our results highlight that INT preserve reactivity to stimuli in states of unconsciousness like sleep, enhancing our understanding of unconscious brain dynamics and their reactivity to stimuli.

Decreased intrinsic neural timescales in obsessive compulsive disorder and two distinct subtypes revealed by heterogeneity through discriminative analysis

BACKGROUND

OCD is featured as the destruction of information storage and processing. The cognition of neurobiological and clinical heterogeneity is in suspense and poorly studied.

METHODS

Ninety-nine patients and matched HCs(n = 104) were recruited and underwent resting-state functional MRI scans. We applied INT to evaluate altered local neural dynamics representing the ability of information integration. Moreover, considering OCD was a highly heterogeneous disorder, we investigated putative OCD subtypes from INT using a novel semi-supervised machine learning, named HYDRA.

RESULTS

Compared with HCs, patients with OCD showed decreased INTs in extensive brain regions, including bilateral cerebellum and precuneus, STG/MTG and PCC, hippocampus in DMN; right IFG/MFG/SFG, SPL and bilateral angular gyrus in CEN and insula, SMA in SN. Moreover, many other regions involved in visual processing also had disrupted dynamics of local neural organization, consisting of bilateral CUN, LING and fusiform gyrus and occipital lobe. HYDRA divided patients into two distinct neuroanatomical subtypes from INT. Subtype 1 showed decreased INTs in distributed networks, while subtype 2 presented increased in several common regions which were also found to be decreased in subtype 1, such as STG, IPL, postcentral gyrus and left insula, supramarginal gyrus.

CONCLUSION

This study showed distinct abnormalities from the perspective of dynamics of local neural organization in OCD. Such alteration and dimensional approach may provide a new insight into the prior traditional cognition of this disorder and to some extent do favor of more precise diagnosis and treatment response in the future.

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.

Enhanced behavioral performance through interareal gamma and beta synchronization

Cognitive functioning requires coordination between brain areas. Between visual areas, feedforward gamma synchronization improves behavioral performance. Here, we investigate whether similar principles hold across brain regions and frequency bands, using simultaneous electrocorticographic recordings from 15 areas of two macaque monkeys during performance of a selective attention task. Short behavioral reaction times (RTs), suggesting efficient interareal communication, occurred when occipital areas V1, V2, V4, and DP showed gamma synchronization, and fronto-central areas S1, 5, F1, F2, and F4 showed beta synchronization. For both area clusters and corresponding frequency bands, deviations from the typically observed phase relations increased RTs. Across clusters and frequency bands, good phase relations occurred in a correlated manner specifically when they processed the behaviorally relevant stimulus. Furthermore, the fronto-central cluster exerted a beta-band influence onto the occipital cluster whose strength predicted short RTs. These results suggest that local gamma and beta synchronization and their inter-regional coordination jointly improve behavioral performance.

The alteration of cortical microstructure similarity in drug-resistant epilepsy correlated with mTOR pathway genes

BACKGROUND

Drug-resistant epilepsy (DRE) is associated with distributed laminar disruptions due to cytoarchitectonic pathologies, which may be characterized by multimodal MRI approaches such as morphometric similarity networks (MSNs). However, the genetic and histological underpinning of MSN alterations in DRE remains poorly understood, hampering its clinical application.

METHODS

We enrolled 60 patients with DRE and 23 controls, acquiring T1 and diffusion spectrum imaging data with a 3.0T GE SIGNA Premier scanner. Morphometric similarity networks (MSNs) were constructed and analyzed to identify microstructure similarity differences between patients and controls. Subsequently, patient-specific MSN alteration patterns were associated with gene expression using the GAMBA tool, and layer-specific neuronal signature mapping were also applied. During these analyses, sex and age were adjusted as covariates and multiple comparisons corrections were applied when appropriate.

FINDINGS

We observed widespread MSN changes in patients with DRE and identified five distinct MSN alteration patterns. Major patterns presented pattern-specific associations with expressions of epilepsy-related genes, particularly involving the mTOR pathway. Histological analysis confirmed the presence of cortical microstructure changes in areas with MSN alterations and revealed cellular abnormalities matching the aforementioned genetic risks.

INTERPRETATION

Our findings highlight the potential of quantifying laminar-related microstructure integrity using MSN to uncover the cytoarchitectonic changes in the pathophysiology of DRE. This approach may facilitate the identification of genetic and histological underpinnings of MSN alterations in DRE, ultimately aiding in the development of targeted therapeutic strategies.

FUNDINGS

The National Natural Science Foundation of China, the Ministry of Science and Technology of the People's Republic of China, and the Beijing Municipal Health Commission.

The biological role of local and global fMRI BOLD signal variability in human brain organization

Variability drives the organization and behavior of complex systems, including the human brain. Understanding the variability of brain signals is thus necessary to broaden our window into brain function and behavior. Few empirical investigations of macroscale brain signal variability have yet been undertaken, given the difficulty in separating biological sources of variance from artefactual noise. Here, we characterize the temporal variability of the most predominant macroscale brain signal, the fMRI BOLD signal, and systematically investigate its statistical, topographical and neurobiological properties. We contrast fMRI acquisition protocols, and integrate across histology, microstructure, transcriptomics, neurotransmitter receptor and metabolic data, fMRI static connectivity, and empirical and simulated magnetoencephalography data. We show that BOLD signal variability represents a spatially heterogeneous, central property of multi-scale multi-modal brain organization, distinct from noise. Our work establishes the biological relevance of BOLD signal variability and provides a lens on brain stochasticity across spatial and temporal scales.

Improvements in task performance after practice are associated with scale-free dynamics of brain activity

Although practicing a task generally benefits later performance on that same task, there are individual differences in practice effects. One avenue to model such differences comes from research showing that brain networks extract functional advantages from operating in the vicinity of criticality, a state in which brain network activity is more scale-free. We hypothesized that higher scale-free signal from fMRI data, measured with the Hurst exponent (H), indicates closer proximity to critical states. We tested whether individuals with higher H during repeated task performance would show greater practice effects. In Study 1, participants performed a dual-n-back task (DNB) twice during MRI (n = 56). In Study 2, we used two runs of n-back task (NBK) data from the Human Connectome Project sample (n = 599). In Study 3, participants performed a word completion task (CAST) across six runs (n = 44). In all three studies, multivariate analysis was used to test whether higher H was related to greater practice-related performance improvement. Supporting our hypothesis, we found patterns of higher H that reliably correlated with greater performance improvement across participants in all three studies. However, the predictive brain regions were distinct, suggesting that the specific spatial H↑ patterns are not task-general.

Local minimization of prediction errors drives learning of invariant object representations in a generative network model of visual perception

The ventral visual processing hierarchy of the cortex needs to fulfill at least two key functions: perceived objects must be mapped to high-level representations invariantly of the precise viewing conditions, and a generative model must be learned that allows, for instance, to fill in occluded information guided by visual experience. Here, we show how a multilayered predictive coding network can learn to recognize objects from the bottom up and to generate specific representations via a top-down pathway through a single learning rule: the local minimization of prediction errors. Trained on sequences of continuously transformed objects, neurons in the highest network area become tuned to object identity invariant of precise position, comparable to inferotemporal neurons in macaques. Drawing on this, the dynamic properties of invariant object representations reproduce experimentally observed hierarchies of timescales from low to high levels of the ventral processing stream. The predicted faster decorrelation of error-neuron activity compared to representation neurons is of relevance for the experimental search for neural correlates of prediction errors. Lastly, the generative capacity of the network is confirmed by reconstructing specific object images, robust to partial occlusion of the inputs. By learning invariance from temporal continuity within a generative model, the approach generalizes the predictive coding framework to dynamic inputs in a more biologically plausible way than self-supervised networks with non-local error-backpropagation. This was achieved simply by shifting the training paradigm to dynamic inputs, with little change in architecture and learning rule from static input-reconstructing Hebbian predictive coding networks.

Training-dependent gradients of timescales of neural dynamics in the primate prefrontal cortex and their contributions to working memory

Cortical neurons exhibit multiple timescales related to dynamics of spontaneous fluctuations (intrinsic timescales) and response to task events (seasonal timescales) in addition to selectivity to task-relevant signals. These timescales increase systematically across the cortical hierarchy, e.g., from parietal to prefrontal and cingulate cortex, pointing to their role in cortical computations. It is currently unknown whether these timescales depend on training in a specific task and/or are an inherent property of neurons, and whether more fine-grained hierarchies of timescales exist within specific cortical regions. To address these questions, we analyzed single-cell recordings within five subregions of the prefrontal cortex (PFC) of male macaques before and after training on a working-memory task. We found fine-grained but opposite gradients of intrinsic and seasonal timescales that mainly appeared after training. Intrinsic timescales decreased whereas seasonal timescales increased from posterior to anterior subregions within both dorsal and ventral PFC. Moreover, training was accompanied by increases in proportions of neurons that exhibited intrinsic and seasonal timescales. These effects were comparable to the emergence of response selectivity due to training. Finally, task selectivity accompanied opposite neural dynamics such that neurons with task-relevant selectivity exhibited longer intrinsic and shorter seasonal timescales. Notably, neurons with longer intrinsic and shorter seasonal timescales exhibited superior population-level coding, but these advantages extended to the delay period mainly after training. Together, our results provide evidence for plastic, fine-grained gradients of timescales within PFC that can influence both single-cell and population coding, pointing to the importance of these timescales in understanding cognition.

Significance statement

Recent studies have demonstrated that neural responses exhibit dynamics with different timescales that follow a certain order or hierarchy across cortical areas. While the hierarchy of timescales is consistent across different tasks, it is unknown if these timescales emerge only after training or if they represent inherent properties of neurons. To answer this question, we estimated multiple timescales in neural response across five subregions of the monkeys' lateral prefrontal cortex before and after training on a working-memory task. Our results provide evidence for fine-grained gradients related to certain neural dynamics. Moreover, we show that these timescales depend on and can be modulated by training in a cognitive task, and contribute to encoding of task-relevant information at single-cell and population levels.

Controlling brain dynamics: Landscape and transition path for working memory

Understanding the underlying dynamical mechanisms of the brain and controlling it is a crucial issue in brain science. The energy landscape and transition path approach provides a possible route to address these challenges. Here, taking working memory as an example, we quantified its landscape based on a large-scale macaque model. The working memory function is governed by the change of landscape and brain-wide state switching in response to the task demands. The kinetic transition path reveals that information flow follows the direction of hierarchical structure. Importantly, we propose a landscape control approach to manipulate brain state transition by modulating external stimulation or inter-areal connectivity, demonstrating the crucial roles of associative areas, especially prefrontal and parietal cortical areas in working memory performance. Our findings provide new insights into the dynamical mechanism of cognitive function, and the landscape control approach helps to develop therapeutic strategies for brain disorders.

Single-cell spatial transcriptome reveals cell-type organization in the macaque cortex

Elucidating the cellular organization of the cerebral cortex is critical for understanding brain structure and function. Using large-scale single-nucleus RNA sequencing and spatial transcriptomic analysis of 143 macaque cortical regions, we obtained a comprehensive atlas of 264 transcriptome-defined cortical cell types and mapped their spatial distribution across the entire cortex. We characterized the cortical layer and region preferences of glutamatergic, GABAergic, and non-neuronal cell types, as well as regional differences in cell-type composition and neighborhood complexity. Notably, we discovered a relationship between the regional distribution of various cell types and the region's hierarchical level in the visual and somatosensory systems. Cross-species comparison of transcriptomic data from human, macaque, and mouse cortices further revealed primate-specific cell types that are enriched in layer 4, with their marker genes expressed in a region-dependent manner. Our data provide a cellular and molecular basis for understanding the evolution, development, aging, and pathogenesis of the primate brain.

Cytoarchitectonic, receptor distribution and functional connectivity analyses of the macaque frontal lobe

Based on quantitative cyto- and receptor architectonic analyses, we identified 35 prefrontal areas, including novel subdivisions of Walker's areas 10, 9, 8B, and 46. Statistical analysis of receptor densities revealed regional differences in lateral and ventrolateral prefrontal cortex. Indeed, structural and functional organization of subdivisions encompassing areas 46 and 12 demonstrated significant differences in the interareal levels of α2 receptors. Furthermore, multivariate analysis included receptor fingerprints of previously identified 16 motor areas in the same macaque brains and revealed 5 clusters encompassing frontal lobe areas. We used the MRI datasets from the non-human primate data sharing consortium PRIME-DE to perform functional connectivity analyses using the resulting frontal maps as seed regions. In general, rostrally located frontal areas were characterized by bigger fingerprints, that is, higher receptor densities, and stronger regional interconnections. Whereas more caudal areas had smaller fingerprints, but showed a widespread connectivity pattern with distant cortical regions. Taken together, this study provides a comprehensive insight into the molecular structure underlying the functional organization of the cortex and, thus, reconcile the discrepancies between the structural and functional hierarchical organization of the primate frontal lobe. Finally, our data are publicly available via the EBRAINS and BALSA repositories for the entire scientific community.

Temporospatial Nestedness in Consciousness: An Updated Perspective on the Temporospatial Theory of Consciousness

Time and space are fundamental elements that permeate the fabric of nature, and their significance in relation to neural activity and consciousness remains a compelling yet unexplored area of research. The Temporospatial Theory of Consciousness (TTC) provides a framework that links time, space, neural activity, and consciousness, shedding light on the intricate relationships among these dimensions. In this review, I revisit the fundamental concepts and mechanisms proposed by the TTC, with a particular focus on the central concept of temporospatial nestedness. I propose an extension of temporospatial nestedness by incorporating the nested relationship between the temporal circuit and functional geometry of the brain. To further unravel the complexities of temporospatial nestedness, future research directions should emphasize the characterization of functional geometry and the temporal circuit across multiple spatial and temporal scales. Investigating the links between these scales will yield a more comprehensive understanding of how spatial organization and temporal dynamics contribute to conscious states. This integrative approach holds the potential to uncover novel insights into the neural basis of consciousness and reshape our understanding of the world-brain dynamic.