Layer-dependent functional connectivity methods

Pre-frontal cortex guides dimension-reducing transformations in the occipito-ventral pathway for categorization behaviors

To interpret our surroundings, the brain uses a visual categorization process. Current theories and models suggest that this process comprises a hierarchy of different computations that transforms complex, high-dimensional inputs into lower-dimensional representations (i.e., manifolds) in support of multiple categorization behaviors. Here, we tested this hypothesis by analyzing these transformations reflected in dynamic MEG source activity while individual participants actively categorized the same stimuli according to different tasks: face expression, face gender, pedestrian gender, and vehicle type. Results reveal three transformation stages guided by the pre-frontal cortex. At stage 1 (high-dimensional, 50-120 ms), occipital sources represent both task-relevant and task-irrelevant stimulus features; task-relevant features advance into higher ventral/dorsal regions, whereas task-irrelevant features halt at the occipital-temporal junction. At stage 2 (121-150 ms), stimulus feature representations reduce to lower-dimensional manifolds, which then transform into the task-relevant features underlying categorization behavior over stage 3 (161-350 ms). Our findings shed light on how the brain's network mechanisms transform high-dimensional inputs into specific feature manifolds that support multiple categorization behaviors.

Cortical Layer-Dependent Signaling in Cognition: Three Computational Modes of the Canonical Circuit

The cerebral cortex performs computations via numerous six-layer modules. The operational dynamics of these modules were studied primarily in early sensory cortices using bottom-up computation for response selectivity as a model, which has been recently revolutionized by genetic approaches in mice. However, cognitive processes such as recall and imagery require top-down generative computation. The question of whether the layered module operates similarly in top-down generative processing as in bottom-up sensory processing has become testable by advances in the layer identification of recorded neurons in behaving monkeys. This review examines recent advances in laminar signaling in these two computations, using predictive coding computation as a common reference, and shows that each of these computations recruits distinct laminar circuits, particularly in layer 5, depending on the cognitive demands. These findings highlight many open questions, including how different interareal feedback pathways, originating from and terminating at different layers, convey distinct functional signals.

Tensor denoising of quantitative multi-parameter mapping

PURPOSE

Quantitative multi-parameter mapping (MPM) provides maps of physical quantities representing physiologically meaningful tissue characteristics, which allows to investigate microstructure-function relationships reflecting normal or pathologic processes in the brain. However, the achievable spatial resolution and stability of MPM for basic research or clinical applications is severely constrained by SNR limits of the MR acquisition process, resulting in relatively long acquisition times. To increase SNR, we denoise MPM acquisitions using principal component analysis along tensors exploiting the Marchenko-Pastur law (tMPPCA).

METHODS

tMPPCA denoising was applied to three sets of MPM raw data before the quantification of maps of proton density, magnetization transfer saturation, R1, and R2*. The regional SNR gain for high-resolution MPM was investigated as well as reproducibility gains for clinically optimized protocols with moderate and high acceleration factors at different image resolutions.

RESULTS

Substantial noise reduction in raw data was achieved, resulting in reduced noise for quantitative mapping up to sixfold without introducing bias of mean values (below 1%). Scan-rescan fluctuations were reduced up to threefold. Denoising allowed to decrease the voxel volume fourfold at the same scan time or reduce the scan time twofold at same voxel volume without loss of sensitivity.

CONCLUSIONS

tMPPCA denoising can (a) improve of fine spatial and temporal patterns, (b) considerably reduce scan time for clinical applications, or (c) increase resolution to potentially push cutting-edge MPM protocols from the upper to the lower limit of the mesoscopic scale.

RF coil design strategies for improving SNR at the ultrahigh magnetic field of 10.5 Tesla

Purpose

To develop multichannel transmit and receive arrays towards capturing the ultimate-intrinsic-SNR (uiSNR) at 10.5 Tesla (T) and to demonstrate the feasibility and potential of whole-brain, high-resolution human brain imaging at this high field strength.

Methods

A dual row 16-channel self-decoupled transmit (Tx) array was converted to a 16Tx/Rx transceiver using custom transmit/receive switches. A 64-channel receive-only (64Rx) array was built to fit into the 16Tx/Rx array. Electromagnetic modeling and experiments were employed to define safe operation limits of the resulting 16Tx/80Rx array and obtain FDA approval for human use.

Results

The 64Rx array alone captured approximately 50% of the central uiSNR at 10.5T while the identical 7T 64Rx array captured ∼76% of uiSNR at this lower field strength. The 16Tx/80Rx configuration brought the fraction of uiSNR captured at 10.5T to levels comparable to the performance of the 64Rx array at 7T. SNR data obtained at the two field strengths with these arrays displayed dependent increases over a large central region. Whole-brain high resolution T 2 * and T 1 weighted anatomical and gradient-recalled echo EPI BOLD fMRI images were obtained at 10.5T for the first time with such an advanced array, illustrating the promise of >10T fields in studying the human brain.

Conclusion

We demonstrated the ability to approach the uiSNR at 10.5T over the human brain with a novel, high channel count array, achieving large SNR gains over 7T, currently the most commonly employed ultrahigh field platform, and demonstrate high resolution and high contrast anatomical and functional imaging at 10.5T.

Principles of cortical areas and their implications for neuroimaging

Cortical organization should constrain the study of how the brain performs behavior and cognition. A fundamental concept in cortical organization is that of arealization: that the cortex is parceled into discrete areas. In part one of this report, we review how non-human animal studies have illuminated principles of cortical arealization by revealing: (1) what defines a cortical area, (2) how cortical areas are formed, (3) how cortical areas interact with one another, and (4) what "computations" or "functions" areas perform. In part two, we discuss how these principles apply to neuroimaging research. In doing so, we highlight several examples where the commonly accepted interpretation of neuroimaging observations requires assumptions that violate the principles of arealization, including nonstationary areas that move on short time scales, large-scale gradients as organizing features, and cortical areas with singular functionality that perfectly map psychological constructs. Our belief is that principles of neurobiology should strongly guide the nature of computational explanations.

Integration of temporal & spatial properties of dynamic functional connectivity based on two-directional two-dimensional principal component analysis for disease analysis

Dynamic functional connectivity, derived from resting-state functional magnetic resonance imaging (rs-fMRI), has emerged as a crucial instrument for investigating and supporting the diagnosis of neurological disorders. However, prevalent features of dynamic functional connectivity predominantly capture either temporal or spatial properties, such as mean and global efficiency, neglecting the significant information embedded in the fusion of spatial and temporal attributes. In addition, dynamic functional connectivity suffers from the problem of temporal mismatch, i.e., the functional connectivity of different subjects at the same time point cannot be matched. To address these problems, this article introduces a novel feature extraction framework grounded in two-directional two-dimensional principal component analysis. This framework is designed to extract features that integrate both spatial and temporal properties of dynamic functional connectivity. Additionally, we propose to use Fourier transform to extract temporal-invariance properties contained in dynamic functional connectivity. Experimental findings underscore the superior performance of features extracted by this framework in classification experiments compared to features capturing individual properties.

