ORIGINAL RESEARCH article

Front. Mol. Biosci., 15 January 2026

Sec. Molecular Diagnostics and Therapeutics

Volume 12 - 2025 | https://doi.org/10.3389/fmolb.2025.1685377

Exploring the risk stratification of carotid plaque hemodynamic using four-dimensional blood flow technology

  • 1. The First Affiliated Hospital of Sun Yat-sen University, Guangzhou, China

  • 2. Beijing Tsinghua Changgung Hospital, Affiliated to Tsinghua University, Beijing, China

  • 3. IDG-McGowan Institute of Brain Science, Tsinghua University, Beijing, China

  • 4. Center for Biomedical Imaging Research, Department of Biomedical Engineering, School of Medicine, Tsinghua University, Beijing, China

Abstract

Introduction:

Carotid plaque rupture is a major cause of cerebrovascular events. This study explores the integration of hemodynamic parameters with structural biomarkers for improved risk stratification.

Methods:

Fifty-seven patients with moderate-to-severe carotid stenosis underwent 4D-flow magnetic resonance imaging (MRI) and high-resolution MRI. Hemodynamic parameters [wall shear stress (WSS) and velocity] were analyzed using GT-Flow software at upstream, throat, and downstream plaque regions. After comparison of characteristic values between symptomatic and asymptomatic plaques, variables with p < 0.1 were included in the multivariate logistic regression model to identify independent risk factors.

Results:

WSS was significantly higher at plaque throat (0.891 ± 0.422 Pa) and downstream (0.971 ± 0.587 Pa) versus upstream (0.649 ± 0.297 Pa; p < 0.001). Symptomatic plaques showed elevated 3D-WSSmean (1.041 ± 0.418 vs. 0.797 ± 0.402 Pa, p = 0.032), WSS up_max (1.345 ± 0.570 vs. 0.970 ± 0.383 Pa, p = 0.004), and stenosis velocity (31.7 ± 9.9 vs. 25.8 ± 10.3 cm/s, p = 0.036). The thin fibrous cap (TFC, OR = 5.34, p = 0.007) and normalized wall index (NWI, adjusted OR = 59.89, p = 0.029) independently predicted symptomatic plaques. The combined model (NWI + TFC + WSS down_max) predicted cerebral ischemic events within 6 months with an AUC of 0.809 (95% CI: 0.699–0.919).

Conclusion:

Integration of downstream hemodynamic profiling (WSS) with structural biomarkers (TFC and NWI) provides a robust stratification tool for cerebrovascular risk assessment. These quantitative parameters offer potential as molecular diagnostics for plaque vulnerability.

1 Introduction

Previous histopathological studies on atherosclerosis biopsies have revealed that histological subtyping provides limited prognostic information: overall histological categories remain similar between patients with different symptom types, and predictions of plaque occurrence, progression, and rupture based on histological subtyping are not entirely reliable (). Moreover, further pathological research has demonstrated an association between histological features of plaques in symptomatic carotid artery stenosis patients (such as lipid core size or fibrous cap inflammation) and plaque outcomes (; ; ). Particularly insightful is recent research indicating that once hemodynamic parameters are introduced, the prognostic role of plaque structural characteristics diminishes. Quantitative comparisons have shown that carotid arteries with a higher NWI have a 3.472-fold likelihood of stroke association, while elevated WSS increases the risk of stroke by 6.974-fold, highlighting the potentially greater impact of hemodynamic factors represented by WSS on plaque prognosis ().

Currently, the consensus for distinguishing high- and low-risk carotid artery plaques is primarily based on imaging structural features (; ; ). However, the influence of hemodynamic factors on the progression of atherosclerotic disease is increasingly emphasized and significantly correlated with prognosis. 4D-flow MRI provides distinct hemodynamic advantages over traditional modalities. Unlike ultrasound (limited by operator dependency and 2D sampling) or Computed Tomography Angiography (CTA) (providing static anatomical assessment only), this technology enables volumetric quantification of time-resolved 3D flow patterns throughout the cardiac cycle. Crucially, it permits direct calculation of multidirectional wall shear stress vectors—including axial, circumferential, and 3D resultant forces—which cannot be derived from conventional imaging. This capability for dynamic hemodynamic profiling offers unique insights into the biomechanical triggers of plaque vulnerability. Therefore, this study aims to utilize 4D flow technology and high-resolution magnetic resonance imaging (HR-MRI) to comprehensively characterize the structure (imaging biomarkers) and non-structural determinants (hemodynamic factors) of carotid artery plaques, enhancing our understanding of their involvement in stroke occurrence.

Our objective is to explore hemodynamic risk stratification of carotid artery plaques based on differential patient outcomes. Through this approach, we can further develop tailored imaging surveillance and intervention strategies for different populations to reduce the risk of stroke recurrence and other adverse outcomes.

2 Methods

2.1 Study design and participants

