Specific and Nonuniform Brain States during Cold Perception in Mice.

The quest to decode the complex supraspinal mechanisms that integrate cutaneous thermal information in the central system is still ongoing. The dorsal horn of the spinal cord is the first hub that encodes thermal input which is then transmitted to brain regions via the spinothalamic and thalamocortical pathways. So far, our knowledge about the strength of the interplay between the brain regions during thermal processing is limited. To address this question, we imaged the brains of adult awake male mice in resting state using functional ultrasound imaging during plantar exposure to constant and varying temperatures. Our study reveals for the first time the following: (1) a dichotomy in the response of the somatomotor-cingulate cortices and the hypothalamus, which was never described before, due to the lack of appropriate tools to study such regions with both good spatial and temporal resolutions. (2) We infer that cingulate areas may be involved in the affective responses to temperature changes. (3) Colder temperatures (ramped down) reinforce the disconnection between the somatomotor-cingulate and hypothalamus networks. (4) Finally, we also confirm the existence in the mouse brain of a brain mode characterized by low cognitive strength present more frequently at resting neutral temperature. The present study points toward the existence of a common hub between somatomotor and cingulate regions, whereas hypothalamus functions are related to a secondary network.

Covariations between pupil diameter and supplementary eye field activity suggest a role in cognitive effort implementation.

In both human and nonhuman primates (NHP), the medial prefrontal region, defined as the supplementary eye field (SEF), can indirectly influence behavior selection through modulation of the primary selection process in the oculomotor structures. To perform this oculomotor control, SEF integrates multiple cognitive signals such as attention, memory, reward, and error. As changes in pupil responses can assess these cognitive efforts, a better understanding of the precise dynamics by which pupil diameter and medial prefrontal cortex activity interact requires thorough investigations before, during, and after changes in pupil diameter. We tested whether SEF activity is related to pupil dynamics during a mixed pro/antisaccade oculomotor task in 2 macaque monkeys. We used functional ultrasound (fUS) imaging to examine temporal changes in brain activity at the 0.1-s time scale and 0.1-mm spatial resolution concerning behavioral performance and pupil dynamics. By combining the pupil signals and real-time imaging of NHP during cognitive tasks, we were able to infer localized cerebral blood volume (CBV) responses within a restricted part of the dorsomedial prefrontal cortex, referred to as the SEF, an area in which antisaccade preparation activity is also recorded. Inversely, SEF neurovascular activity measured by fUS imaging was found to be a robust predictor of specific variations in pupil diameter over short and long-time scales. Furthermore, we directly manipulated pupil diameter and CBV in the SEF using reward modulations. These results bring a novel understanding of the physiological links between pupil and SEF, but it also raises questions about the role of anterior cingulate cortex (ACC), as CBV variations in the ACC seems to be negligible compared to CBV variations in the SEF.

Altered cortical trigeminal fields excitability by spreading depolarization revealed with in vivo functional ultrasound imaging combined with electrophysiology.

Spreading depolarization (SD), usually termed cortical spreading depression has been proposed as the pathophysiological substrate of migraine aura and as an endogenous trigger of headache pain. The links between neurovascular coupling and cortical craniofacial nociceptive activities modulated by SD were assessed by combining in vivo local field potential (LFPs) recordings in the primary somatosensory cortex (S1) with functional ultrasound (fUS) imaging of S1 and caudal insular (INS) cortices of anesthetized male rats. A single SD wave triggered in the primary visual cortex elicited an ipsilateral, quadriphasic hemodynamic and electrophysiological response in S1 with an early phase consisting of concomitant increases of relative cerebral blood volume (rCBV) and LFPs. A transient hypoperfusion was then correlated with the beginning of the neuronal silence, followed by a strong increase of rCBV while synaptic activities remained inhibited.LFPs and rCBV recovery period was followed by a progressive increase in S1 and INS baseline activities and facilitation of cortical responses evoked by periorbital cutaneous receptive fields stimulation. Sensitization of cortical ophthalmic fields by SD was bilateral, occurred with precise spatiotemporal profiles and was significantly reduced by pre-treatment with a NMDA antagonist. Combined high-resolution assessing of neurovascular coupling and electrophysiological activities has revealed a useful preclinical tool for deciphering central sensitization mechanisms involved in migraine attacks.SIGNIFICANCE STATEMENTA crucial unsolved issue is whether visual aura and migraine headache are parallel or sequential processes. Here we show that a single spreading depolarization (SD) wave triggered from the primary visual cortex is powerful enough to elicit progressive, sustained increases of hemodynamic and sensory responses to percutaneous periorbital noxious stimuli recorded in S1 and Insular ophthalmic fields. Sensitization of cortical ophthalmic fields by SD was bilateral, occurred with precise spatiotemporal profiles and was significantly reduced by pre-treatment with a NMDA antagonist. Combined high-resolution assessing of neurovascular coupling and electrophysiological activities has revealed a useful preclinical tool for deciphering central sensitization mechanisms involved in migraine attacks.

