Geometric constraints on human brain function.

The anatomy of the brain necessarily constrains its function, but precisely how remains unclear. The classical and dominant paradigm in neuroscience is that neuronal dynamics are driven by interactions between discrete, functionally specialized cell populations connected by a complex array of axonal fibres1-3. However, predictions from neural field theory, an established mathematical framework for modelling large-scale brain activity4-6, suggest that the geometry of the brain may represent a more fundamental constraint on dynamics than complex interregional connectivity7,8. Here, we confirm these theoretical predictions by analysing human magnetic resonance imaging data acquired under spontaneous and diverse task-evoked conditions. Specifically, we show that cortical and subcortical activity can be parsimoniously understood as resulting from excitations of fundamental, resonant modes of the brain's geometry (that is, its shape) rather than from modes of complex interregional connectivity, as classically assumed. We then use these geometric modes to show that task-evoked activations across over 10,000 brain maps are not confined to focal areas, as widely believed, but instead excite brain-wide modes with wavelengths spanning over 60 mm. Finally, we confirm predictions that the close link between geometry and function is explained by a dominant role for wave-like activity, showing that wave dynamics can reproduce numerous canonical spatiotemporal properties of spontaneous and evoked recordings. Our findings challenge prevailing views and identify a previously underappreciated role of geometry in shaping function, as predicted by a unifying and physically principled model of brain-wide dynamics.

Individuals with problem gambling and obsessive-compulsive disorder learn through distinct reinforcement mechanisms.

Obsessive-compulsive disorder (OCD) and pathological gambling (PG) are accompanied by deficits in behavioural flexibility. In reinforcement learning, this inflexibility can reflect asymmetric learning from outcomes above and below expectations. In alternative frameworks, it reflects perseveration independent of learning. Here, we examine evidence for asymmetric reward-learning in OCD and PG by leveraging model-based functional magnetic resonance imaging (fMRI). Compared with healthy controls (HC), OCD patients exhibited a lower learning rate for worse-than-expected outcomes, which was associated with the attenuated encoding of negative reward prediction errors in the dorsomedial prefrontal cortex and the dorsal striatum. PG patients showed higher and lower learning rates for better- and worse-than-expected outcomes, respectively, accompanied by higher encoding of positive reward prediction errors in the anterior insula than HC. Perseveration did not differ considerably between the patient groups and HC. These findings elucidate the neural computations of reward-learning that are altered in OCD and PG, providing a potential account of behavioural inflexibility in those mental disorders.

Characterizing soundscapes across diverse ecosystems using a universal acoustic feature set.

Natural habitats are being impacted by human pressures at an alarming rate. Monitoring these ecosystem-level changes often requires labor-intensive surveys that are unable to detect rapid or unanticipated environmental changes. Here we have developed a generalizable, data-driven solution to this challenge using eco-acoustic data. We exploited a convolutional neural network to embed soundscapes from a variety of ecosystems into a common acoustic space. In both supervised and unsupervised modes, this allowed us to accurately quantify variation in habitat quality across space and in biodiversity through time. On the scale of seconds, we learned a typical soundscape model that allowed automatic identification of anomalous sounds in playback experiments, providing a potential route for real-time automated detection of irregular environmental behavior including illegal logging and hunting. Our highly generalizable approach, and the common set of features, will enable scientists to unlock previously hidden insights from acoustic data and offers promise as a backbone technology for global collaborative autonomous ecosystem monitoring efforts.

On the intersection between data quality and dynamical modelling of large-scale fMRI signals.