This study is a single-center prospective observational study, selecting patients confirmed to have carotid artery stenosis through screening with carotid ultrasound, CTA, or Digital Subtraction Angiography (DSA) as the study subjects. All participants met the following criteria: (1) aged >18 years; (2) presence of one or more atherosclerosis risk factors, including hypertension [systolic blood pressure (SBP) ≥140 mmHg, diastolic blood pressure (DBP) ≥90 mmHg, or self-reported use of any antihypertensive medication in the past 2 weeks], diabetes (fasting blood glucose ≥7 mmol/L, non-fasting blood glucose ≥11.1 mmol/L, or use of antidiabetic medication), or hyperlipidemia [total cholesterol (TC) ≥240 mg/dL, low-density lipoprotein cholesterol (LDL-C) ≥160 mg/dL, or high-density lipoprotein cholesterol (HDL-C) <40 mg/dL or lipid-lowering medication use]; and (3) asymptomatic carotid stenosis patients who developed symptoms or new ischemic stroke within 6 months. The exclusion criteria were rigorously applied, leading to the exclusion of 47 patients for the following reasons: (1) significant intracranial artery stenosis (>50%) on MRA (n = 9); (2) ischemic stroke of cardioembolic origin (n = 13); (3) non-atherosclerotic etiology, such as dissection or vasculitis (n = 8); (4) prior treatment with bypass or stenting before enrollment (n = 6); and (5) poor image quality (image quality score = 1) that precluded reliable analysis (n = 11).

Consequently, 57 patients with 81 carotid arteries were included in the final analysis (Figure 2). All patients underwent detailed clinical assessments, including comprehensive physical examinations, lifestyle evaluations, present illness history, and past medical history before enrollment. Patients underwent head MRI, High-Resolution Vessel Wall (HRVW)-MRI of the neck, and 4D flow MRI scans separately. HRVW-MRI collected imaging biomarkers, including stenosis severity, plaque length, lumen area, outer diameter area, normalized wall index (NWI) (plaque burden), eccentric index (EI), lipid-rich necrotic core (LRNC), intraplaque hemorrhage (IPH), thin fibrous cap (TFC), calcification, plaque surface irregularity, and plaque imaging subtype. 4D flow collected WSS (3D and axial), flow, and velocity at different positions of the plaque. The individual outcome was defined as positive if there was a new cerebrovascular event (clinically diagnosed TIA or homolateral brain infarction observed on DWI) within 6 months.

2.2 Imaging protocol

All participants underwent 4D blood flow MRI (independently developed) (), HRVW 3D T1- and T2-weighted imaging (T1WI and T2WI), 3D MERGE, SNAP sequence imaging, and cranial MRI. The scans were performed using a Philips Ingenia CX 3.0T MR scanner, with an 8-channel phased-array carotid plaque coil dedicated for scanning. The coronal 3D MERGE sequence had a TR/TE of shortest/shortest, an Field Of View (FOV) of 200 mm × 159 mm, a slice thickness of 0.8 mm, an interslice gap of −0.4 mm, a voxel size of 0.8 mm × 0.8 mm × 0.8 mm, a flip angle of 6°, fat suppression using the spectral presaturation with inversion recovery (SPIR) technique, and a scan time of 4 min 50 s. The coronal 3D T1 VISTA sequence had a TR/TE of shortest/shortest, an FOV of 250 mm × 181 mm, a slice thickness of 0.6 mm, an interslice gap of −0.3 mm, a voxel size of 0.6 mm × 0.6 mm × 0.6 mm, an echo train length of 37, a flip angle of 90°, fat suppression using the spectral attenuated inversion recovery (SPAIR) technique, and a scan time of 1 min 48 s. The coronal 3D T2 VISTA sequence had a TR/TE of 1,600 m/shortest, an FOV of 250 mm × 160 mm, a slice thickness of 0.6 mm, an interslice gap of 0 mm, a voxel size of 0.6 mm × 0.6 mm × 0.6 mm, an echo train length of 50, a flip angle of 90°, the SPAIR technique, and a scan time of 5 min 17 s. The coronal 3D SNAP sequence had a TR/TE of shortest/shortest, an FOV of 200 mm × 159 mm, a slice thickness of 0.8 mm, an interslice gap of 0 mm, a voxel size of 0.8 mm × 0.8 mm × 0.8 mm, a flip angle of 11°, and a scan time of 2 min 36 s. The axial 3D QFLOW sequence had a TR/TE of shortest/shortest, an FOV of 160 mm × 160 mm, a slice thickness of 1.2 mm, an interslice gap of 0 mm, a voxel size of 1 mm × 1 mm × 1.2 mm, a flip angle of 20°, and a single acquisition using the phase contrast (PC) method for volumetric and time-resolved contrast; velocity encoding (VENC) was set at 90 cm/s to avoid aliasing artifacts, and the data were reconstructed into 15 frames using view sharing. The scan time was 11 min 27 s, depending on each participant’s heart rate.

2.3 Image analysis and post-processing

