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.

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.

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.

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.

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.

Correction: Stability of spontaneous, correlated activity in mouse auditory cortex.

[This corrects the article DOI: 10.1371/journal.pcbi.1007360.].

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.

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. mGluR5 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 mGluR5 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.

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.

Language Recovery after Brain Injury: A Structural Network Control Theory Study.

Aphasia recovery after stroke depends on the condition of the remaining, extralesional brain network. Network control theory (NCT) provides a unique, quantitative approach to assess the interaction between brain networks. In this longitudinal, large-scale, whole-brain connectome study, we evaluated whether controllability measures of language-related regions are associated with treated aphasia recovery. Using probabilistic tractography and controlling for the effects of structural lesions, we reconstructed whole-brain diffusion tensor imaging (DTI) connectomes from 68 individuals (20 female, 48 male) with chronic poststroke aphasia who completed a three-week language therapy. Applying principles of NCT, we computed regional (1) average and (2) modal controllability, which decode the ability of a region to (1) spread control input through the brain network and (2) to facilitate brain state transitions. We tested the relationship between pretreatment controllability measures of 20 language-related left hemisphere regions and improvements in naming six months after language therapy using multiple linear regressions and a parsimonious elastic net regression model with cross-validation. Regional controllability of the inferior frontal gyrus (IFG) pars opercularis, pars orbitalis, and the anterior insula were associated with treatment outcomes independently of baseline aphasia severity, lesion volume, age, education, and network size. Modal controllability of the IFG pars opercularis was the strongest predictor of treated aphasia recovery with cross-validation and outperformed traditional graph theory, lesion load, and demographic measures. Regional NCT measures can reflect the status of the residual language network and its interaction with the remaining brain network, being able to predict language recovery after aphasia treatment.SIGNIFICANCE STATEMENT Predicting and understanding language recovery after brain injury remains a challenging, albeit a fundamental aspect of human neurology and neuroscience. In this study, we applied network control theory (NCT) to fully harness the concept of brain networks as dynamic systems and to evaluate their interaction. We studied 68 stroke survivors with aphasia who underwent imaging and longitudinal behavioral assessments coupled with language therapy. We found that the controllability of the inferior frontal regional network significantly predicted recovery in language production six months after treatment. Importantly, controllability outperformed traditional demographic, lesion, and graph-theoretical measures. Our findings shed light on the neurobiological basis of human language and can be translated into personalized rehabilitation approaches.

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.

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.

Broken detailed balance and entropy production in the human brain.

Living systems break detailed balance at small scales, consuming energy and producing entropy in the environment to perform molecular and cellular functions. However, it remains unclear how broken detailed balance manifests at macroscopic scales and how such dynamics support higher-order biological functions. Here we present a framework to quantify broken detailed balance by measuring entropy production in macroscopic systems. We apply our method to the human brain, an organ whose immense metabolic consumption drives a diverse range of cognitive functions. Using whole-brain imaging data, we demonstrate that the brain nearly obeys detailed balance when at rest, but strongly breaks detailed balance when performing physically and cognitively demanding tasks. Using a dynamic Ising model, we show that these large-scale violations of detailed balance can emerge from fine-scale asymmetries in the interactions between elements, a known feature of neural systems. Together, these results suggest that violations of detailed balance are vital for cognition and provide a general tool for quantifying entropy production in macroscopic systems.

Multimodal electrophysiological analyses reveal that reduced synaptic excitatory neurotransmission underlies seizures in a model of NMDAR antibody-mediated encephalitis.

Seizures are a prominent feature in N-Methyl-D-Aspartate receptor antibody (NMDAR antibody) encephalitis, a distinct neuro-immunological disorder in which specific human autoantibodies bind and crosslink the surface of NMDAR proteins thereby causing internalization and a state of NMDAR hypofunction. To further understand ictogenesis in this disorder, and to test a potential treatment compound, we developed an NMDAR antibody mediated rat seizure model that displays spontaneous epileptiform activity in vivo and in vitro. Using a combination of electrophysiological and dynamic causal modelling techniques we show that, contrary to expectation, reduction of synaptic excitatory, but not inhibitory, neurotransmission underlies the ictal events through alterations in the dynamical behaviour of microcircuits in brain tissue. Moreover, in vitro application of a neurosteroid, pregnenolone sulphate, that upregulates NMDARs, reduced established ictal activity. This proof-of-concept study highlights the complexity of circuit disturbances that may lead to seizures and the potential use of receptor-specific treatments in antibody-mediated seizures and epilepsy.

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.

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.

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.

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.