Frequency-domain patterns in foot-force line-of-action: an emergent property of standing balance control.

A recent line of work suggests that the net behavior of the foot-ground interaction force provides insight into quiet-standing-balance dynamics and control. Through human subject experiments, Boehm et al. found that the relative variations of the center of pressure and force orientation emerge as a distinct pattern in the frequency domain, termed the "intersection-point height." Subsequent empirical and simulation-based studies showed that different control strategies are reflected in the distribution of intersection-point height across frequency. To facilitate understanding of the strengths and limitations of the intersection-point height in describing the dynamics and control of standing, the present work establishes a spectral-based method that also enables derivation of a closed-form estimate of the intersection-point height from any linear model of quiet stance. This new method explained observations from prior work, including how the measure captures aspects of control and physiological noise. The analysis presented herein highlights the utility of the frequency-dependent foot-force dynamics in probing the balance controller and provides a tool for model development and validation to further our understanding of the neuromotor control of natural upright posture in humans.

Improving robustness of 3D multi-shot EPI by structured low-rank reconstruction of segmented CAIPI sampling for fMRI at 7T.

Three-dimensional (3D) encoding methods are increasingly being explored as alternatives to two-dimensional (2D) multi-slice acquisitions in fMRI, particularly in cases where high isotropic resolution is needed. 3D multi-shot EPI acquisition, as the workhorse of 3D fMRI imaging, is susceptible to physiological fluctuations which can induce inter-shot phase variations, and thus reducing the achievable tSNR, negating some of the benefit of 3D encoding. This issue can be particularly problematic at ultra-high fields like 7T, which have more severe off-resonance effects. In this work, we aim to improve the temporal stability of 3D multi-shot EPI at 7T by improving its robustness to inter-shot phase variations. We presented a 3D segmented CAIPI sampling trajectory ("seg-CAIPI") and an improved reconstruction method based on Hankel structured low-rank matrix recovery. Simulation and in-vivo results demonstrate that the combination of the seg-CAIPI sampling scheme and the proposed structured low-rank reconstruction is a promising way to effectively reduce the unwanted temporal variance induced by inter-shot physiological fluctuations, and thus improve the robustness of 3D multi-shot EPI for fMRI.

Unsupervised physiological noise correction of functional magnetic resonance imaging data using phase and magnitude information (PREPAIR).

Of the sources of noise affecting blood oxygen level-dependent functional magnetic resonance imaging (fMRI), respiration and cardiac fluctuations are responsible for the largest part of the variance, particularly at high and ultrahigh field. Existing approaches to removing physiological noise either use external recordings, which can be unwieldy and unreliable, or attempt to identify physiological noise from the magnitude fMRI data. Data-driven approaches are limited by sensitivity, temporal aliasing, and the need for user interaction. In the light of the sensitivity of the phase of the MR signal to local changes in the field stemming from physiological processes, we have developed an unsupervised physiological noise correction method using the information carried in the phase and the magnitude of echo-planar imaging data. Our technique, Physiological Regressor Estimation from Phase and mAgnItude, sub-tR (PREPAIR) derives time series signals sampled at the slice TR from both phase and magnitude images. It allows physiological noise to be captured without aliasing, and efficiently removes other sources of signal fluctuations not related to physiology, prior to regressor estimation. We demonstrate that the physiological signal time courses identified with PREPAIR agree well with those from external devices and retrieve challenging cardiac dynamics. The removal of physiological noise was as effective as that achieved with the most used approach based on external recordings, RETROICOR. In comparison with widely used recording-free physiological noise correction tools-PESTICA and FIX, both performed in unsupervised mode-PREPAIR removed significantly more respiratory and cardiac noise than PESTICA, and achieved a larger increase in temporal signal-to-noise-ratio at both 3 and 7 T.

Ballistocardiogram suppression in concurrent EEG-MRI by dynamic modeling of heartbeats.

