A thalamocortical substrate for integrated information via critical synchronous bursting.

Understanding the neurobiological mechanisms underlying consciousness remains a significant challenge. Recent evidence suggests that the coupling between distal-apical and basal-somatic dendrites in thick-tufted layer 5 pyramidal neurons (L5PN), regulated by the nonspecific-projecting thalamus, is crucial for consciousness. Yet, it is uncertain whether this thalamocortical mechanism can support emergent signatures of consciousness, such as integrated information. To address this question, we constructed a biophysical network of dual-compartment thick-tufted L5PN, with dendrosomatic coupling controlled by thalamic inputs. Our findings demonstrate that integrated information is maximized when nonspecific thalamic inputs drive the system into a regime of time-varying synchronous bursting. Here, the system exhibits variable spiking dynamics with broad pairwise correlations, supporting the enhanced integrated information. Further, the observed peak in integrated information aligns with criticality signatures and empirically observed layer 5 pyramidal bursting rates. These results suggest that the thalamocortical core of the mammalian brain may be evolutionarily configured to optimize effective information processing, providing a potential neuronal mechanism that integrates microscale theories with macroscale signatures of consciousness.

Abnormal higher-order network interactions in Parkinson's disease visual hallucinations.

Visual hallucinations in Parkinson's disease can be viewed from a systems-level perspective, whereby dysfunctional communication between brain networks responsible for perception predisposes a person to hallucinate. To this end, abnormal functional interactions between higher-order and primary sensory networks have been implicated in the pathophysiology of visual hallucinations in Parkinson's disease, however the precise signatures remain to be determined. Dimensionality reduction techniques offer a novel means for simplifying the interpretation of multidimensional brain imaging data, identifying hierarchical patterns in the data that are driven by both within- and between-functional network changes. Here, we applied two complementary non-linear dimensionality reduction techniques-diffusion-map embedding and t-distributed stochastic neighbour embedding (t-SNE)-to resting state functional MRI data, in order to characterize the altered functional hierarchy associated with susceptibility to visual hallucinations. Our study involved 77 people with Parkinson's disease (31 with hallucinations; 46 without hallucinations) and 19 age-matched healthy control subjects. In patients with visual hallucinations, we found compression of the unimodal-heteromodal gradient consistent with increased functional integration between sensory and higher order networks. This was mirrored in a traditional functional connectivity analysis, which showed increased connectivity between the visual and default mode networks in the hallucinating group. Together, these results suggest a route by which higher-order regions may have excessive influence over earlier sensory processes, as proposed by theoretical models of hallucinations across disorders. By contrast, the t-SNE analysis identified distinct alterations in prefrontal regions, suggesting an additional layer of complexity in the functional brain network abnormalities implicated in hallucinations, which was not apparent in traditional functional connectivity analyses. Together, the results confirm abnormal brain organization associated with the hallucinating phenotype in Parkinson's disease and highlight the utility of applying convergent dimensionality reduction techniques to investigate complex clinical symptoms. In addition, the patterns we describe in Parkinson's disease converge with those seen in other conditions, suggesting that reduced hierarchical differentiation across sensory-perceptual systems may be a common transdiagnostic vulnerability in neuropsychiatric disorders with perceptual disturbances.

Brain state kinematics and the trajectory of task performance improvement.

Dimensionality reduction techniques offer a unique perspective on brain state dynamics, in which systems-level activity can be tracked through the engagement of a small number of component trajectories. Used in combination with neuroimaging data collected during the performance of cognitive tasks, these approaches can expose the otherwise latent dimensions upon which the brain reconfigures in order to facilitate cognitive performance. Here, we utilized Principal Component Analysis to transform parcellated BOLD timeseries from an fMRI dataset in which 70 human subjects performed an instruction based visuomotor learning task into orthogonal low-dimensional components. We then used Linear Discriminant Analysis to maximise the mean differences between the low-dimensional signatures of fast-and-slow reaction times and early-and-late learners, while also conserving variance present within these groups. The resultant basis set allowed us to describe meaningful differences between these groups and, importantly, to detail the patterns of brain activity which underpin these differences. Our results demonstrate non-linear interactions between three key brain activation maps with convergent trajectories observed at higher task repetitions consistent with optimization. Furthermore, we show subjects with the greatest reaction time improvements have delayed recruitment of left dorsal and lateral prefrontal cortex, as well as deactivation in parts of the occipital lobe and motor cortex, and that the slowest performers have weaker recruitment of somatosensory association cortex and left ventral visual stream, as well as weaker deactivation in the dorsal lateral prefrontal cortex. Overall our results highlight the utility of a kinematic description of brain states, whereby reformatting data into low-dimensional trajectories sensitive to the subtleties of a task can capture non-linear trends in a tractable manner and permit hypothesis generation at the level of brain states.

