Neuronal timescales are functionally dynamic and shaped by cortical microarchitecture

Neurophysiological signatures of cortical micro-architecture

Systematic spatial variation in micro-architecture is observed across the cortex. These micro-architectural gradients are reflected in neural activity, which can be captured by neurophysiological time-series. How spontaneous neurophysiological dynamics are organized across the cortex and how they arise from heterogeneous cortical micro-architecture remains unknown. Here we extensively profile regional neurophysiological dynamics across the human brain by estimating over 6 800 timeseries features from the resting state magnetoencephalography (MEG) signal. We then map regional time-series profiles to a comprehensive multi-modal, multi-scale atlas of cortical micro-architecture, including microstructure, metabolism, neurotransmitter receptors, cell types and laminar differentiation. We find that the dominant axis of neurophysiological dynamics reflects characteristics of power spectrum density and linear correlation structure of the signal, emphasizing the importance of conventional features of electromagnetic dynamics while identifying additional informative features that have traditionally received less attention. Moreover, spatial variation in neurophysiological dynamics is colocalized with multiple micro-architectural features, including genomic gradients, intracortical myelin, neurotransmitter receptors and transporters, and oxygen and glucose metabolism. Collectively, this work opens new avenues for studying the anatomical basis of neural activity.

Spatial and temporal correlations in neural networks with structured connectivity

Correlated fluctuations in the activity of neural populations reflect the network's dynamics and connectivity. The temporal and spatial dimensions of neural correlations are interdependent. However, prior theoretical work mainly analyzed correlations in either spatial or temporal domains, oblivious to their interplay. We show that the network dynamics and connectivity jointly define the spatiotemporal profile of neural correlations. We derive analytical expressions for pairwise correlations in networks of binary units with spatially arranged connectivity in one and two dimensions. We find that spatial interactions among units generate multiple timescales in auto- and cross-correlations. Each timescale is associated with fluctuations at a particular spatial frequency, making a hierarchical contribution to the correlations. External inputs can modulate the correlation timescales when spatial interactions are nonlinear, and the modulation effect depends on the operating regime of network dynamics. These theoretical results open new ways to relate connectivity and dynamics in cortical networks via measurements of spatiotemporal neural correlations.

Dynamical interactions reconfigure the gradient of cortical timescales

The functional organization of the brain is usually presented with a back-to-front gradient of timescales, reflecting regional specialization with sensory areas (back) processing information faster than associative areas (front), which perform information integration. However, cognitive processes require not only local information processing but also coordinated activity across regions. Using magnetoencephalography recordings, we find that the functional connectivity at the edge level (between two regions) is also characterized by a back-to-front gradient of timescales following that of the regional gradient. Unexpectedly, we demonstrate a reverse front-to-back gradient when nonlocal interactions are prominent. Thus, the timescales are dynamic and can switch between back-to-front and front-to-back patterns.

Asymmetric signaling across the hierarchy of cytoarchitecture within the human connectome

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

Correspondence between gene expression and neurotransmitter receptor and transporter density in the human brain

Neurotransmitter receptors modulate signaling between neurons. Thus, neurotransmitter receptors and transporters play a key role in shaping brain function. Due to the lack of comprehensive neurotransmitter receptor/transporter density datasets, microarray gene expression measuring mRNA transcripts is often used as a proxy for receptor densities. In the present report, we comprehensively test the spatial correlation between gene expression and protein density for a total of 27 neurotransmitter receptors, receptor binding-sites, and transporters across 9 different neurotransmitter systems, using both PET and autoradiography radioligand-based imaging modalities. We find poor spatial correspondences between gene expression and density for all neurotransmitter receptors and transporters except four single-protein metabotropic receptors (5-HT_1A, CB₁, D₂, and MOR). These expression-density associations are related to gene differential stability and can vary between cortical and subcortical structures. Altogether, we recommend using direct measures of receptor and transporter density when relating neurotransmitter systems to brain structure and function.

Robust estimation of 1/f activity improves oscillatory burst detection

Neural oscillations often occur as transient bursts with variable amplitude and frequency dynamics. Quantifying these effects is important for understanding brain-behaviour relationships, especially in continuous datasets. To robustly measure bursts, rhythmical periods of oscillatory activity must be separated from arrhythmical background 1/f activity, which is ubiquitous in electrophysiological recordings. The Better OSCillation (BOSC) framework achieves this by defining a power threshold above the estimated background 1/f activity, combined with a duration threshold. Here we introduce a modification to this approach called fBOSC, which uses a spectral parametrisation tool to accurately model background 1/f activity in neural data. fBOSC (which is openly available as a MATLAB toolbox) is robust to power spectra with oscillatory peaks and can also model non-linear spectra. Through a series of simulations, we show that fBOSC more accurately models the 1/f power spectrum compared with existing methods. fBOSC was especially beneficial where power spectra contained a 'knee' below ~.5-10 Hz, which is typical in neural data. We also found that, unlike other methods, fBOSC was unaffected by oscillatory peaks in the neural power spectrum. Moreover, by robustly modelling background 1/f activity, the sensitivity for detecting oscillatory bursts was standardised across frequencies (e.g., theta- and alpha-bands). Finally, using openly available resting state magnetoencephalography and intracranial electrophysiology datasets, we demonstrate the application of fBOSC for oscillatory burst detection in the theta-band. These simulations and empirical analyses highlight the value of fBOSC in detecting oscillatory bursts, including in datasets that are long and continuous with no distinct experimental trials.

Mapping neurotransmitter systems to the structural and functional organization of the human neocortex

Neurotransmitter receptors support the propagation of signals in the human brain. How receptor systems are situated within macro-scale neuroanatomy and how they shape emergent function remain poorly understood, and there exists no comprehensive atlas of receptors. Here we collate positron emission tomography data from more than 1,200 healthy individuals to construct a whole-brain three-dimensional normative atlas of 19 receptors and transporters across nine different neurotransmitter systems. We found that receptor profiles align with structural connectivity and mediate function, including neurophysiological oscillatory dynamics and resting-state hemodynamic functional connectivity. Using the Neurosynth cognitive atlas, we uncovered a topographic gradient of overlapping receptor distributions that separates extrinsic and intrinsic psychological processes. Finally, we found both expected and novel associations between receptor distributions and cortical abnormality patterns across 13 disorders. We replicated all findings in an independently collected autoradiography dataset. This work demonstrates how chemoarchitecture shapes brain structure and function, providing a new direction for studying multi-scale brain organization.