The ballistocardiogram (BCG), the induced electric potentials by the head motion originating from heartbeats, is a prominent source of noise in electroencephalography (EEG) data during magnetic resonance imaging (MRI). Although methods have been proposed to suppress the BCG artifact, more work considering the variability of cardiac cycles and head motion across time and subjects is needed to provide highly robust correction. Here, a method called "dynamic modeling of heartbeats" (DMH) is proposed to reduce BCG artifacts in EEG data recorded inside an MRI system. The DMH method models BCG artifacts by combining EEG points at time instants with similar dynamics. The modeled BCG artifact is then subtracted from the EEG recording to suppress the BCG artifact. Performance of DMH was tested and specifically compared with the Optimal Basis Set (OBS) method on EEG data recorded inside a 3T MRI system with either no MRI acquisition (Inside-MRI), echo-planar imaging (EPI-EEG), or fast MRI acquisition using simultaneous multi-slice and inverse imaging methods (SMS-InI-EEG). In a steady-state visual evoked response (SSVEP) paradigm, the 15-Hz oscillatory neuronal activity at the visual cortex after DMH processing was about 130% of that achieved by OBS processing for Inside-MRI, SMS-InI-EEG, and EPI-EEG conditions. The DMH method is computationally efficient for suppressing BCG artifacts and in the future may help to improve the quality of EEG data recorded in high-field MRI systems for neuroscientific and clinical applications.

Evaluating denoising strategies in resting-state functional magnetic resonance in traumatic brain injury (EpiBioS4Rx).

Resting-state functional MRI is increasingly used in the clinical setting and is now included in some diagnostic guidelines for severe brain injury patients. However, to ensure high-quality data, one should mitigate fMRI-related noise typical of this population. Therefore, we aimed to evaluate the ability of different preprocessing strategies to mitigate noise-related signal (i.e., in-scanner movement and physiological noise) in functional connectivity (FC) of traumatic brain injury (TBI) patients. We applied nine commonly used denoising strategies, combined into 17 pipelines, to 88 TBI patients from the Epilepsy Bioinformatics Study for Anti-epileptogenic Therapy clinical trial. Pipelines were evaluated by three quality control (QC) metrics across three exclusion regimes based on the participant's head movement profile. While no pipeline eliminated noise effects on FC, some pipelines exhibited relatively high effectiveness depending on the exclusion regime. Once high-motion participants were excluded, the choice of denoising pipeline becomes secondary - although this strategy leads to substantial data loss. Pipelines combining spike regression with physiological regressors were the best performers, whereas pipelines that used automated data-driven methods performed comparatively worse. In this study, we report the first large-scale evaluation of denoising pipelines aimed at reducing noise-related FC in a clinical population known to be highly susceptible to in-scanner motion and significant anatomical abnormalities. If resting-state functional magnetic resonance is to be a successful clinical technique, it is crucial that procedures mitigating the effect of noise be systematically evaluated in the most challenging populations, such as TBI datasets.

Correlated noise in brain magnetic resonance elastography.

PURPOSE: Magnetic resonance elastography (MRE) uses phase-contrast MRI to generate mechanical property maps of the in vivo brain through imaging of tissue deformation from induced mechanical vibration. The mechanical property estimation process in MRE can be susceptible to noise from physiological and mechanical sources encoded in the phase, which is expected to be highly correlated. This correlated noise has yet to be characterized in brain MRE, and its effects on mechanical property estimates computed using inversion algorithms are undetermined. METHODS: To characterize the effects of signal noise in MRE, we conducted 3 experiments quantifying (1) physiomechanical sources of signal noise, (2) physiological noise because of cardiac-induced movement, and (3) impact of correlated noise on mechanical property estimates. We use a correlation length metric to estimate the extent that correlated signal persists in MRE images and demonstrate the effect of correlated noise on property estimates through simulations. RESULTS: We found that both physiological noise and vibration noise were greater than image noise and were spatially correlated across all subjects. Added physiological and vibration noise to simulated data resulted in property maps with higher error than equivalent levels of Gaussian noise. CONCLUSION: Our work provides the foundation to understand contributors to brain MRE data quality and provides recommendations for future work to correct for signal noise in MRE.

Heart and bladder detection and segmentation on FDG PET/CT by deep learning.