Core and matrix thalamic sub-populations relate to spatio-temporal cortical connectivity gradients.

Recent neuroimaging experiments have defined low-dimensional gradients of functional connectivity in the cerebral cortex that subserve a spectrum of capacities that span from sensation to cognition. Despite well-known anatomical connections to the cortex, the subcortical areas that support cortical functional organization have been relatively overlooked. One such structure is the thalamus, which maintains extensive anatomical and functional connections with the cerebral cortex across the cortical mantle. The thalamus has a heterogeneous cytoarchitecture, with at least two distinct cell classes that send differential projections to the cortex: granular-projecting 'Core' cells and supragranular-projecting 'Matrix' cells. Here we use high-resolution 7T resting-state fMRI data and the relative amount of two calcium-binding proteins, parvalbumin and calbindin, to infer the relative distribution of these two cell-types (Core and Matrix, respectively) in the thalamus. First, we demonstrate that thalamocortical connectivity recapitulates large-scale, low-dimensional connectivity gradients within the cerebral cortex. Next, we show that diffusely-projecting Matrix regions preferentially correlate with cortical regions with longer intrinsic fMRI timescales. We then show that the Core-Matrix architecture of the thalamus is important for understanding network topology in a manner that supports dynamic integration of signals distributed across the brain. Finally, we replicate our main results in a distinct 3T resting-state fMRI dataset. Linking molecular and functional neuroimaging data, our findings highlight the importance of the thalamic organization for understanding low-dimensional gradients of cortical connectivity.

Quantitative theory of deep brain stimulation of the subthalamic nucleus for the suppression of pathological rhythms in Parkinson's disease.

Deep brain stimulation (DBS) of the subthalamic nucleus (STN) is modeled to explore the mechanisms of this effective, but poorly understood, treatment for motor symptoms of drug-refractory Parkinson's disease and dystonia. First, a neural field model of the corticothalamic-basal ganglia (CTBG) system is developed that reproduces key clinical features of Parkinson's disease, including its characteristic 4-8 Hz and 13-30 Hz electrophysiological signatures. Deep brain stimulation of the STN is then modeled and shown to suppress the pathological 13-30 Hz (beta) activity for physiologically realistic and optimized stimulus parameters. This supports the idea that suppression of abnormally coherent activity in the CTBG system is a major factor in DBS therapy for Parkinson's disease, by permitting normal dynamics to resume. At high stimulus intensities, nonlinear effects in the target population mediate wave-wave interactions between resonant beta activity and the stimulus pulse train, leading to complex spectral structure that shows remarkable similarity to that seen in steady-state evoked potential experiments.

Recruiting neural field theory for data augmentation in a motor imagery brain-computer interface.

We introduce a novel approach to training data augmentation in brain-computer interfaces (BCIs) using neural field theory (NFT) applied to EEG data from motor imagery tasks. BCIs often suffer from limited accuracy due to a limited amount of training data. To address this, we leveraged a corticothalamic NFT model to generate artificial EEG time series as supplemental training data. We employed the BCI competition IV '2a' dataset to evaluate this augmentation technique. For each individual, we fitted the model to common spatial patterns of each motor imagery class, jittered the fitted parameters, and generated time series for data augmentation. Our method led to significant accuracy improvements of over 2% in classifying the "total power" feature, but not in the case of the "Higuchi fractal dimension" feature. This suggests that the fit NFT model may more favorably represent one feature than the other. These findings pave the way for further exploration of NFT-based data augmentation, highlighting the benefits of biophysically accurate artificial data.

Causal evidence for cholinergic stabilization of attractor landscape dynamics.