Submillimeter fMRI Acquisition Techniques for Detection of Laminar and Columnar Level Brain Activation

Since the first demonstration in the early 1990s, functional MRI (fMRI) has emerged as one of the most powerful, noninvasive neuroimaging tools to probe brain functions. Subsequently, fMRI techniques have advanced remarkably, enabling the acquisition of functional signals with a submillimeter voxel size. This innovation has opened the possibility of investigating subcortical neural activities with respect to the cortical depths or cortical columns. For this purpose, numerous previous works have endeavored to design suitable functional contrast mechanisms and dedicated imaging techniques. Depending on the choice of the functional contrast, functional signals can be detected with high sensitivity or with improved spatial specificity to the actual activation site, and the pertaining issues have been discussed in a number of earlier works. This review paper primarily aims to provide an overview of the subcortical fMRI techniques that allow the acquisition of functional signals with a submillimeter resolution. Here, the advantages and disadvantages of the imaging techniques will be described and compared. We also summarize supplementary imaging techniques that assist in the analysis of the subcortical brain activation for more accurate mapping with reduced geometric deformation. This review suggests that there is no single universally accepted method as the gold standard for subcortical fMRI. Instead, the functional contrast and the corresponding readout imaging technique should be carefully determined depending on the purpose of the study. Due to the technical limitations of current fMRI techniques, most subcortical fMRI studies have only targeted partial brain regions. As a future prospect, the spatiotemporal resolution of fMRI will be pushed to satisfy the community's need for a deeper understanding of whole-brain functions and the underlying connectivity in order to achieve the ultimate goal of a time-resolved and layer-specific spatial scale. EVIDENCE LEVEL: 1 TECHNICAL EFFICACY: Stage 1.

Does the brain behave like a (complex) network? I. Dynamics

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network structure does not entail that the brain actually works as a network. Asking whether the brain behaves as a network means asking whether network properties count. From the viewpoint of neurophysiology and, possibly, of brain physics, the most substantial issues a network structure may be instrumental in addressing relate to the influence of network properties on brain dynamics and to whether these properties ultimately explain some aspects of brain function. Here, we address the dynamical implications of complex network, examining which aspects and scales of brain activity may be understood to genuinely behave as a network. To do so, we first define the meaning of networkness, and analyse some of its implications. We then examine ways in which brain anatomy and dynamics can be endowed with a network structure and discuss possible ways in which network structure may be shown to represent a genuine organisational principle of brain activity, rather than just a convenient description of its anatomy and dynamics.

Multilayer Network Analysis across Cortical Depths in Resting-State 7T fMRI

In graph theory, "multilayer networks" represent systems involving several interconnected topological levels. A neuroscience example is the hierarchy of connections between different cortical depths or "lamina". This hierarchy is becoming non-invasively accessible in humans using ultra-high-resolution functional MRI (fMRI). Here, we applied multilayer graph theory to examine functional connectivity across different cortical depths in humans, using 7T fMRI (1-mm3 voxels; 30 participants). Blood oxygenation level dependent (BOLD) signals were derived from five depths between the white matter and pial surface. We then compared networks where the inter-regional connections were limited to a single cortical depth only ("layer-by-layer matrices") to those considering all possible connections between regions and cortical depths ("multilayer matrix"). We utilized global and local graph theory features that quantitatively characterize network attributes such as network composition, nodal centrality, path-based measures, and hub segregation. Detecting functional differences between cortical depths was improved using multilayer connectomics compared to the layer-by-layer versions. Superficial aspects of the cortex dominated information transfer and deeper aspects clustering. These differences were largest in frontotemporal and limbic brain regions. fMRI functional connectivity across different cortical depths may contain neurophysiologically relevant information. Multilayer connectomics could provide a methodological framework for studies on how information flows across this hierarchy.

Strength of predicted information content in the brain biases decision behavior

Prediction-for-perception theories suggest that the brain predicts incoming stimuli to facilitate their categorization.1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17 However, it remains unknown what the information contents of these predictions are, which hinders mechanistic explanations. This is because typical approaches cast predictions as an underconstrained contrast between two categories18,19,20,21,22,23,24-e.g., faces versus cars, which could lead to predictions of features specific to faces or cars, or features from both categories. Here, to pinpoint the information contents of predictions and thus their mechanistic processing in the brain, we identified the features that enable two different categorical perceptions of the same stimuli. We then trained multivariate classifiers to discern, from dynamic MEG brain responses, the features tied to each perception. With an auditory cueing design, we reveal where, when, and how the brain reactivates visual category features (versus the typical category contrast) before the stimulus is shown. We demonstrate that the predictions of category features have a more direct influence (bias) on subsequent decision behavior in participants than the typical category contrast. Specifically, these predictions are more precisely localized in the brain (lateralized), are more specifically driven by the auditory cues, and their reactivation strength before a stimulus presentation exerts a greater bias on how the individual participant later categorizes this stimulus. By characterizing the specific information contents that the brain predicts and then processes, our findings provide new insights into the brain's mechanisms of prediction for perception.

Next-generation MRI scanner designed for ultra-high-resolution human brain imaging at 7 Tesla

To increase granularity in human neuroimaging science, we designed and built a next-generation 7 Tesla magnetic resonance imaging scanner to reach ultra-high resolution by implementing several advances in hardware. To improve spatial encoding and increase the image signal-to-noise ratio, we developed a head-only asymmetric gradient coil (200 mT m-1, 900 T m-1s-1) with an additional third layer of windings. We integrated a 128-channel receiver system with 64- and 96-channel receiver coil arrays to boost signal in the cerebral cortex while reducing g-factor noise to enable higher accelerations. A 16-channel transmit system reduced power deposition and improved image uniformity. The scanner routinely performs functional imaging studies at 0.35-0.45 mm isotropic spatial resolution to reveal cortical layer functional activity, achieves high angular resolution in diffusion imaging and reduces acquisition time for both functional and structural imaging.

A selection and targeting framework of cortical locations for line-scanning fMRI

