Temporal sequences of brain activity at rest are constrained by white matter structure and modulated by cognitive demands

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.

Evolutionarily conserved fMRI network dynamics in the mouse, macaque, and human brain

Evolutionarily relevant networks have been previously described in several mammalian species using time-averaged analyses of fMRI time-series. However, fMRI network activity is highly dynamic and continually evolves over timescales of seconds. Whether the dynamic organization of resting-state fMRI network activity is conserved across mammalian species remains unclear. Using frame-wise clustering of fMRI time-series, we find that intrinsic fMRI network dynamics in awake male macaques and humans is characterized by recurrent transitions between a set of 4 dominant, neuroanatomically homologous fMRI coactivation modes (C-modes), three of which are also plausibly represented in the male rodent brain. Importantly, in all species C-modes exhibit species-invariant dynamic features, including preferred occurrence at specific phases of fMRI global signal fluctuations, and a state transition structure compatible with infraslow coupled oscillator dynamics. Moreover, dominant C-mode occurrence reconstitutes the static organization of the fMRI connectome in all species, and is predictive of ranking of corresponding fMRI connectivity gradients. These results reveal a set of species-invariant principles underlying the dynamic organization of fMRI networks in mammalian species, and offer novel opportunities to relate fMRI network findings across the phylogenetic tree.

Altered Structural Connectivity and Functional Brain Dynamics in Individuals With Heavy Alcohol Use Elucidated via Network Control Theory

BACKGROUND

Heavy alcohol use and its associated conditions, such as alcohol use disorder, impact millions of individuals worldwide. While our understanding of the neurobiological correlates of alcohol use has evolved substantially, we still lack models that incorporate whole-brain neuroanatomical, functional, and pharmacological information under one framework.

METHODS

Here, we utilized diffusion and functional magnetic resonance imaging to investigate alterations to brain dynamics in 130 individuals with a high amount of current alcohol use. We compared these alcohol-using individuals to 308 individuals with minimal use of any substances.

RESULTS

We found that individuals with heavy alcohol use had less dynamic and complex brain activity, and through leveraging network control theory, had increased control energy to complete transitions between activation states. Furthermore, using separately acquired positron emission tomography data, we deployed an in silico evaluation demonstrating that decreased D2 receptor levels, as found previously in individuals with alcohol use disorder, may relate to our observed findings.

CONCLUSIONS

This work demonstrates that whole-brain, multimodal imaging information can be combined under a network control framework to identify and evaluate neurobiological correlates and mechanisms of heavy alcohol use.

Dynamical structure-function correlations provide robust and generalizable signatures of consciousness in humans

Resting-state functional magnetic resonance imaging evolves through a repertoire of functional connectivity patterns which might reflect ongoing cognition, as well as the contents of conscious awareness. We investigated whether the dynamic exploration of these states can provide robust and generalizable markers for the state of consciousness in human participants, across loss of consciousness induced by general anaesthesia or slow wave sleep. By clustering transient states of functional connectivity, we demonstrated that brain activity during unconsciousness is dominated by a recurrent pattern primarily mediated by structural connectivity and with a reduced capacity to transition to other patterns. Our results provide evidence supporting the pronounced differences between conscious and unconscious brain states in terms of whole-brain dynamics; in particular, the maintenance of rich brain dynamics measured by entropy is a critical aspect of conscious awareness. Collectively, our results may have significant implications for our understanding of consciousness and the neural basis of human awareness, as well as for the discovery of robust signatures of consciousness that are generalizable among different brain conditions.

Adaptive Strategies of Single-Sided Deaf Cochlear-Implant Users Revealed Through Resting State Activity: an Auditory PET Brain Imaging Study

Brain plasticity refers to the brain's ability to reorganize its structure or function in response to experiences, learning, and environmental influences. This phenomenon is particularly significant in individuals with deafness, as the brain adapts to compensate for the lack of auditory stimulation. The aim of this study is to investigate whether cochlear implantation can restore a normal pattern of brain activation following auditory stimulation in cases of asymmetric hearing loss. We used a PET-scan technique to assess brain activity after cochlear implantation, specifically during an auditory voice/non-voice discrimination task. The results indicated a nearly normal pattern of brain activity during the auditory discrimination task, except for increased activation in areas related to attentional processes compared to controls. Additionally, brain activity at rest showed significant changes in implanted participants, including cross modal visuo-auditory processing. Therefore, cochlear implants can restore the brain's activation pattern through long-term adaptive adjustments in intrinsic brain activity.

Sex-specific differences in brain activity dynamics of youth with a family history of substance use disorder

An individual's risk of substance use disorder (SUD) is shaped by a complex interplay of potent biosocial factors. Current neurodevelopmental models posit vulnerability to SUD in youth is due to an overreactive reward system and reduced inhibitory control. Having a family history of SUD is a particularly strong risk factor, yet few studies have explored its impact on brain function and structure prior to substance exposure. Herein, we utilized a network control theory approach to quantify sex-specific differences in brain activity dynamics in youth with and without a family history of SUD, drawn from a large cohort of substance-naïve youth from the Adolescent Brain Cognitive Development Study. We summarize brain dynamics by calculating transition energy, which probes the ease with which a whole brain, region or network drives the brain towards a specific spatial pattern of activation (i.e., brain state). Our findings reveal that a family history of SUD is associated with alterations in the brain's dynamics wherein: i) independent of sex, certain regions' transition energies are higher in those with a family history of SUD and ii) there exist sex-specific differences in SUD family history groups at multiple levels of transition energy (global, network, and regional). Family history-by-sex effects reveal that energetic demand is increased in females with a family history of SUD and decreased in males with a family history of SUD, compared to their same-sex counterparts with no SUD family history. Specifically, we localize these effects to higher energetic demands of the default mode network in females with a family history of SUD and lower energetic demands of attention networks in males with a family history of SUD. These results suggest a family history of SUD may increase reward saliency in males and decrease efficiency of top-down inhibitory control in females. This work could be used to inform personalized intervention strategies that may target differing cognitive mechanisms that predispose individuals to the development of SUD.

Contributions of network structure, chemoarchitecture and diagnostic categories to transitions between cognitive topographies

The mechanisms linking the brain's network structure to cognitively relevant activation patterns remain largely unknown. Here, by leveraging principles of network control, we show how the architecture of the human connectome shapes transitions between 123 experimentally defined cognitive activation maps (cognitive topographies) from the NeuroSynth meta-analytic database. Specifically, we systematically integrated large-scale multimodal neuroimaging data from functional magnetic resonance imaging, diffusion tractography, cortical morphometry and positron emission tomography to simulate how anatomically guided transitions between cognitive states can be reshaped by neurotransmitter engagement or by changes in cortical thickness. Our model incorporates neurotransmitter-receptor density maps (18 receptors and transporters) and maps of cortical thickness pertaining to a wide range of mental health, neurodegenerative, psychiatric and neurodevelopmental diagnostic categories (17,000 patients and 22,000 controls). The results provide a comprehensive look-up table charting how brain network organization and chemoarchitecture interact to manifest different cognitive topographies, and establish a principled foundation for the systematic identification of ways to promote selective transitions between cognitive topographies.