neuromaps: structural and functional interpretation of brain maps

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

Advances in human intracranial electroencephalography research, guidelines and good practices

Since the second-half of the twentieth century, intracranial electroencephalography (iEEG), including both electrocorticography (ECoG) and stereo-electroencephalography (sEEG), has provided an intimate view into the human brain. At the interface between fundamental research and the clinic, iEEG provides both high temporal resolution and high spatial specificity but comes with constraints, such as the individual's tailored sparsity of electrode sampling. Over the years, researchers in neuroscience developed their practices to make the most of the iEEG approach. Here we offer a critical review of iEEG research practices in a didactic framework for newcomers, as well addressing issues encountered by proficient researchers. The scope is threefold: (i) review common practices in iEEG research, (ii) suggest potential guidelines for working with iEEG data and answer frequently asked questions based on the most widespread practices, and (iii) based on current neurophysiological knowledge and methodologies, pave the way to good practice standards in iEEG research. The organization of this paper follows the steps of iEEG data processing. The first section contextualizes iEEG data collection. The second section focuses on localization of intracranial electrodes. The third section highlights the main pre-processing steps. The fourth section presents iEEG signal analysis methods. The fifth section discusses statistical approaches. The sixth section draws some unique perspectives on iEEG research. Finally, to ensure a consistent nomenclature throughout the manuscript and to align with other guidelines, e.g., Brain Imaging Data Structure (BIDS) and the OHBM Committee on Best Practices in Data Analysis and Sharing (COBIDAS), we provide a glossary to disambiguate terms related to iEEG research.

Thalamocortical contribution to flexible learning in neural systems

Animal brains evolved to optimize behavior in dynamic environments, flexibly selecting actions that maximize future rewards in different contexts. A large body of experimental work indicates that such optimization changes the wiring of neural circuits, appropriately mapping environmental input onto behavioral outputs. A major unsolved scientific question is how optimal wiring adjustments, which must target the connections responsible for rewards, can be accomplished when the relation between sensory inputs, action taken, and environmental context with rewards is ambiguous. The credit assignment problem can be categorized into context-independent structural credit assignment and context-dependent continual learning. In this perspective, we survey prior approaches to these two problems and advance the notion that the brain's specialized neural architectures provide efficient solutions. Within this framework, the thalamus with its cortical and basal ganglia interactions serves as a systems-level solution to credit assignment. Specifically, we propose that thalamocortical interaction is the locus of meta-learning where the thalamus provides cortical control functions that parametrize the cortical activity association space. By selecting among these control functions, the basal ganglia hierarchically guide thalamocortical plasticity across two timescales to enable meta-learning. The faster timescale establishes contextual associations to enable behavioral flexibility, while the slower one enables generalization to new contexts.

Time-resolved parameterization of aperiodic and periodic brain activity

Macroscopic neural dynamics comprise both aperiodic and periodic signal components. Recent advances in parameterizing neural power spectra offer practical tools for evaluating these features separately. Although neural signals vary dynamically and express non-stationarity in relation to ongoing behaviour and perception, current methods yield static spectral decompositions. Here, we introduce Spectral Parameterization Resolved in Time (SPRiNT) as a novel method for decomposing complex neural dynamics into periodic and aperiodic spectral elements in a time-resolved manner. First, we demonstrate, with naturalistic synthetic data, SPRiNT's capacity to reliably recover time-varying spectral features. We emphasize SPRiNT's specific strengths compared to other time-frequency parameterization approaches based on wavelets. Second, we use SPRiNT to illustrate how aperiodic spectral features fluctuate across time in empirical resting-state EEG data (n=178) and relate the observed changes in aperiodic parameters over time to participants' demographics and behaviour. Lastly, we use SPRiNT to demonstrate how aperiodic dynamics relate to movement behaviour in intracranial recordings in rodents. We foresee SPRiNT responding to growing neuroscientific interests in the parameterization of time-varying neural power spectra and advancing the quantitation of complex neural dynamics at the natural time scales of behaviour.

Sources of Variation in the Spectral Slope of the Sleep EEG

The 1/f spectral slope of the electroencephalogram (EEG) estimated in the γ frequency range has been proposed as an arousal marker that differentiates wake, nonrapid eye movement (NREM) sleep, and rapid eye movement (REM) sleep. Here, we sought to replicate and extend these findings in a large sample, providing a comprehensive characterization of how slope changes with age, sex, and its test-retest reliability as well as potential confounds that could affect the slope estimation. We used 10,255 whole-night polysomnograms (PSGs) from the National Sleep Research Resource (NSRR). All preprocessing steps were performed using an open-source Luna package and the spectral slope was estimated by fitting log-log linear regression models on the absolute power from 30 to 45 Hz separately for wake, NREM, and REM stages. We confirmed that the mean spectral slope grows steeper going from wake to NREM to REM sleep. We found that the choice of mastoid referencing scheme modulated the extent to which electromyogenic, or electrocardiographic artifacts were likely to bias 30- to 45-Hz slope estimates, as well as other sources of technical, device-specific bias. Nonetheless, within individuals, slope estimates were relatively stable over time. Both cross-sectionally and longitudinal, slopes tended to become shallower with increasing age, particularly for REM sleep; males tended to show flatter slopes than females across all states. Our findings support that spectral slope can be a valuable arousal marker for both clinical and research endeavors but also underscore the importance of considering interindividual variation and multiple methodological aspects related to its estimation.

Empirical transmit field bias correction of T1w/T2w myelin maps

T1-weighted divided by T2-weighted (T1w/T2w) myelin maps were initially developed for neuroanatomical analyses such as identifying cortical areas, but they are increasingly used in statistical comparisons across individuals and groups with other variables of interest. Existing T1w/T2w myelin maps contain radiofrequency transmit field (B1+) biases, which may be correlated with these variables of interest, leading to potentially spurious results. Here we propose two empirical methods for correcting these transmit field biases using either explicit measures of the transmit field or alternatively a 'pseudo-transmit' approach that is highly correlated with the transmit field at 3T. We find that the resulting corrected T1w/T2w myelin maps are both better neuroanatomical measures (e.g., for use in cross-species comparisons), and more appropriate for statistical comparisons of relative T1w/T2w differences across individuals and groups (e.g., sex, age, or body-mass-index) within a consistently acquired study at 3T. We recommend that investigators who use the T1w/T2w approach for mapping cortical myelin use these B1+ transmit field corrected myelin maps going forward.

