Useful Metrics to Facilitate Clinical Application of Resting-state fMRI
Tie-Qiang Li1, Guochun Fu2, and Peter Aspelin3

1Department of Medical Physics, Karolinska University Hospital, Stockholm, Sweden, 2Department of CLINTEC, Karolinska Institute, Stockholm, Sweden, 3Department of Radiology, Karolinska University Hospital, Stockholm, Sweden

Synopsis

Resting-state fMRI is useful for studying functional networks in the human brain and the abnormalities associated with various neuropsychiatric disorders. However, a lack of quantitative metrics closely associated with underlying neurophysiological characteristics has made its translation to clinical settings difficult. In this study, we propose two voxel-based metrics, namely the functional connection counter index (CCI) and connection strength index (CSI). These metrics depicts high contrast between different brain tissues and can be used for quantitative data-driven analysis to detect changes in functional connectivity related to neuropathology.

PURPOSE

Resting-state fMRI is a popular way to study functional networks in the human brain and the abnormalities associated with various neuropsychiatric disorders. However, a lack of quantitative metrics closely associated with underlying neurophysiology characteristics has made its translation to clinical settings difficult1,2. In this study, we propose two voxel-based metrics derived from resting state fMRI measurements as potential biomarkers for sensitive detection of pathological changes resulting from neurological disorders.

MATERIALS AND METHODS

The study sample consisted of 27 patients Alzheimer’s disease (AD) (m/f =12/15, age=73±6) and 20 healthy controls (m/f =9/11, age=72±7). The resting-state fMRI data acquisition was conducted using a whole-body 3T clinical MRI scanner (Magnetom Trio, Siemens) equipped with a 32-channel phased-array receiving head coil. The protocol included conventional clinical MRI scans and a session of 10min long resting-state fMRI using a 2D GRE EPI technique. The main acquisition parameters included the following: TE/TR 35/1600 ms, flip angle=90°, 26 slices of 4 mm thick, in-plane resolution of 3x3mm2, iPAT=2, and 400 dynamic timeframes3. The resting-state fMRI datasets underwent a preprocessing pipeline based on AFNI1,2. The main steps included temporal de-spiking, rigid body image registration for motion correction, spatial normalization to the standard MNI template, nuisance signal removal by voxel-wise regression using 14 regressors based on the motion correction parameters, the average signal of the ventricles, baseline detrending removal up to the third order polynomial, low-pass filtering at 0.08Hz, and local Gaussian smoothing up to FWHM=5mm using an eroded gray matter mask. For each voxel inside the brain, we compute the Pearson’s cross-correlation coefficients of the resting-state fMRI time course with that of every other voxel inside the brain. To conduct quantitative data-drive analysis (QDA) of the resting-state fMRI data we computed the following voxel-wise metrics: (1) Connectivity counter index (CCI) which is defined as the number of voxel pairs involving the current voxel in question with the Pearson’s cross-correlation coefficients above a given threshold. (2) Connectivity strength index (CSI) which is defined as the non-zero mean value of the Pearson’s cross-correlation coefficients above a given threshold for all voxel pairs involving the current voxel in question. To optimize the threshold we systematically computed CCI and CSI by setting the threshold 0.2-0.4 with an incrementing step of 0.01. We performed t-test with derived CCI and CSI images to detect possible differences between AD and healthy controls. The statistical significance of t-test results was assessed using a family-wise error rate at p<0.05.

RESULTS

As shown in Figure 1, the derived CCI and CSI metrics depicts high contrast between different brain tissues. The functional connection hubs can be directly identified in the CCI images as the bright brain areas showing higher number of functional connections. The CSI images have similar but less tissue contrast. The brain regions with significantly reduced functional connectivity in AD patients included posterior cingulate, precuneus, inferior parietal lobule and prefrontal cortex as detected by t-test of the CCI and CSI metrics (see Figure 2). This is very consistent with the results reported by previous studies using MRI and other brain imaging modalities3. The optimal threshold value to derive CCI and CSI is about 0.3 as indicated by the maximum detected brain volumes (Figure 3).

DISCUSSION

In addition to the AD data shown here, we have used these metrics to investigate functional connectivity changes in patients with mild trauma brain injuries4, migraine5,6, and psychosis. With a robust pre-processing pipeline for the resting-state fMRI raw data, the derived statistical contrast from CCI and CSI metrics are relatively stable. We have used a consistent data acquisition protocol. It is apparent that we have to investigate further how the metrics are influenced by the data acquisition protocol. We also need to test their robustness in other clinical applications. At higher spatial resolution, the demand for computation capacity can become an issue.

CONCLUSION

The proposed CCI and CSI metrics are relatively stable with good sensitivity to functional connectivity changes related to neuropathology. The metrics can be interpreted straightforward as the local functional connectivity density and strength, respectively. Most importantly, these metrics can be used for quantitative data-driven analysis (e.g., t-test or ANVOA) for direct comparison between subjects groups and over different time periods. Furthermore, they can also be used for studying correlations with other biomarkers and behavior measures.

Acknowledgements

This study was supported by research grants from the Swedish Research Council and ALF Medicine in Stockholm Province. The authors also want to acknowledge the support from Karolinska Institute and Karolinska University Hospital.

References

1) Wang Y, Li TQ. Analysis of whole-brain resting-state FMRI data using hierarchical clustering approach. PloS one 8:e76315 (2013). 2) Wang Y, Li TQ. Dimensionality of ICA in resting-state fMRI investigated by feature optimized classification of independent components with SVM. Front Hum Neurosci 9:259 (2015). 3) X. Li; T.Q. LI, N. ANDREASEN, M.K.WIBERG, E. WESTMAN, L.O. WAHLUND, Ratio of Ass42/P-tau(181p) in CSF is associated with aberrant default mode network in AD, Scientific Reports, 3:1339 (2013). 4) L.E. Nordin, M.C. Möller, P.Julinc, A. Bartfaic, F. Hashim, T.Q. Li, Post mTBI Fatigue is associated with abnormal brain functional connectivity, Scientific Reports, in press (2015). 5) J.C. Juto and R. Hallin, Kinetic oscillation stimulation as treatment of acute migraine. A randomized, controlled pilot study. Headache 55:11 (2015). 6) T.Q. Li, Y. Wang and J.E. Juto, fMRI Study of Migraine Treatment with Kinetic Oscillation Stimulation, Neuroimage Clin, revised (2015).

Figures

Figure 1: (a) Typical CCI images in gray scale for a single subject derived with a correlation coefficient threshold=0.3, (b) the corresponding CCI images obtained with the same correlation coefficient threshold=0.3.

Figure 2. The mask in red color overlaid onto the MNI template showing brain regions with significant reduced CSI (FWER, p<0.05) in AD patients. The CSI data derived with the correlation coefficient threshold=0.3.

Figure 3. The detected brain volumes with significant reduced CCI and CSI in AD patients (FWER, p<0.05) as a function of the employed correlation coefficient threshold.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
3820