Sleep Deprivation Effects on Brain State Dynamics Are Associated With Dopamine D2 Receptor Availability Via Network Control Theory

BACKGROUND

Sleep deprivation (SD) negatively affects brain function. Most brain imaging studies have investigated the effects of SD on static brain function. SD effects on functional brain dynamics and their relationship with molecular changes remain relatively unexplored.

METHODS

We used functional magnetic resonance imaging to examine resting-brain state dynamics after one night of SD compared with rested wakefulness (N = 41) and assessed the association of brain state dynamics with striatal brain dopamine D2 receptor availability measured by positron emission tomography [11C]raclopride using network control theory.

RESULTS

SD reduced dwell time and persistence probabilities, with the strongest effects in two brain states, one characterized by high default mode network and low dorsal attention network activity and the other by high frontoparietal network and low somatomotor network activity. Using network control theory, we showed that after SD, there was an overall increase in the control energy required for brain state transitions, with effects varying for different brain state transitions. Control energy requirement was negatively associated with transition probabilities under SD and restful wakefulness and accounted for SD-induced changes in transition probabilities. Alteration in the energy landscape was associated with SD-induced changes in striatal D2 receptor distribution.

CONCLUSIONS

These findings demonstrate altered occurrence of internally and externally oriented brain states following acute SD and suggest an association with energy requirements for brain state transitions modulated by striatal D2 receptors.

A network control theory pipeline for studying the dynamics of the structural connectome

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains the dynamics of a system. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter the dynamics of a system in a desired way. An interesting development for NCT in the neuroscience field is its application to study behavior and mental health symptoms. To date, NCT has been validated to study different aspects of the human structural connectome. NCT outputs can be monitored throughout developmental stages to study the effects of connectome topology on neural dynamics and, separately, to test the coherence of empirical datasets with brain function and stimulation. Here, we provide a comprehensive pipeline for applying NCT to structural connectomes by following two procedures. The main procedure focuses on computing the control energy associated with the transitions between specific neural activity states. The second procedure focuses on computing average controllability, which indexes nodes' general capacity to control the dynamics of the system. We provide recommendations for comparing NCT outputs against null network models, and we further support this approach with a Python-based software package called 'network control theory for python'. The procedures in this protocol are appropriate for users with a background in network neuroscience and experience in dynamical systems theory.

Imbalanced temporal states of cortical blood-oxygen-level-dependent signal variability during rest in episodic migraine

BACKGROUND

Migraine has been associated with functional brain changes including altered connectivity and activity both during and between headache attacks. Recent studies established that the variability of the blood-oxygen-level-dependent (BOLD) signal is an important attribute of brain activity, which has so far been understudied in migraine. In this study, we investigate how time-varying measures of BOLD variability change interictally in episodic migraine patients.

METHODS

Two independent resting state functional MRI datasets acquired on 3T (discovery cohort) and 1.5T MRI scanners (replication cohort) including 99 episodic migraine patients (n3T = 42, n1.5T=57) and 78 healthy controls (n3T = 46, n1.5T=32) were analyzed in this cross-sectional study. A framework using time-varying measures of BOLD variability was applied to derive BOLD variability states. Descriptors of BOLD variability states such as dwell time and fractional occupancy were calculated, then compared between migraine patients and healthy controls using Mann-Whitney U-tests. Spearman's rank correlation was calculated to test associations with clinical parameters.

RESULTS

Resting-state activity was characterized by states of high and low BOLD signal variability. Migraine patients in the discovery cohort spent more time in the low variability state (mean dwell time: p = 0.014, median dwell time: p = 0.022, maximum dwell time: p = 0.013, fractional occupancy: p = 0.013) and less time in the high variability state (mean dwell time: p = 0.021, median dwell time: p = 0.021, maximum dwell time: p = 0.025, fractional occupancy: p = 0.013). Higher uptime of the low variability state was associated with greater disability as measured by MIDAS scores (maximum dwell time: R = 0.45, p = 0.007; fractional occupancy: R = 0.36, p = 0.035). Similar results were observed in the replication cohort.

CONCLUSION

Episodic migraine patients spend more time in a state of low BOLD variability during rest in headache-free periods, which is associated with greater disability. BOLD variability states show potential as a replicable functional imaging marker in episodic migraine.

Alterations in Dynamic Functional Connectivity in Patients with Cerebral Small Vessel Disease

Cerebral small vessel disease (CSVD) is a common disease that seriously endangers people's health, and is easily overlooked by both patients and clinicians due to its near-silent onset. Dynamic functional connectivity (DFC) is a new concept focusing on the dynamic features and patterns of brain networks that represents a powerful tool for gaining novel insight into neurological diseases. To assess alterations in DFC in CSVD patients, and the correlation of DFC with cognitive function. We enrolled 35 CSVD patients and 31 normal control subjects (NC). Resting-state functional MRI (rs-fMRI) with a sliding-window approach and k-means clustering based on independent component analysis (ICA) was used to evaluate DFC. The temporal properties of fractional windows and the mean dwell time in each state, as well as the number of transitions between each pair of DFC states, were calculated. Additionally, we assessed the functional connectivity (FC) strength of the dynamic states and the associations of altered neuroimaging measures with cognitive performance. A dynamic analysis of all included subjects suggested four distinct functional connectivity states. Compared with the NC group, the CSVD group had more fractional windows and longer mean dwell times in state 4 characterized by sparse FC both inter-network and intra-networks. Additionally, the CSVD group had a reduced number of windows and shorter mean dwell times compared to the NC group in state 3 characterized by highly positive FC between the somatomotor and visual networks, and negative FC in the basal ganglia and somatomotor and visual networks. The number of transitions between state 2 and state 3 and between state 3 and state 4 was significantly reduced in the CSVD group compared to the NC group. Moreover, there was a significant difference in the FC strength between the two groups, and the altered temporal properties of DFC were significantly related to cognitive performance. Our study indicated that CSVD is characterized by altered temporal properties in DFC that may be sensitive neuroimaging biomarkers for early disease identification. Further study of DFC alterations could help us to better understand the progressive dysfunction of networks in CSVD patients.

Network enrichment significance testing in brain-phenotype association studies

Functional networks often guide our interpretation of spatial maps of brain-phenotype associations. However, methods for assessing enrichment of associations within networks of interest have varied in terms of both scientific rigor and underlying assumptions. While some approaches have relied on subjective interpretations, others have made unrealistic assumptions about spatial properties of imaging data, leading to inflated false positive rates. We seek to address this gap in existing methodology by borrowing insight from a method widely used in genetics research for testing enrichment of associations between a set of genes and a phenotype of interest. We propose network enrichment significance testing (NEST), a flexible framework for testing the specificity of brain-phenotype associations to functional networks or other sub-regions of the brain. We apply NEST to study enrichment of associations with structural and functional brain imaging data from a large-scale neurodevelopmental cohort study.

Electroencephalography Spectral-power Volatility Predicts Problem-solving Outcomes

