An MRI-Compatible Hydrodynamic Simulator of Cerebrospinal Fluid Motion in the Cervical Spine

Cerebro-spinal flow pattern in the cervical subarachnoid space of healthy volunteers: Influence of the spinal cord morphology

INTRODUCTION

Toward further cerebro-spinal flow quantification in clinical practice, this study aims at assessing the variations in the cerebro spinal fluid flow pattern associated with change in the morphology of the subarachnoid space of the cervical canal of healthy humans by developing a computational fluid dynamics model.

METHODS

3D T2-space MRI sequence images of the cervical spine were used to segment 11 cervical subarachnoid space. Model validation (time-step, mesh size, size and number of boundary layers, influences of parted inflow and inflow continuous velocity) was performed a 40-year-old patient-specific model. Simulations were performed using computational fluid dynamics approach simulating transient flow (Sparlart-Almaras turbulence model) with a mesh size of 0.6, 6 boundary layers of 0.05 mm, a time step of 20 ms simulated on 15 cycles. Distributions of components velocity and WSS were respectively analyzed within the subarachnoid space (intervertebral et intravertebral levels) and on dura and pia maters.

RESULTS

Mean values cerebro spinal fluid velocity in specific local slices of the canal range between 0.07 and 0.17 m.s-1 and 0.1 and 0.3 m.s-1 for maximum values. Maximum wall shear stress values vary between 0.1 and 0.5 Pa with higher value at the middle of the cervical spine on pia mater and at the lower part of the cervical spine on dura mater. Intra and inter-individual variations of the wall shear stress were highlighted significant correlation gwith compression ratio (r = 0.76), occupation ratio and cross section area of the spinal cord.

CONCLUSION

The inter-individual variability in term of subarachnoid canal morphology and spinal cord position influence the cerebro-spinal flow pattern, highlighting the significance of canal morphology investigation before surgery.

Four-dimensional flow MRI for quantitative assessment of cerebrospinal fluid dynamics: Status and opportunities

Neurological disorders can manifest with altered neurofluid dynamics in different compartments of the central nervous system. These include alterations in cerebral blood flow, cerebrospinal fluid (CSF) flow, and tissue biomechanics. Noninvasive quantitative assessment of neurofluid flow and tissue motion is feasible with phase contrast magnetic resonance imaging (PC MRI). While two-dimensional (2D) PC MRI is routinely utilized in research and clinical settings to assess flow dynamics through a single imaging slice, comprehensive neurofluid dynamic assessment can be limited or impractical. Recently, four-dimensional (4D) flow MRI (or time-resolved three-dimensional PC with three-directional velocity encoding) has emerged as a powerful extension of 2D PC, allowing for large volumetric coverage of fluid velocities at high spatiotemporal resolution within clinically reasonable scan times. Yet, most 4D flow studies have focused on blood flow imaging. Characterizing CSF flow dynamics with 4D flow (i.e., 4D CSF flow) is of high interest to understand normal brain and spine physiology, but also to study neurological disorders such as dysfunctional brain metabolite waste clearance, where CSF dynamics appear to play an important role. However, 4D CSF flow imaging is challenged by the long T1 time of CSF and slower velocities compared with blood flow, which can result in longer scan times from low flip angles and extended motion-sensitive gradients, hindering clinical adoption. In this work, we review the state of 4D CSF flow MRI including challenges, novel solutions from current research and ongoing needs, examples of clinical and research applications, and discuss an outlook on the future of 4D CSF flow.

Spinal cord motion and CSF flow in the cervical spine of 70 healthy participants

