Effects of phase regression on high-resolution functional MRI of the primary visual cortex

Reliability of task-based fMRI in the dorsal horn of the human spinal cord

The application of functional magnetic resonance imaging (fMRI) to the human spinal cord is still a relatively small field of research and faces many challenges. Here we aimed to probe the limitations of task-based spinal fMRI at 3T by investigating the reliability of spinal cord blood oxygen level dependent (BOLD) responses to repeated nociceptive stimulation across two consecutive days in 40 healthy volunteers. We assessed the test-retest reliability of subjective ratings, autonomic responses, and spinal cord BOLD responses to short heat pain stimuli (1s duration) using the intraclass correlation coefficient (ICC). At the group level, we observed robust autonomic responses as well as spatially specific spinal cord BOLD responses at the expected location, but no spatial overlap in BOLD response patterns across days. While autonomic indicators of pain processing showed good-to-excellent reliability, both β-estimates and z-scores of task-related BOLD responses showed poor reliability across days in the target region (gray matter of the ipsilateral dorsal horn). When taking into account the sensitivity of gradient-echo echo planar imaging (GE-EPI) to draining vein signals by including the venous plexus in the analysis, we observed BOLD responses with fair reliability across days. Taken together, these results demonstrate that heat pain stimuli as short as one second are able to evoke a robust and spatially specific BOLD response, which is however strongly variable within participants across time, resulting in low reliability in the dorsal horn gray matter. Further improvements in data acquisition and analysis techniques are thus necessary before event-related spinal cord fMRI as used here can be reliably employed in longitudinal designs or clinical settings.

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.

Predicting the macrovascular contribution to resting-state fMRI functional connectivity at 3 Tesla: A model-informed approach

Macrovascular biases have been a long-standing challenge for fMRI, limiting its ability to detect spatially specific neural activity. Recent experimental studies, including our own (Huck et al., 2023; Zhong et al., 2023), found substantial resting-state macrovascular BOLD fMRI contributions from large veins and arteries, extending into the perivascular tissue at 3 T and 7 T. The objective of this study is to demonstrate the feasibility of predicting, using a biophysical model, the experimental resting-state BOLD fluctuation amplitude (RSFA) and associated functional connectivity (FC) values at 3 Tesla. We investigated the feasibility of both 2D and 3D infinite-cylinder models as well as macrovascular anatomical networks (mVANs) derived from angiograms. Our results demonstrate that: 1) with the availability of mVANs, it is feasible to model macrovascular BOLD FC using both the mVAN-based model and 3D infinite-cylinder models, though the former performed better; 2) biophysical modelling can accurately predict the BOLD pairwise correlation near to large veins (with R 2 ranging from 0.53 to 0.93 across different subjects), but not near to large arteries; 3) compared with FC, biophysical modelling provided less accurate predictions for RSFA; 4) modelling of perivascular BOLD connectivity was feasible at close distances from veins (with R 2 ranging from 0.08 to 0.57), but not arteries, with performance deteriorating with increasing distance. While our current study demonstrates the feasibility of simulating macrovascular BOLD in the resting state, our methodology may also apply to understanding task-based BOLD. Furthermore, these results suggest the possibility of correcting for macrovascular bias in resting-state fMRI and other types of fMRI using biophysical modelling based on vascular anatomy.

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.

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.

Effects of 12-week escitalopram treatment on resting-state functional connectivity of large-scale brain networks in major depressive disorder

