Quantitative Histogram Analysis on Intracranial Atherosclerotic Plaques: A High-Resolution Magnetic Resonance Imaging Study.

BACKGROUND AND
PURPOSE: Intracranial atherosclerosis is one of the main causes of stroke, and high-resolution magnetic resonance imaging provides useful imaging biomarkers related to the risk of ischemic events. This study aims to evaluate differences in histogram features between culprit and nonculprit intracranial atherosclerosis using high-resolution magnetic resonance imaging.
METHODS: Two hundred forty-seven patients with intracranial atherosclerosis who underwent high-resolution magnetic resonance imaging sequentially between January 2015 and December 2016 were recruited. Quantitative features, including stenosis, plaque burden, minimum luminal area, intraplaque hemorrhage, enhancement ratio, and dispersion of signal intensity (coefficient of variation), were analyzed based on T2-, T1-, and contrast-enhanced T1-weighted images. Step-wise regression analysis was used to identify key determinates differentiating culprit and nonculprit plaques and to calculate the odds ratios (ORs) with 95% CIs.
RESULTS: In total, 190 plaques were identified, of which 88 plaques (37 culprit and 51 nonculprit) were located in the middle cerebral artery and 102 (57 culprit and 45 nonculprit) in the basilar artery. Nearly 90% of culprit lesions had a degree of luminal stenosis of <70%. Multiple logistic regression analyses showed that intraplaque hemorrhage (OR, 16.294 [95% CI, 1.043-254.632]; P=0.047), minimum luminal area (OR, 1.468 [95% CI, 1.032-2.087]; P=0.033), and coefficient of variation (OR, 13.425 [95% CI, 3.987-45.204]; P<0.001) were 3 significant features in defining culprit plaques in middle cerebral artery. The enhancement ratio (OR, 9.476 [95% CI, 1.256-71.464]; P=0.029), intraplaque hemorrhage (OR, 2.847 [95% CI, 0.971-10.203]; P=0.046), and coefficient of variation (OR, 10.068 [95% CI, 2.820-21.343]; P<0.001) were significantly associated with plaque type in basilar artery. Coefficient of variation was a strong independent predictor in defining plaque type for both middle cerebral artery and basilar artery with sensitivity, specificity, and accuracy being 0.79, 0.80, and 0.80, respectively.
CONCLUSIONS: Features characterized by high-resolution magnetic resonance imaging provided complementary values over luminal stenosis in defined lesion type for intracranial atherosclerosis; the dispersion of signal intensity in histogram analysis was a particularly effective predictive parameter.

Ischemic stroke is one of the leading killers and the first cause of disability in the world.[1] Recently, intracranial atherosclerosis (ICAS) has been recognized as the most common cause of ischemic stroke worldwide, accounting for ≈9% to 33% of all ischemic strokes.[2-4] The incidence rate in the Chinese population is particularly high, where ICAS accounts for 30% to 50% of ischemic strokes.[5] It has been shown that culprit lesions usually have specific structural features, which cannot be captured by angiography methods such as magnetic resonance angiography or computed tomography angiography, such as a relatively large lipid-rich core with an inflamed thin fibrous cap.[6,7] Therefore, visualization of lesion morphological and compositional features is essential for the risk assessment of ICAS.

High-resolution magnetic resonance imaging (hrMRI) is capable of identifying lesion morphology and different atherosclerotic components[8] including vessel wall remodeling, intraplaque hemorrhage (IPH), and plaque enhancement if a contrast medium is provided.[9] These detailed lesion characteristics predict future cerebrovascular events.[10] hrMRI can provide a direct visualization of atherosclerosis not only in large arteries, for example, carotid artery, but also in small intracranial arteries, including middle cerebral artery (MCA) and basilar artery (BA).[11-13] Both MCA and BA have been shown to contain a large proportion of ICAS plaques, for instance in a study by Qiao et al,[13] in which 32.2% ICAS plaques were located in the MCA and 40% in the BA.

Previous studies have mostly focused on lesion morphological and compositional features. The distribution of signal intensity, which is defined as the number of pixels with signal intensity in a certain categorical range, might also provide clinically relevant information, particularly in relation to the brain and cancer.[14] Pilot studies have reported that signal heterogeneity in hrMRI has potential in differentiating stable and unstable lesions.[11,12] The distribution of signal intensity is often visualized by a histogram, a graphical representation of a frequency distribution. It has been reported that parameters derived from histograms on MR images are useful for the prediction of histological type and grade in head and neck malignancy, as well as differentiating prostate cancer from normal tissues.[15] However, the characteristics of signal histograms in ICAS have not been widely investigated. In this retrospective study, a histogram analysis was performed to identify potential quantitative biomarkers for ICAS as visualized by hrMRI, and associated complementary values over lesion morphological and compositional features were also quantified.

Methods

Data Availability

Anonymized data not published within this article can be provided upon request by the corresponding author.

Study Population

Data from 247 consecutive patients with suspected acute stroke/transient ischemic attack who underwent hrMRI between January 2015 and December 2016 in Changhai Hospital (Shanghai, China) was used for analysis. This study was approved by the Institutional Review Board of Changhai Hospital with all patients having provided written informed consent. The inclusion and exclusion criteria are described in the Data Supplement.

MRI Acquisition

The hrMRI was performed in one of two 3T whole body systems (GE Signa 3.0T HDxt, GE Healthcare, Waukesha; and Skyra Siemens Healthcare, Erlangen, Germany). The purpose and protocols of hrMRI scanning are described in the Data Supplement.

Image Analysis

All atherosclerotic plaque images were analyzed by 2 radiologists (one radiology resident with 7 years’ experience and one senior radiologist with >15 years’ experience in the diagnostics of neuroradiology) who were not involved in statistical analyses. If there was marked eccentric or focal wall-thickening with plaque burden >40%,[11] we considered a plaque was present. Each detected plaque was classified as culprit or nonculprit with the consideration of patient clinical presentations and findings in brain magnetic resonance imaging ([diffusion-weighted imaging (DWI) and fluid-attenuated inversion recovery (FLAIR)]). A culprit plaque was identified as a lesion arising on the ipsilateral side to a fresh infarction on the DWI images with accompanying clinical symptoms, while a plaque was considered to be nonculprit when it occurred in asymptomatic patients. If more than one plaque was present in the same vascular territory, the most stenotic lesion was chosen for analysis. Any disagreements were resolved by consensus, having further assessed other image features such as enhancement ratio. The imaging assessments are described in the Data Supplement.

The histogram of a region of interest is a function showing a frequency distribution of signal intensity in that area. Skewness, kurtosis, and coefficient of variation (CV) are commonly used to characterize the statistical information of a histogram. The skewness describes the skew in the shape of the distribution. A skewed distribution histogram is one that is asymmetrical in shape. The kurtosis describes the peak and/or flatness of the curve peak; a more acute peak has higher kurtosis, and a more flattened peak has lower kurtosis. The CV describes the dispersion of signal intensity.[16] These 3 parameters are defined below.

Where x is the signal intensity value of each pixel within the region of interest; and are the mean and standard deviation (SD) of x, respectively; n is the pixel number. The histogram features were analyzed by 2 experienced reviewers using Image J (National Institutes of Health, Bethesda, MD).[17] The same region of interest was used for all sequences for the same patient and was adjusted accordingly if misalignment was found between differently weighted images because of motion.

Statistical Analysis

Statistical analyses were performed using SPSS 24.0 (IBM). Normality testing was performed to assess the variable distribution. Univariate analysis (t-test or Mann-Whitney U test where appropriate for the comparison of continuous variables, and χ2 for the categorical variables) was used to assess the association of parameter with plaque type (culprit or nonculprit). Multivariate logistic regression analysis was subsequently performed, which included variables with P<0.10 in univariate tests. The odds ratios (ORs) with 95% CIs were calculated by a logistic regression model. The diagnostic performance was described using receiver operating characteristic curves and area under curve (AUC) values. The reproducibility of continuous variables was evaluated using the intraclass coefficient with a 2-way random-effects model, and the Kappa value was determined for the categorical variables.

Results

A total of 247 patients were recruited and data from 57 patients were excluded because of intracranial aneurysm (n=22), dissection (n=11), ≥30% stenosis in extracranial carotid arteries (n=10), Moya-Moya disease (n=7), vasculitis (n=4), and poor image quality (n=3). One hundred ninety patients with ICAS (94 culprit and 96 nonculprit) were included in the final analysis. Eighty-eight lesions were located at MCA (Figure 1A and 1B) and 102 at BA (Figure 1C and 1D). Detailed patient demographics are listed in Table 1.

Figure 1.

High-resolution magnetic resonance imaging (hrMRI) images and histogram showing atherosclerotic plaques in middle cerebral artery (MCA) and basilar artery (BA). A, A 64-year-old man with acute ischemic stroke on the left cerebrum. TOF-magnetic resonance angiography (MRA) demonstrated the stenosis located on the left MCA, and DWI showed a sheet acute infarcts. hrMRI images including T2W, T1W, and CE-T1W at the most stenotic site visualized the plaque and the quantitative histogram parameters based on T1W were shown on the right. B, A 59-year-old woman without any neuro symptom. TOF-MRA demonstrated the stenosis on the right MCA, and DWI showed no acute infarct. hrMRI images including T2W, T1W, and CE-T1W at the most stenosis site visualized the plaque and the quantitative parameters based on T1W were on the right. C, A 56-year-old man with acute ischemic stroke. TOF-MRA demonstrated the stenosis in BA, and DWI showed an acute infarction within the left brain stem. hrMRI images including T2W, T1W, and CE-T1W at the most stenotic site visualized the plaque and the quantitative histogram parameters based on T1W were shown on the right. D, A 62-year-old man without any neuro symptom. TOF-MRA demonstrated the stenosis in BA, and DWI showed no acute infarct with the posterior circulation territory. hrMRI images including T2W, T1W, and CE-T1W at the most stenosis site visualized the plaque and the quantitative parameters based on T1W were on the right (histogram plot: x-axis represents the count of the pixel and y-axis represents the intensity of the plaque).

Table 1.

Patient Demographics

hrMRI Characteristics of ICAS

The hrMRI characteristics of the intracranial atherosclerotic plaques are summarized in Table 2. It can be seen that the stenosis category was not associated with lesion type for either MCA or BA, as nearly 90% of culprit lesions had a degree of luminal stenosis of <70%. hrMRI identified that IPH (P=0.015), minimum luminal area (MLA; P=0.001), and plaque burden (PB; P=0.011) were associated with culprit plaques at MCA. IPH (P=0.009) and enhancement ratio in contrast enhanced T1-weighted (CE-T1W; P=0.006) were the only 2 significant features differentiating culprit and nonculprit plaques at BA.

Table 2.

Comparison of hrMRI Characteristics Between Culprit and Nonculprit Plaques

When data from MCA and BA was pooled, it showed that 3 hrMRI characteristics, IPH (P<0.001), MLA (P<0.001), and PB (P=0.059), were associated with plaque type. Patient demographics including hypertension, diabetes mellitus, and hyperlipidemia, stenosis category, and enhancement ratio in CE-T1W were not significantly different between culprit and nonculprit plaques (P>0.10).

Histogram Analysis for ICAS

Histogram characteristics of the intracranial atherosclerotic plaques are summarized in Table 3 with representative cases shown in Figure 1. Univariate analysis showed that SD (P<0.001), minimum value (P<0.001), maximum value (P=0.037), and CV (P<0.001) were significantly different between culprit and nonculprit lesions in the MCA on T1W images. However, no significant difference was observed in both T2W and CE-T1W images. Similarly, minimum value (P=0.003) and CV (P=0.001) were associated with plaque type in the BA on T1W images whereas no significant difference was found in features on T2W and CE-T1W images.

