The influence of coughing on cerebrospinal fluid pressure in an in vitro syringomyelia model with spinal subarachnoid space stenosis

Water and brain function: effects of hydration status on neurostimulation with transcranial magnetic stimulation

Neurostimulation/neurorecording are tools to study, diagnose, and treat neurological/psychiatric conditions. Both techniques depend on volume conduction between scalp and excitable brain tissue. Here, we examine how neurostimulation with transcranial magnetic stimulation (TMS) is affected by hydration status, a physiological variable that can influence the volume of fluid spaces/cells, excitability, and cellular/global brain functioning. Normal healthy adult participants (32, 9 males) had common motor TMS measures taken in a repeated-measures design from dehydrated (12-h overnight fast/thirst) and rehydrated (identical dehydration protocol followed by rehydration with 1 L water in 1 h) testing days. The target region was left primary motor cortex hand area. Response at the target muscle was recorded with electromyography. Urinalysis confirmed hydration status. Motor hotspot shifted in half of participants. Motor threshold decreased in rehydration, indicating increased excitability. Even after redosing/relocalizing TMS to the new threshold/hotspot, rehydration still showed evidence of increased excitability: recruitment curve measures generally shifted upward and the glutamate-dependent paired-pulse protocol, short intracortical facilitation (SICF), was increased. Short intracortical inhibition (SICI), long intracortical inhibition (LICI), long intracortical facilitation (LICF), and cortical silent period (CSP) were relatively unaffected. The hydration perturbations were mild/subclinical based on the magnitude/speed and urinalysis. Motor TMS measures showed evidence of expected physiological changes of osmotic challenges. Rehydration showed signs of macroscopic and microscopic volume changes including decreased scalp-cortex distance (brain closer to stimulator) and astrocyte swelling-induced glutamate release. Hydration may be a source of variability affecting any techniques dependent on brain volumes/volume conduction. These concepts are important for researchers/clinicians using such techniques or dealing with the wide variety of disease processes involving water balance.NEW & NOTEWORTHY Hydration status can affect brain volumes and excitability, which should affect techniques dependent on electrical volume conduction, including neurostimulation/recording. We test the previously unknown effects of hydration on neurostimulation with TMS and briefly review relevant physiology of hydration. Rehydration showed lower motor threshold, shifted motor hotspot, and generally larger responses even after compensating for threshold/hotspot changes. This is important for clinical and research applications of neurostimulation/neurorecording and the many clinical disorders related to water balance.

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.

The prevalence of redundant nerve roots in standing positional MRI decreases by half in supine and almost to zero in flexed seated position: a retrospective cross-sectional cohort study

PURPOSE

This retrospective cross-sectional cohort study investigated the influence of posture on lordosis (LL), length of the spinal canal (LSC), anteroposterior diameter (APD L1-L5), dural cross-sectional area (DCSA) of the lumbar spinal canal, and the prevalence of redundant nerve roots (RNR) using positional magnetic resonance imaging (MRI) (0.6 T).

METHODS

Sixty-eight patients with single-level degenerative central lumbar spinal stenosis (cLSS) presenting with RNR in the standing position (STA) were also investigated in supine (SUP) or neutral seated (SIT) and flexed seated (FLEX) positions. Additionally, 45 patients complaining of back pain and without MRI evidence of LSS were evaluated. Statistical significance was set at p < 0.05.

RESULTS

Controls (A) and patients with cLSS (B) were comparable in terms of mean age (p = 0.88) and sex (p = 0.22). The progressive transition from STA to FLEX led to a comparable decrease in LL (p = 0.97), an increase in LSC (p = 0.80), and an increase in APD L1-L5 (p = 0.78). The APD of the stenotic level increased disproportionally between the different postures, up to 67% in FLEX compared to 29% in adjacent non-stenotic levels (p < 0.001). Therefore, the prevalence of RNR decreased to 49, 26, and 4% in SUP, SIT, and FLEX, respectively.

CONCLUSION

The prevalence of RNR in standing position was underestimated by half in supine position. Body postures modified LL, LSC, and APD similarly in patients and controls. Stenotic levels compensated for insufficient intraspinal volume with a disproportionate enlargement when switching from the STA to FLEX.

Immediate impact of yogic breathing on pulsatile cerebrospinal fluid dynamics

