Is the local functional connectivity anisotropy (LFCA) in white matter caused by neuro-electric activity?  An examination of potential confounds for orientation-dependent LFCA in fMRI
Michael J. Tobia1, David Gallagher1, Rahul Dewal1, Sebastien Rupprecht1, Prasanna Karunanayaka1, and Qing X. Yang1

1Radiology, Penn State Hershey, Hershey, PA, United States

Synopsis

Phantom experiments showed that fluctuating electric current is sufficient to generate local functional connectivity anisotropy (LFCA), and that effects of motion, such as Lorentz forces, cannot explain the alignment of eigenvectors through neighboring voxels or B0 orientation-dependence. In conclusion, anisotropic correlations of fMRI time series may arise from an alternative non-BOLD contrast mechanism, potentially related to an electric current effect on B0.

Introduction/Purpose

Typically, T2*-weighted GE-EPI is sensitive to the blood-oxygenation-level-dependent (BOLD) effect, which is coupled to gray matter (GM) neuronal activity in myriad complex ways [1], and bears little relation to signals in white matter (WM). Yet WM signals yield functional characteristics including low frequency fluctuations and anisotropic temporal correlations [2] that are orientation-dependent with respect to the direction of B0 [3], which suggested a possible hypothesis that local functional connectivity anisotropy (LFCA) could be related to electric current transmitting along WM fiber pathways. We tested this hypothesis in principle with phantom studies and controlled variables that may confound the interpretation, such as correlation bias between through-plane and in-plane voxel-pairs, motion via Lorentz forces, and electric current intensity dependence. This research addressed these confounds and the sensitivity of LFCA method in determining the mechanism that gives rise to B0 orientation-dependent anisotropic local temporal correlations in WM resting state fMRI time series.

Method

T2*-weighted fMRI (TE/TR 31/2000 ms) were recorded from a spherical distilled water-filled phantom with either 1 mA or .1 mA DC fluctuating at approximately .05 Hz and oriented either parallel or perpendicular to B0. Data were preprocessed according to typical resting state fMRI protocol. Temporal correlations were computed in a neighborhood radius of 2 voxels, and an isotropic tensor model was fitted to the local functional connectivity matrix. Motion correction parameters were calculated using an industry standard realignment algorithm. Detrended motion parameters were regressed out of the phantom data and anisotropy and Eigen values from the LFCA tensor were compared before and after. In addition, phantom data were recorded with the readout gradient in the axial (Bz), sagittal (Bx) or coronal (By) planes, and with the electric current oriented parallel to either Bz or Bx.

Results

Both current intensities produced paradigm-correlated motion in the phantom when the current was oriented perpendicular to B0, but not when oriented parallel to B0 (figure 1). The 1 mA generated detectable motion in 5 of 6 directions (roll, pitch, yaw, dS & dL), with displacement of approximately 0.02 mm rotation, and .03 mm of translational movement. The 0.1 mA also generated detectable rotational motion in the roll and the yaw (approximately .02 degrees). These results clearly show the effect of the Lorentz force. For both current intensities, the motion was significantly negatively correlated with the raw (detrended) signal from a sample voxel, and all positive correlations were non-significant. This shows that when the current is ‘active’ the fMRI signal drops out and global displacement of the volume is detected. After correction of the motion, however, the relation between the paradigm time course and fMRI signal became positive, and anisotropy (and 1st Eigen values) remained highly significant (figure 2). In fact, the 1st Eigen value was larger for the 0.1 mA after correction of the global motion caused the Lorentz force. The temporal correlations between the nearest voxel-pairs are indeed stronger for in-plane directions than between slices, which biased eigenvector alignment. Importantly, such bias, along with LFCA, was not observed when the electric current was oriented parallel to B0, which replicated our previous in vivo findings. Finally, the weaker (0.1 mA) current produced a weaker global movement, however, it accounted for as much of the signal fluctuation as the stronger current.

Discussion

Local anisotropic temporal correlations in a phantom are related to electric current after controlling for several confounds. Fluctuating electric current caused detectable global motion of the phantom, which contributed to but cannot fully account for local correlation anisotropy, nor the corresponding Eigen values and their eigenvector alignment. The B0 orientation-dependence of local temporal correlations cannot be explained by the slice selection orientation of the scan and thus, can be ruled out as a possible cause of the absence of LFCA in the B0 direction. Despite causing a weak signal change, fMRI is sensitive to detect the fluctuation associated with electric current. Together, these results support the hypothesis that orientation-dependent anisotropic local temporal correlations in vivo are related to neuroelectric current in the WM.

Conclusion

Fluctuating electric current is a most likely cause for the correlated signal changes through neighboring voxels that appear as orientation-dependent local functional connectivity anisotropy in white matter withT2*-weighted fMRI.

Acknowledgements

No acknowledgement found.

References

1. Logothetis NK & Wandell BA (2004). Interpreting the BOLD signal. Annual Review of Phsyiology, 66:735-69.

2. Ding et al (2015). Visualizing functional pathways in the human brian using correlation tensors and magnetic resonance imaging. Magnetic Resonance Imaging.

3. Tobia et al (2014). Anisotropy of local functional connectivity (LFC) in resting state fMRI time series: what does it say about the fMRI signal? International Society for Magnetic Resonance in Medicine Annual Meeting 2015. Toronto, ON, CA.

Figures

Global motion correlated to the blocked design DC 'activation'. Gray bars indicate the DC was on, and white bars indicate DC was off. Motion was related to this 'activation' structure when DC was perpendicular to B0 (left, center), but not when parallel to B0 (right).

Effects of motion correction on LFCA. Panel A shows anisotropy maps calculated before and after motion correction for DC intensities of 1 mA and 0.1 mA. Panel B shows the effect of motion correction on measured anisotropy. Panel C shows the effect of motion correction on Eigen values from the LFCA tensor.



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