Simultaneous cortical, subcortical, and brainstem mapping of sensory activation.

Nonpainful tactile sensory stimuli are processed in the cortex, subcortex, and brainstem. Recent functional magnetic resonance imaging studies have highlighted the value of whole-brain, systems-level investigation for examining sensory processing. However, whole-brain functional magnetic resonance imaging studies are uncommon, in part due to challenges with signal to noise when studying the brainstem. Furthermore, differentiation of small sensory brainstem structures such as the cuneate and gracile nuclei necessitates high-resolution imaging. To address this gap in systems-level sensory investigation, we employed a whole-brain, multi-echo functional magnetic resonance imaging acquisition at 3T with multi-echo independent component analysis denoising and brainstem-specific modeling to enable detection of activation across the entire sensory system. In healthy participants, we examined patterns of activity in response to nonpainful brushing of the right hand, left hand, and right foot (n = 10 per location), and found the expected lateralization, with distinct cortical and subcortical responses for upper and lower limb stimulation. At the brainstem level, we differentiated the adjacent cuneate and gracile nuclei, corresponding to hand and foot stimulation respectively. Our findings demonstrate that simultaneous cortical, subcortical, and brainstem mapping at 3T could be a key tool to understand the sensory system in both healthy individuals and clinical cohorts with sensory deficits.

Hemodynamic timing in resting-state and breathing-task BOLD fMRI.

The blood flow response to a vasoactive stimulus demonstrates regional heterogeneity across both the healthy brain and in cerebrovascular pathology. The timing of a regional hemodynamic response is emerging as an important biomarker of cerebrovascular dysfunction, as well as a confound within fMRI analyses. Previous research demonstrated that hemodynamic timing is more robustly characterized when a larger systemic vascular response is evoked by a breathing challenge, compared to when only spontaneous fluctuations in vascular physiology are present (i.e., in resting-state data). However, it is not clear whether hemodynamic delays in these two conditions are physiologically interchangeable, and how methodological signal-to-noise factors may limit their agreement. To address this, we generated whole-brain maps of hemodynamic delays in nine healthy adults. We assessed the agreement of voxel-wise gray matter (GM) hemodynamic delays between two conditions: resting-state and breath-holding. We found that delay values demonstrated poor agreement when considering all GM voxels, but increasingly greater agreement when limiting analyses to voxels showing strong correlation with the GM mean time-series. Voxels showing the strongest agreement with the GM mean time-series were primarily located near large venous vessels, however these voxels explain some, but not all, of the observed agreement in timing. Increasing the degree of spatial smoothing of the fMRI data enhanced the correlation between individual voxel time-series and the GM mean time-series. These results suggest that signal-to-noise factors may be limiting the accuracy of voxel-wise timing estimates and hence their agreement between the two data segments. In conclusion, caution must be taken when using voxel-wise delay estimates from resting-state and breathing-task data interchangeably, and additional work is needed to evaluate their relative sensitivity and specificity to aspects of vascular physiology and pathology.

Comparing end-tidal CO2, respiration volume per time (RVT), and average gray matter signal for mapping cerebrovascular reactivity amplitude and delay with breath-hold task BOLD fMRI.