PURPOSE: Positron emission tomography (PET)/ computed tomography (CT) has been extensively used to quantify metabolically active tumors in various oncology indications. However, FDG-PET/CT often encounters false positives in tumor detection due to 18fluorodeoxyglucose (FDG) accumulation from the heart and bladder that often exhibit similar FDG uptake as tumors. Thus, it is necessary to eliminate this source of physiological noise. Major challenges for this task include: (1) large inter-patient variability in the appearance for the heart and bladder. (2) The size and shape of bladder or heart may appear different on PET and CT. (3) Tumors can be very close or connected to the heart or bladder. APPROACH: A deep learning based approach is proposed to segment the heart and bladder on whole body PET/CT automatically. Two 3D U-Nets were developed separately to segment the heart and bladder, where each network receives the PET and CT as a multi-modal input. Data sets were obtained from retrospective clinical trials and include 575 PET/CT for heart segmentation and 538 for bladder segmentation. RESULTS: The models were evaluated on a test set from an independent trial and achieved a Dice Similarity Coefficient (DSC) of 0.96 for heart segmentation and 0.95 for bladder segmentation, Average Surface Distance (ASD) of 0.44 mm on heart and 0.90 mm on bladder. CONCLUSIONS: This methodology could be a valuable component to the FDG-PET/CT data processing chain by removing FDG physiological noise associated with heart and/or bladder accumulation prior to image analysis by manual, semi- or automated tumor analysis methods.

A retrospective physiological noise correction method for oscillating steady-state imaging.

PURPOSE: Oscillating steady-state imaging (OSSI) is an SNR-efficient steady-state sequence with T2∗ sensitivity suitable for FMRI. Due to the frequency sensitivity of the signal, respiration- and drift-induced field changes can create unwanted signal fluctuations. This study aims to address this issue by developing retrospective signal correction methods that utilize OSSI signal properties to denoise task-based OSSI FMRI experiments. METHODS: A retrospective denoising approach was developed that leverages the unique signal properties of OSSI to perform denoising without a manually specified noise region of interest and works with both voxel timecourses (oscillating steady-state correction [OSSCOR]) or FID timecourses (F-OSSCOR). Simulations were performed to estimate the number of principal components optimal for denoising. In vivo experiments at 3 T field strength were conducted to compare the performance of proposed methods against a standard principal component analysis-based method, measured using mean t score within an region of interest, number of activations, and mean temporal SNR. RESULTS: Correction using OSSCOR was significantly better than the standard method in all metrics. Correction using F-OSSCOR was not significantly different from the standard method using an equal number of principal components. Increasing the number of OSSCOR principal components decreased activation strength and increased the number of suspected false positives. However, increasing the number of principal components in F-OSSCOR increased activation strength with little to no increase in false activation. CONCLUSION: Both OSSCOR and F-OSSCOR substantially reduce physiological noise components and increase temporal SNR, improving the functional results of task-based OSSI functional experiments. F-OSSCOR demonstrates a proof of concept utilization of coil-localized FID signal information for physiological noise correction.

FMRI based on transition-band balanced SSFP in comparison with EPI on a high-performance 0.55 T scanner.

PURPOSE: Low-field (<1 tesla) MRI scanners allow more widespread diagnostic use for a range of cardiac, musculoskeletal, and neurological applications. However, the feasibility of performing robust fMRI at low field has yet to be fully demonstrated. To address this gap, we investigated task-based fMRI using a highly sensitive transition-band balanced steady-state free precession approach and standard EPI on a 0.55 tesla scanner equipped with modern high-performance gradient coils and a receive array. METHODS: TR and flip-angle of transition-band steady-state free precession were optimized for 0.55 tesla by simulations. Static shimming was employed to compensate for concomitant field effects. Visual task-based fMRI data were acquired from 8 healthy volunteers. For comparison, standard EPI data were also acquired with TE = T2∗ . Retrospective image-based correction for physiological effects (RETROICOR) was used to quantify physiological noise effects. RESULTS: Activation was robustly detected using both methods in a 4-min scan time. Transition-band steady-state free precession was found to be sensitive to interference from subtle spatial and temporal (field drift, respiration) variations in the magnetic field, counteracting potential advantages of the reduced magnetic susceptibility effects compared to its utilization at high field. These adverse effects could be partially remedied with static shimming and postprocessing approaches. Standard EPI proved more robust against the sources of interference. CONCLUSION: BOLD contrast is sufficiently large at 0.55 tesla for robust detection of brain activation and may be employed to broaden the spectrum of applications of low-field MRI. Standard EPI outperforms transition-band steady-state free precession in terms of signal stability.