4D flow MRI data were imported into GTFlow software (version 3.2, GyroTools) for preprocessing and parameter calculation. The analysis includes the following steps. (1) Preprocessing: preprocessing steps include automatic offset correction, signal aliasing correction, flow direction correction, and full dynamic preview images to identify and exclude poor-quality images. (2) Segmentation: the target vessel is segmented parallel to the centerline and tracked along the vessel. (3) Calculation: the analysis plane is placed perpendicular to the narrowest part of the carotid artery centerline and dynamically adjusted within 20 time points to measure the velocity (maximum and average), flow (maximum and average), axial WSS (maximum and average), and circumferential WSS (maximum and average) (Figure 1). WSS is calculated according to the method described by . 3D WSS reflects the total WSS along the plane tangent to the local vessel surface and is decomposed into axial and circumferential components. Axial and circumferential WSSs represent the WSS along the flow direction and the vessel circumference. Original image data from MRI and 4D flow were independently measured by separate radiologists and neurologists in a blinded manner and averaged. The average post-processing time per patient using GTFlow software was approximately 30–35 min.

FIGURE 1

To assess the reproducibility of the measurements, inter-observer reliability was evaluated between the radiologist and the neurologist. Intraclass correlation coefficients (ICCs) were calculated using a two-way random-effects model for absolute agreement for continuous variables (WSS and structural parameters). For categorical variables (e.g., plaque classification), Cohen’s kappa (κ) statistic was used. The reliability analysis demonstrated excellent agreement between the two readers: the ICC for 3D WSS mean was 0.92 (95% CI: 0.87–0.95), for NWI was 0.89 (95% CI: 0.82–0.93), and for TFC thickness was 0.91 (95% CI: 0.85–0.94). The agreement for plaque classification (types I–VI) was substantial (κ = 0.78). These results indicate high consistency in image interpretation and parameter quantification.

2.4 Statistical analysis

All data were processed using SPSS software. Continuous variables are expressed as the mean ± standard deviation or median (interquartile range), and categorical variables are expressed as the ““frequency (percentage).” For continuous variables, independent sample t-tests or Mann–Whitney U tests were used to compare the data between two groups, and categorical variables were analyzed using χ2 tests or Fisher’s exact tests. All statistical analyses were performed using SPSS (version 25.0.0 IBM). Independent sample t-tests were used to compare quantitative data between groups (such as stenosis severity), and chi-square tests were used to compare categorical variables between groups. WSS and velocity values were found to follow a normal distribution based on normality tests. Chi-square tests were used to compare categorical data between groups. Logistic regression analysis was conducted separately for the two defined outcomes. Independent sample t-tests were used to compare WSS (maximum and average values) at different locations within the plaque and wall shear stress at rupture and non-rupture locations. The predictive efficacy of WSS in the plaque outcome and individual prognosis was evaluated using the ROC curve. Binary regression analysis was conducted for high-risk imaging features using high and low WSS as outcome variables. All hypothesis tests are two-sided, and a p-value less than 0.05 is considered statistically significant.

2.5 Sample size calculation

PASS-15 calculation: Taking WSS as the reference, a two-sided test is required with α = 0.05, that is, the confidence level (1 − alpha) = 0.95; according to the literature standard, the standard deviation = 0.08 (), and the allowable errors are taken as 0.05 and 0.02, respectively.

When δ = 0.02, N = 62; when δ = 0.05, N = 10.

This study included 57 cases with 81 carotid arteries, which had sufficient statistical power (≥98%) to address this issue.

3 Result

3.1 Demographic baseline characteristics of subjects

The flowchart of this study is shown in Figure 2. This study included a cohort of elderly participants with an average age of 65 years. The predominant comorbidity was hypertension, observed in 68.4% of subjects, and the population was predominantly male (78.9%). A significant majority exhibited unstable plaques, with 87.7% classified as type IV–VI. The mean degree of carotid artery stenosis was 69.9% (Table 1).

FIGURE 2

TABLE 1

CharacteristicFrequency, percentage (mean ± SD)CharacteristicFrequency, percentage (mean ± SD)
Cohort characteristic (n = 57)Plaque histological classification (n = 81)
Age65.52 ± 6.241–II0, 0
Sex (male = 1)30, 78.9III7, 12.3
Hypertension26, 68.4IV–V24, 42.1
Diabetes12, 31.6VI26, 45.6
Statin history25, 65.8Plaque imaging features (n = 81)
Total cholesterol (mmol/L)3.63 ± 0.88IPH27, 47.4
Triglyceride1.23 ± 0.50Calcification28, 49.1
HDL-C1.15 ± 0.28Irregular40, 70.2
LDL-C2.06 ± 0.77TRFC34, 59.6
Smoke24, 63.2Lipid core area (mm2)8.59 ± 12.81
Drink22, 57.9LRNC37, 64.9
Plaque characteristic (n = 81)Proportion of lipid core area13.99% ± 0.12
Symptomatic plaque22, 38.6Hemodynamic parameters (n = 81)
Stenosis69.90% ± 0.143D WSS mean0.89 ± 0.42
Length (mm)20.47 ± 9.19Axi WSS mean0.42 ± 0.33
Lumen area21.64 ± 16.50Velocity28.08 ± 10.45
Cross-sectional area79.85 ± 44.41Flow4.49 ± 2.28
Wall area58.20 ± 30.52WSS up0.65 ± 0.30
Eccentric index0.65 ± 0.15
NWI0.74 ± 0.10WSS down0.97 ± 0.59

Demographic characteristics of study subjects at baseline.

3.2 Periodic changes in carotid artery wall shear stress—characterization by 4D-flow

WSS undergoes periodic changes within a cardiac cycle, reaching its peak during the systolic phase and being at its lowest during the diastolic phase. Axial WSS changes in synchrony with 3DWSS (Figure 3).

FIGURE 3

3.3 Comparison of shear stress upstream, at the stenosis throat, and downstream of carotid artery plaques

WSS measurements were conducted upstream, at the stenosis throat, and downstream of the carotid artery plaques. It was found that for both the average and maximum wall shear stress, the average and maximum wall shear stress at the stenosis throat and downstream of the plaque were significantly higher than those upstream of the plaque (all p < 0.05), while there was no significant difference in the average and maximum WSSs between the stenosis throat and downstream of the plaque (p = 0.795; p = 0.677). Although the differences were not significant, both the average and maximum WSSs at the downstream of the plaque (3DWSSmean 0.971 ± 0.587 N/m2 vs. 3DWSSmax 1.639 ± 0.879 N/m2) were higher than those at the stenosis throat of the plaque (3DWSSmean 0.891 ± 0.422 N/m2 vs. 3DWSSmax 1.492 ± 0.879 N/m2) (Table 2; Figure 4).

TABLE 2

The measured indicatorsLocation in plaques (I) (J)Mean difference (I-J)Standard errorSignificance95% confidence intervalFSig.
3D mean-WSSPlaque throatPlaque up0.2420.0680.0020.0760.4087.8570.001
(0.891 ± 0.422)Plaque down−0.0790.0960.795−0.3120.153
Plaque upstreamPlaque throat−0.2420.0680.002−0.4080.076
(0.649 ± 0.297)Plaque down−0.3220.0870.001−0.5340.109
Plaque downstreamPlaque throat0.0790.0960.795−0.1530.312
(0.971 ± 0.587)Plaque up0.3210.0870.001−0.1090.534
3D max-WSSPlaque throatPlaque up0.3770.1080.0020.1140.6398.6580.000
(1.492 ± 0.650)Plaque down−0.1460.1440.677−0.4980.205
Plaque upstreamPlaque throat−0.3770.1080.002−0.640−0.114
(1.115 ± 0.495)Plaque down−0.5230.1340.001−0.849−0.198
Plaque downstreamPlaque throat0.1460.1450.677−0.2050.498
(1.639 ± 0.879)Plaque up0.5230.1340.0010.1980.849

Comparison of 3D mean and maximum wall shear stress at different locations of carotid plaques.

Plaque up represents the plaque upstream; plaque down represents the plaque downstream.

FIGURE 4

3.4 Comparison of inter-group characteristic values of symptomatic plaques based on the clinical definition of new ischemic events

Compared with the non-symptomatic plaque group, the symptomatic plaque group showed significant differences in imaging markers (Table 3): stenosis degree, lumen area, and NWI (plaque burden) (p < 0.05). However, when comparing with the imaging markers, more wall shear stress-related hemodynamic markers, such as 3D WSS mean (p = 0.032), 3D WSS max (p = 0.014), Axi WSS max (p = 0.079), upstream average WSS (p = 0.010), upstream maximum WSS (p = 0.004), and stenosis velocity (p = 0.036) showed significant differences between symptomatic and non-symptomatic high-risk plaques based on the definition of clinical ischemic events. Variables with p < 0.1 were included in the multivariate logistic regression model.

TABLE 3

CharacteristicIschemia group
(n = 32)
Non-ischemia group (n = 49)Sig (bilateral)Mean differenceStandard error difference95% CI
Lower limitHigher limit
Plaque characteristic
 Stenosis0.756 ± 0.1470.662 ± 0.1150.009−0.0940.035−0.164−0.025
 Length (mm)20.781 ± 12.02120.277 ± 7.0520.842−0.5052.522−5.5594.549
 Lumen area (cm2)16.350 ± 14.99724.968 ± 16.7310.0548.6194.378−0.15517.392
 CSA (cm2)75.259 ± 46.51382.729 ± 43.4730.5417.46912.150−16.88131.820
 Wall area (cm2)58.909 ± 33.80157.760 ± 28.7750.891−1.1498.378−17.93815.640
 Eccentric index0.679 ± 0.1300.632 ± 0.1620.256−0.0470.041−0.1290.035
 NWI0.794 ± 0.0980.706 ± 0.0860.001−0.0880.025−0.137−0.038
 Lipid core area (mm2)10.763 ± 8.1057.221 ± 7.9290.314−3.5423.485−10.5253.441
 Proportion of lipid core area0.141 ± 0.1200.139 ± 0.1280.968−0.0010.034−0.0690.067
