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.

Mode-based morphometry: A multiscale approach to mapping human neuroanatomy.

Voxel-based morphometry (VBM) and surface-based morphometry (SBM) are two widely used neuroimaging techniques for investigating brain anatomy. These techniques rely on statistical inferences at individual points (voxels or vertices), clusters of points, or a priori regions-of-interest. They are powerful tools for describing brain anatomy, but offer little insights into the generative processes that shape a particular set of findings. Moreover, they are restricted to a single spatial resolution scale, precluding the opportunity to distinguish anatomical variations that are expressed across multiple scales. Drawing on concepts from classical physics, here we develop an approach, called mode-based morphometry (MBM), that can describe any empirical map of anatomical variations in terms of the fundamental, resonant modes-eigenmodes-of brain anatomy, each tied to a specific spatial scale. Hence, MBM naturally yields a multiscale characterization of the empirical map, affording new opportunities for investigating the spatial frequency content of neuroanatomical variability. Using simulated and empirical data, we show that the validity and reliability of MBM are either comparable or superior to classical vertex-based SBM for capturing differences in cortical thickness maps between two experimental groups. Our approach thus offers a robust, accurate, and informative method for characterizing empirical maps of neuroanatomical variability that can be directly linked to a generative physical process.

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.

Cortical depth-dependent modeling of visual hemodynamic responses.

A physiologically based three-dimensional (3D) hemodynamic model is developed to predict the experimentally observed blood oxygen level dependent (BOLD) responses versus the cortical depth induced by visual stimuli. Prior 2D approximations are relaxed in order to analyze 3D blood flow dynamics as a function of cortical depth. Comparison of the predictions with experimental data for evoked stimuli demonstrates that the full 3D model performs at least as well as previous approaches while remaining parsimonious. In particular, the 3D model requires significantly fewer assumptions and model parameters than previous models such that there is no longer need to define depth-specific parameter values for spatial spreading, peak amplitude, and hemodynamic velocity.

Two Spatially Distinct Posterior Alpha Sources Fulfill Different Functional Roles in Attention.

Directing attention helps to extract relevant information and suppress distracters. Alpha brain oscillations (8-12 Hz) are crucial for this process, with power decreases facilitating processing of important information and power increases inhibiting brain regions processing irrelevant information. Evidence for this phenomenon arises from visual attention studies (Worden et al., 2000); however, the effect also exists in other modalities, including the somatosensory system (Haegens et al., 2011) and intersensory attention tasks (Foxe and Snyder, 2011). We investigated in human participants (10 females, 10 males) the role of alpha oscillations in focused (0/100%) versus divided (40/60%) attention, both across modalities (visual/somatosensory; Experiment 1) and within the same modality (visual domain: across hemifields; Experiment 2) while recording EEG over 128 scalp electrodes. In Experiment 1, participants divided their attention between visual and somatosensory modality to determine the temporal/spatial frequency of a target stimulus (vibrotactile stimulus/Gabor grating). In Experiment 2, participants divided attention between two visual hemifields to identify the orientation of a Gabor grating. In both experiments, prestimulus alpha power in visual areas decreased linearly with increasing attention to visual stimuli. In contrast, prestimulus alpha power in parietal areas was lower when attention was divided between modalities/hemifields compared with focused attention. These results suggest there are two alpha sources, one of which reflects the "visual spotlight of attention" and the other reflects attentional effort. To our knowledge, this is the first study to show that attention recruits two spatially distinct alpha sources in occipital and parietal brain regions, acting simultaneously but serving different functions in attention.SIGNIFICANCE STATEMENT Attention to one spatial location/sensory modality leads to power changes of alpha oscillations (∼10 Hz) with decreased power over regions processing relevant information and power increases to actively inhibit areas processing "to-be-ignored" information. Here, we used detailed source modeling to investigate EEG data recorded during separate unimodal (visual) and multimodal (visual and somatosensory) attention tasks. Participants either focused their attention on one modality/spatial location or directed it to both. We show for the first time two distinct alpha sources are active simultaneously but play different roles. A sensory (visual) alpha source was linearly modulated by attention representing the "visual spotlight of attention." By contrast, a parietal alpha source was modulated by attentional effort, showing lowest alpha power when attention was divided.

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.