Cerebrospinal fluid (CSF), a clear fluid bathing the central nervous system (CNS), undergoes pulsatile movements. Together with interstitial fluid, CSF plays a critical role for the removal of waste products from the brain, and maintenance of the CNS health. As such, understanding the mechanisms driving CSF movement is of high scientific and clinical impact. Since pulsatile CSF dynamics is sensitive and synchronous to respiratory movements, we are interested in identifying potential integrative therapies such as yogic breathing to regulate CSF dynamics, which has not been reported before. Here, we investigated the pre-intervention baseline data from our ongoing randomized controlled trial, and examined the impact of four yogic breathing patterns: (i) slow, (ii) deep abdominal, (iii) deep diaphragmatic, and (iv) deep chest breathing with the last three together forming a yogic breathing called three-part breath. We utilized our previously established non-invasive real-time phase contrast magnetic resonance imaging approach using a 3T MRI instrument, computed and tested differences in single voxel CSF velocities (instantaneous, respiratory, cardiac 1st and 2nd harmonics) at the level of foramen magnum during spontaneous versus yogic breathing. In examinations of 18 healthy participants (eight females, ten males; mean age 34.9 ± 14 (SD) years; age range: 18-61 years), we observed immediate increase in cranially-directed velocities of instantaneous-CSF 16-28% and respiratory-CSF 60-118% during four breathing patterns compared to spontaneous breathing, with the greatest changes during deep abdominal breathing (28%, p = 0.0008, and 118%, p = 0.0001, respectively). Cardiac pulsation was the primary source of pulsatile CSF motion except during deep abdominal breathing, when there was a comparable contribution of respiratory and cardiac 1st harmonic power [0.59 ± 0.78], suggesting respiration can be the primary regulator of CSF depending on the individual differences in breathing techniques. Further work is needed to investigate the impact of sustained training yogic breathing on pulsatile CSF dynamics for CNS health.

Testing and validation of reciprocating positive displacement pump for benchtop pulsating flow model of cerebrospinal fluid production and other physiologic systems

BACKGROUND

The flow of physiologic fluids through organs and organs systems is an integral component of their function. The complex fluid dynamics in many organ systems are still not completely understood, and in-vivo measurements of flow rates and pressure provide a testament to the complexity of each flow system. Variability in in-vivo measurements and the lack of control over flow characteristics leave a lot to be desired for testing and evaluation of current modes of treatments as well as future innovations. In-vitro models are particularly ideal for studying neurological conditions such as hydrocephalus due to their complex pathophysiology and interactions with therapeutic measures. The following aims to present the reciprocating positive displacement pump, capable of inducing pulsating flow of a defined volume at a controlled beat rate and amplitude. While the other fluidic applications of the pump are currently under investigation, this study was focused on simulating the pulsating cerebrospinal fluid production across profiles with varying parameters.

METHODS

Pumps were manufactured using 3D printed and injection molded parts. The pumps were powered by an Arduino-based board and proprietary software that controls the linear motion of the pumps to achieve the specified output rate at the desired pulsation rate and amplitude. A range of 0.01 [Formula: see text] to 0.7 [Formula: see text] was tested to evaluate the versatility of the pumps. The accuracy and precision of the pumps' output were evaluated by obtaining a total of 150 one-minute weight measurements of degassed deionized water per output rate across 15 pump channels. In addition, nine experiments were performed to evaluate the pumps' control over pulsation rate and amplitude.

RESULTS

Volumetric analysis of a total of 1200 readings determined that the pumps achieved the target output volume rate with a mean absolute error of -0.001034283 [Formula: see text] across the specified domain. It was also determined that the pumps can maintain pulsatile flow at a user-specified beat rate and amplitude.

CONCLUSION

The validation of this reciprocating positive displacement pump system allows for the future validation of novel designs to components used to treat hydrocephalus and other physiologic models involving pulsatile flow. Based on the promising results of these experiments at simulating pulsatile CSF flow, a benchtop model of human CSF production and distribution could be achieved through the incorporation of a chamber system and a compliance component.

Hydrodynamic comparison of shunt and endoscopic third ventriculostomy in adult hydrocephalus using in vitro models and fluid-structure interaction simulation

OBJECTIVE

The comparison of the efficiency of shunt placement and endoscopic third ventriculostomy (ETV) in treating of adult hydrocephalus patients with various intensities and different obstruction intensities in the aqueduct of Sylvius (AS).

METHODS

In vitro models with separated ventricles were simulated and implemented for modeling shunt and ETV surgeries in one healthy subject and hydrocephalus patients with various intensities, as well as three different obstruction intensities in AS and under two cerebrospinal fluid (CSF) dynamic conditions. The fluid-structure interaction simulation was also carried out to validate in vitro results.

RESULTS

The efficiency of both methods in reducing the maximum CSF pressure in the subarachnoid space (MCPS) decreased by an increase in the patient's intensities. Contrary to shunting, the efficiency of ETV in reducing MCPS demonstrated a decline (8.3-16.4%) by an increase in obstruction levels in AS. Based on the findings, shunt efficiency in decreasing MCPS in patients with low intensity was more remarkable compared to ETV. However, ETV was more efficient than shunt in the patient with intracranial hypertension. Further, shunt placement and ETV led to a significant reduction in the amplitude of CSF pressure in the SAS (ACPS) in patients with sneezing, coughing, Valsalva maneuver, and exercising effects in contrast to other patients. Moreover, ACPS reduction was not related to the intensity of the disease in both treatment methods. In contrast to shunt, an increase in the obstruction level in AS led to a reduction in ACPS in ETV in both CSF dynamic conditions.

CONCLUSIONS

The noises from irregular disorders increased the discharging of CSF after shunt placement, and activities such as sneezing, coughing, Valsalva maneuvers, and exercising increased the risk of shunt overdrainage by 10.4~47.8%, especially in the patient with intracranial hypertension. Based on the proposed in vitro ETV and shunt models, an increase of head compliance was higher in ETV compared to the shunt. Eventually, an increase in the obstruction level of AS after ETV led to a decline in head compliance in contrast to shunt.

