Impact of white matter hyperintensities on structural connectivity and cognition in cognitively intact ADNI participants

We used indirect brain mapping with virtual lesion tractography to test the hypothesis that the extent of white matter tract disconnection due to white matter hyperintensities (WMH) is associated with corresponding tract-specific cognitive performance decrements. To estimate tract disconnection, WMH masks were extracted from FLAIR MRI data of 481 cognitively intact participants in the Alzheimer's Disease Neuroimaging Initiative (ADNI) and used as regions of avoidance for fiber tracking in diffusion MRI data from 50 healthy young participants from the Human Connectome Project. Estimated tract disconnection in the right inferior fronto-occipital fasciculus, right frontal aslant tract, and right superior longitudinal fasciculus mediated the effects of WMH volume on executive function. Estimated tract disconnection in the left uncinate fasciculus mediated the effects of WMH volume on memory and in the right frontal aslant tract on language. In a subset of ADNI control participants with amyloid data, positive status increased the probability of periventricular WMH and moderated the relationship between WMH burden and tract disconnection in executive function performance.

Measuring neuronal avalanches to inform brain-computer interfaces

Large-scale interactions among multiple brain regions manifest as bursts of activations called neuronal avalanches, which reconfigure according to the task at hand and, hence, might constitute natural candidates to design brain-computer interfaces (BCIs). To test this hypothesis, we used source-reconstructed magneto/electroencephalography during resting state and a motor imagery task performed within a BCI protocol. To track the probability that an avalanche would spread across any two regions, we built an avalanche transition matrix (ATM) and demonstrated that the edges whose transition probabilities significantly differed between conditions hinged selectively on premotor regions in all subjects. Furthermore, we showed that the topology of the ATMs allows task-decoding above the current gold standard. Hence, our results suggest that neuronal avalanches might capture interpretable differences between tasks that can be used to inform brain-computer interfaces.

XCP-D: A Robust Pipeline for the post-processing of fMRI data

Functional neuroimaging is an essential tool for neuroscience research. Pre-processing pipelines produce standardized, minimally pre-processed data to support a range of potential analyses. However, post-processing is not similarly standardized. While several options for post-processing exist, they tend not to support output from disparate pre-processing pipelines, may have limited documentation, and may not follow BIDS best practices. Here we present XCP-D, which presents a solution to these issues. XCP-D is a collaborative effort between PennLINC at the University of Pennsylvania and the DCAN lab at the University at Minnesota. XCP-D uses an open development model on GitHub and incorporates continuous integration testing; it is distributed as a Docker container or Singularity image. XCP-D generates denoised BOLD images and functional derivatives from resting-state data in either NifTI or CIFTI files, following pre-processing with fMRIPrep, HCP, and ABCD-BIDS pipelines. Even prior to its official release, XCP-D has been downloaded >3,000 times from DockerHub. Together, XCP-D facilitates robust, scalable, and reproducible post-processing of fMRI data.

Individual differences in delay discounting are associated with dorsal prefrontal cortex connectivity in children, adolescents, and adults

Delay discounting is a measure of impulsive choice relevant in adolescence as it predicts many real-life outcomes, including obesity and academic achievement. However, resting-state functional networks underlying individual differences in delay discounting during youth remain incompletely described. Here we investigate the association between multivariate patterns of functional connectivity and individual differences in impulsive choice in a large sample of children, adolescents, and adults. A total of 293 participants (9-23 years) completed a delay discounting task and underwent 3T resting-state fMRI. A connectome-wide analysis using multivariate distance-based matrix regression was used to examine whole-brain relationships between delay discounting and functional connectivity. These analyses revealed that individual differences in delay discounting were associated with patterns of connectivity emanating from the left dorsal prefrontal cortex, a default mode network hub. Greater delay discounting was associated with greater functional connectivity between the dorsal prefrontal cortex and other default mode network regions, but reduced connectivity with regions in the dorsal and ventral attention networks. These results suggest delay discounting in children, adolescents, and adults is associated with individual differences in relationships both within the default mode network and between the default mode and networks involved in attentional and cognitive control.

Microscale Neuronal Activity Collectively Drives Chaotic and Inflexible Dynamics at the Macroscale in Seizures

Neuronal activity propagates through the network during seizures, engaging brain dynamics at multiple scales. Such propagating events can be described through the avalanches framework, which can relate spatiotemporal activity at the microscale with global network properties. Interestingly, propagating avalanches in healthy networks are indicative of critical dynamics, where the network is organised to a phase transition, which optimises certain computational properties. Some have hypothesised that the pathological brain dynamics of epileptic seizures are an emergent property of microscale neuronal networks collectively driving the brain away from criticality. Demonstrating this would provide a unifying mechanism linking microscale spatiotemporal activity with emergent brain dysfunction during seizures. Here, we investigated the effect of drug-induced seizures on critical avalanche dynamics, using in vivo whole-brain 2-photon imaging of GCaMP6s larval zebrafish (males and females) at single neuron resolution. We demonstrate that single neuron activity across the whole brain exhibits a loss of critical statistics during seizures, suggesting that microscale activity collectively drives macroscale dynamics away from criticality. We also construct spiking network models at the scale of the larval zebrafish brain, to demonstrate that only densely connected networks can drive brain-wide seizure dynamics away from criticality. Importantly, such dense networks also disrupt the optimal computational capacities of critical networks, leading to chaotic dynamics, impaired network response properties and sticky states, thus helping to explain functional impairments during seizures. This study bridges the gap between microscale neuronal activity and emergent macroscale dynamics and cognitive dysfunction during seizures.Significance StatementEpileptic seizures are debilitating and impair normal brain function. It is unclear how the coordinated behaviour of neurons collectively impairs brain function during seizures. To investigate this we perform fluorescence microscopy in larval zebrafish, which allows for the recording of whole-brain activity at single-neuron resolution. Using techniques from physics, we show that neuronal activity during seizures drives the brain away from criticality, a regime that enables both high and low activity states, into an inflexible regime that drives high activity states. Importantly, this change is caused by more connections in the network, which we show disrupts the ability of the brain to respond appropriately to its environment. Therefore, we identify key neuronal network mechanisms driving seizures and concurrent cognitive dysfunction.

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.

Sex differences in the functional topography of association networks in youth

