Inter-operator Reliability of Magnetic Resonance Image-Based Computational Fluid Dynamics Prediction of Cerebrospinal Fluid Motion in the Cervical Spine

Numerical analysis of the effect of Syringomyelia on cerebrospinal fluid dynamics

Spinal cord enlargement (SCE) includes conditions such as Syringomyelia, tumors, and tumor-like cases of demyelination, edema, or inflammation. These conditions involve fluid-filled cysts, known as syrinx, or masses of tissue, referred to as tumors, which cause increased pressure within the spinal cord (SC) and obstruct cerebrospinal fluid (CSF) circulation. To assess the impact of SCE location and diameter, we constructed fifteen computational SC models, each featuring a SCE placed in one of five probable locations with 20 %, 40 %, and 60 % stenosis. Our objective was to investigate how the location, diameter, and length of the SCE influence CSF velocity pattern and to identify the most critical location in the SC associated with this condition. The results indicated a velocity increase of 0.5 cm/(s) near models with 60 % stenosis. Importantly, SCE located from T1 to T5 exhibit a more pronounced reduction, exceeding 6.5, in the Womersley number. Our finding suggests that this region is the most vulnerable for SCE formation due to its significant impact on fluid circulation. The identification of specific locations within the SC associated with heightened risk can contribute to an improved understanding, treatment and management of SCE.

A computational fluid dynamics study to assess the impact of coughing on cerebrospinal fluid dynamics in Chiari type 1 malformation

Chiari type 1 malformation is a neurological disorder characterized by an obstruction of the cerebrospinal fluid (CSF) circulation between the brain (intracranial) and spinal cord (spinal) compartments. Actions such as coughing might evoke spinal cord complications in patients with Chiari type 1 malformation, but the underlying mechanisms are not well understood. More insight into the impact of the obstruction on local and overall CSF dynamics can help reveal these mechanisms. Therefore, our previously developed computational fluid dynamics framework was used to establish a subject-specific model of the intracranial and upper spinal CSF space of a healthy control. In this model, we emulated a single cough and introduced porous zones to model a posterior (OBS-1), mild (OBS-2), and severe posterior-anterior (OBS-3) obstruction. OBS-1 and OBS-2 induced minor changes to the overall CSF pressures, while OBS-3 caused significantly larger changes with a decoupling between the intracranial and spinal compartment. Coughing led to a peak in overall CSF pressure. During this peak, pressure differences between the lateral ventricles and the spinal compartment were locally amplified for all degrees of obstruction. These results emphasize the effects of coughing and indicate that severe levels of obstruction lead to distinct changes in intracranial pressure.

Translational approaches for brain delivery of biologics via cerebrospinal fluid

Delivery of biologics via cerebrospinal fluid (CSF) has demonstrated potential to access the tissues of the central nervous system (CNS) by circumventing the blood‐brain barrier and blood‐CSF barrier. Developing an effective CSF drug delivery strategy requires optimization of multiple parameters, including choice of CSF access point, delivery device technology, and delivery kinetics to achieve effective therapeutic concentrations in the target brain region, whereas also considering the biologic modality, mechanism of action, disease indication, and patient population. This review discusses key preclinical and clinical examples of CSF delivery for different biologic modalities (antibodies, nucleic acid‐based therapeutics, and gene therapy) to the brain via CSF or CNS access routes (intracerebroventricular, intrathecal‐cisterna magna, intrathecal‐lumbar, intraparenchymal, and intranasal), including the use of novel device technologies. This review also discusses quantitative models of CSF flow that provide insight into the effect of fluid dynamics in CSF on drug delivery and CNS distribution. Such models can facilitate delivery device design and pharmacokinetic/pharmacodynamic translation from preclinical species to humans in order to optimize CSF drug delivery to brain regions of interest.

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.

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 Common Data Elements for Use in Chiari Malformation Type I Clinical Research: An NIH/NINDS Project

