Evaluating increases in sensitivity from NORDIC for diverse fMRI acquisition strategies

Terra Incognita - Contributions of the Olivo-Cerebellar System to Autism Spectrum Disorder

The inferior olivary nuclei (ION), a key component of the olivo-cerebellar system, remain understudied in autism spectrum disorder (ASD) research despite evident cerebellar involvement. This perspective piece aims to elucidate the critical role of the ION in cerebellar microcircuitry and its potential implications in ASD pathophysiology. We review the olivo- cerebellar system and the structural and functional alterations of the ION in autism, highlighting findings from neuroanatomical, neuroimaging, and behavioral studies. ION disruptions, although underexplored, may have a major role in the symptomatology of ASD, particularly higher-order cognitive abilities, rapid stimuli processing, and motor coordination. Specifically, we highlight how anomalies in olivary neuron morphology and olivo-cerebellar connectivity patterns may underlie deficits in temporal processing and motor learning observed in ASD. Furthermore, we discuss the challenges in brainstem imaging and recent advancements in ultra-high field (UHF) 7 T MRI technology, as standard neuroimaging techniques. As these neuroimaging techniques continue to evolve, further investigation of the functional territories of the ION holds promise for providing essential understandings into ASD processes and may lead to pioneering therapeutic options targeting this crucial brainstem area.

The human brainstem's red nucleus was upgraded to support goal-directed action

The red nucleus, a large brainstem structure, coordinates limb movement for locomotion in quadrupedal animals. In humans, its pattern of anatomical connectivity differs from that of quadrupeds, suggesting a different purpose. Here, we apply our most advanced resting-state functional connectivity based precision functional mapping in highly sampled individuals (n = 5), resting-state functional connectivity in large group-averaged datasets (combined n ~ 45,000), and task based analysis of reward, motor, and action related contrasts from group-averaged datasets (n > 1000) and meta-analyses (n > 14,000 studies) to precisely examine red nucleus function. Notably, red nucleus functional connectivity with motor-effector networks (somatomotor hand, foot, and mouth) is minimal. Instead, connectivity is strongest to the action-mode and salience networks, which are important for action/cognitive control and reward/motivated behavior. Consistent with this, the red nucleus responds to motor planning more than to actual movement, while also responding to rewards. Our results suggest the human red nucleus implements goal-directed behavior by integrating behavioral valence and action plans instead of serving a pure motor-effector function.

Mapping curvature domains in human V4 using CBV-sensitive layer-fMRI at 3T

Introduction

A full understanding of how we see our world remains a fundamental research question in vision neuroscience. While topographic profiling has allowed us to identify different visual areas, the exact functional characteristics and organization of areas up in the visual hierarchy (beyond V1 & V2) is still debated. It is hypothesized that visual area V4 represents a vital intermediate stage of processing spatial and curvature information preceding object recognition. Advancements in magnetic resonance imaging hardware and acquisition techniques (e.g., non-BOLD functional MRI) now permits the capture of cortical layer-specific functional properties and organization of the human brain (including the visual system) at high precision.

Methods

Here, we use functional cerebral blood volume measures to study the modularity in how responses to contours (curvature) are organized within area V4 of the human brain. To achieve this at 3 Tesla (a clinically relevant field strength) we utilize optimized high-resolution 3D-Echo Planar Imaging (EPI) Vascular Space Occupancy (VASO) measurements.

Results

Data here provide the first evidence of curvature domains in human V4 that are consistent with previous findings from non-human primates. We show that VASO and BOLD tSNR maps for functional imaging align with high field equivalents, with robust time series of changes to visual stimuli measured across the visual cortex. V4 curvature preference maps for VASO show strong modular organization compared to BOLD imaging contrast. It is noted that BOLD has a much lower sensitivity (due to known venous vasculature weightings) and specificity to stimulus contrast. We show evidence that curvature domains persist across the cortical depth. The work advances our understanding of the role of mid-level area V4 in human processing of curvature and shape features.