Depth-resolved functional magnetic resonance imaging (fMRI) is an emerging field growing in popularity given the potential of separating signals from different computational processes in cerebral cortex. Conventional acquisition schemes suffer from low spatial and temporal resolutions. Line-scanning methods allow depth-resolved fMRI by sacrificing spatial coverage to sample blood oxygenated level-dependent (BOLD) responses at ultra-high temporal and spatial resolution. For neuroscience applications, it is critical to be able to place the line accurately to (1) sample the right neural population and (2) target that neural population with tailored stimuli or tasks. To this end, we devised a multi-session framework where a target cortical location is selected based on anatomical and functional properties. The line is then positioned according to this information in a separate second session, and we tailor the experiment to focus on the target location. Anatomically, the precision of the line placement was confirmed by projecting a nominal representation of the acquired line back onto the surface. Functional estimates of neural selectivities in the line, as quantified by a visual population-receptive field model, resembled the target selectivities well for most subjects. This functional precision was quantified in detail by estimating the distance between the visual field location of the targeted vertex and the location in visual cortex (V1) that most closely resembled the line-scanning estimates; this distance was on average ~5.5 mm. Given the dimensions of the line, differences in acquisition, session, and stimulus design, this validates that line-scanning can be used to probe local neural sensitivities across sessions. In summary, we present an accurate framework for line-scanning MRI; we believe such a framework is required to harness the full potential of line-scanning and maximize its utility. Furthermore, this approach bridges canonical fMRI experiments with electrophysiological experiments, which in turn allows novel avenues for studying human physiology non-invasively.

Functional neuroimaging as a catalyst for integrated neuroscience

Functional magnetic resonance imaging (fMRI) enables non-invasive access to the awake, behaving human brain. By tracking whole-brain signals across a diverse range of cognitive and behavioural states or mapping differences associated with specific traits or clinical conditions, fMRI has advanced our understanding of brain function and its links to both normal and atypical behaviour. Despite this headway, progress in human cognitive neuroscience that uses fMRI has been relatively isolated from rapid advances in other subdomains of neuroscience, which themselves are also somewhat siloed from one another. In this Perspective, we argue that fMRI is well-placed to integrate the diverse subfields of systems, cognitive, computational and clinical neuroscience. We first summarize the strengths and weaknesses of fMRI as an imaging tool, then highlight examples of studies that have successfully used fMRI in each subdomain of neuroscience. We then provide a roadmap for the future advances that will be needed to realize this integrative vision. In this way, we hope to demonstrate how fMRI can help usher in a new era of interdisciplinary coherence in neuroscience.

Enabling brain-wide mapping of directed functional connectivity at 3T via layer-dependent fMRI with draining-vein suppression

Layer-dependent functional magnetic resonance imaging (fMRI) offers a compelling avenue for investigating directed functional connectivity (FC). To construct a comprehensive map of brain-wide directed FC, several technical criteria must be met, including sub-mm spatial resolution, adequate temporal resolution, functional sensitivity, global brain coverage, and high spatial specificity. Although gradient echo (GE)-based echo planar imaging (EPI) is commonly used for rapid fMRI acquisition, it faces significant challenges due to the draining-vein effect, particularly when utilizing blood oxygen level-dependent (BOLD) contrast. In this study, we mitigated this effect by incorporating velocity-nulling (VN) gradients into a GE-BOLD fMRI sequence, opting for a 3T magnetic field strength over 7T. We also integrated several advanced techniques, such as simultaneous multi-slice (SMS) acceleration and NORDIC denoising, to enhance temporal resolution, spatial coverage, and signal sensitivity. Collectively, the VN fMRI method exhibited notable spatial specificity, as evidenced by the identification of double-peak activation patterns within the primary motor cortex (M1) during a finger-tapping task. Additionally, the technique demonstrated BOLD sensitivity in the lateral geniculate nucleus (LGN). Furthermore, our VN fMRI technique displayed superior robustness when compared to conventional fMRI approaches across participants. Our findings of directed FC elucidate several layer-specific functional relationships between different brain regions and align closely with existing literature. Given the widespread availability of 3T scanners, this technical advancement has the potential for significant impact across multiple domains of neuroscience research.

Insights from personalized models of brain and behavior for identifying biomarkers in psychiatry

A main goal in translational neuroscience is to identify neural correlates of psychopathology ("biomarkers") that can be used to facilitate diagnosis, prognosis, and treatment. This goal has led to substantial research into how psychopathology symptoms relate to large-scale brain systems. However, these efforts have not yet resulted in practical biomarkers used in clinical practice. One reason for this underwhelming progress may be that many study designs focus on increasing sample size instead of collecting additional data within each individual. This focus limits the reliability and predictive validity of brain and behavioral measures in any one person. As biomarkers exist at the level of individuals, an increased focus on validating them within individuals is warranted. We argue that personalized models, estimated from extensive data collection within individuals, can address these concerns. We review evidence from two, thus far separate, lines of research on personalized models of (1) psychopathology symptoms and (2) fMRI measures of brain networks. We close by proposing approaches uniting personalized models across both domains to improve biomarker research.

Functional and microstructural plasticity following social and interoceptive mental training

The human brain supports social cognitive functions, including Theory of Mind, empathy, and compassion, through its intrinsic hierarchical organization. However, it remains unclear how the learning and refinement of social skills shapes brain function and structure. We studied if different types of social mental training induce changes in cortical function and microstructure, investigating 332 healthy adults (197 women, 20-55 years) with repeated multimodal neuroimaging and behavioral testing. Our neuroimaging approach examined longitudinal changes in cortical functional gradients and myelin-sensitive T1 relaxometry, two complementary measures of cortical hierarchical organization. We observed marked changes in intrinsic cortical function and microstructure, which varied as a function of social training content. In particular, cortical function and microstructure changed as a result of attention-mindfulness and socio-cognitive training in regions functionally associated with attention and interoception, including insular and parietal cortices. Conversely, socio-affective and socio-cognitive training resulted in differential microstructural changes in regions classically implicated in interoceptive and emotional processing, including insular and orbitofrontal areas, but did not result in functional reorganization. Notably, longitudinal changes in cortical function and microstructure predicted behavioral change in attention, compassion and perspective-taking. Our work demonstrates functional and microstructural plasticity after the training of social-interoceptive functions, and illustrates the bidirectional relationship between brain organisation and human social skills.

Cortical depth-dependent human fMRI of resting-state networks using EPIK

Introduction

Recent laminar-fMRI studies have substantially improved understanding of the evoked cortical responses in multiple sub-systems; in contrast, the laminar component of resting-state networks spread over the whole brain has been less studied due to technical limitations. Animal research strongly suggests that the supragranular layers of the cortex play a critical role in maintaining communication within the default mode network (DMN); however, whether this is true in this and other human cortical networks remains unclear.

Methods

Here, we used EPIK, which offers unprecedented coverage at sub-millimeter resolution, to investigate cortical broad resting-state dynamics with depth specificity in healthy volunteers.

Results

Our results suggest that human DMN connectivity is primarily supported by intermediate and superficial layers of the cortex, and furthermore, the preferred cortical depth used for communication can vary from one network to another. In addition, the laminar connectivity profile of some networks showed a tendency to change upon engagement in a motor task. In line with these connectivity changes, we observed that the amplitude of the low-frequency-fluctuations (ALFF), as well as the regional homogeneity (ReHo), exhibited a different laminar slope when subjects were either performing a task or were in a resting state (less variation among laminae, i.e., lower slope, during task performance compared to rest).

Discussion

The identification of varied laminar profiles concerning network connectivity, ALFF, and ReHo, observed across two brain states (task vs. rest) has major implications for the characterization of network-related diseases and suggests the potential diagnostic value of laminar fMRI in psychiatric disorders, e.g., to differentiate the cortical dynamics associated with disease stages linked, or not linked, to behavioral changes. The evaluation of laminar-fMRI across the brain encompasses computational challenges; nonetheless, it enables the investigation of a new dimension of the human neocortex, which may be key to understanding neurological disorders from a novel perspective.

Differential Laminar Activation Dissociates Encoding and Retrieval in the Human Medial and Lateral Entorhinal Cortex

The hierarchically organized structures of the medial temporal lobe are critically important for episodic memory function. Accumulating evidence suggests dissociable information processing pathways are maintained throughout these structures including in the medial and lateral entorhinal cortex. Cortical layers provide an additional dimension of dissociation as the primary input to the hippocampus derives from layer 2 neurons in the entorhinal cortex, whereas the deeper layers primarily receive output from the hippocampus. Here, novel high-resolution T2-prepared functional MRI methods were successfully used to mitigate susceptibility artifacts typically affecting MRI signals in this region providing uniform sensitivity across the medial and lateral entorhinal cortex. During the performance of a memory task, healthy human subjects (age 25-33 years, mean age 28.2 ± 3.3 years, 4 female) showed differential functional activation in the superficial and deep layers of the entorhinal cortex associated with task-related encoding and retrieval conditions, respectively. The methods provided here offer an approach to probe layer-specific activation in normal cognition and conditions contributing to memory impairment.SIGNIFICANCE STATEMENT This study provides new evidence for differential neuronal activation in the superficial versus deep layers of the entorhinal cortex associated with encoding and retrieval memory processes, respectively, in cognitively normal adults. The study further shows that this dissociation can be observed in both the medial and the lateral entorhinal cortex. The study was achieved by using a novel functional MRI method allowing us to measure robust functional MRI signals in both the medial and lateral entorhinal cortex that was not possible in previous studies. The methodology established here in healthy human subjects lays a solid foundation for subsequent studies investigating layer-specific and region-specific changes in the entorhinal cortex associated with memory impairment in various conditions such as Alzheimer's disease.

Relating neural oscillations to laminar fMRI connectivity in visual cortex

Laminar functional magnetic resonance imaging (fMRI) holds the potential to study connectivity at the laminar level in humans. Here we analyze simultaneously recorded electroencephalography (EEG) and high-resolution fMRI data to investigate how EEG power modulations, induced by a task with an attentional component, relate to changes in fMRI laminar connectivity between and within brain regions in visual cortex. Our results indicate that our task-induced decrease in beta power relates to an increase in deep-to-deep layer coupling between regions and to an increase in deep/middle-to-superficial layer connectivity within brain regions. The attention-related alpha power decrease predominantly relates to reduced connectivity between deep and superficial layers within brain regions, since, unlike beta power, alpha power was found to be positively correlated to connectivity. We observed no strong relation between laminar connectivity and gamma band oscillations. These results indicate that especially beta band, and to a lesser extent, alpha band oscillations relate to laminar-specific fMRI connectivity. The differential effects for alpha and beta bands indicate that they relate to different feedback-related neural processes that are differentially expressed in intra-region laminar fMRI-based connectivity.

Cerebral blood volume sensitive layer-fMRI in the human auditory cortex at 7T: Challenges and capabilities

The development of ultra high field fMRI signal readout strategies and contrasts has led to the possibility of imaging the human brain in vivo and non-invasively at increasingly higher spatial resolutions of cortical layers and columns. One emergent layer-fMRI acquisition method with increasing popularity is the cerebral blood volume sensitive sequence named vascular space occupancy (VASO). This approach has been shown to be mostly sensitive to locally-specific changes of laminar microvasculature, without unwanted biases of trans-laminar draining veins. Until now, however, VASO has not been applied in the technically challenging cortical area of the auditory cortex. Here, we describe the main challenges we encountered when developing a VASO protocol for auditory neuroscientific applications and the solutions we have adopted. With the resulting protocol, we present preliminary results of laminar responses to sounds and as a proof of concept for future investigations, we map the topographic representation of frequency preference (tonotopy) in the auditory cortex.

Towards functional spin-echo BOLD line-scanning in humans at 7T

OBJECTIVE

Neurons cluster into sub-millimeter spatial structures and neural activity occurs at millisecond resolutions; hence, ultimately, high spatial and high temporal resolutions are required for functional MRI. In this work, we implemented a spin-echo line-scanning (SELINE) sequence to use in high spatial and temporal resolution fMRI.

MATERIALS AND METHODS

A line is formed by simply rotating the spin-echo refocusing gradient to a plane perpendicular to the excited slice and by removing the phase-encoding gradient. This technique promises a combination of high spatial and temporal resolution (250 μm, 500 ms) and microvascular specificity of functional responses. We compared SELINE data to a corresponding gradient-echo version (GELINE).

RESULTS

We demonstrate that SELINE showed much-improved line selection (i.e. a sharper line profile) compared to GELINE, albeit at the cost of a significant drop in functional sensitivity.

DISCUSSION

This low functional sensitivity needs to be addressed before SELINE can be applied for neuroscientific purposes.

Layer-specific, retinotopically-diffuse modulation in human visual cortex in response to viewing emotionally expressive faces

Viewing faces that are perceived as emotionally expressive evokes enhanced neural responses in multiple brain regions, a phenomenon thought to depend critically on the amygdala. This emotion-related modulation is evident even in primary visual cortex (V1), providing a potential neural substrate by which emotionally salient stimuli can affect perception. How does emotional valence information, computed in the amygdala, reach V1? Here we use high-resolution functional MRI to investigate the layer profile and retinotopic distribution of neural activity specific to emotional facial expressions. Across three experiments, human participants viewed centrally presented face stimuli varying in emotional expression and performed a gender judgment task. We found that facial valence sensitivity was evident only in superficial cortical layers and was not restricted to the retinotopic location of the stimuli, consistent with diffuse feedback-like projections from the amygdala. Together, our results provide a feedback mechanism by which the amygdala directly modulates activity at the earliest stage of visual processing.

Long-term morphological and functional dynamics of human stem cell-derived neuronal networks on high-density micro-electrode arrays