Temporal variability is a fundamental property of brain processes and is functionally important to human cognition. This study examined how fluctuations in neural oscillatory activity are related to problem-solving performance as one example of how temporal variability affects high-level cognition. We used volatility to assess step-by-step fluctuations of EEG spectral power while individuals attempted to solve word-association puzzles. Inspired by recent results with hidden-state modeling, we tested the hypothesis that spectral-power volatility is directly associated with problem-solving outcomes. As predicted, volatility was lower during trials solved with insight compared with those solved analytically. Moreover, volatility during prestimulus preparation for problem-solving predicted solving outcomes, including solving success and solving time. These novel findings were replicated in a separate data set from an anagram-solving task, suggesting that less-rapid transitions between neural oscillatory synchronization and desynchronization predict better solving performance and are conducive to solving with insight for these types of problems. Thus, volatility can be a valuable index of cognition-related brain dynamics.

Disrupted brain state dynamics in opioid and alcohol use disorder: attenuation by nicotine use

Substance use disorder (SUD) is a chronic relapsing disorder with long-lasting changes in brain intrinsic networks. While most research to date has focused on static functional connectivity, less is known about the effect of chronic drug use on dynamics of brain networks. Here we investigated brain state dynamics in individuals with opioid use (OUD) and alcohol use disorder (AUD) and assessed how concomitant nicotine use, which is frequent among individuals with OUD and AUD, affects brain dynamics. Resting-state functional magnetic resonance imaging data of 27 OUD, 107 AUD, and 137 healthy participants were included in the analyses. To identify recurrent brain states and their dynamics, we applied a data-driven clustering approach that determines brain states at a single time frame. We found that OUD and AUD non-smokers displayed similar changes in brain state dynamics including decreased fractional occupancy or dwell time in default mode network (DMN)-dominated brain states and increased appearance rate in visual network (VIS)-dominated brain states, which were also reflected in transition probabilities of related brain states. Interestingly, co-use of nicotine affected brain states in an opposite manner by lowering VIS-dominated and enhancing DMN-dominated brain states in both OUD and AUD participants. Our finding revealed a similar pattern of brain state dynamics in OUD and AUD participants that differed from controls, with an opposite effect for nicotine use suggesting distinct effects of various drugs on brain state dynamics. Different strategies for treating SUD may need to be implemented based on patterns of co-morbid drug use.

Controlling the human connectome with spatially diffuse input signals

The human brain is never at "rest"; its activity is constantly fluctuating over time, transitioning from one brain state-a whole-brain pattern of activity-to another. Network control theory offers a framework for understanding the effort - energy - associated with these transitions. One branch of control theory that is especially useful in this context is "optimal control", in which input signals are used to selectively drive the brain into a target state. Typically, these inputs are introduced independently to the nodes of the network (each input signal is associated with exactly one node). Though convenient, this input strategy ignores the continuity of cerebral cortex - geometrically, each region is connected to its spatial neighbors, allowing control signals, both exogenous and endogenous, to spread from their foci to nearby regions. Additionally, the spatial specificity of brain stimulation techniques is limited, such that the effects of a perturbation are measurable in tissue surrounding the stimulation site. Here, we adapt the network control model so that input signals have a spatial extent that decays exponentially from the input site. We show that this more realistic strategy takes advantage of spatial dependencies in structural connectivity and activity to reduce the energy (effort) associated with brain state transitions. We further leverage these dependencies to explore near-optimal control strategies such that, on a per-transition basis, the number of input signals required for a given control task is reduced, in some cases by two orders of magnitude. This approximation yields network-wide maps of input site density, which we compare to an existing database of functional, metabolic, genetic, and neurochemical maps, finding a close correspondence. Ultimately, not only do we propose a more efficient framework that is also more adherent to well-established brain organizational principles, but we also posit neurobiologically grounded bases for optimal control.

Neuropsychology of chronic back pain managed with long-term opioid use

Chronic pain is commonly treated with long-term opioids, but the neuropsychological outcomes associated with stable long-duration opioid use remain unclear. Here, we contrasted the psychological profiles, brain activity, and brain structure of 70 chronic back pain patients on opioids (CBP+O, average opioid exposure 6.2 years) with 70 patients managing their pain without opioids. CBP+O exhibited moderately worse psychological profiles and small differences in brain morphology. However, CBP+O had starkly different spontaneous brain activity, dominated by increased mesocorticolimbic and decreased dorsolateral-prefrontal activity, even after controlling for pain intensity and duration. These differences strongly reflected cortical opioid and serotonin receptor densities and mapped to two antagonistic resting-state circuits. The circuits' dynamics were explained by mesocorticolimbic activity and reflected negative affect. We reassessed a sub-group of CBP+O after they briefly abstained from taking opioids. Network dynamics, but not spontaneous activity, reflected exacerbated signs of withdrawal. Our results have implications for the management and tapering of opioids in chronic pain.

Modular subgraphs in large-scale connectomes underpin spontaneous co-fluctuation events in mouse and human brains

Previous studies have adopted an edge-centric framework to study fine-scale network dynamics in human fMRI. To date, however, no studies have applied this framework to data collected from model organisms. Here, we analyze structural and functional imaging data from lightly anesthetized mice through an edge-centric lens. We find evidence of "bursty" dynamics and events - brief periods of high-amplitude network connectivity. Further, we show that on a per-frame basis events best explain static FC and can be divided into a series of hierarchically-related clusters. The co-fluctuation patterns associated with each cluster centroid link distinct anatomical areas and largely adhere to the boundaries of algorithmically detected functional brain systems. We then investigate the anatomical connectivity undergirding high-amplitude co-fluctuation patterns. We find that events induce modular bipartitions of the anatomical network of inter-areal axonal projections. Finally, we replicate these same findings in a human imaging dataset. In summary, this report recapitulates in a model organism many of the same phenomena observed in previously edge-centric analyses of human imaging data. However, unlike human subjects, the murine nervous system is amenable to invasive experimental perturbations. Thus, this study sets the stage for future investigation into the causal origins of fine-scale brain dynamics and high-amplitude co-fluctuations. Moreover, the cross-species consistency of the reported findings enhances the likelihood of future translation.

Altered structural connectivity and functional brain dynamics in individuals with heavy alcohol use

Heavy alcohol use and its associated conditions, such as alcohol use disorder (AUD), impact millions of individuals worldwide. While our understanding of the neurobiological correlates of AUD has evolved substantially, we still lack models incorporating whole-brain neuroanatomical, functional, and pharmacological information under one framework. Here, we utilize diffusion and functional magnetic resonance imaging to investigate alterations to brain dynamics in N = 130 individuals with a high amount of current alcohol use. We compared these alcohol using individuals to N = 308 individuals with minimal use of any substances. We find that individuals with heavy alcohol use had less dynamic and complex brain activity, and through leveraging network control theory, had increased control energy to complete transitions between activation states. Further, using separately acquired positron emission tomography (PET) data, we deploy an in silico evaluation demonstrating that decreased D2 receptor levels, as found previously in individuals with AUD, may relate to our observed findings. This work demonstrates that whole-brain, multimodal imaging information can be combined under a network control framework to identify and evaluate neurobiological correlates and mechanisms of AUD.

Hierarchical organization of spontaneous co-fluctuations in densely sampled individuals using fMRI

Edge time series decompose functional connectivity into its framewise contributions. Previous studies have focused on characterizing the properties of high-amplitude frames (time points when the global co-fluctuation amplitude takes on its largest value), including their cluster structure. Less is known about middle- and low-amplitude co-fluctuations (peaks in co-fluctuation time series but of lower amplitude). Here, we directly address those questions, using data from two dense-sampling studies: the MyConnectome project and Midnight Scan Club. We develop a hierarchical clustering algorithm to group peak co-fluctuations of all magnitudes into nested and multiscale clusters based on their pairwise concordance. At a coarse scale, we find evidence of three large clusters that, collectively, engage virtually all canonical brain systems. At finer scales, however, each cluster is dissolved, giving way to increasingly refined patterns of co-fluctuations involving specific sets of brain systems. We also find an increase in global co-fluctuation magnitude with hierarchical scale. Finally, we comment on the amount of data needed to estimate co-fluctuation pattern clusters and implications for brain-behavior studies. Collectively, the findings reported here fill several gaps in current knowledge concerning the heterogeneity and richness of co-fluctuation patterns as estimated with edge time series while providing some practical guidance for future studies.

Using network control theory to study the dynamics of the structural connectome

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains dynamics. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter dynamics in a desired way. We have extensively developed and validated the application of NCT to the human structural connectome. Through these efforts, we have studied (i) how different aspects of connectome topology affect neural dynamics, (ii) whether NCT outputs cohere with empirical data on brain function and stimulation, and (iii) how NCT outputs vary across development and correlate with behavior and mental health symptoms. In this protocol, we introduce a framework for applying NCT to structural connectomes following two main pathways. Our primary pathway focuses on computing the control energy associated with transitioning between specific neural activity states. Our second pathway focuses on computing average controllability, which indexes nodes' general capacity to control dynamics. We also provide recommendations for comparing NCT outputs against null network models. Finally, we support this protocol with a Python-based software package called network control theory for python (nctpy).

Alprazolam modulates persistence energy during emotion processing in first-degree relatives of individuals with schizophrenia: a network control study

Schizophrenia is marked by deficits in facial affect processing associated with abnormalities in GABAergic circuitry, deficits also found in first-degree relatives. Facial affect processing involves a distributed network of brain regions including limbic regions like amygdala and visual processing areas like fusiform cortex. Pharmacological modulation of GABAergic circuitry using benzodiazepines like alprazolam can be useful for studying this facial affect processing network and associated GABAergic abnormalities in schizophrenia. Here, we use pharmacological modulation and computational modeling to study the contribution of GABAergic abnormalities toward emotion processing deficits in schizophrenia. Specifically, we apply principles from network control theory to model persistence energy - the control energy required to maintain brain activation states - during emotion identification and recall tasks, with and without administration of alprazolam, in a sample of first-degree relatives and healthy controls. Here, persistence energy quantifies the magnitude of theoretical external inputs during the task. We find that alprazolam increases persistence energy in relatives but not in controls during threatening face processing, suggesting a compensatory mechanism given the relative absence of behavioral abnormalities in this sample of unaffected relatives. Further, we demonstrate that regions in the fusiform and occipital cortices are important for facilitating state transitions during facial affect processing. Finally, we uncover spatial relationships (i) between regional variation in differential control energy (alprazolam versus placebo) and (ii) both serotonin and dopamine neurotransmitter systems, indicating that alprazolam may exert its effects by altering neuromodulatory systems. Together, these findings provide a new perspective on the distributed emotion processing network and the effect of GABAergic modulation on this network, in addition to identifying an association between schizophrenia risk and abnormal GABAergic effects on persistence energy during threat processing.

The Cellular Underpinnings of the Human Cortical Connectome

The functional properties of the human brain arise, in part, from the vast assortment of cell types that pattern the cortex. The cortical sheet can be broadly divided into distinct networks, which are further embedded into processing streams, or gradients, that extend from unimodal systems through higher-order association territories. Here, using transcriptional data from the Allen Human Brain Atlas, we demonstrate that imputed cell type distributions are spatially coupled to the functional organization of cortex, as estimated through fMRI. Cortical cellular profiles follow the macro-scale organization of the functional gradients as well as the associated large-scale networks. Distinct cellular fingerprints were evident across networks, and a classifier trained on post-mortem cell-type distributions was able to predict the functional network allegiance of cortical tissue samples. These data indicate that the in vivo organization of the cortical sheet is reflected in the spatial variability of its cellular composition.

Larger lesion volume in people with multiple sclerosis is associated with increased transition energies between brain states and decreased entropy of brain activity

Quantifying the relationship between the brain's functional activity patterns and its structural backbone is crucial when relating the severity of brain pathology to disability in multiple sclerosis (MS). Network control theory (NCT) characterizes the brain's energetic landscape using the structural connectome and patterns of brain activity over time. We applied NCT to investigate brain-state dynamics and energy landscapes in controls and people with MS (pwMS). We also computed entropy of brain activity and investigated its association with the dynamic landscape's transition energy and lesion volume. Brain states were identified by clustering regional brain activity vectors, and NCT was applied to compute the energy required to transition between these brain states. We found that entropy was negatively correlated with lesion volume and transition energy, and that larger transition energies were associated with pwMS with disability. This work supports the notion that shifts in the pattern of brain activity in pwMS without disability results in decreased transition energies compared to controls, but, as this shift evolves over the disease, transition energies increase beyond controls and disability occurs. Our results provide the first evidence in pwMS that larger lesion volumes result in greater transition energy between brain states and decreased entropy of brain activity.

Effect of LSD and music on the time-varying brain dynamics

RATIONALE

Psychedelics are getting closer to being widely used in clinical treatment. Music is known as a key element of psychedelic-assisted therapy due to its psychological effects, specifically on the emotion, meaning-making, and sensory processing. However, there is still a lack of understanding in how psychedelics influence brain activity in experimental settings involving music listening.

OBJECTIVES

The main goal of our research was to investigate the effect of music, as a part of "setting," on the brain states dynamics after lysergic acid diethylamide (LSD) intake.

METHODS

We used an open dataset, where a group of 15 participants underwent two functional MRI scanning sessions under LSD and placebo influence. Every scanning session contained three runs: two resting-state runs separated by one run with music listening. We applied K-Means clustering to identify the repetitive patterns of brain activity, so-called brain states. For further analysis, we calculated states' dwell time, fractional occupancy and transition probability.

RESULTS

The interaction effect of music and psychedelics led to change in the time-varying brain activity of the task-positive state. LSD, regardless of the music, affected the dynamics of the state of combined activity of DMN, SOM, and VIS networks. Crucially, we observed that the music itself could potentially have a long-term influence on the resting-state, in particular on states involving task-positive networks.

CONCLUSIONS

This study indicates that music, as a crucial element of "setting," can potentially have an influence on the subject's resting-state during psychedelic experience. Further studies should replicate these results on a larger sample size.