There is substantial evidence that neuromodulatory systems critically influence brain state dynamics; however, most work has been purely descriptive. Here, we quantify, using data combining local inactivation of the basal forebrain with simultaneous measurement of resting-state fMRI activity in the macaque, the causal role of long-range cholinergic input to the stabilization of brain states in the cerebral cortex. Local inactivation of the nucleus basalis of Meynert (nbM) leads to a decrease in the energy barriers required for an fMRI state transition in cortical ongoing activity. Moreover, the inactivation of particular nbM sub-regions predominantly affects information transfer in cortical regions known to receive direct anatomical projections. We demonstrate these results in a simple neurodynamical model of cholinergic impact on neuronal firing rates and slow hyperpolarizing adaptation currents. We conclude that the cholinergic system plays a critical role in stabilizing macroscale brain state dynamics.

A sensorimotor-association axis of thalamocortical connection development.

Human cortical development follows a sensorimotor-to-association sequence during childhood and adolescence1-6. The brain's capacity to enact this sequence over decades indicates that it relies on intrinsic mechanisms to regulate inter-regional differences in the timing of cortical maturation, yet regulators of human developmental chronology are not well understood. Given evidence from animal models that thalamic axons modulate windows of cortical plasticity7-12, here we evaluate the overarching hypothesis that structural connections between the thalamus and cortex help to coordinate cortical maturational heterochronicity during youth. We first introduce, cortically annotate, and anatomically validate a new atlas of human thalamocortical connections using diffusion tractography. By applying this atlas to three independent youth datasets (ages 8-23 years; total N = 2,676), we reproducibly demonstrate that thalamocortical connections develop along a maturational gradient that aligns with the cortex's sensorimotor-association axis. Associative cortical regions with thalamic connections that take longest to mature exhibit protracted expression of neurochemical, structural, and functional markers indicative of higher circuit plasticity as well as heightened environmental sensitivity. This work highlights a central role for the thalamus in the orchestration of hierarchically organized and environmentally sensitive windows of cortical developmental malleability.

Neuronal connected burst cascades bridge macroscale adaptive signatures across arousal states.

The human brain displays a rich repertoire of states that emerge from the microscopic interactions of cortical and subcortical neurons. Difficulties inherent within large-scale simultaneous neuronal recording limit our ability to link biophysical processes at the microscale to emergent macroscopic brain states. Here we introduce a microscale biophysical network model of layer-5 pyramidal neurons that display graded coarse-sampled dynamics matching those observed in macroscale electrophysiological recordings from macaques and humans. We invert our model to identify the neuronal spike and burst dynamics that differentiate unconscious, dreaming, and awake arousal states and provide insights into their functional signatures. We further show that neuromodulatory arousal can mediate different modes of neuronal dynamics around a low-dimensional energy landscape, which in turn changes the response of the model to external stimuli. Our results highlight the promise of multiscale modelling to bridge theories of consciousness across spatiotemporal scales.

The non-specific matrix thalamus facilitates the cortical information processing modes relevant for conscious awareness.

The neurobiological mechanisms of arousal and anesthesia remain poorly understood. Recent evidence highlights the key role of interactions between the cerebral cortex and the diffusely projecting matrix thalamic nuclei. Here, we interrogate these processes in a whole-brain corticothalamic neural mass model endowed with targeted and diffusely projecting thalamocortical nuclei inferred from empirical data. This model captures key features seen in propofol anesthesia, including diminished network integration, lowered state diversity, impaired susceptibility to perturbation, and decreased corticocortical coherence. Collectively, these signatures reflect a suppression of information transfer across the cerebral cortex. We recover these signatures of conscious arousal by selectively stimulating the matrix thalamus, recapitulating empirical results in macaque, as well as wake-like information processing states that reflect the thalamic modulation of large-scale cortical attractor dynamics. Our results highlight the role of matrix thalamocortical projections in shaping many features of complex cortical dynamics to facilitate the unique communication states supporting conscious awareness.

Parallel processing relies on a distributed, low-dimensional cortico-cerebellar architecture.