Impact

Knowledge of how the functional architecture and hierarchical integration of local contours (curvature) contribute to formation of shapes can inform computational models of object recognition. Techniques described here allow for quantification of individual differences in functional architecture of mid-level visual areas to help drive a better understanding of how changes in functional brain organization relate to difference in visual perception.

NORDIC denoising on VASO data

The use of submillimeter resolution functional magnetic resonance imaging (fMRI) is increasing in popularity due to the prospect of studying human brain activation non-invasively at the scale of cortical layers and columns. This method, known as laminar fMRI, is inherently signal-to-noise ratio (SNR)-limited, especially at lower field strengths, with the dominant noise source being of thermal origin. Furthermore, laminar fMRI is challenged with signal displacements due to draining vein effects in conventional gradient-echo blood oxygen level-dependent (BOLD) imaging contrasts. fMRI contrasts such as cerebral blood volume (CBV)-sensitive vascular space occupancy (VASO) sequences have the potential to mitigate draining vein effects. However, VASO comes along with another reduction in detection sensitivity. NOise Reduction with DIstribution Corrected (NORDIC) PCA (principal component analysis) is a denoising technique specifically aimed at suppressing thermal noise, which has proven useful for increasing the SNR of high-resolution functional data. While NORDIC has been examined for BOLD acquisitions, its application to VASO data has been limited, which was the focus of the present study. We present a preliminary analysis to evaluate NORDIC's capability to suppress thermal noise while preserving the VASO signal across a wide parameter space at 3T. For the data presented here, with a proper set of parameters, NORDIC reduced thermal noise with minimal bias on the underlying signal and preserved spatial resolution. Denoising performance was found to vary with different implementation strategies and parameter choices, for which we provide recommendations. We conclude that when applied properly, NORDIC has the potential to overcome the sensitivity limitations of laminar-specific VASO fMRI. Since very few groups currently have 3T VASO data, by sharing our analysis and code, we can compile and compare the effects of NORDIC across a broader range of acquisition parameters and study designs. Such a communal effort will help develop robust recommendations that will increase the utility of laminar fMRI at lower field strengths.

Principles of intensive human neuroimaging

The rise of large, publicly shared functional magnetic resonance imaging (fMRI) data sets in human neuroscience has focused on acquiring either a few hours of data on many individuals ('wide' fMRI) or many hours of data on a few individuals ('deep' fMRI). In this opinion article, we highlight an emerging approach within deep fMRI, which we refer to as 'intensive' fMRI: one that strives for extensive sampling of cognitive phenomena to support computational modeling and detailed investigation of brain function at the single voxel level. We discuss the fundamental principles, trade-offs, and practical considerations of intensive fMRI. We also emphasize that intensive fMRI does not simply mean collecting more data: it requires careful design of experiments to enable a rich hypothesis space, optimizing data quality, and strategically curating public resources to maximize community impact.

Precision Functional Mapping to Advance Developmental Psychiatry Research

Many psychiatric conditions have their roots in early development. Individual differences in prenatal brain function (which is influenced by a combination of genetic risk and the prenatal environment) likely interact with individual differences in postnatal experience, resulting in substantial variation in brain functional organization and development in infancy. Neuroimaging has been a powerful tool for understanding typical and atypical brain function and holds promise for uncovering the neurodevelopmental basis of psychiatric illness; however, its clinical utility has been relatively limited thus far. A substantial challenge in this endeavor is the traditional approach of averaging brain data across groups despite individuals varying in their brain organization, which likely obscures important clinically relevant individual variation. Precision functional mapping (PFM) is a neuroimaging technique that allows the capture of individual-specific and highly reliable functional brain properties. Here, we discuss how PFM, through its focus on individuals, has provided novel insights for understanding brain organization across the life span and its promise in elucidating the neural basis of psychiatric disorders. We first summarize the extant literature on PFM in normative populations, followed by its limited utilization in studying psychiatric conditions in adults. We conclude by discussing the potential for infant PFM in advancing developmental precision psychiatry applications, given that many psychiatric disorders start during early infancy and are associated with changes in individual-specific functional neuroanatomy. By exploring the intersection of PFM, development, and psychiatric research, this article underscores the importance of individualized approaches in unraveling the complexities of brain function and improving clinical outcomes across development.

