A hierarchy of intrinsic timescales across primate cortex

Hierarchical consciousness: the Nested Observer Windows model

Foremost in our experience is the intuition that we possess a unified conscious experience. However, many observations run counter to this intuition: we experience paralyzing indecision when faced with two appealing behavioral choices, we simultaneously hold contradictory beliefs, and the content of our thought is often characterized by an internal debate. Here, we propose the Nested Observer Windows (NOW) Model, a framework for hierarchical consciousness wherein information processed across many spatiotemporal scales of the brain feeds into subjective experience. The model likens the mind to a hierarchy of nested mosaic tiles-where an image is composed of mosaic tiles, and each of these tiles is itself an image composed of mosaic tiles. Unitary consciousness exists at the apex of this nested hierarchy where perceptual constructs become fully integrated and complex behaviors are initiated via abstract commands. We define an observer window as a spatially and temporally constrained system within which information is integrated, e.g. in functional brain regions and neurons. Three principles from the signal analysis of electrical activity describe the nested hierarchy and generate testable predictions. First, nested observer windows disseminate information across spatiotemporal scales with cross-frequency coupling. Second, observer windows are characterized by a high degree of internal synchrony (with zero phase lag). Third, observer windows at the same spatiotemporal level share information with each other through coherence (with non-zero phase lag). The theoretical framework of the NOW Model accounts for a wide range of subjective experiences and a novel approach for integrating prominent theories of consciousness.

Primate thalamic nuclei select abstract rules and shape prefrontal dynamics

Flexible behavior depends on abstract rules to generalize beyond specific instances, and outcome monitoring to adjust actions. Cortical circuits are posited to read out rules from high-dimensional representations of task-relevant variables in prefrontal cortex (PFC). We instead hypothesized that converging inputs from PFC, directly or via basal ganglia (BG), enable primate-specific thalamus to select rules. To test this, we simultaneously measured spiking activity across PFC and two connected thalamic nuclei of monkeys applying rules. Abstract rule information first appeared in the ventroanterior thalamus (VA) - the main thalamic hub between BG and PFC. The mediodorsal thalamus (MD) also represented rule information before PFC, which persisted after rule cues were removed, to help maintain activation of relevant posterior PFC cell ensembles. MD, a major recipient of midbrain dopamine input, was first to represent information about behavioral outcomes. This persisted after the trial (also in PFC). A PFC-BG-thalamus model reproduced key findings, and thalamic-lesion modeling disrupted PFC rule representations. These results suggest a revised view of the neural basis of flexible behavior in primates, featuring a central role for thalamus in selecting high-level cognitive information from PFC and implementing post-error behavioral adjustments, and of the functional organization of PFC along its anterior-posterior dimension.

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.

Different extrapolation of moving object locations in perception, smooth pursuit, and saccades

The ability to accurately perceive and track moving objects is crucial for many everyday activities. In this study, we use a "double-drift stimulus" to explore the processing of visual motion signals that underlie perception, pursuit, and saccade responses to a moving object. Participants were presented with peripheral moving apertures filled with noise that either drifted orthogonally to the aperture's direction or had no net motion. Participants were asked to saccade to and track these targets with their gaze as soon as they appeared and then to report their direction. In the trials with internal motion, the target disappeared at saccade onset so that the first 100 ms of the postsaccadic pursuit response was driven uniquely by peripheral information gathered before saccade onset. This provided independent measures of perceptual, pursuit, and saccadic responses to the double-drift stimulus on a trial-by-trial basis. Our analysis revealed systematic differences between saccadic responses, on one hand, and perceptual and pursuit responses, on the other. These differences are unlikely to be caused by differences in the processing of motion signals because both saccades and pursuits seem to rely on shared target position and velocity information. We conclude that our results are instead due to a difference in how the processing mechanisms underlying perception, pursuit, and saccades combine motor signals with target position. These findings advance our understanding of the mechanisms underlying dissociation in visual processing between perception and eye movements.

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

BACKGROUND

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

PURPOSE

To assess the intrinsic timescale and GMV in PD.

STUDY TYPE

Prospective.

POPULATION

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

FIELD STRENGTH/SEQUENCE

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

ASSESSMENT

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

STATISTICAL TEST

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

RESULTS

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

CONCLUSION

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

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 1.

Bayesian inference is facilitated by modular neural networks with different time scales

Various animals, including humans, have been suggested to perform Bayesian inferences to handle noisy, time-varying external information. In performing Bayesian inference by the brain, the prior distribution must be acquired and represented by sampling noisy external inputs. However, the mechanism by which neural activities represent such distributions has not yet been elucidated. Our findings reveal that networks with modular structures, composed of fast and slow modules, are adept at representing this prior distribution, enabling more accurate Bayesian inferences. Specifically, the modular network that consists of a main module connected with input and output layers and a sub-module with slower neural activity connected only with the main module outperformed networks with uniform time scales. Prior information was represented specifically by the slow sub-module, which could integrate observed signals over an appropriate period and represent input means and variances. Accordingly, the neural network could effectively predict the time-varying inputs. Furthermore, by training the time scales of neurons starting from networks with uniform time scales and without modular structure, the above slow-fast modular network structure and the division of roles in which prior knowledge is selectively represented in the slow sub-modules spontaneously emerged. These results explain how the prior distribution for Bayesian inference is represented in the brain, provide insight into the relevance of modular structure with time scale hierarchy to information processing, and elucidate the significance of brain areas with slower time scales.

Curiosity: primate neural circuits for novelty and information seeking

For many years, neuroscientists have investigated the behavioural, computational and neurobiological mechanisms that support value-based decisions, revealing how humans and animals make choices to obtain rewards. However, many decisions are influenced by factors other than the value of physical rewards or second-order reinforcers (such as money). For instance, animals (including humans) frequently explore novel objects that have no intrinsic value solely because they are novel and they exhibit the desire to gain information to reduce their uncertainties about the future, even if this information cannot lead to reward or assist them in accomplishing upcoming tasks. In this Review, I discuss how circuits in the primate brain responsible for detecting, predicting and assessing novelty and uncertainty regulate behaviour and give rise to these behavioural components of curiosity. I also briefly discuss how curiosity-related behaviours arise during postnatal development and point out some important reasons for the persistence of curiosity across generations.

Dynamic fading memory and expectancy effects in the monkey primary visual cortex

In order to investigate the involvement of the primary visual cortex (V1) in working memory (WM), parallel, multisite recordings of multi-unit activity were obtained from monkey V1 while the animals performed a delayed match-to-sample (DMS) task. During the delay period, V1 population firing rate vectors maintained a lingering trace of the sample stimulus that could be reactivated by intervening impulse stimuli that enhanced neuronal firing. This fading trace of the sample did not require active engagement of the monkeys in the DMS task and likely reflects the intrinsic dynamics of recurrent cortical networks in lower visual areas. This renders an active, attention-dependent involvement of V1 in the maintenance of WM contents unlikely. By contrast, population responses to the test stimulus depended on the probabilistic contingencies between sample and test stimuli. Responses to tests that matched expectations were reduced which agrees with concepts of predictive coding.

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.

Optimal decoding of neural dynamics occurs at mesoscale spatial and temporal resolutions

Introduction

Understanding the neural code has been one of the central aims of neuroscience research for decades. Spikes are commonly referred to as the units of information transfer, but multi-unit activity (MUA) recordings are routinely analyzed in aggregate forms such as binned spike counts, peri-stimulus time histograms, firing rates, or population codes. Various forms of averaging also occur in the brain, from the spatial averaging of spikes within dendritic trees to their temporal averaging through synaptic dynamics. However, how these forms of averaging are related to each other or to the spatial and temporal units of information representation within the neural code has remained poorly understood.

Materials and methods

In this work we developed NeuroPixelHD, a symbolic hyperdimensional model of MUA, and used it to decode the spatial location and identity of static images shown to n = 9 mice in the Allen Institute Visual Coding-NeuroPixels dataset from large-scale MUA recordings. We parametrically varied the spatial and temporal resolutions of the MUA data provided to the model, and compared its resulting decoding accuracy.

Results

For almost all subjects, we found 125ms temporal resolution to maximize decoding accuracy for both the spatial location of Gabor patches (81 classes for patches presented over a 9×9 grid) as well as the identity of natural images (118 classes corresponding to 118 images) across the whole brain. This optimal temporal resolution nevertheless varied greatly between different regions, followed a sensory-associate hierarchy, and was significantly modulated by the central frequency of theta-band oscillations across different regions. Spatially, the optimal resolution was at either of two mesoscale levels for almost all mice: the area level, where the spiking activity of all neurons within each brain area are combined, and the population level, where neuronal spikes within each area are combined across fast spiking (putatively inhibitory) and regular spiking (putatively excitatory) neurons, respectively. We also observed an expected interplay between optimal spatial and temporal resolutions, whereby increasing the amount of averaging across one dimension (space or time) decreases the amount of averaging that is optimal across the other dimension, and vice versa.