Cerebrovascular reactivity (CVR), defined as the cerebral blood flow response to a vasoactive stimulus, is an imaging biomarker with demonstrated utility in a range of diseases and in typical development and aging processes. A robust and widely implemented method to map CVR involves using a breath-hold task during a BOLD fMRI scan. Recording end-tidal CO2 (PETCO2) changes during the breath-hold task is recommended to be used as a reference signal for modeling CVR amplitude in standard units (%BOLD/mmHg) and CVR delay in seconds. However, obtaining reliable PETCO2 recordings requires equipment and task compliance that may not be achievable in all settings. To address this challenge, we investigated two alternative reference signals to map CVR amplitude and delay in a lagged general linear model (lagged-GLM) framework: respiration volume per time (RVT) and average gray matter BOLD response (GM-BOLD). In 8 healthy adults with multiple scan sessions, we compare spatial agreement of CVR maps from RVT and GM-BOLD to those generated with PETCO2. We define a threshold to determine whether a PETCO2 recording has "sufficient" quality for CVR mapping and perform these comparisons in 16 datasets with sufficient PETCO2 and 6 datasets with insufficient PETCO2. When PETCO2 quality is sufficient, both RVT and GM-BOLD produce CVR amplitude maps that are nearly identical to those from PETCO2 (after accounting for differences in scale), with the caveat they are not in standard units to facilitate between-group comparisons. CVR delays are comparable to PETCO2 with an RVT regressor but may be underestimated with the average GM-BOLD regressor. Importantly, when PETCO2 quality is insufficient, RVT and GM-BOLD CVR recover reasonable CVR amplitude and delay maps, provided the participant attempted the breath-hold task. Therefore, our framework offers a solution for achieving high quality CVR maps in both retrospective and prospective studies where sufficient PETCO2 recordings are not available and especially in populations where obtaining reliable measurements is a known challenge (e.g., children). Our results have the potential to improve the accessibility of CVR mapping and to increase the prevalence of this promising metric of vascular health.

Spatial distribution of hand-grasp motor task activity in spinal cord functional magnetic resonance imaging.

Upper extremity motor paradigms during spinal cord functional magnetic resonance imaging (fMRI) can provide insight into the functional organization of the cord. Hand-grasping is an important daily function with clinical significance, but previous studies of similar squeezing movements have not reported consistent areas of activity and are limited by sample size and simplistic analysis methods. Here, we study spinal cord fMRI activation using a unimanual isometric hand-grasping task that is calibrated to participant maximum voluntary contraction (MVC). Two task modeling methods were considered: (1) a task regressor derived from an idealized block design (Ideal) and (2) a task regressor based on the recorded force trace normalized to individual MVC (%MVC). Across these two methods, group motor activity was highly lateralized to the hemicord ipsilateral to the side of the task. Activation spanned C5-C8 and was primarily localized to the C7 spinal cord segment. Specific differences in spatial distribution are also observed, such as an increase in C8 and dorsal cord activity when using the %MVC regressor. Furthermore, we explored the impact of data quantity and spatial smoothing on sensitivity to hand-grasp motor task activation. This analysis shows a large increase in number of active voxels associated with the number of fMRI runs, sample size, and spatial smoothing, demonstrating the impact of experimental design choices on motor activation.

A practical modification to a resting state fMRI protocol for improved characterization of cerebrovascular function.

Cerebrovascular reactivity (CVR), defined here as the Blood Oxygenation Level Dependent (BOLD) response to a CO2 pressure change, is a useful metric of cerebrovascular function. Both the amplitude and the timing (hemodynamic lag) of the CVR response can bring insight into the nature of a cerebrovascular pathology and aid in understanding noise confounds when using functional Magnetic Resonance Imaging (fMRI) to study neural activity. This research assessed a practical modification to a typical resting-state fMRI protocol, to improve the characterization of cerebrovascular function. In 9 healthy subjects, we modelled CVR and lag in three resting-state data segments, and in data segments which added a 2-3 minute breathing task to the start of a resting-state segment. Two different breathing tasks were used to induce fluctuations in arterial CO2 pressure: a breath-hold task to induce hypercapnia (CO2 increase) and a cued deep breathing task to induce hypocapnia (CO2 decrease). Our analysis produced voxel-wise estimates of the amplitude (CVR) and timing (lag) of the BOLD-fMRI response to CO2 by systematically shifting the CO2 regressor in time to optimize the model fit. This optimization inherently increases gray matter CVR values and fit statistics. The inclusion of a simple breathing task, compared to a resting-state scan only, increases the number of voxels in the brain that have a significant relationship between CO2 and BOLD-fMRI signals, and improves our confidence in the plausibility of voxel-wise CVR and hemodynamic lag estimates. We demonstrate the clinical utility and feasibility of this protocol in an incidental finding of Moyamoya disease, and explore the possibilities and challenges of using this protocol in younger populations. This hybrid protocol has direct applications for CVR mapping in both research and clinical settings and wider applications for fMRI denoising and interpretation.

ICA-based denoising strategies in breath-hold induced cerebrovascular reactivity mapping with multi echo BOLD fMRI.

Performing a BOLD functional MRI (fMRI) acquisition during breath-hold (BH) tasks is a non-invasive, robust method to estimate cerebrovascular reactivity (CVR). However, movement and breathing-related artefacts caused by the BH can substantially hinder CVR estimates due to their high temporal collinearity with the effect of interest, and attention has to be paid when choosing which analysis model should be applied to the data. In this study, we evaluate the performance of multiple analysis strategies based on lagged general linear models applied on multi-echo BOLD fMRI data, acquired in ten subjects performing a BH task during ten sessions, to obtain subject-specific CVR and haemodynamic lag estimates. The evaluated approaches range from conventional regression models, i.e. including drifts and motion timecourses as nuisance regressors, applied on single-echo or optimally-combined data, to more complex models including regressors obtained from multi-echo independent component analysis with different grades of orthogonalization in order to preserve the effect of interest, i.e. the CVR. We compare these models in terms of their ability to make signal intensity changes independent from motion, as well as the reliability as measured by voxelwise intraclass correlation coefficients of both CVR and lag maps over time. Our results reveal that a conservative independent component analysis model applied on the optimally-combined multi-echo fMRI signal offers the largest reduction of motion-related effects in the signal, while yielding reliable CVR amplitude and lag estimates, although a conventional regression model applied on the optimally-combined data results in similar estimates. This work demonstrates the usefulness of multi-echo based fMRI acquisitions and independent component analysis denoising for precision mapping of CVR in single subjects based on BH paradigms, fostering its potential as a clinically-viable neuroimaging tool for individual patients. It also proves that the way in which data-driven regressors should be incorporated in the analysis model is not straight-forward due to their complex interaction with the BH-induced BOLD response.

Machine learning-based prediction of MRI-induced power absorption in the tissue in patients with simplified deep brain stimulation lead models.

Interaction of an active electronic implant such as a deep brain stimulation (DBS) system and MRI RF fields can induce excessive tissue heating, limiting MRI accessibility. Efforts to quantify RF heating mostly rely on electromagnetic (EM) simulations to assess individualized specific absorption rate (SAR), but such simulations require extensive computational resources. Here, we investigate if a predictive model using machine learning (ML) can predict the local SAR in the tissue around tips of implanted leads from the distribution of the tangential component of the MRI incident electric field, Etan. A dataset of 260 unique patient-derived and artificial DBS lead trajectories was constructed, and the 1 g-averaged SAR, 1gSARmax, at the lead-tip during 1.5 T MRI was determined by EM simulations. Etan values along each lead's trajectory and the simulated SAR values were used to train and test the ML algorithm. The resulting predictions of the ML algorithm indicated that the distribution of Etan could effectively predict 1gSARmax at the DBS lead-tip (R = 0.82). Our results indicate that ML has the potential to provide a fast method for predicting MR-induced power absorption in the tissue around tips of implanted leads such as those in active electronic medical devices.

Vascular physiology drives functional brain networks.

We present the first evidence for vascular regulation driving fMRI signals in specific functional brain networks. Using concurrent neuronal and vascular stimuli, we collected 30 BOLD fMRI datasets in 10 healthy individuals: a working memory task, flashing checkerboard stimulus, and CO2 inhalation challenge were delivered in concurrent but orthogonal paradigms. The resulting imaging data were averaged together and decomposed using independent component analysis, and three "neuronal networks" were identified as demonstrating maximum temporal correlation with the neuronal stimulus paradigms: Default Mode Network, Task Positive Network, and Visual Network. For each of these, we observed a second network component with high spatial overlap. Using dual regression in the original 30 datasets, we extracted the time-series associated with these network pairs and calculated the percent of variance explained by the neuronal or vascular stimuli using a normalized R2 parameter. In each pairing, one network was dominated by the appropriate neuronal stimulus, and the other was dominated by the vascular stimulus as represented by the end-tidal CO2 time-series recorded in each scan. We acquired a second dataset in 8 of the original participants, where no CO2 challenge was delivered and CO2 levels fluctuated naturally with breathing variations. Although splitting of functional networks was not robust in these data, performing dual regression with the network maps from the original analysis in this new dataset successfully replicated our observations. Thus, in addition to responding to localized metabolic changes, the brain's vasculature may be regulated in a coordinated manner that mimics (and potentially supports) specific functional brain networks. Multi-modal imaging and advances in fMRI acquisition and analysis could facilitate further study of the dual nature of functional brain networks. It will be critical to understand network-specific vascular function, and the behavior of a coupled vascular-neural network, in future studies of brain pathology.

Multiparametric measurement of cerebral physiology using calibrated fMRI.

The ultimate goal of calibrated fMRI is the quantitative imaging of oxygen metabolism (CMRO2), and this has been the focus of numerous methods and approaches. However, one underappreciated aspect of this quest is that in the drive to measure CMRO2, many other physiological parameters of interest are often acquired along the way. This can significantly increase the value of the dataset, providing greater information that is clinically relevant, or detail that can disambiguate the cause of signal variations. This can also be somewhat of a double-edged sword: calibrated fMRI experiments combine multiple parameters into a physiological model that requires multiple steps, thereby providing more opportunity for error propagation and increasing the noise and error of the final derived values. As with all measurements, there is a trade-off between imaging time, spatial resolution, coverage, and accuracy. In this review, we provide a brief overview of the benefits and pitfalls of extracting multiparametric measurements of cerebral physiology through calibrated fMRI experiments.

Changes in arterial cerebral blood volume during lower body negative pressure measured with MRI.

Cerebral Autoregulation (CA), defined as the ability of the cerebral vasculature to maintain stable levels of blood flow despite changes in systemic blood pressure, is a critical factor in neurophysiological health. Magnetic resonance imaging (MRI) is a powerful technique for investigating cerebrovascular function, offering high spatial resolution and wide fields of view (FOV), yet it is relatively underutilized as a tool for assessment of CA. The aim of this study was to demonstrate the potential of using MRI to measure changes in cerebrovascular resistance in response to lower body negative pressure (LBNP). A Pulsed Arterial Spin Labeling (PASL) approach with short inversion times (TI) was used to estimate cerebral arterial blood volume (CBVa) in eight healthy subjects at baseline and -40mmHg LBNP. We estimated group mean CBVa values of 3.13 ± 1.00 and 2.70 ± 0.38 for baseline and lbnp respectively, which were the result of a differential change in CBVa during -40mmHg LBNP that was dependent on baseline CBVa. These data suggest that the PASL CBVa estimates are sensitive to the complex cerebrovascular response that occurs during the moderate orthostatic challenge delivered by LBNP, which we speculatively propose may involve differential changes in vascular tone within different segments of the arterial vasculature. These novel data provide invaluable insight into the mechanisms that regulate perfusion of the brain, and establishes the use of MRI as a tool for studying CA in more detail.

Ultra-high-field arterial spin labelling MRI for non-contrast assessment of cortical lesion perfusion in multiple sclerosis.

OBJECTIVES: To assess the feasibility of using an optimised ultra-high-field high-spatial-resolution low-distortion arterial spin labelling (ASL) MRI acquisition to measure focal haemodynamic pathology in cortical lesions (CLs) in multiple sclerosis (MS). METHODS: Twelve MS patients (eight female, mean age 50 years; range 35-64 years) gave informed consent and were scanned on a 7 Tesla Philips Achieva scanner. Perfusion data were collected at multiple post-labelling delay times using a single-slice flow-sensitive alternating inversion recovery ASL protocol with a balanced steady-state free precession readout scheme. CLs were identified using a high-resolution Phase-Sensitive Inversion Recovery (PSIR) scan. Significant differences in perfusion within CLs compared to immediately surrounding normal appearing grey matter (NAGMlocal) and total cortical normal appearing grey matter (NAGMcortical) were assessed using paired t-tests. RESULTS: Forty CLs were identified in PSIR scans that overlapped with the ASL acquisition coverage. After excluding lesions due to small size or intravascular contamination, 27 lesions were eligible for analysis. Mean perfusion was 40 ± 25 ml/100 g/min in CLs, 53 ± 12 ml/100 g/min in NAGMlocal, and 53 ± 8 ml/100 g/min in NAGMcortical. CL perfusion was significantly reduced by 23 ± 9% (mean ± SE, p = 0.013) and 26 ± 9% (p = 0.006) relative to NAGMlocal and NAGMcortical perfusion, respectively. CONCLUSION: This is the first ASL MRI study quantifying CL perfusion in MS at 7 Tesla, demonstrating that an optimised ASL acquisition is sensitive to focal haemodynamic pathology previously observed using dynamic susceptibility contrast MRI. ASL requires no exogenous contrast agent, making it a more appropriate tool to monitor longitudinal perfusion changes in MS, providing a new window to study lesion development. KEY POINTS: • Perfusion can be quantified within cortical lesions in multiple sclerosis using an optimised high spatial resolution arterial spin Labelling MRI acquisition at ultra-high-field. • The majority of cortical lesions assessed using arterial spin labelling are hypo-perfused compared to normal appearing grey matter, in agreement with dynamic susceptibility contrast MRI literature. • Arterial spin labelling MRI, which does not involve the injection of a contrast agent, is a safe and appropriate technique for repeat scanning of an individual patient.

Arterial CO2 Fluctuations Modulate Neuronal Rhythmicity: Implications for MEG and fMRI Studies of Resting-State Networks.

UNLABELLED: A fast emerging technique for studying human resting state networks (RSNs) is based on spontaneous temporal fluctuations in neuronal oscillatory power, as measured by magnetoencephalography. However, it has been demonstrated recently that this power is sensitive to modulations in arterial CO2 concentration. Arterial CO2 can be modulated by natural fluctuations in breathing pattern, as might typically occur during the acquisition of an RSN experiment. Here, we demonstrate for the first time the fine-scale dependence of neuronal oscillatory power on arterial CO2 concentration, showing that reductions in alpha, beta, and gamma power are observed with even very mild levels of hypercapnia (increased arterial CO2). We use a graded hypercapnia paradigm and participant feedback to rule out a sensory cause, suggesting a predominantly physiological origin. Furthermore, we demonstrate that natural fluctuations in arterial CO2, without administration of inspired CO2, are of a sufficient level to influence neuronal oscillatory power significantly in the delta-, alpha-, beta-, and gamma-frequency bands. A more thorough understanding of the relationship between physiological factors and cortical rhythmicity is required. In light of these findings, existing results, paradigms, and analysis techniques for the study of resting-state brain data should be revisited. SIGNIFICANCE STATEMENT: In this study, we show for the first time that neuronal oscillatory power is intimately linked to arterial CO2 concentration down to the fine-scale modulations that occur during spontaneous breathing. We extend these results to demonstrate a correlation between neuronal oscillatory power and spontaneous arterial CO2 fluctuations in awake humans at rest. This work identifies a need for studies investigating resting-state networks in the human brain to measure and account for the impact of spontaneous changes in arterial CO2 on the neuronal signals of interest. Changes in breathing pattern that are time locked to task performance could also lead to confounding effects on neuronal oscillatory power when considering the electrophysiological response to functional stimulation.

Potential pitfalls when denoising resting state fMRI data using nuisance regression.

In resting state fMRI, it is necessary to remove signal variance associated with noise sources, leaving cleaned fMRI time-series that more accurately reflect the underlying intrinsic brain fluctuations of interest. This is commonly achieved through nuisance regression, in which the fit is calculated of a noise model of head motion and physiological processes to the fMRI data in a General Linear Model, and the "cleaned" residuals of this fit are used in further analysis. We examine the statistical assumptions and requirements of the General Linear Model, and whether these are met during nuisance regression of resting state fMRI data. Using toy examples and real data we show how pre-whitening, temporal filtering and temporal shifting of regressors impact model fit. Based on our own observations, existing literature, and statistical theory, we make the following recommendations when employing nuisance regression: pre-whitening should be applied to achieve valid statistical inference of the noise model fit parameters; temporal filtering should be incorporated into the noise model to best account for changes in degrees of freedom; temporal shifting of regressors, although merited, should be achieved via optimisation and validation of a single temporal shift. We encourage all readers to make simple, practical changes to their fMRI denoising pipeline, and to regularly assess the appropriateness of the noise model used. By negotiating the potential pitfalls described in this paper, and by clearly reporting the details of nuisance regression in future manuscripts, we hope that the field will achieve more accurate and precise noise models for cleaning the resting state fMRI time-series.

The absolute CBF response to activation is preserved during elevated perfusion: Implications for neurovascular coupling measures.

Functional magnetic resonance imaging (fMRI) techniques in which the blood oxygenation level dependent (BOLD) and cerebral blood flow (CBF) response to a neural stimulus are measured, can be used to estimate the fractional increase in the cerebral metabolic rate of oxygen consumption (CMRO2) that accompanies evoked neural activity. A measure of neurovascular coupling is obtained from the ratio of fractional CBF and CMRO2 responses, defined as n, with the implicit assumption that relative rather than absolute changes in CBF and CMRO2 adequately characterise the flow-metabolism response to neural activity. The coupling parameter n is important in terms of its effect on the BOLD response, and as potential insight into the flow-metabolism relationship in both normal and pathological brain function. In 10 healthy human subjects, BOLD and CBF responses were measured to test the effect of baseline perfusion (modulated by a hypercapnia challenge) on the coupling parameter n during graded visual stimulation. A dual-echo pulsed arterial spin labelling (PASL) sequence provided absolute quantification of CBF in baseline and active states as well as relative BOLD signal changes, which were used to estimate CMRO2 responses to the graded visual stimulus. The absolute CBF response to the visual stimuli were constant across different baseline CBF levels, meaning the fractional CBF responses were reduced at the hyperperfused baseline state. For the graded visual stimuli, values of n were significantly reduced during hypercapnia induced hyperperfusion. Assuming the evoked neural responses to the visual stimuli are the same for both baseline CBF states, this result has implications for fMRI studies that aim to measure neurovascular coupling using relative changes in CBF. The coupling parameter n is sensitive to baseline CBF, which would confound its interpretation in fMRI studies where there may be significant differences in baseline perfusion between groups. The absolute change in CBF, as opposed to the change relative to baseline, may more closely match the underlying increase in neural activity in response to a stimulus.

Is fMRI "noise" really noise? Resting state nuisance regressors remove variance with network structure.

Noise correction is a critical step towards accurate mapping of resting state BOLD fMRI connectivity. Noise sources related to head motion or physiology are typically modelled by nuisance regressors, and a generalised linear model is applied to regress out the associated signal variance. In this study, we use independent component analysis (ICA) to characterise the data variance typically discarded in this pre-processing stage in a cohort of 12 healthy volunteers. The signal variance removed by 24, 12, 6, or only 3 head motion parameters demonstrated network structure typically associated with functional connectivity, and certain networks were discernable in the variance extracted by as few as 2 physiologic regressors. Simulated nuisance regressors, unrelated to the true data noise, also removed variance with network structure, indicating that any group of regressors that randomly sample variance may remove highly structured "signal" as well as "noise." Furthermore, to support this we demonstrate that random sampling of the original data variance continues to exhibit robust network structure, even when as few as 10% of the original volumes are considered. Finally, we examine the diminishing returns of increasing the number of nuisance regressors used in pre-processing, showing that excessive use of motion regressors may do little better than chance in removing variance within a functional network. It remains an open challenge to understand the balance between the benefits and confounds of noise correction using nuisance regressors.

Early anti-correlated BOLD signal changes of physiologic origin.

Negative BOLD signals that are synchronous with resting state fluctuations have been observed in large vessels in the cortical sulci and surrounding the ventricles. In this study, we investigated the origin of these negative BOLD signals by applying a Cued Deep Breathing (CDB) task to create transient hypocapnia and a resultant global fMRI signal decrease. We hypothesized that a global stimulus would amplify the effect in large vessels and that using a global negative (vasoconstrictive) stimulus would test whether these voxels exhibit either inherently negative or simply anti-correlated BOLD responses. Significantly anti-correlated, but positive, BOLD signal changes during respiratory challenges were identified in voxels primarily located near edges of brain spaces containing CSF. These positive BOLD responses occurred earlier than the negative CDB response across most of gray matter voxels. These findings confirm earlier suggestions that in some brain regions, local, fractional changes in CSF volume may overwhelm BOLD-related signal changes, leading to signal anti-correlation. We show that regions with CDB anti-correlated signals coincide with most, but not all, of the regions with negative BOLD signal changes observed during a visual and motor stimulus task. Thus, the addition of a physiological challenge to fMRI experiments can help identify which negative BOLD signals are passive physiological anti-correlations and which may have a putative neuronal origin.

Simultaneous cortical, subcortical, and brainstem mapping of sensory activation.

Non-painful tactile sensory stimuli are processed in the cortex, subcortex, and brainstem. Recent functional magnetic resonance imaging (fMRI) studies have highlighted the value of whole-brain, systems-level investigation for examining pain processing. However, whole-brain fMRI studies are uncommon, in part due to challenges with signal to noise when studying the brainstem. Furthermore, the differentiation of small sensory brainstem structures such as the cuneate and gracile nuclei necessitates high resolution imaging. To address this gap in systems-level sensory investigation, we employed a whole-brain, multi-echo fMRI acquisition at 3T with multi-echo independent component analysis (ME-ICA) denoising and brainstem-specific modeling to enable detection of activation across the entire sensory system. In healthy participants, we examined patterns of activity in response to non-painful brushing of the right hand, left hand, and right foot, and found the expected lateralization, with distinct cortical and subcortical responses for upper and lower limb stimulation. At the brainstem level, we were able to differentiate the small, adjacent cuneate and gracile nuclei, corresponding to hand and foot stimulation respectively. Our findings demonstrate that simultaneous cortical, subcortical, and brainstem mapping at 3T could be a key tool to understand the sensory system in both healthy individuals and clinical cohorts with sensory deficits.

Macrovascular blood flow and microvascular cerebrovascular reactivity are regionally coupled in adolescence.