Multi-echo Acquisition and Thermal Denoising Advances Precision Functional Imaging

The characterization of individual functional brain organization with Precision Functional Mapping has provided important insights in recent years in adults. However, little is known about the ontogeny of inter-individual differences in brain functional organization during human development. Precise characterization of systems organization during periods of high plasticity is likely to be essential for discoveries promoting lifelong health. Obtaining precision fMRI data during development has unique challenges that highlight the importance of establishing new methods to improve data acquisition, processing, and analysis. Here, we investigate two methods that can facilitate attaining this goal: multi-echo (ME) data acquisition and thermal noise removal with Noise Reduction with Distribution Corrected (NORDIC) principal component analysis. We applied these methods to precision fMRI data from adults, children, and newborn infants. In adults, both ME acquisitions and NORDIC increased temporal signal to noise ratio (tSNR) as well as the split-half reliability of functional connectivity matrices, with the combination helping more than either technique alone. The benefits of NORDIC denoising replicated in both our developmental samples. ME acquisitions revealed longer and more variable T2* relaxation times across the brain in infants relative to older children and adults, leading to major differences in the echo weighting for optimally combining ME data. This result suggests ME acquisitions may be a promising tool for optimizing developmental fMRI, albeit application in infants needs further investigation. The present work showcases methodological advances that improve Precision Functional Mapping in adults and developmental populations and, at the same time, highlights the need for further improvements in infant specific fMRI.

NOise Reduction with DIstribution Corrected (NORDIC) principal component analysis improves brain activity detection across rodent and human functional MRI contexts

NOise Reduction with DIstribution Corrected (NORDIC) principal component analysis (PCA) has been shown to selectively suppress thermal noise and improve the temporal signal-to-noise ratio (tSNR) in human functional magnetic resonance imaging (fMRI). However, the feasibility to improve data quality for rodent fMRI using NORDIC PCA remains uncertain. NORDIC PCA may also be particularly beneficial for improving topological brain mapping, as conventional mapping requires precise spatiotemporal signals from large datasets (ideally ~1 hour acquisition) for individual representations. In this study, we evaluated the effects of NORDIC PCA compared with "Standard" processing in various rodent fMRI contexts that range from task-evoked optogenetic fMRI to resting-state fMRI. We also evaluated the effects of NORDIC PCA on human resting-state and retinotopic mapping fMRI via population receptive field (pRF) modeling. In rodent optogenetic fMRI, apart from doubling the tSNR, NORDIC PCA resulted in a larger number of activated voxels and a significant decrease in the variance of evoked brain responses without altering brain morphology. In rodent resting-state fMRI, we found that NORDIC PCA induced a nearly threefold increase in tSNR and preserved task-free relative cerebrovascular reactivity (rCVR) across cortical depth. NORDIC PCA further improved the detection of TGN020-induced aquaporin-4 inhibition on rCVR compared with Standard processing without NORDIC PCA. NORDIC PCA also increased the tSNR for both human resting-state and pRF fMRI, and for the latter also increased activation cluster sizes while retaining retinotopic organization. This suggests that NORDIC PCA preserves the spatiotemporal precision of fMRI signals needed for pRF analysis, and effectively captures small activity changes with high sensitivity. Taken together, these results broadly demonstrate the value of NORDIC PCA for the enhanced detection of neural dynamics across various rodent and human fMRI contexts. This can in turn play an important role in improving fMRI image quality and sensitivity for translational and preclinical neuroimaging research.

Precision fMRI and cluster-failure in the individual brain

Advances in neuroimaging acquisition protocols and denoising techniques, along with increasing magnetic field strengths, have dramatically improved the temporal signal-to-noise ratio (tSNR) in functional magnetic resonance imaging (fMRI). This permits spatial resolution with submillimeter voxel sizes and ultrahigh temporal resolution and opens a route toward performing precision fMRI in the brains of individuals. Yet ultrahigh spatial and temporal resolution comes at a cost: it reduces tSNR and, therefore, the sensitivity to the blood oxygen level-dependent (BOLD) effect and other functional contrasts across the brain. Here we investigate the potential of various smoothing filters to improve BOLD sensitivity while preserving the spatial accuracy of activated clusters in single-subject analysis. We introduce adaptive-weight smoothing with optimized metrics (AWSOM), which addresses this challenge extremely well. AWSOM employs a local inference approach that is as sensitive as cluster-corrected inference of data smoothed with large Gaussian kernels, but it preserves spatial details across multiple tSNR levels. This is essential for examining whole-brain fMRI data because tSNR varies across the entire brain, depending on the distance of a brain region from the receiver coil, the type of setup, acquisition protocol, preprocessing, and resolution. We found that cluster correction in single subjects results in inflated family-wise error and false positive rates. AWSOM effectively suppresses false positives while remaining sensitive even to small clusters of activated voxels. Furthermore, it preserves signal integrity, that is, the relative activation strength of significant voxels, making it a valuable asset for a wide range of fMRI applications. Here we demonstrate these features and make AWSOM freely available to the research community for download.

Leveraging Multi-Echo EPI to Enhance BOLD Sensitivity in Task-based Olfactory fMRI

Functional magnetic resonance imaging (fMRI) using blood-oxygenation-level-dependent (BOLD) contrast relies on gradient echo echo-planar imaging (GE-EPI) to quantify dynamic susceptibility changes associated with the hemodynamic response to neural activity. However, acquiring BOLD fMRI in human olfactory regions is particularly challenging due to their proximity to the sinuses where large susceptibility gradients induce magnetic field distortions. BOLD fMRI of the human olfactory system is further complicated by respiratory artifacts that are highly correlated with event onsets in olfactory tasks. Multi-echo EPI (ME-EPI) acquires gradient echo data at multiple echo times (TEs) during a single acquisition and can leverage signal evolution over the multiple echo times to enhance BOLD sensitivity and reduce artifactual signal contributions. In the current study, we developed a ME-EPI acquisition protocol for olfactory task-based fMRI and demonstrated significant improvement in BOLD signal sensitivity over conventional single-echo EPI (1E-EPI). The observed improvement arose from both an increase in BOLD signal changes through a T 2 * -weighted echo combination and a reduction in non-BOLD artifacts through the application of the Multi-Echo Independent Components Analysis (ME-ICA) denoising method. This study represents one of the first direct comparisons between 1E-EPI and ME-EPI in high-susceptibility regions and provides compelling evidence in favor of using ME-EPI for future task-based fMRI studies.

Evaluating the effect of denoising submillimeter auditory fMRI data with NORDIC

Functional magnetic resonance imaging (fMRI) has emerged as an essential tool for exploring human brain function. Submillimeter fMRI, in particular, has emerged as a tool to study mesoscopic computations. The inherently low signal-to-noise ratio (SNR) at submillimeter resolutions warrants the use of denoising approaches tailored at reducing thermal noise-the dominant contributing noise component in high-resolution fMRI. NOise Reduction with DIstribution Corrected Principal Component Analysis (NORDIC PCA) is one of such approaches, and has been benchmarked against other approaches in several applications. Here, we investigate the effects that two versions of NORDIC denoising have on auditory submillimeter data. While investigating auditory functional responses poses unique challenges, we anticipated NORDIC to have a positive impact on the data on the basis of previous applications. Our results show that NORDIC denoising improves the detection sensitivity and the reliability of estimates in submillimeter auditory fMRI data. These effects can be explained by the reduction of the noise-induced signal variability. However, we did observe a reduction in the average response amplitude (percent signal change) within regions of interest, which may suggest that a portion of the signal of interest, which could not be distinguished from general i.i.d. noise, was also removed. We conclude that, while evaluating the effects of the signal reduction induced by NORDIC may be necessary for each application, using NORDIC in high-resolution auditory fMRI studies may be advantageous because of the large reduction in variability of the estimated responses.

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.