The management of Chiari I malformation (CMI) is controversial because treatment methods vary and treatment decisions rest on incomplete understanding of its complex symptom patterns, etiologies, and natural history. Validity of studies that attempt to compare treatment of CMI has been limited because of variable terminology and methods used to describe study subjects. The goal of this project was to standardize terminology and methods by developing a comprehensive set of Common Data Elements (CDEs), data definitions, case report forms (CRFs), and outcome measure recommendations for use in CMI clinical research, as part of the CDE project at the National Institute of Neurological Disorders and Stroke (NINDS) of the US National Institutes of Health. A working group, comprising over 30 experts, developed and identified CDEs, template CRFs, data dictionaries, and guidelines to aid investigators starting and conducting CMI clinical research studies. The recommendations were compiled, internally reviewed, and posted online for external public comment. In October 2016, version 1.0 of the CMI CDE recommendations became available on the NINDS CDE website. The recommendations span these domains: Core Demographics/Epidemiology; Presentation/Symptoms; Co-Morbidities/Genetics; Imaging; Treatment; and Outcome. Widespread use of CDEs could facilitate CMI clinical research trial design, data sharing, retrospective analyses, and consistent data sharing between CMI investigators around the world. Updating of CDEs will be necessary to keep them relevant and applicable to evolving research goals for understanding CMI and its treatment.

Characterization of intrathecal cerebrospinal fluid geometry and dynamics in cynomolgus monkeys (macaca fascicularis) by magnetic resonance imaging

Recent advancements have been made toward understanding the diagnostic and therapeutic potential of cerebrospinal fluid (CSF) and related hydrodynamics. Increased understanding of CSF dynamics may lead to improved detection of central nervous system (CNS) diseases and optimized delivery of CSF based CNS therapeutics, with many proposed therapeutics hoping to successfully treat or cure debilitating neurological conditions. Before significant strides can be made toward the research and development of interventions designed for human use, additional research must be carried out with representative subjects such as non-human primates (NHP). This study presents a geometric and hydrodynamic characterization of CSF in eight cynomolgus monkeys (Macaca fascicularis) at baseline and two-week follow-up.

Results showed that CSF flow along the entire spine was laminar with a Reynolds number ranging up to 80 and average Womersley number ranging from 4.1–7.7. Maximum CSF flow rate occurred ~25 mm caudal to the foramen magnum. Peak CSF flow rate ranged from 0.3–0.6 ml/s at the C3-C4 level. Geometric analysis indicated that average intrathecal CSF volume below the foramen magnum was 7.4 ml. The average surface area of the spinal cord and dura was 44.7 and 66.7 cm² respectively. Subarachnoid space cross-sectional area and hydraulic diameter ranged from 7–75 mm² and 2–3.7 mm, respectively. Stroke volume had the greatest value of 0.14 ml at an axial location corresponding to C3-C4.

Cardiac-Related Spinal Cord Tissue Motion at the Foramen Magnum is Increased in Patients with Type I Chiari Malformation and Decreases Postdecompression Surgery

OBJECTIVE

Type 1 Chiari malformation (CM-I) is a craniospinal disorder historically defined by cerebellar tonsillar position greater than 3-5 mm below the foramen magnum (FM). This definition has come under question because quantitative measurements of cerebellar herniation do not always correspond with symptom severity. Researchers have proposed several additional radiographic diagnostic criteria based on dynamic motion of fluids and/or tissues. The present study objective was to determine if cardiac-related craniocaudal spinal cord tissue displacement is an accurate indicator of the presence of CM-I and if tissue displacement is altered with decompression.

METHODS

A cohort of 20 symptomatic patients underwent decompression surgery. Fifteen healthy volunteers were recruited for comparison with the CM-I group. Axial phase-contrast magnetic resonance imaging (PC-MRI) measurements were collected before and after surgery at the FM with cranial-caudal velocity encoding and 20 frames per cardiac cycle with retrospective reconstruction. Spinal cord motion (SCM) at the FM was quantified based on the peak-to-peak integral of average spinal cord velocity.

RESULTS

Tissue motion for the presurgical group was significantly greater than controls (P = 0.0009). Motion decreased after surgery (P = 0.058) with an effect size of -0.151 mm and a standard error of 0.066 mm. Postoperatively, no statistical difference from controls in bulk displacement at the FM was found (P = 0.200) after post hoc testing using the Tukey adjustment for multiple comparisons.

CONCLUSIONS

These results support SCM measurement by PC-MRI as a possible noninvasive radiographic diagnostic for CM-I. Dynamic measurement of SCM provides unique diagnostic information about CM-I alongside static quantification of tonsillar position and other intracranial morphometrics.