Hemodynamic parameter
 3D WSS min0.580 ± 0.2960.472 ± 0.2760.968−0.0010.034−0.0690.067
 3D WSS mean1.041 ± 0.4180.797 ± 0.4020.032−0.2440.111−0.467−0.022
 3D WSS max1.754 ± 0.6631.327 ± 0.5930.014−0.4270.169−0.765−0.088
 Axi WSS min0.104 ± 0.3950.018 ± 0.2570.369−0.0860.095−0.2790.107
 Axi WSS mean0.497 ± 0.3980.367 ± 0.2660.144−0.1300.088−0.3070.046
 Axi WSS max1.105 ± 0.6810.836 ± 0.4560.079−0.2690.150−0.5700.032
 WSS up mean0.774 ± 0.3450.570 ± 0.2350.010−0.2040.077−0.358−0.050
 WSS up max1.345 ± 0.5700.970 ± 0.3830.004−0.3750.126−0.628−0.123
 WSS down mean1.067 ± 0.5840.910 ± 0.5890.329−0.1570.160−0.4770.163
 WSS down max1.813 ± 0.8491.529 ± 0.8930.240−0.2830.238−0.7610.194
 Velocity31.714 ± 9.90925.803 ± 10.2580.036−5.9102.755−11.431−0.388
 Vel up24.967 ± 8.81622.013 ± 6.7550.159−2.9542.070−7.1031.194
 Vel down32.267 ± 11.20830.229 ± 13.1560.550−2.0373.387−8.8244.750
 Flow4.175 ± 2.4914.688 ± 2.1570.4140.5130.623−0.7361.762
 Flow up5.101 ± 2.5475.333 ± 1.8530.6930.2320.583−0.9381.401
 Flow down3.395 ± 2.0403.932 ± 2.1120.3480.5370.567−0.6001.674

Comparison of characteristic values of symptomatic plaques based on the clinical definition of new ischemic events (continuous variables).

CSA, cross-sectional area; 3D WSS min/mean/max represent the minimum, average, and maximum wall shear stress resultant forces, respectively; Axi WSS min/mean/max represent the minimum, average, and maximum axial wall shear stresses, respectively; WSS up represents the upstream wall shear stress of the plaque; WSS down represents the downstream wall shear stress of the plaque; Velocity, Vel up, and Vel down represent the flow velocity at the stenosis site, upstream of the plaque, and downstream of the plaque, respectively; Flow, Flow up, and Flow down represent the flow volume at the stenosis site, upstream of the plaque, and downstream of the plaque, respectively. Bold values indicate statistical significance (p < 0.05).

In the categorical variable data between the symptomatic and asymptomatic plaque groups (Table 4), there were significant differences in drinking history (p = 0.075), calcification (p = 0.082), thin fibrous cap (p = 0.007), and plaque imaging classification (p = 0.016) between the groups. Among the calculated OR values, the baseline data’s TFC (OR = 5.344, 95% CI = 1.499–19) was the most significant risk factor.

TABLE 4

CharacteristicSig. (two-tailed)OR95% CI of OR
Lower limitHigher limit
Cohort characteristic
 Sex0.3352.1920.5229.201
 Hypertension0.7381.2220.3763.973
 Diabetes0.5381.4290.4584.453
 Smoke0.2352.0000.6326.333
 Drink0.0752.8240.8958.906
 Statin history0.4650.6620.2182.008
Plaque imaging feature
 IPH0.1602.1670.7316.419
 Calcification0.0822.6250.8737.893
 Irregular0.7381.2220.3763.973
 TRFC0.0075.3441.49919.052
 LRNC0.6821.2660.4093.917
 Classification0.016---

Comparison of eigenvalues between symptomatic and asymptomatic plaque groups (chi-square test).

Bold values indicate statistical significance (p < 0.1).

There was a nearly significant difference in the imaging classification of plaques (I–VI) between symptomatic plaques and asymptomatic plaques (one-way ANOVA, p = 0.054) (Table 5). The proportion of unstable plaques (IV–V and VI) in the high-risk plaque group was significantly higher than that of stable plaques (I–II) (p = 0.003 vs. p < 0.001). However, there was no significant difference in the proportion between hemorrhagic plaques and non-hemorrhagic plaques among unstable plaques (p = 0.766).

TABLE 5

Plaque histological classificationMean difference (I-J)Standard errorSignificance95% confidence intervalFSig.
I–IIIV–V−0.3750.1010.003−0.63−0.123.0810.054
VI−0.5000.1000.000−0.76−0.24
IV–VI–II0.3750.1010.0030.120.63
VI−0.1250.1420.766−0.480.23
VII–II0.5000.1000.0000.240.76