Discussion

Our findings corroborate existing empirical practices of spatiotemporal binning and averaging in MUA data analysis, and provide a rigorous computational framework for optimizing the level of such aggregations. Our findings can also synthesize these empirical practices with existing knowledge of the various sources of biological averaging in the brain into a new theory of neural information processing in which the unit of information varies dynamically based on neuronal signal and noise correlations across space and time.

Unsupervised learning of mid-level visual representations

Recently, a confluence between trends in neuroscience and machine learning has brought a renewed focus on unsupervised learning, where sensory processing systems learn to exploit the statistical structure of their inputs in the absence of explicit training targets or rewards. Sophisticated experimental approaches have enabled the investigation of the influence of sensory experience on neural self-organization and its synaptic bases. Meanwhile, novel algorithms for unsupervised and self-supervised learning have become increasingly popular both as inspiration for theories of the brain, particularly for the function of intermediate visual cortical areas, and as building blocks of real-world learning machines. Here we review some of these recent developments, placing them in historical context and highlighting some research lines that promise exciting breakthroughs in the near future.

Neural tracking of continuous acoustics: properties, speech-specificity and open questions

Human speech is a particularly relevant acoustic stimulus for our species, due to its role of information transmission during communication. Speech is inherently a dynamic signal, and a recent line of research focused on neural activity following the temporal structure of speech. We review findings that characterise neural dynamics in the processing of continuous acoustics and that allow us to compare these dynamics with temporal aspects in human speech. We highlight properties and constraints that both neural and speech dynamics have, suggesting that auditory neural systems are optimised to process human speech. We then discuss the speech-specificity of neural dynamics and their potential mechanistic origins and summarise open questions in the field.

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.

Intrinsic timescales and predictive allostatic interoception in brain health and disease

The cognitive neuroscience of brain diseases faces challenges in understanding the complex relationship between brain structure and function, the heterogeneity of brain phenotypes, and the lack of dimensional and transnosological explanations. This perspective offers a framework combining the predictive coding theory of allostatic interoceptive overload (PAIO) and the intrinsic neural timescales (INT) theory to provide a more dynamic understanding of brain health in psychiatry and neurology. PAIO integrates allostasis and interoception to assess the interaction between internal patterns and environmental stressors, while INT shows that different brain regions operate on different intrinsic timescales. The allostatic overload can be understood as a failure of INT, which involves a breakdown of proper temporal integration and segregation. This can lead to dimensional disbalances between exteroceptive/interoceptive inputs across brain and whole-body levels (cardiometabolic, cardiovascular, inflammatory, immune). This approach offers new insights, presenting novel perspectives on brain spatiotemporal hierarchies and interactions. By integrating these theories, the paper opens innovative paths for studying brain health dynamics, which can inform future research in brain health and disease.

Individual variability in functional connectivity of human auditory cortex

Individual variability in functional connectivity underlies individual differences in cognition and behaviors, yet its association with functional specialization in the auditory cortex remains elusive. Using resting-state functional magnetic resonance imaging data from the Human Connectome Project, this study was designed to investigate the spatial distribution of auditory cortex individual variability in its whole-brain functional network architecture. An inherent hierarchical axis of the variability was discerned, which radiates from the medial to lateral orientation, with the left auditory cortex demonstrating more pronounced variations than the right. This variability exhibited a significant correlation with the variations in structural and functional metrics in the auditory cortex. Four auditory cortex subregions, which were identified from a clustering analysis based on this variability, exhibited unique connectional fingerprints and cognitive maps, with certain subregions showing specificity to speech perception functional activation. Moreover, the lateralization of the connectional fingerprint exhibited a U-shaped trajectory across the subregions. These findings emphasize the role of individual variability in functional connectivity in understanding cortical functional organization, as well as in revealing its association with functional specialization from the activation, connectome, and cognition perspectives.

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.

Inferring neural dynamics of memory during naturalistic social communication