Functional imaging evidence for task-induced deactivation and disconnection of a major default mode network hub in the mouse brain.

The default mode network (DMN) has been defined in functional brain imaging studies as a set of highly connected brain areas, which are active during wakeful rest and inactivated during task-based stimulation. DMN function is characteristically impaired in major neuropsychiatric diseases, emphasizing its interest for translational research. However, in the mouse, a major preclinical rodent model, there is still no functional imaging evidence supporting DMN deactivation and deconnection during high-demanding cognitive/sensory tasks. Here we have developed functional ultrasound (fUS) imaging to properly visualize both activation levels and functional connectivity patterns, in head-restrained awake and behaving mice, and investigated their modulation during a sensory-task, whisker stimulation. We identified reproducible and highly symmetric resting-state networks, with overall connectivity strength directly proportional to the wakefulness level of the animal. We show that unilateral whisker stimulation leads to the expected activation of the contralateral barrel cortex in lightly sedated mice, while interhemispheric inhibition reduces activity in the ipsilateral barrel cortex. Whisker stimulation also leads to elevated bilateral connectivity in the hippocampus. Importantly, in addition to functional changes in these major hubs of tactile information processing, whisker stimulation during genuine awake resting-state periods leads to highly specific reductions both in activation and interhemispheric correlation within the restrosplenial cortex, a major hub of the DMN. These results validate an imaging technique for the study of activation and connectivity in the lightly sedated awake mouse brain and provide evidence supporting an evolutionary preserved function of the DMN, putatively improving translational relevance of preclinical models of neuropsychiatric diseases.

Functional ultrasound imaging of deep visual cortex in awake nonhuman primates.

Deep regions of the brain are not easily accessible to investigation at the mesoscale level in awake animals or humans. We have recently developed a functional ultrasound (fUS) technique that enables imaging hemodynamic responses to visual tasks. Using fUS imaging on two awake nonhuman primates performing a passive fixation task, we constructed retinotopic maps at depth in the visual cortex (V1, V2, and V3) in the calcarine and lunate sulci. The maps could be acquired in a single-hour session with relatively few presentations of the stimuli. The spatial resolution of the technology is illustrated by mapping patterns similar to ocular dominance (OD) columns within superficial and deep layers of the primary visual cortex. These acquisitions using fUS suggested that OD selectivity is mostly present in layer IV but with extensions into layers II/III and V. This imaging technology provides a new mesoscale approach to the mapping of brain activity at high spatiotemporal resolution in awake subjects within the whole depth of the cortex.

4D functional ultrasound imaging of whole-brain activity in rodents.

We extended the capabilities of functional ultrasound to whole-brain four-dimensional (4D) neuroimaging. Our multiplane-wave transmission scheme on matrix arrays at thousands of frames per second provides volumetric recordings of cerebral blood volume changes at high spatiotemporal resolution. We illustrated the approach in rats while providing multiple sensory stimuli, for 4D functional connectivity and during instantaneous tracking of epileptiform events.

Pharmaco-fUS: Quantification of pharmacologically-induced dynamic changes in brain perfusion and connectivity by functional ultrasound imaging in awake mice.

There is a critical need for reliable quantitative biomarkers to assess functional brain alterations in mouse models of neuropsychiatric diseases, but current imaging methods measuring drug effects through the neurovascular coupling, face issues including poor sensitivity, drug-induced changes in global brain perfusion and the effects of anesthesia. Here we demonstrate the proof-of-concept of a minimally-invasive fUS imaging approach to detect the acute cholinergic modulatory effects of Scopolamine (ScoP) on functional brain connectivity in awake and behaving mice, through the intact skull. A machine-learning algorithm constructed an ad-hoc pharmacological score from the ScoP-induced changes in connectivity patterns of five mice. The discrimination model shows important ScoP-induced increase of the hippocampo-cortical connectivity. The pharmacological score led to robust discrimination of ScoP treatment from baseline in an independent dataset and showed, in another independent group, dose-dependent specific effects of central cholinergic modulation of functional connectivity, independent from global brain perfusion changes. In conclusion, we introduce pharmaco-fUS as a simple, robust, specific and sensitive modality to monitor drug effects on perfusion and functional connectivity in the awake mouse brain.