Comprehensive electrophysiological characterizations of human induced pluripotent stem cell (hiPSC)-derived neuronal networks are essential to determine to what extent these in vitro models recapitulate the functional features of in vivo neuronal circuits. High-density micro-electrode arrays (HD-MEAs) offer non-invasive recording with the best spatial and temporal resolution possible to date. For 3 months, we tracked the morphology and activity features of developing networks derived from a transgenic hiPSC line in which neurogenesis is inducible by neurogenic transcription factor overexpression. Our morphological data revealed large-scale structural changes from homogeneously distributed neurons in the first month to the formation of neuronal clusters over time. This led to a constant shift in position of neuronal cells and clusters on HD-MEAs and corresponding changes in spatial distribution of the network activity maps. Network activity appeared as scarce action potentials (APs), evolved as local bursts with longer duration and changed to network-wide synchronized bursts with higher frequencies but shorter duration over time, resembling the emerging burst features found in the developing human brain. Instantaneous firing rate data indicated that the fraction of fast spiking neurons (150–600 Hz) increases sharply after 63 days post induction (dpi). Inhibition of glutamatergic synapses erased burst features from network activity profiles and confirmed the presence of mature excitatory neurotransmission. The application of GABAergic receptor antagonists profoundly changed the bursting profile of the network at 120 dpi. This indicated a GABAergic switch from excitatory to inhibitory neurotransmission during circuit development and maturation. Our results suggested that an emerging GABAergic system at older culture ages is involved in regulating spontaneous network bursts. In conclusion, our data showed that long-term and continuous microscopy and electrophysiology readouts are crucial for a meaningful characterization of morphological and functional maturation in stem cell-derived human networks. Most importantly, assessing the level and duration of functional maturation is key to subject these human neuronal circuits on HD-MEAs for basic and biomedical applications.

It's about time: Linking dynamical systems with human neuroimaging to understand the brain

Most human neuroscience research to date has focused on statistical approaches that describe stationary patterns of localized neural activity or blood flow. While these patterns are often interpreted in light of dynamic, information-processing concepts, the static, local, and inferential nature of the statistical approach makes it challenging to directly link neuroimaging results to plausible underlying neural mechanisms. Here, we argue that dynamical systems theory provides the crucial mechanistic framework for characterizing both the brain's time-varying quality and its partial stability in the face of perturbations, and hence, that this perspective can have a profound impact on the interpretation of human neuroimaging results and their relationship with behavior. After briefly reviewing some key terminology, we identify three key ways in which neuroimaging analyses can embrace a dynamical systems perspective: by shifting from a local to a more global perspective, by focusing on dynamics instead of static snapshots of neural activity, and by embracing modeling approaches that map neural dynamics using "forward" models. Through this approach, we envisage ample opportunities for neuroimaging researchers to enrich their understanding of the dynamic neural mechanisms that support a wide array of brain functions, both in health and in the setting of psychopathology.

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.

What is the relationship between functional network architectures inferred from electromagnetic and haemodynamic data? This study shows that superposition of electromagnetic networks at canonical frequency bands manifests as highly structured patterns of haemodynamic functional connectivity in the human brain, reflecting systematic variation in cytoarchitecture.

Toward an integrative neurovascular framework for studying brain networks

Brain functional connectivity based on the measure of blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signals has become one of the most widely used measurements in human neuroimaging. However, the nature of the functional networks revealed by BOLD fMRI can be ambiguous, as highlighted by a recent series of experiments that have suggested that typical resting-state networks can be replicated from purely vascular or physiologically driven BOLD signals. After going through a brief review of the key concepts of brain network analysis, we explore how the vascular and neuronal systems interact to give rise to the brain functional networks measured with BOLD fMRI. This leads us to emphasize a view of the vascular network not only as a confounding element in fMRI but also as a functionally relevant system that is entangled with the neuronal network. To study the vascular and neuronal underpinnings of BOLD functional connectivity, we consider a combination of methodological avenues based on multiscale and multimodal optical imaging in mice, used in combination with computational models that allow the integration of vascular information to explain functional connectivity.

The promise of layer-specific neuroimaging for testing predictive coding theories of psychosis

Predictive coding potentially provides an explanatory model for understanding the neurocognitive mechanisms of psychosis. It proposes that cognitive processes, such as perception and inference, are implemented by a hierarchical system, with the influence of each level being a function of the estimated precision of beliefs at that level. However, predictive coding models of psychosis are insufficiently constrained-any phenomenon can be explained in multiple ways by postulating different changes to precision at different levels of processing. One reason for the lack of constraint in these models is that the core processes are thought to be implemented by the function of specific cortical layers, and the technology to measure layer specific neural activity in humans has until recently been lacking. As a result, our ability to constrain the models with empirical data has been limited. In this review we provide a brief overview of predictive processing models of psychosis and then describe the potential for newly developed, layer specific neuroimaging techniques to test and thus constrain these models. We conclude by discussing the most promising avenues for this research as well as the technical and conceptual challenges which may limit its application.

Pre-processing of Sub-millimeter GE-BOLD fMRI Data for Laminar Applications

Over the past 30 years, brain function has primarily been evaluated non-invasively using functional magnetic resonance imaging (fMRI) with gradient-echo (GE) sequences to measure blood-oxygen-level-dependent (BOLD) signals. Despite the multiple advantages of GE sequences, e.g., higher signal-to-noise ratio, faster acquisitions, etc., their relatively inferior spatial localization compromises the routine use of GE-BOLD in laminar applications. Here, in an attempt to rescue the benefits of GE sequences, we evaluated the effect of existing pre-processing methods on the spatial localization of signals obtained with EPIK, a GE sequence that affords voxel volumes of 0.25 mm3 with near whole-brain coverage. The methods assessed here apply to both task and resting-state fMRI data assuming the availability of reconstructed magnitude and phase images.

Simultaneous pure T2 and varying T2'-weighted BOLD fMRI using Echo Planar Time-resolved Imaging for mapping cortical-depth dependent responses

Spin-echo (SE) BOLD fMRI has high microvascular specificity, and thus provides a more reliable means to localize neural activity compared to conventional gradient-echo BOLD fMRI. However, the most common SE BOLD acquisition method, SE-EPI, is known to suffer from T2' contrast contamination with undesirable draining vein bias. To address this, in this study, we extended a recently developed distortion/blurring-free multi-shot EPI technique, Echo-Planar Time-resolved Imaging (EPTI), to cortical-depth dependent SE-fMRI at 7T to test whether it could provide purer SE BOLD contrast with minimal T2' contamination for improved neuronal specificity. From the same acquisition, the time-resolved feature of EPTI also provides a series of asymmetric SE (ASE) images with varying T2' weightings, and enables extraction of data equivalent to conventional SE EPI with different echo train lengths (ETLs). This allows us to systematically examine how T2'-contribution affects different SE acquisition strategies using a single dataset. A low-rank spatiotemporal subspace reconstruction was implemented for the SE-EPTI acquisition, which incorporates corrections for both shot-to-shot phase variations and dynamic B0 drifts. SE-EPTI was used in a visual task fMRI experiment to demonstrate that i) the pure SE image provided by EPTI results in the highest microvascular specificity; ii) the ASE EPTI series, with a graded introduction of T2' weightings at time points farther away from the pure SE, show a gradual sensitivity increase along with increasing draining vein bias; iii) the longer ETL seen in conventional SE EPI acquisitions will induce more draining vein bias. Consistent results were observed across multiple subjects, demonstrating the robustness of the proposed technique for SE-BOLD fMRI with high specificity.

Mapping the human lateral geniculate nucleus and its cytoarchitectonic subdivisions using quantitative MRI

The human lateral geniculate nucleus (LGN) of the visual thalamus is a key subcortical processing site for visual information analysis. Due to its small size and deep location within the brain, a non-invasive characterization of the LGN and its microstructurally distinct magnocellular (M) and parvocellular (P) subdivisions in humans is challenging. Here, we investigated whether structural quantitative MRI (qMRI) methods that are sensitive to underlying microstructural tissue features enable MR-based mapping of human LGN M and P subdivisions. We employed high-resolution 7 Tesla in-vivo qMRI in N = 27 participants and ultra-high resolution 7 Tesla qMRI of a post-mortem human LGN specimen. We found that a quantitative assessment of the LGN and its subdivisions is possible based on microstructure-informed qMRI contrast alone. In both the in-vivo and post-mortem qMRI data, we identified two components of shorter and longer longitudinal relaxation time (T₁) within the LGN that coincided with the known anatomical locations of a dorsal P and a ventral M subdivision, respectively. Through ground-truth histological validation, we further showed that the microstructural MRI contrast within the LGN pertains to cyto- and myeloarchitectonic tissue differences between its subdivisions. These differences were based on cell and myelin density, but not on iron content. Our qMRI-based mapping strategy paves the way for an in-depth understanding of LGN function and microstructure in humans. It further enables investigations into the selective contributions of LGN subdivisions to human behavior in health and disease.

Simulating Local Deformations in the Human Cortex Due to Blood Flow-Induced Changes in Mechanical Tissue Properties: Impact on Functional Magnetic Resonance Imaging

Investigating human brain tissue is challenging due to the complexity and the manifold interactions between structures across different scales. Increasing evidence suggests that brain function and microstructural features including biomechanical features are related. More importantly, the relationship between tissue mechanics and its influence on brain imaging results remains poorly understood. As an important example, the study of the brain tissue response to blood flow could have important theoretical and experimental consequences for functional magnetic resonance imaging (fMRI) at high spatial resolutions. Computational simulations, using realistic mechanical models can predict and characterize the brain tissue behavior and give us insights into the consequent potential biases or limitations of in vivo, high-resolution fMRI. In this manuscript, we used a two dimensional biomechanical simulation of an exemplary human gyrus to investigate the relationship between mechanical tissue properties and the respective changes induced by focal blood flow changes. The model is based on the changes in the brain’s stiffness and volume due to the vasodilation evoked by neural activity. Modeling an exemplary gyrus from a brain atlas we assessed the influence of different potential mechanisms: (i) a local increase in tissue stiffness (at the level of a single anatomical layer), (ii) an increase in local volume, and (iii) a combination of both effects. Our simulation results showed considerable tissue displacement because of these temporary changes in mechanical properties. We found that the local volume increase causes more deformation and consequently higher displacement of the gyrus. These displacements introduced considerable artifacts in our simulated fMRI measurements. Our results underline the necessity to consider and characterize the tissue displacement which could be responsible for fMRI artifacts.

Towards new concepts for a biological neuroscience of consciousness

In the search for a sound model of consciousness, we aim at introducing new concepts: closure, compositionality, biobranes and autobranes. This is important to overcome reductionism and to bring life back into the neuroscience of consciousness. Using these definitions, we conjecture that consciousness co-arises with the non-trivial composition of biological closure in the form of biobranes and autobranes: conscious processes generate closed activity at various levels and are, in turn, themselves, supported by biobranes and autobranes. This approach leads to a non-reductionist biological and simultaneously phenomenological theory of conscious experience, giving new perspectives for a science of consciousness. Future works will implement experimental definitions and computational simulations to characterize these dynamical biobranes interacting.

Naturalistic Stimuli: A Paradigm for Multi-Scale Functional Characterization of the Human Brain

Movies, audio stories, and virtual reality are increasingly used as stimuli for functional brain imaging. Such naturalistic paradigms are in sharp contrast to the tradition of experimental reductionism in neuroscience research. Being complex, dynamic, and diverse, naturalistic stimuli set up a more ecologically relevant condition and induce highly reproducible brain responses across a wide range of spatiotemporal scales. Here, we review recent technical advances and scientific findings on imaging the brain under naturalistic stimuli. Then we elaborate on the premise of using naturalistic paradigms for multi-scale, multi-modal, and high-throughput functional characterization of the human brain. We further highlight the growing potential of using deep learning models to infer neural information processing from brain responses to naturalistic stimuli. Lastly, we advocate large-scale collaborations to combine brain imaging and recording data across experiments, subjects, and labs that use the same set of naturalistic stimuli.

LayNii: A software suite for layer-fMRI

High-resolution fMRI in the sub-millimeter regime allows researchers to resolve brain activity across cortical layers and columns non-invasively. While these high-resolution data make it possible to address novel questions of directional information flow within and across brain circuits, the corresponding data analyses are challenged by MRI artifacts, including image blurring, image distortions, low SNR, and restricted coverage. These challenges often result in insufficient spatial accuracy of conventional analysis pipelines. Here we introduce a new software suite that is specifically designed for layer-specific functional MRI: LayNii. This toolbox is a collection of command-line executable programs written in C/C++ and is distributed opensource and as pre-compiled binaries for Linux, Windows, and macOS. LayNii is designed for layer-fMRI data that suffer from SNR and coverage constraints and thus cannot be straightforwardly analyzed in alternative software packages. Some of the most popular programs of LayNii contain 'layerification' and columnarization in the native voxel space of functional data as well as many other layer-fMRI specific analysis tasks: layer-specific smoothing, model-based vein mitigation of GE-BOLD data, quality assessment of artifact dominated sub-millimeter fMRI, as well as analyses of VASO data.

Linking cortical circuit models to human cognition with laminar fMRI

Laboratory animal research has provided significant knowledge into the function of cortical circuits at the laminar level, which has yet to be fully leveraged towards insights about human brain function on a similar spatiotemporal scale. The use of functional magnetic resonance imaging (fMRI) in conjunction with neural models provides new opportunities to gain important insights from current knowledge. During the last five years, human studies have demonstrated the value of high-resolution fMRI to study laminar-specific activity in the human brain. This is mostly performed at ultra-high-field strengths (≥ 7 T) and is known as laminar fMRI. Advancements in laminar fMRI are beginning to open new possibilities for studying questions in basic cognitive neuroscience. In this paper, we first review recent methodological advances in laminar fMRI and describe recent human laminar fMRI studies. Then, we discuss how the use of laminar fMRI can help bridge the gap between cortical circuit models and human cognition.

Age-related changes in cerebrovascular health and their effects on neural function and cognition: A comprehensive review

The process of aging includes changes in cellular biology that affect local interactions between cells and their environments and eventually propagate to systemic levels. In the brain, where neurons critically depend on an efficient and dynamic supply of oxygen and glucose, age-related changes in the complex interaction between the brain parenchyma and the cerebrovasculature have effects on health and functioning that negatively impact cognition and play a role in pathology. Thus, cerebrovascular health is considered one of the main mechanisms by which a healthy lifestyle, such as habitual cardiorespiratory exercise and a healthful diet, could lead to improved cognitive outcomes with aging. This review aims at detailing how the physiology of the cerebral vascular system changes with age and how these changes lead to differential trajectories of cognitive maintenance or decline. This provides a framework for generating specific mechanistic hypotheses about the efficacy of proposed interventions and lifestyle covariates that contribute to enhanced cognitive well-being. Finally, we discuss the methodological implications of age-related changes in the cerebral vasculature for human cognitive neuroscience research and propose directions for future experiments aimed at investigating age-related changes in the relationship between physiology and cognitive mechanisms.

ULTRAHIGH FIELD and ULTRAHIGH RESOLUTION fMRI

Functional magnetic resonance imaging (fMRI) has become one of the most powerful tools for investigating the human brain. Ultrahigh magnetic field (UHF) of 7 Tesla has played a critical role in enabling higher resolution and more accurate (relative to the neuronal activity) functional maps. However, even with these gains, the fMRI approach is challenged relative to the spatial scale over which brain function is organized. Therefore, going forward, significant advances in fMRI are still needed. Such advances will predominantly come from magnetic fields significantly higher than 7 Tesla, which is the most commonly used UHF platform today, and additional technologies that will include developments in pulse sequences, image reconstruction, noise suppression, and image analysis in order to further enhance and augment the gains than can be realized by going to higher magnetic fields.

Validating layer-specific VASO across species

Cerebral blood volume (CBV) has been shown to be a robust and important physiological parameter for quantitative interpretation of functional (f)MRI, capable of delivering highly localized mapping of neural activity. Indeed, with recent advances in ultra-high-field (≥7T) MRI hardware and associated sequence libraries, it has become possible to capture non-invasive CBV weighted fMRI signals across cortical layers. One of the most widely used approaches to achieve this (in humans) is through vascular-space-occupancy (VASO) fMRI. Unfortunately, the exact contrast mechanisms of layer-dependent VASO fMRI have not been validated for human fMRI and thus interpretation of such data is confounded. Here we validate the signal source of layer-dependent SS-SI VASO fMRI using multi-modal imaging in a rat model in response to neuronal activation (somatosensory cortex) and respiratory challenge (hypercapnia). In particular VASO derived CBV measures are directly compared to concurrent measures of total haemoglobin changes from high resolution intrinsic optical imaging spectroscopy (OIS). Quantified cortical layer profiling is demonstrated to be in agreement between VASO and contrast enhanced fMRI (using monocrystalline iron oxide nanoparticles, MION). Responses show high spatial localisation to layers of cortical processing independent of confounding large draining veins which can hamper BOLD fMRI studies, (depending on slice positioning). Thus, a cross species comparison is enabled using VASO as a common measure. We find increased VASO based CBV reactivity (3.1 ± 1.2 fold increase) in humans compared to rats. Together, our findings confirm that the VASO contrast is indeed a reliable estimate of layer-specific CBV changes. This validation study increases the neuronal interpretability of human layer-dependent VASO fMRI as an appropriate method in neuroscience application studies, in which the presence of large draining intracortical and pial veins limits neuroscientific inference with BOLD fMRI.

Laminar fMRI using T2-prepared multi-echo FLASH

Functional magnetic resonance imaging (fMRI) using blood oxygenation level dependent (BOLD) contrast at a sub-millimeter scale is a promising technique to probe neural activity at the level of cortical layers. While gradient echo (GRE) BOLD sequences exhibit the highest sensitivity, their signal is confounded by unspecific extravascular (EV) and intravascular (IV) effects of large intracortical ascending veins and pial veins leading to a downstream blurring effect of local signal changes. In contrast, spin echo (SE) fMRI promises higher specificity towards signal changes near the microvascular compartment. However, the T₂-weighted signal is typically sampled with a gradient echo readout imposing additional T₂'-weighting. In this work, we used a T₂-prepared (T₂-prep) sequence with short GRE readouts to investigate its capability to acquire laminar fMRI data during a visual task in humans at 7 T. By varying the T₂-prep echo time (TE_prep) and acquiring multiple gradient echoes (TE_GRE) per excitation, we studied the specificity of the sequence and the influence of possible confounding contributions to the shape of laminar fMRI profiles. By fitting and extrapolating the multi-echo GRE data to a TE_GRE = 0 ms condition, we show for the first time laminar profiles free of T₂'-pollution, confined to gray matter. This finding is independent of TE_prep, except for the shortest one (31 ms) where hints of a remaining intravascular component can be seen. For TE_GRE > 0 ms a prominent peak at the pial surface is observed that increases with longer TE_GRE and dominates the shape of the profiles independent of the amount of T₂-weighting. Simulations show that the peak at the pial surface is a result of static EV dephasing around pial vessels in CSF visible in GM due to partial voluming. Additionally, another, weaker, static dephasing effect is observed throughout all layers of the cortex, which is particularly obvious in the data with shortest T₂-prep echo time. Our simulations show that this cannot be explained by intravascular dephasing but that it is likely caused by extravascular effects of the intracortical and pial veins. We conclude that even for TE_GRE as short as 2.3 ms, the T₂'-weighting added to the T₂-weighting is enough to dramatically affect the laminar specificity of the BOLD signal change. However, the bulk of this corruption stems from CSF partial volume effects which can in principle be addressed by increasing the spatial resolution of the acquisition.

The nonhuman primate neuroimaging and neuroanatomy project

Multi-modal neuroimaging projects such as the Human Connectome Project (HCP) and UK Biobank are advancing our understanding of human brain architecture, function, connectivity, and their variability across individuals using high-quality non-invasive data from many subjects. Such efforts depend upon the accuracy of non-invasive brain imaging measures. However, 'ground truth' validation of connectivity using invasive tracers is not feasible in humans. Studies using nonhuman primates (NHPs) enable comparisons between invasive and non-invasive measures, including exploration of how "functional connectivity" from fMRI and "tractographic connectivity" from diffusion MRI compare with long-distance connections measured using tract tracing. Our NonHuman Primate Neuroimaging & Neuroanatomy Project (NHP_NNP) is an international effort (6 laboratories in 5 countries) to: (i) acquire and analyze high-quality multi-modal brain imaging data of macaque and marmoset monkeys using protocols and methods adapted from the HCP; (ii) acquire quantitative invasive tract-tracing data for cortical and subcortical projections to cortical areas; and (iii) map the distributions of different brain cell types with immunocytochemical stains to better define brain areal boundaries. We are acquiring high-resolution structural, functional, and diffusion MRI data together with behavioral measures from over 100 individual macaques and marmosets in order to generate non-invasive measures of brain architecture such as myelin and cortical thickness maps, as well as functional and diffusion tractography-based connectomes. We are using classical and next-generation anatomical tracers to generate quantitative connectivity maps based on brain-wide counting of labeled cortical and subcortical neurons, providing ground truth measures of connectivity. Advanced statistical modeling techniques address the consistency of both kinds of data across individuals, allowing comparison of tracer-based and non-invasive MRI-based connectivity measures. We aim to develop improved cortical and subcortical areal atlases by combining histological and imaging methods. Finally, we are collecting genetic and sociality-associated behavioral data in all animals in an effort to understand how genetic variation shapes the connectome and behavior.

Multimodal EEG-fMRI: advancing insight into large-scale human brain dynamics

Advances in the acquisition and analysis of functional magnetic resonance imaging (fMRI) data are revealing increasingly rich spatiotemporal structure across the human brain. Nonetheless, uncertainty surrounding the origins of fMRI hemodynamic signals, and in the link between large-scale fMRI patterns and ongoing functional states, presently limits the neurobiological conclusions one can draw from fMRI alone. Electroencephalography (EEG) provides complementary information about neural electrical activity and state change, and simultaneously acquiring EEG together with fMRI presents unique opportunities for studying large-scale brain activity and gaining more information from fMRI itself. Here, we discuss recent progress in the use of concurrent EEG-fMRI to enrich the investigation of neural and physiological states and clarify the origins of fMRI hemodynamic signals. Throughout, we outline perspectives on future directions and open challenges.

Cortical hierarchy, dual counterstream architecture and the importance of top-down generative networks

Hierarchy is a major organizational principle of the cortex and underscores modern computational theories of cortical function. The local microcircuit amplifies long-distance inter-areal input, which show distance-dependent changes in their laminar profiles. Statistical modeling of these changes in laminar profiles demonstrates that inputs from multiple hierarchical levels to their target areas show remarkable consistency, allowing the construction of a cortical hierarchy based on a principle of hierarchical distance. The statistical modeling that is applied to structure can also be applied to laminar differences in the oscillatory coherence between areas thereby determining a functional hierarchy of the cortex. Close examination of the anatomy of inter-areal connectivity reveals a dual counterstream architecture with well-defined distance-dependent feedback and feedforward pathways in both the supra- and infragranular layers, suggesting a multiplicity of feedback pathways with well-defined functional properties. These findings are consistent with feedback connections providing a generative network involved in a wide range of cognitive functions. A dynamical model constrained by connectivity data sheds insight into the experimentally observed signatures of frequency-dependent Granger causality for feedforward versus feedback signaling. Concerted experiments capitalizing on recent technical advances and combining tract-tracing, high-resolution fMRI, optogenetics and mathematical modeling hold the promise of a much improved understanding of lamina-constrained mechanisms of neural computation and cognition. However, because inter-areal interactions involve cortical layers that have been the target of important evolutionary changes in the primate lineage, these investigations will need to include human and non-human primate comparisons.

Comparing hand movement rate dependence of cerebral blood volume and BOLD responses at 7T

Functional magnetic resonance imaging (fMRI) based on the Blood Oxygenation Level Dependent (BOLD) contrast takes advantage of the coupling between neuronal activity and the hemodynamics to allow a non-invasive localisation of the neuronal activity. In general, fMRI experiments assume a linear relationship between neuronal activation and the observed hemodynamics. However, the relationship between BOLD responses, neuronal activity, and behaviour are often nonlinear. In addition, the nonlinearity between BOLD responses and behaviour may be related to neuronal process rather than a neurovascular uncoupling. Further, part of the nonlinearity may be driven by vascular nonlinearity effects in particular from large vessel contributions. fMRI based on cerebral blood volume (CBV), promises a higher microvascular specificity, potentially without vascular nonlinearity effects and reduced contamination of the large draining vessels compared to BOLD. In this study, we aimed to investigate differences in BOLD and VASO-CBV signal changes during a hand movement task over a broad range of movement rates. We used a double readout 3D-EPI sequence at 7T to simultaneously measure VASO-CBV and BOLD responses in the sensorimotor cortex. The measured BOLD and VASO-CBV responses increased very similarly in a nonlinear fashion, plateauing for movement rates larger than 1 Hz. Our findings show a tight relationship between BOLD and VASO-CBV responses, indicating that the overall interplay of CBV and BOLD responses are similar for the assessed range of movement rates. These results suggest that the observed nonlinearity of neuronal origin is already present in VASO-CBV measurements, and consequently shows relatively unchanged BOLD responses.

Forging a path to mesoscopic imaging success with ultra-high field functional magnetic resonance imaging

Functional magnetic resonance imaging (fMRI) studies with ultra-high field (UHF, 7+ Tesla) technology enable the acquisition of high-resolution images. In this work, we discuss recent achievements in UHF fMRI at the mesoscopic scale, on the order of cortical columns and layers, and examine approaches to addressing common challenges. As researchers push to smaller and smaller voxel sizes, acquisition and analysis decisions have greater potential to degrade spatial accuracy, and UHF fMRI data must be carefully interpreted. We consider the impact of acquisition decisions on the spatial specificity of the MR signal with a representative dataset with 0.8 mm isotropic resolution. We illustrate the trade-offs in contrast with noise ratio and spatial specificity of different acquisition techniques and show that acquisition blurring can increase the effective voxel size by as much as 50% in some dimensions. We further describe how different sources of degradations to spatial resolution in functional data may be characterized. Finally, we emphasize that progress in UHF fMRI depends not only on scientific discovery and technical advancement, but also on informal discussions and documentation of challenges researchers face and overcome in pursuit of their goals. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

New acquisition techniques and their prospects for the achievable resolution of fMRI

This work reviews recent advances in technologies for functional magnetic resonance imaging (fMRI) of the human brain and highlights the push for higher functional specificity based on increased spatial resolution and specific MR contrasts to reveal previously undetectable functional properties of small-scale cortical structures. We discuss how the combination of MR hardware, advanced acquisition techniques and various MR contrast mechanisms have enabled recent progress in functional neuroimaging. However, these advanced fMRI practices have only been applied to a handful of neuroscience questions to date, with the majority of the neuroscience community still using conventional imaging techniques. We thus discuss upcoming challenges and possibilities for fMRI technology development in human neuroscience. We hope that readers interested in functional brain imaging acquire an understanding of current and novel developments and potential future applications, even if they don't have a background in MR physics or engineering. We summarize the capabilities of standard fMRI acquisition schemes with pointers to relevant literature and comprehensive reviews and introduce more recent developments.