Time-resolved network control analysis links reduced control energy under DMT with the serotonin 2a receptor, signal diversity, and subjective experience

Psychedelics offer a profound window into the functioning of the human brain and mind through their robust acute effects on perception, subjective experience, and brain activity patterns. In recent work using a receptor-informed network control theory framework, we demonstrated that the serotonergic psychedelics lysergic acid diethylamide (LSD) and psilocybin flatten the brain's control energy landscape in a manner that covaries with more dynamic and entropic brain activity. Contrary to LSD and psilocybin, whose effects last for hours, the serotonergic psychedelic N,N-dimethyltryptamine (DMT) rapidly induces a profoundly immersive altered state of consciousness lasting less than 20 minutes, allowing for the entirety of the drug experience to be captured during a single resting-state fMRI scan. Using network control theory, which quantifies the amount of input necessary to drive transitions between functional brain states, we integrate brain structure and function to map the energy trajectories of 14 individuals undergoing fMRI during DMT and placebo. Consistent with previous work, we find that global control energy is reduced following injection with DMT compared to placebo. We additionally show longitudinal trajectories of global control energy correlate with longitudinal trajectories of EEG signal diversity (a measure of entropy) and subjective ratings of drug intensity. We interrogate these same relationships on a regional level and find that the spatial patterns of DMT's effects on these metrics are correlated with serotonin 2a receptor density (obtained from separately acquired PET data). Using receptor distribution and pharmacokinetic information, we were able to successfully recapitulate the effects of DMT on global control energy trajectories, demonstrating a proof-of-concept for the use of control models in predicting pharmacological intervention effects on brain dynamics.

Transitions between cognitive topographies: contributions of network structure, neuromodulation, and disease

Patterns of neural activity underlie human cognition. Transitions between these patterns are orchestrated by the brain's network architecture. What are the mechanisms linking network structure to cognitively relevant activation patterns? Here we implement principles of network control to investigate how the architecture of the human connectome shapes transitions between 123 experimentally defined cognitive activation maps (cognitive topographies) from the NeuroSynth meta-analytic engine. We also systematically incorporate neurotransmitter receptor density maps (18 receptors and transporters) and disease-related cortical abnormality maps (11 neurodegenerative, psychiatric and neurodevelopmental diseases; N = 17 000 patients, N = 22 000 controls). Integrating large-scale multimodal neuroimaging data from functional MRI, diffusion tractography, cortical morphometry, and positron emission tomography, we simulate how anatomically-guided transitions between cognitive states can be reshaped by pharmacological or pathological perturbation. Our results provide a comprehensive look-up table charting how brain network organisation and chemoarchitecture interact to manifest different cognitive topographies. This computational framework establishes a principled foundation for systematically identifying novel ways to promote selective transitions between desired cognitive topographies.

Task-based differences in brain state dynamics and their relation to cognitive ability

Transient patterns of interregional connectivity form and dissipate in response to varying cognitive demands. Yet, it is not clear how different cognitive demands influence brain state dynamics, and whether these dynamics relate to general cognitive ability. Here, using functional magnetic resonance imaging (fMRI) data, we characterised shared, recurrent, global brain states in 187 participants across the working memory, emotion, language, and relation tasks from the Human Connectome Project. Brain states were determined using Leading Eigenvector Dynamics Analysis (LEiDA). In addition to the LEiDA-based metrics of brain state lifetimes and probabilities, we also computed information-theoretic measures of Block Decomposition Method of complexity, Lempel-Ziv complexity and transition entropy. Information theoretic metrics are notable in their ability to compute relationships amongst sequences of states over time, compared to lifetime and probability, which capture the behaviour of each state in isolation. We then related task-based brain state metrics to fluid intelligence. We observed that brain states exhibited stable topology across a range of numbers of clusters (K = 2:15). Most metrics of brain state dynamics, including state lifetime, probability, and all information theoretic metrics, reliably differed between tasks. However, relationships between state dynamic metrics and cognitive abilities varied according to the task, the metric, and the value of K, indicating that there are contextual relationships between task-dependant state dynamics and trait cognitive ability. This study provides evidence that the brain reconfigures across time in response to cognitive demands, and that there are contextual, rather than generalisable, relationships amongst task, state dynamics, and cognitive ability.

Toward Best Practices for Imaging Transcriptomics of the Human Brain

Modern brainwide transcriptional atlases provide unprecedented opportunities for investigating the molecular correlates of brain organization, as quantified using noninvasive neuroimaging. However, integrating neuroimaging data with transcriptomic measures is not straightforward, and careful consideration is required to make valid inferences. In this article, we review recent work exploring how various methodological choices affect 3 main phases of imaging transcriptomic analyses, including 1) processing of transcriptional atlas data; 2) relating transcriptional measures to independently derived neuroimaging phenotypes; and 3) evaluating the functional implications of identified associations through gene enrichment analyses. Our aim is to facilitate the development of standardized and reproducible approaches for this rapidly growing field. We identify sources of methodological variability, key choices that can affect findings, and considerations for mitigating false positive and/or spurious results. Finally, we provide an overview of freely available open-source toolboxes implementing current best-practice procedures across all 3 analysis phases.

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.

Mindful attention promotes control of brain network dynamics for self-regulation and discontinues the past from the present

Mindful attention is characterized by acknowledging the present experience as a transient mental event. Early stages of mindfulness practice may require greater neural effort for later efficiency. Early effort may self-regulate behavior and focalize the present, but this understanding lacks a computational explanation. Here we used network control theory as a model of how external control inputs-operationalizing effort-distribute changes in neural activity evoked during mindful attention across the white matter network. We hypothesized that individuals with greater network controllability, thereby efficiently distributing control inputs, effectively self-regulate behavior. We further hypothesized that brain regions that utilize greater control input exhibit shorter intrinsic timescales of neural activity. Shorter timescales characterize quickly discontinuing past processing to focalize the present. We tested these hypotheses in a randomized controlled study that primed participants to either mindfully respond or naturally react to alcohol cues during fMRI and administered text reminders and measurements of alcohol consumption during 4 wk postscan. We found that participants with greater network controllability moderated alcohol consumption. Mindful regulation of alcohol cues, compared to one's own natural reactions, reduced craving, but craving did not differ from the baseline group. Mindful regulation of alcohol cues, compared to the natural reactions of the baseline group, involved more-effortful control of neural dynamics across cognitive control and attention subnetworks. This effort persisted in the natural reactions of the mindful group compared to the baseline group. More-effortful neural states had shorter timescales than less effortful states, offering an explanation for how mindful attention promotes being present.

Asymmetric signaling across the hierarchy of cytoarchitecture within the human connectome

Cortical variations in cytoarchitecture form a sensory-fugal axis that shapes regional profiles of extrinsic connectivity and is thought to guide signal propagation and integration across the cortical hierarchy. While neuroimaging work has shown that this axis constrains local properties of the human connectome, it remains unclear whether it also shapes the asymmetric signaling that arises from higher-order topology. Here, we used network control theory to examine the amount of energy required to propagate dynamics across the sensory-fugal axis. Our results revealed an asymmetry in this energy, indicating that bottom-up transitions were easier to complete compared to top-down. Supporting analyses demonstrated that asymmetries were underpinned by a connectome topology that is wired to support efficient bottom-up signaling. Lastly, we found that asymmetries correlated with differences in communicability and intrinsic neuronal time scales and lessened throughout youth. Our results show that cortical variation in cytoarchitecture may guide the formation of macroscopic connectome topology.