Large-scale dynamics of the brain are routinely modelled using systems of nonlinear dynamical equations that describe the evolution of population-level activity, with distinct neural populations often coupled according to an empirically measured structural connectivity matrix. This modelling approach has been used to generate insights into the neural underpinnings of spontaneous brain dynamics, as recorded with techniques such as resting state functional MRI (fMRI). In fMRI, researchers have many degrees of freedom in the way that they can process the data and recent evidence indicates that the choice of pre-processing steps can have a major effect on empirical estimates of functional connectivity. However, the potential influence of such variations on modelling results are seldom considered. Here we show, using three popular whole-brain dynamical models, that different choices during fMRI preprocessing can dramatically affect model fits and interpretations of findings. Critically, we show that the ability of these models to accurately capture patterns in fMRI dynamics is mostly driven by the degree to which they fit global signals rather than interesting sources of coordinated neural dynamics. We show that widespread deflections can arise from simple global synchronisation. We introduce a simple two-parameter model that captures these fluctuations and performs just as well as more complex, multi-parameter biophysical models. From our combined analyses of data and simulations, we describe benchmarks to evaluate model fit and validity. Although most models are not resilient to denoising, we show that relaxing the approximation of homogeneous neural populations by more explicitly modelling inter-regional effective connectivity can improve model accuracy at the expense of increased model complexity. Our results suggest that many complex biophysical models may be fitting relatively trivial properties of the data, and underscore a need for tighter integration between data quality assurance and model development.

Multimodal gradients across mouse cortex.

The primate cerebral cortex displays a hierarchy that extends from primary sensorimotor to association areas, supporting increasingly integrated function underpinned by a gradient of heterogeneity in the brain's microcircuits. The extent to which these hierarchical gradients are unique to primate or may reflect a conserved mammalian principle of brain organization remains unknown. Here we report the topographic similarity of large-scale gradients in cytoarchitecture, gene expression, interneuron cell densities, and long-range axonal connectivity, which vary from primary sensory to prefrontal areas of mouse cortex, highlighting an underappreciated spatial dimension of mouse cortical specialization. Using the T1-weighted:T2-weighted (T1w:T2w) magnetic resonance imaging map as a common spatial reference for comparison across species, we report interspecies agreement in a range of large-scale cortical gradients, including a significant correspondence between gene transcriptional maps in mouse cortex with their human orthologs in human cortex, as well as notable interspecies differences. Our results support the view of systematic structural variation across cortical areas as a core organizational principle that may underlie hierarchical specialization in mammalian brains.

The correspondence between EMG and EEG measures of changes in cortical excitability following transcranial magnetic stimulation.

KEY POINTS: TMS is commonly used to study excitatory/inhibitory neurotransmission in cortical circuits. Changes in cortical excitability following TMS are typically measured from hand (using EMG; limited to motor cortex) or scalp (using EEG); however, it is unclear whether these two measures represent the same activity when assessing motor cortex. We found that TMS-EMG and TMS-EEG measures of motor cortex excitability are differentially affected by sensory confounds at different time points, masking any actual relationship between them in the time domain. In the frequency domain, local high-frequency oscillations in EEG recordings were minimally confounded by sensory artefacts and demonstrated strong correlations with EMG measures of cortical excitability across time, regardless of TMS intensity or waveform. Therefore, despite the effects of sensory artefacts, the two measures of motor cortex excitability share a response component, suggesting that they index a similar cortical activity and perhaps the same neuronal population. ABSTRACT: Transcranial magnetic stimulation (TMS) is a powerful tool for investigating cortical circuits. Changes in cortical excitability following TMS are typically assessed by measuring changes in either conditioned motor-evoked potentials (MEPs) following paired-pulse TMS over motor cortex or evoked potentials measured with electroencephalography following single-pulse TMS (TEPs). However, it is unclear whether these two measures of cortical excitability index the same cortical response. Twenty-four healthy participants received local and interhemispheric paired-pulse TMS over motor cortex with eight inter-pulse intervals, sub- and suprathreshold conditioning intensities, and two different pulse waveforms, while MEPs were recorded from a hand muscle. TEPs were also recorded in response to single-pulse TMS using the conditioning pulse alone. The relationships between TEPs and conditioned-MEPs were evaluated using metrics sensitive to both their magnitude at each time point and their overall shape across time. The impacts of undesired sensory potentials resulting from TMS pulse and muscle contractions were also assessed on both measures. Both conditioned-MEPs and TEPs were sensitive to re-afferent somatosensory activity following motor-evoked responses, but over different post-stimulus time points. Moreover, the amplitude of low-frequency oscillations in TEPs was strongly correlated with the sensory potentials, whereas early and local high-frequency responses showed minimal relationships. Accordingly, conditioned-MEPs did not correlate with TEPs in the time domain but showed high shape similarity with the amplitude of high-frequency oscillations in TEPs. Therefore, despite the effects of sensory confounds, the TEP and MEP measures share a response component, suggesting that they index a similar cortical response and perhaps the same neuronal populations.