Edge-centric analysis of time-varying functional brain networks with applications in autism spectrum disorder

The interaction between brain regions changes over time, which can be characterized using time-varying functional connectivity (tvFC). The common approach to estimate tvFC uses sliding windows and offers limited temporal resolution. An alternative method is to use the recently proposed edge-centric approach, which enables the tracking of moment-to-moment changes in co-fluctuation patterns between pairs of brain regions. Here, we first examined the dynamic features of edge time series and compared them to those in the sliding window tvFC (sw-tvFC). Then, we used edge time series to compare subjects with autism spectrum disorder (ASD) and healthy controls (CN). Our results indicate that relative to sw-tvFC, edge time series captured rapid and bursty network-level fluctuations that synchronize across subjects during movie-watching. The results from the second part of the study suggested that the magnitude of peak amplitude in the collective co-fluctuations of brain regions (estimated as root sum square (RSS) of edge time series) is similar in CN and ASD. However, the trough-to-trough duration in RSS signal is greater in ASD, compared to CN. Furthermore, an edge-wise comparison of high-amplitude co-fluctuations showed that the within-network edges exhibited greater magnitude fluctuations in CN. Our findings suggest that high-amplitude co-fluctuations captured by edge time series provide details about the disruption of functional brain dynamics that could potentially be used in developing new biomarkers of mental disorders.

A personalizable autonomous neural mass model of epileptic seizures

Work in the last two decades has shown that neural mass models (NMM) can realistically reproduce and explain epileptic seizure transitions as recorded by electrophysiological methods (EEG, SEEG). In previous work, advances were achieved by increasing excitation and heuristically varying network inhibitory coupling parameters in the models. Based on these early studies, we provide a laminar NMM capable of realistically reproducing the electrical activity recorded by SEEG in the epileptogenic zone during interictal to ictal states. With the exception of the external noise input into the pyramidal cell population, the model dynamics are autonomous. By setting the system at a point close to bifurcation, seizure-like transitions are generated, including pre-ictal spikes, low voltage fast activity, and ictal rhythmic activity. A novel element in the model is a physiologically motivated algorithm for chloride dynamics: the gain of GABAergic post-synaptic potentials is modulated by the pathological accumulation of chloride in pyramidal cells due to high inhibitory input and/or dysfunctional chloride transport. In addition, in order to simulate SEEG signals for comparison with real seizure recordings, the NMM is embedded first in a layered model of the neocortex and then in a realistic physical model. We compare modeling results with data from four epilepsy patient cases. By including key pathophysiological mechanisms, the proposed framework captures succinctly the electrophysiological phenomenology observed in ictal states, paving the way for robust personalization methods based on NMMs.

Whole-brain neuronal MCT2 lactate transporter expression links metabolism to human brain structure and function

Brain activity is constrained by local availability of chemical energy, which is generated through compartmentalized metabolic processes. By analyzing data of whole human brain gene expression, we characterize the spatial distribution of seven glucose and monocarboxylate membrane transporters that mediate astrocyte–neuron lactate shuttle transfer of energy. We found that the gene coding for neuronal MCT2 is the only gene enriched in cerebral cortex where its abundance is inversely correlated with cortical thickness. Coexpression network analysis revealed that MCT2 was the only gene participating in an organized gene cluster enriched in K+ dynamics. Indeed, the expression of KATP subunits, which mediate lactate increases with spiking activity, is spatially coupled to MCT2 distribution. Notably, MCT2 expression correlated with fluorodeoxyglucose positron emission tomography task-dependent glucose utilization. Finally, the MCT2 messenger RNA gradient closely overlaps with functional MRI brain regions associated with attention, arousal, and stress. Our results highlight neuronal MCT2 lactate transporter as a key component of the cross-talk between astrocytes and neurons and a link between metabolism, cortical structure, and state-dependent brain function.

Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex

Whole-brain neural communication is typically estimated from statistical associations among electromagnetic or haemodynamic time-series. The relationship between functional network architectures recovered from these 2 types of neural activity remains unknown. Here, we map electromagnetic networks (measured using magnetoencephalography (MEG)) to haemodynamic networks (measured using functional magnetic resonance imaging (fMRI)). We find that the relationship between the 2 modalities is regionally heterogeneous and systematically follows the cortical hierarchy, with close correspondence in unimodal cortex and poor correspondence in transmodal cortex. Comparison with the BigBrain histological atlas reveals that electromagnetic-haemodynamic coupling is driven by laminar differentiation and neuron density, suggesting that the mapping between the 2 modalities can be explained by cytoarchitectural variation. Importantly, haemodynamic connectivity cannot be explained by electromagnetic activity in a single frequency band, but rather arises from the mixing of multiple neurophysiological rhythms. Correspondence between the two is largely driven by MEG functional connectivity at the beta (15 to 29 Hz) frequency band. Collectively, these findings demonstrate highly organized but only partly overlapping patterns of connectivity in MEG and fMRI functional networks, opening fundamentally new avenues for studying the relationship between cortical microarchitecture and multimodal connectivity patterns.

Closing the mechanistic gap: the value of microarchitecture in understanding cognitive networks

Cognitive neuroscience aims to provide biologically relevant accounts of cognition. Contemporary research linking spatial patterns of neural activity to psychological constructs describes 'where' hypothesised functions occur, but not 'how' these regions contribute to cognition. Technological, empirical, and conceptual advances allow this mechanistic gap to be closed by embedding patterns of functional activity in macro- and microscale descriptions of brain organisation. Recent work on the default mode network (DMN) and the multiple demand network (MDN), for example, highlights a microarchitectural landscape that may explain how activity in these networks integrates varied information, thus providing an anatomical foundation that will help to explain how these networks contribute to many different cognitive states. This perspective highlights emerging insights into how microarchitecture can constrain network accounts of human cognition.

Multiple timescales of sensory-evidence accumulation across the dorsal cortex

