Brain States and Transitions: Insights from Computational Neuroscience

Resting-state dynamic functional connectivity in major depressive disorder: A systematic review

As a novel measure, dynamic functional connectivity (dFC) provides insight into the dynamic nature of brain networks and their interactions in resting-state, surpassing traditional static functional connectivity in pathological conditions such as depression. Since a comprehensive review is still lacking, we then reviewed forty-five eligible papers to explore pathological mechanisms of major depressive disorder (MDD) from perspectives including abnormal brain regions and functional networks, brain state, topological properties, relevant recognition, along with longitudinal studies. Though inconsistencies could be found, common findings are: (1) From different perspectives based on dFC, default-mode network (DMN) with its subregions exhibited a close relation to the pathological mechanism of MDD. (2) With a corrupted integrity within large-scale functional networks and imbalance between them, longer fraction time in a relatively weakly-connected state may be a possible property of MDD concerning its relation with DMN. Abnormal transition frequencies between states were correlated to the severity of MDD. (3) Including dynamic properties in topological network metrics enhanced recognition effect. In all, this review summarized its use for clinical diagnosis and treatment, elucidating the non-stationary of MDD patients' aberrant brain activity in the absence of stimuli and bringing new views into its underlying neuro mechanism.

The model of the brain as a complex system: Interactions of physical, neural and mental states with neurocognitive functions

The isolated approaching of physical, neural and mental states and the binary classification into stable traits and fluctuating states previously lead to a limited understanding concerning underlying processes and possibilities to explain, measure and regulate neural and mental performance along with the interaction of mental states and neurocognitive traits. In this article these states are integrated by i) differentiating the model of the brain as a complex, self-organizing system, ii) showing possibilities to measure this model, iii) offering a classification of mental states and iv) presenting a holistic operationalization of state regulations and trait trainings to enhance neural and mental high-performance on a macro- and micro scale. This model integrates current findings from the theory of constructed emotions, the theory of thousand brains and complex systems theory and yields several testable hypotheses to provide an integrated reference frame for future research and applied target points to regulate and enhance performance.

Inducing a meditative state by artificial perturbations: A mechanistic understanding of brain dynamics underlying meditation

Contemplative neuroscience has increasingly explored meditation using neuroimaging. However, the brain mechanisms underlying meditation remain elusive. Here, we implemented a mechanistic framework to explore the spatiotemporal dynamics of expert meditators during meditation and rest, and controls during rest. We first applied a model-free approach by defining a probabilistic metastable substate (PMS) space for each condition, consisting of different probabilities of occurrence from a repertoire of dynamic patterns. Moreover, we implemented a model-based approach by adjusting the PMS of each condition to a whole-brain model, which enabled us to explore in silico perturbations to transition from resting-state to meditation and vice versa. Consequently, we assessed the sensitivity of different brain areas regarding their perturbability and their mechanistic local-global effects. Overall, our work reveals distinct whole-brain dynamics in meditation compared to rest, and how transitions can be induced with localized artificial perturbations. It motivates future work regarding meditation as a practice in health and as a potential therapy for brain disorders.

In vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here, we noninvasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the gamma-aminobutyric acid (GABA) agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in the association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 y old) and Asian (7.2 to 7.9 y old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

Brain states as wave-like motifs

There is ample evidence of wave-like activity in the brain at multiple scales and levels. This emerging literature supports the broader adoption of a wave perspective of brain activity. Specifically, a brain state can be described as a set of recurring, sequential patterns of propagating brain activity, namely a wave. We examine a collective body of experimental work investigating wave-like properties. Based on these works, we consider brain states as waves using a scale-agnostic framework across time and space. Emphasis is placed on the sequentiality and periodicity associated with brain activity. We conclude by discussing the implications, prospects, and experimental opportunities of this framework.

The Thermodynamics of Mind

To not only survive, but also thrive, the brain must efficiently orchestrate distributed computation across space and time. This requires hierarchical organisation facilitating fast information transfer and processing at the lowest possible metabolic cost. Quantifying brain hierarchy is difficult but can be estimated from the asymmetry of information flow. Thermodynamics has successfully characterised hierarchy in many other complex systems. Here, we propose the 'Thermodynamics of Mind' framework as a natural way to quantify hierarchical brain orchestration and its underlying mechanisms. This has already provided novel insights into the orchestration of hierarchy in brain states including movie watching, where the hierarchy of the brain is flatter than during rest. Overall, this framework holds great promise for revealing the orchestration of cognition.

Insomnia disorder characterized by probabilistic metastable substates using blood-oxygenation-level-dependent (BOLD) phase signals

PURPOSE

From a clinical point of view, how to force a transition from insomnia brain state to healthy brain state by external driven stimulation is of great interest. This needs to define brain state of insomnia disorder as metastable substates. The current study was to identify recurrent substates of insomnia disorder in terms of probability of occurrence, lifetime, and alternation profiles by using leading eigenvector dynamics analysis (LEiDA) method.

METHODS

We enrolled 32 patients with insomnia disorder and 30 healthy subjects. We firstly obtained the BOLD phase coherence matrix from Hilbert transform of BOLD signals and then extracted all the leading eigenvectors from the BOLD phase coherence matrix for all subjects across all time points. Lastly, we clustered the leading eigenvectors using a k-means clustering algorithm to find the probabilistic metastable substates (PMS) and calculate the probability of occurrence and associated lifetime for substates.

RESULTS

The resulting 3 clusters were optimal for brain state of insomnia disorder and healthy brain state, respectively. The occurred probabilities of the PMS were significantly different between the patients with insomnia disorder and healthy subjects, with 0.51 versus 0.44 for PMS-1 (p < 0.001), 0.25 versus 0.27 for PMS-2 (p = 0.051), and 0.24 versus 0.29 for PMS-3 (p < 0.001), as well as the lifetime (in TR) of 36.65 versus 33.15 for PMS-1 (p = 0.068), 14.36 versus 15.43 for PMS-2 (p = 0.117), and 14.80 versus 16.34 for PMS-3 (p = 0.042). The values of the diagonal of the transition matrix were much higher than the probabilities of switching states, indicating the metastable nature of substates.

CONCLUSION

The resulted probabilistic metastable substates hint the characteristic brain dynamics of insomnia disorder. The results may lay a foundation to help determine how to force a transition from insomnia brain state to healthy brain state by external driven stimulation.

Innovation at the Intersection: Emerging Translational Research in Neurology and Psychiatry

Translational research in neurological and psychiatric diseases is a rapidly advancing field that promises to redefine our approach to these complex conditions [...].

Multi-timescale neuromodulation strategy for closed-loop deep brain stimulation in Parkinson's disease

Objective.Beta triggered closed-loop deep brain stimulation (DBS) shows great potential for improving the efficacy while reducing side effect for Parkinson's disease. However, there remain great challenges due to the dynamics and stochasticity of neural activities. In this study, we aimed to tune the amplitude of beta oscillations with different time scales taking into account influence of inherent variations in the basal ganglia-thalamus-cortical circuit.Approach. A dynamic basal ganglia-thalamus-cortical mean-field model was established to emulate the medication rhythm. Then, a dynamic target model was designed to embody the multi-timescale dynamic of beta power with milliseconds, seconds and minutes. Moreover, we proposed a closed-loop DBS strategy based on a proportional-integral-differential (PID) controller with the dynamic control target. In addition, the bounds of stimulation amplitude increments and different parameters of the dynamic target were considered to meet the clinical constraints. The performance of the proposed closed-loop strategy, including beta power modulation accuracy, mean stimulation amplitude, and stimulation variation were calculated to determine the PID parameters and evaluate neuromodulation performance in the computational dynamic mean-field model.Main results. The Results show that the dynamic basal ganglia-thalamus-cortical mean-field model simulated the medication rhythm with the fasted and the slowest rate. The dynamic control target reflected the temporal variation in beta power from milliseconds to minutes. With the proposed closed-loop strategy, the beta power tracked the dynamic target with a smoother stimulation sequence compared with closed-loop DBS with the constant target. Furthermore, the beta power could be modulated to track the control target under different long-term targets, modulation strengths, and bounds of the stimulation increment.Significance. This work provides a new method of closed-loop DBS for multi-timescale beta power modulation with clinical constraints.

Virtual brain twins: from basic neuroscience to clinical use

Virtual brain twins are personalized, generative and adaptive brain models based on data from an individual's brain for scientific and clinical use. After a description of the key elements of virtual brain twins, we present the standard model for personalized whole-brain network models. The personalization is accomplished using a subject's brain imaging data by three means: (1) assemble cortical and subcortical areas in the subject-specific brain space; (2) directly map connectivity into the brain models, which can be generalized to other parameters; and (3) estimate relevant parameters through model inversion, typically using probabilistic machine learning. We present the use of personalized whole-brain network models in healthy ageing and five clinical diseases: epilepsy, Alzheimer's disease, multiple sclerosis, Parkinson's disease and psychiatric disorders. Specifically, we introduce spatial masks for relevant parameters and demonstrate their use based on the physiological and pathophysiological hypotheses. Finally, we pinpoint the key challenges and future directions.

Re-awakening the brain: Forcing transitions in disorders of consciousness by external in silico perturbation

A fundamental challenge in neuroscience is accurately defining brain states and predicting how and where to perturb the brain to force a transition. Here, we investigated resting-state fMRI data of patients suffering from disorders of consciousness (DoC) after coma (minimally conscious and unresponsive wakefulness states) and healthy controls. We applied model-free and model-based approaches to help elucidate the underlying brain mechanisms of patients with DoC. The model-free approach allowed us to characterize brain states in DoC and healthy controls as a probabilistic metastable substate (PMS) space. The PMS of each group was defined by a repertoire of unique patterns (i.e., metastable substates) with different probabilities of occurrence. In the model-based approach, we adjusted the PMS of each DoC group to a causal whole-brain model. This allowed us to explore optimal strategies for promoting transitions by applying off-line in silico probing. Furthermore, this approach enabled us to evaluate the impact of local perturbations in terms of their global effects and sensitivity to stimulation, which is a model-based biomarker providing a deeper understanding of the mechanisms underlying DoC. Our results show that transitions were obtained in a synchronous protocol, in which the somatomotor network, thalamus, precuneus and insula were the most sensitive areas to perturbation. This motivates further work to continue understanding brain function and treatments of disorders of consciousness.

Exploring Motor Network Connectivity in State-Dependent Transcranial Magnetic Stimulation: A Proof-of-Concept Study

State-dependent non-invasive brain stimulation (NIBS) informed by electroencephalography (EEG) has contributed to the understanding of NIBS inter-subject and inter-session variability. While these approaches focus on local EEG characteristics, it is acknowledged that the brain exhibits an intrinsic long-range dynamic organization in networks. This proof-of-concept study explores whether EEG connectivity of the primary motor cortex (M1) in the pre-stimulation period aligns with the Motor Network (MN) and how the MN state affects responses to the transcranial magnetic stimulation (TMS) of M1. One thousand suprathreshold TMS pulses were delivered to the left M1 in eight subjects at rest, with simultaneous EEG. Motor-evoked potentials (MEPs) were measured from the right hand. The source space functional connectivity of the left M1 to the whole brain was assessed using the imaginary part of the phase locking value at the frequency of the sensorimotor μ-rhythm in a 1 s window before the pulse. Group-level connectivity revealed functional links between the left M1, left supplementary motor area, and right M1. Also, pulses delivered at high MN connectivity states result in a greater MEP amplitude compared to low connectivity states. At the single-subject level, this relation is more highly expressed in subjects that feature an overall high cortico-spinal excitability. In conclusion, this study paves the way for MN connectivity-based NIBS.

Release the Krakencoder: A unified brain connectome translation and fusion tool

Brain connectivity can be estimated in many ways, depending on modality and processing strategy. Here we present the Krakencoder, a joint connectome mapping tool that simultaneously, bidirectionally translates between structural (SC) and functional connectivity (FC), and across different atlases and processing choices via a common latent representation. These mappings demonstrate unprecedented accuracy and individual-level identifiability; the mapping between SC and FC has identifiability 42-54% higher than existing models. The Krakencoder combines all connectome flavors via a shared low-dimensional latent space. This "fusion" representation i) better reflects familial relatedness, ii) preserves age- and sex-relevant information and iii) enhances cognition-relevant information. The Krakencoder can be applied without retraining to new, out-of-age-distribution data while still preserving inter-individual differences in the connectome predictions and familial relationships in the latent representations. The Krakencoder is a significant leap forward in capturing the relationship between multi-modal brain connectomes in an individualized, behaviorally- and demographically-relevant way.

Information decomposition and the informational architecture of the brain