Memory processes in complex behaviors like social communication require forming representations of the past that grow with time. The neural mechanisms that support such continually growing memory remain unknown. We address this gap in the context of fly courtship, a natural social behavior involving the production and perception of long, complex song sequences. To study female memory for male song history in unrestrained courtship, we present 'Natural Continuation' (NC)-a general, simulation-based model comparison procedure to evaluate candidate neural codes for complex stimuli using naturalistic behavioral data. Applying NC to fly courtship revealed strong evidence for an adaptive population mechanism for how female auditory neural dynamics could convert long song histories into a rich mnemonic format. Song temporal patterning is continually transformed by heterogeneous nonlinear adaptation dynamics, then integrated into persistent activity, enabling common neural mechanisms to retain continuously unfolding information over long periods and yielding state-of-the-art predictions of female courtship behavior. At a population level this coding model produces multi-dimensional advection-diffusion-like responses that separate songs over a continuum of timescales and can be linearly transformed into flexible output signals, illustrating its potential to create a generic, scalable mnemonic format for extended input signals poised to drive complex behavioral responses. This work thus shows how naturalistic behavior can directly inform neural population coding models, revealing here a novel process for memory formation.

Mixed Selectivity Coding of Content-Temporal Detail by Dorsomedial Posterior Parietal Neurons

The dorsomedial posterior parietal cortex (dmPPC) is part of a higher-cognition network implicated in elaborate processes underpinning memory formation, recollection, episode reconstruction, and temporal information processing. Neural coding for complex episodic processing is however under-documented. Here, we recorded extracellular neural activities from three male rhesus macaques (Macaca mulatta) and revealed a set of neural codes of "neuroethogram" in the primate parietal cortex. Analyzing neural responses in macaque dmPPC to naturalistic videos, we discovered several groups of neurons that are sensitive to different categories of ethogram items, low-level sensory features, and saccadic eye movement. We also discovered that the processing of category and feature information by these neurons is sustained by the accumulation of temporal information over a long timescale of up to 30 s, corroborating its reported long temporal receptive windows. We performed an additional behavioral experiment with additional two male rhesus macaques and found that saccade-related activities could not account for the mixed neuronal responses elicited by the video stimuli. We further observed monkeys' scan paths and gaze consistency are modulated by video content. Taken altogether, these neural findings explain how dmPPC weaves fabrics of ongoing experiences together in real time. The high dimensionality of neural representations should motivate us to shift the focus of attention from pure selectivity neurons to mixed selectivity neurons, especially in increasingly complex naturalistic task designs.

A week in the life of the human brain: stable states punctuated by chaotic transitions

Many important neurocognitive states, such as performing natural activities and fluctuations of arousal, shift over minutes-to-hours in the real-world. We harnessed 3-12 days of continuous multi-electrode intracranial recordings in twenty humans during natural behavior (socializing, using digital devices, sleeping, etc.) to study real-world neurodynamics. Applying deep learning with dynamical systems approaches revealed that brain networks formed consistent stable states that predicted behavior and physiology. Changes in behavior were associated with bursts of rapid neural fluctuations where brain networks chaotically explored many configurations before settling into new states. These trajectories traversed an hourglass-shaped structure anchored around a set of networks that slowly tracked levels of outward awareness related to wake-sleep stages, and a central attractor corresponding to default mode network activation. These findings indicate ways our brains use rapid, chaotic transitions that coalesce into neurocognitive states slowly fluctuating around a stabilizing central equilibrium to balance flexibility and stability during real-world behavior.

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.

A role for the serotonin 2A receptor in the expansion and functioning of human transmodal cortex

Integrating independent but converging lines of research on brain function and neurodevelopment across scales, this article proposes that serotonin 2A receptor (5-HT2AR) signalling is an evolutionary and developmental driver and potent modulator of the macroscale functional organization of the human cerebral cortex. A wealth of evidence indicates that the anatomical and functional organization of the cortex follows a unimodal-to-transmodal gradient. Situated at the apex of this processing hierarchy-where it plays a central role in the integrative processes underpinning complex, human-defining cognition-the transmodal cortex has disproportionately expanded across human development and evolution. Notably, the adult human transmodal cortex is especially rich in 5-HT2AR expression and recent evidence suggests that, during early brain development, 5-HT2AR signalling on neural progenitor cells stimulates their proliferation-a critical process for evolutionarily-relevant cortical expansion. Drawing on multimodal neuroimaging and cross-species investigations, we argue that, by contributing to the expansion of the human cortex and being prevalent at the apex of its hierarchy in the adult brain, 5-HT2AR signalling plays a major role in both human cortical expansion and functioning. Owing to its unique excitatory and downstream cellular effects, neuronal 5-HT2AR agonism promotes neuroplasticity, learning and cognitive and psychological flexibility in a context-(hyper)sensitive manner with therapeutic potential. Overall, we delineate a dual role of 5-HT2ARs in enabling both the expansion and modulation of the human transmodal cortex.

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.