A characteristic feature of human cognition is our ability to 'multi-task'-performing two or more tasks in parallel-particularly when one task is well learned. How the brain supports this capacity remains poorly understood. Most past studies have focussed on identifying the areas of the brain-typically the dorsolateral prefrontal cortex-that are required to navigate information-processing bottlenecks. In contrast, we take a systems neuroscience approach to test the hypothesis that the capacity to conduct effective parallel processing relies on a distributed architecture that interconnects the cerebral cortex with the cerebellum. The latter structure contains over half of the neurons in the adult human brain and is well suited to support the fast, effective, dynamic sequences required to perform tasks relatively automatically. By delegating stereotyped within-task computations to the cerebellum, the cerebral cortex can be freed up to focus on the more challenging aspects of performing the tasks in parallel. To test this hypothesis, we analysed task-based fMRI data from 50 participants who performed a task in which they either balanced an avatar on a screen (balance), performed serial-7 subtractions (calculation) or performed both in parallel (dual task). Using a set of approaches that include dimensionality reduction, structure-function coupling, and time-varying functional connectivity, we provide robust evidence in support of our hypothesis. We conclude that distributed interactions between the cerebral cortex and cerebellum are crucially involved in parallel processing in the human brain.

The role of the locus coeruleus in shaping adaptive cortical melodies.

Neural dynamics are shaped and constrained by the projections of a small nucleus in the pons: the noradrenergic locus coeruleus (LC). Much like a bow to the brain's violin, activity in the LC lacks content specificity, but instead dynamically shapes the excitability and receptivity of neurons across the brain. In this review, we explain how the style of the bowing technique, which is analogous to different firing modes in the LC, affects distinct activity patterns in the rest of the brain. Through this analogical lens, we provide intuitive insights into how the complex activity of the LC acts to coordinate adaptive neural dynamics.

The music of the hemispheres: Cortical eigenmodes as a physical basis for large-scale brain activity and connectivity patterns.

Neuroscience has had access to high-resolution recordings of large-scale cortical activity and structure for decades, but still lacks a generally adopted basis to analyze and interrelate results from different individuals and experiments. Here it is argued that the natural oscillatory modes of the cortex-cortical eigenmodes-provide a physically preferred framework for systematic comparisons across experimental conditions and imaging modalities. In this framework, eigenmodes are analogous to notes of a musical instrument, while commonly used statistical patterns parallel frequently played chords. This intuitive perspective avoids problems that often arise in neuroimaging analyses, and connects to underlying mechanisms of brain activity. We envisage this approach will lead to novel insights into whole-brain function, both in existing and prospective datasets, and facilitate a unification of empirical findings across presently disparate analysis paradigms and measurement modalities.

It's about time: Linking dynamical systems with human neuroimaging to understand the brain.

Most human neuroscience research to date has focused on statistical approaches that describe stationary patterns of localized neural activity or blood flow. While these patterns are often interpreted in light of dynamic, information-processing concepts, the static, local, and inferential nature of the statistical approach makes it challenging to directly link neuroimaging results to plausible underlying neural mechanisms. Here, we argue that dynamical systems theory provides the crucial mechanistic framework for characterizing both the brain's time-varying quality and its partial stability in the face of perturbations, and hence, that this perspective can have a profound impact on the interpretation of human neuroimaging results and their relationship with behavior. After briefly reviewing some key terminology, we identify three key ways in which neuroimaging analyses can embrace a dynamical systems perspective: by shifting from a local to a more global perspective, by focusing on dynamics instead of static snapshots of neural activity, and by embracing modeling approaches that map neural dynamics using "forward" models. Through this approach, we envisage ample opportunities for neuroimaging researchers to enrich their understanding of the dynamic neural mechanisms that support a wide array of brain functions, both in health and in the setting of psychopathology.

The ascending arousal system shapes neural dynamics to mediate awareness of cognitive states.

Models of cognitive function typically focus on the cerebral cortex and hence overlook functional links to subcortical structures. This view does not consider the role of the highly-conserved ascending arousal system's role and the computational capacities it provides the brain. We test the hypothesis that the ascending arousal system modulates cortical neural gain to alter the low-dimensional energy landscape of cortical dynamics. Here we use spontaneous functional magnetic resonance imaging data to study phasic bursts in both locus coeruleus and basal forebrain, demonstrating precise time-locked relationships between brainstem activity, low-dimensional energy landscapes, network topology, and spatiotemporal travelling waves. We extend our analysis to a cohort of experienced meditators and demonstrate locus coeruleus-mediated network dynamics were associated with internal shifts in conscious awareness. Together, these results present a view of brain organization that highlights the ascending arousal system's role in shaping both the dynamics of the cerebral cortex and conscious awareness.