To explain how the brain orchestrates information-processing for cognition, we must understand information itself. Importantly, information is not a monolithic entity. Information decomposition techniques provide a way to split information into its constituent elements: unique, redundant, and synergistic information. We review how disentangling synergistic and redundant interactions is redefining our understanding of integrative brain function and its neural organisation. To explain how the brain navigates the trade-offs between redundancy and synergy, we review converging evidence integrating the structural, molecular, and functional underpinnings of synergy and redundancy; their roles in cognition and computation; and how they might arise over evolution and development. Overall, disentangling synergistic and redundant information provides a guiding principle for understanding the informational architecture of the brain and cognition.

The evolution of whole-brain turbulent dynamics during recovery from traumatic brain injury

It has been previously shown that traumatic brain injury (TBI) is associated with reductions in metastability in large-scale networks in resting-state fMRI (rsfMRI). However, little is known about how TBI affects the local level of synchronization and how this evolves during the recovery trajectory. Here, we applied a novel turbulent dynamics framework to investigate whole-brain dynamics using an rsfMRI dataset from a cohort of moderate to severe TBI patients and healthy controls (HCs). We first examined how several measures related to turbulent dynamics differ between HCs and TBI patients at 3, 6, and 12 months post-injury. We found a significant reduction in these empirical measures after TBI, with the largest change at 6 months post-injury. Next, we built a Hopf whole-brain model with coupled oscillators and conducted in silico perturbations to investigate the mechanistic principles underlying the reduced turbulent dynamics found in the empirical data. A simulated attack was used to account for the effect of focal lesions. This revealed a shift to lower coupling parameters in the TBI dataset and, critically, decreased susceptibility and information-encoding capability. These findings confirm the potential of the turbulent framework to characterize longitudinal changes in whole-brain dynamics and in the reactivity to external perturbations after TBI.