Where the genome meets the connectome: Understanding how genes shape human brain connectivity.

The integration of modern neuroimaging methods with genetically informative designs and data can shed light on the molecular mechanisms underlying the structural and functional organization of the human connectome. Here, we review studies that have investigated the genetic basis of human brain network structure and function through three complementary frameworks: (1) the quantification of phenotypic heritability through classical twin designs; (2) the identification of specific DNA variants linked to phenotypic variation through association and related studies; and (3) the analysis of correlations between spatial variations in imaging phenotypes and gene expression profiles through the integration of neuroimaging and transcriptional atlas data. We consider the basic foundations, strengths, limitations, and discoveries associated with each approach. We present converging evidence to indicate that anatomical connectivity is under stronger genetic influence than functional connectivity and that genetic influences are not uniformly distributed throughout the brain, with phenotypic variation in certain regions and connections being under stronger genetic control than others. We also consider how the combination of imaging and genetics can be used to understand the ways in which genes may drive brain dysfunction in different clinical disorders.

Scaling of gene transcriptional gradients with brain size across mouse development.

The structure of the adult brain is the result of complex physical mechanisms acting in three-dimensional space through development. Consequently, the brain's spatial embedding plays a key role in its organization, including the gradient-like patterning of gene expression that encodes the molecular underpinning of functional specialization. However, we do not yet understand how changes in brain shape and size that occur across development influence the brain's transcriptional architecture. Here we investigate the spatial embedding of transcriptional patterns of over 1800 genes across seven time points through mouse-brain development using data from the Allen Developing Mouse Brain Atlas. We find that transcriptional similarity decreases exponentially with separation distance across all developmental time points, with a correlation length scale that follows a power-law scaling relationship with a linear dimension of brain size. This scaling suggests that the mouse brain achieves a characteristic balance between local molecular similarity (homogeneous gene expression within a specialized brain area) and longer-range diversity (between functionally specialized brain areas) throughout its development. Extrapolating this mouse developmental scaling relationship to the human cortex yields a prediction consistent with the value measured from microarray data. We introduce a simple model of brain growth as spatially autocorrelated gene-expression gradients that expand through development, which captures key features of the mouse developmental data. Complementing the well-known exponential distance rule for structural connectivity, our findings characterize an analogous exponential distance rule for transcriptional gradients that scales across mouse brain development, providing new understanding of spatial constraints on the brain's molecular patterning.

Cortical Excitation:Inhibition Imbalance Causes Abnormal Brain Network Dynamics as Observed in Neurodevelopmental Disorders.

Abnormal brain development manifests itself at different spatial scales. However, whether abnormalities at the cellular level can be diagnosed from network activity measured with functional magnetic resonance imaging (fMRI) is largely unknown, yet of high clinical relevance. Here a putative mechanism reported in neurodevelopmental disorders, that is, excitation-to-inhibition ratio (E:I), was chemogenetically increased within cortical microcircuits of the mouse brain and measured via fMRI. Increased E:I caused a significant "reduction" of long-range connectivity, irrespective of whether excitatory neurons were facilitated or inhibitory Parvalbumin (PV) interneurons were suppressed. Training a classifier on fMRI signals, we were able to accurately classify cortical areas exhibiting increased E:I. This classifier was validated in an independent cohort of Fmr1y/- knockout mice, a model for autism with well-documented loss of parvalbumin neurons and chronic alterations of E:I. Our findings demonstrate a promising novel approach towards inferring microcircuit abnormalities from macroscopic fMRI measurements.