Comparison of the imaging classification of plaques between the symptomatic plaque group and the asymptomatic plaque group.

Bold values indicate statistical significance (p < 0.05).

3.5 Binary logistic regression analysis of symptomatic plaques based on the clinical definition of new ischemic events

By comparing the above inter-group indicators, parameters with p < 0.1 were selected [stenosis degree (p = 0.009), lumen area (p = 0.054), NWI (p = 0.001), 3D WSS mean (p = 0.032), 3D WSS max (p = 0.014), Axi WSS max (p = 0.079), WSS up mean (p = 0.010), WSS up max (p = 0.004), velocity (p = 0.036), calcification, plaque imaging classification, and TFC], and the imaging parameters of interest IPH and WSS down max were included in the regression analysis. The WSS parameters at the same site of the plaque were only included in the logistic regression model as the maximum value.

Binary stepwise regression showed that only the plaque imaging markers—NWI (p = 0.043) and TFC (p = 0.021)—were independent risk factors for symptomatic high-risk plaques based on clinical ischemic events (Table 6).

TABLE 6

CharacteristicBStandard errorSignificanceOR95% CI of OR
Plaque radiologic characteristic
 Plaque classification1.0290.9880.2982.7990.40319.428
 Stenosis−0.4944.5510.9140.6100.00045.468
 Lumen area0.0640.0450.1531.0660.9771.164
 NWI22.02310.8630.04336.2922.08164.000
 IPH−1.6301.3310.2210.1960.0142.659
 Calcification−0.5280.9270.5690.5900.0963.627
 TRFC2.2620.9780.0219.6001.41165.332
Hemodynamic index
 3DWSS mean0.0764.0010.9851.0790.0002743.951
 3DWSS max0.5792.2760.7991.7840.021154.361
 WSS up mean−1.9094.8680.6950.1480.0002064.533
 WSS up max2.0732.6420.4337.9510.0451410.768
 WSS down mean−2.2322.3350.3390.1070.00110.437
 WSS down max0.6441.2860.3091.9040.15323.663
 Axi WSS max1.0441.0450.3182.8420.36722.032
 Velocity narrow0.0860.0840.3021.0900.9251.284

Symptomatic plaque defined by clinical new ischemic events—binary logistic regression analysis.

Velocity narrowing indicates the flow velocity at the smallest lumen area. The remaining hemodynamic parameters are the same as those in Table 3.

Combining the inter-group comparison and regression analysis, if the indicators of asymptomatic patients, including stenosis degree, lumen area, NWI, 3D WSS mean, 3D WSS max, Axi WSS max, WSS up mean, WSS up max, changes in the flow velocity at the stenosis site, and unstable plaques (types IV–VI), are abnormal, especially when NWI or TFC occurs, it is necessary to be vigilant about the progression of the plaque and the recurrence of ischemic events.

3.6 Diagnosis of multicollinearity

As shown in Table 7, after reducing the independent variables included in the model, none of the variables had a tolerance (Tol) of less than 0.1 or a variance inflation factor (VIF) greater than 10, indicating that multicollinearity did not exist. A robustness analysis was conducted by correcting for confounding factors using logistic regression.

TABLE 7

CharacteristicUnstandardized coefficientStandardization coefficienttSig.Collinearity statistics
BStandard errorBetaToleranceVIF
Plaque classification0.0650.1330.0920.4880.6280.3832.613
Stenosis0.2770.6640.0760.4180.6780.4032.484
Wall area0.0060.0060.2071.0580.2960.3502.856
NWI2.1421.1590.4351.8490.0710.2434.117
IPH−0.1300.172−0.134−0.7590.4520.4322.313
Calcification−0.0540.134−0.055−0.4010.6900.7091.410
TFC0.2870.1290.2892.2290.0310.7991.251
3DWSSmax0.1410.2140.1860.6580.5140.1685.961
WSS up max0.2380.1550.2391.5300.1330.5511.816
WSS down max−0.0990.105−0.176−0.9420.3520.3832.611
Axi WSS max0.0420.1560.0480.2680.7900.4202.382
Velocity0.0040.0100.0770.3650.7170.3043.285

Diagnosis of multicollinearity in binary logistic regression variables of symptomatic plaques.

3.7 Robustness analysis—adjustment for confounding factors

After correcting for confounding factors (degree of stenosis, IPH, calcification, and imaging classification of plaques), the binary logistic regression analysis once again confirmed that the imaging markers of plaques, namely, NWI (p = 0.029) and TFC (p = 0.020), were independent risk factors for symptomatic high-risk plaques according to the definition of clinical ischemic events (Table 8).

TABLE 8

