Abstract
BACKGROUND AND PURPOSE: Chronic hemodynamic impairment in high-grade carotid occlusive disease is thought to cause microstructural abnormalities that might be subclinical or lead to subtle symptoms including cognitive impairment. Quantitative MR imaging allows assessing pathologic structural changes beyond macroscopically visible tissue damage. In this study, high-resolution quantitative T2 mapping combined with DSC-based PWI was used to investigate quantitative T2 changes as a potential marker of microstructural damage in relation to hemodynamic impairment in patients with unilateral high-grade carotid occlusive disease.
MATERIALS AND METHODS: Eighteen patients with unilateral high-grade ICA or MCA stenosis/occlusion were included in the study. T2 values and deconvolved perfusion parameters, including relative CBF, relative CBV, and the relative CBF/relative CBV ratio as a potential indicator of local cerebral perfusion pressure, were determined within areas with delayed TTP and compared with values from contralateral unaffected areas after segmentation of normal-appearing hypoperfused WM and cortical regions. Hemispheric asymmetry indices were calculated for all parameters.
RESULTS: Quantitative T2 was significantly prolonged (P < .01) in hypoperfused tissue and correlated significantly (P < .01) with TTP delay and relative CBF/relative CBV reduction in WM. Significant correlations (P < .001) between TTP delay and the relative CBF/relative CBV ratio were found both in WM and in cortical areas.
CONCLUSIONS: Quantitative T2 can be used as a marker of microstructural tissue damage even in normal-appearing GM and WM within a vascular territory affected by high-grade carotid occlusive disease. Furthermore, the extent of damage correlates with the degree of hemodynamic failure measured by DSC perfusion parameters.
ABBREVIATIONS:
- AI
- asymmetry index
- CPP
- cerebral perfusion pressure
- q
- quantitative
- rCBF
- relative CBF
- rCBV
- relative CBV
In patients with ICA and MCA occlusive disease, chronic hemodynamic compromise is associated with an increased risk of ischemic stroke1⇓⇓–4 because it may lead to hemodynamic infarction and it presumably also predisposes the patient to embolic phenomena.5 Apart from focal neurologic symptoms due to acute ischemia, chronic hemodynamic impairment in unilateral high-grade ICA stenosis or occlusion may be associated with measurable global cognitive impairment,6,7 even in patients without visible ischemic lesions on conventional MR images (“silent”).7 Pathophysiologically, these clinical findings seem to be related to microstructural changes, including diffuse demyelination, the loss of axons, and gliosis.6,8⇓⇓⇓–12
On conventional MR imaging, ipsilateral atrophy and signal increase on T2WI due to increased water content and gliotic changes have been described in severe chronic hypoperfusion.13 However, in many patients, no abnormalities are visible on conventional T2WI, though microstructural damage associated with stenosis-related chronic hypoperfusion might be reasonably assumed. Compared with conventional MR imaging, quantitative (q)MR imaging techniques are more sensitive to microstructural and metabolic tissue changes.14,15 In patients with diffuse chronic hypoxia (ischemic leukoaraiosis), quantitative MR imaging including qT2 revealed significant T2 prolongation in the normal-appearing WM, suggesting microstructural damage.14 Therefore, qT2 might detect subtle changes, including microstructural damage in the territory of a chronically stenotic cerebral artery that is not visible on conventional MR imaging, and might be a promising imaging biomarker for microstructural changes related to chronic hypoperfusion.
Because these microstructural changes possibly depend on the degree of hemodynamic failure and cerebral perfusion pressure (CPP) impairment,8 it might be interesting to investigate the correlation between potential qT2 changes in the dependent territory and hemodynamic measures, including the CPP. By means of PET, the CBF/CBV ratio has been validated as an indicator of local CPP,16 and a negative correlation with DSC-based TTP delay has been demonstrated in patients with unilateral ICA or MCA stenosis/occlusion.17 Consequently, if this relationship can be reproduced, DSC-based relative (r)CBF/rCBV might be a suitable index for local CPP as well.
In this study, an advanced high-resolution qMRI technique combined with DSC perfusion imaging was used to explore qT2 changes in hypoperfused tissue in patients with high-grade ICA and MCA occlusive disease. The purpose was to evaluate the following: 1) whether qT2 mapping displays a significant increase of T2 values in chronically hypoperfused noninfarcted and normal-appearing tissue as a sign of microstructural damage, and 2) whether these changes of qT2 are related to measures of hemodynamic impairment, including autoregulatory capacity and CPP.
Materials and Methods
Patients
Eighteen consecutive patients, some having been described previously,18 with unilateral symptomatic or asymptomatic high-grade stenosis or occlusion of the ICA or MCA were included. Inclusion criteria were the following: 1) Doppler sonography evidence of unilateral, >70% (NASCET criteria) high-grade extracranial ICA stenosis or ICA occlusion; or 2) Doppler sonography or MR angiographic evidence of a high-grade (>50%) unilateral intracranial ICA or proximal MCA stenosis/occlusion. Patients with relevant bilateral stenoses (contralateral stenosis of >50%) were not considered for inclusion in the study. The study was approved by the local institutional review board, and written informed consent was obtained from all patients before study enrollment.
MR Imaging Protocol
MR imaging measurements were performed on a 3T whole-body scanner (Magnetom Trio; Siemens, Erlangen, Germany) using the body coil of the scanner for radiofrequency transmission and an 8-channel phased array head coil for signal reception. The MR imaging examination included quantitative T2 mapping as well as DWI and PWI, MRA, and conventional T1WI and T2WI.
Anatomic imaging was based on a T1-weighted MPRAGE sequence with the following imaging parameters: TR/TE/TI = 2250/2.6/900 ms, FOV = 256 × 224 × 144 mm3, whole-brain coverage, isotropic spatial resolution = 1 mm, 22% oversampling in the slice-encoding (left-right) direction to avoid aliasing, bandwidth = 200 Hz/pixel, excitation angle = 9°, duration = 7 minutes 23 seconds. This dataset was used for tissue segmentation.
Quantitative T2 mapping was based on a fast spin-echo sequence with an echo-train length of 11 echoes per excitation, an echo spacing of 17.1 ms, and the following imaging parameters: 50 axial slices with 2-mm slice thickness, no interslice gap, TR = 10 seconds, bandwidth = 100 Hz/pixel, 180° refocusing pulses, matrix size = 192 × 132 (readout × phase encoding), FOV = 240 × 165 mm2, and in-plane resolution = 1.25 × 1.25 mm2. For quantitative T2-mapping, 5 datasets were acquired with different TE values (17, 86, 103, 120, 188 ms), keeping all other acquisition parameters constant. The total duration was 11 minutes 50 seconds.
PWI was based on a gradient-echo EPI sequence with the following imaging parameters: TE = 35 ms, TR = 1500 ms, flip angle = 90°, FOV = 192 × 192 mm, matrix = 64 × 64, slice thickness = 4 mm, number of slices = 16, voxel size =2.0 × 2.0 × 4.0 mm3, and acquisition time = 1 minute 15 seconds. The contrast agent (0.1 mmol/kg of Gd-DTPA, Magnevist; Bayer HealthCare Pharmaceuticals, Wayne, New Jersey) was injected into an antecubital vein using a power injector at a rate of 5 mL/s followed by a flush with 10 mL of saline.
Image Postprocessing and Analysis
Perfusion-weighted MR raw images were processed on a pixel-by-pixel basis to generate maps of the nondeconvolved TTP, rCBF, and rCBV. To determine the shape of the arterial input function, 2 experienced neuroradiologic readers manually selected 5–10 pixels in consensus over the proximal MCA in the unaffected hemisphere. For calculation of rCBF and rCBV maps, we used the model-independent (nonparametric) singular-value decomposition method described by Ostergaard et al.19
Further image postprocessing and analysis were performed with FMRIB Software Library tools (FSL, Version 5.0; https://fsl.fmrib.ox.ac.uk/fsl). All images were skull-stripped before coregistration. Quantitative T2 maps were linearly coregistered to the T1-weighted images. Because PWI parameter maps do not provide sufficient anatomic contrast for coregistration, we used the following: The first volume of the PWI time-series, which showed adequate anatomic contrast for coregistration, was extracted and coregistered to the T1-weighted image, saving the coregistration matrix. This matrix was applied to coregister TTP, rCBF, and rCBV parameter maps to the T1-weighted image. For coregistration of DWI maps, the first diffusion-weighted image (b=0, purely T2-weighted) was coregistered to the T1-weighted image. The coregistration matrix of this step was used for later coregistration of the diffusion-weighted image with b=1000 s/mm2 to the T1-weighted image.
T1-weighted images were segmented with the FMRIB Automated Segmentation Tool (FAST; https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FAST),20 yielding partial volume estimate maps for white matter, gray matter, and CSF. A lower partial volume estimate threshold of 0.95 was applied to both the WM and the GM maps.21 For segmentation of the cerebral cortex, voxels of noncortical structures, including subcortical lesions or (enlarged) perivascular spaces misclassified as GM, were removed from the GM maps. Both the WM and cortex maps were binarized to receive tissue masks for the qT2 and PWI parameter maps (Fig 1, upper row, and Fig 2A).
Maps of rCBF and rCBV were processed voxelwise to calculate maps of the rCBF/rCBV ratio as possible indicators of local CPP.16,22 The binary WM and cortex masks were applied to the qT2 and PWI parameter maps to generate WM and cortex maps for each parameter (Fig 1, lower row). After the masking procedure, macroscopically visible lesions (including DWI lesions) were automatically removed from the respective parameter maps (On-line Figure). The second T2-weighted raw image (TE = 86 ms) from each patient was thoroughly inspected and manually corrected if necessary to exclude residual visible lesions or to avoid removing parenchyma beyond the edges of the lesions.
3D ROIs were defined on the basis of areas with TTP delay in the affected hemisphere. For WM, an elliptic periventricular ROI was placed in the unaffected hemisphere. The mean TTP plus 2 SDs was defined as the upper limit of normoperfusion and was used as a lower threshold for TTP maps. On the basis of this threshold, different degrees of TTP delay were determined by gradually increasing the threshold in intervals of 2 seconds (eg, 0–2 and 2–4 seconds, up to the maximum TTP delay) (Fig 1, upper row). For the cortical ribbon, the delineation of areas with different degrees of TTP delay was not deemed feasible. Therefore, only cortical areas were included that depended, with an extremely high likelihood (considering all available imaging data including CTA/MRA and DSA to assess a potentially adjusted watershed area, the presence of communicating arteries, and so forth), on the stenosed/occluded artery and should have been affected by resulting perfusion abnormalities (Fig 2). Parameter values from WM areas with different ranges of TTP delay and from cortical areas were compared with values from corresponding contralateral unaffected areas by mirroring the ROIs to the contralateral side. To determine relative changes within hypoperfused tissue, we calculated a hemispheric asymmetry index (AI) for each parameter using the following equation23,24: where Ri and Rc represent the ROI values in the ipsilateral and contralateral hemispheres, respectively. For a comparative morphometric measure, the volumes of affected WM and hypoperfused cortex (including visible lesions) were determined on the segmented T1-weighted images and compared with tissue volumes on the contralateral side. Before volumetric analysis, the tissue maps were linearly coregistered to standard space (Montreal Neurological Institute 152, 1 mm; http://www.bic.mni.mcgill.ca/ServicesAtlases/) using 6 df to minimize interpolation errors.
Statistical Analysis
Statistical analysis was performed using SPSS 22 (IBM, Armonk, New York). Because several parameters were not normally distributed (Kolmogorov-Smirnov test), we only used nonparametric statistical testing. Parameter values between corresponding ROIs (ipsilateral/contralateral hemispheres) were compared with the 2-sided Wilcoxon signed rank test. For evaluation of a parameter correlation, the 2-sided Spearman signed rank test was used. A P < .05 was considered statistically significant. Correction for multiple statistical testing was done by performing false discovery rate correction for all parameters. Data are given as mean ± SD unless indicated otherwise.
Results
The mean age of the patients (15 men, 3 women) was 57.7 ± 14.8 years (range, 28–83 years). Sites of stenosis and occlusion were the extracranial ICA (n = 8), the proximal MCA (M1 segment) (n = 7), and the intracranial ICA (n = 3). The right hemisphere was affected in 9 patients, while the other 9 patients had a left-sided pathology. Seven patients were recently symptomatic (stroke or TIA within the last 30 days before MR imaging), 6 had been formerly symptomatic (several years before MR imaging examination), and 5 patients had a completely asymptomatic vascular pathology. Demographic and clinical baseline characteristics for all patients are summarized in the Table. Four patients showed small (<10 cm3) ischemic lesions on DWI, which were removed from the parameter maps before further analysis as described above.
qT2 and Perfusion Parameters within Perfusion-Restricted Normal-Appearing WM
In normal-appearing WM with any TTP delay (>0 seconds, entire TTP lesion) ipsilateral to the vessel pathology, T2 values were significantly increased compared with corresponding contralateral normoperfused areas (110.08 ± 19.27 ms versus 106.39 ± 15.71 ms, AI = 3.07% ± 4.1%, P = .008). Perfusion parameters also showed significant changes within the entire TTP-delayed area with significant decreases of rCBF (AI = −15.57% ± 12.38%, P = .001) and an increase of rCBV (AI = 5.8% ± 9.68%, P = .041), resulting in a significant decrease of the rCBF/rCBV ratio (AI = −22.11% ± 11.94%, P < .001) compared with contralateral unaffected areas. These changes were coherent for the respective parameters at different degrees of perfusion delay, except for a low TTP delay of 0–2 seconds and a severe TTP delay of 6–8 seconds (On-line Table). Asymmetry indices increased for qT2 and decreased for the rCBF/rCBV ratio with the increasing TTP delay (On-line Table). AI for qT2 showed a strong positive correlation with the TTP delay (r = 0.595, P < .001) and a strong negative correlation with the AI for rCBF/rCBV (r = −0.352, P = .002), with the latter showing a negative correlation with the TTP delay (r = −0.570, P < .001) (Fig 3A–C).
qT2 and Perfusion Parameters within TTP-Delayed Cortical Areas
Similarly, in the TTP-delayed cortex, qT2 was significantly elevated in comparison with the corresponding contralateral cortical areas of the unaffected hemisphere (148.78 ± 37.46 ms versus 142.58 ± 31.18 ms, AI = 3.61% ± 5.02%, P = .007). Analogous to the WM ROIs, significant decreases of rCBF (AI = −7.35% ± 12.88%, P = .022) and significant increases of rCBV (AI = 11.24% ± 45.53%, P = .013) were detected, resulting in significant decreases of the rCBF/rCBV ratio (AI = −15.17% ± 13.41%, P < .001), compared with unaffected cortical areas (On-line Table). In contrast to WM, relative changes of T2 values in hypoperfused cortical areas did not correlate significantly with changes of perfusion parameters (AI qT2 versus TTP delay: r = 0.311, P = .21; AI qT2 versus AI rCBF/rCBV: r = −0.292, P = .256), while again a strong significant negative signed correlation was found between the rCBF/rCBV AI and TTP delay (r = −0.806, P < .001) (Fig 3D).
Volumetric Analysis
Volumes of affected WM and hypoperfused cortex did not differ significantly from volumes of the contralateral unaffected tissue (269.33 ± 25.5 cm3 versus 270.81 ± 27.29 cm3, P = .372 and 79.54 ± 20.5 cm3 versus 79.12 ± 20.16 cm3, P = .948). False discovery rate correction was performed for the respective parameters, and the corrected level of significance was P = .033.
Discussion
This study revealed significantly increased qT2 values in perfusion-delayed normal-appearing WM and in cortical areas in patients with unilateral high-grade carotid occlusive disease compared with corresponding contralateral areas. In addition, these changes correlated with increasing hemodynamic failure as indicated by TTP delay, suggesting silent hypoperfusion-related microstructural tissue damage. Furthermore, the rCBF/rCBV ratio seems to be a promising parameter to assess CPP in chronic cerebral hypoperfusion.
In patients with high-grade carotid stenosis—even asymptomatic—qT2 showed a significant increase within normal-appearing GM and WM. This is in line with the literature findings, in which experimental studies revealed structural damage, including demyelination and axonal loss, in chronically hypoperfused brain parenchyma—in fact, the extent of structural damage correlated with the degree of perfusion impairment.8 In addition, glial activation has been demonstrated in experimentally induced chronic cerebral hypoperfusion.8 Gliosis in chronically hypoperfused tissue, which can be expected to cause prolongation of qT2 due to increased water content, may be caused by subtotal ischemic necrosis not leading to focal softening or cavitation (so-called “selective neuronal loss”).25,26 Hence, selective neuronal death may be one of the underlying causes of the significant prolongation of qT2 in normal-appearing hypoperfused tissue in this study. Furthermore, the qT2 increase in normal-appearing depending WM correlated significantly with the perfusion abnormalities depicted as perfusion delay on TTP maps, suggesting that the extent of tissue damage is related to the degree of perfusion delay. Most interesting, qT2 prolongation within hypoperfused tissue also showed a strong negative correlation with changes of the rCBF/rCBV ratio, which, for its part, correlated significantly with TTP delay (Fig 3A–C). By means of PET imaging, the CBF/CBV ratio has been validated as an index of local CPP,16 which itself is a function of mean arterial blood pressure. Because obtaining rCBF and rCBV via DSC MR imaging (a bolus-tracking technique) is subject to various sources of imprecision and is technically fundamentally different from measuring perfusion with PET, MR imaging–derived rCBF/rCBV might not necessarily have the same significance. However, because TTP has been shown to correlate with PET-derived CBF/CBV,17 the finding of a significant relationship between the rCBF/rCBV ratio and TTP delay in this study supports the hypothesis that rCBF/rCBV obtained with DSC MR imaging could also reflect the local CPP. The significant correlation between qT2 and TTP as well as the rCBF/rCBV ratio suggests that pathologic alterations of hypoperfused WM are mainly related to the degree of local CPP compromise, which determines the restriction of microvascular perfusion and oxygen supply and promotes microstructural tissue damage, even in normal-appearing tissue.
In cortical areas, the same changes of qT2 and perfusion parameters as in WM were detected, again indicating reduced perfusion pressure (On-line Table). Whereas a strong negative correlation between TTP delay and rCBF/rCBV was also found in the hypoperfused cortex, no significant correlation of qT2 increase with perfusion measures was detected. This difference between the cortex and WM might be explained by the anatomic structure of the respective supplying arterial vasculature. Because the perforating arteries supplying the deep WM are hardly collateralized, the extent of ischemic tissue damage depends on the severity of hemodynamic compromise.27,28 In contrast, a variable degree of leptomeningeal collateralization compensating for decreased perfusion can be expected in chronically hypoperfused cortex.29 This might explain why the magnitudes of both TTP delay and rCBF/rCBV reduction within affected cortical areas are lower than in WM (On-line Table and Fig 3). On the other hand, the occurrence of embolism or microembolism, which is not necessarily directly related to perfusion abnormalities, might be more predominant as a mechanism of tissue damage in the cortex.30
That volumetric quantification of normal-appearing WM and cortex did not reveal significant volume reduction of hypoperfused tissue excludes relevant partial volume effects from CSF, with a resulting increase of tissue water as the reason for the increase of T2 values observed in this study. The different underlying causes of microstructural changes in hypoperfused tissue, which may include a variety of pathologies, cannot be differentiated with qT2 mapping. Because microstructural changes in carotid occlusive disease have been shown to regress partly after stent placement6 or carotid endarterectomy,11,31 the underlying mechanisms might also comprise (partially) reversible conditions, though the pathophysiologic correlate for this finding remains currently unclear.
In summary, the results of this study suggest the presence of microstructural alterations that are not assessable using conventional MR imaging within chronically hypoperfused tissue in patients with carotid occlusive disease and exhausted autoregulation. In addition, our results might support the hypothesis that the DSC-derived rCBF/rCBV ratio could be used as an index of local CPP, comparable with TTP delay, as an appropriate indicator of hemodynamic failure in chronic cerebrovascular disease.
Limitations
This study has several limitations. First, perfusion values determined with DSC techniques only provide representative relative values, which must be compared with values from reference regions. Second, because this is not a longitudinal study, the temporal evolution of qT2 values and perfusion parameters and their possible response to reperfusion procedures was not evaluated. Third, although prolongation of qT2 has been described in chronic cerebral hypoperfusion and would be consistent with microstructural damage manifesting itself as gliosis and neuronal loss, a formal histologic validation of qT2 as a marker of microstructural tissue damage due to hypoperfusion is lacking. This should be pursued and further investigated in future studies. Additionally, the relationship between pathologic qT2 increases and clinical correlates such as cognitive performance must be investigated for the relevance in clinical practice. Finally, an investigation of qT2 changes in larger patient collectives with well-defined risk profiles and different statuses of reperfusion as well as a PET-based validation of MR imaging–based rCBF/rCBV as an index of CPP would be of interest.
Conclusions
Quantitative T2 mapping seems to detect microstructural damage in chronically hypoperfused but normal-appearing cortex and WM, which correlates with the degree of perfusion pressure impairment due to carotid occlusive disease. The rCBF/rCBV ratio and the nondeconvolved TTP determined with DSC MR imaging seem to be suitable for detecting reduced CPP and hemodynamic failure in these patients. After further research to elucidate the clinical correlates of these parameters and their response to reperfusion, these techniques could help to identify patients at risk of progressive cognitive decline who would potentially benefit from revascularization.
Footnotes
Disclosures: Alexander Seiler—UNRELATED: Grants/Grants Pending: Goethe University Frankfurt, Faculty of Medicine.* Ralf Deichmann—RELATED: Grant: support to the Brain Imaging Center Frankfurt, Comments: Bundesministerium für Bildung und Forschung (Germany) DLR 01GO0203, Deutsche Forschungsgemeinschaft (Germany) ZA 233/1–1*; UNRELATED: Consultancy: Wellcome Trust Centre for Neuroimaging at UCL, London, UK, Comments: consultant for MR imaging scanner procurement; Employment: University of Frankfurt; Grants/Grants Pending: Land of Hesse, Germany: LOEWE research grant “Neuronale Koordination Forschungsschwerpunkt Frankfurt,” 2011–2013 Deutsche Forschungsgemeinschaft: “Imaging of damage and repair in MS patients—Impact of new lesions on network activity, Teilprojekt B5” 2012–2016, Comments: public funding*; Royalties: published book chapters: 1) R. Deichmann, “Principles of MRI and Functional MRI.” In: fMRI Techniques and Protocols. M. Filippi, ed. Humana Press, Springer, Heidelberg, 2009; 2) R. Deichmann, U. Noeth, N. Weiskopf, “The Basics of Functional Magnetic Resonance Imaging.” In: EEG-fMRI. C. Mulert, L. Lemieux, eds. Springer, Heidelberg, 2010 (so far, no royalties paid); Travel/Accommodations/Meeting Expenses Unrelated to Activities Listed, Comments: occasionally travel expenses for invited lectures given at nonprofit organizations. Ulrike Nöth—RELATED: Grant: support to Brain Imaging Center Frankfurt, Comments: Bundesministerium für Bildung und Forschung (Germany) DLR 01GO0203, Deutsche Forschungsgemeinschaft (Germany) ZA 233/1–1.* Waltraud Pfeilschifter—UNRELATED: Consultancy: Sanofi-Aventis; Grants/Grants Pending: Stryker Neurovascular, Boehringer Ingelheim, Medtronic*; Payment for Lectures Including Service on Speakers Bureaus: Stryker Neurovascular, Boehringer Ingelheim. *Money paid to the institution.
Ralf Deichmann and Ulrike Nöth were supported by the Brain Imaging Center Bundesministerium für Bildung und Forschung (Germany) DLR 01GO0203, Deutsche Forschungsgemeinschaft (Germany) ZA 233/1-1.
Indicates open access to non-subscribers at www.ajnr.org
References
- Received November 28, 2017.
- Accepted after revision March 15, 2018.
- © 2018 by American Journal of Neuroradiology