Super-Resolution Ultrasound Imaging Reveals Spatiotemporal Vascular and Glial Remodelling Following Spinal Cord Injury

Background Spinal cord injury (SCI) induces complex secondary pathophysiological cascades, including microvascular disruption, glial activation, and chronic neuroinflammation. These vascular changes are dynamic and spatially heterogeneous across acute, subacute, and chronic phases. However, limitations in in vivo imaging resolution have impeded precise, longitudinal evaluation of spinal microcirculation. Objective To quantify post-SCI alterations in spinal cord vessel density, flow speed, and morphology over time using ultrasound localization microscopy (ULM), and to correlate these changes with chronic astrocytic reactivity assessed by GFAP expression. Methods Eleven Sprague-Dawley rats underwent T9–T11 laminectomy; nine received a standardized 200 kdyn contusion using the Infinite Horizons Impactor and were assigned to acute (24h), subacute (7d), and chronic (21d) groups (n = 3/group), while two served as uninjured controls. ULM was performed using a Vevo 3100 high-frequency ultrasound system with a 40 MHz transducer following intravenous Definity® microbubble injection. Vessel density (VD) and flow speed (FS) were quantified at the lesion epicenter and adjacent rostral/caudal regions. Arterial and venous compartments were distinguished by flow direction. Histological sections were stained for glial fibrillary acidic protein (GFAP), and %Area was quantified across phases. One-way ANOVA and Tukey’s test evaluated between-group differences; Pearson correlation identified phase-specific relationships between vascular metrics and chronic GFAP expression. Results ULM provided super-resolved imaging of spinal microvasculature at 5–10 µm resolution. Venous VD significantly declined from 0.92 ± 0.18 (acute) to 0.42 ± 0.23 (subacute) and 0.11 ± 0.07 (chronic, p < 0.001), particularly in dorsal regions. Arterial VD remained stable across timepoints (p > 0.05), while FS demonstrated transient reductions in both arterial and venous systems during the subacute phase. Chronic venous FS partially recovered by Day 21. GFAP-positive area increased significantly in the chronic group (3.3 ± 0.87%) compared to acute (1.52 ± 0.42%) and subacute (1.40 ± 0.05%) (p < 0.05), indicating delayed astrocyte reactivity. Correlation analysis revealed significant inverse associations between chronic GFAP %Area and early ventral venous density (r = − 1.0, p < 0.05), as well as positive correlations with chronic ventral arterial FS (r = 1.0, p < 0.01), suggesting early microvascular deficits may predict long-term gliosis. Conclusion This study demonstrates the feasibility and power of ULM to quantify dynamic microvascular remodeling in SCI with unprecedented resolution. Venous rarefaction precedes chronic astrocytosis, while arterial flow speed partially recovers over time. These findings highlight key spatiotemporal biomarkers of injury progression and support the use of ULM as a noninvasive platform for monitoring SCI pathophysiology and guiding therapeutic timing.