Pulsatile spinal cord and CSF velocities related to the cardiac cycle can be depicted by phase-contrast MRI. Among patients with spontaneous intracranial hypotension, we have recently described relevant differences compared with healthy controls in segment C2/C3. The method might be a promising tool to solve clinical and diagnostic ambiguities. Therefore, it is important to understand the physiological range and the effects of clinical and anatomical parameters in healthy volunteers. Within a prospective study, 3D T2-weighted MRI for spinal canal anatomy and cardiac-gated phase-contrast MRI adapted to CSF flow and spinal cord motion for time-resolved velocity data and derivatives were performed in 70 participants (age 20-79 years) in segments C2/C3 and C5/C6. Correlations were analyzed by multiple linear regression models; p < 0.01 was required to assume a significant impact of clinical or anatomical data quantified by the regression coefficient B. Data showed that in C2/C3, the CSF and spinal cord craniocaudal velocity ranges were 4.5 ± 0.9 and 0.55 ± 0.15 cm/s; the total displacements were 1.1 ± 0.3 and 0.07 ± 0.02 cm, respectively. The craniocaudal range of the CSF flow rate was 8.6 ± 2.4 mL/s; the CSF stroke volume was 2.1 ± 0.7 mL. In C5/C5, physiological narrowing of the spinal canal caused higher CSF velocity ranges and lower stroke volume (C5/C6 B = +1.64 cm/s, p < 0.001; B = -0.4 mL, p = 0.002, respectively). Aging correlated to lower spinal cord motion (e.g., B = -0.01 cm per 10 years of aging, p < 0.001). Increased diastolic blood pressure was associated with lower spinal cord motion and CSF flow parameters (e.g., C2/C3 CSF stroke volume B = -0.3 mL per 10 mmHg, p < 0.001). Males showed higher CSF flow and spinal cord motion (e.g., CSF stroke volume B = +0.5 mL, p < 0.001; total displacement spinal cord B = +0.016 cm, p = 0.002). We therefore propose to stratify data for age and sex and to adjust for diastolic blood pressure and segmental narrowing in future clinical studies.

Cerebrospinal Fluid Hydrodynamics in Chiari I Malformation and Syringomyelia: Modeling Pathophysiology

Anatomic MRI, MRI flow studies, and intraoperative ultrasonography demonstrate that the Chiari I malformation obstructs CSF pathways at the foramen magnum and prevents normal CSF movement through the foramen magnum. Impaired CSF displacement across the foramen magnum during the cardiac cycle increases pulsatile hindbrain motion, pressure transmission to the spinal subarachnoid space, and the amplitude of CSF subarachnoid pressure waves driving CSF into the spinal cord. Central canal septations in adults prevent syrinx formation by CSF directly transmitting its pressure wave from the fourth ventricle to the central canal.

Quality Control for 4D Flow MR Imaging

In recent years, 4D flow MRI has become increasingly important in clinical applications for the blood vessels in the whole body, heart, and cerebrospinal fluid. 4D flow MRI has advantages over 2D cine phase-contrast (PC) MRI in that any targeted area of interest can be analyzed post-hoc, but there are some factors to be considered, such as ensuring measurement accuracy, a long imaging time and post-processing complexity, and interobserver variability.Due to the partial volume phenomenon caused by low spatial and temporal resolutions, the accuracy of flow measurement in 4D flow MRI is reduced. For spatial resolution, it is recommended to include at least four voxels in the vessel of interest, and if possible, six voxels. In large vessels such as the aorta, large voxels can be secured and SNR can be maintained, but in small cerebral vessels, SNR is reduced, resulting in reduced accuracy. A temporal resolution of less than 40 ms is recommended. The velocity-to-noise ratio (VNR) of low-velocity blood flow is low, resulting in poor measurement accuracy. The use of dual velocity encoding (VENC) or multi-VENC is recommended to avoid velocity wrap around and to increase VNR. In order to maintain sufficient spatio-temporal resolution, a longer imaging time is required, leading to potential patient movement during examination and a corresponding decrease in measurement accuracy.For the clinical application of new technologies, including various acceleration techniques, in vitro and in vivo accuracy verification based on existing accuracy-validated 2D cine PC MRI and 4D flow MRI, as well as accuracy verification on the conservation of mass' principle, should be performed, and intraobserver repeatability, interobserver reproducibility, and test-retest reproducibility should be checked.