Addressing challenges of high spatial resolution UHF fMRI for group analysis of higher-order cognitive tasks: An inter-sensory task directing attention between visual and somatosensory domains.

Functional MRI at ultra-high field (UHF, ≥7 T) provides significant increases in BOLD contrast-to-noise ratio (CNR) compared with conventional field strength (3 T), and has been exploited for reduced field-of-view, high spatial resolution mapping of primary sensory areas. Applying these high spatial resolution methods to investigate whole brain functional responses to higher-order cognitive tasks leads to a number of challenges, in particular how to perform robust group-level statistical analyses. This study addresses these challenges using an inter-sensory cognitive task which modulates top-down attention at graded levels between the visual and somatosensory domains. At the individual level, highly focal functional activation to the task and task difficulty (modulated by attention levels) were detectable due to the high CNR at UHF. However, to assess group level effects, both anatomical and functional variability must be considered during analysis. We demonstrate the importance of surface over volume normalisation and the requirement of no spatial smoothing when assessing highly focal activity. Using novel group analysis on anatomically parcellated brain regions, we show that in higher cognitive areas (parietal and dorsal-lateral-prefrontal cortex) fMRI responses to graded attention levels were modulated quadratically, whilst in visual cortex and VIP, responses were modulated linearly. These group fMRI responses were not seen clearly using conventional second-level GLM analyses, illustrating the limitations of a conventional approach when investigating such focal responses in higher cognitive regions which are more anatomically variable. The approaches demonstrated here complement other advanced analysis methods such as multivariate pattern analysis, allowing UHF to be fully exploited in cognitive neuroscience.

The spatiotemporal hemodynamic response function for depth-dependent functional imaging of human cortex.

The gray matter of human cortex is characterized by depth-dependent differences in neuronal activity and connections (Shipp, 2007) as well as in the associated vasculature (Duvernoy et al., 1981). The resolution limit of functional magnetic resonance imaging (fMRI) measurements is now below a millimeter, promising the non-invasive measurement of these properties in awake and behaving humans (Muckli et al., 2015; Olman et al., 2012; Ress et al., 2007). To advance this endeavor, we present a detailed spatiotemporal hemodynamic response function (HRF) reconstructed through the use of high-resolution, submillimeter fMRI. We decomposed the HRF into directions tangential and perpendicular to the cortical surface and found that key spatial properties of the HRF change significantly with depth from the cortical surface. Notably, we found that the spatial spread of the HRF increases linearly from 4.8mm at the gray/white matter boundary to 6.6mm near the cortical surface. Using a hemodynamic model, we posit that this effect can be explained by the depth profile of the cortical vasculature, and as such, must be taken into account to properly estimate the underlying neuronal responses at different cortical depths.

A multiscale characterization of cortical shape asymmetries in early psychosis.