In this study, the effects of antidepressants on large-scale brain networks and the neural basis of individual differences in response were explored. A total of 41 patients with major depressive disorder (MDD) and 42 matched healthy controls (HCs) were scanned by resting-state functional magnetic resonance imaging separately at baseline and after a 12-week follow-up. The patients with MDD received escitalopram for 12 weeks. After treatment, patients were classified into those with MDD in remission [MDDr, endpoint 17-item Hamilton Depression Rating Scale (HAMD) total score ≤7] and those in nonremission (MDDnr). The human Brainnetome Atlas was used to define large-scale networks and compute within- and between-network resting-state functional connectivity (rsFC). Results showed the decreased subcortical network (SCN)-ventral attention network (VAN) connectivity at baseline increased in patients with MDD after 12-week treatment, and it was comparable with that of HCs. This change was only observed in patients with MDDr. However, the decreased within-network rsFC in SCN and default mode network (DMN) persisted in all patients with MDD, including those with MDDr and MDDnr, after treatment. The strength of SCN-VAN connectivity at baseline was significantly negatively correlated with the reduction rate of HAMD score in all patients with MDD. Thus, SCN-VAN connectivity may be an antidepressant target associated with depressive state changes and a predictor of treatment response to serotonin reuptake inhibitors. The within-network rsFC in SCN and DMN may reflect a trait-like abnormality in MDD. These findings provide further insights into the mechanism of antidepressants and their individual differences in response. The trial name is "Appropriate technology study of MDD diagnosis and treatment based on objective indicators and measurement" (URL: http://www.chictr.org.cn/showproj.aspx?proj=21377; registration number: ChiCTR-OOC-17012566).

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.

Laminar functional magnetic resonance imaging in vision research

Magnetic resonance imaging (MRI) scanners at ultra-high magnetic fields have become available to use in humans, thus enabling researchers to investigate the human brain in detail. By increasing the spatial resolution, ultra-high field MR allows both structural and functional characterization of cortical layers. Techniques that can differentiate cortical layers, such as histological studies and electrode-based measurements have made critical contributions to the understanding of brain function, but these techniques are invasive and thus mainly available in animal models. There are likely to be differences in the organization of circuits between humans and even our closest evolutionary neighbors. Thus research on the human brain is essential. Ultra-high field MRI can observe differences between cortical layers, but is non-invasive and can be used in humans. Extensive previous literature has shown that neuronal connections between brain areas that transmit feedback and feedforward information terminate in different layers of the cortex. Layer-specific functional MRI (fMRI) allows the identification of layer-specific hemodynamic responses, distinguishing feedback and feedforward pathways. This capability has been particularly important for understanding visual processing, as it has allowed researchers to test hypotheses concerning feedback and feedforward information in visual cortical areas. In this review, we provide a general overview of successful ultra-high field MRI applications in vision research as examples of future research.

Correcting for Superficial Bias in 7T Gradient Echo fMRI

The arrival of submillimeter ultra high-field fMRI makes it possible to compare activation profiles across cortical layers. However, the blood oxygenation level dependent (BOLD) signal measured by gradient echo (GE) fMRI is biased toward superficial layers of the cortex, which is a serious confound for laminar analysis. Several univariate and multivariate analysis methods have been proposed to correct this bias. We compare these methods using computational simulations of 7T fMRI data from regions of interest (ROI) during a visual attention paradigm. We also tested the methods on a pilot dataset of human 7T fMRI data. The simulations show that two methods–the ratio of ROI means across conditions and a novel application of Deming regression–offer the most robust correction for superficial bias. Deming regression has the additional advantage that it does not require that the conditions differ in their mean activation over voxels within an ROI. When applied to the pilot dataset, we observed strikingly different layer profiles when different attention metrics were used, but were unable to discern any differences in laminar attention across layers when Deming regression or ROI ratio was applied. Our simulations demonstrates that accurate correction of superficial bias is crucial to avoid drawing erroneous conclusions from laminar analyses of GE fMRI data, and this is affirmed by the results from our pilot 7T fMRI data.

Receiver phase alignment using fitted SVD derived sensitivities from routine prescans

Magnetic resonance imaging radio frequency arrays are composed of multiple receive coils that have their signals combined to form an image. Combination requires an estimate of the radio frequency coil sensitivities to align signal phases and prevent destructive interference. At lower fields this can be accomplished using a uniform physical reference coil. However, at higher fields, uniform volume coils are lacking and, when available, suffer from regions of low receive sensitivity that result in poor sensitivity estimation and combination. Several approaches exist that do not require a physical reference coil but require manual intervention, specific prescans, or must be completed post-acquisition. This makes these methods impractical for large multi-volume datasets such as those collected for novel types of functional MRI or quantitative susceptibility mapping, where magnitude and phase are important. This pilot study proposes a fitted SVD method which utilizes existing combination methods to create a phase sensitive combination method targeted at large multi-volume datasets. This method uses any multi-image prescan to calculate the relative receive sensitivities using voxel-wise singular value decomposition. These relative sensitivities are fitted to the solid harmonics using an iterative least squares fitting algorithm. Fits of the relative sensitivities are used to align the phases of the receive coils and improve combination in subsequent acquisitions during the imaging session. This method is compared against existing approaches in the human brain at 7 Tesla by examining the combined data for the presence of singularities and changes in phase signal-to-noise ratio. Two additional applications of the method are also explored, using the fitted SVD method in an asymmetrical coil and in a case with subject motion. The fitted SVD method produces singularity-free images and recovers between 95-100% of the phase signal-to-noise ratio depending on the prescan data resolution. Using solid harmonic fitting to interpolate singular value decomposition derived receive sensitivities from existing prescans allows the fitted SVD method to be used on all acquisitions within a session without increasing exam duration. Our fitted SVD method is able to combine imaging datasets accurately without supervision during online reconstruction.