Search

The study of brain network interactions during naturalistic stimuli facilitates a deeper understanding of human brain function. To estimate large-scale brain networks evoked with naturalistic stimuli, a tensor component analysis (TCA) based framework was used to characterize shared spatio-temporal patterns across subjects in a purely data-driven manner. In this framework, a third-order tensor is constructed from the timeseries extracted from all brain regions from a given parcellation, for all participants, with modes of the tensor corresponding to spatial distribution, time series and participants. TCA then reveals spatially and temporally shared components, i.e., evoked networks with the naturalistic stimuli, their time courses of activity and subject loadings of each component. To enhance the reproducibility of the estimation with the adaptive TCA algorithm, a novel spectral clustering method, tensor spectral clustering, was proposed and applied to evaluate the stability of the TCA algorithm. We demonstrated the effectiveness of the proposed framework via simulations and real fMRI data collected during a motor task with a traditional fMRI study design. We also applied the proposed framework to fMRI data collected during passive movie watching to illustrate how reproducible brain networks are evoked by naturalistic movie viewing., Competing Interests: Declaration of Competing Interest The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022. Published by Elsevier Inc.)

Pooling magnetic resonance imaging (MRI) data across research studies, or utilizing shared data from imaging repositories, presents exceptional opportunities to advance and enhance reproducibility of neuroscience research. However, scanner confounds hinder pooling data collected on different scanners or across software and hardware upgrades on the same scanner, even when all acquisition protocols are harmonized. These confounds reduce power and can lead to spurious findings. Unfortunately, methods to address this problem are scant. In this study, we propose a novel denoising approach that implements a data-driven linked independent component analysis (LICA) to identify scanner-related effects for removal from multimodal MRI to denoise scanner effects. We utilized multi-study data to test our proposed method that were collected on a single 3T scanner, pre- and post-software and major hardware upgrades and using different acquisition parameters. Our proposed denoising method shows a greater reduction of scanner-related variance compared with standard GLM confound regression or ICA-based single-modality denoising. Although we did not test it here, for combining data across different scanners, LICA should prove even better at identifying scanner effects as between-scanner variability is generally much larger than within-scanner variability. Our method has great promise for denoising scanner effects in multi-study and in large-scale multi-site studies that may be confounded by scanner differences., (Copyright © 2019. Published by Elsevier Inc.)

Cluster-size tests (CSTs) based on random field theory (RFT) are commonly adopted to identify significant differences in brain images. However, the use of RFT in CSTs rests on the assumption of uniform smoothness (stationarity). When images are non-stationary, CSTs based on RFT will likely lead to increased false positives in smooth regions and reduced power in rough regions. An adjustment to the cluster size according to the local smoothness at each voxel has been proposed for the standard test based on RFT to address non-stationarity, however, this technique requires images with a large degree of spatial smoothing, large degrees of freedom and high intensity thresholding. Recently, we proposed a voxelation-corrected 3D CST based on Gaussian random field theory that does not place constraints on the degree of spatial smoothness. However, this approach is only applicable to stationary images, requiring further modification to enable use for non-stationary images. In this study, we present modifications of this method to develop a voxelation-corrected non-stationary 3D CST based on RFT. Both simulated and real data were used to compare the voxelation-corrected non-stationary CST to the standard cluster-size adjusted non-stationary CST based on RFT and the voxelation-corrected stationary CST. We found that voxelation-corrected stationary CST is liberal for non-stationary images and the voxelation-corrected non-stationary CST performs better than cluster-size adjusted non-stationary CST based on RFT under low smoothness, low intensity threshold and low degrees of freedom., (Published by Elsevier Inc.)