In vitro evaluation of cerebrospinal fluid velocity measurement in type I Chiari malformation: repeatability, reproducibility, and agreement using 2D phase contrast and 4D flow MRI

BACKGROUND

Phase contrast magnetic resonance imaging, PC MRI, is a valuable tool allowing for non-invasive quantification of CSF dynamics, but has lacked adoption in clinical practice for Chiari malformation diagnostics. To improve these diagnostic practices, a better understanding of PC MRI based measurement agreement, repeatability, and reproducibility of CSF dynamics is needed.

METHODS

An anatomically realistic in vitro subject specific model of a Chiari malformation patient was scanned three times at five different scanning centers using 2D PC MRI and 4D Flow techniques to quantify intra-scanner repeatability, inter-scanner reproducibility, and agreement between imaging modalities. Peak systolic CSF velocities were measured at nine axial planes using 2D PC MRI, which were then compared to 4D Flow peak systolic velocity measurements extracted at those exact axial positions along the model.

RESULTS

Comparison of measurement results showed good overall agreement of CSF velocity detection between 2D PC MRI and 4D Flow (p = 0.86), fair intra-scanner repeatability (confidence intervals ± 1.5 cm/s), and poor inter-scanner reproducibility. On average, 4D Flow measurements had a larger variability than 2D PC MRI measurements (standard deviations 1.83 and 1.04 cm/s, respectively).

CONCLUSION

Agreement, repeatability, and reproducibility of 2D PC MRI and 4D Flow detection of peak CSF velocities was quantified using a patient-specific in vitro model of Chiari malformation. In combination, the greatest factor leading to measurement inconsistency was determined to be a lack of reproducibility between different MRI centers. Overall, these findings may help lead to better understanding for application of 2D PC MRI and 4D Flow techniques as diagnostic tools for CSF dynamics quantification in Chiari malformation and related diseases.

Correlation of a new hydrodynamic index with other effective indexes in Chiari I malformation patients with different associations

This study aimed to find a new CSF hydrodynamic index to assess Chiari type I malformation (CM-I) patients’ conditions and examine the relationship of this new index with morphometric and volumetric changes in these patients and their clinical symptoms. To this end, 58 CM-I patients in four groups and 20 healthy subjects underwent PC-MRI. Ten morphometric and three volumetric parameters were calculated. The CSF hydrodynamic parameters were also analyzed through computational fluid dynamic (CFD) simulation. The maximum CSF pressure was identified as a new hydrodynamic parameter to assess the CM-I patients’ conditions. This parameter was similar in patients with the same symptoms regardless of the group to which they belonged. The result showed a weak correlation between the maximum CSF pressure and the morphometric parameters in the patients. Among the volumetric parameters, PCF volume had the highest correlation with the maximum CSF pressure, which its value being higher in patients with CM-I/SM/scoliosis (R² = 65.6%, P = 0.0022) than in the other patients. PCF volume was the more relevant volumetric parameter to assess the patients’ symptoms. The values of PCF volume were greater in patients that headache symptom was more obvious than other symptoms, as compared to the other patients.

Non-invasive MRI quantification of cerebrospinal fluid dynamics in amyotrophic lateral sclerosis patients

Background

Developing novel therapeutic agents to treat amyotrophic lateral sclerosis (ALS) has been difficult due to multifactorial pathophysiologic processes at work. Intrathecal drug administration shows promise due to close proximity of cerebrospinal fluid (CSF) to affected tissues. Development of effective intrathecal pharmaceuticals will rely on accurate models of how drugs are dispersed in the CSF. Therefore, a method to quantify these dynamics and a characterization of differences across disease states is needed.

Methods

Complete intrathecal 3D CSF geometry and CSF flow velocities at six axial locations in the spinal canal were collected by T2-weighted and phase-contrast MRI, respectively. Scans were completed for eight people with ALS and ten healthy controls. Manual segmentation of the spinal subarachnoid space was performed and coupled with an interpolated model of CSF flow within the spinal canal. Geometric and hydrodynamic parameters were then generated at 1 mm slice intervals along the entire spine. Temporal analysis of the waveform spectral content and feature points was also completed.