Prior work has shown that there is substantial interindividual variation in the spatial distribution of functional networks across the cerebral cortex, or functional topography. However, it remains unknown whether there are sex differences in the topography of individualized networks in youth. Here, we leveraged an advanced machine learning method (sparsity-regularized non-negative matrix factorization) to define individualized functional networks in 693 youth (ages 8 to 23 y) who underwent functional MRI as part of the Philadelphia Neurodevelopmental Cohort. Multivariate pattern analysis using support vector machines classified participant sex based on functional topography with 82.9% accuracy (P < 0.0001). Brain regions most effective in classifying participant sex belonged to association networks, including the ventral attention, default mode, and frontoparietal networks. Mass univariate analyses using generalized additive models with penalized splines provided convergent results. Furthermore, transcriptomic data from the Allen Human Brain Atlas revealed that sex differences in multivariate patterns of functional topography were spatially correlated with the expression of genes on the X chromosome. These results highlight the role of sex as a biological variable in shaping functional topography.

Mobile footprinting: linking individual distinctiveness in mobility patterns to mood, sleep, and brain functional connectivity

Mapping individual differences in behavior is fundamental to personalized neuroscience, but quantifying complex behavior in real world settings remains a challenge. While mobility patterns captured by smartphones have increasingly been linked to a range of psychiatric symptoms, existing research has not specifically examined whether individuals have person-specific mobility patterns. We collected over 3000 days of mobility data from a sample of 41 adolescents and young adults (age 17-30 years, 28 female) with affective instability. We extracted summary mobility metrics from GPS and accelerometer data and used their covariance structures to identify individuals and calculated the individual identification accuracy-i.e., their "footprint distinctiveness". We found that statistical patterns of smartphone-based mobility features represented unique "footprints" that allow individual identification (p < 0.001). Critically, mobility footprints exhibited varying levels of person-specific distinctiveness (4-99%), which was associated with age and sex. Furthermore, reduced individual footprint distinctiveness was associated with instability in affect (p < 0.05) and circadian patterns (p < 0.05) as measured by environmental momentary assessment. Finally, brain functional connectivity, especially those in the somatomotor network, was linked to individual differences in mobility patterns (p < 0.05). Together, these results suggest that real-world mobility patterns may provide individual-specific signatures relevant for studies of development, sleep, and psychopathology.

Model-based stationarity filtering of long-term memory data applied to resting-state blood-oxygen-level-dependent signal

Resting-state blood-oxygen-level-dependent (BOLD) signal acquired through functional magnetic resonance imaging is a proxy of neural activity and a key mechanism for assessing neurological conditions. Therefore, practical tools to filter out artefacts that can compromise the assessment are required. On the one hand, a variety of tailored methods to preprocess the data to deal with identified sources of noise (e.g., head motion, heart beating, and breathing, just to mention a few) are in place. But, on the other hand, there might be unknown sources of unstructured noise present in the data. Therefore, to mitigate the effects of such unstructured noises, we propose a model-based filter that explores the statistical properties of the underlying signal (i.e., long-term memory). Specifically, we consider autoregressive fractional integrative process filters. Remarkably, we provide evidence that such processes can model the signals at different regions of interest to attain stationarity. Furthermore, we use a principled analysis where a ground-truth signal with statistical properties similar to the BOLD signal under the injection of noise is retrieved using the proposed filters. Next, we considered preprocessed (i.e., the identified sources of noise removed) resting-state BOLD data of 98 subjects from the Human Connectome Project. Our results demonstrate that the proposed filters decrease the power in the higher frequencies. However, unlike the low-pass filters, the proposed filters do not remove all high-frequency information, instead they preserve process-related higher frequency information. Additionally, we considered four different metrics (power spectrum, functional connectivity using the Pearson’s correlation, coherence, and eigenbrains) to infer the impact of such filter. We provided evidence that whereas the first three keep most of the features of interest from a neuroscience perspective unchanged, the latter exhibits some variations that could be due to the sporadic activity filtered out.

Dissociable multi-scale patterns of development in personalized brain networks

The brain is organized into networks at multiple resolutions, or scales, yet studies of functional network development typically focus on a single scale. Here, we derive personalized functional networks across 29 scales in a large sample of youths (n = 693, ages 8-23 years) to identify multi-scale patterns of network re-organization related to neurocognitive development. We found that developmental shifts in inter-network coupling reflect and strengthen a functional hierarchy of cortical organization. Furthermore, we observed that scale-dependent effects were present in lower-order, unimodal networks, but not higher-order, transmodal networks. Finally, we found that network maturation had clear behavioral relevance: the development of coupling in unimodal and transmodal networks are dissociably related to the emergence of executive function. These results suggest that the development of functional brain networks align with and refine a hierarchy linked to cognition.

Daily sensation-seeking and urgency in young adults: Examining associations with alcohol use and self-defined risky behaviors

PURPOSE

To examine the dynamic expression of sensation-seeking and urgency in daily life and the implications for alcohol use and risk-taking during young adulthood.

METHODS

Daily diary surveys were administered to young adults (n = 77) aged 18-25 years every evening for 21 days to assess day's sensation-seeking, urgency, risk-taking, and alcohol use.

RESULTS

Days of higher than usual sensation-seeking are also days of higher than usual risk-taking and are more likely to be alcohol use days than days of lower than usual sensation-seeking. Day's urgency was not associated with day's alcohol use or risk-taking. We extracted 10 themes from self-reports of the day's riskiest behavior: transportation (29.9%), social (22.8%), recreation (17.4%), work (14.8%), school (13.5%), food (9.5%), sleep (9.2%), substance use (5.8%), other (5.2%), and jaywalking (1.5%), and 14.6% of self-reported risky behaviors were considered threatening to safety, health, or wellbeing.

CONCLUSIONS

Risks taken during daily life have mostly positive outcomes and a minority represent threats to safety, health, and wellbeing. Risk-taking and alcohol use in young adult's daily lives is more likely to be driven by the desire to experience novel and exciting experiences than by rash action.

Daily fluctuations in young children's persistence