Assessment of biopsy-proven liver fibrosis by two-dimensional shear wave elastography: An individual patient data-based meta-analysis.

Two-dimensional shear wave elastography (2D-SWE) has proven to be efficient for the evaluation of liver fibrosis in small to moderate-sized clinical trials. We aimed at running a larger-scale meta-analysis of individual data. Centers which have worked with Aixplorer ultrasound equipment were contacted to share their data. Retrospective statistical analysis used direct and paired receiver operating characteristic and area under the receiver operating characteristic curve (AUROC) analyses, accounting for random effects. Data on both 2D-SWE and liver biopsy were available for 1,134 patients from 13 sites, as well as on successful transient elastography in 665 patients. Most patients had chronic hepatitis C (n = 379), hepatitis B (n = 400), or nonalcoholic fatty liver disease (n = 156). AUROCs of 2D-SWE in patients with hepatitis C, hepatitis B, and nonalcoholic fatty liver disease were 86.3%, 90.6%, and 85.5% for diagnosing significant fibrosis and 92.9%, 95.5%, and 91.7% for diagnosing cirrhosis, respectively. The AUROC of 2D-SWE was 0.022-0.084 (95% confidence interval) larger than the AUROC of transient elastography for diagnosing significant fibrosis (P = 0.001) and 0.003-0.034 for diagnosing cirrhosis (P = 0.022) in all patients. This difference was strongest in hepatitis B patients. CONCLUSION: 2D-SWE has good to excellent performance for the noninvasive staging of liver fibrosis in patients with hepatitis B; further prospective studies are needed for head-to-head comparison between 2D-SWE and other imaging modalities to establish disease-specific appropriate cutoff points for assessment of fibrosis stage. (Hepatology 2018;67:260-272).

EEG and functional ultrasound imaging in mobile rats.

We developed an integrated experimental framework that extends the brain exploration capabilities of functional ultrasound imaging to awake and mobile rats. In addition to acquiring hemodynamic data, this method further allows parallel access to electroencephalography (EEG) recordings of neuronal activity. We illustrate this approach with two proofs of concept: a behavioral study on theta rhythm activation in a maze running task and a disease-related study on spontaneous epileptic seizures.

3D functional ultrasound imaging of the cerebral visual system in rodents.

3D functional imaging of the whole brain activity during visual task is a challenging task in rodents due to the complex tri-dimensional shape of involved brain regions and the fine spatial and temporal resolutions required to reveal the visual tract. By coupling functional ultrasound (fUS) imaging with a translational motorized stage and an episodic visual stimulation device, we managed to accurately map and to recover the activity of the visual cortices, the Superior Colliculus (SC) and the Lateral Geniculate Nuclei (LGN) in 3D. Cerebral Blood Volume (CBV) responses during visual stimuli were found to be highly correlated with the visual stimulus time profile in visual cortices (r=0.6), SC (r=0.7) and LGN (r=0.7). These responses were found dependent on flickering frequency and contrast, and optimal stimulus parameters for largest CBV increases were obtained. In particular, increasing the flickering frequency higher than 7Hz revealed a decrease of visual cortices response while the SC response was preserved. Finally, cross-correlation between CBV signals exhibited significant delays (d=0.35s +/-0.1s) between blood volume response in SC and visual cortices in response to our visual stimulus. These results emphasize the interest of fUS imaging as a whole brain neuroimaging modality for brain vision studies in rodent models.

4D microvascular imaging based on ultrafast Doppler tomography.