Towards personalized precision functional mapping in infancy

The precise network topology of functional brain systems is highly specific to individuals and undergoes dramatic changes during critical periods of development. Large amounts of high-quality resting state data are required to investigate these individual differences, but are difficult to obtain in early infancy. Using the template matching method, we generated a set of infant network templates to use as priors for individualized functional resting-state network mapping in two independent neonatal datasets with extended acquisition of resting-state functional MRI (fMRI) data. We show that template matching detects all major adult resting-state networks in individual infants and that the topology of these resting-state network maps is individual-specific. Interestingly, there was no plateau in within-subject network map similarity with up to 25 minutes of resting-state data, suggesting that the amount and/or quality of infant data required to achieve stable or high-precision network maps is higher than adults. These findings are a critical step towards personalized precision functional brain mapping in infants, which opens new avenues for clinical applicability of resting-state fMRI and potential for robust prediction of how early functional connectivity patterns relate to subsequent behavioral phenotypes and health outcomes.

Evaluating the effect of denoising submillimeter auditory fMRI data with NORDIC

Functional magnetic resonance imaging (fMRI) has emerged as an essential tool for exploring human brain function. Submillimeter fMRI, in particular, has emerged as a tool to study mesoscopic computations. The inherently low signal-to-noise ratio (SNR) at submillimeter resolutions warrants the use of denoising approaches tailored at reducing thermal noise - the dominant contributing noise component in high resolution fMRI. NORDIC PCA is one of such approaches, and has been benchmarked against other approaches in several applications. Here, we investigate the effects that two versions of NORDIC denoising have on auditory submillimeter data. As investigating auditory functional responses poses unique challenges, we anticipated that the benefit of this technique would be especially pronounced. Our results show that NORDIC denoising improves the detection sensitivity and the reliability of estimates in submillimeter auditory fMRI data. These effects can be explained by the reduction of the noise-induced signal variability. However, we also observed a reduction in the average response amplitude (percent signal), which may suggest that a small amount of signal was also removed. We conclude that, while evaluating the effects of the signal reduction induced by NORDIC may be necessary for each application, using NORDIC in high resolution auditory fMRI studies may be advantageous because of the large reduction in variability of the estimated responses.

The brainstem's red nucleus was evolutionarily upgraded to support goal-directed action

The red nucleus is a large brainstem structure that coordinates limb movement for locomotion in quadrupedal animals (Basile et al., 2021). The humans red nucleus has a different pattern of anatomical connectivity compared to quadrupeds, suggesting a unique purpose (Hatschek, 1907). Previously the function of the human red nucleus remained unclear at least partly due to methodological limitations with brainstem functional neuroimaging (Sclocco et al., 2018). Here, we used our most advanced resting-state functional connectivity (RSFC) based precision functional mapping (PFM) in highly sampled individuals (n = 5) and large group-averaged datasets (combined N ~ 45,000), to precisely examine red nucleus functional connectivity. Notably, red nucleus functional connectivity to motor-effector networks (somatomotor hand, foot, and mouth) was minimal. Instead, red nucleus functional connectivity along the central sulcus was specific to regions of the recently discovered somato-cognitive action network (SCAN; (Gordon et al., 2023)). Outside of primary motor cortex, red nucleus connectivity was strongest to the cingulo-opercular (CON) and salience networks, involved in action/cognitive control (Dosenbach et al., 2007; Newbold et al., 2021) and reward/motivated behavior (Seeley, 2019), respectively. Functional connectivity to these two networks was organized into discrete dorsal-medial and ventral-lateral zones. Red nucleus functional connectivity to the thalamus recapitulated known structural connectivity of the dento-rubral thalamic tract (DRTT) and could prove clinically useful in functionally targeting the ventral intermediate (VIM) nucleus. In total, our results indicate that far from being a 'motor' structure, the red nucleus is better understood as a brainstem nucleus for implementing goal-directed behavior, integrating behavioral valence and action plans.