In-vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here we non-invasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically-plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the GABA-agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 years old) and Asian (7.2 to 7.9 years old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

Chronic full-band recordings with graphene microtransistors as neural interfaces for discrimination of brain states

Brain states such as sleep, anesthesia, wakefulness, or coma are characterized by specific patterns of cortical activity dynamics, from local circuits to full-brain emergent properties. We previously demonstrated that full-spectrum signals, including the infraslow component (DC, direct current-coupled), can be recorded acutely in multiple sites using flexible arrays of graphene solution-gated field-effect transistors (gSGFETs). Here, we performed chronic implantation of 16-channel gSGFET arrays over the rat cerebral cortex and recorded full-band neuronal activity with two objectives: (1) to test the long-term stability of implanted devices; and (2) to investigate full-band activity during the transition across different levels of anesthesia. First, we demonstrate it is possible to record full-band signals with stability, fidelity, and spatiotemporal resolution for up to 5.5 months using chronic epicortical gSGFET implants. Second, brain states generated by progressive variation of levels of anesthesia could be identified as traditionally using the high-pass filtered (AC, alternating current-coupled) spectrogram: from synchronous slow oscillations in deep anesthesia through to asynchronous activity in the awake state. However, the DC signal introduced a highly significant improvement for brain-state discrimination: the DC band provided an almost linear information prediction of the depth of anesthesia, with about 85% precision, using a trained algorithm. This prediction rose to about 95% precision when the full-band (AC + DC) spectrogram was taken into account. We conclude that recording infraslow activity using gSGFET interfaces is superior for the identification of brain states, and further supports the preclinical and clinical use of graphene neural interfaces for long-term recordings of cortical activity.

Brain dynamics predictive of response to psilocybin for treatment-resistant depression

Psilocybin therapy for depression has started to show promise, yet the underlying causal mechanisms are not currently known. Here, we leveraged the differential outcome in responders and non-responders to psilocybin (10 and 25 mg, 7 days apart) therapy for depression-to gain new insights into regions and networks implicated in the restoration of healthy brain dynamics. We used large-scale brain modelling to fit the spatiotemporal brain dynamics at rest in both responders and non-responders before treatment. Dynamic sensitivity analysis of systematic perturbation of these models enabled us to identify specific brain regions implicated in a transition from a depressive brain state to a healthy one. Binarizing the sample into treatment responders (>50% reduction in depressive symptoms) versus non-responders enabled us to identify a subset of regions implicated in this change. Interestingly, these regions correlate with in vivo density maps of serotonin receptors 5-hydroxytryptamine 2a and 5-hydroxytryptamine 1a, which psilocin, the active metabolite of psilocybin, has an appreciable affinity for, and where it acts as a full-to-partial agonist. Serotonergic transmission has long been associated with depression, and our findings provide causal mechanistic evidence for the role of brain regions in the recovery from depression via psilocybin.

The connectivity of the human frontal pole cortex, and a theory of its involvement in exploit versus explore

The frontal pole is implicated in humans in whether to exploit resources versus explore alternatives. Effective connectivity, functional connectivity, and tractography were measured between six human frontal pole regions and for comparison 13 dorsolateral and dorsal prefrontal cortex regions, and the 360 cortical regions in the Human Connectome Project Multi-modal-parcellation atlas in 171 HCP participants. The frontal pole regions have effective connectivity with Dorsolateral Prefrontal Cortex regions, the Dorsal Prefrontal Cortex, both implicated in working memory; and with the orbitofrontal and anterior cingulate cortex reward/non-reward system. There is also connectivity with temporal lobe, inferior parietal, and posterior cingulate regions. Given this new connectivity evidence, and evidence from activations and damage, it is proposed that the frontal pole cortex contains autoassociation attractor networks that are normally stable in a short-term memory state, and maintain stability in the other prefrontal networks during stable exploitation of goals and strategies. However, if an input from the orbitofrontal or anterior cingulate cortex that expected reward, non-reward, or punishment is received, this destabilizes the frontal pole and thereby other prefrontal networks to enable exploration of competing alternative goals and strategies. The frontal pole connectivity with reward systems may be key in exploit versus explore.

Intrinsic neural timescales relate to the dynamics of infraslow neural waves

The human brain is a highly dynamic organ that operates across a variety of timescales, the intrinsic neural timescales (INT). In addition to the INT, the neural waves featured by its phase-related processes including their cycles with peak/trough and rise/fall play a key role in shaping the brain's neural activity. However, the relationship between the brain's ongoing wave dynamics and INT remains yet unclear. In this study, we utilized functional magnetic resonance imaging (fMRI) rest and task data from the Human Connectome Project (HCP) to investigate the relationship of infraslow wave dynamics [as measured in terms of speed by changes in its peak frequency (PF)] with INT. Our findings reveal that: (i) the speed of phase dynamics (PF) is associated with distinct parts of the ongoing phase cycles, namely higher PF in peak/trough and lower PF in rise/fall; (ii) there exists a negative correlation between phase dynamics (PF) and INT such that slower PF relates to longer INT; (iii) exposure to a movie alters both PF and INT across the different phase cycles, yet their negative correlation remains intact. Collectively, our results demonstrate that INT relates to infraslow phase dynamics during both rest and task states.

Depression recurrence is accompanied by longer periods in default mode and more frequent attentional and reward processing dynamic brain-states during resting-state activity

Recurrence in major depressive disorder (MDD) is common, but neurobiological models capturing vulnerability for recurrences are scarce. Disturbances in multiple resting-state networks have been linked to MDD, but most approaches focus on stable (vs. dynamic) network characteristics. We investigated how the brain's dynamical repertoire changes after patients transition from remission to recurrence of a new depressive episode. Sixty two drug-free, MDD-patients with ≥2 episodes underwent a baseline resting-state fMRI scan when in remission. Over 30-months follow-up, 11 patients with a recurrence and 17 matched-remitted MDD-patients without a recurrence underwent a second fMRI scan. Recurrent patterns of functional connectivity were characterized by applying Leading Eigenvector Dynamics Analysis (LEiDA). Differences between baseline and follow-up were identified for the 11 non-remitted patients, while data from the 17 matched-remitted patients was used as a validation dataset. After the transition into a depressive state, basal ganglia-anterior cingulate cortex (ACC) and visuo-attentional networks were detected significantly more often, whereas default mode network activity was found to have a longer duration. Additionally, the fMRI signal in the basal ganglia-ACC areas underlying the reward network, were significantly less synchronized with the rest of the brain after recurrence (compared to a state of remission). No significant changes were observed in the matched-remitted patients who were scanned twice while in remission. These findings characterize changes that may be associated with the transition from remission to recurrence and provide initial evidence of altered dynamical exploration of the brain's repertoire of functional networks when a recurrent depressive episode occurs.

Violations of the fluctuation-dissipation theorem reveal distinct nonequilibrium dynamics of brain states

The brain is a nonequilibrium system whose dynamics change in different brain states, such as wakefulness and deep sleep. Thermodynamics provides the tools for revealing these nonequilibrium dynamics. We used violations of the fluctuation-dissipation theorem to describe the hierarchy of nonequilibrium dynamics associated with different brain states. Together with a whole-brain model fitted to empirical human neuroimaging data, and deriving the appropriate analytical expressions, we were able to capture the deviation from equilibrium in different brain states that arises from asymmetric interactions and hierarchical organization.

Cortical Synchrony and Information Flow during Transition from Wakefulness to Light Non-Rapid Eye Movement Sleep

Sleep is a highly stereotyped phenomenon, requiring robust spatiotemporal coordination of neural activity. Understanding how the brain coordinates neural activity with sleep onset can provide insights into the physiological functions subserved by sleep and the pathologic phenomena associated with sleep onset. We quantified whole-brain network changes in synchrony and information flow during the transition from wakefulness to light non-rapid eye movement (NREM) sleep, using MEG imaging in a convenient sample of 14 healthy human participants (11 female; mean 63.4 years [SD 11.8 years]). We furthermore performed computational modeling to infer excitatory and inhibitory properties of local neural activity. The transition from wakefulness to light NREM was identified to be encoded in spatially and temporally specific patterns of long-range synchrony. Within the delta band, there was a global increase in connectivity from wakefulness to light NREM, which was highest in frontoparietal regions. Within the theta band, there was an increase in connectivity in fronto-parieto-occipital regions and a decrease in temporal regions from wakefulness to Stage 1 sleep. Patterns of information flow revealed that mesial frontal regions receive hierarchically organized inputs from broad cortical regions upon sleep onset, including direct inflow from occipital regions and indirect inflow via parieto-temporal regions within the delta frequency band. Finally, biophysical neural mass modeling demonstrated changes in the anterior-to-posterior distribution of cortical excitation-to-inhibition with increased excitation-to-inhibition model parameters in anterior regions in light NREM compared with wakefulness. Together, these findings uncover whole-brain corticocortical structure and the orchestration of local and long-range, frequency-specific cortical interactions in the sleep-wake transition.SIGNIFICANCE STATEMENT Our work uncovers spatiotemporal cortical structure of neural synchrony and information flow upon the transition from wakefulness to light non-rapid eye movement sleep. Mesial frontal regions were identified to receive hierarchically organized inputs from broad cortical regions, including both direct inputs from occipital regions and indirect inputs via the parieto-temporal regions within the delta frequency range. Biophysical neural mass modeling revealed a spatially heterogeneous, anterior-posterior distribution of cortical excitation-to-inhibition. Our findings shed light on the orchestration of local and long-range cortical neural structure that is fundamental to sleep onset, and support an emerging view of cortically driven regulation of sleep homeostasis.

Brain state transitions primarily impact the spontaneous rate of slow-firing neurons

The spontaneous firing of neurons is modulated by brain state. Here, we examine how such modulation impacts the overall distribution of firing rates in neuronal populations of neocortical, hippocampal, and thalamic areas across natural and pharmacologically driven brain state transitions. We report that across all the examined combinations of brain area and state transition category, the structure of rate modulation is similar, with almost all fast-firing neurons experiencing proportionally weak modulation, while slow-firing neurons exhibit high inter-neuron variability in the modulation magnitude, leading to a stronger modulation on average. We further demonstrate that this modulation structure is linked to the left-skewed distribution of firing rates on the logarithmic scale and is recapitulated by bivariate log-gamma, but not Gaussian, distributions. Our findings indicate that a preconfigured log-rate distribution with rigid fast-firing neurons and a long left tail of malleable slow-firing neurons is a generic property of forebrain neuronal circuits.

Distinct brain state dynamics of native and second language processing during narrative listening in late bilinguals

The process of complex cognition, which includes language processing, is dynamic in nature and involves various network modes or cognitive modes. This dynamic process can be manifested by a set of brain states and transitions between them. Previous neuroimaging studies have shed light on how bilingual brains support native language (L1) and second language (L2) through a shared network. However, the mechanism through which this shared brain network enables L1 and L2 processing remains unknown. This study examined this issue by testing the hypothesis that L1 and L2 processing is associated with distinct brain state dynamics in terms of brain state integration and transition flexibility. A group of late Chinese-English bilinguals was scanned using functional magnetic resonance imaging (fMRI) while listening to eight short narratives in Chinese (L1) and English (L2). Brain state dynamics were modeled using the leading eigenvector dynamic analysis framework. The results show that L1 processing involves more integrated states and frequent transitions between integrated and segregated states, while L2 processing involves more segregated states and fewer transitions. Our work provides insight into the dynamic process of narrative listening comprehension in late bilinguals and sheds new light on the neural representation of language processing and related disorders.

Characterizing brain stage-dependent pupil dynamics based on lateral hypothalamic activity

Pupil dynamics presents varied correlation features with brain activity under different vigilant levels. The modulation of brain dynamic stages can arise from the lateral hypothalamus (LH), where diverse neuronal cell types contribute to arousal regulation in opposite directions via the anterior cingulate cortex (ACC). However, the relationship of the LH and pupil dynamics has seldom been investigated. Here, we performed local field potential (LFP) recordings at the LH and ACC, and whole-brain fMRI with simultaneous fiber photometry Ca2+ recording in the ACC, to evaluate their correlation with brain state-dependent pupil dynamics. Both LFP and functional magnetic resonance imaging (fMRI) data showed various correlations to pupil dynamics across trials that span negative, null, and positive correlation values, demonstrating brain state-dependent coupling features. Our results indicate that the correlation of pupil dynamics with ACC LFP and whole-brain fMRI signals depends on LH activity, suggesting a role of the latter in brain dynamic stage regulation.

Hierarchical organization of the human ventral visual streams revealed with magnetoencephalography

The hierarchical organization between 25 ventral stream visual cortical regions and 180 cortical regions was measured with magnetoencephalography using the Human Connectome Project Multimodal Parcellation atlas in 83 Human Connectome Project participants performing a visual memory task. The aim was to reveal the hierarchical organization using a whole-brain model based on generative effective connectivity with this fast neuroimaging method. V1-V4 formed a first group of interconnected regions. Especially V4 had connectivity to a ventrolateral visual stream: V8, the fusiform face cortex, and posterior inferior temporal cortex PIT. These regions in turn had effectivity connectivity to inferior temporal cortex visual regions TE2p and TE1p. TE2p and TE1p then have connectivity to anterior temporal lobe regions TE1a, TE1m, TE2a, and TGv, which are multimodal. In a ventromedial visual stream, V1-V4 connect to ventromedial regions VMV1-3 and VVC. VMV1-3 and VVC connect to the medial parahippocampal gyrus PHA1-3, which, with the VMV regions, include the parahippocampal scene area. The medial parahippocampal PHA1-3 regions have connectivity to the hippocampal system regions the perirhinal cortex, entorhinal cortex, and hippocampus. These effective connectivities of two ventral visual cortical streams measured with magnetoencephalography provide support to the hierarchical organization of brain systems measured with fMRI, and new evidence on directionality.

Establishing brain states in neuroimaging data

The definition of a brain state remains elusive, with varying interpretations across different sub-fields of neuroscience-from the level of wakefulness in anaesthesia, to activity of individual neurons, voltage in EEG, and blood flow in fMRI. This lack of consensus presents a significant challenge to the development of accurate models of neural dynamics. However, at the foundation of dynamical systems theory lies a definition of what constitutes the 'state' of a system-i.e., a specification of the system's future. Here, we propose to adopt this definition to establish brain states in neuroimaging timeseries by applying Dynamic Causal Modelling (DCM) to low-dimensional embedding of resting and task condition fMRI data. We find that ~90% of subjects in resting conditions are better described by first-order models, whereas ~55% of subjects in task conditions are better described by second-order models. Our work calls into question the status quo of using first-order equations almost exclusively within computational neuroscience and provides a new way of establishing brain states, as well as their associated phase space representations, in neuroimaging datasets.

Cross-attractor modeling of resting-state functional connectivity in psychiatric disorders

Resting-state functional connectivity (RSFC) is altered across various psychiatric disorders. Brain network modeling (BNM) has the potential to reveal the neurobiological underpinnings of such abnormalities by dynamically modeling the structure-function relationship and examining biologically relevant parameters after fitting the models with real data. Although innovative BNM approaches have been developed, two main issues need to be further addressed. First, previous BNM approaches are primarily limited to simulating noise-driven dynamics near a chosen attractor (or a stable brain state). An alternative approach is to examine multi(or cross)-attractor dynamics, which can be used to better capture non-stationarity and switching between states in the resting brain. Second, previous BNM work is limited to characterizing one disorder at a time. Given the large degree of co-morbidity across psychiatric disorders, comparing BNMs across disorders might provide a novel avenue to generate insights regarding the dynamical features that are common across (vs. specific to) disorders. Here, we address these issues by (1) examining the layout of the attractor repertoire over the entire multi-attractor landscape using a recently developed cross-attractor BNM approach; and (2) characterizing and comparing multiple disorders (schizophrenia, bipolar, and ADHD) with healthy controls using an openly available and moderately large multimodal dataset from the UCLA Consortium for Neuropsychiatric Phenomics. Both global and local differences were observed across disorders. Specifically, the global coupling between regions was significantly decreased in schizophrenia patients relative to healthy controls. At the same time, the ratio between local excitation and inhibition was significantly higher in the schizophrenia group than the ADHD group. In line with these results, the schizophrenia group had the lowest switching costs (energy gaps) across groups for several networks including the default mode network. Paired comparison also showed that schizophrenia patients had significantly lower energy gaps than healthy controls for the somatomotor and visual networks. Overall, this study provides preliminary evidence supporting transdiagnostic multi-attractor BNM approaches to better understand psychiatric disorders' pathophysiology.

Controlling brain dynamics: Landscape and transition path for working memory

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

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.

The Mystery of Mental Integrity: Clarifying Its Relevance to Neurotechnologies

The concept of mental integrity is currently a significant topic in discussions concerning the regulation of neurotechnologies. Technologies such as deep brain stimulation and brain-computer interfaces are believed to pose a unique threat to mental integrity, and some authors have advocated for a legal right to protect it. Despite this, there remains uncertainty about what mental integrity entails and why it is important. Various interpretations of the concept have been proposed, but the literature on the subject is inconclusive. Here we consider a number of possible interpretations and argue that the most plausible one concerns neurotechnologies that bypass one's reasoning capacities, and do so specifically in ways that reliably lead to alienation from one's mental states. This narrows the scope of what constitutes a threat to mental integrity and offers a more precise role for the concept to play in the ethical evaluation of neurotechnologies.

Computational modeling allows unsupervised classification of epileptic brain states across species

Current advances in epilepsy treatment aim to personalize and responsively adjust treatment parameters to overcome patient heterogeneity in treatment efficiency. For tailoring treatment to the individual and the current brain state, tools are required that help to identify the patient- and time-point-specific parameters of epilepsy. Computational modeling has long proven its utility in gaining mechanistic insight. Recently, the technique has been introduced as a diagnostic tool to predict individual treatment outcomes. In this article, the Wendling model, an established computational model of epilepsy dynamics, is used to automatically classify epileptic brain states in intracranial EEG from patients (n = 4) and local field potential recordings from in vitro rat data (high-potassium model of epilepsy, n = 3). Five-second signal segments are classified to four types of brain state in epilepsy (interictal, preonset, onset, ictal) by comparing a vector of signal features for each data segment to four prototypical feature vectors obtained by Wendling model simulations. The classification result is validated against expert visual assessment. Model-driven brain state classification achieved a classification performance significantly above chance level (mean sensitivity 0.99 on model data, 0.77 on rat data, 0.56 on human data in a four-way classification task). Model-driven prototypes showed similarity with data-driven prototypes, which we obtained from real data for rats and humans. Our results indicate similar electrophysiological patterns of epileptic states in the human brain and the animal model that are well-reproduced by the computational model, and captured by a key set of signal features, enabling fully automated and unsupervised brain state classification in epilepsy.

Computational modelling in disorders of consciousness: Closing the gap towards personalised models for restoring consciousness

Disorders of consciousness are complex conditions characterised by persistent loss of responsiveness due to brain injury. They present diagnostic challenges and limited options for treatment, and highlight the urgent need for a more thorough understanding of how human consciousness arises from coordinated neural activity. The increasing availability of multimodal neuroimaging data has given rise to a wide range of clinically- and scientifically-motivated modelling efforts, seeking to improve data-driven stratification of patients, to identify causal mechanisms for patient pathophysiology and loss of consciousness more broadly, and to develop simulations as a means of testing in silico potential treatment avenues to restore consciousness. As a dedicated Working Group of clinicians and neuroscientists of the international Curing Coma Campaign, here we provide our framework and vision to understand the diverse statistical and generative computational modelling approaches that are being employed in this fast-growing field. We identify the gaps that exist between the current state-of-the-art in statistical and biophysical computational modelling in human neuroscience, and the aspirational goal of a mature field of modelling disorders of consciousness; which might drive improved treatments and outcomes in the clinic. Finally, we make several recommendations for how the field as a whole can work together to address these challenges.

In vivo mapping of pharmacologically induced functional reorganization onto the human brain's neurotransmitter landscape

To understand how pharmacological interventions can exert their powerful effects on brain function, we need to understand how they engage the brain's rich neurotransmitter landscape. Here, we bridge microscale molecular chemoarchitecture and pharmacologically induced macroscale functional reorganization, by relating the regional distribution of 19 neurotransmitter receptors and transporters obtained from positron emission tomography, and the regional changes in functional magnetic resonance imaging connectivity induced by 10 different mind-altering drugs: propofol, sevoflurane, ketamine, lysergic acid diethylamide (LSD), psilocybin, N,N-Dimethyltryptamine (DMT), ayahuasca, 3,4-methylenedioxymethamphetamine (MDMA), modafinil, and methylphenidate. Our results reveal a many-to-many mapping between psychoactive drugs' effects on brain function and multiple neurotransmitter systems. The effects of both anesthetics and psychedelics on brain function are organized along hierarchical gradients of brain structure and function. Last, we show that regional co-susceptibility to pharmacological interventions recapitulates co-susceptibility to disorder-induced structural alterations. Collectively, these results highlight rich statistical patterns relating molecular chemoarchitecture and drug-induced reorganization of the brain's functional architecture.

Low-dimensional organization of global brain states of reduced consciousness

Brain states are frequently represented using a unidimensional scale measuring the richness of subjective experience (level of consciousness). This description assumes a mapping between the high-dimensional space of whole-brain configurations and the trajectories of brain states associated with changes in consciousness, yet this mapping and its properties remain unclear. We combine whole-brain modeling, data augmentation, and deep learning for dimensionality reduction to determine a mapping representing states of consciousness in a low-dimensional space, where distances parallel similarities between states. An orderly trajectory from wakefulness to patients with brain injury is revealed in a latent space whose coordinates represent metrics related to functional modularity and structure-function coupling, increasing alongside loss of consciousness. Finally, we investigate the effects of model perturbations, providing geometrical interpretation for the stability and reversibility of states. We conclude that conscious awareness depends on functional patterns encoded as a low-dimensional trajectory within the vast space of brain configurations.

Auditory cortical connectivity in humans

To understand auditory cortical processing, the effective connectivity between 15 auditory cortical regions and 360 cortical regions was measured in 171 Human Connectome Project participants, and complemented with functional connectivity and diffusion tractography. 1. A hierarchy of auditory cortical processing was identified from Core regions (including A1) to Belt regions LBelt, MBelt, and 52; then to PBelt; and then to HCP A4. 2. A4 has connectivity to anterior temporal lobe TA2, and to HCP A5, which connects to dorsal-bank superior temporal sulcus (STS) regions STGa, STSda, and STSdp. These STS regions also receive visual inputs about moving faces and objects, which are combined with auditory information to help implement multimodal object identification, such as who is speaking, and what is being said. Consistent with this being a "what" ventral auditory stream, these STS regions then have effective connectivity to TPOJ1, STV, PSL, TGv, TGd, and PGi, which are language-related semantic regions connecting to Broca's area, especially BA45. 3. A4 and A5 also have effective connectivity to MT and MST, which connect to superior parietal regions forming a dorsal auditory "where" stream involved in actions in space. Connections of PBelt, A4, and A5 with BA44 may form a language-related dorsal stream.

From Neurotransmitters to Networks: Transcending Organisational Hierarchies with Molecular-informed Functional Imaging

The human brain exhibits complex interactions across micro, meso-, and macro-scale organisational principles. Recent synergistic multi-modal approaches have begun to link micro-scale information to systems level dynamics, transcending organisational hierarchies and offering novel perspectives into the brain's function and dysfunction. Specifically, the distribution of micro-scale properties (such as receptor density or gene expression) can be mapped onto macro-scale measures from functional MRI to provide novel neurobiological insights. Methodological approaches to enrich functional imaging analyses with molecular information are rapidly evolving, with several streams of research having developed relatively independently, each offering unique potential to explore the trans-hierarchical functioning of the brain. Here, we address the three principal streams of research - spatial correlation, molecular-enriched network, and in-silico whole brain modelling analyses - to provide a critical overview of the different sources of molecular information, how this information can be utilised within analyses of fMRI data, the merits and pitfalls of each methodology, and, through the use of key examples, highlight their promise to shed new light on key domains of neuroscientific inquiry.

Reduced emergent character of neural dynamics in patients with a disrupted connectome

High-level brain functions are widely believed to emerge from the orchestrated activity of multiple neural systems. However, lacking a formal definition and practical quantification of emergence for experimental data, neuroscientists have been unable to empirically test this long-standing conjecture. Here we investigate this fundamental question by leveraging a recently proposed framework known as "Integrated Information Decomposition," which establishes a principled information-theoretic approach to operationalise and quantify emergence in dynamical systems - including the human brain. By analysing functional MRI data, our results show that the emergent and hierarchical character of neural dynamics is significantly diminished in chronically unresponsive patients suffering from severe brain injury. At a functional level, we demonstrate that emergence capacity is positively correlated with the extent of hierarchical organisation in brain activity. Furthermore, by combining computational approaches from network control theory and whole-brain biophysical modelling, we show that the reduced capacity for emergent and hierarchical dynamics in severely brain-injured patients can be mechanistically explained by disruptions in the patients' structural connectome. Overall, our results suggest that chronic unresponsiveness resulting from severe brain injury may be related to structural impairment of the fundamental neural infrastructures required for brain dynamics to support emergence.

The human posterior parietal cortex: effective connectome, and its relation to function

The effective connectivity between 21 regions in the human posterior parietal cortex, and 360 cortical regions was measured in 171 Human Connectome Project (HCP) participants using the HCP atlas, and complemented with functional connectivity and diffusion tractography. Intraparietal areas LIP, VIP, MIP, and AIP have connectivity from early cortical visual regions, and to visuomotor regions such as the frontal eye fields, consistent with functions in eye saccades and tracking. Five superior parietal area 7 regions receive from similar areas and from the intraparietal areas, but also receive somatosensory inputs and connect with premotor areas including area 6, consistent with functions in performing actions to reach for, grasp, and manipulate objects. In the anterior inferior parietal cortex, PFop, PFt, and PFcm are mainly somatosensory, and PF in addition receives visuo-motor and visual object information, and is implicated in multimodal shape and body image representations. In the posterior inferior parietal cortex, PFm and PGs combine visuo-motor, visual object, and reward input and connect with the hippocampal system. PGi in addition provides a route to motion-related superior temporal sulcus regions involved in social interactions. PGp has connectivity with intraparietal regions involved in coordinate transforms and may be involved in idiothetic update of hippocampal visual scene representations.

Reduced resilience of brain state transitions in anti-N-methyl-D-aspartate receptor encephalitis

Patients with anti-N-methyl-aspartate receptor (NMDA) receptor encephalitis suffer from a severe neuropsychiatric syndrome, yet most patients show no abnormalities in routine magnetic resonance imaging. In contrast, advanced neuroimaging studies have consistently identified disrupted functional connectivity in these patients, with recent work suggesting increased volatility of functional state dynamics. Here, we investigate these network dynamics through the spatiotemporal trajectory of meta-state transitions, yielding a time-resolved account of brain state exploration in anti-NMDA receptor encephalitis. To this end, resting-state functional magnetic resonance imaging data were acquired in 73 patients with anti-NMDA receptor encephalitis and 73 age- and sex-matched healthy controls. Time-resolved functional connectivity was clustered into brain meta-states, giving rise to a time-resolved transition network graph with states as nodes and transitions between brain meta-states as weighted, directed edges. Network topology, robustness and transition cost of these transition networks were compared between groups. Transition networks of patients showed significantly lower local efficiency (t = -2.41, p_FDR  = .029), lower robustness (t = -2.01, p_FDR  = .048) and higher leap size (t = 2.18, p_FDR  = .037) compared with controls. Furthermore, the ratio of within-to-between module transitions and state similarity was significantly lower in patients. Importantly, alterations of brain state transitions correlated with disease severity. Together, these findings reveal systematic alterations of transition networks in patients, suggesting that anti-NMDA receptor encephalitis is characterized by reduced stability of brain state transitions and that this reduced resilience of transition networks plays a clinically relevant role in the manifestation of the disease.

The human posterior cingulate, retrosplenial, and medial parietal cortex effective connectome, and implications for memory and navigation

The human posterior cingulate, retrosplenial, and medial parietal cortex are involved in memory and navigation. The functional anatomy underlying these cognitive functions was investigated by measuring the effective connectivity of these Posterior Cingulate Division (PCD) regions in the Human Connectome Project-MMP1 atlas in 171 HCP participants, and complemented with functional connectivity and diffusion tractography. First, the postero-ventral parts of the PCD (31pd, 31pv, 7m, d23ab, and v23ab) have effective connectivity with the temporal pole, inferior temporal visual cortex, cortex in the superior temporal sulcus implicated in auditory and semantic processing, with the reward-related vmPFC and pregenual anterior cingulate cortex, with the inferior parietal cortex, and with the hippocampal system. This connectivity implicates it in hippocampal episodic memory, providing routes for "what," reward and semantic schema-related information to access the hippocampus. Second, the antero-dorsal parts of the PCD (especially 31a and 23d, PCV, and also RSC) have connectivity with early visual cortical areas including those that represent spatial scenes, with the superior parietal cortex, with the pregenual anterior cingulate cortex, and with the hippocampal system. This connectivity implicates it in the "where" component for hippocampal episodic memory and for spatial navigation. The dorsal-transitional-visual (DVT) and ProStriate regions where the retrosplenial scene area is located have connectivity from early visual cortical areas to the parahippocampal scene area, providing a ventromedial route for spatial scene information to reach the hippocampus. These connectivities provide important routes for "what," reward, and "where" scene-related information for human hippocampal episodic memory and navigation. The midcingulate cortex provides a route from the anterior dorsal parts of the PCD and the supracallosal part of the anterior cingulate cortex to premotor regions.

Progressive alteration of dynamic functional connectivity patterns in subcortical ischemic vascular cognitive impairment patients

Alterations in the temporal evolution of brain states in the process of cognitive impairment aggravation due to subcortical ischemic vascular disease (SIVD) is not understood. The dynamic functional connectivity was investigated to identify the abnormal temporal properties of brain states associated with cognitive impairment caused by SIVD. Eighteen patients with subcortical ischemic vascular cognitive impairment with no dementia (SIVCIND), 19 dementia patients (SIVaD) and 26 normal controls were enrolled. We found that the occupancy rate and mean lifetime of brain states were associated with cognitive performance. SIVCIND had a higher occupancy rate and longer mean lifetime in weakly connected states than normal controls. SIVaD had similar but more extensive changes in the temporal properties of brain states. In addition, switching from weakly connected states to more strongly connected states was more difficult in SIVCIND and SIVaD patients than in normal controls, especially in SIVaD patients. The results revealed that not only the transition to but also maintenance in strongly connected states became increasingly difficult when SIVD-related cognitive impairment progressed into a more severe stage.

Distributed harmonic patterns of structure-function dependence orchestrate human consciousness

A central question in neuroscience is how consciousness arises from the dynamic interplay of brain structure and function. Here we decompose functional MRI signals from pathological and pharmacologically-induced perturbations of consciousness into distributed patterns of structure-function dependence across scales: the harmonic modes of the human structural connectome. We show that structure-function coupling is a generalisable indicator of consciousness that is under bi-directional neuromodulatory control. We find increased structure-function coupling across scales during loss of consciousness, whether due to anaesthesia or brain injury, capable of discriminating between behaviourally indistinguishable sub-categories of brain-injured patients, tracking the presence of covert consciousness. The opposite harmonic signature characterises the altered state induced by LSD or ketamine, reflecting psychedelic-induced decoupling of brain function from structure and correlating with physiological and subjective scores. Overall, connectome harmonic decomposition reveals how neuromodulation and the network architecture of the human connectome jointly shape consciousness and distributed functional activation across scales.

Toward naturalistic neuroscience: Mechanisms underlying the flattening of brain hierarchy in movie-watching compared to rest and task

Identifying the functional specialization of the brain has moved from using cognitive tasks and resting state to using ecological relevant, naturalistic movies. We leveraged a large-scale neuroimaging dataset to directly investigate the hierarchical reorganization of functional brain activity when watching naturalistic films compared to performing seven cognitive tasks and resting. A thermodynamics-inspired whole-brain model paradigm revealed the generative underlying mechanisms for changing the balance in causal interactions between brain regions in different conditions. Paradoxically, the hierarchy is flatter for movie-watching, and the level of nonreversibility is significantly smaller in comparison to both rest and tasks, where the latter in turn have the highest levels of hierarchy and nonreversibility. The underlying mechanisms were revealed by the model-based generative effective connectivity (GEC). Naturalistic films could therefore provide a fast and convenient way to measure important changes in GEC (integrating functional and anatomical connectivity) found in, for example, neuropsychiatric disorders. Overall, this study demonstrates the benefits of moving toward a more naturalistic neuroscience.

Optimal Control Costs of Brain State Transitions in Linear Stochastic Systems

The brain is a system that performs numerous functions by controlling its states. Quantifying the cost of this control is essential as it reveals how the brain can be controlled based on the minimization of the control cost, and which brain regions are most important to the optimal control of transitions. Despite its great potential, the current control paradigm in neuroscience uses a deterministic framework and is therefore unable to consider stochasticity, severely limiting its application to neural data. Here, to resolve this limitation, we propose a novel framework for the evaluation of control costs based on a linear stochastic model. Following our previous work, we quantified the optimal control cost as the minimal Kullback-Leibler divergence between the uncontrolled and controlled processes. In the linear model, we established an analytical expression for minimal cost and showed that we can decompose it into the cost for controlling the mean and covariance of brain activity. To evaluate the utility of our novel framework, we examined the significant brain regions in the optimal control of transitions from the resting state to seven cognitive task states in human whole-brain imaging data of either sex. We found that, in realizing the different transitions, the lower visual areas commonly played a significant role in controlling the means, while the posterior cingulate cortex commonly played a significant role in controlling the covariances.SIGNIFICANCE STATEMENT The brain performs many cognitive functions by controlling its states. Quantifying the cost of this control is essential as it reveals how the brain can be optimally controlled in terms of the cost, and which brain regions are most important to the optimal control of transitions. Here, we built a novel framework to quantify control cost that takes account of stochasticity of neural activity, which is ignored in previous studies. We established the analytical expression of the stochastic control cost, which enables us to compute the cost in high-dimensional neural data. We identified the significant brain regions for the optimal control in cognitive tasks in human whole-brain imaging data.

Human amygdala compared to orbitofrontal cortex connectivity, and emotion

The amygdala and orbitofrontal cortex have been implicated in emotion. To understand these regions better in humans, their effective connectivity with 360 cortical regions was measured in 171 humans from the Human Connectome Project, and complemented with functional connectivity and diffusion tractography. The human amygdala has effective connectivity from few cortical regions compared to the orbitofrontal cortex: primarily from auditory cortex A5 and the related superior temporal gyrus and temporal pole regions; the piriform (olfactory) cortex; the lateral orbitofrontal cortex 47m; somatosensory cortex; the hippocampus, entorhinal cortex, perirhinal cortex, and parahippocampal TF; and from the cholinergic nucleus basalis. The amygdala has effective connectivity to the hippocampus, entorhinal and perirhinal cortex; to the temporal pole; and to the lateral orbitofrontal cortex. The orbitofrontal cortex has effective connectivity from gustatory, olfactory, and temporal visual, auditory and pole cortex, and to the pregenual anterior and posterior cingulate cortex, hippocampal system, and prefrontal cortex, and provides for rewards and punishers to be used in reported emotions, and memory and navigation to goals. Given the paucity of amygdalo-neocortical connectivity in humans, it is proposed that the human amygdala is involved primarily in autonomic and conditioned responses via brainstem connectivity, rather than in reported (declarative) emotion.

The human orbitofrontal cortex, vmPFC, and anterior cingulate cortex effective connectome: emotion, memory, and action

The human orbitofrontal cortex, ventromedial prefrontal cortex (vmPFC), and anterior cingulate cortex are involved in reward processing and thereby in emotion but are also implicated in episodic memory. To understand these regions better, the effective connectivity between 360 cortical regions and 24 subcortical regions was measured in 172 humans from the Human Connectome Project and complemented with functional connectivity and diffusion tractography. The orbitofrontal cortex has effective connectivity from gustatory, olfactory, and temporal visual, auditory, and pole cortical areas. The orbitofrontal cortex has connectivity to the pregenual anterior and posterior cingulate cortex and hippocampal system and provides for rewards to be used in memory and navigation to goals. The orbitofrontal and pregenual anterior cortex have connectivity to the supracallosal anterior cingulate cortex, which projects to midcingulate and other premotor cortical areas and provides for action-outcome learning including limb withdrawal or flight or fight to aversive and nonreward stimuli. The lateral orbitofrontal cortex has outputs to language systems in the inferior frontal gyrus. The medial orbitofrontal cortex connects to the nucleus basalis of Meynert and the pregenual cingulate to the septum, and damage to these cortical regions may contribute to memory impairments by disrupting cholinergic influences on the neocortex and hippocampus.