Protosequences in human cortical organoids model intrinsic states in the developing cortex

Neuronal firing sequences are thought to be the basic building blocks of neural coding and information broadcasting within the brain. However, when sequences emerge during neurodevelopment remains unknown. We demonstrate that structured firing sequences are present in spontaneous activity of human brain organoids and ex vivo neonatal brain slices from the murine somatosensory cortex. We observed a balance between temporally rigid and flexible firing patterns that are emergent phenomena in human brain organoids and early postnatal murine somatosensory cortex, but not in primary dissociated cortical cultures. Our findings suggest that temporal sequences do not arise in an experience-dependent manner, but are rather constrained by an innate preconfigured architecture established during neurogenesis. These findings highlight the potential for brain organoids to further explore how exogenous inputs can be used to refine neuronal circuits and enable new studies into the genetic mechanisms that govern assembly of functional circuitry during early human brain development.

The neuronal implementation of representational geometry in primate prefrontal cortex

Modern neuroscience has seen the rise of a population-doctrine that represents cognitive variables using geometrical structures in activity space. Representational geometry does not, however, account for how individual neurons implement these representations. Leveraging the principle of sparse coding, we present a framework to dissect representational geometry into biologically interpretable components that retain links to single neurons. Applied to extracellular recordings from the primate prefrontal cortex in a working memory task with interference, the identified components revealed disentangled and sequential memory representations including the recovery of memory content after distraction, signals hidden to conventional analyses. Each component was contributed by small subpopulations of neurons with distinct spiking properties and response dynamics. Modeling showed that such sparse implementations are supported by recurrently connected circuits as in prefrontal cortex. The perspective of neuronal implementation links representational geometries to their cellular constituents, providing mechanistic insights into how neural systems encode and process information.

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.

Evidence for increased parallel information transmission in human brain networks compared to macaques and male mice

Brain communication, defined as information transmission through white-matter connections, is at the foundation of the brain's computational capacities that subtend almost all aspects of behavior: from sensory perception shared across mammalian species, to complex cognitive functions in humans. How did communication strategies in macroscale brain networks adapt across evolution to accomplish increasingly complex functions? By applying a graph- and information-theory approach to assess information-related pathways in male mouse, macaque and human brains, we show a brain communication gap between selective information transmission in non-human mammals, where brain regions share information through single polysynaptic pathways, and parallel information transmission in humans, where regions share information through multiple parallel pathways. In humans, parallel transmission acts as a major connector between unimodal and transmodal systems. The layout of information-related pathways is unique to individuals across different mammalian species, pointing at the individual-level specificity of information routing architecture. Our work provides evidence that different communication patterns are tied to the evolution of mammalian brain networks.

Complexity and 1/f slope jointly reflect brain states

Characterization of brain states is essential for understanding its functioning in the absence of external stimuli. Brain states differ on their balance between excitation and inhibition, and on the diversity of their activity patterns. These can be respectively indexed by 1/f slope and Lempel-Ziv complexity (LZc). However, whether and how these two brain state properties relate remain elusive. Here we analyzed the relation between 1/f slope and LZc with two in-silico approaches and in both rat EEG and monkey ECoG data. We contrasted resting state with propofol anesthesia, which directly modulates the excitation-inhibition balance. We found convergent results among simulated and empirical data, showing a strong, inverse and non trivial monotonic relation between 1/f slope and complexity, consistent at both ECoG and EEG scales. We hypothesize that differentially entropic regimes could underlie the link between the excitation-inhibition balance and the vastness of the repertoire of brain systems.

Towards a biologically annotated brain connectome

The brain is a network of interleaved neural circuits. In modern connectomics, brain connectivity is typically encoded as a network of nodes and edges, abstracting away the rich biological detail of local neuronal populations. Yet biological annotations for network nodes - such as gene expression, cytoarchitecture, neurotransmitter receptors or intrinsic dynamics - can be readily measured and overlaid on network models. Here we review how connectomes can be represented and analysed as annotated networks. Annotated connectomes allow us to reconceptualize architectural features of networks and to relate the connection patterns of brain regions to their underlying biology. Emerging work demonstrates that annotated connectomes help to make more veridical models of brain network formation, neural dynamics and disease propagation. Finally, annotations can be used to infer entirely new inter-regional relationships and to construct new types of network that complement existing connectome representations. In summary, biologically annotated connectomes offer a compelling way to study neural wiring in concert with local biological features.

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.