Identifying and removing widespread signal deflections from fMRI data: Rethinking the global signal regression problem.

One of the most controversial procedures in the analysis of resting-state functional magnetic resonance imaging (rsfMRI) data is global signal regression (GSR): the removal, via linear regression, of the mean signal averaged over the entire brain. On one hand, the global mean signal contains variance associated with respiratory, scanner-, and motion-related artifacts, and its removal via GSR improves various quality-control metrics, enhances the anatomical specificity of functional-connectivity patterns, and can increase the behavioral variance explained by such patterns. On the other hand, GSR alters the distribution of regional signal correlations in the brain, can induce artifactual anticorrelations, may remove real neural signal, and can distort case-control comparisons of functional-connectivity measures. Global signal fluctuations can be identified visually from a matrix of colour-coded signal intensities, called a carpet plot, in which rows represent voxels and columns represent time. Prior to GSR, large, periodic bands of coherent signal changes that affect most of the brain are often apparent; after GSR, these apparently global changes are greatly diminished. Here, using three independent datasets, we show that reordering carpet plots to emphasize cluster structure in the data reveals a greater diversity of spatially widespread signal deflections (WSDs) than previously thought. Their precise form varies across time and participants, and GSR is only effective in removing specific kinds of WSDs. We present an alternative, iterative correction method called Diffuse Cluster Estimation and Regression (DiCER), that identifies representative signals associated with large clusters of coherent voxels. DiCER is more effective than GSR at removing diverse WSDs as visualized in carpet plots, reduces correlations between functional connectivity and head-motion estimates, reduces inter-individual variability in global correlation structure, and results in comparable or improved identification of canonical functional-connectivity networks. Using task fMRI data across 47 contrasts from 7 tasks in the Human Connectome Project, we also present evidence that DiCER is more successful than GSR in preserving the spatial structure of expected task-related activation patterns. Our findings indicate that care must be exercised when examining WSDs (and their possible removal) in rsfMRI data, and that DiCER is a viable alternative to GSR for removing anatomically widespread and temporally coherent signals. All code for implementing DiCER and replicating our results is available at https://github.com/BMHLab/DiCER.

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

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

A practical guide to linking brain-wide gene expression and neuroimaging data.

The recent availability of comprehensive, brain-wide gene expression atlases such as the Allen Human Brain Atlas (AHBA) has opened new opportunities for understanding how spatial variations on molecular scale relate to the macroscopic neuroimaging phenotypes. A rapidly growing body of literature is demonstrating relationships between gene expression and diverse properties of brain structure and function, but approaches for combining expression atlas data with neuroimaging are highly inconsistent, with substantial variations in how the expression data are processed. The degree to which these methodological variations affect findings is unclear. Here, we outline a seven-step analysis pipeline for relating brain-wide transcriptomic and neuroimaging data and compare how different processing choices influence the resulting data. We suggest that studies using the AHBA should work towards a unified data processing pipeline to ensure consistent and reproducible results in this burgeoning field.

Hub connectivity, neuronal diversity, and gene expression in the Caenorhabditis elegans connectome.

Studies of nervous system connectivity, in a wide variety of species and at different scales of resolution, have identified several highly conserved motifs of network organization. One such motif is a heterogeneous distribution of connectivity across neural elements, such that some elements act as highly connected and functionally important network hubs. These brain network hubs are also densely interconnected, forming a so-called rich club. Recent work in mouse has identified a distinctive transcriptional signature of neural hubs, characterized by tightly coupled expression of oxidative metabolism genes, with similar genes characterizing macroscale inter-modular hub regions of the human cortex. Here, we sought to determine whether hubs of the neuronal C. elegans connectome also show tightly coupled gene expression. Using open data on the chemical and electrical connectivity of 279 C. elegans neurons, and binary gene expression data for each neuron across 948 genes, we computed a correlated gene expression score for each pair of neurons, providing a measure of their gene expression similarity. We demonstrate that connections between hub neurons are the most similar in their gene expression while connections between nonhubs are the least similar. Genes with the greatest contribution to this effect are involved in glutamatergic and cholinergic signaling, and other communication processes. We further show that coupled expression between hub neurons cannot be explained by their neuronal subtype (i.e., sensory, motor, or interneuron), separation distance, chemically secreted neurotransmitter, birth time, pairwise lineage distance, or their topological module affiliation. Instead, this coupling is intrinsically linked to the identity of most hubs as command interneurons, a specific class of interneurons that regulates locomotion. Our results suggest that neural hubs may possess a distinctive transcriptional signature, preserved across scales and species, that is related to the involvement of hubs in regulating the higher-order behaviors of a given organism.

A transcriptional signature of hub connectivity in the mouse connectome.

Connectivity is not distributed evenly throughout the brain. Instead, it is concentrated on a small number of highly connected neural elements that act as network hubs. Across different species and measurement scales, these hubs show dense interconnectivity, forming a core or "rich club" that integrates information across anatomically distributed neural systems. Here, we show that projections between connectivity hubs of the mouse brain are both central (i.e., they play an important role in neural communication) and costly (i.e., they extend over long anatomical distances) aspects of network organization that carry a distinctive genetic signature. Analyzing the neuronal connectivity of 213 brain regions and the transcriptional coupling, across 17,642 genes, between each pair of regions, we find that coupling is highest for pairs of connected hubs, intermediate for links between hubs and nonhubs, and lowest for connected pairs of nonhubs. The high transcriptional coupling associated with hub connectivity is driven by genes regulating the oxidative synthesis and metabolism of ATP--the primary energetic currency of neuronal communication. This genetic signature contrasts that identified for neuronal connectivity in general, which is driven by genes regulating neuronal, synaptic, and axonal structure and function. Our findings establish a direct link between molecular function and the large-scale topology of neuronal connectivity, showing that brain hubs display a tight coordination of gene expression, often over long anatomical distances, that is intimately related to the metabolic requirements of these highly active network elements.

An evaluation of the efficacy, reliability, and sensitivity of motion correction strategies for resting-state functional MRI.

Estimates of functional connectivity derived from resting-state functional magnetic resonance imaging (rs-fMRI) are sensitive to artefacts caused by in-scanner head motion. This susceptibility has motivated the development of numerous denoising methods designed to mitigate motion-related artefacts. Here, we compare popular retrospective rs-fMRI denoising methods, such as regression of head motion parameters and mean white matter (WM) and cerebrospinal fluid (CSF) (with and without expansion terms), aCompCor, volume censoring (e.g., scrubbing and spike regression), global signal regression and ICA-AROMA, combined into 19 different pipelines. These pipelines were evaluated across five different quality control benchmarks in four independent datasets associated with varying levels of motion. Pipelines were benchmarked by examining the residual relationship between in-scanner movement and functional connectivity after denoising; the effect of distance on this residual relationship; whole-brain differences in functional connectivity between high- and low-motion healthy controls (HC); the temporal degrees of freedom lost during denoising; and the test-retest reliability of functional connectivity estimates. We also compared the sensitivity of each pipeline to clinical differences in functional connectivity in independent samples of people with schizophrenia and obsessive-compulsive disorder. Our results indicate that (1) simple linear regression of regional fMRI time series against head motion parameters and WM/CSF signals (with or without expansion terms) is not sufficient to remove head motion artefacts; (2) aCompCor pipelines may only be viable in low-motion data; (3) volume censoring performs well at minimising motion-related artefact but a major benefit of this approach derives from the exclusion of high-motion individuals; (4) while not as effective as volume censoring, ICA-AROMA performed well across our benchmarks for relatively low cost in terms of data loss; (5) the addition of global signal regression improved the performance of nearly all pipelines on most benchmarks, but exacerbated the distance-dependence of correlations between motion and functional connectivity; and (6) group comparisons in functional connectivity between healthy controls and schizophrenia patients are highly dependent on preprocessing strategy. We offer some recommendations for best practice and outline simple analyses to facilitate transparent reporting of the degree to which a given set of findings may be affected by motion-related artefact.