Cortical areas seem to form a hierarchy of intrinsic timescales, but the relevance of this organization for cognitive behavior remains unknown. In particular, decisions requiring the gradual accrual of sensory evidence over time recruit widespread areas across this hierarchy. Here, we tested the hypothesis that this recruitment is related to the intrinsic integration timescales of these widespread areas. We trained mice to accumulate evidence over seconds while navigating in virtual reality and optogenetically silenced the activity of many cortical areas during different brief trial epochs. We found that the inactivation of all tested areas affected the evidence-accumulation computation. Specifically, we observed distinct changes in the weighting of sensory evidence occurring during and before silencing, such that frontal inactivations led to stronger deficits on long timescales than posterior cortical ones. Inactivation of a subset of frontal areas also led to moderate effects on behavioral processes beyond evidence accumulation. Moreover, large-scale cortical Ca2+ activity during task performance displayed different temporal integration windows. Our findings suggest that the intrinsic timescale hierarchy of distributed cortical areas is an important component of evidence-accumulation mechanisms.

Methodological considerations for studying neural oscillations

Neural oscillations are ubiquitous across recording methodologies and species, broadly associated with cognitive tasks, and amenable to computational modelling that investigates neural circuit generating mechanisms and neural population dynamics. Because of this, neural oscillations offer an exciting potential opportunity for linking theory, physiology and mechanisms of cognition. However, despite their prevalence, there are many concerns-new and old-about how our analysis assumptions are violated by known properties of field potential data. For investigations of neural oscillations to be properly interpreted, and ultimately developed into mechanistic theories, it is necessary to carefully consider the underlying assumptions of the methods we employ. Here, we discuss seven methodological considerations for analysing neural oscillations. The considerations are to (1) verify the presence of oscillations, as they may be absent; (2) validate oscillation band definitions, to address variable peak frequencies; (3) account for concurrent non-oscillatory aperiodic activity, which might otherwise confound measures; measure and account for (4) temporal variability and (5) waveform shape of neural oscillations, which are often bursty and/or nonsinusoidal, potentially leading to spurious results; (6) separate spatially overlapping rhythms, which may interfere with each other; and (7) consider the required signal-to-noise ratio for obtaining reliable estimates. For each topic, we provide relevant examples, demonstrate potential errors of interpretation, and offer suggestions to address these issues. We primarily focus on univariate measures, such as power and phase estimates, though we discuss how these issues can propagate to multivariate measures. These considerations and recommendations offer a helpful guide for measuring and interpreting neural oscillations.

Temporal continuity of self: Long autocorrelation windows mediate self-specificity

The self is characterized by an intrinsic temporal component consisting in continuity across time. On the neural level, this temporal continuity manifests in the brain's intrinsic neural timescales (INT) that can be measured by the autocorrelation window (ACW). Recent EEG studies reveal a relationship between resting state ACW and self-consciousness. However, it remains unclear whether ACW exhibits different degrees of task-related changes during self-specific compared to non-self-specific activities. To this end, participants in our study initially recorded an eight-minute autobiographical narrative. Following a resting-state session, participants were presented with their own narrative and the narrative of a stranger while undergoing concurrent EEG recording. Behaviorally, subjects evaluated both of the narratives and indicated their perceptions of positivity or negativity on a moment-to-moment basis by positioning a cursor relative to the center of the computer screen. Our results indicate: (a) greater spatial extension and velocity in the behavioral cursor movement during the self narrative assessment compared to the non-self narrative assessment; and (b) longer neural ACWs in response to the self- compared to the non-self narrative and rest. These findings demonstrate the importance of longer temporal windows in neural activity measured by ACWs for self-specificity. More broadly, the results highlight the relevance of temporal continuity for the self on the neural level. Such temporal continuity may, correspondingly, also manifest on the psychological level as a "common currency" between brain and self.

Spectral parameterization for studying neurodevelopment: How and why

A growing body of literature suggests that the explicit parameterization of neural power spectra is important for the appropriate physiological interpretation of periodic and aperiodic electroencephalogram (EEG) activity. In this paper, we discuss why parameterization is an imperative step for developmental cognitive neuroscientists interested in cognition and behavior across the lifespan, as well as how parameterization can be readily accomplished with an automated spectral parameterization ("specparam") algorithm (Donoghue et al., 2020a). We provide annotated code for power spectral parameterization, via specparam, in Jupyter Notebook and R Studio. We then apply this algorithm to EEG data in childhood (N = 60; Mage = 9.97, SD = 0.95) to illustrate its utility for developmental cognitive neuroscientists. Ultimately, the explicit parameterization of EEG power spectra may help us refine our understanding of how dynamic neural communication contributes to normative and aberrant cognition across the lifespan. Data and annotated analysis code for this manuscript are available on GitHub as a supplement to the open-access specparam toolbox.

Cortical Cross-Frequency Coupling Is Affected by in utero Exposure to Antidepressant Medication

Up to five percent of human infants are exposed to maternal antidepressant medication by serotonin reuptake inhibitors (SRI) during pregnancy, yet the SRI effects on infants’ early neurodevelopment are not fully understood. Here, we studied how maternal SRI medication affects cortical frequency-specific and cross-frequency interactions estimated, respectively, by phase-phase correlations (PPC) and phase-amplitude coupling (PAC) in electroencephalographic (EEG) recordings. We examined the cortical activity in infants after fetal exposure to SRIs relative to a control group of infants without medical history of any kind. Our findings show that the sleep-related dynamics of PPC networks are selectively affected by in utero SRI exposure, however, those alterations do not correlate to later neurocognitive development as tested by neuropsychological evaluation at two years of age. In turn, phase-amplitude coupling was found to be suppressed in SRI infants across multiple distributed cortical regions and these effects were linked to their neurocognitive outcomes. Our results are compatible with the overall notion that in utero drug exposures may cause subtle, yet measurable changes in the brain structure and function. Our present findings are based on the measures of local and inter-areal neuronal interactions in the cortex which can be readily used across species, as well as between different scales of inspection: from the whole animals to in vitro preparations. Therefore, this work opens a framework to explore the cellular and molecular mechanisms underlying neurodevelopmental SRI effects at all translational levels.

Intrinsic timescales as an organizational principle of neural processing across the whole rhesus macaque brain