Reliability of resting-state functional connectivity in the human spinal cord: Assessing the impact of distinct noise sources

The investigation of spontaneous fluctuations of the blood-oxygen-level-dependent (BOLD) signal has recently been extended from the brain to the spinal cord, where it has stimulated interest from a clinical perspective. A number of resting-state functional magnetic resonance imaging (fMRI) studies have demonstrated robust functional connectivity between the time series of BOLD fluctuations in bilateral dorsal horns and between those in bilateral ventral horns, in line with the functional neuroanatomy of the spinal cord. A necessary step prior to extension to clinical studies is assessing the reliability of such resting-state signals, which we aimed to do here in a group of 45 healthy young adults at the clinically prevalent field strength of 3T. When investigating connectivity in the entire cervical spinal cord, we observed fair to good reliability for dorsal-dorsal and ventral-ventral connectivity, whereas reliability was poor for within- and between-hemicord dorsal-ventral connectivity. Considering how prone spinal cord fMRI is to noise, we extensively investigated the impact of distinct noise sources and made two crucial observations: removal of physiological noise led to a reduction in functional connectivity strength and reliability - due to the removal of stable and participant-specific noise patterns - whereas removal of thermal noise considerably increased the detectability of functional connectivity without a clear influence on reliability. Finally, we also assessed connectivity within spinal cord segments and observed that while the pattern of connectivity was similar to that of whole cervical cord, reliability at the level of single segments was consistently poor. Taken together, our results demonstrate the presence of reliable resting-state functional connectivity in the human spinal cord even after thoroughly accounting for physiological and thermal noise, but at the same time urge caution if focal changes in connectivity (e.g. due to segmental lesions) are to be studied, especially in a longitudinal manner.

MP-PCA denoising of fMRI time-series data can lead to artificial activation "spreading"

MP-PCA denoising has become the method of choice for denoising MRI data since it provides an objective threshold to separate the signal components from unwanted thermal noise components. In rodents, thermal noise in the coils is an important source of noise that can reduce the accuracy of activation mapping in fMRI. Further confounding this problem, vendor data often contains zero-filling and other post-processing steps that may violate MP-PCA assumptions. Here, we develop an approach to denoise vendor data and assess activation "spreading" caused by MP-PCA denoising in rodent task-based fMRI data. Data was obtained from N = 3 mice using conventional multislice and ultrafast acquisitions (1 s and 50 ms temporal resolution, respectively), during visual stimulation. MP-PCA denoising produced SNR gains of 64% and 39% and Fourier Spectral Amplitude (FSA) increases in BOLD maps of 9% and 7% for multislice and ultrafast data, respectively, when using a small [2 2] denoising window. Larger windows provided higher SNR and FSA gains with increased spatial extent of activation that may or may not represent real activation. Simulations showed that MP-PCA denoising can incur activation "spreading" with increased false positive rate and smoother functional maps due to local "bleeding" of principal components, and that the optimal denoising window for improved specificity of functional mapping, based on Dice score calculations, depends on the data's tSNR and functional CNR. This "spreading" effect applies also to another recently proposed low-rank denoising method (NORDIC), although to a lesser degree. Our results bode well for enhancing spatial and/or temporal resolution in future fMRI work, while taking into account the sensitivity/specificity trade-offs of low-rank denoising methods.