Table 3.

Histogram Features of Atherosclerotic Plaques on MCA and BA

When data from both MCA and BA were pooled, area (P=0.002), SD (P<0.001), minimum value (P<0.001), maximum value (P=0.014), and CV (P<0.001) were all significantly different between culprit and nonculprit lesions on T1W images and no significant difference was found in both T2W and CE-T1W images. No significant difference was found in term of either skewness or kurtosis.

Univariate and Multivariate Logistic Regression Analysis

Potential predictive parameters listed in Tables 2 and 3 with P<0.1 were included in the regression analyses (Table 4). The univariate logistic regression analysis shown that luminal stenosis, IPH, MLA, PB, histogram parameters, including SD, minimum value, maximum value and CV, were significantly different between culprit and nonculprit plaques in the MCA, and enhancement ratios in CE-T1, IPH, MLA, area, minimum value and CV were significantly different between culprit and nonculprit lesions in the BA. Multiple logistic regression analyses showed that IPH (OR, 16.294v[95% CI, 1.043–254.632]; P=0.047), MLA (OR, 1.468 [95% CI, 1.032–2.087]; P=0.033), and CV (OR, 13.425 [95% CI, 3.987–45.204]; P<0.001) were 3 significant features associated with plaque type in the MCA. The AUC values of IPH, MLA, and CV were 0.571, 0.712, and 0.800, respectively (Figure 2A). The combination of these 3 features improved the AUC to 0.861 with sensitivity and specificity being 0.83 and 0.80, respectively. Except for IPH, MLA- and CV-defined AUC were both significantly higher than the one based on luminal stenosis (AUC=0.609).

Table 4.

Logistic Regression of the Significant Characteristics on Intracranial Atherosclerotic Plaques

Figure 2.

Receiver operating characteristic curves to differentiate culprit and nonculprit plaques, and odds ratio (OR) and 95% CI on the basis of multivariable logistic regression for each parameter. A, For middle cerebral artery (MCA); (B) for basilar artery (BA); (C) pooled data from both MCA and BA; and (D) OR and 95% CI on the basis of multivariable logistic regression.

The enhancement ratio in CE-T1W (OR, 9.476 [95% CI, 1.256–71.464]; P=0.029), IPH (OR, 2.847 [95% CI, 0.971–10.203]; P=0.046), and CV (OR, 10.068 [95% CI, 2.820–21.343]; P<0.001) were significantly associated with the plaque type for lesions in BA. Their AUC of receiver operating characteristic curves were 0.643, 0.605, and 0.708, respectively (Figure 2B). The combination of these three features improved the AUC to 0.804 with sensitivity and specificity being 0.79 and 0.77, respectively. In BA, compared with stenosis (AUC=0.546), IPH and enhancement did not improve AUC significantly, while CV did (P=0.033).

When data from both MCA and BA were pooled for the multiple logistic regression analysis (Figure 2C), the location where plaque was present (OR, 0.352 [95% CI, 0.134–0.925]; P=0.034), IPH (OR, 5.851 [95% CI, 1.624–21.085]; P=0.007), and CV (OR, 7.875; 95% CI, 2.754–22.515; P<0.001) were significant factors in differentiating plaque type with AUC of receiver operating characteristic curves being 0.569, 0.596, and 0.766, respectively (Figure 2D). When these 3 features were combined, the AUC was improved to 0.831 with sensitivity and specificity being 0.76 and 0.77, respectively. However, it was not significantly higher than the one of CV (P=0.362) whose sensitivity was 0.79, the specificity 0.80, and the accuracy 0.80 (Table I in the Data Supplement).

Reproducibility of Measurements

The reliability of hrMRI for the identification of morphological and compositional features in ICAS was assessed through the intra- and inter-observer agreement analysis, as listed in Table II in the Data Supplement.