Hierarchical temporal dynamics are a fundamental computational property of the brain; however, there are no whole brain, noninvasive investigations into timescales of neural processing in animal models. To that end, we used the spatial resolution and sensitivity of ultrahigh field functional magnetic resonance imaging (fMRI) performed at 10.5 T to probe timescales across the whole macaque brain. We uncovered within-species consistency between timescales estimated from fMRI and electrophysiology. Crucially, we extended existing electrophysiological hierarchies to whole-brain topographies. Our results validate the complementary use of hemodynamic and electrophysiological intrinsic timescales, establishing a basis for future translational work. Further, with these results in hand, we were able to show that one facet of the high-dimensional functional connectivity (FC) topography of any region in the brain is closely related to hierarchical temporal dynamics. We demonstrated that intrinsic timescales are organized along spatial gradients that closely match FC gradient topographies across the whole brain. We conclude that intrinsic timescales are a unifying organizational principle of neural processing across the whole brain.

A flexible Bayesian framework for unbiased estimation of timescales

Timescales characterize the pace of change for many dynamic processes in nature. Timescales are usually estimated by fitting the exponential decay of data autocorrelation in the time or frequency domain. Here we show that this standard procedure often fails to recover the correct timescales due to a statistical bias arising from the finite sample size. We develop an alternative approach which estimates timescales by fitting the sample autocorrelation or power spectrum with a generative model based on a mixture of Ornstein-Uhlenbeck (OU) processes using adaptive approximate Bayesian computations (aABC). Our method accounts for finite sample size and noise in data and returns a posterior distribution of timescales that quantifies the estimation uncertainty and can be used for model selection. We demonstrate the accuracy of our method on synthetic data and illustrate its application to recordings from primate cortex. We provide a customizable Python package implementing our framework with different generative models suitable for diverse applications.

Power signatures of habenular neuronal signals in patients with bipolar or unipolar depressive disorders correlate with their disease severity

The habenula is an epithalamic structure implicated in negative reward mechanisms and plays a downstream modulatory role in regulation of dopaminergic and serotonergic functions. Human and animal studies show its hyperactivity in depression which is curtailed by the antidepressant response of ketamine. Deep brain stimulation of habenula (DBS) for major depression have also shown promising results. However, direct neuronal activity of habenula in human studies have rarely been reported. Here, in a cross-sectional design, we acquired both spontaneous resting state and emotional task-induced neuronal recordings from habenula from treatment resistant depressed patients undergoing DBS surgery. We first characterise the aperiodic component (1/f slope) of the power spectrum, interpreted to signify excitation-inhibition balance, in resting and task state. This aperiodicity for left habenula correlated between rest and task and which was significantly positively correlated with depression severity. Time-frequency responses to the emotional picture viewing task show condition differences in beta and gamma frequencies for left habenula and alpha for right habenula. Notably, alpha activity for right habenula was negatively correlated with depression severity. Overall, from direct habenular recordings, we thus show findings convergent with depression models of aberrant excitatory glutamatergic output of the habenula driving inhibition of monoaminergic systems.

Hierarchical timescales in the neocortex: Mathematical mechanism and biological insights

A cardinal feature of the neocortex is the progressive increase of the spatial receptive fields along the cortical hierarchy. Recently, theoretical and experimental findings have shown that the temporal response windows also gradually enlarge, so that early sensory neural circuits operate on short timescales whereas higher-association areas are capable of integrating information over a long period of time. While an increased receptive field is accounted for by spatial summation of inputs from neurons in an upstream area, the emergence of timescale hierarchy cannot be readily explained, especially given the dense interareal cortical connectivity known in the modern connectome. To uncover the required neurobiological properties, we carried out a rigorous analysis of an anatomically based large-scale cortex model of macaque monkeys. Using a perturbation method, we show that the segregation of disparate timescales is defined in terms of the localization of eigenvectors of the connectivity matrix, which depends on three circuit properties: 1) a macroscopic gradient of synaptic excitation, 2) distinct electrophysiological properties between excitatory and inhibitory neuronal populations, and 3) a detailed balance between long-range excitatory inputs and local inhibitory inputs for each area-to-area pathway. Our work thus provides a quantitative understanding of the mechanism underlying the emergence of timescale hierarchy in large-scale primate cortical networks.

Intrinsic neural timescales: temporal integration and segregation

We are continuously bombarded by external inputs of various timescales from the environment. How does the brain process this multitude of timescales? Recent resting state studies show a hierarchy of intrinsic neural timescales (INT) with a shorter duration in unimodal regions (e.g., visual cortex and auditory cortex) and a longer duration in transmodal regions (e.g., default mode network). This unimodal-transmodal hierarchy is present across acquisition modalities [electroencephalogram (EEG)/magnetoencephalogram (MEG) and fMRI] and can be found in different species and during a variety of different task states. Together, this suggests that the hierarchy of INT is central to the temporal integration (combining successive stimuli) and segregation (separating successive stimuli) of external inputs from the environment, leading to temporal segmentation and prediction in perception and cognition.

Resting EEG spectral slopes are associated with age-related differences in information processing speed

BACKGROUND

Previous research has shown the slope of the EEG power spectrum differentiates between older and younger adults in various experimental cognitive tasks. We extend that work, assessing the relation between the EEG power spectrum and performance on the Repeatable Battery for the Assessment of Neuropsychological Status (RBANS).

METHODS

Twenty-one younger and twenty-three older adults completed the RBANS with EEG data collected at rest. Using spectral parameterization, we tested the mediating effect of the spectral slope on differences in subsequent cognitive task performance.

RESULTS

Older adults performed reliably worse on the RBANS overall, and on the Attention and Delayed Memory domains specifically. However, evidence of mediation was only found for the Coding subtest.

CONCLUSIONS

The slope of the resting EEG power spectrum mediated age-related differences in cognition, but only in a task requiring speeded processing. Mediation was not statistically significant for delayed memory, even though age-related differences were present.

Age-Related Trajectories of Brain Structure–Function Coupling in Female Roller Derby Athletes