Robustness of convolutional neural networks to physiological electrocardiogram noise.

The electrocardiogram (ECG) is a widespread diagnostic tool in healthcare and supports the diagnosis of cardiovascular disorders. Deep learning methods are a successful and popular technique to detect indications of disorders from an ECG signal. However, there are open questions around the robustness of these methods to various factors, including physiological ECG noise. In this study, we generate clean and noisy versions of an ECG dataset before applying symmetric projection attractor reconstruction (SPAR) and scalogram image transformations. A convolutional neural network is used to classify these image transforms. For the clean ECG dataset, F1 scores for SPAR attractor and scalogram transforms were 0.70 and 0.79, respectively. Scores decreased by less than 0.05 for the noisy ECG datasets. Notably, when the network trained on clean data was used to classify the noisy datasets, performance decreases of up to 0.18 in F1 scores were seen. However, when the network trained on the noisy data was used to classify the clean dataset, the decrease was less than 0.05. We conclude that physiological ECG noise impacts classification using deep learning methods and careful consideration should be given to the inclusion of noisy ECG signals in the training data when developing supervised networks for ECG classification. This article is part of the theme issue 'Advanced computation in cardiovascular physiology: new challenges and opportunities'.

Model-based physiological noise removal in fast fMRI.

Recent improvements in the speed and sensitivity of fMRI acquisition techniques suggest that fast fMRI can be used to detect and precisely localize sub-second neural dynamics. This enhanced temporal resolution has enormous potential for neuroscientists. However, physiological noise poses a major challenge for the analysis of fast fMRI data. Physiological noise scales with sensitivity, and its autocorrelation structure is altered in rapidly sampled data, suggesting that new approaches are needed for physiological noise removal in fast fMRI. Existing strategies either rely on external physiological recordings, which can be noisy or difficult to collect, or employ data-driven approaches which make assumptions that may not hold true in fast fMRI. We created a statistical model of harmonic regression with autoregressive noise (HRAN) to estimate and remove cardiac and respiratory noise from the fMRI signal directly. This technique exploits the fact that cardiac and respiratory noise signals are fully sampled (rather than aliasing) when imaging at fast rates, allowing us to track and model physiology over time without requiring external physiological measurements. We then created a joint model of neural hemodynamics, and physiological and autocorrelated noise to more accurately remove noise. We first verified that HRAN accurately estimates cardiac and respiratory dynamics and that our model demonstrates goodness-of-fit in fast fMRI data. In task-driven data, we then demonstrated that HRAN is able to remove physiological noise while leaving the neural signal intact, thereby increasing detection of task-driven voxels. Finally, we established that in both simulations and fast fMRI data HRAN is able to improve statistical inferences as compared with gold-standard physiological noise removal techniques. In conclusion, we created a tool that harnesses the novel information in fast fMRI to remove physiological noise, enabling broader use of the technology to study human brain function.

Functional engagement of white matter in resting-state brain networks.

The topological characteristics of functional networks, derived from measurements of resting-state connectivity in gray matter (GM), are associated with individual cognitive abilities or specific dysfunctions. However, blood oxygen level-dependent (BOLD) signals in white matter (WM) are usually ignored or even regressed out as nuisance factors in the data analyses that underlie network models. Recent studies have demonstrated reliable detection of WM BOLD signals and imply these reflect associated neural activities. Here we evaluate quantitatively the contributions of individual WM voxels to the identification of functional networks, which we term their engagement (or conceptually, their importance). We quantify the engagement by measuring the reductions of connectivity, produced by ignoring the signal fluctuations within each WM voxel, with respect to both the entire network (global) or a single GM node (local). We observed highly reproducible spatial distributions of global engagement maps, as well as a trend toward increased relevance of deep WM voxels at delayed times. Local engagement maps exhibit homogeneous spatial distributions with respect to internal nodes that constitute a well-recognized sub-functional network, but inhomogeneous distributions with respect to other nodes. WM voxels show distinct distributions of engagement depending on their anatomical locations. These findings demonstrate the important role of WM in network modeling, thus supporting the need for changes of conventional views that WM signal variations represent only physiological noise.