Computational models link cellular mechanisms of neuromodulation to large-scale neural dynamics.

Decades of neurobiological research have disclosed the diverse manners in which the response properties of neurons are dynamically modulated to support adaptive cognitive functions. This neuromodulation is achieved through alterations in the biophysical properties of the neuron. However, changes in cognitive function do not arise directly from the modulation of individual neurons, but are mediated by population dynamics in mesoscopic neural ensembles. Understanding this multiscale mapping is an important but nontrivial issue. Here, we bridge these different levels of description by showing how computational models parametrically map classic neuromodulatory processes onto systems-level models of neural activity. The ensuing critical balance of systems-level activity supports perception and action, although our knowledge of this mapping remains incomplete. In this way, quantitative models that link microscale neuronal neuromodulation to systems-level brain function highlight gaps in knowledge and suggest new directions for integrating theoretical and experimental work.

Diffuse neural coupling mediates complex network dynamics through the formation of quasi-critical brain states.

The biological mechanisms that allow the brain to balance flexibility and integration remain poorly understood. A potential solution may lie in a unique aspect of neurobiology, which is that numerous brain systems contain diffuse synaptic connectivity. Here, we demonstrate that increasing diffuse cortical coupling within a validated biophysical corticothalamic model traverses the system through a quasi-critical regime in which spatial heterogeneities in input noise support transient critical dynamics in distributed subregions. The presence of quasi-critical states coincides with known signatures of complex, adaptive brain network dynamics. Finally, we demonstrate the presence of similar dynamic signatures in empirical whole-brain human neuroimaging data. Together, our results establish that modulating the balance between local and diffuse synaptic coupling in a thalamocortical model subtends the emergence of quasi-critical brain states that act to flexibly transition the brain between unique modes of information processing.

The Low-Dimensional Neural Architecture of Cognitive Complexity Is Related to Activity in Medial Thalamic Nuclei.

Cognitive activity emerges from large-scale neuronal dynamics that are constrained to a low-dimensional manifold. How this low-dimensional manifold scales with cognitive complexity, and which brain regions regulate this process, are not well understood. We addressed this issue by analyzing sub-second high-field fMRI data acquired during performance of a task that systematically varied the complexity of cognitive reasoning. We show that task performance reconfigures the low-dimensional manifold and that deviations from these patterns relate to performance errors. We further demonstrate that individual differences in thalamic activity relate to reconfigurations of the low-dimensional architecture during task engagement.

Suppression of Parkinsonian Beta Oscillations by Deep Brain Stimulation: Determination of Effective Protocols.

A neural field model of the corticothalamic-basal ganglia system is developed that describes enhanced beta activity within subthalamic and pallidal circuits in Parkinson's disease (PD) via system resonances. A model of deep brain stimulation (DBS) of typical clinical targets, the subthalamic nucleus (STN) and globus pallidus internus (GPi), is added and studied for several distinct stimulation protocols that are used for treatment of the motor symptoms of PD and that reduce pathological beta band activity (13-30 Hz) in the corticothalamic-basal ganglia network. The resulting impact of DBS on enhanced beta activity in the STN and GPi, as well as cortico-subthalamic and cortico-pallidal coherence, are studied. Both STN-DBS and GPi-DBS are found to be effective for suppressing peak STN and GPi power in the beta band, with GPi-DBS being slightly more effective in both the STN and the GPi for all stimulus protocols tested. The largest decrease in cortico-STN coherence is observed during STN-DBS, whereas GPi-DBS is most effective for reducing cortico-GPi coherence. A reduction of the pathologically large STN connection strengths that define the parkinsonian state results in enhanced 6 Hz activity and could thus represent a compensatory mechanism that has the side effect of driving parkinsonian tremor-like oscillations. This model provides a method for systematically testing effective DBS protocols that agrees with experimental and clinical findings. Furthermore, the model suggests GPi-DBS and STN-DBS have distinct impacts on elevated synchronization between the basal ganglia and motor cortex in PD.