Cluster-size tests (CST) based on random field theory have been widely adopted in fMRI data analysis to detect brain activation. However, most existing approaches can be used appropriately only when the image is highly smoothed in the spatial domain. Unfortunately, spatial smoothing degrades spatial specificity. Recently, a threshold-free cluster enhancement technique was proposed which does not require spatial smoothing, but this method can be used only for group level analysis. Advances in imaging technology now yield high quality high spatial resolution imaging data in single subjects and an inference approach that retains the benefits of greater spatial resolution is called for. In this work, we present a new CST with a correction for voxelation to address this problem. The theoretical formulation of the new approach based on Gaussian random fields is developed to estimate statistical significance using 3D statistical parametric maps without assuming spatial smoothness. Simulated phantom and resting-state fMRI experimental data are then used to compare the voxelation-corrected procedure to the widely used standard random field theory. Unlike standard random field theory approaches, which require heavy spatial smoothing, the new approach has a higher sensitivity for localizing activation regions without the requirement of spatial smoothness., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Independent component analysis (ICA) is widely used in resting state functional connectivity studies. ICA is a data-driven method, which uses no a priori anatomical or functional assumptions. However, as a result, it still relies on the user to distinguish the independent components (ICs) corresponding to neuronal activation, peripherally originating signals (without directly attributable neuronal origin, such as respiration, cardiac pulsation and Mayer wave), and acquisition artifacts. In this concurrent near infrared spectroscopy (NIRS)/functional MRI (fMRI) resting state study, we developed a method to systematically and quantitatively identify the ICs that show strong contributions from signals originating in the periphery. We applied group ICA (MELODIC from FSL) to the resting state data of 10 healthy participants. The systemic low frequency oscillation (LFO) detected simultaneously at each participant's fingertip by NIRS was used as a regressor to correlate with every subject-specific IC time course. The ICs that had high correlation with the systemic LFO were those closely associated with previously described sensorimotor, visual, and auditory networks. The ICs associated with the default mode and frontoparietal networks were less affected by the peripheral signals. The consistency and reproducibility of the results were evaluated using bootstrapping. This result demonstrates that systemic, low frequency oscillations in hemodynamic properties overlay the time courses of many spatial patterns identified in ICA analyses, which complicates the detection and interpretation of connectivity in these regions of the brain., (Copyright © 2013 Elsevier Inc. All rights reserved.)

Networks of brain regions having synchronized fluctuations of the blood oxygen level-dependent functional magnetic resonance imaging (BOLD fMRI) time-series at rest, or "resting state networks" (RSNs), are emerging as a basis for understanding intrinsic brain activity. RSNs are topographically consistent with activity-related networks subserving sensory, motor, and cognitive processes, and studying their spontaneous fluctuations following acute drug challenge may provide a way to understand better the neuroanatomical substrates of drug action. The present within-subject double-blind study used BOLD fMRI at 3T to investigate the functional networks influenced by the non-benzodiazepine hypnotic zolpidem (Ambien). Zolpidem is a positive modulator of γ-aminobutyric acid(A) (GABA(A)) receptors, and engenders sedative effects that may be explained in part by how it modulates intrinsic brain activity. Healthy participants (n=12) underwent fMRI scanning 45 min after acute oral administration of zolpidem (0, 5, 10, or 20mg), and changes in BOLD signal were measured while participants gazed at a static fixation point (i.e., at rest). Data were analyzed using group independent component analysis (ICA) with dual regression and results indicated that compared to placebo, the highest dose of zolpidem increased functional connectivity within a number of sensory, motor, and limbic networks. These results are consistent with previous studies showing an increase in functional connectivity at rest following administration of the positive GABA(A) receptor modulators midazolam and alcohol, and suggest that investigating how zolpidem modulates intrinsic brain activity may have implications for understanding the etiology of its powerful sedative effects., (Copyright © 2013 Elsevier Inc. All rights reserved.)