Neuropeptide Signaling is Required to Implement a Line Attractor Encoding a Persistent Internal Behavioral State

Internal states drive survival behaviors, but their neural implementation is not well understood. Recently we identified a line attractor in the ventromedial hypothalamus (VMH) that represents an internal state of aggressiveness. Line attractors can be implemented by recurrent connectivity and/or neuromodulatory signaling, but evidence for the latter is scant. Here we show that neuropeptidergic signaling is necessary for line attractor dynamics in this system, using a novel approach that integrates cell type-specific, anatomically restricted CRISPR/Cas9-based gene editing with microendoscopic calcium imaging. Co-disruption of receptors for oxytocin and vasopressin in adult VMH Esr1 + neurons that control aggression suppressed attack, reduced persistent neural activity and eliminated line attractor dynamics, while only modestly impacting neural activity and sex- or behavior-tuning. These data identify a requisite role for neuropeptidergic signaling in implementing a behaviorally relevant line attractor. Our approach should facilitate mechanistic studies in neuroscience that bridge different levels of biological function and abstraction.

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.

Hierarchical network structure as the source of hierarchical dynamics (power-law frequency spectra) in living and non-living systems: How state-trait continua (body plans, personalities) emerge from first principles in biophysics

Living systems are hierarchical control systems that display a small world network structure. In such structures, many smaller clusters are nested within fewer larger ones, producing a fractal-like structure with a 'power-law' cluster size distribution (a mereology). Just like their structure, the dynamics of living systems shows fractal-like qualities: the timeseries of inner message passing and overt behavior contain high frequencies or 'states' (treble) that are nested within lower frequencies or 'traits' (bass), producing a power-law frequency spectrum that is known as a 'state-trait continuum' in the behavioral sciences. Here, we argue that the power-law dynamics of living systems results from their power-law network structure: organisms 'vertically encode' the deep spatiotemporal structure of their (anticipated) environments, to the effect that many small clusters near the base of the hierarchy produce high frequency signal changes and fewer larger clusters at its top produce ultra-low frequencies. Such ultra-low frequencies exert a tonic regulatory pressure that produces morphological as well as behavioral traits (i.e., body plans and personalities). Nested-modular structure causes higher frequencies to be embedded within lower frequencies, producing a power-law state-trait continuum. At the heart of such dynamics lies the need for efficient energy dissipation through networks of coupled oscillators, which also governs the dynamics of non-living systems (e.q., earthquakes, stock market fluctuations). Since hierarchical structure produces hierarchical dynamics, the development and collapse of hierarchical structure (e.g., during maturation and disease) should leave specific traces in system dynamics (shifts in lower frequencies, i.e. morphological and behavioral traits) that may serve as early warning signs to system failure. The applications of this idea range from (bio)physics and phylogenesis to ontogenesis and clinical medicine.

Architecture of the brain's visual system enhances network stability and performance through layers, delays, and feedback

In the visual system of primates, image information propagates across successive cortical areas, and there is also local feedback within an area and long-range feedback across areas. Recent findings suggest that the resulting temporal dynamics of neural activity are crucial in several vision tasks. In contrast, artificial neural network models of vision are typically feedforward and do not capitalize on the benefits of temporal dynamics, partly due to concerns about stability and computational costs. In this study, we focus on recurrent networks with feedback connections for visual tasks with static input corresponding to a single fixation. We demonstrate mathematically that a network's dynamics can be stabilized by four key features of biological networks: layer-ordered structure, temporal delays between layers, longer distance feedback across layers, and nonlinear neuronal responses. Conversely, when feedback has a fixed distance, one can omit delays in feedforward connections to achieve more efficient artificial implementations. We also evaluated the effect of feedback connections on object detection and classification performance using standard benchmarks, specifically the COCO and CIFAR10 datasets. Our findings indicate that feedback connections improved the detection of small objects, and classification performance became more robust to noise. We found that performance increased with the temporal dynamics, not unlike what is observed in core vision of primates. These results suggest that delays and layered organization are crucial features for stability and performance in both biological and artificial recurrent neural 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.

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.

Simple and complex cells revisited: toward a selectivity-invariance model of object recognition