Improved physiological noise regression in fNIRS: A multimodal extension of the General Linear Model using temporally embedded Canonical Correlation Analysis.

For the robust estimation of evoked brain activity from functional Near-Infrared Spectroscopy (fNIRS) signals, it is crucial to reduce nuisance signals from systemic physiology and motion. The current best practice incorporates short-separation (SS) fNIRS measurements as regressors in a General Linear Model (GLM). However, several challenging signal characteristics such as non-instantaneous and non-constant coupling are not yet addressed by this approach and additional auxiliary signals are not optimally exploited. We have recently introduced a new methodological framework for the unsupervised multivariate analysis of fNIRS signals using Blind Source Separation (BSS) methods. Building onto the framework, in this manuscript we show how to incorporate the advantages of regularized temporally embedded Canonical Correlation Analysis (tCCA) into the supervised GLM. This approach allows flexible integration of any number of auxiliary modalities and signals. We provide guidance for the selection of optimal parameters and auxiliary signals for the proposed GLM extension. Its performance in the recovery of evoked HRFs is then evaluated using both simulated ground truth data and real experimental data and compared with the GLM with short-separation regression. Our results show that the GLM with tCCA significantly improves upon the current best practice, yielding significantly better results across all applied metrics: Correlation (HbO max. +45%), Root Mean Squared Error (HbO max. -55%), F-Score (HbO up to 3.25-fold) and p-value as well as power spectral density of the noise floor. The proposed method can be incorporated into the GLM in an easily applicable way that flexibly combines any available auxiliary signals into optimal nuisance regressors. This work has potential significance both for conventional neuroscientific fNIRS experiments as well as for emerging applications of fNIRS in everyday environments, medicine and BCI, where high Contrast to Noise Ratio is of importance for single trial analysis.

Extraction of the cardiac waveform from simultaneous multislice fMRI data using slice sorted averaging and a deep learning reconstruction filter.

Cardiac signal contamination has long confounded the analysis of blood-oxygenation-level-dependent (BOLD) functional magnetic resonance imaging (fMRI). Cardiac pulsation results in significant BOLD signal changes, especially in and around blood vessels. Until the advent of simultaneous multislice echo-planar imaging (EPI) acquisition, the time resolution of whole brain EPI was insufficient to avoid cardiac aliasing (and acquisitions with repetition times (TRs) under 400-500 ms are still uncommon). As a result, direct detection and removal of the cardiac signal with spectral filters is generally not possible. Modelling methods have been developed to mitigate cardiac contamination, and recently developed techniques permit the visualization of cardiac signal propagation through the brain in undersampled data (e.g., TRs > 1s), which is useful in its own right for finding blood vessels. However, both of these techniques require data from which to estimate cardiac phase, which is generally not available for the data in many large databases of existing imaging data, and even now is not routinely recorded in many fMRI experiments. Here we present a method to estimate the cardiac waveform directly from a multislice fMRI dataset, without additional physiological measurements, such as plethysmograms. The pervasive spatial extent and temporal structure of the cardiac contamination signal across the brain offers an opportunity to exploit the nature of multislice imaging to extract this signal from the fMRI data itself. While any particular slice is recorded at the TR of the imaging experiment, slices are recorded much more quickly - typically from 10 to 20 Hz - sufficiently fast to fully sample the cardiac signal. Using the fairly permissive assumptions that the cardiac signal is a) pseudoperiodic b) somewhat coherent within any given slice, and c) is similarly shaped throughout the brain, we can extract a good estimate of the cardiac phase as a function of time from fMRI data alone. If we make further assumptions about the shape and consistency of cardiac waveforms, we can develop a deep learning filter to greatly improve our estimate of the cardiac waveform.

Identification of physiological response functions to correct for fluctuations in resting-state fMRI related to heart rate and respiration.