Results

Comparison of ALS and control groups revealed a reduction in CSF flow magnitude and increased flow propagation velocities in the ALS cohort. Other differences in spectral harmonic content and geometric comparisons may support an overall decrease in intrathecal compliance in the ALS group. Notably, there was a high degree of variability between cases, with one ALS patient displaying nearly zero CSF flow along the entire spinal canal.

Conclusion

While our sample size limits statistical confidence about the differences observed in this study, it was possible to measure and quantify inter-individual and cohort variability in a non-invasive manner. Our study also shows the potential for MRI based measurements of CSF geometry and flow to provide information about  the hydrodynamic environment of the spinal subarachnoid space. These dynamics may be studied further to understand the behavior of CSF solute transport in healthy and diseased states.

Development of the MRI Flow Phantom System Focused on Low Speed Flows in Fluid Machinery

At present, magnetic resonance imaging (MRI), which is one of the noninvasive diagnostic imaging techniques for medical use, has been applied to industrial product inspection. In this study, a hydraulic model test system, called the “flow phantom system,” which can generate steady flow and oscillating flow for time-resolved MRI, was developed to aid optimum design of fluid machinery. Simple flow phantoms are tested under laminar flow or oscillating flow to evaluate the validity and effectiveness of the proposed system. In this article, the test results of double cylindrical pipe flow, elbow pipe flow, cylindrical valve flow, and jet pump flow, which are often seen in fluid machines, are tested using the 2D time-spatial labeling inversion pulse (Time-SLIP) method, which can track a labeled water mass and visualize it using two-dimensional images. MRI-detected flow patterns were compared with particle image velocimetry (PIV) or numerical simulation.

A 3D subject-specific model of the spinal subarachnoid space with anatomically realistic ventral and dorsal spinal cord nerve rootlets

BACKGROUND

The spinal subarachnoid space (SSS) has a complex 3D fluid-filled geometry with multiple levels of anatomic complexity, the most salient features being the spinal cord and dorsal and ventral nerve rootlets. An accurate anthropomorphic representation of these features is needed for development of in vitro and numerical models of cerebrospinal fluid (CSF) dynamics that can be used to inform and optimize CSF-based therapeutics.

METHODS

A subject-specific 3D model of the SSS was constructed based on high-resolution anatomic MRI. An expert operator completed manual segmentation of the CSF space with detailed consideration of the anatomy. 31 pairs of semi-idealized dorsal and ventral nerve rootlets (NR) were added to the model based on anatomic reference to the magnetic resonance (MR) imaging and cadaveric measurements in the literature. Key design criteria for each NR pair included the radicular line, descending angle, number of NR, attachment location along the spinal cord and exit through the dura mater. Model simplification and smoothing was performed to produce a final model with minimum vertices while maintaining minimum error between the original segmentation and final design. Final model geometry and hydrodynamics were characterized in terms of axial distribution of Reynolds number, Womersley number, hydraulic diameter, cross-sectional area and perimeter.

RESULTS

The final model had a total of 139,901 vertices with a total CSF volume within the SSS of 97.3 cm³. Volume of the dura mater, spinal cord and NR was 123.1, 19.9 and 5.8 cm³. Surface area of these features was 318.52, 112.2 and 232.1 cm² respectively. Maximum Reynolds number was 174.9 and average Womersley number was 9.6, likely indicating presence of a laminar inertia-dominated oscillatory CSF flow field.

CONCLUSIONS

This study details an anatomically realistic anthropomorphic 3D model of the SSS based on high-resolution MR imaging of a healthy human adult female. The model is provided for re-use under the Creative Commons Attribution-ShareAlike 4.0 International license (CC BY-SA 4.0) and can be used as a tool for development of in vitro and numerical models of CSF dynamics for design and optimization of intrathecal therapeutics.