Discussion

This study demonstrated that compared with the traditional assessment criterion (luminal stenosis), hrMRI-defined features provided a significant added value in differentiating culprit and nonculprit lesions in both the MCA and BA. Moreover, compared with conventional hrMRI-defined parameters, including IPH, MLA, and PB, histogram analysis further improved the accuracy for the lesion type differentiation. Histogram-defined signal dispersion, CV, was a strong predictor for the differentiation of lesion type for both MCA and BA. Compared with CV, the combination of other parameters, including IPH, MLA, and enhancement ratio did not improve the predictive power significantly.

Assessing the difference of hrMRI-defined plaque characteristics between culprit and nonculprit lesions is important, and it is also important to quantify the power of specific lesion characteristics in differentiating culprit and nonculprit lesions. Most previous studies were designed to identify higher risk hrMRI-defined features associated with patient symptoms. It has been suggested that compared with nonculprit plaques, culprit plaques are thicker with a larger wall area, volume, and greater PB and higher prevalence of expansive remodeling[3,18] and symptomatic patients tended to have a lesion with an irregular lumen surface.[18] Qiao et al[13] reported that although positive remodeling was marginally associated with downstream stroke presence, arterial remodeling may provide insight into stroke risk. Several articles have reported that the imaging intensity features such as IPH and plaque enhancement had a significant association with recent infarction in the territory of the culprit plaque in ICAS.[19] In our study, 9/37 (24.3%) of culprit plaques and 16/51 (31.4%) of nonculprit plaques showed enhancement in MCA (P=0.633), where enhancement is defined by a signal intensity increase of at least 15% compared with the precontrast image. In BA, 34/57 (59.6%) of culprit plaques and 15/45 (33.3%) of nonculprit plaques showed enhancement (P=0.01). The results obtained indicate that enhancement ratio was predictive for BA, which agreed with a previous study which found that enhancement was independently associated with culprit plaques in BA,[12] whereas the enhancement ratio was not significantly different for lesions in MCA. A previous study also pointed out that more culprit lesions showed enhancement than nonculprit ones in MCA. However, adding enhancement did not significantly improve the diagnostic accuracy.[11]

Histogram analysis has previously been performed with coronary computed tomography angiology and carotid B-mode sonography.[17,20,21] Obtained results demonstrated that histogram-defined parameters could discriminate symptomatic from asymptomatic plaques.[22,23] This study has demonstrated its capacity in differentiating lesion type for ICAS. As shown in Table 3, CV was capable of differentiating culprit and nonculprit lesions only in T1W images, not in T2W or CE-T1W images. Compared with nonculprit lesions, culprit lesions had a larger CV as a consequence of a larger difference between the minimum and maximum intensity value and/or smaller mean intensity value in T1W images. Although the exact mechanism is unclear, it may be related to atherosclerotic inclusions in the lesion. According to in vivo studies with atherosclerotic plaques from the carotid artery[24] and ex vivo studies with lesions from MCA,[6] hrMRI provides contrasts which could differentiate various atherosclerotic components in the lesion. For instance, hyperintensity in T1W images may be associated with high risk recent IPH or lipid core.[25] Because of their deep location and small size, it is hard to segment the components of intracranial atherosclerotic plaque. However, it could be rational to hypothesize that the wider signal intensity dispersion in T1W images indicates more types of components that might have been well known, for example, lipid, IPH, calcium and fibrous tissue, or which may not have been widely reported previously. For example, a recent study showed that the main chemical structure of advanced glycation end-products was a better candidate molecular imaging probe to identify vulnerable plaques on T1W than T2W images.[26]

Results obtained from this study suggest that TOF and T1W sequences might be sufficient in characterising the risk of intracranial atherosclerotic plaques based on the following facts: (1) both stenosis and MLA can be quantified by either TOF or T1W image; (2) the combination of TOF and T1W sequence can reliably identify IPH[27]; and (3) histogram-defined CV on T1W image can differentiate lesion type with high sensitivity and specificity. Although enhancement ratio is predictive for BA, excluding this factor only reduced AUC to 0.720. The imaging protocol with only TOF and T1W sequence is therefore beneficial for patients who are contraindicated with contrast agents. However, future follow-up studies are necessary to validate the predictive power of this scanning protocol in predicting subsequent ischemic events.

This study confirmed that IPH was associated with lesion type in both MCA and BA (Tables 2 and 4). Xu et al[19] reported that hrMRI-defined IPH was observed more often in culprit MCA plaques than in nonculprit plaques. A recent study showed that IPH was the strongest independent marker of symptomatic status in patents with atherosclerotic plaque in BA.[28] IPH is mostly attributed to fragile and leaky neovascularities with endothelial disruption and large local deformation.[29] It is one of the key factors leading to plaque progression[30] and therefore clinical symptoms.[31] It is clear that IPH has a long-term effect on plaque progression. However, the duration until the clinical event remains unclear as patients with lesions with IPH may remain asymptomatic for as long as 54 months.[31] Early identification of IPH might prove invaluable in optimizing management to minimize the risk of future adverse events, but future studies are needed to decide when aggressive treatments are truly required.

This study revealed that lesions in MCA and BA had common and distinct features (Table 4). MLA was predictive for MCA, but not for BA, and enhancement ratio was predictive for BA, but not for MCA. Multivariate regression analysis with data from both MCA and BA showed that location (MCA or BA) was an independent factor implying the difference between lesions from these 2 locations (Figure 2C). Such location-dependent features have been reported previously—compared with anterior circulation plaques, posterior circulation plaques had a larger plaque burden, higher remodeling ratio, and more often exhibited positive remodeling.[13] Postmortem studies showed that plaques in BA tended to be concentric as opposed to eccentric as in the anterior circulation with higher incidence of IPH and tended to have a higher proportion of internal elastic lamina and less elastin.[32] Such location-dependent features might be because of the differences in local arterial configuration and therefore the local hemodynamic environment.

Despite interesting findings, there are limitations to this study. First, this study is limited to Chinese populations. Second, our analysis was performed using 2-dimensional (2D) imaging data—its application is limited by low spatial resolution in the slice-select direction, long required acquisition times, and the difficulty of positioning 2D slices in one scan to capture multiple intracranial vessels with varying orientations. It is possible that the use of 3D imaging data could better characterize plaque features for improved histogram analysis. Moreover, limited by the imaging resolution, only lesions in large arteries, for example, MCA and BA, are visible as results most of studies are focused on large intracranial arteries. However, lesions in smaller arteries which cannot be pictured by current hrMRI techniques might also cause infarction. The development on this aspect is urgently needed. Finally, the results obtained from the retrospective analysis in this study need to be further confirmed by prospective studies including different ethnic populations to determine whether this model can be used for diagnosis of the vulnerability of atherosclerotic plaques by TOF and T1W images alone.

Acknowledgments

Drs Shi and Li contributed to the study concept and design, analysis and interpretation of data, drafting/revising the manuscript for content, and statistical analysis. Dr Zhao was involved in analysis and interpretation of data. Dr Peng contributed to the study concept and design and manuscript revision. Z. Meddings participated in manuscript revision. Drs Jiang and Liu performed analysis and interpretation of data. Dr Teng took part in the study concept and design, analysis, and interpretation of data. Dr Lu contributed to the study concept and design and revising the manuscript for content.

Sources of Funding

This study was supported in part by the National Natural Science Foundation of China (No. 81670396, No. 81270413, and No. 31470910), the Youth Foundation of Changhai Hospital (2018QNB009) and British Heart Foundation (PG/18/14/33562).

Disclosures

Dr Zhongzhao Teng is the chief scientist in Tenoke Ltd., Cambridge, UK and Nanjing Jingsan Medical Science and Technology, Ltd., Nanjing, China. Other authors do not have any conflict of interest to declare.