Cerebrovascular imaging assessments are particularly challenging in adolescent cohorts, where not all modalities are appropriate, and rapid brain maturation alters hemodynamics at both macro- and microvascular scales. In a preliminary sample of healthy adolescents (n=12, 8-25 years), we investigated relationships between 4D flow MRI-derived blood velocity and blood flow in bilateral anterior, middle, and posterior cerebral arteries and BOLD cerebrovascular reactivity in associated vascular territories. As hypothesized, higher velocities in large arteries are associated with an earlier response to a vasodilatory stimulus (cerebrovascular reactivity delay) in the downstream territory. Higher blood flow through these arteries is associated with a larger BOLD response to a vasodilatory stimulus (cerebrovascular reactivity amplitude) in the associated territory. These trends are consistent in a case study of adult moyamoya disease. In our small adolescent cohort, macrovascular-microvascular relationships for velocity/delay and flow/CVR change with age, though underlying mechanisms are unclear. Our work emphasizes the need to better characterize this key stage of human brain development, when cerebrovascular hemodynamics are changing, and standard imaging methods offer limited insight into these processes. We provide important normative data for future comparisons in pathology, where combining macro- and microvascular assessments may better help us prevent, stratify, and treat cerebrovascular disease.

MRI mapping of hemodynamics in the human spinal cord.

Impaired spinal cord vascular function contributes to numerous neurological pathologies, making it important to be able to noninvasively characterize these changes. Here, we propose a functional magnetic resonance imaging (fMRI)-based method to map spinal cord vascular reactivity (SCVR). We used a hypercapnic breath-holding task, monitored with end-tidal CO2 (PETCO2), to evoke a systemic vasodilatory response during concurrent blood oxygenation level-dependent (BOLD) fMRI. SCVR amplitude and hemodynamic delay were mapped at the group level in 27 healthy participants as proof-of-concept of the approach, and then in two highly-sampled participants to probe feasibility/stability of individual SCVR mapping. Across the group and the highly-sampled individuals, a strong ventral SCVR amplitude was initially observed without accounting for local regional variation in the timing of the vasodilatory response. Shifted breathing traces (PETCO2) were used to account for temporal differences in the vasodilatory response across the spinal cord, producing maps of SCVR delay. These delay maps reveal an earlier ventral and later dorsal response and demonstrate distinct gray matter regions concordant with territories of arterial supply. The SCVR fMRI methods described here enable robust mapping of spatiotemporal hemodynamic properties of the human spinal cord. This noninvasive approach has exciting potential to provide early insight into pathology-driven vascular changes in the cord, which may precede and predict future irreversible tissue damage and guide the treatment of several neurological pathologies involving the spine.

Removing motion and physiological artifacts from intrinsic BOLD fluctuations using short echo data.

Differing noise variance across study populations has been shown to cause artifactual group differences in functional connectivity measures. In this study, we investigate the use of short echo time functional MRI data to correct for these noise sources in blood oxygenation level dependent (BOLD)-weighted time series. A dual-echo sequence was used to simultaneously acquire data at both a short (TE=3.3 ms) and a BOLD-weighted (TE=35 ms) echo time. This approach is effectively "free," using dead-time in the pulse sequence to collect an additional echo without affecting overall scan time or temporal resolution. The proposed correction method uses voxelwise regression of the short TE data from the BOLD-weighted data to remove noise variance. In addition to a typical resting state scan, non-compliant behavior associated with patient groups was simulated via increased head motion or physiological fluctuations in 10 subjects. Short TE data showed significant correlations with the traditional motion-related and physiological noise regressors used in current connectivity analyses. Following traditional preprocessing, the extent of significant additional variance explained by the short TE data regressors was significantly correlated with the average head motion across the scan in the resting data (r(2)=0.93, p<0.0001). The reduction in data variance following the inclusion of short TE regressors was also correlated with scan head motion (r(2)=0.48, p=0.027). Task-related data were used to demonstrate the effects of the short TE correction on BOLD activation time series with known temporal structure; the size and strength of the activation were significantly decreased, but it is not clear whether this reflects BOLD contamination in the short TE data or correlated changes in blood volume. Finally, functional connectivity maps of the default mode network were constructed using a seed correlation approach. The effects of short TE correction and low-pass filtering on the resulting correlations maps were compared. Results suggest that short TE correction more accurately differentiates artifactual correlations from the correlations of interest in conditions of amplified noise.