4D ultrasound microvascular imaging was demonstrated by applying ultrafast Doppler tomography (UFD-T) to the imaging of brain hemodynamics in rodents. In vivo real-time imaging of the rat brain was performed using ultrasonic plane wave transmissions at very high frame rates (18,000 frames per second). Such ultrafast frame rates allow for highly sensitive and wide-field-of-view 2D Doppler imaging of blood vessels far beyond conventional ultrasonography. Voxel anisotropy (100 µm × 100 µm × 500 µm) was corrected for by using a tomographic approach, which consisted of ultrafast acquisitions repeated for different imaging plane orientations over multiple cardiac cycles. UFT-D allows for 4D dynamic microvascular imaging of deep-seated vasculature (up to 20 mm) with a very high 4D resolution (respectively 100 µm × 100 µm × 100 µm and 10 ms) and high sensitivity to flow in small vessels (>1 mm/s) for a whole-brain imaging technique without requiring any contrast agent. 4D ultrasound microvascular imaging in vivo could become a valuable tool for the study of brain hemodynamics, such as cerebral flow autoregulation or vascular remodeling after ischemic stroke recovery, and, more generally, tumor vasculature response to therapeutic treatment.

Investigating liver stiffness and viscosity for fibrosis, steatosis and activity staging using shear wave elastography.

BACKGROUND & AIMS: Quantitative shear wave elastography was shown to be an effective tool for the non-invasive diagnosis and staging of chronic liver diseases. The liver shear modulus, estimated from the propagation velocity of shear waves, is correlated to the degree of fibrosis and can therefore be used for the non-invasive staging of fibrosis. METHODS: We performed a clinical prospective study in a total of 120 patients with various chronic liver diseases to compare the accuracy of supersonic shear imaging (SSI), a technique based on acoustic radiation and ultrafast ultrasound imaging, to 1D transient elastography (FibroScan) for the staging and grading of fibrosis as assessed by liver biopsy. Since shear wave propagation spectroscopy can also provide additional mechanical information on soft tissues, such as viscosity, we also investigated those new mechanical parameters as possible predictors of fibrosis, steatosis, and disease activity. RESULTS: SSI was successfully performed in 98.3% of patients and it was shown to be as accurate as FibroScan for the staging of fibrosis both for the whole population (N=120) and for the subgroup with viral hepatitis (n=70) (AUC=0.85 [0.77-0.96] and 0.89 [0.81-0.97] for significant fibrosis, AUC=0.90 [0.83-0.97] and 0.87 [0.75-0.98] for cirrhosis, with respect to SSI [n=68/70] and FibroScan [n=66/68]). Viscosity could also be used to stage the degree of fibrosis (AUC=0.76 [0.64-0.87] for significant fibrosis and AUC=0.87 [0.74-0.99] for cirrhosis), for the subgroup of patients with viral hepatitis (n=67/70) but was a poor predictor of disease activity and steatosis levels. CONCLUSIONS: Supersonic shear imaging is a robust technique for the staging of liver fibrosis. Liver viscosity was found to be correlated with fibrosis but not to steatosis or disease activity.

Curved toroidal row column addressed transducer for 3D ultrafast ultrasound imaging.

3D Imaging of the human heart at high frame rate is of major interest for various clinical applications. Electronic complexity and cost has prevented the dissemination of 3D ultrafast imaging into the clinic. Row column addressed (RCA) transducers provide volumetric imaging at ultrafast frame rate by using a low electronic channel count, but current models are ill-suited for transthoracic cardiac imaging due to field-of-view limitations. In this study, we proposed a mechanically curved RCA with an aperture adapted for transthoracic cardiac imaging (24 × 16 mm²). The RCA has a toroidal curved surface of 96 elements along columns (curvature radius rC = 4.47 cm) and 64 elements along rows (curvature radius rR = 3 cm). We implemented delay and sum beamforming with an analytical calculation of the propagation of a toroidal wave which was validated using simulations (Field II). The imaging performance was evaluated on a calibrated phantom. Experimental 3D imaging was achieved up to 12 cm deep with a total angular aperture of 30° for both lateral dimensions. The Contrast-to-Noise ratio increased by 12 dB from 2 to 128 virtual sources. Then, 3D Ultrasound Localization Microscopy (ULM) was characterized in a sub-wavelength tube diameter. Finally, 3D ULM was demonstrated on a perfused ex-vivo swine heart to image the coronary microcirculation.

Dynamics of cerebral blood volume during and after middle cerebral artery occlusion in rats - Comparison between ultrafast ultrasound and dynamic susceptibility contrast-enhanced MRI measurements.