CharacteristicBStandard errorSig.OR95% CI of the OR value
Plaque radiologic characteristic
 Plaque classification−20.35812424.5220.9990.0000.000-
 Stenosis−1.3934.6320.7640.2480.000217.168
 Lumen area0.0720.0460.1191.0750.9821.178
 NWI24.81611.5930.03259.8908.119441.000
 IPH1.1951.3230.3663.3040.24744.184
 Calcification0.7980.9570.4052.2210.34014.501
 TFC−2.3601.0130.0200.0940.0130.687
Hemodynamic index
 3D WSS max0.8191.3830.5542.2690.15134.122
 WSS up max1.2881.1450.2613.6270.38434.233
 WSS down max−0.6390.6280.3090.5280.1541.807
 Axi WSS max1.0291.0320.3192.7970.37021.141
 Velocity0.0310.0690.6561.0310.9011.181

Binary logistic regression analysis for correcting the symptoms of plaques affected by confounding factors.

3.8 Predictive efficacy of carotid plaque wall shear stress and independent risk factors for new-onset stroke

In the previous study, we found that the WSS down max was an independent risk factor for plaque rupture. Therefore, we incorporated this hemodynamic factor into the prediction model. The predictive ability of the downstream max value of WSS for new-onset stroke was AUC = 0.622, with a cutoff value of 1.302. NWI (AUC = 0.740) and TFC (AUC = 0.681) also showed good predictive efficacy. The combined index (NWI + TFC + downstream max of WSS) has the highest predictive ability (AUC = 0.809), significantly higher than downstream max of WSS and TFC (p < 0.05), suggesting that the combined screening of histology and hemodynamics is more capable of understanding the prognosis of individuals with carotid artery stenosis (Figure 5, Table 9).

FIGURE 5

TABLE 9

Test variableAUCpStandard errorAsymptotic sig.Asymptotic 95% CIYouden index maxCutoff value
Lower limitUpper limit
WSS down max0.6220.0350.0770.1230.4720.7720.2871.302
NWI0.7400.2130.0700.0020.6040.8770.4430.823
TFC0.6810.0200.0720.0230.5390.8220.3610.5
Joint indicator0.809-0.05600.6990.9190.4860.225

Predictive efficacy of WSS and independent risk factors for new-onset stroke.

4 Discussion

This study conducted detailed measurements of WSS upstream, at the throat region (the area with the smallest lumen area), and the downstream of the carotid artery plaques and discovered some interesting phenomena. The results showed that both the average WSS and the maximum WSS at the laryngeal region of the plaque and downstream of the plaque were significantly higher than those upstream of the plaque (p < 0.05). Notably, although not statistically significant, WSS at the downstream was numerically higher than that at the throat.

These findings differ from the research results on coronary artery plaques (; ). In the coronary arteries, the downstream of the plaque has a lower WSS than the upstream due to blood flow disturbance and energy loss (). This difference between carotid artery plaques and coronary artery plaques may be due to the unique anatomical structure of the carotid artery. The corner at the bifurcation of the carotid artery may cause vortices to form downstream of the plaque, which could be one of the reasons for the higher WSS downstream of the plaque in this study. This continuous spiral flow pattern has been reported by , and it has been shown to be beneficial in patients with ICA stenosis ≤10% as this “zigzag” of the carotid artery can help protect the ICA from atherosclerosis (). Dilba characterized the regional distribution of WSS based on the baseline MRI contour and flow on symptomatic carotid arteries and also found that high WSS was concentrated in the downstream region of the plaque, which was consistent with the speculated vortex-shaped region.

To verify this hypothesis, this study conducted further analysis by reconstructing streamlines and pathlines (Figure 4). The analysis results supported our assumption that the special geometric shape at the carotid artery bifurcation did indeed form vortices inside the lumen downstream of the plaque (as indicated by the red arrows), which might be the main reason for the higher WSS downstream of the plaque. In addition to the above explanation, other factors may also affect the WSS distribution upstream and downstream of the carotid artery plaque. For example, the morphology and composition of the carotid artery plaque, the degree of vascular remodeling, and local hemodynamic characteristics could all affect the distribution of WSS (14). The significance of this result lies in the fact that it reveals the uniqueness of the WSS distribution upstream and downstream of the carotid artery plaque compared to that of coronary artery plaques, which is of great significance for understanding the progression of carotid atherosclerosis and clinical management. The downstream area of the carotid artery plaque may have higher WSS due to the helical blood flow pattern, which plays a protective role in mild atherosclerosis and makes it more prone to plaque rupture in moderate to severe carotid artery stenosis. The higher wall shear stress downstream of the plaque is associated with excessive expression of MMP-9 (). This mechanotransduction pathway indicates that quantification of wall shear stress can serve as a non-invasive alternative indicator for molecular fragility biomarkers. Our research results indicate that the degree of stenosis, lumen area, NWI, TFC, 3D WSS mean/max, Axi WSS max, WSS up mean, WSS up max, changes in flow velocity at the stenosis site, and the occurrence of unstable plaques (types IV–VI) are all risk factors for symptomatic plaques that require vigilance. In particular, NWI and TFC have been confirmed as independent risk factors through regression analysis. Their identification as independent predictors aligns with their well-established roles in plaque instability and thromboembolism. These risk factors may affect the stability of the plaque and subsequently influence the occurrence of ischemic events.