Nonuniform Moving Boundary Method for Computational Fluid Dynamics Simulation of Intrathecal Cerebrospinal Flow Distribution in a Cynomolgus Monkey

A detailed quantification and understanding of cerebrospinal fluid (CSF) dynamics may improve detection and treatment of central nervous system (CNS) diseases and help optimize CSF system-based delivery of CNS therapeutics. This study presents a computational fluid dynamics (CFD) model that utilizes a nonuniform moving boundary approach to accurately reproduce the nonuniform distribution of CSF flow along the spinal subarachnoid space (SAS) of a single cynomolgus monkey. A magnetic resonance imaging (MRI) protocol was developed and applied to quantify subject-specific CSF space geometry and flow and define the CFD domain and boundary conditions. An algorithm was implemented to reproduce the axial distribution of unsteady CSF flow by nonuniform deformation of the dura surface. Results showed that maximum difference between the MRI measurements and CFD simulation of CSF flow rates was <3.6%. CSF flow along the entire spine was laminar with a peak Reynolds number of ∼150 and average Womersley number of ∼5.4. Maximum CSF flow rate was present at the C4-C5 vertebral level. Deformation of the dura ranged up to a maximum of 134 μm. Geometric analysis indicated that total spinal CSF space volume was ∼8.7 ml. Average hydraulic diameter, wetted perimeter, and SAS area were 2.9 mm, 37.3 mm and 27.24 mm2, respectively. CSF pulse wave velocity (PWV) along the spine was quantified to be 1.2 m/s.

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.

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

GOAL

Develop and test an MRI-compatible hydrodynamic simulator of cerebrospinal fluid (CSF) motion in the cervical spinal subarachnoid space. Four anatomically realistic subject-specific models were created based on a 22-year-old healthy volunteer and a five-year-old patient diagnosed with Chiari I malformation.

METHODS

The in vitro models were based on manual segmentation of high-resolution T2-weighted MRI of the cervical spine. Anatomically realistic dorsal and ventral spinal cord nerve rootlets (NR) were added. Models were three dimensional (3-D) printed by stereolithography with 50-μm layer thickness. A computer controlled pump system was used to replicate the shape of the subject specific in vivo CSF flow measured by phase-contrast MRI. Each model was then scanned by T2-weighted and 4-D phase contrast MRI (4D flow).

RESULTS

Cross-sectional area, wetted perimeter, and hydraulic diameter were quantified for each model. The oscillatory CSF velocity field (flow jets near NR, velocity profile shape, and magnitude) had similar characteristics to previously reported studies in the literature measured by in vivo MRI.

CONCLUSION

This study describes the first MRI-compatible hydrodynamic simulator of CSF motion in the cervical spine with anatomically realistic NR. NR were found to impact CSF velocity profiles to a great degree.

SIGNIFICANCE

CSF hydrodynamics are thought to be altered in craniospinal disorders such as Chiari I malformation. MRI scanning techniques and protocols can be used to quantify CSF flow alterations in disease states. The provided in vitro models can be used to test the reliability of these protocols across MRI scanner manufacturers and machines.

Accuracy of 4D Flow Measurement of Cerebrospinal Fluid Dynamics in the Cervical Spine: An In Vitro Verification Against Numerical Simulation

Abnormal alterations in cerebrospinal fluid (CSF) flow are thought to play an important role in pathophysiology of various craniospinal disorders such as hydrocephalus and Chiari malformation. Three directional phase contrast MRI (4D Flow) has been proposed as one method for quantification of the CSF dynamics in healthy and disease states, but prior to further implementation of this technique, its accuracy in measuring CSF velocity magnitude and distribution must be evaluated. In this study, an MR-compatible experimental platform was developed based on an anatomically detailed 3D printed model of the cervical subarachnoid space and subject specific flow boundary conditions. Accuracy of 4D Flow measurements was assessed by comparison of CSF velocities obtained within the in vitro model with the numerically predicted velocities calculated from a spatially averaged computational fluid dynamics (CFD) model based on the same geometry and flow boundary conditions. Good agreement was observed between CFD and 4D Flow in terms of spatial distribution and peak magnitude of through-plane velocities with an average difference of 7.5 and 10.6% for peak systolic and diastolic velocities, respectively. Regression analysis showed lower accuracy of 4D Flow measurement at the timeframes corresponding to low CSF flow rate and poor correlation between CFD and 4D Flow in-plane velocities.