Contact and collision sports are believed to accelerate brain aging. Postmortem studies of the human brain have implicated tau deposition in and around the perivascular space as a biomarker of an as yet poorly understood neurodegenerative process. Relatively little is known about the effects that collision sport participation has on the age-related trajectories of macroscale brain structure and function, particularly in female athletes. Diffusion MRI and resting-state functional MRI were obtained from female collision sport athletes (n = 19 roller derby (RD) players; 23–45 years old) and female control participants (n = 14; 20–49 years old) to quantify structural coupling (S_C) and decoupling (S_D). The novel and interesting finding is that RD athletes, but not controls, exhibited increasing S_C with age in two association networks: the frontoparietal network, important for cognitive control, and default-mode network, a task-negative network (permuted p = 0.0006). Age-related increases in S_C were also observed in sensorimotor networks (RD, controls) and age-related increases in S_D were observed in association networks (controls) (permuted p ≤ 0.0001). These distinct patterns suggest that competing in RD results in compressed neuronal timescales in critical networks as a function of age and encourages the broader study of female athlete brains across the lifespan.

The diversity and specificity of functional connectivity across spatial and temporal scales

Macroscopic neuroimaging modalities in humans have revealed the organization of brain-wide activity into distributed functional networks that re-organize according to behavioral demands. However, the inherent coarse-graining of macroscopic measurements conceals the diversity and specificity in responses and connectivity of many individual neurons contained in each local region. New invasive approaches in animals enable recording and manipulating neural activity at meso- and microscale resolution, with cell-type specificity and temporal precision down to milliseconds. Determining how brain-wide activity patterns emerge from interactions across spatial and temporal scales will allow us to identify the key circuit mechanisms contributing to global brain states and how the dynamic activity of these states enables adaptive behavior.

A mechanism for inter-areal coherence through communication based on connectivity and oscillatory power

Inter-areal coherence between field potentials is a widespread phenomenon in cortex. Coherence has been hypothesized to reflect phase-synchronization between oscillators and flexibly gate communication according to behavioral and cognitive demands. We reveal an alternative mechanism where coherence is not the cause but the consequence of communication and naturally emerges because spiking activity in a sending area causes post-synaptic potentials both in the same and in other areas. Consequently, coherence depends in a lawful manner on power and phase-locking in the sender and connectivity. Changes in oscillatory power explained prominent changes in fronto-parietal and LGN-V1 coherence across behavioral conditions. Optogenetic experiments and excitatory-inhibitory network simulations identified afferent synaptic inputs rather than spiking entrainment as the principal determinant of coherence. These findings suggest that unique spectral profiles of different brain areas automatically give rise to large-scale coherence patterns that follow anatomical connectivity and continuously reconfigure as a function of behavior and cognition.

Standardizing workflows in imaging transcriptomics with the abagen toolbox

Gene expression fundamentally shapes the structural and functional architecture of the human brain. Open-access transcriptomic datasets like the Allen Human Brain Atlas provide an unprecedented ability to examine these mechanisms in vivo; however, a lack of standardization across research groups has given rise to myriad processing pipelines for using these data. Here, we develop the abagen toolbox, an open-access software package for working with transcriptomic data, and use it to examine how methodological variability influences the outcomes of research using the Allen Human Brain Atlas. Applying three prototypical analyses to the outputs of 750,000 unique processing pipelines, we find that choice of pipeline has a large impact on research findings, with parameters commonly varied in the literature influencing correlations between derived gene expression and other imaging phenotypes by as much as ρ ≥ 1.0. Our results further reveal an ordering of parameter importance, with processing steps that influence gene normalization yielding the greatest impact on downstream statistical inferences and conclusions. The presented work and the development of the abagen toolbox lay the foundation for more standardized and systematic research in imaging transcriptomics, and will help to advance future understanding of the influence of gene expression in the human brain.

Sensory-motor cortices shape functional connectivity dynamics in the human brain

Large-scale biophysical circuit models provide mechanistic insights into the micro-scale and macro-scale properties of brain organization that shape complex patterns of spontaneous brain activity. We developed a spatially heterogeneous large-scale dynamical circuit model that allowed for variation in local synaptic properties across the human cortex. Here we show that parameterizing local circuit properties with both anatomical and functional gradients generates more realistic static and dynamic resting-state functional connectivity (FC). Furthermore, empirical and simulated FC dynamics demonstrates remarkably similar sharp transitions in FC patterns, suggesting the existence of multiple attractors. Time-varying regional fMRI amplitude may track multi-stability in FC dynamics. Causal manipulation of the large-scale circuit model suggests that sensory-motor regions are a driver of FC dynamics. Finally, the spatial distribution of sensory-motor drivers matches the principal gradient of gene expression that encompasses certain interneuron classes, suggesting that heterogeneity in excitation-inhibition balance might shape multi-stability in FC dynamics.

Revealing the Relevant Spatiotemporal Scale Underlying Whole-Brain Dynamics

The brain rapidly processes and adapts to new information by dynamically transitioning between whole-brain functional networks. In this whole-brain modeling study we investigate the relevance of spatiotemporal scale in whole-brain functional networks. This is achieved through estimating brain parcellations at different spatial scales (100-900 regions) and time series at different temporal scales (from milliseconds to seconds) generated by a whole-brain model fitted to fMRI data. We quantify the richness of the dynamic repertoire at each spatiotemporal scale by computing the entropy of transitions between whole-brain functional networks. The results show that the optimal relevant spatial scale is around 300 regions and a temporal scale of around 150 ms. Overall, this study provides much needed evidence for the relevant spatiotemporal scales and recommendations for analyses of brain dynamics.

When makes you unique: Temporality of the human brain fingerprint

The extraction of “fingerprints” from human brain connectivity data has become a new frontier in neuroscience. However, the time scales of human brain identifiability are still largely unexplored. We here investigate the dynamics of brain fingerprints along two complementary axes: (i) What is the optimal time scale at which brain fingerprints integrate information and (ii) when best identification happens. Using dynamic identifiability, we show that the best identification emerges at longer time scales; however, short transient “bursts of identifiability,” associated with neuronal activity, persist even when looking at shorter functional interactions. Furthermore, we report evidence that different parts of connectome fingerprints relate to different time scales, i.e., more visual-somatomotor at short temporal windows and more frontoparietal-DMN driven at increasing temporal windows. Last, different cognitive functions appear to be meta-analytically implicated in dynamic fingerprints across time scales. We hope that this investigation will advance our understanding of what makes our brains unique.

Timescales of Cognition in the Brain

We live in a world that changes on many timescales. To learn and make decisions appropriately, the human brain has evolved to integrate various types of information, such as sensory evidence and reward feedback, on multiple timescales. This is reflected in cortical hierarchies of timescales consisting of heterogeneous neuronal activities and expression of genes related to neurotransmitters critical for learning. We review the recent findings on how timescales of sensory and reward integration are affected by the temporal properties of sensory and reward signals in the environment. Despite existing evidence linking behavioral and neuronal timescales, future studies must examine how neural computations at multiple timescales are adjusted and combined to influence behavior flexibly.

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.

Atypical genomic cortical patterning in autism with poor early language outcome

Cortical regionalization develops via genomic patterning along anterior-posterior (A-P) and dorsal-ventral (D-V) gradients. Here, we find that normative A-P and D-V genomic patterning of cortical surface area (SA) and thickness (CT), present in typically developing and autistic toddlers with good early language outcome, is absent in autistic toddlers with poor early language outcome. Autistic toddlers with poor early language outcome are instead specifically characterized by a secondary and independent genomic patterning effect on CT. Genes involved in these effects can be traced back to midgestational A-P and D-V gene expression gradients and different prenatal cell types (e.g., progenitor cells and excitatory neurons), are functionally important for vocal learning and human-specific evolution, and are prominent in prenatal coexpression networks enriched for high-penetrance autism risk genes. Autism with poor early language outcome may be explained by atypical genomic cortical patterning starting in prenatal development, which may detrimentally affect later regional functional specialization and circuit formation.

Profiling Transcription Initiation in Peripheral Leukocytes Reveals Severity-Associated Cis-Regulatory Elements in Critical COVID-19

The contribution of transcription factors (TFs) and gene regulatory programs in the immune response to COVID-19 and their relationship to disease outcome is not fully understood. Analysis of genome-wide changes in transcription at both promoter-proximal and distal cis-regulatory DNA elements, collectively termed the 'active cistrome,' offers an unbiased assessment of TF activity identifying key pathways regulated in homeostasis or disease. Here, we profiled the active cistrome from peripheral leukocytes of critically ill COVID-19 patients to identify major regulatory programs and their dynamics during SARS-CoV-2 associated acute respiratory distress syndrome (ARDS). We identified TF motifs that track the severity of COVID- 19 lung injury, disease resolution, and outcome. We used unbiased clustering to reveal distinct cistrome subsets delineating the regulation of pathways, cell types, and the combinatorial activity of TFs. We found critical roles for regulatory networks driven by stimulus and lineage determining TFs, showing that STAT and E2F/MYB regulatory programs targeting myeloid cells are activated in patients with poor disease outcomes and associated with single nucleotide genetic variants implicated in COVID-19 susceptibility. Integration with single-cell RNA-seq found that STAT and E2F/MYB activation converged in specific neutrophils subset found in patients with severe disease. Collectively we demonstrate that cistrome analysis facilitates insight into disease mechanisms and provides an unbiased approach to evaluate global changes in transcription factor activity and stratify patient disease severity.

The brain and its time: intrinsic neural timescales are key for input processing

We process and integrate multiple timescales into one meaningful whole. Recent evidence suggests that the brain displays a complex multiscale temporal organization. Different regions exhibit different timescales as described by the concept of intrinsic neural timescales (INT); however, their function and neural mechanisms remains unclear. We review recent literature on INT and propose that they are key for input processing. Specifically, they are shared across different species, i.e., input sharing. This suggests a role of INT in encoding inputs through matching the inputs' stochastics with the ongoing temporal statistics of the brain's neural activity, i.e., input encoding. Following simulation and empirical data, we point out input integration versus segregation and input sampling as key temporal mechanisms of input processing. This deeply grounds the brain within its environmental and evolutionary context. It carries major implications in understanding mental features and psychiatric disorders, as well as going beyond the brain in integrating timescales into artificial intelligence.

Frontal circuit specialisations for decision making

There is widespread consensus that distributed circuits across prefrontal and anterior cingulate cortex (PFC/ACC) are critical for reward-based decision making. The circuit specialisations of these areas in primates were likely shaped by their foraging niche, in which decision making is typically sequential, attention-guided and temporally extended. Here, I argue that in humans and other primates, PFC/ACC circuits are functionally specialised in two ways. First, microcircuits found across PFC/ACC are highly recurrent in nature and have synaptic properties that support persistent activity across temporally extended cognitive tasks. These properties provide the basis of a computational account of time-varying neural activity within PFC/ACC as a decision is being made. Second, the macrocircuit connections (to other brain areas) differ between distinct PFC/ACC cytoarchitectonic subregions. This variation in macrocircuit connections explains why PFC/ACC subregions make unique contributions to reward-based decision tasks and how these contributions are shaped by attention. They predict dissociable neural representations to emerge in orbitofrontal, anterior cingulate and dorsolateral prefrontal cortex during sequential attention-guided choice, as recently confirmed in neurophysiological recordings.

Wiring of higher-order cortical areas: Spatiotemporal development of cortical hierarchy

A hierarchical development of cortical areas was suggested over a century ago, but the diversity and complexity of cortical hierarchy properties have so far prevented a formal demonstration. The aim of this review is to clarify the similarities and differences in the developmental processes underlying cortical development of primary and higher-order areas. We start by recapitulating the historical and recent advances underlying the biological principle of cortical hierarchy in adults. We then revisit the arguments for a hierarchical maturation of cortical areas, and further integrate the principles of cortical areas specification during embryonic and postnatal development. We highlight how the dramatic expansion in cortical size might have contributed to the increased number of association areas sustaining cognitive complexification in evolution. Finally, we summarize the recent observations of an alteration of cortical hierarchy in neuropsychiatric disorders and discuss their potential developmental origins.

Neuronal Sequence Models for Bayesian Online Inference