The independent predictive values of NWI and TFC can be explained by their central roles in plaque vulnerability. A high NWI signifies substantial plaque burden, which is often associated with inflammatory activity and a larger lipid core, increasing susceptibility to rupture under hemodynamic stress (; ; ). Concurrently, TFC represents a critical structural failure, where a thin or ruptured fibrous cap directly facilitates thrombosis and the release of emboli, leading to cerebral ischemia (; ). Although WSS was not an independent factor in the final model, hemodynamic forces may indirectly contribute by influencing the morphological and compositional development of these high-risk features ().

The predictive power of the downstream maximum of WSS (WSS down max) for new-onset stroke showed moderate strength (AUC = 0.622), with a cutoff value of 1.302. In addition, plaque morphological indicators such as TFC and NWI also showed good predictive efficacy. The combined model (NWI + TFC + WSS down_max) achieved significantly enhanced predictive ability (AUC = 0.809), indicating that combined screening of histology and hemodynamics can more comprehensively assess the prognosis of patients with carotid artery stenosis. Our findings, enabled by 4D-flow MRI, show that WSS-based hemodynamic assessment predicts the stroke risk as effectively as structural biomarkers. This highlights its potential as a non-invasive tool for evaluating plaque vulnerability based on direct biomechanical profiling. It is recommended that in type IV–VI plaques with large plaque load, thin fibrous caps, and calcification as observed on carotid magnetic resonance imaging, mechanical indicators related to WSS be measured as a routine screening method for stroke risk. The combined model (NWI + TFC + WSS down max) achieved a clinical-grade predictive value (AUC = 0.809; cutoff = 0.225), suggesting its utility for identifying high-risk plaques requiring intensive intervention. However, implementation barriers remain: 4D-flow acquisition time (11 min 27 s) currently limits emergency applications, and standardized WSS quantification protocols are still evolving (). Future integration of accelerated 4D-flow sequences with deep learning reconstruction could enable point-of-care hemodynamic risk stratification (). This hemodynamic-structural paradigm shift enables the precision identification of vulnerable plaques through three synergistic biomarkers: (1) NWI quantifying plaque burden, (2) TFC indicating cap integrity, and (3) downstream WSS mapping biomechanical stress (Figure 4G). Validation in multicenter cohorts established standardized risk thresholds for clinical deployment.

This study has limitations inherent to its cross-sectional design, which precludes the establishment of causal relationships between plaque characteristics and outcomes. Furthermore, our analysis focused on the culprit plaque and did not account for the potential influence of total plaque burden or collateral circulation status on overall stroke risk. Studies have shown a significant association between the total number of plaques and the recurrence of acute stroke events (). However, the role of the number of plaques in plaque rupture has been proven to be limited (). Future longitudinal studies tracking the evolution of hemodynamic and structural parameters are warranted.

5 Conclusions

Our findings establish 4D-flow hemodynamic profiling as a transformative diagnostic tool for cerebrovascular risk stratification. The integration of downstream wall shear stress quantification with structural biomarkers (TFC and NWI) creates a novel paradigm where biomechanical forces directly inform molecular vulnerability assessment. This synergistic approach enables the precision identification of high-risk plaques warranting targeted interventions, advancing the frontier of molecular diagnostics in cerebrovascular disorders.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by the Ethics Committee of Beijing Tsinghua Changgung Hospital (approval no. [22028-0-02]). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

S-yZ: Data curation, Investigation, Project administration, Writing – original draft, Writing – review and editing. H-lZ: Data curation, Investigation, Software, Validation, Writing – original draft. X-wS: Data curation, Supervision, Validation, Writing – original draft, Writing – review and editing. JW: Data curation, Supervision, Validation, Visualization, Writing – review and editing. RL: Formal analysis, Investigation, Supervision, Visualization, Writing – review and editing.

Funding

The authors declare that no financial support was received for the research and/or publication of this article.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2025.1685377/full#supplementary-material

References

Summary

Keywords

4D-flow magnetic resonance imaging, high-resolution magnetic resonance imaging, carotid artery plaque, hemodynamics, risk stratification

Citation

Zhang S, Zhao H, Song X, Wu J and Li R (2026) Exploring the risk stratification of carotid plaque hemodynamic using four-dimensional blood flow technology. Front. Mol. Biosci. 12:1685377. doi: 10.3389/fmolb.2025.1685377

Received

13 August 2025

Revised

13 November 2025

Accepted

14 November 2025

Published

15 January 2026

Volume

12 - 2025

Edited by

Lijun Wang, Naval Medical University, China

Reviewed by

Xuewei Xie, Capital Medical University, China

Dongwei Dai, Fudan University, China

Updates

Copyright

*Correspondence: Xiao-wei Song, ; Jian Wu, ; Rui Li,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics