The Nature of the Task Influences Intrinsic Connectivity Networks: An Exploratory fMRI Study in Healthy Subjects

Parcellation-based tractographic modeling of the salience network through meta-analysis

BACKGROUND

The salience network (SN) is a transitory mediator between active and passive states of mind. Multiple cortical areas, including the opercular, insular, and cingulate cortices have been linked in this processing, though knowledge of network connectivity has been devoid of structural specificity.

OBJECTIVE

The current study sought to create an anatomically specific connectivity model of the neural substrates involved in the salience network.

METHODS

A literature search of PubMed and BrainMap Sleuth was conducted for resting-state and task-based fMRI studies relevant to the salience network according to PRISMA guidelines. Publicly available meta-analytic software was utilized to extract relevant fMRI data for the creation of an activation likelihood estimation (ALE) map and relevant parcellations from the human connectome project overlapping with the ALE data were identified for inclusion in our SN model. DSI-based fiber tractography was then performed on publicaly available data from healthy subjects to determine the structural connections between cortical parcellations comprising the network.

RESULTS

Nine cortical regions were found to comprise the salience network: areas AVI (anterior ventral insula), MI (middle insula), FOP4 (frontal operculum 4), FOP5 (frontal operculum 5), a24pr (anterior 24 prime), a32pr (anterior 32 prime), p32pr (posterior 32 prime), and SCEF (supplementary and cingulate eye field), and 46. The frontal aslant tract was found to connect the opercular-insular cluster to the middle cingulate clusters of the network, while mostly short U-fibers connected adjacent nodes of the network.

CONCLUSION

Here we provide an anatomically specific connectivity model of the neural substrates involved in the salience network. These results may serve as an empiric basis for clinical translation in this region and for future study which seeks to expand our understanding of how specific neural substrates are involved in salience processing and guide subsequent human behavior.

Examining the Influence of Spatial Smoothing on Spatiotemporal Features of Intrinsic Connectivity Networks at Low ICA Model Order

Using a relatively high model order of independent component analysis (ICA with 75 ICs) of functional magnetic resonance imaging (fMRI) data, we have reported a clear effect of spatial smoothing Gaussian kernel size on spatiotemporal properties of intrinsic connectivity networks (ICNs). However, many if not the majority of ICA fMRI studies are usually performed at low model order, e.g., 20-IC decomposition, as such low order is generally enough to extract the few networks of interest such as the default-mode network (DMN). The aim of this study is to investigate if we can replicate the spatial smoothing effects on spatiotemporal features of ICNs at low ICA model order. Same resting state fMRI data that we used with 75-IC analysis were used here. Spatial smoothing using an isotropic Gaussian filter kernel with full width at half maximum (FWHM) of 4, 8, and 12 mm was applied during preprocessing. ICNs were identified from 20-IC decomposition and evaluated in terms of three primary features: spatial map intensity, functional network connectivity (FNC), and power spectra. The results identified similar effects of spatial smoothing on spatial map intensities and power spectra at p < 0.01, false discovery rate (FDR) corrected for multiple comparisons. Reduced spatial smoothing kernel size resulted in decreased spatial map intensities as well as a generally decreased low-frequency power (0.01 - 0.10 Hz) but increased high-frequency power (0.15 - 0.25 Hz). FNC, however, did not show a uniform change in correlation values with the size of smoothing kernel. Notably, FNC between DMNs decreased but FNC between central executive and visual networks increased with an increase in smoothing kernel size. These preliminary findings confirm spatial smoothing influences ICN features regardless of model order. The discussion focuses on differences between observed changes at low and high ICA model orders.

The Influence of Spatial Smoothing Kernel Size on the Whole-brain Dynamic Functional Network Connectivity and Meta-state Parameters

In functional magnetic resonance imaging (fMRI), spatial smoothing procedure is generally a stable step in the preprocessing stream. Previous research (including ours) suggested dependency of the static functional connectivity on the size of the spatial smoothing kernel size. But its impact on the time-varying patterns of functional connectivity has not been investigated. Here, we sought to identify the effects of spatial smoothing on brain dynamics by performing dynamic functional network connectivity (dFNC) and meta-state analysis, a unique approach capable of examining a higher-dimensional temporal dynamism of whole-brain functional connectivity. Gaussian smoothing kernel with different widths at half of the maximum of the height of the Gaussian (4, 8, and 12 mm FWHM) were used during preprocessing prior to the group independent component analysis (ICA) with a relatively high model order of 75. dFNC was conducted using the sliding-time window approach and k-means clustering algorithm. Meta-state dynamics method was performed by reducing the number of windowed FNC correlations using principal components analysis (PCA), temporal and spatial ICA and k-means. Results revealed robust effects of spatial smoothing on the connectivity dynamics of several network pairs including a variety of cognitive/attention networks in a connectivity state with the highest occurrence (FDR corrected-p < 0.01). Meta-state analyses indicated significant changes in meta-state metrics including the number of meta-states, meta-state changes, meta-state span, and the total distance. These changes were particularly pronounced when we compared resting state data smoothed with 8 vs. 12 mm FWHM. Our preliminary findings give insights into the effects of spatial smoothing kernel size on the dynamics of functional connectivity and its consequences on meta-state parameters. It also provides further indication of the importance of evaluating variance associated with preprocessing steps on analysis outcomes.