Various imaging and electrophysiological studies in a number of different species and brain regions have revealed that neuronal dynamics associated with diverse behavioral patterns and cognitive tasks take on a sequence-like structure, even when encoding stationary concepts. These neuronal sequences are characterized by robust and reproducible spatiotemporal activation patterns. This suggests that the role of neuronal sequences may be much more fundamental for brain function than is commonly believed. Furthermore, the idea that the brain is not simply a passive observer but an active predictor of its sensory input, is supported by an enormous amount of evidence in fields as diverse as human ethology and physiology, besides neuroscience. Hence, a central aspect of this review is to illustrate how neuronal sequences can be understood as critical for probabilistic predictive information processing, and what dynamical principles can be used as generators of neuronal sequences. Moreover, since different lines of evidence from neuroscience and computational modeling suggest that the brain is organized in a functional hierarchy of time scales, we will also review how models based on sequence-generating principles can be embedded in such a hierarchy, to form a generative model for recognition and prediction of sensory input. We shortly introduce the Bayesian brain hypothesis as a prominent mathematical description of how online, i.e., fast, recognition, and predictions may be computed by the brain. Finally, we briefly discuss some recent advances in machine learning, where spatiotemporally structured methods (akin to neuronal sequences) and hierarchical networks have independently been developed for a wide range of tasks. We conclude that the investigation of specific dynamical and structural principles of sequential brain activity not only helps us understand how the brain processes information and generates predictions, but also informs us about neuroscientific principles potentially useful for designing more efficient artificial neuronal networks for machine learning tasks.

Comparing spatial null models for brain maps

Technological and data sharing advances have led to a proliferation of high-resolution structural and functional maps of the brain. Modern neuroimaging research increasingly depends on identifying correspondences between the topographies of these maps; however, most standard methods for statistical inference fail to account for their spatial properties. Recently, multiple methods have been developed to generate null distributions that preserve the spatial autocorrelation of brain maps and yield more accurate statistical estimates. Here, we comprehensively assess the performance of ten published null frameworks in statistical analyses of neuroimaging data. To test the efficacy of these frameworks in situations with a known ground truth, we first apply them to a series of controlled simulations and examine the impact of data resolution and spatial autocorrelation on their family-wise error rates. Next, we use each framework with two empirical neuroimaging datasets, investigating their performance when testing (1) the correspondence between brain maps (e.g., correlating two activation maps) and (2) the spatial distribution of a feature within a partition (e.g., quantifying the specificity of an activation map within an intrinsic functional network). Finally, we investigate how differences in the implementation of these null models may impact their performance. In agreement with previous reports, we find that naive null models that do not preserve spatial autocorrelation consistently yield elevated false positive rates and unrealistically liberal statistical estimates. While spatially-constrained null models yielded more realistic, conservative estimates, even these frameworks suffer from inflated false positive rates and variable performance across analyses. Throughout our results, we observe minimal impact of parcellation and resolution on null model performance. Altogether, our findings highlight the need for continued development of statistically-rigorous methods for comparing brain maps. The present report provides a harmonised framework for benchmarking and comparing future advancements.

Inverted central auditory hierarchies for encoding local intervals and global temporal patterns

In sensory systems, representational features of increasing complexity emerge at successive stages of processing. In the mammalian auditory pathway, the clearest change from brainstem to cortex is defined by what is lost, not by what is gained, in that high-fidelity temporal coding becomes increasingly restricted to slower acoustic modulation rates.^1,2 Here, we explore the idea that sluggish temporal processing is more than just an inability for fast processing, but instead reflects an emergent specialization for encoding sound features that unfold on very slow timescales.^3,4 We performed simultaneous single unit ensemble recordings from three hierarchical stages of auditory processing in awake mice - the inferior colliculus (IC), medial geniculate body of the thalamus (MGB) and primary auditory cortex (A1). As expected, temporal coding of brief local intervals (0.001 - 0.1 s) separating consecutive noise bursts was robust in the IC and declined across MGB and A1. By contrast, slowly developing (∼1 s period) global rhythmic patterns of inter-burst interval sequences strongly modulated A1 spiking, were weakly captured by MGB neurons, and not at all by IC neurons. Shifts in stimulus regularity were not represented by changes in A1 spike rates, but rather in how the spikes were arranged in time. These findings show that low-level auditory neurons with fast timescales encode isolated sound features but not the longer gestalt, while the extended timescales in higher-level areas can facilitate sensitivity to slower contextual changes in the sensory environment.

Topographic gradients of intrinsic dynamics across neocortex

The intrinsic dynamics of neuronal populations are shaped by both microscale attributes and macroscale connectome architecture. Here we comprehensively characterize the rich temporal patterns of neural activity throughout the human brain. Applying massive temporal feature extraction to regional haemodynamic activity, we systematically estimate over 6000 statistical properties of individual brain regions' time-series across the neocortex. We identify two robust spatial gradients of intrinsic dynamics, one spanning a ventromedial-dorsolateral axis and dominated by measures of signal autocorrelation, and the other spanning a unimodal-transmodal axis and dominated by measures of dynamic range. These gradients reflect spatial patterns of gene expression, intracortical myelin and cortical thickness, as well as structural and functional network embedding. Importantly, these gradients are correlated with patterns of meta-analytic functional activation, differentiating cognitive versus affective processing and sensory versus higher-order cognitive processing. Altogether, these findings demonstrate a link between microscale and macroscale architecture, intrinsic dynamics, and cognition.

Timescales of spontaneous fMRI fluctuations relate to structural connectivity in the brain

Intrinsic timescales of activity fluctuations vary hierarchically across the brain. This variation reflects a broad gradient of functional specialization in information storage and processing, with integrative association areas displaying slower timescales that are thought to reflect longer temporal processing windows. The organization of timescales is associated with cognitive function, distinctive between individuals, and disrupted in disease, but we do not yet understand how the temporal properties of activity dynamics are shaped by the brain's underlying structural connectivity network. Using resting-state fMRI and diffusion MRI data from 100 healthy individuals from the Human Connectome Project, here we show that the timescale of resting-state fMRI dynamics increases with structural connectivity strength, matching recent results in the mouse brain. Our results hold at the level of individuals, are robust to parcellation schemes, and are conserved across a range of different timescale- related statistics. We establish a comprehensive BOLD dynamical signature of structural connectivity strength by comparing over 6,000 time series features, highlighting a range of new temporal features for characterizing BOLD dynamics, including measures of stationarity and symbolic motif frequencies. Our findings indicate a conserved property of mouse and human brain organization in which a brain region's spontaneous activity fluctuations are closely related to their surrounding structural scaffold.