Ex-vivo quantification of ovine pia arachnoid complex biomechanical properties under uniaxial tension

Background

The pia arachnoid complex (PAC) is a cerebrospinal fluid-filled tissue conglomerate that surrounds the brain and spinal cord. Pia mater adheres directly to the surface of the brain while the arachnoid mater adheres to the deep surface of the dura mater. Collagen fibers, known as subarachnoid trabeculae (SAT) fibers, and microvascular structure lie intermediately to the pia and arachnoid meninges. Due to its structural role, alterations to the biomechanical properties of the PAC may change surface stress loading in traumatic brain injury (TBI) caused by sub-concussive hits. The aim of this study was to quantify the mechanical and morphological properties of ovine PAC.

Methods

Ovine brain samples (n = 10) were removed from the skull and tissue was harvested within 30 min post-mortem. To access the PAC, ovine skulls were split medially from the occipital region down the nasal bone on the superior and inferior aspects of the skull. A template was used to remove arachnoid samples from the left and right sides of the frontal and occipital regions of the brain. 10 ex-vivo samples were tested with uniaxial tension at 2 mm s⁻¹, average strain rate of 0.59 s⁻¹, until failure at < 5 h post extraction. The force and displacement data were acquired at 100 Hz. PAC tissue collagen fiber microstructure was characterized using second-harmonic generation (SHG) imaging on a subset of n = 4 stained tissue samples. To differentiate transverse blood vessels from SAT by visualization of cell nuclei and endothelial cells, samples were stained with DAPI and anti-von Willebrand Factor, respectively. The Mooney-Rivlin model for average stress–strain curve fit was used to model PAC material properties.

Results

The elastic modulus, ultimate stress, and ultimate strain were found to be 7.7 ± 3.0, 2.7 ± 0.76 MPa, and 0.60 ± 0.13, respectively. No statistical significance was found across brain dissection locations in terms of biomechanical properties. SHG images were post-processed to obtain average SAT fiber intersection density, concentration, porosity, tortuosity, segment length, orientation, radial counts, and diameter as 0.23, 26.14, 73.86%, 1.07 ± 0.28, 17.33 ± 15.25 µm, 84.66 ± 49.18°, 8.15%, 3.46 ± 1.62 µm, respectively.

Conclusion

For the sizes, strain, and strain rates tested, our results suggest that ovine PAC mechanical behavior is isotropic, and that the Mooney-Rivlin model is an appropriate curve-fitting constitutive equation for obtaining material parameters of PAC tissues.

Intracranial pressure in the American Alligator (Alligator mississippiensis): reptilian meninges and orthostatic gradients

The cranial meninges of reptiles differ from the more widely studied mammalian pattern in that the intraventricular and subarachnoid spaces are, at least partially, isolated. This study was undertaken to investigate the bulk flow of cerebrospinal fluid, and the resulting changes in intracranial pressure, in a common reptilian species. Intracranial pressure was measured using ocular ultrasonography and by surgically implanting pressure cannulae into the cranial subarachnoid space. The system was then challenged by: rotating the animal to create orthostatic gradients, perturbation of the vascular system, administration of epinephrine, and cephalic cutaneous heating. Pressure changes determined from the implanted catheters and through quantification of the optic nerve sheath were highly correlated and showed a significant linear relationship with orthostatic gradients. The catheter pressure responses were phasic, with an initial rapid response followed by a much slower response; each phase accounted for roughly half of the total pressure change. No significant relationship was found between intracranial pressure and either heart rate or blood flow. The focal application of heat and the administration of epinephrine both increased intracranial pressure, the latter influence being particularly pronounced.

Enhanced in vitro model of the CSF dynamics

BACKGROUND

Fluid dynamics of the craniospinal system are complex and still not completely understood. In vivo flow and pressure measurements of the cerebrospinal fluid (CSF) are limited. Whereas in silico modeling can be an adequate pathway for parameter studies, in vitro modeling of the craniospinal system is essential for testing and evaluation of therapeutic measures associated with innovative implants relating to, for example, normal pressure hydrocephalus and other fluid disorders. Previously-reported in vitro models focused on the investigation of only one hypothesis of the fluid dynamics rather than developing a modular set-up to allow changes in focus of the investigation. The aim of this study is to present an enhanced and validated in vitro model of the CSF system which enables the future embedding of implants, the validation of in silico models or phase-contrast magnetic resonance imaging (PC-MRI) measurements and a variety of sensitivity analyses regarding pathological behavior, such as reduced CSF compliances, higher resistances or altered blood dynamics.

METHODS

The in vitro model consists of a ventricular system which is connected via the aqueduct to the cranial and spinal subarachnoid spaces. Two compliance chambers are integrated to cushion the arteriovenous blood flow generated by a cam plate unit enabling the modeling of patient specific flow dynamics. The CSF dynamics are monitored using three cranial pressure sensors and a spinal ultrasound flow meter. Measurements of the in vitro spinal flow were compared to cervical flow data recorded with PC-MRI from nine healthy young volunteers, and pressure measurements were compared to the literature values reported for intracranial pressure (ICP) to validate the newly developed in vitro model.

RESULTS

The maximum spinal CSF flow recorded in the in vitro simulation was 133.60 ml/min in the caudal direction and 68.01 ml/min in the cranial direction, whereas the PC-MRI flow data of the subjects showed 122.82 ml/min in the caudal and 77.86 ml/min in the cranial direction. In addition, the mean ICP (in vitro) was 12.68 mmHg and the pressure wave amplitude, 4.86 mmHg, which is in the physiological range.

CONCLUSIONS

The in vitro pressure values were in the physiological range. The amplitudes of the flow results were in good agreement with PC-MRI data of young and healthy volunteers. However, the maximum cranial flow in the in vitro model occurred earlier than in the PC-MRI data, which might be due to a lack of an in vitro dynamic compliance. Implementing dynamic compliances and related sensitivity analyses are major aspects of our ongoing research.

Computer simulation of syringomyelia in dogs

Background

Syringomyelia is a pathological condition in which fluid-filled cavities (syringes) form and expand in the spinal cord. Syringomyelia is often linked with obstruction of the craniocervical junction and a Chiari malformation, which is similar in both humans and animals. Some brachycephalic toy breed dogs such as Cavalier King Charles Spaniels (CKCS) are particularly predisposed. The exact mechanism of the formation of syringomyelia is undetermined and consequently with the lack of clinical explanation, engineers and mathematicians have resorted to computer models to identify possible physical mechanisms that can lead to syringes. We developed a computer model of the spinal cavity of a CKCS suffering from a large syrinx. The model was excited at the cranial end to simulate the movement of the cerebrospinal fluid (CSF) and the spinal cord due to the shift of blood volume in the cranium related to the cardiac cycle. To simulate the normal condition, the movement was prescribed to the CSF. To simulate the pathological condition, the movement of CSF was blocked.

Results

For normal conditions the pressure in the SAS was approximately 400 Pa and the same applied to all stress components in the spinal cord. The stress was uniformly distributed along the length of the spinal cord. When the blockage between the cranial and spinal CSF spaces forced the cord to move with the cardiac cycle, shear and axial normal stresses in the cord increased significantly. The sites where the elevated stress was most pronounced coincided with the axial locations where the syringes typically form, but they were at the perimeter rather than in the central portion of the cord. This elevated stress originated from the bending of the cord at the locations where its curvature was high.

Conclusions

The results suggest that it is possible that repetitive stressing of the spinal cord caused by its exaggerated movement could be a cause for the formation of initial syringes. Further consideration of factors such as cord tethering and the difference in mechanical properties of white and grey matter is needed to fully explore this possibility.

Electronic supplementary material

The online version of this article (10.1186/s12917-018-1410-7) contains supplementary material, which is available to authorized users.

Vertical diplopia and oscillopsia due to midbrain keyhole aqueduct syndrome associated with severe cough

Purpose

Midline structural defects in the neural axis can give rise to neuro-ophthalmic symptoms. We report a rare case of keyhole aqueduct syndrome presenting after two years of severe cough due to gastroesophageal reflux disease.

Observations

A 58-year-old woman with a 2-year history of daily, severe cough presented to the neuro-ophthalmology clinic with progressive diplopia and oscillopsia. Examination revealed a 1-2 Hz down-beating nystagmus in primary gaze that worsened with left, right, and down gazes. Gaze evoked nystagmus and mild paresis were also seen with up gaze. There was an incomitant left hypertropia due to skew deviation that worsened with right and up gazes and improved with down gaze. She also had a right-sided ptosis and a 3 mm anisocoria not due to cranial nerve 3 paresis or Horner's syndrome. Brain magnetic resonance imaging showed a 1.5 mm × 11.7 mm × 6 mm midline cleft in the ventral midbrain communicating with the cerebral aqueduct, consistent with keyhole aqueduct syndrome. Her nystagmus and diplopia improved with oral acetazolamide treatment, at high doses of 2500-3000 mg per day.

Conclusions and importance

We report the first case of midbrain keyhole aqueduct syndrome with ocular motor and other neuro-ophthalmic manifestations associated with severe cough. Although her cough was effectively treated and intracranial pressure measurement was normal, her ophthalmic symptoms continued to progress, which is common in previous cases reported. Treatment with acetazolamide led to significant improvement, supporting the use of acetazolamide in this rare condition.

Translaminar pressure in Caucasian normal tension glaucoma patients

PURPOSE

The aim of this study was to examine the translaminar pressure difference (TLP) in Caucasian patients with normal tension glaucoma (NTG) and its possible impact onto the pathogenesis of NTG.

METHODS

Retrospective analysis of medical records of patients with open-angle glaucoma (OAG) in the period from 2005 to 2015 from the Ophthalmology Department, Cantonal Hospital Aarau, Switzerland. A total of 67 eyes of 38 patients (mean age 68.6 ± 11.3 years, 21 women and 17 men) fulfilled the diagnostic criteria of progressive NTG and underwent lumbar puncture (LP) during computer-assisted cisternography (CT - cisternography). The intraocular pressure (IOP) and lumbar cerebrospinal fluid pressure (CSF-p) were analysed and the TLP calculated. The TLP was compared with the mean defect (MD) of visual fields. Statistical analysis was performed with the one and two-tailed paired and unpaired t-test and the non-parametric Spearman correlation test.

RESULTS

The mean lumbar opening CSF-p measured 11.6 ± 3.7 mmHg. The mean IOP in the right eye measured 14.7 ± 2.4 mmHg, in the left eye 14.7 ± 2.5 mmHg. The calculated mean TLP was 3.0 ± 4.2 mmHg in the right and 3.3 ± 4.3 mmHg in the left eye. There was no significant correlation between TLP and the MD of visual fields in both eyes.

CONCLUSIONS

This study did not confirm either a lower lumbar CSF-p or increased TLP compared to previous retrospective and prospective studies. As cerebrospinal fluid (CSF) flow is not homogenous throughout all CSF spaces and CSF-p and IOP fluctuate, the current view on TLP needs modifications to improve its validity.

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.

Strain rate dependent behavior of the porcine spinal cord under transverse dynamic compression

The accurate description of the mechanical properties of spinal cord tissue benefits to clinical evaluation of spinal cord injuries and is a required input for analysis tools such as finite element models. Unfortunately, available data in the literature generally relate mechanical properties of the spinal cord under quasi-static loading conditions, which is not adapted to the study of traumatic behavior, as neurological tissue adopts a viscoelastic behavior. Thus, the objective of this study is to describe mechanical properties of the spinal cord up to mechanical damage, under dynamic loading conditions. A total of 192 porcine cervical to lumbar spinal cord samples were compressed in a transverse direction. Loading conditions included ramp tests at 0.5, 5 or 50 s−1 and cyclic loading at 1, 10 or 20 Hz. Results showed that spinal cord behavior was significantly influenced by strain rate. Mechanical damage occurred at 0.64, 0.68 and 0.73 strains for 0.5, 5 or 50 s−1 loadings, respectively. Variations of behavior between the tested strain rates were explained by cyclic loading results, which revealed behavior more or less viscous depending on strain rate. Also, a parameter (stress multiplication factor) was introduced to allow transcription of a stress–strain behavior curve to different strain rates. This factor was described and was significantly different for cervical, thoracic and lumbar vertebral heights, and for the strain rates evaluated in this study.

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

For the first time, inter-operator dependence of MRI based computational fluid dynamics (CFD) modeling of cerebrospinal fluid (CSF) in the cervical spinal subarachnoid space (SSS) is evaluated. In vivo MRI flow measurements and anatomy MRI images were obtained at the cervico-medullary junction of a healthy subject and a Chiari I malformation patient. 3D anatomies of the SSS were reconstructed by manual segmentation by four independent operators for both cases. CFD results were compared at nine axial locations along the SSS in terms of hydrodynamic and geometric parameters. Intraclass correlation (ICC) assessed the inter-operator agreement for each parameter over the axial locations and coefficient of variance (CV) compared the percentage of variance for each parameter between the operators. Greater operator dependence was found for the patient (0.19 < ICC < 0.99) near the craniovertebral junction compared to the healthy subject (ICC > 0.78). For the healthy subject, hydraulic diameter and Womersley number had the least variance (CV = ~2%). For the patient, peak diastolic velocity and Reynolds number had the smallest variance (CV = ~3%). These results show a high degree of inter-operator reliability for MRI-based CFD simulations of CSF flow in the cervical spine for healthy subjects and a lower degree of reliability for patients with Type I Chiari malformation.

Spontaneous Improvement of Chiari I Malformation and Syringomyelia in a Patient With Cystic Fibrosis

BACKGROUND AND IMPORTANCE:

Syringomyelia is highly associated with Chiari I malformation, but the pathophysiologic mechanism of syrinx formation and its relation to downward cerebellar tonsillar displacement remains elusive. Cough, Valsalva maneuver, and other physiological strains transiently exacerbate the clinical symptoms of these conditions, exert profound effects on the flow dynamics across the craniospinal junction, and are thought to play an important role in the pathogenesis of syringomyelia.

CLINICAL PRESENTATION:

We report the case of a patient with cystic fibrosis who presented during an exacerbation of bronchiectasis and was found to have a Chiari I malformation with associated syringomyelia. Eight months later, when the patient had returned to baseline pulmonary status, repeat imaging showed interval improvement in both the size of the syrinx and descent of cerebellar tonsils.

CONCLUSION:

This rare case of spontaneous improvement of syringomyelia and Chiari I malformation attributable to relief from chronic cough offers interesting insight into the mechanism of these disorders.

Poro-elastic modeling of Syringomyelia - a systematic study of the effects of pia mater, central canal, median fissure, white and gray matter on pressure wave propagation and fluid movement within the cervical spinal cord

Syringomyelia, fluid-filled cavities within the spinal cord, occurs frequently in association with a Chiari I malformation and produces some of its most severe neurological symptoms. The exact mechanism causing syringomyelia remains unknown. Since syringomyelia occurs frequently in association with obstructed cerebrospinal fluid (CSF) flow, it has been hypothesized that syrinx formation is mechanically driven. In this study we model the spinal cord tissue either as a poro-elastic medium or as a solid linear elastic medium, and simulate the propagation of pressure waves through an anatomically plausible 3D geometry, with boundary conditions based on in vivo CSF pressure measurements. Then various anatomic and tissue properties are modified, resulting in a total of 11 variations of the model that are compared. The results show that an open segment of the central canal and a stiff pia (relative to the cord) both increase the radial pressure gradients and enhance interstitial fluid flow in the central canal. The anterior median fissure, anisotropic permeability of the white matter, and Poisson ratio play minor roles.

Syringomyelia in an older patient

We describe the case of an 80-year-old man who presented with lower limb upper motor neurone weakness and spinothalamic tract sensory deficit secondary to previously undiagnosed syringomyelia. The case highlights the need for methodical history, examination and investigation in elderly patients to achieve diagnostic accuracy.

The impact of spinal cord nerve roots and denticulate ligaments on cerebrospinal fluid dynamics in the cervical spine

Cerebrospinal fluid (CSF) dynamics in the spinal subarachnoid space (SSS) have been thought to play an important pathophysiological role in syringomyelia, Chiari I malformation (CM), and a role in intrathecal drug delivery. Yet, the impact that fine anatomical structures, including nerve roots and denticulate ligaments (NRDL), have on SSS CSF dynamics is not clear. In the present study we assessed the impact of NRDL on CSF dynamics in the cervical SSS. The 3D geometry of the cervical SSS was reconstructed based on manual segmentation of MRI images of a healthy volunteer and a patient with CM. Idealized NRDL were designed and added to each of the geometries based on in vivo measurments in the literature and confirmation by a neuroanatomist. CFD simulations were performed for the healthy and patient case with and without NRDL included. Our results showed that the NRDL had an important impact on CSF dynamics in terms of velocity field and flow patterns. However, pressure distribution was not altered greatly although the NRDL cases required a larger pressure gradient to maintain the same flow. Also, the NRDL did not alter CSF dynamics to a great degree in the SSS from the foramen magnum to the C1 level for the healthy subject and CM patient with mild tonsillar herniation (∼6 mm). Overall, the NRDL increased fluid mixing phenomena and resulted in a more complex flow field. Comparison of the streamlines of CSF flow revealed that the presence of NRDL lead to the formation of vortical structures and remarkably increased the local mixing of the CSF throughout the SSS.

Effect of the central canal in the spinal cord on fluid movement within the cord

Computational studies are used to demonstrate the effect of oscillating CSF flow on pressures within the spinal cord. We tested the hypothesis that the central canal in the spinal cord affects spinal cord pressure gradients resulting from oscillatory CSF flow. Two computational models of the spinal cord were created with the same dimensions. Model 1 had a homogeneous porous structure. Model 2 had the same structure with the addition of a central fluid filled space, representing the central canal of the cord. We simulated oscillatory flow in the fluid space using standard computational fluid dynamics tools. For all phases of the CSF flow cycle and for specific projections through the model we calculated pressure gradients and fluid movement in the cord models. Pressures in the models varied through the flow cycle. Model 1 had linearly varying pressure along its long axis that varied with the cycle and had no pressure gradients across the cord. Model 2 had nonlinear varying pressure along its long axis varying with the time in the cycle and had transient centrifugal and centripetal pressure gradients with a central fluid space. The radial pressures varied linearly with distance from the fluid space. Centrifugal and centripetal pressure gradients resulted in radially directed fluid flow in the cord. The central canal within the spinal cord alters the pressure fields occurring during oscillatory CSF flow and creates centrifugal and centripetal fluid flux in the cord.

An update on the pathogenesis of syringomyelia secondary to Chiari-like malformations in dogs

Syringomyelia (SM) is a spinal cord disease that can cause neuropathic pain in dogs. The pathogenesis of SM secondary to Chiari-like malformation (CM) has been the focus of intense research in recent years. The gulf in our understanding of CM/SM in dogs relative to the analogous human condition has progressively narrowed. CM is primarily a disease of abnormal geometric morphometry affecting the caudal cranial fossa and the brain parenchyma contained within it. This review describes how advanced imaging techniques have revealed a series of morphometric abnormalities associated with CM/SM. The series is presented in a logical order to help describe the pathogenesis of CM and the subsequent formation of syringes, with particular reference to the concepts of craniospinal compliance and cerebrospinal fluid pulse pressure timing.

Computational fluid dynamics modelling of cerebrospinal fluid pressure in Chiari malformation and syringomyelia

The pathogenesis of syringomyelia in association with Chiari malformation (CM) is unclear. Studies of patients with CM have shown alterations in the CSF velocity profile and these could contribute to syrinx development or enlargement. Few studies have considered the fluid mechanics of CM patients with and without syringomyelia separately. Three subject-specific CFD models were developed for a normal participant, a CM patient with syringomyelia and a CM patient without syringomyelia. Model geometries, CSF flow rate data and CSF velocity validation data were collected from MRI scans of the 3 subjects. The predicted peak CSF pressure was compared for the 3 models. An extension of the study performed geometry and flow substitution to investigate the relative effects of anatomy and CSF flow profile on resulting spinal CSF pressure. Based on 50 monitoring locations for each of the models, the CM models had significantly higher magnitude (p<0.01) peak CSF pressure compared with normal. When using the same CSF input flow waveform, changing the upper spinal geometry changed the magnitude of the CSF pressure gradient, and when using the same upper spinal geometry, changing the input flow waveform changed the timing of the peak pressure. This study may assist in understanding syringomyelia mechanisms and relative effects of CSF velocity profile and spinal geometry on CSF pressure.

Comparison of 4D phase-contrast MRI flow measurements to computational fluid dynamics simulations of cerebrospinal fluid motion in the cervical spine

Cerebrospinal fluid (CSF) dynamics in the cervical spinal subarachnoid space (SSS) have been thought to be important to help diagnose and assess craniospinal disorders such as Chiari I malformation (CM). In this study we obtained time-resolved three directional velocity encoded phase-contrast MRI (4D PC MRI) in three healthy volunteers and four CM patients and compared the 4D PC MRI measurements to subject-specific 3D computational fluid dynamics (CFD) simulations. The CFD simulations considered the geometry to be rigid-walled and did not include small anatomical structures such as nerve roots, denticulate ligaments and arachnoid trabeculae. Results were compared at nine axial planes along the cervical SSS in terms of peak CSF velocities in both the cranial and caudal direction and visual interpretation of thru-plane velocity profiles. 4D PC MRI peak CSF velocities were consistently greater than the CFD peak velocities and these differences were more pronounced in CM patients than in healthy subjects. In the upper cervical SSS of CM patients the 4D PC MRI quantified stronger fluid jets than the CFD. Visual interpretation of the 4D PC MRI thru-plane velocity profiles showed greater pulsatile movement of CSF in the anterior SSS in comparison to the posterior and reduction in local CSF velocities near nerve roots. CFD velocity profiles were relatively uniform around the spinal cord for all subjects. This study represents the first comparison of 4D PC MRI measurements to CFD of CSF flow in the cervical SSS. The results highlight the utility of 4D PC MRI for evaluation of complex CSF dynamics and the need for improvement of CFD methodology. Future studies are needed to investigate whether integration of fine anatomical structures and gross motion of the brain and/or spinal cord into the computational model will lead to a better agreement between the two techniques.

Syrinx fluid transport: modeling pressure-wave-induced flux across the spinal pial membrane

Syrinxes are fluid-filled cavities of the spinal cord that characterize syringomyelia, a disease involving neurological damage. Their formation and expansion is poorly understood, which has hindered successful treatment. Syrinx cavities are hydraulically connected with the spinal subarachnoid space (SSS) enveloping the spinal cord via the cord interstitium and the network of perivascular spaces (PVSs), which surround blood vessels penetrating the pial membrane that is adherent to the cord surface. Since the spinal canal supports pressure wave propagation, it has been hypothesized that wave-induced fluid exchange across the pial membrane may play a role in syrinx filling. To investigate this conjecture a pair of one-dimensional (1-d) analytical models were developed from classical elastic tube theory coupled with Darcy's law for either perivascular or interstitial flow. The results show that transpial flux serves as a mechanism for damping pressure waves by alleviating hoop stress in the pial membrane. The timescale ratio over which viscous and inertial forces compete was explicitly determined, which predicts that dilated PVS, SSS flow obstructions, and a stiffer and thicker pial membrane-all associated with syringomyelia-will increase transpial flux and retard wave travel. It was also revealed that the propagation of a pressure wave is aided by a less-permeable pial membrane and, in contrast, by a more-permeable spinal cord. This is the first modeling of the spinal canal to include both pressure-wave propagation along the spinal axis and a pathway for fluid to enter and leave the cord, which provides an analytical foundation from which to approach the full poroelastic problem.

A one-dimensional model of the spinal cerebrospinal-fluid compartment

Modeling of the cerebrospinal fluid (CSF) system in the spine is strongly motivated by the need to understand the origins of pathological conditions such as the emergence and growth of fluid-filled cysts in the spinal cord. In this study, a one-dimensional (1D) approximation for the flow in elastic conduits was used to formulate a model of the spinal CSF compartment. The modeling was based around a coaxial geometry in which the inner elastic cylinder represented the spinal cord, middle elastic tube represented the dura, and the outermost tube represented the vertebral column. The fluid-filled annuli between the cord and dura, and the dura and vertebral column, represented the subarachnoid and epidural spaces, respectively. The system of governing equations was constructed by applying a 1D form of mass and momentum conservation to all segments of the model. The developed 1D model was used to simulate CSF pulse excited by pressure disturbances in the subarachnoid and epidural spaces. The results were compared to those obtained from an equivalent two-dimensional finite element (FE) model which was implemented using a commercial software package. The analysis of linearized governing equations revealed the existence of three types of waves, of which the two slower waves can be clearly related to the wave modes identified in previous similar studies. The third, much faster, wave emanates directly from the vertebral column and has little effect on the deformation of the spinal cord. The results obtained from the 1D model and its FE counterpart were found to be in good general agreement even when sharp spatial gradients of the spinal cord stiffness were included; both models predicted large radial displacements of the cord at the location of an initial cyst. This study suggests that 1D modeling, which is computationally inexpensive and amenable to coupling with the models of the cranial CSF system, should be a useful approach for the analysis of some aspects of the CSF dynamics in the spine. The simulation of the CSF pulse excited by a pressure disturbance in the epidural space, points to the possibility that regions of the spinal cord with abnormally low stiffness may be prone to experiencing large strains due to coughing and sneezing.

A coupled hydrodynamic model of the cardiovascular and cerebrospinal fluid system

Coupling of the cardiovascular and cerebrospinal fluid (CSF) system is considered to be important to understand the pathophysiology of cerebrovascular and craniospinal disease and intrathecal drug delivery. A coupled cardiovascular and CSF system model was designed to examine the relation of spinal cord (SC) blood flow (SCBF) and CSF pulsations along the spinal subarachnoid space (SSS). A one-dimensional (1-D) cardiovascular tree model was constructed including a simplified SC arterial network. Connection between the cardiovascular and CSF system was accomplished by a transfer function based on in vivo measurements of CSF and cerebral blood flow. A 1-D tube model of the SSS was constructed based on in vivo measurements in the literature. Pressure and flow throughout the cardiovascular and CSF system were determined for different values of craniospinal compliance. SCBF results indicated that the cervical, thoracic, and lumbar SC each had a signature waveform shape. The cerebral blood flow to CSF transfer function reproduced an in vivo-like CSF flow waveform. The 1-D tube model of the SSS resulted in a distribution of CSF pressure and flow and a wave speed that were similar to those in vivo. The SCBF to CSF pulse delay was found to vary a great degree along the spine depending on craniospinal compliance and vascular anatomy. The properties and anatomy of the SC arterial network and SSS were found to have an important impact on pressure and flow and perivascular fluid movement to the SC. Overall, the coupled model provides predictions about the flow and pressure environment in the SC and SSS. More detailed measurements are needed to fully validate the model.

Cerebrospinal fluid hydrodynamics in type I Chiari malformation

PURPOSE

The objective of this study was to review past studies that have used engineering analysis to examine cerebrospinal fluid hydrodynamics in cranial and spinal subarachnoid spaces in both healthy humans and those affected by type I Chiari malformation.

METHODS

A PubMed search of literature pertaining to cerebrospinal fluid hydrodynamics was performed with a particular focus on those that utilized methods such as computational fluid dynamics or experimental flow modeling.

DISCUSSION

From the engineer's perspective, type I Chiari malformation is an abnormal geometry of the cerebellum that causes increased resistance to cerebrospinal fluid flow between the intracranial and spinal subarachnoid space. As such, understanding the hydrodynamics of cerebrospinal fluid in the craniospinal subarachnoid space has long been thought to be important in the diagnosis and management of type I Chiari malformation. Hydrodynamic quantification of cerebrospinal fluid motion in the subarachnoid space may better reflect the pathophysiology of the disorder and serve as a prognostic indicator in conjunction with geometric magnetic resonance measurements that are currently used clinically. This review discusses the results of studies that have sought to quantify the hydrodynamics of cerebrospinal fluid motion using computational and experimental modeling and critiques the methods by which the results were obtained.

CONCLUSION

Researchers have found differences in cerebrospinal fluid velocities and pressures in type I Chiari malformation patients compared to healthy subjects. However, further research is necessary to determine the causal relationship between changes to hydrodynamic parameters such as cerebrospinal fluid velocity, pressure, resistance to flow, and craniospinal compliance and clinical aspects such as neurological symptoms, radiological evidence of severity, and surgical success.

A lumped-parameter model of the cerebrospinal system for investigating arterial-driven flow in posttraumatic syringomyelia

Fluid transport in syringomyelia has remained enigmatic ever since the disease was first identified some three centuries ago. However, accumulating evidence in the last decade from animal studies implicates arterial pulsations in syrinx formation. In particular, it has been suggested that a phase difference between the pressure pulse in the spinal subarachnoid space and the perivascular spaces, due to a pathologically disturbed cerebrospinal fluid (CSF) or blood supply, could result in a net influx of CSF into the spinal cord (SC). A lumped-parameter model is developed of the cerebrospinal system to investigate this conjecture. It is found that although this phase-lag mechanism may operate, it requires the SC to have an intrinsic storage capacity due to the collapsibility of the contained venous reservoir. This net flux is associated with a higher mean pressure in the SC than the SSS which is maintained in the periodic steady state. According to our simulations the mechanical perturbations of arachnoiditis exacerbate the phase-lag effect, which may be partially alleviated by the presence of a posttraumatic syrinx and more completely by a syringo-subarachnoid shunt.