Confounding noise in BOLD fMRI data arises primarily from fluctuations in blood flow and oxygenation due to cardiac and respiratory effects, spontaneous low frequency oscillations (LFO) in arterial pressure, and non-task related neural activity. Cardiac noise is particularly problematic, as the low sampling frequency of BOLD fMRI ensures that these effects are aliased in recorded data. Various methods have been proposed to estimate the noise signal through measurement and transformation of the cardiac and respiratory waveforms (e.g. RETROICOR and respiration volume per time (RVT)) and model-free estimation of noise variance through examination of spatial and temporal patterns. We have previously demonstrated that by applying a voxel-specific time delay to concurrently acquired near infrared spectroscopy (NIRS) data, we can generate regressors that reflect systemic blood flow and oxygenation fluctuations effects. Here, we apply this method to the task of removing physiological noise from BOLD data. We compare the efficacy of noise removal using various sets of noise regressors generated from NIRS data, and also compare the noise removal to RETROICOR+RVT. We compare the results of resting state analyses using the original and noise filtered data, and we evaluate the bias for the different noise filtration methods by computing null distributions from the resting data and comparing them with the expected theoretical distributions. Using the best set of processing choices, six NIRS-generated regressors with voxel-specific time delays explain a median of 10.5% of the variance throughout the brain, with the highest reductions being seen in gray matter. By comparison, the nine RETROICOR+RVT regressors together explain a median of 6.8% of the variance in the BOLD data. Detection of resting state networks was enhanced with NIRS denoising, and there were no appreciable differences in the bias of the different techniques. Physiological noise regressors generated using Regressor Interpolation at Progressive Time Delays (RIPTiDe) offer an effective method for efficiently removing hemodynamic noise from BOLD data., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Several studies have demonstrated age-related regional differences in the magnitude of the BOLD signal using task-based fMRI. It has been suggested that functional changes reflect either compensatory or de-differentiation mechanisms, both of which assume response to a specific stimulus. Here, we have tested whether ageing affects both task-based and resting brain function, and the extent to which functional changes are mediated by reductions in grey matter (GM) volume. Two groups, of 22 healthy younger and 22 older volunteers, underwent an imaging protocol involving structural and functional MRI, both during a memory task and at rest. The two groups had similar socio-demographical characteristics and cognitive performance. Image analysis revealed both structural and functional differences. Increased BOLD signal in older relative to younger volunteers was mainly observed in the frontal lobes, both during the task and at rest. Functional changes in the frontal lobes were largely located in brain regions spared from GM loss, and adding GM covariates to the fMRI analysis did not significantly alter the group differences. Our results are consistent with the suggestion that, during normal ageing, the brain responds to neuronal loss by fine-tuning connections between spared neurons. Longitudinal studies will be necessary to fully test this hypothesis., (Copyright © 2011 Elsevier Inc. All rights reserved.)

The purpose of this report is to implement novel modifications to overcome the limitations of an existing algorithm for estimating the local statistical noise in a positron emission tomography (PET) image without performing repeated measures. The original algorithm is based on a modification of the filtered back-projection algorithm that allows the variance to be estimated using only a single sinogram. In addition, the effects of photon absorption, random coincidences, radioactive decay, and detector nonuniformity are taken into account. However, there are some limitations when applying this method with modern scanners. In particular, it is common practice to interleave the projections in the sinogram (to increase the sampling rate along each projection) and to perform an interpolation when actually back-projecting to reconstruct the images. Both of these procedures introduce covariance among the elements of the projections, which is cumbersome and impractical to deal with using the existing technique for creating a variance image. An alternative image reconstruction scheme that is shown to be equivalent to image reconstruction using traditional filtered back-projection greatly simplifies the estimation of the variance image. The proposed methods were tested by Monte Carlo simulations and by using repeated scans of a uniform phantom filled with F-18. Results demonstrate that the proposed methods are very rigorous and stable when compared to calculations of the local variance using either repeated measures with a large number of measurements, or region-of-interest estimates of the variance, assuming homogeneous variance structure. In addition, strategies for extending the proposed technique are discussed that would permit the estimation of the variance due to measurement error of a pixel in a brain map from both single subjects and pooled group data.

A technique for calculating the uncertainty in the location of an activation site in a PET image, without performing repeated measures, is presented. With the development of new fMRI methods for measuring cerebral hemodynamics, demonstration of the efficacy of these techniques will be critical to establish clinical utility. Comparisons with PET are a powerful tool for validating these new fMRI techniques. In addition to the fact that PET techniques are well-established methods for making physiological measurements in vivo, PET methods are also free of the geometric distortions and nonuniform signal-to-noise artifacts (due to signal dropout) common in fMRI techniques. Comparisons reported previously have been limited by the large number of trials acquired in single-subject fMRI studies and the small number of trials in a PET study (due to the radiation dose to the patient or the interscan delays for tracer decay). Our method calculates both the center of mass (CM) of a predefined region of interest and the uncertainty in the location of the CM using the preimage PET data (sinograms). Results of phantom studies demonstrate that our method is an unbiased measurement equivalent to that of repeated measures with a large number of images. Extension of this technique to estimate the uncertainty in the location of an activation site in a PET statistical parametric map will permit precise rigorous comparisons of PET and fMRI methods in single subjects without the constraints imposed by the relatively small number of PET measurements.