Functional magnetic resonance imaging (fMRI) is widely viewed as the gold standard for studying brain function due to its high spatial resolution and non-invasive nature. However, it is well established that changes in breathing patterns and heart rate strongly influence the blood oxygen-level dependent (BOLD) fMRI signal and this, in turn, can have considerable effects on fMRI studies, particularly resting-state studies. The dynamic effects of physiological processes are often quantified by using convolution models along with simultaneously recorded physiological data. In this context, physiological response function (PRF) curves (cardiac and respiratory response functions), which are convolved with the corresponding physiological fluctuations, are commonly employed. While it has often been suggested that the PRF curves may be region- or subject-specific, it is still an open question whether this is the case. In the present study, we propose a novel framework for the robust estimation of PRF curves and use this framework to rigorously examine the implications of using population-, subject-, session- and scan-specific PRF curves. The proposed framework was tested on resting-state fMRI and physiological data from the Human Connectome Project. Our results suggest that PRF curves vary significantly across subjects and, to a lesser extent, across sessions from the same subject. These differences can be partly attributed to physiological variables such as the mean and variance of the heart rate during the scan. The proposed methodological framework can be used to obtain robust scan-specific PRF curves from data records with duration longer than 5 min, exhibiting significantly improved performance compared to previously defined canonical cardiac and respiration response functions. Besides removing physiological confounds from the BOLD signal, accurate modeling of subject- (or session-/scan-) specific PRF curves is of importance in studies that involve populations with altered vascular responses, such as aging subjects.

Dual regression physiological modeling of resting-state EPI power spectra: Effects of healthy aging.

Aging and disease-related changes in the arteriovasculature have been linked to elevated levels of cardiac cycle-induced pulsatility in the cerebral microcirculation. Functional magnetic resonance imaging (fMRI), acquired fast enough to unalias the cardiac frequency contributions, can be used to study these physiological signals in the brain. Here, we propose an iterative dual regression analysis in the frequency domain to model single voxel power spectra of echo planar imaging (EPI) data using external recordings of the cardiac and respiratory cycles as input. We further show that a data-driven variant, without external physiological traces, produces comparable results. We use this framework to map and quantify cardiac and respiratory contributions in healthy aging. We found a significant increase in the spatial extent of cardiac modulated white matter voxels with age, whereas the overall strength of cardiac-related EPI power did not show an age effect.

Intracortical smoothing of small-voxel fMRI data can provide increased detection power without spatial resolution losses compared to conventional large-voxel fMRI data.

Continued improvement in MRI acquisition technology has made functional MRI (fMRI) with small isotropic voxel sizes down to 1 mm and below more commonly available. Although many conventional fMRI studies seek to investigate regional patterns of cortical activation for which conventional voxel sizes of 3 mm and larger provide sufficient spatial resolution, smaller voxels can help avoid contamination from adjacent white matter (WM) and cerebrospinal fluid (CSF), and thereby increase the specificity of fMRI to signal changes within the gray matter. Unfortunately, temporal signal-to-noise ratio (tSNR), a metric of fMRI sensitivity, is reduced in high-resolution acquisitions, which offsets the benefits of small voxels. Here we introduce a framework that combines small, isotropic fMRI voxels acquired at 7 T field strength with a novel anatomically-informed, surface mesh-navigated spatial smoothing that can provide both higher detection power and higher resolution than conventional voxel sizes. Our smoothing approach uses a family of intracortical surface meshes and allows for kernels of various shapes and sizes, including curved 3D kernels that adapt to and track the cortical folding pattern. Our goal is to restrict smoothing to the cortical gray matter ribbon and avoid noise contamination from CSF and signal dilution from WM via partial volume effects. We found that the intracortical kernel that maximizes tSNR does not maximize percent signal change (ΔS/S), and therefore the kernel configuration that optimizes detection power cannot be determined from tSNR considerations alone. However, several kernel configurations provided a favorable balance between boosting tSNR and ΔS/S, and allowed a 1.1-mm isotropic fMRI acquisition to have higher performance after smoothing (in terms of both detection power and spatial resolution) compared to an unsmoothed 3.0-mm isotropic fMRI acquisition. Overall, the results of this study support the strategy of acquiring voxels smaller than the cortical thickness, even for studies not requiring high spatial resolution, and smoothing them down within the cortical ribbon with a kernel of an appropriate shape to achieve the best performance-thus decoupling the choice of fMRI voxel size from the spatial resolution requirements of the particular study. The improvement of this new intracortical smoothing approach over conventional surface-based smoothing is expected to be modest for conventional resolutions, however the improvement is expected to increase with higher resolutions. This framework can also be applied to anatomically-informed intracortical smoothing of higher-resolution data (e.g. along columns and layers) in studies with prior information about the spatial structure of activation.

