Effective Connectivity Measured with Layer-Dependent Resting-State Blood Volume fMRI in Humans
Laurentius Huber1, Daniel A Handwerker1, Javier Gonzalez-Castillo1, David C Jangraw1, Maria Guidi2, Dimo Ivanov3, Benedikt A Poser3, Jozien Goense4, and Peter A Bandettini1

1Section of Functional Imaging Methods, National Institute of Mental Health, Bethesda, MD, United States, 2Human Cognitive and Brain Sciences, Max Planck Institute, Leipzig, Germany, 3MBIC, Maastricht University, Maastricht, Netherlands, 4Institute of Neuroscience and Psychology, University of Glasgow, Glasgow, United Kingdom

Synopsis

Measurements of layer-dependent cortical activity provide insight on how feedforward/feedback functional connectivity affects a given cortical area. Here, we simultaneously measure layer-dependent changes in resting-state BOLD and CBV with VASO. We demonstrate that the superior specificity of CBV fMRI reveals layer-dependent resting-state activity better than GE-BOLD fMRI and gives indications of effective connectivity in the human sensory-motor system. In particular, superficial and deeper layers in M1 show different connectivity patterns than those associated with the middle layer, likely driven by input from S1. Our data show that the middle layer in M1 correlates with contralateral M1, while it anti-correlates with contralateral S1.

Purpose

Cortical layer-dependent high-resolution fMRI can help understand the effects of afferent and efferent functional connectivity on brain activity at different cortical layers [1]. Because of known differences in input-output characteristics for anatomically defined cortical layers, we hypothesize that individual layers should show different resting-state signal fluctuations [2]. Here, we report on a new fMRI acquisition and analysis method that might be able to measure layer-dependent resting-state activity fluctuations. In addition, we show effective functional connectivity of layer-dependent resting-state signal fluctuations in the human sensory-motor system.

Methods

Experiments in four volunteers were performed on a 7T Siemens scanner with a 32-channel NOVA Medical head coil. We used a multi-contrast VASO-BOLD sequence with slices positioned to be nearly perpendicular to the cortex to capture M1 and S1 of both hemispheres without aliasing artifacts in the A-P direction (Fig. 1A). Sequence parameters (similar to [3]), TI1/TI2/TR = 1.1/2.6/3.0 s, 7 slices with GRAPPA-2 and FLEET [4] without partial-Fourier imaging, rectangular imaging matrix with 140 phase encoding steps and 70 data points. Two functional scans were collected per subject: (1) 12 min long rest, (2) 12 min long unilateral finger-tapping (alternating 30 s rest and 30 s act). Physiological noise correction was performed with RETROICOR based on acquired respiration and cardiac traces. Segmentation and cortical layering algorithms (based on the equi-volume approach [5]) were implemented in C++ and applied directly on the EPI images from the functional scans (Fig. 1). In order to obtain smooth surfaces, EPI data were rescaled by expanding k-space by a factor of 4 with zero-filling (blue curved arrow in Fig. 1B-C). This allows us to model 10-15 layers across the cortical depth (4 mm in M1; color edges in Fig 1D). See also Fig. 2I for an example of smooth surfaces in rescaled EPI space.

Results

Functional VASO EPI time-series provide sufficient T1-weighted contrast to apply surface/layering algorithms [5] directly on the functional data (Fig. 1D), obviating the need for coregistration and distortion-correction. tSNR of VASO (18) is about 48% of BOLD tSNR (38). However, CBV-weighted functional activity (purple in Fig. 1F/H) shows higher specificity to GM compared to BOLD (orange in Fig. 1E/G) for finger-tapping and specific ICA components (manually selected based on its position in M1). Placing seeds in individual layers of M1 shows distinct layer-specific correlation with numerous brain regions within the FOV (Fig. 1I). Layer-dependent correlation matrices of right-hand M1 with itself (Fig. 2A/E), as well as with ipsilateral S1 (B/F), contralateral M1 (C/G) and contralateral S1 (D/H), confirm that CBV-based fMRI reveals layer-dependent features, which are blurred and shifted towards the cortical surface in BOLD-fMRI.

Discussion

The higher localization specificity of CBV-fMRI compared to BOLD-fMRI is consistent with previous task-based VASO studies [1,3]. Specific input into superficial or deeper cortical layers of M1 is different between the individual brain areas. While input from primary somato-sensory areas correlates most strongly with upper layers in M1 (yellow arrows in Fig. 1I), the input from premotor regions is mostly correlated with deeper layers (blue arrows in Fig. 1I). The cross-correlation matrices (Fig. 2A-G) show layer-dependent features in M1, but not in S1. The presence of layer-dependent features in M1 but not in S1 suggests that the nominal in-plane resolution of 0.9 mm is sufficient to reveal layer-dependent resting-state activity in thick cortices (e.g. M1, 4 mm), but may not be sufficient for thinner cortices (e.g. S1, 1.8-2 mm). However, it might also be the case the individual layers in S1 are more synchronized, compared to M1. The fact that S1 projects its input into M1 mostly in upper layers in M1 (black ellipse in Fig. 1I and 2F) with a secondary peak in deeper layers (blue ellipse in Fig. 2F) is consistent with MEMRI tracer studies in rats [6].

Conclusion

We were able to show layer-dependent resting-state signal features with CBV-based fMRI in humans. The higher specificity of CBV fMRI reveals layer-dependent signal variations more clearly than BOLD. We could identify that both M1 areas have the same source of resting-state signal fluctuations in their middle layer (green lines in Fig. 2), which is anticorrelated with the primary sensory areas (blue lines in Fig. 2).

Acknowledgements

Supported by the NIMH Intramural Research Program. MG was supported by the FP7 Marie-Curie-Actions of the European Commission FP7-PEOPLE-2012-ITN-316716.

References

[1] Huber et al., NeuroImage, 2015, 107:23-33; [2] Polimeni et al., ISMRM, 2010,#353; [3] Huber et al., NeuroImage, 2015, doi:10.1016/j.neuroimage.2015.10.082; [4] Polimeni et al., MRM, 2015, doi:10.1002/mrm.25628; [5] Waehnert et al., NeuroImage, 2014, 93:210-220; [6] Yu et al., Nature Methods, 2014, 11:55-58.

Figures

(A-D) T1-weighting in functional EPI data is used for the definition of cortical layers. (E-H) Task driven and resting-state activity maps. CBV fMRI is more specific but has a lower CNR than BOLD fMRI. (I) Map of maximum correlation to M1 layers. Colors in (I) represent the layer to which they correlate most strongly (independent of the sign). Blue refers to stronger correlation to deep layers of M1-seed and yellow refers to stronger correlation to superficial layers to M1-seed.

Layer-dependent cross correlation matrices between M1 and S1 cortices for BOLD (A-D) and CBV (E-H) resting-state signal fluctuations. There is a layer-dependent signature visible in CBV resting-state (indicated with dashed lines). (I) depicts the ROIS and the approximate cortical depths.



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