This paper presents a theoretical perspective on modeling ventral stream processing by revisiting the computational abstraction of simple and complex cells. In parallel to David Marr's vision theory, we organize the new perspective into three levels. At the computational level, we abstract simple and complex cells into space partitioning and composition in a topological space based on the redundancy exploitation hypothesis of Horace Barlow. At the algorithmic level, we present a hierarchical extension of sparse coding by exploiting the manifold constraint in high-dimensional space (i.e., the blessing of dimensionality). The resulting over-parameterized models for object recognition differ from existing hierarchical models by disentangling the objectives of selectivity and invariance computation. It is possible to interpret our hierarchical construction as a computational implementation of cortically local subspace untangling for object recognition and face representation, which are closely related to exemplar-based and axis-based coding in the medial temporal lobe. At the implementation level, we briefly discuss two possible implementations based on asymmetric sparse autoencoders and divergent spiking neural networks.

Evaluation of the Hierarchical Correspondence between the Human Brain and Artificial Neural Networks: A Review

Artificial neural networks (ANNs) that are heavily inspired by the human brain now achieve human-level performance across multiple task domains. ANNs have thus drawn attention in neuroscience, raising the possibility of providing a framework for understanding the information encoded in the human brain. However, the correspondence between ANNs and the brain cannot be measured directly. They differ in outputs and substrates, neurons vastly outnumber their ANN analogs (i.e., nodes), and the key algorithm responsible for most of modern ANN training (i.e., backpropagation) is likely absent from the brain. Neuroscientists have thus taken a variety of approaches to examine the similarity between the brain and ANNs at multiple levels of their information hierarchy. This review provides an overview of the currently available approaches and their limitations for evaluating brain-ANN correspondence.

Neuromodulation of striatal D1 cells shapes BOLD fluctuations in anatomically connected thalamic and cortical regions

Understanding how the brain's macroscale dynamics are shaped by underlying microscale mechanisms is a key problem in neuroscience. In animal models, we can now investigate this relationship in unprecedented detail by directly manipulating cellular-level properties while measuring the whole-brain response using resting-state fMRI. Here, we focused on understanding how blood-oxygen-level-dependent (BOLD) dynamics, measured within a structurally well-defined striato-thalamo-cortical circuit in mice, are shaped by chemogenetically exciting or inhibiting D1 medium spiny neurons (MSNs) of the right dorsomedial caudate putamen (CPdm). We characterize changes in both the BOLD dynamics of individual cortical and subcortical brain areas, and patterns of inter-regional coupling (functional connectivity) between pairs of areas. Using a classification approach based on a large and diverse set of time-series properties, we found that CPdm neuromodulation alters BOLD dynamics within thalamic subregions that project back to dorsomedial striatum. In the cortex, changes in local dynamics were strongest in unimodal regions (which process information from a single sensory modality) and weakened along a hierarchical gradient towards transmodal regions. In contrast, a decrease in functional connectivity was observed only for cortico-striatal connections after D1 excitation. Our results show that targeted cellular-level manipulations affect local BOLD dynamics at the macroscale, such as by making BOLD dynamics more predictable over time by increasing its self-correlation structure. This contributes to ongoing attempts to understand the influence of structure-function relationships in shaping inter-regional communication at subcortical and cortical levels.

A phylogenetically-conserved axis of thalamocortical connectivity in the human brain

The thalamus enables key sensory, motor, emotive, and cognitive processes via connections to the cortex. These projection patterns are traditionally considered to originate from discrete thalamic nuclei, however recent work showing gradients of molecular and connectivity features in the thalamus suggests the organisation of thalamocortical connections occurs along a continuous dimension. By performing a joint decomposition of densely sampled gene expression and non-invasive diffusion tractography in the adult human thalamus, we define a principal axis of genetic and connectomic variation along a medial-lateral thalamic gradient. Projections along this axis correspond to an anterior-posterior cortical pattern and are aligned with electrophysiological properties of the cortex. The medial-lateral axis demonstrates phylogenetic conservation, reflects transitions in neuronal subtypes, and shows associations with neurodevelopment and common brain disorders. This study provides evidence for a supra-nuclear axis of thalamocortical organisation characterised by a graded transition in molecular properties and anatomical connectivity.

Neurophysiological signatures of cortical micro-architecture