On the analysis of rapidly sampled fMRI data.

Recent advances in parallel imaging and simultaneous multi-slice techniques have permitted whole-brain fMRI acquisitions at sub-second sampling intervals, without significantly sacrificing the spatial coverage and resolution. Apart from probing brain function at finer temporal scales, faster sampling rates may potentially lead to enhanced functional sensitivity, owing possibly to both cleaner neural representations (due to less aliased physiological noise) and additional statistical benefits (due to more degrees of freedom for a fixed scan duration). Accompanying these intriguing aspects of fast acquisitions, however, confusion has also arisen regarding (1) how to preprocess/analyze these fast fMRI data, and (2) what exactly is the extent of benefits with fast acquisitions, i.e., how fast is fast enough for a specific research aim? The first question is motivated by the altered spectral distribution and noise characteristics at short sampling intervals, while the second question seeks to reconcile the complicated trade-offs between the functional contrast-to-noise ratio and the effective degrees of freedom. Although there have been recent efforts to empirically approach different aspects of these two questions, in this work we discuss, from a theoretical perspective accompanied by some illustrative, proof-of-concept experimental in vivo human fMRI data, a few considerations that are rarely mentioned, yet are important for both preprocessing and optimizing statistical inferences for studies that employ acquisitions with sub-second sampling intervals. Several summary recommendations include concerns regarding advisability of relying on low-pass filtering to de-noise physiological contributions, employment of statistical models with sufficient complexity to account for the substantially increased serial correlation, and cautions regarding using rapid sampling to enhance functional sensitivity given that different analysis models may associate with distinct trade-offs between contrast-to-noise ratios and the effective degrees of freedom. As an example, we demonstrate that as TR shortens, the intrinsic differences in how noise is accommodated in general linear models and Pearson correlation analyses (assuming Gaussian distributed stochastic signals and noise) can result in quite different outcomes, either gaining or losing statistical power.

Cluster failure revisited: Impact of first level design and physiological noise on cluster false positive rates.

Methodological research rarely generates a broad interest, yet our work on the validity of cluster inference methods for functional magnetic resonance imaging (fMRI) created intense discussion on both the minutia of our approach and its implications for the discipline. In the present work, we take on various critiques of our work and further explore the limitations of our original work. We address issues about the particular event-related designs we used, considering multiple event types and randomization of events between subjects. We consider the lack of validity found with one-sample permutation (sign flipping) tests, investigating a number of approaches to improve the false positive control of this widely used procedure. We found that the combination of a two-sided test and cleaning the data using ICA FIX resulted in nominal false positive rates for all data sets, meaning that data cleaning is not only important for resting state fMRI, but also for task fMRI. Finally, we discuss the implications of our work on the fMRI literature as a whole, estimating that at least 10% of the fMRI studies have used the most problematic cluster inference method (p = .01 cluster defining threshold), and how individual studies can be interpreted in light of our findings. These additional results underscore our original conclusions, on the importance of data sharing and thorough evaluation of statistical methods on realistic null data.

Influence of background suppression and retrospective realignment on free-breathing renal perfusion measurement using pseudo-continuous ASL.

PURPOSE: To assess the influence of background suppression and retrospective realignment on physiological noise and image quality in free-breathing renal pseudo-continuous arterial spin labeling (pCASL). METHODS: Ten subjects were scanned at 3T with a pCASL prepared single-slice coronal acquisition through the kidneys under free breathing. Multiple acquisitions were performed with various levels of residual background signal based on optimization of pulse timings to achieve specific background suppression levels (<2%, <5%, <10%, <20%). A retrospective non-rigid motion-correction strategy was also implemented. RESULTS: Decreasing level of residual background signal was associated with higher temporal SNR. The retrospective motion-correction provided an additional but not statistically significant improvement in tSNR. The highest image quality was obtained with the lowest level of residual background signal accompanied by the retrospective motion-correction, although no significant difference in quantitative renal blood-flow could be observed. CONCLUSIONS: Renal perfusion measurement with ASL under free breathing is feasible and robust against physiological noise when using strong background suppression strategies. Finally, retrospective motion-correction further improves image quality but cannot replace background suppression.