Tomographic perfusion imaging techniques are integral to translational stroke research paradigms that advance our understanding of the disease. Functional ultrasound (fUS) is an emerging technique that informs on cerebral blood volume (CBV) through ultrasensitive Doppler and flow velocity (CBFv) through ultrafast localization microscopy. It is not known how experimental results compare with a classical CBV-probing technique such as dynamic susceptibility contrast-enhanced perfusion MRI (DSC-MRI). To that end, we assessed hemodynamics based on uUS (n = 6) or DSC-MRI (n = 7) before, during and up to three hours after 90-minute filament-induced middle cerebral artery occlusion (MCAO) in rats. Recanalization was followed by a brief hyperperfusion response, after which CBV and CBFv temporarily normalized but progressively declined after one hour in the lesion territory. DSC-MRI data corroborated the incomplete restoration of CBV after recanalization, which may have been caused by the free-breathing anesthetic regimen. During occlusion, MCAO-induced hypoperfusion was more discrepant between either technique, likely attributable to artefactual signal mechanisms related to slow flow, and processing algorithms employed for either technique. In vivo uUS- and DSC-MRI-derived measures of CBV enable serial whole-brain assessment of post-stroke hemodynamics, but readouts from both techniques need to be interpreted cautiously in situations of very low blood flow.

Decoding motor plans using a closed-loop ultrasonic brain-machine interface.

Brain-machine interfaces (BMIs) enable people living with chronic paralysis to control computers, robots and more with nothing but thought. Existing BMIs have trade-offs across invasiveness, performance, spatial coverage and spatiotemporal resolution. Functional ultrasound (fUS) neuroimaging is an emerging technology that balances these attributes and may complement existing BMI recording technologies. In this study, we use fUS to demonstrate a successful implementation of a closed-loop ultrasonic BMI. We streamed fUS data from the posterior parietal cortex of two rhesus macaque monkeys while they performed eye and hand movements. After training, the monkeys controlled up to eight movement directions using the BMI. We also developed a method for pretraining the BMI using data from previous sessions. This enabled immediate control on subsequent days, even those that occurred months apart, without requiring extensive recalibration. These findings establish the feasibility of ultrasonic BMIs, paving the way for a new class of less-invasive (epidural) interfaces that generalize across extended time periods and promise to restore function to people with neurological impairments.

Leptomeningeal collaterals regulate reperfusion in ischemic stroke and rescue the brain from futile recanalization.

Recanalization is the mainstay of ischemic stroke treatment. However, even with timely clot removal, many stroke patients recover poorly. Leptomeningeal collaterals (LMCs) are pial anastomotic vessels with yet-unknown functions. We applied laser speckle imaging, ultrafast ultrasound, and two-photon microscopy in a thrombin-based mouse model of stroke and fibrinolytic treatment to show that LMCs maintain cerebral autoregulation and allow for gradual reperfusion, resulting in small infarcts. In mice with poor LMCs, distal arterial segments collapse, and deleterious hyperemia causes hemorrhage and mortality after recanalization. In silico analyses confirm the relevance of LMCs for preserving perfusion in the ischemic region. Accordingly, in stroke patients with poor collaterals undergoing thrombectomy, rapid reperfusion resulted in hemorrhagic transformation and unfavorable recovery. Thus, we identify LMCs as key components regulating reperfusion and preventing futile recanalization after stroke. Future therapeutic interventions should aim to enhance collateral function, allowing for beneficial reperfusion after stroke.

PerceptFlow: Real-Time Ultrafast Doppler Image Enhancement Using Deep Convolutional Neural Network and Perceptual Loss.

Ultrafast ultrasound is an emerging imaging modality derived from standard medical ultrasound. It allows for a high spatial resolution of 100 µm and a temporal resolution in the millisecond range with techniques such as ultrafast Doppler imaging. Ultrafast Doppler imaging has become a priceless tool for neuroscience, especially for visualizing functional vascular structures and navigating the brain in real time. Yet, the quality of a Doppler image strongly depends on experimental conditions and is easily subject to artifacts and deterioration, especially with transcranial imaging, which often comes at the cost of higher noise and lower sensitivity to small blood vessels. A common solution to better visualize brain vasculature is either accumulating more information, integrating the image over several seconds or using standard filter-based enhancement techniques, which often over-smooth the image, thus failing both to preserve sharp details and to improve our perception of the vasculature. In this study we propose combining the standard Doppler accumulation process with a real-time enhancement strategy, based on deep-learning techniques, using perceptual loss (PerceptFlow). With our perceptual approach, we bypass the need for long integration times to enhance Doppler images. We applied and evaluated our proposed method on transcranial Doppler images of mouse brains, outperforming state-of-the-art filters. We found that, in comparison to standard filters such as the Gaussian filter (GF) and block-matching and 3-D filtering (BM3D), PerceptFlow was capable of reducing background noise with a significant increase in contrast and contrast-to-noise ratio, as well as better preserving details without compromising spatial resolution.

Co-variations of cerebral blood volume and single neurons discharge during resting state and visual cognitive tasks in non-human primates.

To better understand how the brain allows primates to perform various sets of tasks, the ability to simultaneously record neural activity at multiple spatiotemporal scales is challenging but necessary. However, the contribution of single-unit activities (SUAs) to neurovascular activity remains to be fully understood. Here, we combine functional ultrasound imaging of cerebral blood volume (CBV) and SUA recordings in visual and fronto-medial cortices of behaving macaques. We show that SUA provides a significant estimate of the neurovascular response below the typical fMRI spatial resolution of 2mm3. Furthermore, our results also show that SUAs and CBV activities are statistically uncorrelated during the resting state but correlate during tasks. These results have important implications for interpreting functional imaging findings while one constructs inferences of SUA during resting state or tasks.

In vivo whole brain microvascular imaging in mice using transcranial 3D Ultrasound Localization Microscopy.

BACKGROUND: Non-invasive high-resolution imaging of the cerebral vascular anatomy and function is key for the study of intracranial aneurysms, stenosis, arteriovenous malformations, and stroke, but also neurological pathologies, such as degenerative diseases. Direct visualization of the microvascular networks in the whole brain remains however challenging in vivo. METHODS: In this work, we performed 3D ultrafast ultrasound localization microscopy (ULM) using a 2D ultrasound matrix array and mapped the whole-brain microvasculature and flow at microscopic resolution in C57Bl6 mice in vivo. FINDINGS: We demonstrated that the mouse brain vasculature can be imaged directly through the intact skull at a spatial resolution of 20 µm and over the whole brain depth and at high temporal resolution (750 volumes.s-1). Individual microbubbles were tracked to estimate the flow velocities that ranged from 2 mm.s-1 in arterioles and venules up to 100 mm.s-1 in large vessels. The vascular maps were registered automatically with the Allen atlas in order to extract quantitative vascular parameters such as local flow rates and velocities in regions of interest. INTERPRETATION: We show the potential of 3D ULM to provide new insights into whole-brain vascular flow in mice models at unprecedented vascular scale for an in vivo technique. This technology is highly translational and has the potential to become a major tool for the clinical investigation of the cerebral microcirculation. FUNDING: This study was supported by the European Research Council under the European Union's Seventh Framework Program (FP/2007-2013) / ERC Grant Agreement n° 311025 and by the Fondation Bettencourt-Schueller under the program "Physics for Medicine". We acknowledge the ART (Technological Research Accelerator) biomedical ultrasound program of INSERM.

Ultrafast Doppler imaging and ultrasound localization microscopy reveal the complexity of vascular rearrangement in chronic spinal lesion.

Acute spinal cord injury (SCI) leads to severe damage to the microvascular network. The process of spontaneous repair is accompanied by formation of new blood vessels; their functionality, however, presumably very important for functional recovery, has never been clearly established, as most studies so far used fixed tissues. Here, combining ultrafast Doppler imaging and ultrasound localization microscopy (ULM) on the same animals, we proceeded at a detailed analysis of structural and functional vascular alterations associated with the establishment of chronic SCI, both at macroscopic and microscopic scales. Using a standardized animal model of SCI, our results demonstrate striking hemodynamic alterations in several subparts of the spinal cord: a reduced blood velocity in the lesion site, and an asymmetrical hypoperfusion caudal but not rostral to the lesion. In addition, the worsening of many evaluated parameters at later time points suggests that the neoformed vascular network is not yet fully operational, and reveals ULM as an efficient in vivo readout for spinal cord vascular alterations. Finally, we show statistical correlations between the diverse biomarkers of vascular dysfunction and SCI severity. The imaging modality developed here will allow evaluating recovery of vascular function over time in pre-clinical models of SCI. Also, used on SCI patients in combination with other quantitative markers of neural tissue damage, it may help classifying lesion severity and predict possible treatment outcomes in patients.