Predicting Depression From Language-Based Emotion Dynamics: Longitudinal Analysis of Facebook and Twitter Status Updates.

BACKGROUND: Frequent expression of negative emotion words on social media has been linked to depression. However, metrics have relied on average values, not dynamic measures of emotional volatility. OBJECTIVE: The aim of this study was to report on the associations between depression severity and the variability (time-unstructured) and instability (time-structured) in emotion word expression on Facebook and Twitter across status updates. METHODS: Status updates and depression severity ratings of 29 Facebook users and 49 Twitter users were collected through the app MoodPrism. The average proportion of positive and negative emotion words used, within-person variability, and instability were computed. RESULTS: Negative emotion word instability was a significant predictor of greater depression severity on Facebook (rs(29)=.44, P=.02, 95% CI 0.09-0.69), even after controlling for the average proportion of negative emotion words used (partial rs(26)=.51, P=.006) and within-person variability (partial rs(26)=.49, P=.009). A different pattern emerged on Twitter where greater negative emotion word variability indicated lower depression severity (rs(49)=-.34, P=.01, 95% CI -0.58 to 0.09). Differences between Facebook and Twitter users in their emotion word patterns and psychological characteristics were also explored. CONCLUSIONS: The findings suggest that negative emotion word instability may be a simple yet sensitive measure of time-structured variability, useful when screening for depression through social media, though its usefulness may depend on the social media platform.

Structural connectome topology relates to regional BOLD signal dynamics in the mouse brain.

Brain dynamics are thought to unfold on a network determined by the pattern of axonal connections linking pairs of neuronal elements; the so-called connectome. Prior work has indicated that structural brain connectivity constrains pairwise correlations of brain dynamics ("functional connectivity"), but it is not known whether inter-regional axonal connectivity is related to the intrinsic dynamics of individual brain areas. Here we investigate this relationship using a weighted, directed mesoscale mouse connectome from the Allen Mouse Brain Connectivity Atlas and resting state functional MRI (rs-fMRI) time-series data measured in 184 brain regions in eighteen anesthetized mice. For each brain region, we measured degree, betweenness, and clustering coefficient from weighted and unweighted, and directed and undirected versions of the connectome. We then characterized the univariate rs-fMRI dynamics in each brain region by computing 6930 time-series properties using the time-series analysis toolbox, hctsa. After correcting for regional volume variations, strong and robust correlations between structural connectivity properties and rs-fMRI dynamics were found only when edge weights were accounted for, and were associated with variations in the autocorrelation properties of the rs-fMRI signal. The strongest relationships were found for weighted in-degree, which was positively correlated to the autocorrelation of fMRI time series at time lag τ = 34 s (partial Spearman correlation ρ=0.58), as well as a range of related measures such as relative high frequency power (f > 0.4 Hz: ρ=-0.43). Our results indicate that the topology of inter-regional axonal connections of the mouse brain is closely related to intrinsic, spontaneous dynamics such that regions with a greater aggregate strength of incoming projections display longer timescales of activity fluctuations.

Developmental Changes in Brain Network Hub Connectivity in Late Adolescence.