Children's behavior changes from day to day, but the factors that contribute to its variability are understudied. We developed a novel repeated measures paradigm to study children's persistence by capitalizing on a task that children complete every day: toothbrushing (N = 81; 48% female; 36-47 months; 80% white, 14% Multiracial, 10% Hispanic, 2% Asian, 1% Black; 1195 observations collected between January 2019 and March 2020). Children brushed longer on days when their parents used more praise (d = .23) and less instruction (d = -.22). Sensitivity to mood, sleep, and parent stress varied across children, suggesting that identifying the factors that shape an individual child's persistence could lead to personalized interventions.

A multilayer network model of neuron-astrocyte populations in vitro reveals mGluR5 inhibition is protective following traumatic injury

Astrocytes communicate bidirectionally with neurons, enhancing synaptic plasticity and promoting the synchronization of neuronal microcircuits. Despite recent advances in understanding neuron-astrocyte signaling, little is known about astrocytic modulation of neuronal activity at the population level, particularly in disease or following injury. We used high-speed calcium imaging of mixed cortical cultures in vitro to determine how population activity changes after disruption of glutamatergic signaling and mechanical injury. We constructed a multilayer network model of neuron-astrocyte connectivity, which captured distinct topology and response behavior from single-cell-type networks. mGluR₅ inhibition decreased neuronal activity, but did not on its own disrupt functional connectivity or network topology. In contrast, injury increased the strength, clustering, and efficiency of neuronal but not astrocytic networks, an effect that was not observed in networks pretreated with mGluR₅ inhibition. Comparison of spatial and functional connectivity revealed that functional connectivity is largely independent of spatial proximity at the microscale, but mechanical injury increased the spatial-functional correlation. Finally, we found that astrocyte segments of the same cell often belong to separate functional communities based on neuronal connectivity, suggesting that astrocyte segments function as independent entities. Our findings demonstrate the utility of multilayer network models for characterizing the multiscale connectivity of two distinct but functionally dependent cell populations.

Astrocytes communicate bidirectionally with neurons, enhancing synaptic plasticity and promoting the synchronization of neuronal microcircuits. We constructed a multilayer network model of neuron-astrocyte connectivity based on calcium activity in mixed cortical cultures, and used this model to evaluate the effect of glutamatergic inhibition and mechanical injury on network topology. We found that injury increased the strength, clustering, and efficiency of neuronal but not astrocytic networks, an effect that was not observed in injured networks pretreated with a glutamate receptor antagonist. Our findings demonstrate the utility of multilayer network models for characterizing the multiscale connectivity of two distinct but functionally dependent cell populations.

Functional brain network community structure in childhood: Unfinished territories and fuzzy boundaries

Adult cortex is organized into distributed functional communities. Yet, little is known about community architecture of children’s brains. Here, we uncovered the community structure of cortex in childhood using fMRI data from 670 children aged 9–11 years (48% female, replication sample n = 544, 56% female) from the Adolescent Brain and Cognitive Development study. We first applied a data-driven community detection approach to cluster cortical regions into communities, then employed a generative model-based approach called the weighted stochastic block model to further probe community interactions. Children showed similar community structure to adults, as defined by Yeo and colleagues in 2011, in early-developing sensory and motor communities, but differences emerged in transmodal areas. Children have more cortical territory in the limbic community, which is involved in emotion processing, than adults. Regions in association cortex interact more flexibly across communities, creating uncertainty for the model-based assignment algorithm, and perhaps reflecting cortical boundaries that are not yet solidified. Uncertainty was highest for cingulo-opercular areas involved in flexible deployment of cognitive control. Activation and deactivation patterns during a working memory task showed that both the data-driven approach and a set of adult communities statistically capture functional organization in middle childhood. Collectively, our findings suggest that community boundaries are not solidified by middle childhood.

Neurodevelopment of the association cortices: Patterns, mechanisms, and implications for psychopathology

The human brain undergoes a prolonged period of cortical development that spans multiple decades. During childhood and adolescence, cortical development progresses from lower-order, primary and unimodal cortices with sensory and motor functions to higher-order, transmodal association cortices subserving executive, socioemotional, and mentalizing functions. The spatiotemporal patterning of cortical maturation thus proceeds in a hierarchical manner, conforming to an evolutionarily rooted, sensorimotor-to-association axis of cortical organization. This developmental program has been characterized by data derived from multimodal human neuroimaging and is linked to the hierarchical unfolding of plasticity-related neurobiological events. Critically, this developmental program serves to enhance feature variation between lower-order and higher-order regions, thus endowing the brain's association cortices with unique functional properties. However, accumulating evidence suggests that protracted plasticity within late-maturing association cortices, which represents a defining feature of the human developmental program, also confers risk for diverse developmental psychopathologies.

Network Controllability in Transmodal Cortex Predicts Positive Psychosis Spectrum Symptoms

BACKGROUND

The psychosis spectrum (PS) is associated with structural dysconnectivity concentrated in transmodal cortex. However, understanding of this pathophysiology has been limited by an overreliance on examining direct interregional connectivity. Using network control theory, we measured variation in both direct and indirect connectivity to a region to gain new insights into the pathophysiology of the PS.

METHODS

We used psychosis symptom data and structural connectivity in 1068 individuals from the Philadelphia Neurodevelopmental Cohort. Applying a network control theory metric called average controllability, we estimated each brain region's capacity to leverage its direct and indirect structural connections to control linear brain dynamics. Using nonlinear regression, we determined the accuracy with which average controllability could predict PS symptoms in out-of-sample testing. We also examined the predictive performance of regional strength, which indexes only direct connections to a region, as well as several graph-theoretic measures of centrality that index indirect connectivity. Finally, we assessed how the prediction performance for PS symptoms varied over the functional hierarchy spanning unimodal to transmodal cortex.

RESULTS

Average controllability outperformed all other connectivity features at predicting positive PS symptoms and was the only feature to yield above-chance predictive performance. Improved prediction for average controllability was concentrated in transmodal cortex, whereas prediction performance for strength was uniform across the cortex, suggesting that indexing indirect connections through average controllability is crucial in association cortex.

CONCLUSIONS

Examining interindividual variation in direct and indirect structural connections to transmodal cortex is crucial for accurate prediction of positive PS symptoms.

Synchronization of coupled Kuramoto oscillators under resource constraints