Psychosis has often been linked to abnormal cortical asymmetry, but prior results have been inconsistent. Here, we applied a novel spectral shape analysis to characterize cortical shape asymmetries in patients with early psychosis across different spatial scales. We used the Human Connectome Project for Early Psychosis dataset (aged 16-35), comprising 56 healthy controls (37 males, 19 females) and 112 patients with early psychosis (68 males, 44 females). We quantified shape variations of each hemisphere over different spatial frequencies and applied a general linear model to compare differences between healthy controls and patients with early psychosis. We further used canonical correlation analysis to examine associations between shape asymmetries and clinical symptoms. Cortical shape asymmetries, spanning wavelengths from about 22 to 75 mm, were significantly different between healthy controls and patients with early psychosis (Cohen's d = 0.28-0.51), with patients showing greater asymmetry in cortical shape than controls. A single canonical mode linked the asymmetry measures to symptoms (canonical correlation analysis r = 0.45), such that higher cortical asymmetry was correlated with more severe excitement symptoms and less severe emotional distress. Significant group differences in the asymmetries of traditional morphological measures of cortical thickness, surface area, and gyrification, at either global or regional levels, were not identified. Cortical shape asymmetries are more sensitive than other morphological asymmetries in capturing abnormalities in patients with early psychosis. These abnormalities are expressed at coarse spatial scales and are correlated with specific symptom domains.

Hemodynamic traveling waves in human visual cortex.

Functional MRI (fMRI) experiments rely on precise characterization of the blood oxygen level dependent (BOLD) signal. As the spatial resolution of fMRI reaches the sub-millimeter range, the need for quantitative modelling of spatiotemporal properties of this hemodynamic signal has become pressing. Here, we find that a detailed physiologically-based model of spatiotemporal BOLD responses predicts traveling waves with velocities and spatial ranges in empirically observable ranges. Two measurable parameters, related to physiology, characterize these waves: wave velocity and damping rate. To test these predictions, high-resolution fMRI data are acquired from subjects viewing discrete visual stimuli. Predictions and experiment show strong agreement, in particular confirming BOLD waves propagating for at least 5-10 mm across the cortical surface at speeds of 2-12 mm s-1. These observations enable fundamentally new approaches to fMRI analysis, crucial for fMRI data acquired at high spatial resolution.

Mode-based morphometry: A multiscale approach to mapping human neuroanatomy.

Voxel-based morphometry (VBM) and surface-based morphometry (SBM) are two widely used neuroimaging techniques for investigating brain anatomy. These techniques rely on statistical inferences at individual points (voxels or vertices), clusters of points, or a priori regions-of-interest. They are powerful tools for describing brain anatomy, but offer little insights into the generative processes that shape a particular set of findings. Moreover, they are restricted to a single spatial resolution scale, precluding the opportunity to distinguish anatomical variations that are expressed across multiple scales. Drawing on concepts from classical physics, here we develop an approach, called mode-based morphometry (MBM), that can describe any empirical map of anatomical variations in terms of the fundamental, resonant modes--eigenmodes--of brain anatomy, each tied to a specific spatial scale. Hence, MBM naturally yields a multiscale characterization of the empirical map, affording new opportunities for investigating the spatial frequency content of neuroanatomical variability. Using simulated and empirical data, we show that the validity and reliability of MBM are either comparable or superior to classical vertex-based SBM for capturing differences in cortical thickness maps between two experimental groups. Our approach thus offers a robust, accurate, and informative method for characterizing empirical maps of neuroanatomical variability that can be directly linked to a generative physical process.

Disruptions of Hierarchical Cortical Organization in Early Psychosis and Schizophrenia.

BACKGROUND: The cerebral cortex is organized hierarchically along an axis that spans unimodal sensorimotor to transmodal association areas. This hierarchy is often characterized using low-dimensional embeddings, termed gradients, of interregional functional coupling estimates measured with resting-state functional magnetic resonance imaging. Such analyses may offer insights into the pathophysiology of schizophrenia, which has been frequently linked to dysfunctional interactions between association and sensorimotor areas. METHODS: To examine disruptions of hierarchical cortical function across distinct stages of psychosis, we applied diffusion map embedding to 2 independent functional magnetic resonance imaging datasets: one comprising 114 patients with early psychosis and 48 control participants, and the other comprising 50 patients with established schizophrenia and 121 control participants. Then, we analyzed the primary sensorimotor-to-association and secondary visual-to-sensorimotor gradients of each participant in both datasets. RESULTS: There were no significant differences in regional gradient scores between patients with early psychosis and control participants. Patients with established schizophrenia showed significant differences in the secondary, but not primary, gradient compared with control participants. Gradient differences in schizophrenia were characterized by lower within-network dispersion in the dorsal attention (false discovery rate [FDR]-corrected p [pFDR] < .001), visual (pFDR = .003), frontoparietal (pFDR = .018), and limbic (pFDR = .020) networks and lower between-network dispersion between the visual network and other networks (pFDR < .001). CONCLUSIONS: These findings indicate that differences in cortical hierarchical function occur along the secondary visual-to-sensorimotor axis rather than the primary sensorimotor-to-association axis as previously thought. The absence of differences in early psychosis suggests that visual-sensorimotor abnormalities may emerge as the illness progresses.

The music of the hemispheres: Cortical eigenmodes as a physical basis for large-scale brain activity and connectivity patterns.

Neuroscience has had access to high-resolution recordings of large-scale cortical activity and structure for decades, but still lacks a generally adopted basis to analyze and interrelate results from different individuals and experiments. Here it is argued that the natural oscillatory modes of the cortex-cortical eigenmodes-provide a physically preferred framework for systematic comparisons across experimental conditions and imaging modalities. In this framework, eigenmodes are analogous to notes of a musical instrument, while commonly used statistical patterns parallel frequently played chords. This intuitive perspective avoids problems that often arise in neuroimaging analyses, and connects to underlying mechanisms of brain activity. We envisage this approach will lead to novel insights into whole-brain function, both in existing and prospective datasets, and facilitate a unification of empirical findings across presently disparate analysis paradigms and measurement modalities.

The individuality of shape asymmetries of the human cerebral cortex.

Asymmetries of the cerebral cortex are found across diverse phyla and are particularly pronounced in humans, with important implications for brain function and disease. However, many prior studies have confounded asymmetries due to size with those due to shape. Here, we introduce a novel approach to characterize asymmetries of the whole cortical shape, independent of size, across different spatial frequencies using magnetic resonance imaging data in three independent datasets. We find that cortical shape asymmetry is highly individualized and robust, akin to a cortical fingerprint, and identifies individuals more accurately than size-based descriptors, such as cortical thickness and surface area, or measures of inter-regional functional coupling of brain activity. Individual identifiability is optimal at coarse spatial scales (~37 mm wavelength), and shape asymmetries show scale-specific associations with sex and cognition, but not handedness. While unihemispheric cortical shape shows significant heritability at coarse scales (~65 mm wavelength), shape asymmetries are determined primarily by subject-specific environmental effects. Thus, coarse-scale shape asymmetries are highly personalized, sexually dimorphic, linked to individual differences in cognition, and are primarily driven by stochastic environmental influences.

Determination of Dynamic Brain Connectivity via Spectral Analysis.

Spectral analysis based on neural field theory is used to analyze dynamic connectivity via methods based on the physical eigenmodes that are the building blocks of brain dynamics. These approaches integrate over space instead of averaging over time and thereby greatly reduce or remove the temporal averaging effects, windowing artifacts, and noise at fine spatial scales that have bedeviled the analysis of dynamical functional connectivity (FC). The dependences of FC on dynamics at various timescales, and on windowing, are clarified and the results are demonstrated on simple test cases, demonstrating how modes provide directly interpretable insights that can be related to brain structure and function. It is shown that FC is dynamic even when the brain structure and effective connectivity are fixed, and that the observed patterns of FC are dominated by relatively few eigenmodes. Common artifacts introduced by statistical analyses that do not incorporate the physical nature of the brain are discussed and it is shown that these are avoided by spectral analysis using eigenmodes. Unlike most published artificially discretized "resting state networks" and other statistically-derived patterns, eigenmodes overlap, with every mode extending across the whole brain and every region participating in every mode-just like the vibrations that give rise to notes of a musical instrument. Despite this, modes are independent and do not interact in the linear limit. It is argued that for many purposes the intrinsic limitations of covariance-based FC instead favor the alternative of tracking eigenmode coefficients vs. time, which provide a compact representation that is directly related to biophysical brain dynamics.

Dynamical consequences of regional heterogeneity in the brain's transcriptional landscape.

Brain regions vary in their molecular and cellular composition, but how this heterogeneity shapes neuronal dynamics is unclear. Here, we investigate the dynamical consequences of regional heterogeneity using a biophysical model of whole-brain functional magnetic resonance imaging (MRI) dynamics in humans. We show that models in which transcriptional variations in excitatory and inhibitory receptor (E:I) gene expression constrain regional heterogeneity more accurately reproduce the spatiotemporal structure of empirical functional connectivity estimates than do models constrained by global gene expression profiles or MRI-derived estimates of myeloarchitecture. We further show that regional transcriptional heterogeneity is essential for yielding both ignition-like dynamics, which are thought to support conscious processing, and a wide variance of regional-activity time scales, which supports a broad dynamical range. We thus identify a key role for E:I heterogeneity in generating complex neuronal dynamics and demonstrate the viability of using transcriptomic data to constrain models of large-scale brain function.