Uncovering the biological basis of control energy: Structural and metabolic correlates of energy inefficiency in temporal lobe epilepsy

Network control theory is increasingly used to profile the brain's energy landscape via simulations of neural dynamics. This approach estimates the control energy required to simulate the activation of brain circuits based on structural connectome measured using diffusion magnetic resonance imaging, thereby quantifying those circuits' energetic efficiency. The biological basis of control energy, however, remains unknown, hampering its further application. To fill this gap, investigating temporal lobe epilepsy as a lesion model, we show that patients require higher control energy to activate the limbic network than healthy volunteers, especially ipsilateral to the seizure focus. The energetic imbalance between ipsilateral and contralateral temporolimbic regions is tracked by asymmetric patterns of glucose metabolism measured using positron emission tomography, which, in turn, may be selectively explained by asymmetric gray matter loss as evidenced in the hippocampus. Our investigation provides the first theoretical framework unifying gray matter integrity, metabolism, and energetic generation of neural dynamics.

neuromaps: structural and functional interpretation of brain maps

Imaging technologies are increasingly used to generate high-resolution reference maps of brain structure and function. Comparing experimentally generated maps to these reference maps facilitates cross-disciplinary scientific discovery. Although recent data sharing initiatives increase the accessibility of brain maps, data are often shared in disparate coordinate systems, precluding systematic and accurate comparisons. Here we introduce neuromaps, a toolbox for accessing, transforming and analyzing structural and functional brain annotations. We implement functionalities for generating high-quality transformations between four standard coordinate systems. The toolbox includes curated reference maps and biological ontologies of the human brain, such as molecular, microstructural, electrophysiological, developmental and functional ontologies. Robust quantitative assessment of map-to-map similarity is enabled via a suite of spatial autocorrelation-preserving null models. neuromaps combines open-access data with transparent functionality for standardizing and comparing brain maps, providing a systematic workflow for comprehensive structural and functional annotation enrichment analysis of the human brain.

Cofluctuation analysis reveals aberrant default mode network patterns in adolescents and youths with autism spectrum disorder

Resting-state functional connectivity (rsFC) approaches provide informative estimates of the functional architecture of the brain, and recently-proposed cofluctuation analysis temporally unwraps FC at every moment in time, providing refined information for quantifying brain dynamics. As a brain network disorder, autism spectrum disorder (ASD) was characterized by substantial alteration in FC, but the contribution of moment-to-moment-activity cofluctuations to the overall dysfunctional connectivity pattern in ASD remains poorly understood. Here, we used the cofluctuation approach to explore the underlying dynamic properties of FC in ASD, using a large multisite resting-state functional magnetic resonance imaging (rs-fMRI) dataset (ASD = 354, typically developing controls [TD] = 446). Our results verified that the networks estimated using high-amplitude frames were highly correlated with the traditional rsFC. Moreover, these frames showed higher average amplitudes in participants with ASD than those in the TD group. Principal component analysis was performed on the activity patterns in these frames and aggregated over all subjects. The first principal component (PC1) corresponds to the default mode network (DMN), and the PC1 coefficients were greater in participants with ASD than those in the TD group. Additionally, increased ASD symptom severity was associated with the increased coefficients, which may result in excessive internally oriented cognition and social cognition deficits in individuals with ASD. Our finding highlights the utility of cofluctuation approaches in prevalent neurodevelopmental disorders and verifies that the aberrant contribution of DMN to rsFC may underline the symptomatology in adolescents and youths with ASD.

Receptor-informed network control theory links LSD and psilocybin to a flattening of the brain's control energy landscape

Psychedelics including lysergic acid diethylamide (LSD) and psilocybin temporarily alter subjective experience through their neurochemical effects. Serotonin 2a (5-HT2a) receptor agonism by these compounds is associated with more diverse (entropic) brain activity. We postulate that this increase in entropy may arise in part from a flattening of the brain's control energy landscape, which can be observed using network control theory to quantify the energy required to transition between recurrent brain states. Using brain states derived from existing functional magnetic resonance imaging (fMRI) datasets, we show that LSD and psilocybin reduce control energy required for brain state transitions compared to placebo. Furthermore, across individuals, reduction in control energy correlates with more frequent state transitions and increased entropy of brain state dynamics. Through network control analysis that incorporates the spatial distribution of 5-HT2a receptors (obtained from publicly available positron emission tomography (PET) data under non-drug conditions), we demonstrate an association between the 5-HT2a receptor and reduced control energy. Our findings provide evidence that 5-HT2a receptor agonist compounds allow for more facile state transitions and more temporally diverse brain activity. More broadly, we demonstrate that receptor-informed network control theory can model the impact of neuropharmacological manipulation on brain activity dynamics.

Equalization of Brain State Occupancy Accompanies Cognitive Impairment in Cerebral Small Vessel Disease

BACKGROUND

Cognitive impairment is a hallmark of cerebral small vessel disease (cSVD). Functional magnetic resonance imaging has highlighted connections between patterns of brain activity and variability in behavior. We aimed to characterize the associations between imaging markers of cSVD, dynamic connectivity, and cognitive impairment.

METHODS

We obtained magnetic resonance imaging and clinical data from the population-based Hamburg City Health Study. cSVD was quantified by white matter hyperintensities and peak-width of skeletonized mean diffusivity (PSMD). Resting-state blood oxygen level-dependent signals were clustered into discrete brain states, for which fractional occupancies (%) and dwell times (seconds) were computed. Cognition in multiple domains was assessed using validated tests. Regression analysis was used to quantify associations between white matter damage, spatial coactivation patterns, and cognitive function.

RESULTS

Data were available for 979 participants (ages 45-74 years, median white matter hyperintensity volume 0.96 mL). Clustering identified five brain states with the most time spent in states characterized by activation (+) or suppression (-) of the default mode network (DMN) (fractional occupancy: DMN+ = 25.1 ± 7.2%, DMN- = 25.5 ± 7.2%). Every 4.7-fold increase in white matter hyperintensity volume was associated with a 0.95-times reduction of the odds of occupying DMN+ or DMN-. Time spent in DMN-related brain states was associated with executive function.

CONCLUSIONS

Associations between white matter damage, whole-brain spatial coactivation patterns, and cognition suggest equalization of time spent in different brain states as a marker for cSVD-associated cognitive decline. Reduced gradients between brain states in association with brain damage and cognitive impairment reflect the dedifferentiation hypothesis of neurocognitive aging in a network-theoretical context.

Frontoparietal network activation is associated with motor recovery in ischemic stroke patients

Strokes cause lesions that damage brain tissue, disrupt normal brain activity patterns and can lead to impairments in motor function. Although modulation of cortical activity is central to stimulation-based rehabilitative therapies, aberrant and adaptive patterns of brain activity after stroke have not yet been fully characterized. Here, we apply a brain dynamics analysis approach to study longitudinal brain activity patterns in individuals with ischemic pontine stroke. We first found 4 commonly occurring brain states largely characterized by high amplitude activations in the visual, frontoparietal, default mode, and motor networks. Stroke subjects spent less time in the frontoparietal state compared to controls. For individuals with dominant-hand CST damage, more time spent in the frontoparietal state from 1 week to 3-6 months post-stroke was associated with better motor recovery over the same time period, an association which was independent of baseline impairment. Furthermore, the amount of time spent in brain states was linked empirically to functional connectivity. This work suggests that when the dominant-hand CST is compromised in stroke, resting state configurations may include increased activation of the frontoparietal network, which may facilitate compensatory neural pathways that support recovery of motor function when traditional motor circuits of the dominant-hemisphere are compromised.

External drivers of BOLD signal's non-stationarity

A fundamental challenge in neuroscience is to uncover the principles governing how the brain interacts with the external environment. However, assumptions about external stimuli fundamentally constrain current computational models. We show in silico that unknown external stimulation can produce error in the estimated linear time-invariant dynamical system. To address these limitations, we propose an approach to retrieve the external (unknown) input parameters and demonstrate that the estimated system parameters during external input quiescence uncover spatiotemporal profiles of external inputs over external stimulation periods more accurately. Finally, we unveil the expected (and unexpected) sensory and task-related extra-cortical input profiles using functional magnetic resonance imaging data acquired from 96 subjects (Human Connectome Project) during the resting-state and task scans. This dynamical systems model of the brain offers information on the structure and dimensionality of the BOLD signal's external drivers and shines a light on the likely external sources contributing to the BOLD signal's non-stationarity. Our findings show the role of exogenous inputs in the BOLD dynamics and highlight the importance of accounting for external inputs to unravel the brain's time-varying functional dynamics.

Changes in dynamic transitions between integrated and segregated states underlie visual hallucinations in Parkinson's disease

Hallucinations are a core feature of psychosis and common in Parkinson's. Their transient, unexpected nature suggests a change in dynamic brain states, but underlying causes are unknown. Here, we examine temporal dynamics and underlying structural connectivity in Parkinson's-hallucinations using a combination of functional and structural MRI, network control theory, neurotransmitter density and genetic analyses. We show that Parkinson's-hallucinators spent more time in a predominantly Segregated functional state with fewer between-state transitions. The transition from integrated-to-segregated state had lower energy cost in Parkinson's-hallucinators; and was therefore potentially preferable. The regional energy needed for this transition was correlated with regional neurotransmitter density and gene expression for serotoninergic, GABAergic, noradrenergic and cholinergic, but not dopaminergic, receptors. We show how the combination of neurochemistry and brain structure jointly shape functional brain dynamics leading to hallucinations and highlight potential therapeutic targets by linking these changes to neurotransmitter systems involved in early sensory and complex visual processing.

Predicting Spatio-Temporal Human Brain Response Using fMRI

The transformation and transmission of brain stimuli reflect the dynamical brain activity in space and time. Compared with functional magnetic resonance imaging (fMRI), magneto- or electroencephalography (M/EEG) fast couples to the neural activity through generated magnetic fields. However, the MEG signal is inhomogeneous throughout the whole brain, which is affected by the signal-to-noise ratio, the sensors' location and distance. Current non-invasive neuroimaging modalities such as fMRI and M/EEG excel high resolution in space or time but not in both. To solve the main limitations of current technique for brain activity recording, we propose a novel recurrent memory optimization approach to predict the internal behavioral states in space and time. The proposed method uses Optimal Polynomial Projections to capture the long temporal history with robust online compression. The training process takes the pairs of fMRI and MEG data as inputs and predicts the recurrent brain states through the Siamese network. In the testing process, the framework only uses fMRI data to generate the corresponding neural response in space and time. The experimental results with Human connectome project (HCP) show that the predicted signal could reflect the neural activity with high spatial resolution as fMRI and high temporal resolution as MEG signal. The experimental results demonstrate for the first time that the proposed method is able to predict the brain response in both milliseconds and millimeters using only fMRI signal.

The expanding horizons of network neuroscience: From description to prediction and control

The field of network neuroscience has emerged as a natural framework for the study of the brain and has been increasingly applied across divergent problems in neuroscience. From a disciplinary perspective, network neuroscience originally emerged as a formal integration of graph theory (from mathematics) and neuroscience (from biology). This early integration afforded marked utility in describing the interconnected nature of neural units, both structurally and functionally, and underscored the relevance of that interconnection for cognition and behavior. But since its inception, the field has not remained static in its methodological composition. Instead, it has grown to use increasingly advanced graph-theoretic tools and to bring in several other disciplinary perspectives-including machine learning and systems engineering-that have proven complementary. In doing so, the problem space amenable to the discipline has expanded markedly. In this review, we discuss three distinct flavors of investigation in state-of-the-art network neuroscience: (i) descriptive network neuroscience, (ii) predictive network neuroscience, and (iii) a perturbative network neuroscience that draws on recent advances in network control theory. In considering each area, we provide a brief summary of the approaches, discuss the nature of the insights obtained, and highlight future directions.

Development of Thalamocortical Structural Connectivity in Typically Developing and Psychosis Spectrum Youths

BACKGROUND

Thalamocortical white matter connectivity is disrupted in psychosis and is hypothesized to play a role in its etiology and associated cognitive impairment. Attenuated cognitive symptoms often begin in adolescence, during a critical phase of white matter and cognitive development. However, little is known about the development of thalamocortical white matter connectivity and its association with cognition.

METHODS

This study characterized effects of age, sex, psychosis symptomatology, and cognition in thalamocortical networks in a large sample of youths (N = 1144, ages 8-22 years, 46% male) from the Philadelphia Neurodevelopmental Cohort, which included 316 typically developing youths, 330 youths on the psychosis spectrum, and 498 youths with other psychopathology. Probabilistic tractography was used to quantify percent total connectivity between the thalamus and six cortical regions and assess microstructural properties (i.e., fractional anisotropy) of thalamocortical white matter tracts.

RESULTS

Overall, percent total connectivity of the thalamus was weakly associated with age and was not associated with psychopathology or cognition. In contrast, fractional anisotropy of all thalamocortical tracts increased significantly with age, was generally higher in males than females, and was lowest in youths on the psychosis spectrum. Fractional anisotropy of tracts linking the thalamus to prefrontal and posterior parietal cortices was related to better cognitive function across subjects.

CONCLUSIONS

By characterizing the pattern of typical development and alterations in those at risk for psychotic disorders, this study provides a foundation for further conceptualization of thalamocortical white matter microstructure as a marker of neurodevelopment supporting cognition and an important risk marker for psychosis.

Null models in network neuroscience

Recent advances in imaging and tracing technology provide increasingly detailed reconstructions of brain connectomes. Concomitant analytic advances enable rigorous identification and quantification of functionally important features of brain network architecture. Null models are a flexible tool to statistically benchmark the presence or magnitude of features of interest, by selectively preserving specific architectural properties of brain networks while systematically randomizing others. Here we describe the logic, implementation and interpretation of null models of connectomes. We introduce randomization and generative approaches to constructing null networks, and outline a taxonomy of network methods for statistical inference. We highlight the spectrum of null models - from liberal models that control few network properties, to conservative models that recapitulate multiple properties of empirical networks - that allow us to operationalize and test detailed hypotheses about the structure and function of brain networks. We review emerging scenarios for the application of null models in network neuroscience, including for spatially embedded networks, annotated networks and correlation-derived networks. Finally, we consider the limits of null models, as well as outstanding questions for the field.

Comparison of Resting-State Functional MRI Methods for Characterizing Brain Dynamics

Resting-state functional MRI (fMRI) exhibits time-varying patterns of functional connectivity. Several different analysis approaches have been developed for examining these resting-state dynamics including sliding window connectivity (SWC), phase synchrony (PS), co-activation pattern (CAP), and quasi-periodic patterns (QPP). Each of these approaches can be used to generate patterns of activity or inter-areal coordination which vary across time. The individual frames can then be clustered to produce temporal groupings commonly referred to as “brain states.” Several recent publications have investigated brain state alterations in clinical populations, typically using a single method for quantifying frame-wise functional connectivity. This study directly compares the results of k-means clustering in conjunction with three of these resting-state dynamics methods (SWC, CAP, and PS) and quantifies the brain state dynamics across several metrics using high resolution data from the human connectome project. Additionally, these three dynamics methods are compared by examining how the brain state characterizations vary during the repeated sequences of brain states identified by a fourth dynamic analysis method, QPP. The results indicate that the SWC, PS, and CAP methods differ in the clusters and trajectories they produce. A clear illustration of these differences is given by how each one results in a very different clustering profile for the 24s sequences explicitly identified by the QPP algorithm. PS clustering is sensitive to QPPs with the mid-point of most QPP sequences grouped into the same single cluster. CAPs are also highly sensitive to QPPs, separating each phase of the QPP sequences into different sets of clusters. SWC (60s window) is less sensitive to QPPs. While the QPPs are slightly more likely to occur during specific SWC clusters, the SWC clustering does not vary during the 24s QPP sequences, the goal of this work is to improve both the practical and theoretical understanding of different resting-state dynamics methods, thereby enabling investigators to better conceptualize and implement these tools for characterizing functional brain networks.

Control theory illustrates the energy efficiency in the dynamic reconfiguration of functional connectivity

The brain's functional connectivity fluctuates over time instead of remaining steady in a stationary mode even during the resting state. This fluctuation establishes the dynamical functional connectivity that transitions in a non-random order between multiple modes. Yet it remains unexplored how the transition facilitates the entire brain network as a dynamical system and what utility this mechanism for dynamic reconfiguration can bring over the widely used graph theoretical measurements. To address these questions, we propose to conduct an energetic analysis of functional brain networks using resting-state fMRI and behavioral measurements from the Human Connectome Project. Through comparing the state transition energy under distinct adjacent matrices, we justify that dynamic functional connectivity leads to 60% less energy cost to support the resting state dynamics than static connectivity when driving the transition through default mode network. Moreover, we demonstrate that combining graph theoretical measurements and our energy-based control measurements as the feature vector can provide complementary prediction power for the behavioral scores (Combination vs. Control: t = 9.41, p = 1.64e-13; Combination vs. Graph: t = 4.92, p = 3.81e-6). Our approach integrates statistical inference and dynamical system inspection towards understanding brain networks.

Simulated visual hallucinations in virtual reality enhance cognitive flexibility

Historically, psychedelic drugs are known to modulate cognitive flexibility, a central aspect of cognition permitting adaptation to changing environmental demands. Despite proof suggesting phenomenological similarities between artificially-induced and actual psychedelic altered perception, experimental evidence is still lacking about whether the former is also able to modulate cognitive flexibility. To address this, we measure participants' cognitive flexibility through behavioral tasks after the exposure to virtual reality panoramic videos and their hallucinatory-like counterparts generated by the DeepDream algorithm. Results show that the estimated semantic network has a flexible structure when preceded by altered videos. Crucially, following the simulated psychedelic exposure, individuals also show an attenuated contribution of the automatic process and chaotic dynamics underlying the decision process. This suggests that simulated altered perceptual phenomenology enhances cognitive flexibility, presumably due to a reorganization in the cognitive dynamics that facilitates the exploration of uncommon decision strategies and inhibits automated choices.

Unique spatiotemporal fMRI dynamics in the awake mouse brain

Human imaging studies have shown that spontaneous brain activity exhibits stereotypic spatiotemporal reorganization in awake, conscious conditions with respect to minimally conscious states. However, whether and how this phenomenon can be generalized to lower mammalian species remains unclear. Leveraging a robust protocol for resting-state fMRI (rsfMRI) mapping in non-anesthetized, head-fixed mice, we investigated functional network topography and dynamic structure of spontaneous brain activity in wakeful animals. We found that rsfMRI networks in the awake state, while anatomically comparable to those observed under anesthesia, are topologically configured to maximize interregional communication, departing from the underlying community structure of the mouse axonal connectome. We further report that rsfMRI activity in wakeful animals exhibits unique spatiotemporal dynamics characterized by a state-dependent, dominant occurrence of coactivation patterns encompassing a prominent participation of arousal-related forebrain nuclei and functional anti-coordination between visual-auditory and polymodal cortical areas. We finally show that rsfMRI dynamics in awake mice exhibits a stereotypical temporal structure, in which state-dominant coactivation patterns are configured as network attractors. These findings suggest that spontaneous brain activity in awake mice is critically shaped by state-specific involvement of basal forebrain arousal systems and document that its dynamic structure recapitulates distinctive, evolutionarily relevant principles that are predictive of conscious states in higher mammalian species.

•

fMRI networks in awake mice depart from underlying anatomical structure

•

fMRI dynamics in wakeful mice is critically shaped by arousal-related nuclei

•

Occurrence and topography of rsfMRI coactivation patterns define conscious states

•

fMRI coactivation dynamics defines a signature of consciousness in the mouse brain

Network controllability mediates the relationship between rigid structure and flexible dynamics

Precisely how the anatomical structure of the brain supports a wide range of complex functions remains a question of marked importance in both basic and clinical neuroscience. Progress has been hampered by the lack of theoretical frameworks explaining how a structural network of relatively rigid interareal connections can produce a diverse repertoire of functional neural dynamics. Here, we address this gap by positing that the brain's structural network architecture determines the set of accessible functional connectivity patterns according to predictions of network control theory. In a large developmental cohort of 823 youths aged 8 to 23 years, we found that the flexibility of a brain region's functional connectivity was positively correlated with the proportion of its structural links extending to different cognitive systems. Notably, this relationship was mediated by nodes' boundary controllability, suggesting that a region's strategic location on the boundaries of modules may underpin the capacity to integrate information across different cognitive processes. Broadly, our study provides a mechanistic framework that illustrates how temporal flexibility observed in functional networks may be mediated by the controllability of the underlying structural connectivity.