The Influence of Spatial Smoothing Kernel Size on the Temporal Features of Intrinsic Connectivity Networks

Spatial smoothing is a common preprocessing step in the analysis of functional magnetic resonance imaging (fMRI) data. However, little is known about the effect of spatial smoothing kernel size on the temporal properties of functional brain networks. This study presents a pilot investigation on the influence of spatial smoothing using independent component analysis (ICA) as a data-driven technique to extract functional networks of brain in the form of intrinsic connectivity networks (ICNs). BOLD resting state fMRI data were collected from 22 healthy subjects on a 3.0 T MRI scanner. 3D spatial smoothing was applied using a Gaussian filter with full width at half maximum (FWHM) kernel sizes of 4 mm, 8 mm, and 12 mm in the preprocessing step. Group ICA with the Infomax algorithm was performed at 75-IC decomposition. Network temporal features including functional network connectivity (FNC) and BOLD power spectra were calculated and compared pairwise using a paired t-test with a false discovery rate (FDR) correction for multiple comparisons. Results revealed robust effects of smoothing kernel size on FNC measures of most ICNs, largely indicating a decrease in inter-network connectivity as the smoothing kernel size decreased. Power spectra analysis showed increased high-frequency power (0.15 - 0.25 Hz) but decreased low-frequency power (0.01 - 0.10 Hz) with a decrease in the smoothing kernel size (corrected p< 0.01). These findings provide a preliminary observation on the effect of spatial smoothing kernel size on the FNC and power spectra.

An ICA Investigation into the Effect of Physiological Noise Correction on Dynamic Functional Network Connectivity and Meta-state Metrics

Physiological fluctuations such as cardiac pulsations (heart rate) and respiratory rhythm (breathing) have been studied in the resting state functional magnetic resonance imaging (rs-fMRI) studies as the potential sources of confounds in functional connectivity. Independent component analysis (ICA) provides a data driven approach to investigate functional connectivity at the network level. However, the effect of physiological noise correction on the dynamic of ICA-derived networks has not yet been studied. The goal of this study was to investigate the effect of retrospective correction of cardiorespiratory artifacts on the time-varying aspects of functional network connectivity. Blood oxygenation-level dependent (BOLD) rs-fMRI data were collected from healthy subjects using a 3.0T MRI scanner. Whole-brain dynamic functional network connectivity (dFNC) was computed using sliding time window correlation, and k-means clustering of windowed correlation matrices. Results showed significant effects of physiological denoising on dFNC between several network pairs in particular the subcortical, and cognitive/attention networks (false discovery rate [FDR]-corrected p < 0.01). Meta-state dynamics further revealed significant changes in the number of unique windows for each subject, number of times each subject changes from one meta-state to other, and sum of L1 distances between successive meta-states. In conclusion, removal of artifacts is important for achieving reliable fMRI results, however a more cautious approach should be adapted in regressing such "noise" in ICA functional connectivity approach. More experiments are needed to investigate impact of denoising on dFNC especially across different datasets.

An ICA Investigation into the Effect of Physiological Noise Correction on Dimensionality and Spatial Maps of Intrinsic Connectivity Networks

Physiological processes such as cardiac pulsations and respiration can induce signal modulations in functional magnetic resonance imaging (fMRI) time series, and confound inferences made about neural processing from analyses of the blood oxygenation level-dependent (BOLD) signals. Retrospective image space correction of physiological noise (RETROICOR) is a widely used approach to reduce physiological signals in data. Independent component analysis (ICA) is a valuable blind source separation method for analyzing brain networks, referred to as intrinsic connectivity networks (ICNs). Previously, we showed that temporal properties of the ICA-derived networks such as spectral power and functional network connectivity could be impacted by RETROICOR corrections. The goal of this study is to investigate the effect of retrospective correction of physiological artifacts on the ICA dimensionality (model order) and intensities of ICN spatial maps. To this aim, brain BOLD fMRI, heartbeat, and respiration were measured in 22 healthy subjects during resting state. ICA dimensionality was estimated using minimum description length (MDL) based on i.i.d. data samples and smoothness FWHM kernel, and entropy-rate based order selection by finite memory length model (ER-FM) and autoregressive model (ER-AR). Differences in spatial maps between the raw and denoised data were compared using the paired t-test and false discovery rate (FDR) thresholding was used to correct for multiple comparisons. Results showed that ICA dimensionality was greater in the raw data compared to the denoised data. Significant differences were found in the intensities of spatial maps for three ICNs: basal ganglia, precuneus, and frontal network. These preliminary results indicate that the retrospective physiological noise correction can induce change in the resting state spatial map intensity related to the within-network connectivity. More research is needed to understand this effect.