Systematic spatial variation in micro-architecture is observed across the cortex. These micro-architectural gradients are reflected in neural activity, which can be captured by neurophysiological time-series. How spontaneous neurophysiological dynamics are organized across the cortex and how they arise from heterogeneous cortical micro-architecture remains unknown. Here we extensively profile regional neurophysiological dynamics across the human brain by estimating over 6800 time-series features from the resting state magnetoencephalography (MEG) signal. We then map regional time-series profiles to a comprehensive multi-modal, multi-scale atlas of cortical micro-architecture, including microstructure, metabolism, neurotransmitter receptors, cell types and laminar differentiation. We find that the dominant axis of neurophysiological dynamics reflects characteristics of power spectrum density and linear correlation structure of the signal, emphasizing the importance of conventional features of electromagnetic dynamics while identifying additional informative features that have traditionally received less attention. Moreover, spatial variation in neurophysiological dynamics is co-localized with multiple micro-architectural features, including gene expression gradients, intracortical myelin, neurotransmitter receptors and transporters, and oxygen and glucose metabolism. Collectively, this work opens new avenues for studying the anatomical basis of neural activity.

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.

Automated customization of large-scale spiking network models to neuronal population activity

Understanding brain function is facilitated by constructing computational models that accurately reproduce aspects of brain activity. Networks of spiking neurons capture the underlying biophysics of neuronal circuits, yet the dependence of their activity on model parameters is notoriously complex. As a result, heuristic methods have been used to configure spiking network models, which can lead to an inability to discover activity regimes complex enough to match large-scale neuronal recordings. Here we propose an automatic procedure, Spiking Network Optimization using Population Statistics (SNOPS), to customize spiking network models that reproduce the population-wide covariability of large-scale neuronal recordings. We first confirmed that SNOPS accurately recovers simulated neural activity statistics. Then, we applied SNOPS to recordings in macaque visual and prefrontal cortices and discovered previously unknown limitations of spiking network models. Taken together, SNOPS can guide the development of network models and thereby enable deeper insight into how networks of neurons give rise to brain function.

Routing states transition during oscillatory bursts and attentional selection

Brain-wide information routing relies on the spatio-temporal dynamics of neural activity, but it remains unclear how routing states emerge at fast spiking timescales and relate to slower activity dynamics during cognitive processes. Here, we show that localized spiking events participate in directional routing states with spiking activity in distant brain areas that dynamically switch or amplify states during oscillatory bursts, attentional selection, and decision-making. Modeling and neural recordings from lateral prefrontal cortex (LPFC), anterior cingulate cortex (ACC), and striatum of nonhuman primates revealed that cross-regional routing states arise within 20 ms following individual neuron spikes, with LPFC spikes leading the activity in ACC and striatum. The baseline routing state amplified during LPFC beta bursts between LPFC and striatum and switched direction during ACC theta/alpha bursts between ACC and LPFC. Selective attention amplified theta-/alpha-band-specific lead ensembles in ACC, while decision-making increased the lead of ACC and LPFC spikes to the striatum.

The neural basis of resting-state fMRI functional connectivity in fronto-limbic circuits revealed by chemogenetic manipulation

Measures of fMRI resting-state functional connectivity (rs-FC) are an essential tool for basic and clinical investigations of fronto-limbic circuits. Understanding the relationship between rs-FC and neural activity in these circuits is therefore vital. Here we introduced inhibitory designer receptors exclusively activated by designer drugs (DREADDs) into the macaque amygdala and activated them with a highly selective and potent DREADD agonist, deschloroclozapine. We evaluated the causal effect of activating the DREADD receptors on rs-FC and neural activity within circuits connecting amygdala and frontal cortex. Interestingly, activating the inhibitory DREADD increased rs-FC between amygdala and ventrolateral prefrontal cortex. Neurophysiological recordings revealed that the DREADD-induced increase in fMRI rs-FC was associated with increased local field potential coherency in the alpha band (6.5-14.5Hz) between amygdala and ventrolateral prefrontal cortex. Thus, our multi-disciplinary approach reveals the specific signature of neuronal activity that underlies rs-FC in fronto-limbic circuits.

Cerebellar contributions across behavioural timescales: a review from the perspective of cerebro-cerebellar interactions

Performing successful adaptive behaviour relies on our ability to process a wide range of temporal intervals with certain precision. Studies on the role of the cerebellum in temporal information processing have adopted the dogma that the cerebellum is involved in sub-second processing. However, emerging evidence shows that the cerebellum might be involved in suprasecond temporal processing as well. Here we review the reciprocal loops between cerebellum and cerebral cortex and provide a theoretical account of cerebro-cerebellar interactions with a focus on how cerebellar output can modulate cerebral processing during learning of complex sequences. Finally, we propose that while the ability of the cerebellum to support millisecond timescales might be intrinsic to cerebellar circuitry, the ability to support supra-second timescales might result from cerebellar interactions with other brain regions, such as the prefrontal cortex.