Abstract
Background Carotid atherosclerosis is a major etiology of stroke. Although intraplaque hemorrhage (IPH) is known to increase stroke risk and plaque burden, its long-term effects on plaque dynamics remain unclear.
Objectives This study aimed to evaluate the long-term impact of IPH on carotid plaque burden progression using deep learning-based segmentation on multi-contrast vessel wall imaging (VWI).
Methods Twenty-eight asymptomatic subjects with carotid atherosclerosis underwent an average of 4.7 ± 0.6 VWI scans over 5.8 ± 1.1 years. Deep learning pipelines segmented the carotid vessel walls and IPH. Bilateral plaque progression was analyzed using generalized estimating equations, and linear mixed-effects models evaluated long-term associations between IPH occurrence, IPH volume (%HV), and plaque burden (%WV) progression.
Results Two subjects with ipsilateral IPH developed new ischemic infarcts during follow- up. IPH was detected in 23/50 of arteries. Of arteries without IPH at baseline, 11/39 developed new IPH that persisted, while 5/11 arteries with baseline IPH exhibited it throughout the study. Bilateral plaque growth was significantly correlated (r = 0.54, p < 0.001), but this symmetry was weakened with IPH presence. The progression rate for arteries without IPH was -0.001 %/year (p = 0.895). However, IPH presence or development at any point was associated with a 2.34% absolute increase in %WV (p < 0.001), and %HV was associated with 0.73% per 2-fold increase over the mean of %HV (p = 0.005).
Conclusions IPH may persist asymptomatically for extended periods. While arteries without IPH demonstrated minimal progression under contemporary treatment, IPH significantly accelerated long-term plaque growth.
Introduction
Stroke remains the fifth leading cause of mortality in the United States, with an 8.4% increase in the age-adjusted death rate from 2011 to 20211. Carotid atherosclerosis is a leading etiology in stroke2. The progression of carotid atheroma, which leads to a reduction in vessel lumen and eventual plaque rupture, is closely related to the risk of cerebrovascular ischemic events3,4.
Understanding the mechanisms driving plaque evolution is essential for improving stroke prevention strategies.
A recent meta-analysis demonstrated the prevalence of intraplaque hemorrhage (IPH) in patients with both symptomatic and asymptomatic carotid stenosis and identified IPH as a stronger predictor of stroke than any known clinical risk factors5. Moreover, IPH is associated with increasing lumen stenosis and rapid vessel wall growth6–9. Despite this, little is known about the long-term behavior of IPH, including its occurrence, resolution, and distribution over extended periods, as most studies have been limited to follow-ups of one to two years7–9. Its long-term impact on plaque progression also remains inadequately explored6.
Traditionally, manual analysis of longitudinal vessel wall imaging MRI (VWI) has limited the ability to study carotid atherosclerosis due to its labor-intensive and expertise-dependent nature. Advances in deep learning-based segmentation now enabled automated, reproducible measurements of the carotid vessel wall10–12 and IPH, facilitating comprehensive studies of plaque dynamics over longer follow-up periods.
This study leverages these tools to analyze carotid atherosclerosis in an asymptomatic cohort over approximately six years. By investigating both systemic and localized trends in plaque evolution, this work provides novel insights into the natural history of carotid atherosclerosis, emphasizing the critical role of IPH in driving plaque burden progression and its implications for stroke prevention.
Methods
Study population
Participants for this study were recruited as part of a prospective cohort study over a ten-year period (2012-2022), aimed at examining the natural history of carotid atherosclerosis through serial VWI. Eligible participants included patients aged 35 years and older with clinically identified atherosclerosis in at least one carotid artery and without cerebrovascular events in the six months prior to enrollment. Study procedures and consent forms were reviewed and approved by the institutional review board, and all participants provided written informed consent prior to enrollment. Detailed information on cohort recruitment has previously been published9.
To monitor the evolution of atherosclerotic plaque, serial VWI scans with consistent parameters (detailed below) were conducted up to five times during the study period. For the analysis of the long-term plaque burden progression and the effects of IPH, we conducted a retrospective analysis of a subset of participants. Inclusion criteria for this retrospective analysis required subjects to have undergone at least three VWIs over a minimum span of five years. Arteries with a history of carotid endarterectomy (CEA) or stenting at study entry were excluded. Scans of arteries performed after CEA or stenting were also excluded.
MR Imaging Protocol
Two large-coverage three-dimensional isotropic VWI sequences were performed using a 3.0T Philips Ingenia CX scanner (Philips Healthcare, Best, The Netherlands): 1) Motion-sensitized driven equilibrium prepared rapid gradient echo sequence (MERGE13) was used to delineate vessel wall morphology; 2) an inversion-recovery gradient echo sequence with phase-sensitive reconstruction (simultaneous non-contrast angiography and intraplaque hemorrhage imaging, SNAP14) was used to detect and measure IPH. Detailed description of both sequences and imaging parameters can be found elsewhere9,15.
Ischemic infarcts assessment
Conventional brain MRI examinations, including diffusion weighted imaging (DWI), susceptibility weighted imaging (SWI), and T2-fluid attenuated inversion recovery (FLAIR), were performed at all VWI time points. A board-certified neuroradiologist directly compared the baseline and follow-up brain MRI images and scored for any incident infarcts (FLAIR and DWI images) or hemorrhages (SWI images). Vascular territory and signs of acuity or associated hemorrhage were noted for the infarcts.
Deep Learning-based Vessel Wall and IPH segmentation
Two deep learning-based imaging analysis pipelines were developed to bilaterally segment the carotid vessel wall and IPH on MERGE and SNAP sequences, respectively. Each scan was segmented independently.
For the segmentation of carotid lumen and outer wall boundary using MERGE16, the location of the left and right common carotid bifurcations was identified automatically. On each side, a 3D bounding box encompassing the bifurcation area as well as the distal and proximal segments of the bifurcation was extracted and then used for a preliminary 3D lumen segmentation. A trained radiologist reviewed the preliminary segmentation to assess image quality and ensure that the lumen from all slices encompassing the bifurcation region were correctly localized. Finally, 2D patches along both carotids were segmented slice-by-slice using a 2.5D convolutional neural network.
To segment IPH from SNAP, a nnUNet17-based algorithm was trained to identify hyperintense signals on the 2D axial reformatted images (with 2mm spacing) surrounding the bilateral bifurcations.
Detailed information about the training dataset and extensive validation studies can be found in the supplementary materials. The external validation datasets encompass histological validation of IPH detection and scan-rescan reproducibility assessments of both vessel wall and IPH segmentation.
For each MERGE image, after obtaining lumen and outer wall segmentations of the same artery, the range of slices with wall thickness greater than 2mm in any scan was included as atherosclerotic plaque. Plaque burden was defined as the percent wall volume (%WV), calculated as wall volume divided by total vessel volume, multiplied by 100.
For each SNAP image per artery per time point, IPH was categorized as present (IPH+) or absent (IPH-). IPH volume was aggregated as the sum of IPH areas per slice multiplied by 2mm, and percent hemorrhage volume (%HV) was defined as IPH volume divided by wall volume, multiplied by 100.
Statistical Analysis
Summary statistics of baseline demographics are presented as mean ± standard deviations (SDs) for continuous variables and counts with percentage for categorical variables.
For each analyzable plaque, a linear model was fitted across all available time points to determine the long-term growth rate of plaque burden18, defined as annualized absolute change in %WV, expressed in units of percent per year. The method for deriving the long-term plaque growth rate is illustrated in Supplementary Figure 1. Generalized estimating equation (GEE)-based linear regression models were used to evaluate associations of the long-term growth rate of each side with the contralateral side, adjusted for potential confounding factors of age, sex and statin use. Effect sizes were summarized using the partial correlation coefficient.
In all arteries, to investigate whether long-term progression rate in plaque burden depends on presence and volume of IPH, we employed a linear mixed-effects model. This model used %WV as the primary outcome variable and included random intercepts to account for within-subject and within-artery correlation. As fixed effects, the model included time as a continuous variable to estimate the growth rate of %WV and account for variation in scan intervals, a binary predictor indicating IPH presence, and %HV to assess the effect of IPH volume. %HV was log2- transformed to reduce right skewness, and further transformed such that it equals zero when IPH is absent, and equals the deviation from the population mean when IPH is present. The model setup was designed to disentangle the effects of the presence and volume of IPH on plaque burden progression. This model was further adjusted for age, sex and statin use.
All statistical analyses were performed using R 4.4.0. A two-sided p-value < 0.05 was considered statistically significant.
Results
Patient Characteristics
A total of 98 subjects were recruited prospectively, and population demographics can be found in the referenced publication9. Among these, 28 subjects underwent at least three VWIs over a span of at least five years and were therefore included in this retrospective analysis. All subjects were under contemporary medical management, and the baseline clinical information of these 28 subjects is presented in Table 1. After three scans were excluded due to poor image quality, each subject had an average of 4.7 ± 0.6 scans within 5.8 ± 1.1 years. Supplementary Figure 2 shows the study flow chart.
Among the 56 carotid arteries within the 28 subjects, six were excluded due to prior history of CEA or having fewer than three available scans due to CEA during the study period. The remaining 50 arteries from 28 subjects were included in artery-level analyses. Of these subjects, there were 22 with bilateral arteries available (44 arteries) for the subject-level analysis where plaque growth on each side was compared.
IPH Presence, Occurrences, and Relation to New Infarcts
Among the 50 carotid arteries from 28 subjects, 23 arteries (46%) from 17 subjects (60.7%) showed presence of IPH on VWI at any time point during the study period. Of the 39 IPH- arteries at the baseline scan, eleven (28.2%) became IPH+ on all subsequent scans after the initial occurrence of IPH, and one artery (2.56%) developed new IPH on the second scan but reverted to IPH- on scans 3 to 5. Of the 11 arteries that were IPH+ at baseline, IPH persisted in five (45.5%) throughout the study, whereas six (54.5%) demonstrated resolution of IPH before the final scan.
During the follow-up, three subjects (10.7%) developed new ischemic infarcts. Of these, one subject had no detectable IPH in either carotid artery at any time point. Another subject consistently exhibited ipsilateral IPH+ and contralateral IPH- across all scans. In the third subject, a new infarct was detected at the third scan, coinciding with the appearance of new ipsilateral IPH. Subsequently, the contralateral side transitioned from IPH- to IPH+ at the fourth scan. Overall, two new infarcts occurred among the 17 subjects with either unilateral or bilateral IPH.
Bilateral Symmetry in Plaque Evolution and the Disruptive Role of IPH
In the 22 subjects with bilateral analyzable arteries, the average %WV growth rate was -0.19 ± 0.74 %/year for the left carotid artery and 0.04 ± 0.54 %/year for the right. The long-term plaque growth rate of one side was moderately correlated (r = 0.54, p < 0.001) with the contralateral side. The significant correlation remained after adjusting for age, sex and statin use (r = 0.43, p = 0.005). Figure 1 provides an example of bilateral plaque burden increases over five scans. The association between plaque progression on the left and right sides is visualized in Figure 2.
The 22 subjects were further divided into two sub-groups based on the presence of IPH. In nine subjects, IPH was not detected at any time point, while 13 subjects had either unilateral or bilateral IPH during the study period. Within the IPH- subgroup, significant associations of plaque progression between one side and its contralateral part were found (r = 0.57, p < 0.001). However, in the group with IPH present, the correlation coefficient decreased, and the association no longer remained significant (r = 0.28, p = 0.187). This suggests that the presence of IPH disrupts the systemic, symmetric trend of bilateral plaque evolution. The correlation of bilateral long-term plaque progression in the two subgroups is visualized in Figure 2.
IPH’s Long-term Impact on Plaque Progression
Building on these findings, we further evaluated both the presence and the volume of IPH on %WV progression using a linear mixed-effects model. In the unadjusted model, the random effects for patients and arteries accounted for individual variability, with SDs of 4.3% and 1.89%, respectively. The SD of the residual term is 3.08%. This suggests considerable variability in %WV progression among patients and between bilateral arteries within the same patient.
The progression of %WV over time in arteries without IPH was minimal, with an average rate of -0.001 %/year (p = 0.895). However, the presence and occurrence of IPH at any point was a significant predictor, associated with a 2.34% absolute increase in %WV (p < 0.001). This indicates that patients with IPH had higher %WV compared to those without IPH and that the development of new IPH was significantly associated with an increase in plaque burden.
Furthermore, %HV was significantly associated with %WV, with a coefficient of 0.73 (p = 0.005). This means that a 1-unit increase in the log2-transformed %HV, equivalent to a doubling of the mean of %HV, resulted in a 0.73% absolute increase in %WV. These findings highlight the significant impact of both the presence and volume of IPH on the progression of carotid plaque burden over time, as illustrated in Figure 3. Figure 4 provides an example of plaque progression over four scans, highlighting the increasing plaque burden alongside the occurrence (scan 3) and subsequent volume expansion (scan 4) of IPH.
When adjusted for age, sex, and statin use, the linear mixed-effects model continued to show strong associations between IPH and %WV progression. The presence of IPH was associated with an adjusted 2.29% absolute increase in %WV (p < 0.001), and the %HV association persisted, with a coefficient of 0.74% (p = 0.004), similar to the unadjusted model. Age, sex, and statin use had minimal impact on %WV progression (p = 0.321, 0.568, and 0.710, respectively).
Discussion
By leveraging two deep learning-based segmentation pipelines on two large-coverage 3D MR vessel wall imaging (VWI) sequences, we comprehensively evaluated carotid plaque features across an average of five scans over six years in an asymptomatic patient cohort. The role of intraplaque hemorrhage (IPH) in long-term plaque progression can be characterized in two critical ways. First, while plaque progression on one side was significantly correlated with the contralateral side, reflecting the systemic nature of atherosclerotic disease, this symmetry was notably weakened in subjects with IPH, suggesting a more localized effect. Second, despite contemporary medical management, the occurrence and volume of IPH were strongly associated with plaque burden increase, whereas plaques without IPH remained stable over the study period. These findings deepen our understanding of the long-term and pivotal role IPH plays in the natural history of carotid atherosclerosis and underscore the importance of early IPH detection for optimized strategies in managing atherosclerosis.
Previous studies have demonstrated carotid plaque symmetry by comparing cross-sectional measurements of plaque burden between bilateral carotid arteries, as observed in ex vivo imaging studies19 and in vivo assessments using ultrasound20 and VWI21. Our study extends this understanding with longitudinal follow-up over a six-year period, which reveals not just static morphological symmetry but also symmetrical trends in bilateral plaque progression. Our results suggest that in patients with bilateral disease, monitoring the more advanced plaque could offer valuable prognostic insights into the less advanced contralateral plaque, as signs of progression in the more advanced plaque may indicate elevated bilateral risk.
Notably, as demonstrated by Li et al.21, IPH exhibits low concordance between left and right carotid plaques, with unilateral or bilateral IPH frequently reported in the current and other studies, underscoring the role of localized factors in its pathogenesis. Studies showed that IPH development, likely driven by the rupture of neo-vessels, is shaped by conditions such as hemodynamic forces22, calcification9, and inflammation23 within the plaque microenvironment. At the same time, systemic factors still significantly contribute to the development of IPH. For instance, the Rotterdam Study24, in a cross-sectional analysis, identified pulse pressure as the most significant predictor of IPH, independent of plaque burden and other cardiovascular risk factors. Similarly, another longitudinal study associated male sex and hypertension with an increase in IPH9. Understanding the interplay between these localized and systemic influences is critical for advancing our understanding of pathophysiology of IPH development and its impact on plaque stability.
IPH is a dynamic process with the potential for growth or resolution over time. Longitudinal follow-up studies have shown that its hyperintense signal on VWI, likely attributable to methemoglobin from hemoglobin breakdown after erythrocyte extravasation, can persist throughout the course of atherosclerotic disease. Previous studies7,25,26 using T1-weighted sequences reported IPH resolution rates of 0%, 0%, and 6% over 17 to 18 months of follow-up. More recently, one study9 found that IPH resolved in only 3 out of 89 scan pairs (3.4%) over an average of 1.8 years, utilizing the SNAP sequence which combines inversion-recovery gradient echo with phase-sensitive acquisition to enhance IPH-to-wall contrast and measurement accuracy27,28. Similarly, our study, using SNAP to detect IPH, demonstrated that all but one newly developed IPH persisted to the study end, around half arteries with pre-existing IPH continued to exhibit it throughout the study period, and hyperintense IPH signals may be observed for as long as six years. This persistence may be explained by ongoing erythrocyte extravasation from leaky neo-vessels29.
IPH is widely recognized as a strong indicator of carotid plaque vulnerability, with its presence linked to an increased risk of future ischemic cerebrovascular events30, independent of symptom status31 or degree of stenosis5,32. In our study, despite more than half of subjects exhibited IPH at some point during follow-up, the incidence of ischemic symptoms over at least five years was low. Histological studies29,33 have supported this observation by identifying neovasculature, a potential source of IPH, predominantly located in the adventitia away from the plaque surface, thus preserving surface integrity. This suggests that assessing plaque surface disruption34,35 could be an important additional metric for stratifying stroke risk, particularly in asymptomatic patients. Nonetheless, IPH has been consistently associated with a rapid increase in plaque burden7,9, and our findings added to the body of evidence by further highlighting the sustained and long-term impact of IPH presence and volume on the progression of atherosclerotic disease. Statin therapy improves cardiovascular outcomes and has demonstrated plaque regression36–38 in multiple imaging studies, by the mechanism of modulating inflammation and decreasing or even depleting the lipid content within plaques39–42. However, the recent AIM-HIGH MRI sub-study42 showed that, despite two years of intensive lipid-lowering therapy, carotid arteries with IPH experienced larger increases in lipid core size and greater reductions in lumen area compared to plaques without hemorrhage, which exhibited stabilized plaque burden and significant reductions in lipid content. This suggests the IPH’s interplay with a potential non-lipid-driven pathway of plaque progression. One study43 identified CD163+ macrophages, induced by IPH, as potential contributors to a positive feedback loop involving increased vascular permeability, inflammation, and further IPH increase. Our findings similarly highlighted the role of IPH in counteracting the systemic stabilizing effects of contemporary medical treatment. While non-IPH plaques remained stable, current therapies may be insufficient to prevent IPH-associated plaque burden increase. Further investigation into the mechanisms behind this non-lipid-driven pathway is essential and could pave the way for novel therapeutic approaches specifically targeting IPH. Overall, there remains a significant gap in longitudinal imaging studies investigating the pathophysiology of carotid atherosclerosis and IPH, as most studies are either short-term7,9,44 or involve long-term clinical follow-up without repeated imaging31,45. To the best of our knowledge, the longest imaging follow-up duration reported is 54 months6. However, this study was limited by a very small sample size and the use of 2D VWI, which provided inadequate coverage of the carotid bifurcation. This limitation was further compounded by challenges in identifying consistent imaging coverage for bilateral arteries and comparing multiple time points. Indeed, Geleri et al. reported that, when using the large-coverage 3D VWI sequence MERGE, only 54% of lesions were fully covered and 25% partially covered by 2D VWI in their analysis of 381 advanced plaques from the CARE-II study46. Recognizing these limitations, our study utilized MERGE, which provided adequate imaging coverage for bilateral and longitudinal analysis.
The emergence of deep learning-based imaging analysis also offers a promising solution to the labor-intensive nature of longitudinal VWI studies, reducing the need for extensive reviewer training and minimizing bias, as well-validated deep learning tools provide consistent and reproducible analyses across time points. While significant efforts have been made in segmenting the vessel wall on VWI10–12, robust algorithms for detailed plaque composition segmentation remain underdeveloped47. In our study, we successfully utilized deep learning to integrate multi- contrast VWI, demonstrating its potential to deliver more detailed and comprehensive assessments of plaque morphology and composition48.
Our study has a few limitations. First, although each patient underwent an average of five scans, the requirement to analyze the long-term evolution of atherosclerosis over more than five years restricted the number of participants in this retrospective analysis. Second, the retrospective nature of this analysis resulted in a relatively healthy and stabilized patient cohort, as individuals with severe symptoms—indicative of plaques experiencing dramatic reductions in stability— were more likely to drop out of the prospective cohort early. Third, despite extensive follow-up, the absolute number of newly detected or resolved hemorrhages remained relatively low, underscoring the need for larger studies to better understand these potentially course-altering events in the plaque natural history. Finally, future research should aim to develop more advanced deep learning algorithms capable of simultaneously segmenting additional plaque components, such as lipid contents and calcification49, to further elucidate the complex interplay between these plaque components and the pathobiology of atherosclerosis.
Conclusion
By leveraging deep learning-based segmentation on multi-contrast VWI, we successfully monitored long-term plaque dynamics across an average of five scans over a follow-up of six years in an asymptomatic cohort. Around half of the arteries with pre-existing IPH continued to exhibit it throughout the study period, while newly developed IPH also persisted, yet caused few new ischemic symptoms. Importantly, IPH disrupted the systemic symmetry of bilateral plaque progression: arteries without IPH demonstrated minimal growth under contemporary medical treatment, whereas the occurrence and volume of IPH were strongly associated with accelerated plaque burden progression.
Perspectives
Competency in medical knowledge
Demonstrates advanced knowledge of carotid atherosclerosis pathophysiology by leveraging deep learning-based segmentation on multi- contrast VWI to reveal that IPH disrupts bilateral plaque symmetry, persists over time, and significantly accelerates plaque burden progression, even under contemporary medical treatment. Competency in patient care: Emphasizes the importance of early detection of IPH using advanced imaging techniques and advocates for personalized management strategies to mitigate its impact on plaque progression and stroke risk.
Translational outlook
The present study advocates the adoption of longitudinal patient monitoring and the integration of deep learning-based imaging tools into clinical practice to enable early assessment of plaque burden and IPH, while advancing the development of IPH- targeted therapeutics to improve patient outcomes in carotid atherosclerosis.
Funding disclosure: This study was funded by National Institutes of Health under grant R01- HL103609, R01-NS125635 and R01-NS127317.
Abbreviations list
- IPH
- intraplaque hemorrhage
- VWI
- vessel wall magnetic resonance imaging
- CEA
- carotid endarterectomy
- MERGE
- motion-sensitized driven equilibrium prepared rapid gradient echo
- SNAP
- simultaneous non-contrast angiography and intraplaque hemorrhage imaging
- %WV
- percent wall volume
- %HV
- percent hemorrhage volume
Data Availability
All data produced in the present study are available upon reasonable request to the authors.