A fundamental understanding of synchronized behavior in multiagent systems can be acquired by studying analytically tractable Kuramoto models. However, such models typically diverge from many real systems whose dynamics evolve under nonnegligible resource constraints. Here we construct a system of coupled Kuramoto oscillators that consume or produce resources as a function of their oscillation frequency. At high coupling, we observe strongly synchronized dynamics, whereas at low coupling, we observe independent oscillator dynamics as expected from standard Kuramoto models. For intermediate coupling, which typically induces a partially synchronized state, we empirically observe that (and theoretically explain why) the system can exist in either: (i) a state in which the order parameter oscillates in time, or (ii) a state in which multiple synchronization states are simultaneously stable. Whether (i) or (ii) occurs depends upon whether the oscillators consume or produce resources, respectively. Relevant for systems as varied as coupled neurons and social groups, our paper lays important groundwork for future efforts to develop quantitative predictions of synchronized dynamics for systems embedded in environments marked by sparse resources.

Quantifying the compressibility of complex networks

Many complex networks depend upon biological entities for their preservation. Such entities, from human cognition to evolution, must first encode and then replicate those networks under marked resource constraints. Networks that survive are those that are amenable to constrained encoding-or, in other words, are compressible. But how compressible is a network? And what features make one network more compressible than another? Here, we answer these questions by modeling networks as information sources before compressing them using rate-distortion theory. Each network yields a unique rate-distortion curve, which specifies the minimal amount of information that remains at a given scale of description. A natural definition then emerges for the compressibility of a network: the amount of information that can be removed via compression, averaged across all scales. Analyzing an array of real and model networks, we demonstrate that compressibility increases with two common network properties: transitivity (or clustering) and degree heterogeneity. These results indicate that hierarchical organization-which is characterized by modular structure and heterogeneous degrees-facilitates compression in complex networks. Generally, our framework sheds light on the interplay between a network's structure and its capacity to be compressed, enabling investigations into the role of compression in shaping real-world networks.

Evaluating the sensitivity of functional connectivity measures to motion artifact in resting-state fMRI data

Functional connectivity (FC) networks are typically inferred from resting-state fMRI data using the Pearson correlation between BOLD time series from pairs of brain regions. However, alternative methods of estimating functional connectivity have not been systematically tested for their sensitivity or robustness to head motion artifact. Here, we evaluate the sensitivity of eight different functional connectivity measures to motion artifact using resting-state data from the Human Connectome Project. We report that FC estimated using full correlation has a relatively high residual distance-dependent relationship with motion compared to partial correlation, coherence, and information theory-based measures, even after implementing rigorous methods for motion artifact mitigation. This disadvantage of full correlation, however, may be offset by higher test-retest reliability, fingerprinting accuracy, and system identifiability. FC estimated by partial correlation offers the best of both worlds, with low sensitivity to motion artifact and intermediate system identifiability, with the caveat of low test-retest reliability and fingerprinting accuracy. We highlight spatial differences in the sub-networks affected by motion with different FC metrics. Further, we report that intra-network edges in the default mode and retrosplenial temporal sub-networks are highly correlated with motion in all FC methods. Our findings indicate that the method of estimating functional connectivity is an important consideration in resting-state fMRI studies and must be chosen carefully based on the parameters of the study.

Electrocorticography and stereo EEG provide distinct measures of brain connectivity: implications for network models

Brain network models derived from graph theory have the potential to guide functional neurosurgery, and to improve rates of post-operative seizure freedom for patients with epilepsy. A barrier to applying these models clinically is that intracranial EEG electrode implantation strategies vary by centre, region and country, from cortical grid & strip electrodes (Electrocorticography), to purely stereotactic depth electrodes (Stereo EEG), to a mixture of both. To determine whether models derived from one type of study are broadly applicable to others, we investigate the differences in brain networks mapped by electrocorticography and stereo EEG in a cohort of patients who underwent surgery for temporal lobe epilepsy and achieved a favourable outcome. We show that networks derived from electrocorticography and stereo EEG define distinct relationships between resected and spared tissue, which may be driven by sampling bias of temporal depth electrodes in patients with predominantly cortical grids. We propose a method of correcting for the effect of internodal distance that is specific to electrode type and explore how additional methods for spatially correcting for sampling bias affect network models. Ultimately, we find that smaller surgical targets tend to have lower connectivity with respect to the surrounding network, challenging notions that abnormal connectivity in the epileptogenic zone is typically high. Our findings suggest that effectively applying computational models to localize epileptic networks requires accounting for the effects of spatial sampling, particularly when analysing both electrocorticography and stereo EEG recordings in the same cohort, and that future network studies of epilepsy surgery should also account for differences in focality between resection and ablation. We propose that these findings are broadly relevant to intracranial EEG network modelling in epilepsy and an important step in translating them clinically into patient care.

Gender bias in academia: A lifetime problem that needs solutions

Despite increased awareness of the lack of gender equity in academia and a growing number of initiatives to address issues of diversity, change is slow, and inequalities remain. A major source of inequity is gender bias, which has a substantial negative impact on the careers, work-life balance, and mental health of underrepresented groups in science. Here, we argue that gender bias is not a single problem but manifests as a collection of distinct issues that impact researchers' lives. We disentangle these facets and propose concrete solutions that can be adopted by individuals, academic institutions, and society.

The modulation of brain network integration and arousal during exploration

There is growing interest in how neuromodulators shape brain networks. Recent neuroimaging studies provide evidence that brainstem arousal systems, such as the locus coeruleus-norepinephrine system (LC-NE), influence functional connectivity and brain network topology, suggesting they have a role in flexibly reconfiguring brain networks in order to adapt behavior and cognition to environmental demands. To date, however, the relationship between brainstem arousal systems and functional connectivity has not been assessed within the context of a task with an established relationship between arousal and behavior, with most prior studies relying on incidental variations in arousal or pharmacological manipulation and static brain networks constructed over long periods of time. These factors have likely contributed to a heterogeneity of effects across studies. To address these issues, we took advantage of the association between LC-NE-linked arousal and exploration to probe the relationships between exploratory choice, arousal—as measured indirectly via pupil diameter—and brain network dynamics. Exploration in a bandit task was associated with a shift toward fewer, more weakly connected modules that were more segregated in terms of connectivity and topology but more integrated with respect to the diversity of cognitive systems represented in each module. Functional connectivity strength decreased, and changes in connectivity were correlated with changes in pupil diameter, in line with the hypothesis that brainstem arousal systems influence the dynamic reorganization of brain networks. More broadly, we argue that carefully aligning dynamic network analyses with task designs can increase the temporal resolution at which behaviorally- and cognitively-relevant modulations can be identified, and offer these results as a proof of concept of this approach.

Predicted disconnectome associated with progressive periventricular white matter ischemia

We used a virtual lesion DTI fiber tracking approach with healthy subject DTI data and simulated periventricular white matter (PVWM) lesion masks to predict the sequence of connectivity changes associated with progressive PVWM ischemia. We found that the optic radiations, inferior fronto-occipital fasciculus, inferior longitudinal fasciculus, corpus callosum, temporopontine tract and fornix were affected in early simulated ischemic injury, and that the connectivity of subcortical, cerebellar, and visual regions were significantly disrupted with increasing simulated lesion severity. The results of this study provide insights into the spatial-temporal changes of the brain structural connectome under progressive PVWM ischemia. The virtual lesion approach provides a meaningful proxy to the spatial-temporal changes of the brain's structural connectome and can be used to further characterize the cognitive sequelae of progressive PVWM ischemia in both normal aging and dementia.

QSIPrep: an integrative platform for preprocessing and reconstructing diffusion MRI data

Diffusion-weighted magnetic resonance imaging (dMRI) is the primary method for noninvasively studying the organization of white matter in the human brain. Here we introduce QSIPrep, an integrative software platform for the processing of diffusion images that is compatible with nearly all dMRI sampling schemes. Drawing on a diverse set of software suites to capitalize on their complementary strengths, QSIPrep facilitates the implementation of best practices for processing of diffusion images.

Response inhibition in adolescents is moderated by brain connectivity and social network structure

The social environment an individual is embedded in influences their ability and motivation to engage self-control processes, but little is known about the neural mechanisms underlying this effect. Many individuals successfully regulate their behavior even when they do not show strong activation in canonical self-control brain regions. Thus, individuals may rely on other resources to compensate, including daily experiences navigating and managing complex social relationships that likely bolster self-control processes. Here, we employed a network neuroscience approach to investigate the role of social context and social brain systems in facilitating self-control in adolescents. We measured brain activation using functional magnetic resonance imaging (fMRI) as 62 adolescents completed a Go/No-Go response inhibition task. We found that self-referential brain systems compensate for weaker activation in executive function brain systems, especially for adolescents with more friends and more communities in their social networks. Collectively, our results indicate a critical role for self-referential brain systems during the developmental trajectory of self-control throughout adolescence.

DotMotif: an open-source tool for connectome subgraph isomorphism search and graph queries

Recent advances in neuroscience have enabled the exploration of brain structure at the level of individual synaptic connections. These connectomics datasets continue to grow in size and complexity; methods to search for and identify interesting graph patterns offer a promising approach to quickly reduce data dimensionality and enable discovery. These graphs are often too large to be analyzed manually, presenting significant barriers to searching for structure and testing hypotheses. We combine graph database and analysis libraries with an easy-to-use neuroscience grammar suitable for rapidly constructing queries and searching for subgraphs and patterns of interest. Our approach abstracts many of the computer science and graph theory challenges associated with nanoscale brain network analysis and allows scientists to quickly conduct research at scale. We demonstrate the utility of these tools by searching for motifs on simulated data and real public connectomics datasets, and we share simple and complex structures relevant to the neuroscience community. We contextualize our findings and provide case studies and software to motivate future neuroscience exploration.

Phase/Amplitude Synchronization of Brain Signals During Motor Imagery BCI Tasks

In the last decade, functional connectivity (FC) has been increasingly adopted based on its ability to capture statistical dependencies between multivariate brain signals. However, the role of FC in the context of brain-computer interface applications is still poorly understood. To address this gap in knowledge, we considered a group of 20 healthy subjects during an EEG-based hand motor imagery (MI) task. We studied two well-established FC estimators, i.e. spectral- and imaginary-coherence, and we investigated how they were modulated by the MI task. We characterized the resulting FC networks by extracting the strength of connectivity of each EEG sensor and we compared the discriminant power with respect to standard power spectrum features. At the group level, results showed that while spectral-coherence based network features were increasing in the sensorimotor areas, those based on imaginary-coherence were significantly decreasing. We demonstrated that this opposite, but complementary, behavior was respectively determined by the increase in amplitude and phase synchronization between the brain signals. At the individual level, we eventually assessed the potential of these network connectivity features in a simple off-line classification scenario. Taken together, our results provide fresh insights into the oscillatory mechanisms subserving brain network changes during MI and offer new perspectives to improve BCI performance.

Brain network dynamics during working memory are modulated by dopamine and diminished in schizophrenia

Dynamical brain state transitions are critical for flexible working memory but the network mechanisms are incompletely understood. Here, we show that working memory performance entails brain-wide switching between activity states using a combination of functional magnetic resonance imaging in healthy controls and individuals with schizophrenia, pharmacological fMRI, genetic analyses and network control theory. The stability of states relates to dopamine D1 receptor gene expression while state transitions are influenced by D2 receptor expression and pharmacological modulation. Individuals with schizophrenia show altered network control properties, including a more diverse energy landscape and decreased stability of working memory representations. Our results demonstrate the relevance of dopamine signaling for the steering of whole-brain network dynamics during working memory and link these processes to schizophrenia pathophysiology.

Working memory requires the brain to switch between cognitive states and activity patterns. Here, the authors show that the steering of these neural network dynamics is influenced by dopamine D1- and D2-receptor function and altered in schizophrenia.

Environmental influences on the pace of brain development

Childhood socio-economic status (SES), a measure of the availability of material and social resources, is one of the strongest predictors of lifelong well-being. Here we review evidence that experiences associated with childhood SES affect not only the outcome but also the pace of brain development. We argue that higher childhood SES is associated with protracted structural brain development and a prolonged trajectory of functional network segregation, ultimately leading to more efficient cortical networks in adulthood. We hypothesize that greater exposure to chronic stress accelerates brain maturation, whereas greater access to novel positive experiences decelerates maturation. We discuss the impact of variation in the pace of brain development on plasticity and learning. We provide a generative theoretical framework to catalyse future basic science and translational research on environmental influences on brain development.

Computational modeling of tau pathology spread reveals patterns of regional vulnerability and the impact of a genetic risk factor

Neuropathological staging studies have suggested that tau pathology spreads through the brain in Alzheimer's disease (AD) and other tauopathies, but it is unclear how neuroanatomical connections, spatial proximity, and regional vulnerability contribute. In this study, we seed tau pathology in the brains of nontransgenic mice with AD tau and quantify pathology development over 9 months in 134 brain regions. Network modeling of pathology progression shows that diffusion through the connectome is the best predictor of tau pathology patterns. Further, deviations from pure neuroanatomical spread are used to estimate regional vulnerability to tau pathology and identify related gene expression patterns. Last, we show that pathology spread is altered in mice harboring a mutation in leucine-rich repeat kinase 2. While tau pathology spread is still constrained by anatomical connectivity in these mice, it spreads preferentially in a retrograde direction. This study provides a framework for understanding neuropathological progression in tauopathies.

Structural control energy of resting-state functional brain states reveals less cost-effective brain dynamics in psychosis vulnerability

How the brain's white-matter anatomy constrains brain activity is an open question that might give insights into the mechanisms that underlie mental disorders such as schizophrenia. Chromosome 22q11.2 deletion syndrome (22q11DS) is a neurodevelopmental disorder with an extremely high risk for psychosis providing a test case to study developmental aspects of schizophrenia. In this study, we used principles from network control theory to probe the implications of aberrant structural connectivity for the brain's functional dynamics in 22q11DS. We retrieved brain states from resting-state functional magnetic resonance images of 78 patients with 22q11DS and 85 healthy controls. Then, we compared them in terms of persistence control energy; that is, the control energy that would be required to persist in each of these states based on individual structural connectivity and a dynamic model. Persistence control energy was altered in a broad pattern of brain states including both energetically more demanding and less demanding brain states in 22q11DS. Further, we found a negative relationship between persistence control energy and resting-state activation time, which suggests that the brain reduces energy by spending less time in energetically demanding brain states. In patients with 22q11DS, this behavior was less pronounced, suggesting a deficiency in the ability to reduce energy through brain activation. In summary, our results provide initial insights into the functional implications of altered structural connectivity in 22q11DS, which might improve our understanding of the mechanisms underlying the disease.

Crystalline shielding mitigates structural rearrangement and localizes memory in jammed systems under oscillatory shear

The nature of yield in amorphous materials under stress has yet to be fully elucidated. In particular, understanding how microscopic rearrangement gives rise to macroscopic structural and rheological signatures in disordered systems is vital for the prediction and characterization of yield and the study of how memory is stored in disordered materials. Here, we investigate the evolution of local structural homogeneity on an individual particle level in amorphous jammed two-dimensional (athermal) systems under oscillatory shear and relate this evolution to rearrangement, memory, and macroscale rheological measurements. We define the structural metric crystalline shielding, and show that it is predictive of rearrangement propensity and structural volatility of individual particles under shear. We use this metric to identify localized regions of the system in which the material's memory of its preparation is preserved. Our results contribute to a growing understanding of how local structure relates to dynamic response and memory in disordered systems.

The Citation Diversity Statement: A Practice of Transparency, A Way of Life

Appending a Citation Diversity Statement to a paper is a simple and effective way to increase awareness about citation bias and help mitigate it. Here, we describe why reducing citation bias is important and how to include a Citation Diversity Statement in your next publication.

Transdiagnostic dimensions of psychopathology explain individuals' unique deviations from normative neurodevelopment in brain structure

Psychopathology is rooted in neurodevelopment. However, clinical and biological heterogeneity, together with a focus on case-control approaches, have made it difficult to link dimensions of psychopathology to abnormalities of neurodevelopment. Here, using the Philadelphia Neurodevelopmental Cohort, we built normative models of cortical volume and tested whether deviations from these models better predicted psychiatric symptoms compared to raw cortical volume. Specifically, drawing on the p-factor hypothesis, we distilled 117 clinical symptom measures into six orthogonal psychopathology dimensions: overall psychopathology, anxious-misery, externalizing disorders, fear, positive psychosis symptoms, and negative psychosis symptoms. We found that multivariate patterns of deviations yielded improved out-of-sample prediction of psychopathology dimensions compared to multivariate patterns of raw cortical volume. We also found that correlations between overall psychopathology and deviations in ventromedial prefrontal, inferior temporal, and dorsal anterior cingulate cortices were stronger than those observed for specific dimensions of psychopathology (e.g., anxious-misery). Notably, these same regions are consistently implicated in a range of putatively distinct disorders. Finally, we performed conventional case-control comparisons of deviations in a group of individuals with depression and a group with attention-deficit hyperactivity disorder (ADHD). We observed spatially overlapping effects between these groups that diminished when controlling for overall psychopathology. Together, our results suggest that modeling cortical brain features as deviations from normative neurodevelopment improves prediction of psychiatric symptoms in out-of-sample testing, and that p-factor models of psychopathology may assist in separating biomarkers that are disorder-general from those that are disorder-specific.

Associations between Neighborhood SES and Functional Brain Network Development

Higher socioeconomic status (SES) in childhood is associated with stronger cognitive abilities, higher academic achievement, and lower incidence of mental illness later in development. While prior work has mapped the associations between neighborhood SES and brain structure, little is known about the relationship between SES and intrinsic neural dynamics. Here, we capitalize upon a large cross-sectional community-based sample (Philadelphia Neurodevelopmental Cohort, ages 8-22 years, n = 1012) to examine associations between age, SES, and functional brain network topology. We characterize this topology using a local measure of network segregation known as the clustering coefficient and find that it accounts for a greater degree of SES-associated variance than mesoscale segregation captured by modularity. High-SES youth displayed stronger positive associations between age and clustering than low-SES youth, and this effect was most pronounced for regions in the limbic, somatomotor, and ventral attention systems. The moderating effect of SES on positive associations between age and clustering was strongest for connections of intermediate length and was consistent with a stronger negative relationship between age and local connectivity in these regions in low-SES youth. Our findings suggest that, in late childhood and adolescence, neighborhood SES is associated with variation in the development of functional network structure in the human brain.

Atypical neural topographies underpin dysfunctional pattern separation in temporal lobe epilepsy

Episodic memory is the ability to accurately remember events from our past. The process of pattern separation is hypothesized to underpin this ability and is defined as the ability to orthogonalize memory traces, to maximize the features that make them unique. Contemporary cognitive neuroscience suggests that pattern separation entails complex interactions between the hippocampus and the neocortex, where specific hippocampal subregions shape neural reinstatement in the neocortex. To test this hypothesis, the current work studied both healthy controls and patients with temporal lobe epilepsy (TLE) who present with hippocampal structural anomalies. In all participants, we measured neural activity using functional magnetic resonance imaging (fMRI) while they retrieved memorized items compared to lure items which share features with the target. Behaviorally, TLE patients were less able to exclude lures than controls, and showed a reduction in pattern separation. To assess the hypothesized relationship between neural patterns in the hippocampus and the neocortex, we identified topographic gradients of intrinsic connectivity along neocortical and hippocampal subfield surfaces and identified the topographic profile of the neural activity accompanying pattern separation. In healthy controls, pattern separation followed a graded pattern of neural activity, both along the hippocampal long axis (and peaked in anterior segments that are more heavily engaged in transmodal processing) and along the neocortical hierarchy running from unimodal to transmodal regions (peaking in transmodal default mode regions). In TLE patients, however, this concordance between task-based functional activations and topographic gradients was markedly reduced. Furthermore, person specific measures of concordance between task-related activity and connectivity gradients in patients and controls related to inter-individual differences in behavioral measures of pattern separation and episodic memory, highlighting the functional relevance of the observed topographic motifs. Our work is consistent with an emerging understanding that successful discrimination between memories with similar features entails a shift in the locus of neural activity away from sensory systems, a pattern that is mirrored along the hippocampal long axis and with respect to neocortical hierarchies. More broadly, our study establishes topographic profiling using intrinsic connectivity gradients captures the functional underpinnings of episodic memory processes in manner that is sensitive to their reorganization in pathology.

BCI learning induces core-periphery reorganization in M/EEG multiplex brain networks

Brain-computer interfaces (BCIs) constitute a promising tool for communication and control. However, mastering non-invasive closed-loop systems remains a learned skill that is difficult to develop for a non-negligible proportion of users. The involved learning process induces neural changes associated with a brain network reorganization that remains poorly understood. To address this inter-subject variability, we adopted a multilayer approach to integrate brain network properties from electroencephalographic (EEG) and magnetoencephalographic (MEG) data resulting from a four-session BCI training program followed by a group of healthy subjects. Our method gives access to the contribution of each layer to multilayer network that tends to be equal with time. We show that regardless the chosen modality, a progressive increase in the integration of somatosensory areas in the α band was paralleled by a decrease of the integration of visual processing and working memory areas in the β band. Notably, only brain network properties in multilayer network correlated with future BCI scores in the α2 band: positively in somatosensory and decision-making related areas and negatively in associative areas. Our findings cast new light on neural processes underlying BCI training. Integrating multimodal brain network properties provides new information that correlates with behavioral performance and could be considered as a potential marker of BCI learning.

Data-driven control of complex networks

Our ability to manipulate the behavior of complex networks depends on the design of efficient control algorithms and, critically, on the availability of an accurate and tractable model of the network dynamics. While the design of control algorithms for network systems has seen notable advances in the past few years, knowledge of the network dynamics is a ubiquitous assumption that is difficult to satisfy in practice. In this paper we overcome this limitation, and develop a data-driven framework to control a complex network optimally and without any knowledge of the network dynamics. Our optimal controls are constructed using a finite set of data, where the unknown network is stimulated with arbitrary and possibly random inputs. Although our controls are provably correct for networks with linear dynamics, we also characterize their performance against noisy data and in the presence of nonlinear dynamics, as they arise in power grid and brain networks.

Neurocognitive and functional heterogeneity in depressed youth

Depression is a common psychiatric illness that often begins in youth, and is sometimes associated with cognitive deficits. However, there is significant variability in cognitive dysfunction, likely reflecting biological heterogeneity. We sought to identify neurocognitive subtypes and their neurofunctional signatures in a large cross-sectional sample of depressed youth. Participants were drawn from the Philadelphia Neurodevelopmental Cohort, including 712 youth with a lifetime history of a major depressive episode and 712 typically developing (TD) youth matched on age and sex. A subset (MDD n = 368, TD n = 200) also completed neuroimaging. Cognition was assessed with the Penn Computerized Neurocognitive Battery. A recently developed semi-supervised machine learning algorithm was used to delineate neurocognitive subtypes. Subtypes were evaluated for differences in both clinical psychopathology and brain activation during an n-back working memory fMRI task. We identified three neurocognitive subtypes in the depressed group. Subtype 1 was high-performing (high accuracy, moderate speed), Subtype 2 was cognitively impaired (low accuracy, slow speed), and Subtype 3 was impulsive (low accuracy, fast speed). While subtypes did not differ in clinical psychopathology, they diverged in their activation profiles in regions critical for executive function, which mirrored differences in cognition. Taken together, these data suggest disparate mechanisms of cognitive vulnerability and resilience in depressed youth, which may inform the identification of biomarkers for prognosis and treatment response.

Pairwise maximum entropy model explains the role of white matter structure in shaping emergent co-activation states

A major challenge in neuroscience is determining a quantitative relationship between the brain's white matter structural connectivity and emergent activity. We seek to uncover the intrinsic relationship among brain regions fundamental to their functional activity by constructing a pairwise maximum entropy model (MEM) of the inter-ictal activation patterns of five patients with medically refractory epilepsy over an average of ~14 hours of band-passed intracranial EEG (iEEG) recordings per patient. We find that the pairwise MEM accurately predicts iEEG electrodes' activation patterns' probability and their pairwise correlations. We demonstrate that the estimated pairwise MEM's interaction weights predict structural connectivity and its strength over several frequencies significantly beyond what is expected based solely on sampled regions' distance in most patients. Together, the pairwise MEM offers a framework for explaining iEEG functional connectivity and provides insight into how the brain's structural connectome gives rise to large-scale activation patterns by promoting co-activation between connected structures.

Time-evolving controllability of effective connectivity networks during seizure progression

Over one third of the estimated 3 million people with epilepsy in the United States are medication resistant. Responsive neurostimulation from chronically implanted electrodes provides a promising treatment alternative to resective surgery. However, determining optimal personalized stimulation parameters, including when and where to intervene to guarantee a positive patient outcome, is a major open challenge. Network neuroscience and control theory offer useful tools that may guide improvements in parameter selection for control of anomalous neural activity. Here we use a method to characterize dynamic controllability across consecutive effective connectivity (EC) networks based on regularized partial correlations between implanted electrodes during the onset, propagation, and termination regimes of 34 seizures. We estimate regularized partial correlation adjacency matrices from 1-s time windows of intracranial electrocorticography recordings using the Graphical Least Absolute Shrinkage and Selection Operator (GLASSO). Average and modal controllability metrics calculated from each resulting EC network track the time-varying controllability of the brain on an evolving landscape of conditionally dependent network interactions. We show that average controllability increases throughout a seizure and is negatively correlated with modal controllability throughout. Our results support the hypothesis that the energy required to drive the brain to a seizure-free state from an ictal state is smallest during seizure onset, yet we find that applying control energy at electrodes in the seizure onset zone may not always be energetically favorable. Our work suggests that a low-complexity model of time-evolving controllability may offer insights for developing and improving control strategies targeting seizure suppression.

Multimodal in vivo recording using transparent graphene microelectrodes illuminates spatiotemporal seizure dynamics at the microscale

Neurological disorders such as epilepsy arise from disrupted brain networks. Our capacity to treat these disorders is limited by our inability to map these networks at sufficient temporal and spatial scales to target interventions. Current best techniques either sample broad areas at low temporal resolution (e.g. calcium imaging) or record from discrete regions at high temporal resolution (e.g. electrophysiology). This limitation hampers our ability to understand and intervene in aberrations of network dynamics. Here we present a technique to map the onset and spatiotemporal spread of acute epileptic seizures in vivo by simultaneously recording high bandwidth microelectrocorticography and calcium fluorescence using transparent graphene microelectrode arrays. We integrate dynamic data features from both modalities using non-negative matrix factorization to identify sequential spatiotemporal patterns of seizure onset and evolution, revealing how the temporal progression of ictal electrophysiology is linked to the spatial evolution of the recruited seizure core. This integrated analysis of multimodal data reveals otherwise hidden state transitions in the spatial and temporal progression of acute seizures. The techniques demonstrated here may enable future targeted therapeutic interventions and novel spatially embedded models of local circuit dynamics during seizure onset and evolution.

Modeling brain, symptom, and behavior in the winds of change

Neuropsychopharmacology addresses pressing questions in the study of three intertwined complex systems: the brain, human behavior, and symptoms of illness. The field seeks to understand the perturbations that impinge upon those systems, either driving greater health or illness. In the pursuit of this aim, investigators often perform analyses that make certain assumptions about the nature of the systems that are being perturbed. Those assumptions can be encoded in powerful computational models that serve to bridge the wide gulf between a descriptive analysis and a formal theory of a system's response. Here we review a set of three such models along a continuum of complexity, moving from a local treatment to a network treatment: one commonly applied form of the general linear model, impulse response models, and network control models. For each, we describe the model's basic form, review its use in the field, and provide a frank assessment of its relative strengths and weaknesses. The discussion naturally motivates future efforts to interlink data analysis, computational modeling, and formal theory. Our goal is to inspire practitioners to consider the assumptions implicit in their analytical approach, align those assumptions to the complexity of the systems under study, and take advantage of exciting recent advances in modeling the relations between perturbations and system function.

Towards precise resting-state fMRI biomarkers in psychiatry: synthesizing developments in transdiagnostic research, dimensional models of psychopathology, and normative neurodevelopment

Searching for biomarkers has been a chief pursuit of the field of psychiatry. Toward this end, studies have catalogued candidate resting-state biomarkers in nearly all forms of mental disorder. However, it is becoming increasingly clear that these biomarkers lack specificity, limiting their capacity to yield clinical impact. We discuss three avenues of research that are overcoming this limitation: (i) the adoption of transdiagnostic research designs, which involve studying and explicitly comparing multiple disorders from distinct diagnostic axes of psychiatry; (ii) dimensional models of psychopathology that map the full spectrum of symptomatology and that cut across traditional disorder boundaries; and (iii) modeling individuals' unique functional connectomes throughout development. We provide a framework for tying these subfields together that draws on tools from machine learning and network science.

Generative network models of altered structural brain connectivity in schizophrenia

Alterations in the structural connectome of schizophrenia patients have been widely characterized, but the mechanisms remain largely unknown. Generative network models have recently been introduced as a tool to test the biological underpinnings of altered brain network formation. We evaluated different generative network models in healthy controls (n=152), schizophrenia patients (n=66), and their unaffected first-degree relatives (n=32), and we identified spatial and topological factors contributing to network formation. We further investigated how these factors relate to cognition and to polygenic risk for schizophrenia. Our data show that among the four tested classes of generative network models, structural brain networks were optimally accounted for by a two-factor model combining spatial constraints and topological neighborhood structure. The same wiring model explained brain network formation across study groups. However, relatives and schizophrenia patients exhibited significantly lower spatial constraints and lower topological facilitation compared to healthy controls. Further exploratory analyses point to potential associations of the model parameter reflecting spatial constraints with the polygenic risk for schizophrenia and cognitive performance. Our results identify spatial constraints and local topological structure as two interrelated mechanisms contributing to regular brain network formation as well as altered connectomes in schizophrenia and healthy individuals at familial risk for schizophrenia. On an exploratory level, our data further point to the potential relevance of spatial constraints for the genetic risk for schizophrenia and general cognitive functioning, thereby encouraging future studies in following up on these observations to gain further insights into the biological basis and behavioral relevance of model parameters.

Network structure of cascading neural systems predicts stimulus propagation and recovery

OBJECTIVE

Many neural systems display spontaneous, spatiotemporal patterns of neural activity that are crucial for information processing. While these cascading patterns presumably arise from the underlying network of synaptic connections between neurons, the precise contribution of the network's local and global connectivity to these patterns and information processing remains largely unknown.

APPROACH

Here, we demonstrate how network structure supports information processing through network dynamics in empirical and simulated spiking neurons using mathematical tools from linear systems theory, network control theory, and information theory.

MAIN RESULTS

In particular, we show that activity, and the information that it contains, travels through cycles in real and simulated networks.

SIGNIFICANCE

Broadly, our results demonstrate how cascading neural networks could contribute to cognitive faculties that require lasting activation of neuronal patterns, such as working memory or attention.