| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
(Stroke. 2006;37:93.)
© 2006 American Heart Association, Inc.
Original Contributions |
From the Imaging Research Laboratories (S.E.C., V.B., B.K.R.) and Blackburn Cardiovascular Genetics Laboratory and Vascular Biology Research Group (R.A.H.), Robarts Research Institute; Department of Medical Biophysics (S.E.C., B.K.R.), University of Western Ontario; and Department of Clinical and Neurological Sciences (V.B., R.R.H.), London Health Sciences Centre, London, Ontario, Canada.
Correspondence to Brian Rutt, Imaging Research Laboratories, Robarts Research Institute, PO Box 5015, 100 Perth Dr, London, Ontario, Canada. E-mail brian.rutt{at}imaging.robarts.ca
| Abstract |
|---|
|
|
|---|
Methods Forty-eight carotid endarterectomy specimen cross-sections had AHA lesion grade determined histopathologically and were concurrently imaged using combinations of 8 MRI contrast weightings in vitro. A maximum likelihood classification algorithm generated MRI "maps" of plaque components, and an AHA lesion grade was assigned correspondingly. Additional analyses compared classification accuracy obtained with a commonly used set of magnetic resonance contrast weightings [proton density (PDw), T1 (T1w), and partial T2 (T2w)] to accuracy obtained with the combination of PDw, T1w, and diffusion-weighted (Dw) contrast.
Results For the 8-contrast combination, the sensitivities for fibrous tissue, necrotic core, calcification, and hemorrhage detection were 83%, 67%, 86%, and 77%, respectively. The corresponding specificities were 81%, 78%, 99%, and 97%. Good agreement (79%) between magnetic resonance and histopathology for AHA classification was achieved. For the PDw, T1w, and Dw combination, the overall classification accuracy was insignificantly different at 78%, whereas the overall classification accuracy using PDw, T1w, and partial T2w contrast weightings was significantly lower at 67%.
Conclusions This study provides proof-of-principle that the composition of atherosclerotic plaques determined by automated classification of high-resolution ex vivo MRI accurately reflects lesion composition defined by histopathological examination.
Key Words: atherosclerosis carotid arteries MRI
| Introduction |
|---|
|
|
|---|
Current imaging studies of atherosclerotic lesions often rely on a human observers interpretation of magnetic resonance (MR) images, making results subject to considerable variability.5 Automated classification, on the other hand, could yield objective and reproducible assessment of atherosclerotic plaque composition based on MR images and therefore, lead to an increased sensitivity to disease progression. Promising early work done in this field was limited to a degree, because either the classifier required a priori knowledge of the number of plaque constituents present or was validated with a small data set.8,9
We have developed a quantitative and reproducible method for extracting relevant information about plaque composition using ex vivo MRI in conjunction with a maximum likelihood classification algorithm. It exhibited a high degree of accuracy when applied to a broad range of plaque types; this validation was accomplished using a large test data set, independent of the training data. In addition, we compared the accuracy of plaque compositional analysis based on the combination of proton-densityweighted (PDw), T1-weighted (T1w), and diffusion-weighted (Dw) images to a commonly used set of contrast weightings: PDw, T1w, and partial T2 weighted (T2w).
| Methods |
|---|
|
|
|---|
After imaging, histology was obtained; well-matched axial MR, µCT, and histological images were selected; and the histological and µCT images were registered to the corresponding MR images. The MR images were classified on a pixel-by-pixel basis into 4 tissue classes: fibrous/loose connective tissue, hemorrhage, necrosis, and calcification, using a maximum likelihood classification algorithm. The accuracy of the classified images was determined by pixel-by-pixel comparison to a pathologists interpretation of the histopathology. AHA grading was derived both conventionally from the histological images and from the classified MR images, and these gradings were compared.
MR and µCT Imaging
MRI and µCT imaging of the endarterectomy specimens were performed using a protocol described previously.8 Image acquisition parameters were as follows: proton density weighted with repetition time (TR)/echo time (TE) 2000/7.8 ms; partial T2w with TR/TE 2000/32 ms; T2w with TR/TE 2000/62 ms; T1w with TR/TE 500/7.8 ms fast spin-echo; fast imaging using steady-state acquisition with TR/TE/flip 5.3/2.5/35; T1w spoiled gradient echo with TR/TE/flip 35/5.9/25; T1w spoiled gradient echo with magnetization transfer, TR/TE/flip 35/5.9/25; and Dw (b=1500 s/mm2) spin-echo with TR/TE 2000/25 ms. The in-plane resolution for all of the sequences was 156x156 µm, and the slice thickness ranged from 200 to 1000 µm.
After MRI, the endarterectomy samples were fixed in formalin and imaged with µCT.8 Because the specimens were decalcified before histological sectioning, the µCT images served as the gold standard for identification of calcified plaque regions.
After imaging, 4-µmthick sections were taken at 0.5-mm intervals along the length of the plaque and stained with Movats pentachrome and hematoxylin/eosin. The histological sections provided the reference standard for the presence of various plaque components other than calcification.
MR Image Processing and Classification
From the 13 carotid endarterectomy specimens, 48 axial cross-sections were selected for analysis, with a range of 3 to 5 slices per specimen. Although there were >48 cross-sections available with MR, µCT, and histological images acquired, these 48 were selected as the minimal set for which all 3 of the image types were available, artifact-free, registerable to one another, and subject to the requirement that within 1 specimen, selected cross-sections were not too similar (ie, different total plaque cross-sectional area, arrangement, and/or number of plaque components). To avoid selection bias, slices were selected before automated classification and pathologist interpretation. Selection based on these criteria maximized the range of plaque and tissue types included in the study while minimizing errors in classification because of registration inaccuracies.
The MR images of the endarterectomy specimens were classified using the maximum likelihood classifier, a commonly used supervised classification algorithm.10 Small regions of interest representing fibrous/loose connective tissue, lipid-rich necrosis, calcification, and hemorrhage were selected from the 12 endarterectomy slices in the training data set. These training regions constituted only 2.5% of the total plaque cross-sectional area and were selected by an expert pathologist to enclose regions of unambiguously defined, uniform tissue on the basis of examination of the histological images. These small training regions were selected from 12 of the 48 cross-sections, spanning and uniformly distributed across 8 of the 13 endarterectomy specimens (the "training set"). The fact that such a small proportion of the total plaque pixels was used for training meant that the entire set of 48 cross-sections (the "test set") could be used to independently validate the classifier (97.5% of the analyzed pixels were different than the training pixels).
Based on the training data, the signal intensity distribution for each plaque component was modeled as a multivariate Gaussian distribution. Initial analysis showed that the pixel percentage and training region distribution we used produced a robust estimate of the underlying intensity distributions (mean intensity and SD, both used for pixel classification with this method). For every pixel in the test set, the maximum likelihood classification algorithm reports the tissue class with the highest probability given the overlapping distributions modeled by the training data. After the classifier was trained, all 48 of the cross-sectional plaques were fully classified.
Validation of Classified MR Images by Comparison With Histopathology
The accuracy of the classifier was determined by comparing the classified image on a pixel-by-pixel basis with µCT and histology, which served as the gold standards for calcification and other plaque components, respectively. The µCT images were registered to the MR images using linear registration (scaling, rotation, and translation only), with the ImageJ plug-in Align3TP software (National Institutes of Health), whereas the digitized histological images were registered to the MR images using the Delaunay triangle nonlinear registration method (ENVI, Research Systems Incorporated) so that shearing and other distortions caused by histological processing could be corrected. Applying a threshold to the µCT images generated "truth" regions of interest for calcification. A pathologists interpretation of the histology served as the truth regions for other plaque components. Using an electronic tablet, the pathologist, blinded to both the classifier results and the appearance of the MR images, manually traced appropriate regions of interest onto the digitized histological images, classifying all of the pixels within the plaque into 1 of the 4 defined tissue classes.
Registration of histology to the MR images permitted the 2 images to be overlaid for direct comparison of the MR-derived classification of each individual pixel to the truth class for that same pixel. Thus, the registration of histology to MR allowed for the sensitivity, specificity, and overall accuracy to be calculated for each plaque component on a pixel-by-pixel basis.8
Evaluation of Sensitivity, Specificity, and Accuracy of MR Classification
The sensitivity, specificity, and overall accuracy for all of the plaque components were determined for the following 3 combinations of MR contrast used in the classification process: (1) all 8 contrast weightings; (2) only PDw, T1w, and Dw contrast weightings; and (3) only PDw, T1w, and partial T2w contrast weightings. This set has often been used in other studies.11,12
The overall accuracies for each of the 3 conditions were compared with a paired t test. For this purpose, data from cross-sections extracted from the same endarterectomy specimen were combined to ensure that the t statistic was calculated using truly independent samples (n=13 endarterectomy specimens), a requirement for the t test.13
Modified AHA Grading System for Use With Classified MR Images
The AHA scientific statement described the characteristic components and pathogenic mechanisms of early7 and advanced6 atherosclerotic lesions. Because MRI is not able to resolve tissue structure at the histological level, the AHA classification criteria were modified for use with noinvasive image data. This was accomplished by combining type I and II lesions, as well as type IV and Va lesions, and using additional criteria as proposed previously.11,14 Table 1 details the specific criteria used to define this modified AHA grading system, which were applied to the MR-generated classification maps. For example, for the classified MR images, a plaque was labeled type VI if either the labeled necrotic core (of size >35 pixels) communicated directly with the lumen at any point (with contact length of at least 2 pixels) or a contiguous region of >150 pixels was labeled hemorrhage.
|
The corresponding histological sections were also classified according to AHA guidelines.6,7 This grading was performed by an expert pathologist using the traced and fully classified histological images. The AHA grade assigned to the classified MR image was compared with the true AHA grade as defined by the pathologist.
| Results |
|---|
|
|
|---|
=0.941) between repeated classifications, separated by 6 weeks. For practical reasons, this study did not involve multiple pathologists, and, therefore, interrater variability was not assessed.
Accuracy of Maximum Likelihood Classification of MR Images
The sensitivities and specificities for each plaque component were calculated for the 3 different conditions (Table 2). A paired t test demonstrated that the overall accuracy obtained with the best set of 3 contrast weightings (PDw, T1w, and Dw) (78%) was not significantly different from that obtained using all 8 of the contrast weightings (78%; P=0.42). In addition, the classifier accuracies achieved with all 8 of the contrasts and with these best 3 contrast weightings were both significantly higher than the accuracy attained with a commonly used combination of 3 contrasts (PDw, T1w, and partial T2w; 67%; P=0.004 and 0.010, respectively).
|
Comparison of AHA Plaque Classification by MR and Histopathology
Our sample of endarterectomy specimens consisted primarily of advanced atherosclerotic lesions (type IV-Va and type VI). Table 3 summarizes the results of AHA grading derived from maximum likelihood classification based on an input of 8 MR contrast weightings and AHA grading based on pathological classification for all 48 of the endarterectomy cross-sections analyzed in this study. The overall agreement was 79% (
=0.711). Figures 1 and 2
illustrate the different categories of AHA classification.
|
|
|
Often, classified images of the same plaque with a different set of contrast weightings would be assigned the same AHA rating, although the pixel-by-pixel accuracies differed. For example, Figure 3 demonstrates a type IV-Va plaque as assigned by histological examination. Maximum likelihood classification with PDw, T1w, and Dw agreed well with the pathologists tracing of the histology (pixel-by-pixel accuracy of 78%) and also correctly graded the plaque as type IV-Va. Classification based on PDw, T1w, and partial T2w also successfully rated this plaque as type IV-Va. However, the spatial distribution and percentage area of the lipid core by MRI was quite different from the pathologists segmentation and yielded a pixel-by-pixel accuracy of only 45%. With these 3 MR contrast weightings, MRI underestimated the size of the lipid-rich necrotic core and, in general, placed the necrotic core farther from the lumen.
|
| Discussion |
|---|
|
|
|---|
In this study, statistical analysis confirmed that the standard 3-contrast technique (PDw, T1w, and partial T2w) does not perform as well as the combination of PDw, T1w, and Dw. This suggests that significant improvement in overall accuracy of plaque characterization may be achieved if practical in vivo diffusion imaging protocols can be developed for high-resolution plaque imaging. Previous studies have identified both T2w and Dw as useful contrasts for identifying a necrotic core.5,14,17 However, these studies did not directly compare T2w and Dw; therefore, the relative advantage of Dw over T2w has not been clear until now.
Although the current study was limited to ex vivo analysis, it would not have been feasible to attempt this in an in vivo setting. The use of endarterectomy specimens allowed for the careful study of human atherosclerotic lesions with MR sequences not commonly used in clinical imaging of atherosclerosis (eg, Dw imaging), elucidating the potential value in additional investigation of these sequences. For example, a recent ex vivo study demonstrated optimal contrast between thrombosis and vessel wall with Dw MRI.18 In vivo plaque imaging will naturally incur decreased resolution, motion artifacts, and greater variability in MR signal characteristics, all of which will complicate the implementation and interpretation of both manual and automated classification methods. Therefore, the use of ex vivo data is a critical first step in the establishment and validation of sophisticated MR pulse sequences and plaque characterization algorithms.
The use of ex vivo specimens also carries with it some fundamental limitations, and extrapolation of our results to the in vivo setting must be approached with caution. Because of the differences in image quality mentioned above, we have not accounted for the limited resolution and signal-to-noise ratio, motion artifact, and other imperfections intrinsic to in vivo MRI. In addition, high-quality Dw images are more difficult to achieve in vivo, especially in the vicinity of the carotid bifurcation where flow characteristics and motion are significant factors.
We previously reported a sensitivity of 83.9% for necrotic core using a minimum-distance-to-means classifier.8 Despite using the maximum likelihood classifier (a more sophisticated algorithm) in the current study, the sensitivity for necrotic core dropped to 67%. We believe this difference is because of the addition of hemorrhage as a tissue type. Hemorrhage and necrotic core have similar MR signal characteristics and, thus, were difficult to separate with accuracy.
Several articles in the field of MR characterization of atherosclerotic plaque have reported values for sensitivity and specificity higher than those described here.12,14,19 However, these studies analyzed atherosclerotic plaques in quadrants, and the simple presence or absence of the plaque components was then assessed with MRI. Such methodology does not evaluate the ability of MR to precisely determine the spatial location and size of plaque components. The pixel-by-pixel method does provide such spatial information and leads to a different and far more stringent definition of MR classification accuracy, with correspondingly lower values for sensitivity, specificity, and overall accuracy.
| Acknowledgments |
|---|
Received July 26, 2005; revision received October 21, 2005; accepted October 25, 2005.
| References |
|---|
|
|
|---|
2. Carr S, Farb A, Pearce WH, Virmani R, Yao JS Atherosclerotic plaque rupture in symptomatic carotid artery stenosis. J Vasc Surg. 1996; 23: 755766.[CrossRef][Medline] [Order article via Infotrieve]
3. Glagov S, Bassiouny HS, Giddens DP, Zarins CK. Pathobiology of plaque modeling and complication. Surg Clin North Am. 1995; 75: 545556.[Medline] [Order article via Infotrieve]
4. Fayad ZA, Fuster V. Clinical imaging of the high-risk or vulnerable atherosclerotic plaque. Circ Res. 2001; 89: 305316.
5. Yuan C, Mitsumori LM, Ferguson MS, Polissar NL, Echelard D, Ortiz G, Small R, Davies JW, Kerwin WS, Hatsukami TS. In vivo accuracy of multispectral magnetic resonance imaging for identifying lipid-rich necrotic cores and intraplaque hemorrhage in advanced human carotid plaques. Circulation. 2001; 104: 20512056.
6. Stary HC, Chandler AB, Dinsmore RE, Fuster V, Glagov S, Insull WJ, Rosenfeld ME, Schwartz CJ, Wagner WD, Wissler RW. A definition of advanced types of atherosclerotic lesions and a histological classification of atherosclerosis. A report from the Committee on Vascular Lesions of the Council on Arteriosclerosis, Am Heart Association. Circulation. 1995; 92: 13551374.
7. Stary HC, Chandler AB, Glagov S, Guyton JR, Insull WJ, Rosenfeld ME, Schaffer SA, Schwartz CJ, Wagner WD, Wissler RW. A definition of initial, fatty streak, and intermediate lesions of atherosclerosis. A report from the Committee on Vascular Lesions of the Council on Arteriosclerosis, American Heart Association. Circulation. 1994; 89: 24622478.
8. Clarke SE, Hammond RR, Mitchell JR, Rutt BK. Quantitative assessement of carotid plaque composition using multi-contrast MRI and registered histology. Magn Reson Med. 2003; 50: 11991208.[CrossRef][Medline] [Order article via Infotrieve]
9. Itskovich VV, Samber DD, Mani V, Aguinaldo JG, Fallon JT, Tang CY, Fuster V, Fayad ZA. Quantification of human atherosclerotic plaques using spatially enhanced cluster analysis of multicontrast-weighted magnetic resonance images. Magn Reson Med. 2004; 52: 515523.[CrossRef][Medline] [Order article via Infotrieve]
10. Richards JA. Remote Sensing Digital Image Analysis: An Introduction. Berlin: Springer-Verlag; 1993.
11. Fayad ZA, Nahar T, Fallon JT, Goldman M, Aguinaldo JG, Badimon JJ, Shinnar M, Chesebro JH, Fuster V. In vivo magnetic resonance evaluation of atherosclerotic plaques in the human thoracic aorta: a comparison with transesophageal echocardiography. Circulation. 2000; 101: 25032509.
12. Itskovich VV, Choudhury RP, Aguinaldo JG, Fallon JT, Omerhodzic S, Fisher EA, Fayad ZA. Characterization of aortic root atherosclerosis in ApoE knockout mice: high-resolution in vivo and ex vivo MRM with histological correlation. Magn Reson Med. 2003; 49: 381385.[CrossRef][Medline] [Order article via Infotrieve]
13. Motulsky H. Intuitive Biostatistics. New York: Oxford University Press; 1995.
14. Cai JM, Hatsukami TS, Ferguson MS, Small R, Polissar NL, Yuan C. Classification of human carotid atherosclerotic lesions with in vivo multicontrast magnetic resonance imaging. Circulation. 2002; 106: 13681373.
15. Adams GJ, Greene J, Vick GW 3rd, Harrist R, Kimball KT, Karmonik C, Ballantyne CM, Insull W Jr, Morrisett JD. Tracking regression and progression of atherosclerosis in human carotid arteries using high-resolution magnetic resonance imaging. Magn Reson Imaging. 2004; 22: 12491258.[CrossRef][Medline] [Order article via Infotrieve]
16. Adame IM, van der Geest RJ, Wasserman BA, Mohamed MA, Reiber JH, Lelieveldt BP. Automatic segmentation and plaque characterization in atherosclerotic carotid artery MR images. MAGMA. 2004; 16: 227234.[CrossRef][Medline] [Order article via Infotrieve]
17. Toussaint JF, Southern JF, Fuster V, Kantor HL. Water diffusion properties of human atherosclerosis and thrombosis measured by pulse field gradient nuclear magnetic resonance. Arterioscler Thromb Vasc Biol. 1997; 17: 542546.
18. Viereck J, Ruberg FL, Qiao Y, Perez AS, Detwiller K, Johnstone M, Hamilton JA. MRI of atherothrombosis associated with plaque rupture. Arterioscler Thromb Vasc Biol. 2005; 25: 240245.
19. Shinnar M, Fallon JT, Wehrli S, Levin M, Dalmacy D, Fayad ZA, Badimon JJ, Harrington M, Harrington E, Fuster V. The diagnostic accuracy of ex vivo MRI for human atherosclerotic plaque characterization. Arterioscler Thromb Vasc Biol. 1999; 19: 27562761.
This article has been cited by other articles:
![]() |
C. Oppenheim, O. Naggara, E. Touze, J.-C. Lacour, E. Schmitt, F. Bonneville, S. Crozier, E. Guegan-Massardier, E. Gerardin, X. Leclerc, et al. High-Resolution MR Imaging of the Cervical Arterial Wall: What the Radiologist Needs to Know RadioGraphics, September 1, 2009; 29(5): 1413 - 1431. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Wintermark, S.S. Jawadi, J.H. Rapp, T. Tihan, E. Tong, D.V. Glidden, S. Abedin, S. Schaeffer, G. Acevedo-Bolton, B. Boudignon, et al. High-Resolution CT Imaging of Carotid Artery Atherosclerotic Plaques AJNR Am. J. Neuroradiol., May 1, 2008; 29(5): 875 - 882. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Qiao, I. Ronen, J. Viereck, F. L. Ruberg, and J. A. Hamilton Identification of Atherosclerotic Lipid Deposits by Diffusion-Weighted Imaging Arterioscler Thromb Vasc Biol, June 1, 2007; 27(6): 1440 - 1446. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Touze, J.-F. Toussaint, J. Coste, E. Schmitt, F. Bonneville, P. Vandermarcq, J.-Y. Gauvrit, F. Douvrin, J.-F. Meder, J.-L. Mas, et al. Reproducibility of High-Resolution MRI for the Identification and the Quantification of Carotid Atherosclerotic Plaque Components: Consequences for Prognosis Studies and Therapeutic Trials Stroke, June 1, 2007; 38(6): 1812 - 1819. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. M. Adame, P. J.H. de Koning, B. P.F. Lelieveldt, B. A. Wasserman, J. H.C. Reiber, and R. J. van der Geest An Integrated Automated Analysis Method for Quantifying Vessel Stenosis and Plaque Burden From Carotid MRI Images: Combined Postprocessing of MRA and Vessel Wall MR Stroke, August 1, 2006; 37(8): 2162 - 2164. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Stroke Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2006 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |