In MR cerebral vasoreactivity (CVR) experiments, responses to the gas stimuli are expected to be with no major time lag. However, it has been shown that different brain regions could respond with different timing during a vasoreactivity. The presence of potential latencies could lead to a misinterpretation of the resulting CVR maps. We attempt to correct CVR maps for physiological latency differences, and propose an alternative way to display both corrected CVR and CVR latency. The data from a Moya Moya patient highlight that even without a strong perfusion alteration, the vasoreactivity is strongly delayed or even completely disrupted.
To investigate the feasibility and interest of this new approach, we examined a 60 years old subject with an occlusion of the left middle cerebral artery (Figure 1A) and a suspected Moyamoya disease. No arguments for a chronic ischemia, no coronary history, no evidences of previous stroke/transient cerebral ischemia or headaches, and no known arterial risk factors. No major sign of asymmetric perfusion was identified. We also display the results for a healthy subject.
Neuroimaging data were collected at 3T (MAGNETOM Skyra, Siemens Healthcare) using a 32-channel head coil. A hypercapnic challenge was performed during the MR scan to measure CVR. The paradigm of inhalation involved a shift between air and carbon-dioxid-enriched gas (BACTAL®, 8% of carbone dioxid). A MAGLIFE system (Maglife, Schiller medical) measured physiologic parameters such as mean arterial blood pressure, heart rate, and arterial oxygen saturation. The corresponding MRI data were acquired using an EPI sequence (262 volumes, 3 gas blocks). A time-of-flight angiography was included to check for vascular abnormalities. Finally, for the patient, a susceptibility based MR perfusion sequence (Gadolinium based) and a diffusion sequence were performed. Tmax maps were estimated from prefusion images using RAPID software (RAPID 4.7.1.1).
The data was processed with MATLAB (R2014a) using the SPM12 toolbox. The CVR-latency map was created first by estimating the delay between mean end tidal CO2 obtained on a standard population and BOLD signal of each voxel based on the maximum value of the cross-correlation. The corrected CVR maps are created using a two step procedure, first temporal realignment is achieved according to the CVR-latency map, then a pearson-correlation is implemented between etCO2 and each voxel timecourse intensity. The corrected CVR map indicates this correlation value. Finally, a RGB map is generated for the corrected CVR, with the color indicating the delay and the brightness the “confidence” in this result (degree of correlation).
[1] Bandettini, P.A.,Wong, E.C., Hinks, R.S., Tikofsky, R.S., Hyde, J.S.,1992. Time course EPI of human brain function during task activation. Magn. Reson. Med. 25, 390–397.
[2] Kwong, K.K., Belliveau, J.W., Chesler, D.A., Goldberg, I.E.,Weisskoff, R.M., Poncelet, B.P., Kennedy, D.N., Hoppel, B.E., Cohen, M.S., Turner, R., 1992. Dynamic magnetic reso- nance imaging of human brain activity during primary sensory stimulation. Proc. Natl. Acad. Sci. U.S.A. 89, 5675–5679.
[3] Logothetis, N.K., Pauls, J., Augath, M., Trinath, T., Oeltermann, A., 2001. Neurophysio- logical investigation of the basis of the fMRI signal. Nature 412, 150–157.
[4] Ogawa, S., Tank, D.W., Menon, R., Ellermann, J.M., Kim, S.G., Merkle, H., Ugurbil, K.,1992. Intrinsic signal changes accompanying sensory stimulation: functional brain map- ping with magnetic resonance imaging. Proc. Natl. Acad. Sci. U.S.A. 89, 5951–5955.
[5] Bandettini, P.A., 1999. The temporal resolution of MRI. Functional MRI 205–220.
[6] Lee, A.T., Glover, G.H., Meyer, C.H., 1995. Discrimination of large venous vessels in time- course spiral blood-oxygen-level-dependent magnetic-resonance functional neuroimaging. Magn. Reson. Med. 33, 745–754
[7] C. Chang, M. E. Thomason, and G. H. Glover, “Mapping and correction of vascular hemodynamic latency in the BOLD signal,” Neuroimage, vol. 43, no. 1, pp. 90–102, 2008.