The human brain undergoes substantial development throughout adolescence and into early adulthood. This maturational process is thought to include the refinement of connectivity between putative connectivity hub regions of the brain, which collectively form a dense core that enhances the functional integration of anatomically distributed, and functionally specialized, neural systems. Here, we used longitudinal diffusion magnetic resonance imaging to characterize changes in connectivity between 80 cortical and subcortical anatomical regions over a 2 year period in 31 adolescents between the ages of 15 and 19 years. Connectome-wide analysis indicated that only a small subset of connections showed evidence of statistically significant developmental change over the study period, with 8% and 6% of connections demonstrating decreased and increased structural connectivity, respectively. Nonetheless, these connections linked 93% and 90% of the 80 regions, respectively, pointing to a selective, yet anatomically distributed pattern of developmental changes that involves most of the brain. Hub regions showed a distinct tendency to be highly connected to each other, indicating robust "rich-club" organization. Moreover, connectivity between hubs was disproportionately influenced by development, such that connectivity between subcortical hubs decreased over time, whereas frontal-subcortical and frontal-parietal hub-hub connectivity increased over time. These findings suggest that late adolescence is characterized by selective, yet significant remodeling of hub-hub connectivity, with the topological organization of hubs shifting emphasis from subcortical hubs in favor of an increasingly prominent role for frontal hub regions.

Extracting interpretable signatures of whole-brain dynamics through systematic comparison.

The brain's complex distributed dynamics are typically quantified using a limited set of manually selected statistical properties, leaving the possibility that alternative dynamical properties may outperform those reported for a given application. Here, we address this limitation by systematically comparing diverse, interpretable features of both intra-regional activity and inter-regional functional coupling from resting-state functional magnetic resonance imaging (rs-fMRI) data, demonstrating our method using case-control comparisons of four neuropsychiatric disorders. Our findings generally support the use of linear time-series analysis techniques for rs-fMRI case-control analyses, while also identifying new ways to quantify informative dynamical fMRI structures. While simple statistical representations of fMRI dynamics performed surprisingly well (e.g., properties within a single brain region), combining intra-regional properties with inter-regional coupling generally improved performance, underscoring the distributed, multifaceted changes to fMRI dynamics in neuropsychiatric disorders. The comprehensive, data-driven method introduced here enables systematic identification and interpretation of quantitative dynamical signatures of multivariate time-series data, with applicability beyond neuroimaging to diverse scientific problems involving complex time-varying systems.

Uncovering a stability signature of brain dynamics associated with meditation experience using massive time-series feature extraction.

Previous research has examined resting electroencephalographic (EEG) data to explore brain activity related to meditation. However, previous research has mostly examined power in different frequency bands. The practical objective of this study was to comprehensively test whether other types of time-series analysis methods are better suited to characterize brain activity related to meditation. To achieve this, we compared >7000 time-series features of the EEG signal to comprehensively characterize brain activity differences in meditators, using many measures that are novel in meditation research. Eyes-closed resting-state EEG data from 49 meditators and 46 non-meditators was decomposed into the top eight principal components (PCs). We extracted 7381 time-series features from each PC and each participant and used them to train classification algorithms to identify meditators. Highly differentiating individual features from successful classifiers were analysed in detail. Only the third PC (which had a central-parietal maximum) showed above-chance classification accuracy (67 %, pFDR = 0.007), for which 405 features significantly distinguished meditators (all pFDR < 0.05). Top-performing features indicated that meditators exhibited more consistent statistical properties across shorter subsegments of their EEG time-series (higher stationarity) and displayed an altered distributional shape of values about the mean. By contrast, classifiers trained with traditional band-power measures did not distinguish the groups (pFDR > 0.05). Our novel analysis approach suggests the key signatures of meditators' brain activity are higher temporal stability and a distribution of time-series values suggestive of longer, larger, or more frequent non-outlying voltage deviations from the mean within the third PC of their EEG data. The higher temporal stability observed in this EEG component might underpin the higher attentional stability associated with meditation. The novel time-series properties identified here have considerable potential for future exploration in meditation research and the analysis of neural dynamics more broadly.