Cerebrospinal fluid dynamics in the human cranial subarachnoid space: an overlooked mediator of cerebral disease. I. Computational model

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.

Large-scale in-silico analysis of CSF dynamics within the subarachnoid space of the optic nerve

BACKGROUND

Impaired cerebrospinal fluid (CSF) dynamics is involved in the pathophysiology of neurodegenerative diseases of the central nervous system and the optic nerve (ON), including Alzheimer's and Parkinson's disease, as well as frontotemporal dementia. The smallness and intricate architecture of the optic nerve subarachnoid space (ONSAS) hamper accurate measurements of CSF dynamics in this space, and effects of geometrical changes due to pathophysiological processes remain unclear. The aim of this study is to investigate CSF dynamics and its response to structural alterations of the ONSAS, from first principles, with supercomputers.

METHODS

Large-scale in-silico investigations were performed by means of computational fluid dynamics (CFD) analysis. High-order direct numerical simulations (DNS) have been carried out on ONSAS geometry at a resolution of 1.625 μm/pixel. Morphological changes on the ONSAS microstructure have been examined in relation to CSF pressure gradient (CSFPG) and wall strain rate, a quantitative proxy for mass transfer of solutes.

RESULTS

A physiological flow speed of 0.5 mm/s is achieved by imposing a hydrostatic pressure gradient of 0.37-0.67 Pa/mm across the ONSAS structure. At constant volumetric rate, the relationship between pressure gradient and CSF-accessible volume is well captured by an exponential curve. The ONSAS microstructure exhibits superior mass transfer compared to other geometrical shapes considered. An ONSAS featuring no microstructure displays a threefold smaller surface area, and a 17-fold decrease in mass transfer rate. Moreover, ONSAS trabeculae seem key players in mass transfer.

CONCLUSIONS

The present analysis suggests that a pressure drop of 0.1-0.2 mmHg over 4 cm is sufficient to steadily drive CSF through the entire subarachnoid space. Despite low hydraulic resistance, great heterogeneity in flow speeds puts certain areas of the ONSAS at risk of stagnation. Alterations of the ONSAS architecture aimed at mimicking pathological conditions highlight direct relationships between CSF volume and drainage capability. Compared to the morphological manipulations considered herein, the original ONSAS architecture seems optimized towards providing maximum mass transfer across a wide range of pressure gradients and volumetric rates, with emphasis on trabecular structures. This might shed light on pathophysiological processes leading to damage associated with insufficient CSF flow in patients with optic nerve compartment syndrome.

A mechanistic pharmacokinetic model for intrathecal administration of antisense oligonucleotides

Intrathecal administration is an important mode for delivering biological agents targeting central nervous system (CNS) diseases. However, current clinical practices lack a sound theorical basis for a quantitative understanding of the variables and conditions that govern the delivery efficiency and specific tissue targeting especially in the brain. This work presents a distributed mechanistic pharmacokinetic model (DMPK) for predictive analysis of intrathecal drug delivery to CNS. The proposed DMPK model captures the spatiotemporal dispersion of antisense oligonucleotides (ASO) along the neuraxis over clinically relevant time scales of days and weeks as a function of infusion, physiological and molecular properties. We demonstrate its prediction capability using biodistribution data of antisense oligonucleotide (ASO) administration in non-human primates. The results are in close agreement with the observed ASO pharmacokinetics in all key compartments of the central nervous system. The model enables determination of optimal injection parameters such as intrathecal infusion volume and duration for maximum ASO delivery to the brain. Our quantitative model-guided analysis is suitable for identifying optimal parameter settings to target specific brain regions with therapeutic drugs such as ASOs.

Effect of cilia-induced surface velocity on cerebrospinal fluid exchange in the lateral ventricles

Ciliary motility disorders are known to cause hydrocephalus. The instantaneous velocity of cerebrospinal fluid (CSF) flow is dominated by artery pulsation, and it remains unclear why ciliary dysfunction results in hydrocephalus. In this study, we investigated the effects of cilia-induced surface velocity on CSF flow using computational fluid dynamics. A geometric model of the human ventricles was constructed using medical imaging data. The CSF produced by the choroid plexus and cilia-induced surface velocity were given as the velocity boundary conditions at the ventricular walls. We developed healthy and reduced cilia motility models based on experimental data of cilia-induced velocity in healthy wild-type and Dpcd-knockout mice. The results indicate that there is almost no difference in intraventricular pressure between healthy and reduced cilia motility models. Additionally, it was found that newly produced CSF from the choroid plexus did not spread to the anterior and inferior horns of the lateral ventricles in the reduced cilia motility model. These findings suggest that a ciliary motility disorder could delay CSF exchange in the anterior and inferior horns of the lateral ventricles.

Intrathecal delivery and its applications in leptomeningeal disease

Intrathecal delivery (IT) of opiates into the cerebrospinal fluid (CSF) for anesthesia and pain relief has been used clinically for decades, but this relatively straightforward approach of bypassing the blood-brain barrier has been underutilized for other indications because of its lack of utility in delivering small lipid-soluble drugs. However, emerging evidence suggests that IT drug delivery be an efficacious strategy for the treatment of cancers in which there is leptomeningeal spread of disease. In this review, we discuss CSF flow dynamics and CSF clearance pathways in the context of intrathecal delivery. We discuss human and animal studies of several new classes of therapeutic agents-cellular, protein, nucleic acid, and nanoparticle-based small molecules-that may benefit from IT delivery. The complexity of the CSF compartment presents several key challenges in predicting biodistribution of IT-delivered drugs. New approaches and strategies are needed that can overcome the high rates of turnover in the CSF to reach specific tissues or cellular targets.

Intrathecal Drug Delivery: Advances and Applications in the Management of Chronic Pain Patient

Advances in our understanding of the biology of spinal systems in organizing and defining the content of exteroceptive information upon which higher centers define the state of the organism and its role in the regulation of somatic and automatic output, defining the motor response of the organism, along with the unique biology and spatial organization of this space, have resulted in an increased focus on therapeutics targeted at this extracranial neuraxial space. Intrathecal (IT) drug delivery systems (IDDS) are well-established as an effective therapeutic approach to patients with chronic non-malignant or malignant pain and as a tool for management of patients with severe spasticity and to deliver therapeutics that address a myriad of spinal pathologies. The risk to benefit ratio of IDD makes it a useful interventional approach. While not without risks, this approach has a significant therapeutic safety margin when employed using drugs with a validated safety profile and by skilled practioners. The present review addresses current advances in our understanding of the biology and dynamics of the intrathecal space, therapeutic platforms, novel therapeutics, delivery technology, issues of safety and rational implementation of its therapy, with a particular emphasis upon the management of pain.

Importance of CSF-based Aβ clearance with age in humans increases with declining efficacy of blood-brain barrier/proteolytic pathways

The kinetics of amyloid beta turnover within human brain is still poorly understood. We previously found a dramatic decline in the turnover of Aβ peptides in normal aging. It was not known if brain interstitial fluid/cerebrospinal fluid (ISF/CSF) fluid exchange, CSF turnover, blood-brain barrier function or proteolysis were affected by aging or the presence of β amyloid plaques. Here, we describe a non-steady state physiological model developed to decouple CSF fluid transport from other processes. Kinetic parameters were estimated using: (1) MRI-derived brain volumes, (2) stable isotope labeling kinetics (SILK) of amyloid-β peptide (Aβ), and (3) lumbar CSF Aβ concentration during SILK. Here we show that changes in blood-brain barrier transport and/or proteolysis were largely responsible for the age-related decline in Aβ turnover rates. CSF-based clearance declined modestly in normal aging but became increasingly important due to the slowing of other processes. The magnitude of CSF-based clearance was also lower than that due to blood-brain barrier function plus proteolysis. These results suggest important roles for blood-brain barrier transport and proteolytic degradation of Aβ in the development Alzheimer’s Disease in humans.

To understand if brain interstitial fluid/cerebrospinal fluid (ISF/CSF) exchange, CSF turnover, blood-brain barrier function or proteolysis were affected by aging or the presence of β amyloid plaques, Elbert et al. develop a non-steady state physiological model using MRI-derived brain volumes, stable isotope labeling kinetics of Aβ, and lumbar CSF Aβ concentration. Their model suggests an important role for blood-brain barrier transport and proteolytic degradation of Aβ in the development Alzheimer’s Disease in humans.

Boundary conditions investigation to improve computer simulation of cerebrospinal fluid dynamics in hydrocephalus patients

Three-D head geometrical models of eight healthy subjects and 11 hydrocephalus patients were built using their CINE phase-contrast MRI data and used for computer simulations under three different inlet/outlet boundary conditions (BCs). The maximum cerebrospinal fluid (CSF) pressure and the ventricular system volume were more effective and accurate than the other parameters in evaluating the patients' conditions. In constant CSF pressure, the computational patient models were 18.5% more sensitive to CSF volume changes in the ventricular system under BC "C". Pulsatile CSF flow rate diagrams were used for inlet and outlet BCs of BC "C". BC "C" was suggested to evaluate the intracranial compliance of the hydrocephalus patients. The results suggested using the computational fluid dynamic (CFD) method and the fully coupled fluid-structure interaction (FSI) method for the CSF dynamic analysis in patients with external and internal hydrocephalus, respectively.

Quantification of Oscillatory Shear Stress from Reciprocating CSF Motion on 4D Flow Imaging

BACKGROUND AND PURPOSE

Oscillatory shear stress could not be directly measured in consideration of direction, although cerebrospinal fluid has repetitive movements synchronized with heartbeat. Our aim was to evaluate the important of oscillatory shear stress in the cerebral aqueduct and foramen magnum in idiopathic normal pressure hydrocephalus by comparing it with wall shear stress and the oscillatory shear index in patients with idiopathic normal pressure hydrocephalus.

MATERIALS AND METHODS

By means of the 4D flow application, oscillatory shear stress, wall shear stress, and the oscillatory shear index were measured in 41 patients with idiopathic normal pressure hydrocephalus, 23 with co-occurrence of idiopathic normal pressure hydrocephalus and Alzheimer-type dementia, and 9 age-matched controls. These shear stress parameters at the cerebral aqueduct were compared with apertures and stroke volumes at the foramen of Magendie and cerebral aqueduct.

RESULTS

Two wall shear stress magnitude peaks during a heartbeat were changed to periodic oscillation by converting oscillatory shear stress. The mean oscillatory shear stress amplitude and time-averaged wall shear stress values at the dorsal and ventral regions of the cerebral aqueduct in the idiopathic normal pressure hydrocephalus groups were significantly higher than those in controls. Furthermore, those at the ventral region of the cerebral aqueduct in the idiopathic normal pressure hydrocephalus group were also significantly higher than those in the co-occurrence of idiopathic normal pressure hydrocephalus with Alzheimer-type dementia group. The oscillatory shear stress amplitude at the dorsal region of the cerebral aqueduct was significantly associated with foramen of Magendie diameters, whereas it was strongly associated with the stroke volume at the upper end of the cerebral aqueduct rather than that at the foramen of Magendie.

CONCLUSIONS

Oscillatory shear stress, which reflects wall shear stress vector changes better than the conventional wall shear stress magnitude and the oscillatory shear index, can be directly measured on 4D flow MR imaging. Oscillatory shear stress at the cerebral aqueduct was considerably higher in patients with idiopathic normal pressure hydrocephalus.

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.

Functional hyperemia drives fluid exchange in the paravascular space

The brain lacks a conventional lymphatic system to remove metabolic waste. It has been proposed that directional fluid movement through the arteriolar paravascular space (PVS) promotes metabolite clearance. We performed simulations to examine if arteriolar pulsations and dilations can drive directional CSF flow in the PVS and found that arteriolar wall movements do not drive directional CSF flow. We propose an alternative method of metabolite clearance from the PVS, namely fluid exchange between the PVS and the subarachnoid space (SAS). In simulations with compliant brain tissue, arteriolar pulsations did not drive appreciable fluid exchange between the PVS and the SAS. However, when the arteriole dilated, as seen during functional hyperemia, there was a marked exchange of fluid. Simulations suggest that functional hyperemia may serve to increase metabolite clearance from the PVS. We measured blood vessels and brain tissue displacement simultaneously in awake, head-fixed mice using two-photon microscopy. These measurements showed that brain deforms in response to pressure changes in PVS, consistent with our simulations. Our results show that the deformability of the brain tissue needs to be accounted for when studying fluid flow and metabolite transport.

Spatial distribution of human arachnoid trabeculae

Traumatic brain injury (TBI) is a common injury modality affecting a diverse patient population. Axonal injury occurs when the brain experiences excessive deformation as a result of head impact. Previous studies have shown that the arachnoid trabeculae (AT) in the subarachnoid space significantly influence the magnitude and distribution of brain deformation during impact. However, the quantity and spatial distribution of cranial AT in humans is unknown. Quantification of these microstructural features will improve understanding of force transfer during TBI, and may be a valuable dataset for microneurosurgical procedures. In this study, we quantify the spatial distribution of cranial AT in seven post-mortem human subjects. Optical coherence tomography (OCT) was used to conduct in situ imaging of AT microstructure across the surface of the human brain. OCT images were segmented to quantify the relative amounts of trabecular structures through a volume fraction (VF) measurement. The average VF for each brain ranged from 22.0% to 29.2%. Across all brains, there was a positive spatial correlation, with VF significantly greater by 12% near the superior aspect of the brain (p < .005), and significantly greater by 5%-10% in the frontal lobes (p < .005). These findings suggest that the distribution of AT between the brain and skull is heterogeneous, region-dependent, and likely contributes to brain deformation patterns. This study is the first to image and quantify human AT across the cerebrum and identify region-dependencies. Incorporation of this spatial heterogeneity may improve the accuracy of computational models of human TBI and enhance understanding of brain dynamics.

Methods to measure, model and manipulate fluid flow in brain

The brain consists of a complex network of cells and matrix that is cushioned and nourished by multiple types of fluids: cerebrospinal fluid, blood, and interstitial fluid. The movement of these fluids through the tissues has recently gained more attention due to implications in Alzheimer's Disease and glioblastoma. Therefore, methods to study these fluid flows are necessary and timely for the current study of neuroscience. Imaging modalities such as magnetic resonance imaging have been used clinically and pre-clinically to image flows in healthy and diseased brains. These measurements have been used to both parameterize and validate models of fluid flow both computational and in vitro. Both of these models can elucidate the changes to fluid flow that occur during disease and can assist in linking the compartments of fluid flow with one another, a difficult challenge experimentally. In vitro models, though in limited use with fluid flow, allow the examination of cellular responses to physiological flow. To determine causation, in vivo methods have been developed to manipulate flow, including both physical and pharmacological manipulations, at each point of fluid movement of origination resulting in exciting findings in the preclinical setting. With new targets, such as the brain-draining lymphatics and glymphatic system, fluid flow and tissue drainage within the brain is an exciting and growing research area. In this review, we discuss the methods that currently exist to examine and test hypotheses related to fluid flow in the brain as we attempt to determine its impact on neural function.

Impact of Neurapheresis System on Intrathecal Cerebrospinal Fluid Dynamics: A Computational Fluid Dynamics Study

It has been hypothesized that early and rapid filtration of blood from cerebrospinal fluid (CSF) in postsubarachnoid hemorrhage patients may reduce hospital stay and related adverse events. In this study, we formulated a subject-specific computational fluid dynamics (CFD) model to parametrically investigate the impact of a novel dual-lumen catheter-based CSF filtration system, the Neurapheresis™ system (Minnetronix Neuro, Inc., St. Paul, MN), on intrathecal CSF dynamics. The operating principle of this system is to remove CSF from one location along the spine (aspiration port), externally filter the CSF routing the retentate to a waste bag, and return permeate (uncontaminated CSF) to another location along the spine (return port). The CFD model allowed parametric simulation of how the Neurapheresis system impacts intrathecal CSF velocities and steady-steady streaming under various Neurapheresis flow settings ranging from 0.5 to 2.0 ml/min and with a constant retentate removal rate of 0.2 ml/min simulation of the Neurapheresis system were compared to a lumbar drain simulation with a typical CSF removal rate setting of 0.2 ml/min. Results showed that the Neurapheresis system at a maximum flow of 2.0 ml/min increased average steady streaming CSF velocity 2× in comparison to lumbar drain (0.190 ± 0.133 versus 0.093 ± 0.107 mm/s, respectively). This affect was localized to the region within the Neurapheresis flow loop. The mean velocities introduced by the flow loop were relatively small in comparison to normal cardiac-induced CSF velocities.

Anthropomorphic Model of Intrathecal Cerebrospinal Fluid Dynamics Within the Spinal Subarachnoid Space: Spinal Cord Nerve Roots Increase Steady-Streaming

Cerebrospinal fluid (CSF) dynamics are thought to play a vital role in central nervous system (CNS) physiology. The objective of this study was to investigate the impact of spinal cord (SC) nerve roots (NR) on CSF dynamics. A subject-specific computational fluid dynamics (CFD) model of the complete spinal subarachnoid space (SSS) with and without anatomically realistic NR and nonuniform moving dura wall deformation was constructed. This CFD model allowed detailed investigation of the impact of NR on CSF velocities that is not possible in vivo using magnetic resonance imaging (MRI) or other noninvasive imaging methods. Results showed that NR altered CSF dynamics in terms of velocity field, steady-streaming, and vortical structures. Vortices occurred in the cervical spine around NR during CSF flow reversal. The magnitude of steady-streaming CSF flow increased with NR, in particular within the cervical spine. This increase was located axially upstream and downstream of NR due to the interface of adjacent vortices that formed around NR.

FSI simulation of CSF hydrodynamic changes in a large population of non-communicating hydrocephalus patients during treatment process with regard to their clinical symptoms

3D fluid-structure interaction modelling was utilized for simulation of 13 normal subjects, 11 non-communicating hydrocephalus (NCH) patients at pre-treatment phase, and 3 patients at five post-treatment phases. Evaluation of ventricles volume and maximum CSF pressure (before shunting) following results validation indicated that these parameters were the most proper hydrodynamic indices and the NCH type doesn't have any significant effect on changes in two indices. The results confirmed an appropriate correlation between these indices although the correlation decreased slightly after the occurrence of disease. NCH raises the intensity of vortex and pulsatility (2.4 times) of CSF flow while the flow remains laminar. On day 18 after shunting, the CSF pressure decreased 81.0% and all clinical symptoms of patients vanished except for headache. Continuing this investigation during the treatment process showed that maximum CSF pressure is the most sensitive parameter to patients' clinical symptoms. Maximum CSF pressure has decreased proportional to the level of decrease in clinical symptoms and has returned close to the pressure range in normal subjects faster than other parameters and simultaneous with disappearance of patients' clinical symptoms (from day 81 after shunting). However, phase lag between flow rate and pressure gradient functions and the degree of CSF pulsatility haven't returned to normal subjects' conditions even 981 days after shunting and NCH has also caused a permanent volume change (of 20.1%) in ventricles. Therefore, patients have experienced a new healthy state in new hydrodynamic conditions after shunting and healing. Increase in patients' intracranial compliance was predicted with a more accurate non-invasive method than previous experimental methods up to more than 981 days after shunting. The changes in hydrodynamic parameters along with clinical reports of patients can help to gain more insight into the pathophysiology of NCH patients.

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.

Subarachnoid Trabeculae: A Comprehensive Review of Their Embryology, Histology, Morphology, and Surgical Significance

INTRODUCTION

Brain is suspended in cerebrospinal fluid (CSF)-filled subarachnoid space by subarachnoid trabeculae (SAT), which are collagen-reinforced columns stretching between the arachnoid and pia maters. Much neuroanatomic research has been focused on the subarachnoid cisterns and arachnoid matter but reported data on the SAT are limited. This study provides a comprehensive review of subarachnoid trabeculae, including their embryology, histology, morphologic variations, and surgical significance.

METHODS

A literature search was conducted with no date restrictions in PubMed, Medline, EMBASE, Wiley Online Library, Cochrane, and Research Gate. Terms for the search included but were not limited to subarachnoid trabeculae, subarachnoid trabecular membrane, arachnoid mater, subarachnoid trabeculae embryology, subarachnoid trabeculae histology, and morphology. Articles with a high likelihood of bias, any study published in nonpopular journals (not indexed in PubMed or MEDLINE), and studies with conflicting data were excluded.

RESULTS

A total of 1113 articles were retrieved. Of these, 110 articles including 19 book chapters, 58 original articles, 31 review articles, and 2 case reports met our inclusion criteria.

CONCLUSIONS

SAT provide mechanical support to neurovascular structures through cell-to-cell interconnections and specific junctions between the pia and arachnoid maters. They vary widely in appearance and configuration among different parts of the brain. The complex network of SAT is inhomogeneous and mainly located in the vicinity of blood vessels. Microsurgical procedures should be performed with great care, and sharp rather than blunt trabecular dissection is recommended because of the close relationship to neurovascular structures. The significance of SAT for cerebrospinal fluid flow and hydrocephalus is to be determined.

Fluid structure interaction model analysis of cerebrospinal fluid circulation in patients with continuous-flow left ventricular assist devices

PURPOSE

The current 1-dimensional fluid structure interaction model (FSI) for understanding cerebrospinal fluid (CSF) circulation requires pulsatility as a precondition and has not been applied to patients with continuous-flow left ventricular assist devices (CF-LVAD) where pulsatility is chronically reduced. Our study aims to characterize the behavior of CSF pressure and flow in patients with CF-LVADs using a computational FSI model.

METHODS

Utilizing the computational FSI model, CSF production in choroid plexuses of the 4 ventricles was specified as a boundary condition for the model. The other source of production from capillary ultrafiltrate spaces was accounted for by the mass conservation equation. The primary CSF absorption sites (i.e., arachnoid granulations) were treated as the outlet boundary conditions. We established a low pulse wave to represent patients with a CF-LVAD.

RESULTS

From the model, low pulse conditions resulted in a reduction in CSF pressure amplitude and velocity though the overall flow rate was unchanged.

CONCLUSIONS

The existing FSI model is not a suitable representation of CSF flow in CF-LVAD patients. More studies are needed to elucidate the role of pulsatility in CSF flow and the compensatory changes in CSF production and absorption that occur in patients with CF-LVADs in whom pulsatility is diminished.

Direct numerical simulation of transitional hydrodynamics of the cerebrospinal fluid in Chiari I malformation: The role of cranio-vertebral junction

Obstruction to the cerebrospinal fluid (CSF) outflow caused by the herniation of cerebellar tonsils as a result of Chiari malformation type I leads to altered CSF hydrodynamics. This contribution explores the minutest characteristics of the CSF hydrodynamics in cervical subarachnoid space (SAS) of a healthy subject and 2 Chiari patients by performing highly resolved direct numerical simulation. The lattice Boltzmann method is used for the simulations because of its scalability on modern supercomputers that allow us to simulate up to approximately 10⁹ cells while resolving the Kolmogorov microscales. The results depict that whereas the complex CSF flow remains largely laminar in the SAS of a healthy subject, constriction of the cranio-vertebral junction in Chiari I patients causes manifold fluctuations in the hydrodynamics of the CSF. These fluctuations resemble a flow that is in a transitional regime rather than laminar or fully developed turbulence. The fluctuations confine near the cranio-vertebral junction and are triggered due to the tonsillar herniation, which perturbs the flow as a result of altered anatomy of the SAS.

The Connected Steady State Model and the Interdependence of the CSF Proteome and CSF Flow Characteristics

Here we show that the hydrodynamic radii-dependent entry of blood proteins into cerebrospinal fluid (CSF) can best be modeled with a diffusional system of consecutive interdependent steady states between barrier-restricted molecular flux and bulk flow of CSF. The connected steady state model fits precisely to experimental results and provides the theoretical backbone to calculate the in-vivo hydrodynamic radii of blood-derived proteins as well as individual barrier characteristics. As the experimental reference set we used a previously published large-scale patient cohort of CSF to serum quotient ratios of immunoglobulins in relation to the respective albumin quotients. We related the inter-individual variances of these quotient relationships to the individual CSF flow time and barrier characteristics. We claim that this new concept allows the diagnosis of inflammatory processes with Reibergrams derived from population-based thresholds to be shifted to individualized judgment, thereby improving diagnostic sensitivity. We further use the source-dependent gradient patterns of proteins in CSF as intrinsic tracers for CSF flow characteristics. We assume that the rostrocaudal gradient of blood-derived proteins is a consequence of CSF bulk flow, whereas the slope of the gradient is a consequence of the unidirectional bulk flow and bidirectional pulsatile flow of CSF. Unlike blood-derived proteins, the influence of CSF flow characteristics on brain-derived proteins in CSF has been insufficiently discussed to date. By critically reviewing existing experimental data and by reassessing their conformity to CSF flow assumptions we conclude that the biomarker potential of brain-derived proteins in CSF can be improved by considering individual subproteomic dynamics of the CSF system.

Barrier dysfunction or drainage reduction: differentiating causes of CSF protein increase

BACKGROUND

Cerebrospinal fluid (CSF) protein analysis is an important element in the diagnostic chain for various central nervous system (CNS) pathologies. Among multiple existing approaches to interpreting measured protein levels, the Reiber diagram is particularly robust with respect to physiologic inter-individual variability, as it uses multiple subject-specific anchoring values. Beyond reliable identification of abnormal protein levels, the Reiber diagram has the potential to elucidate their pathophysiologic origin. In particular, both reduction of CSF drainage from the cranio-spinal space as well as blood-CNS barrier dysfunction have been suggested ρas possible causes of increased concentration of blood-derived proteins. However, there is disagreement on which of the two is the true cause.

METHODS

We designed two computational models to investigate the mechanisms governing protein distribution in the spinal CSF. With a one-dimensional model, we evaluated the distribution of albumin and immunoglobulin G (IgG), accounting for protein transport rates across blood-CNS barriers, CSF dynamics (including both dispersion induced by CSF pulsations and advection by mean CSF flow) and CSF drainage. Dispersion coefficients were determined a priori by computing the axisymmetric three-dimensional CSF dynamics and solute transport in a representative segment of the spinal canal.

RESULTS

Our models reproduce the empirically determined hyperbolic relation between albumin and IgG quotients. They indicate that variation in CSF drainage would yield a linear rather than the expected hyperbolic profile. In contrast, modelled barrier dysfunction reproduces the experimentally observed relation.

CONCLUSIONS

High levels of albumin identified in the Reiber diagram are more likely to originate from a barrier dysfunction than from a reduction in CSF drainage. Our in silico experiments further support the hypothesis of decreasing spinal CSF drainage in rostro-caudal direction and emphasize the physiological importance of pulsation-driven dispersion for the transport of large molecules in the CSF.

A numerical investigation of intrathecal isobaric drug dispersion within the cervical subarachnoid space

Intrathecal drug and gene vector delivery is a procedure to release a solute within the cerebrospinal fluid. This procedure is currently used in clinical practice and shows promise for treatment of several central nervous system pathologies. However, intrathecal delivery protocols and systems are not yet optimized. The aim of this study was to investigate the effects of injection parameters on solute distribution within the cervical subarachnoid space using a numerical platform. We developed a numerical model based on a patient-specific three dimensional geometry of the cervical subarachnoid space with idealized dorsal and ventral nerve roots and denticulate ligament anatomy. We considered the drug as massless particles within the flow field and with similar properties as the CSF, and we analyzed the effects of anatomy, catheter position, angle and injection flow rate on solute distribution within the cerebrospinal fluid by performing a series of numerical simulations. Results were compared quantitatively in terms of drug peak concentration, spread, accumulation rate and appearance instant over 15 seconds following the injection. Results indicated that solute distribution within the cervical spine was altered by all parameters investigated within the time range analyzed following the injection. The presence of spinal cord nerve roots and denticulate ligaments increased drug spread by 60% compared to simulations without these anatomical features. Catheter position and angle were both found to alter spread rate up to 86%, and catheter flow rate altered drug peak concentration up to 78%. The presented numerical platform fills a first gap towards the realization of a tool to parametrically assess and optimize intrathecal drug and gene vector delivery protocols and systems. Further investigation is needed to analyze drug spread over a longer clinically relevant time frame.

Choroidal fissure acts as an overflow device in cerebrospinal fluid drainage: morphological comparison between idiopathic and secondary normal-pressure hydrocephalus

To clarify the pathogenesis of two different types of adult-onset normal-pressure hydrocephalus (NPH), we investigated cerebrospinal fluid distribution on the high-field three-dimensional MRI. The subarachnoid spaces in secondary NPH were smaller than those in the controls, whereas those in idiopathic NPH were of similar size to the controls. In idiopathic NPH, however, the basal cistern and Sylvian fissure were enlarged in concurrence with ventricular enlargement towards the z-direction, but the convexity subarachnoid space was severely diminished. In this article, we provide evidence that the key cause of the disproportionate cerebrospinal fluid distribution in idiopathic NPH is the compensatory direct CSF communication between the inferior horn of the lateral ventricles and the ambient cistern at the choroidal fissure. In contrast, all parts of the subarachnoid spaces were equally and severely decreased in secondary NPH. Blockage of CSF drainage from the subarachnoid spaces could cause the omnidirectional ventricular enlargement in secondary NPH.

Quantitative measurement of cerebral blood volume using velocity-selective pulse trains

PURPOSE

To develop a non-contrast-enhanced MRI method for cerebral blood volume (CBV) mapping using velocity-selective (VS) pulse trains.

METHODS

The new pulse sequence applied velocity-sensitive gradient waveforms in the VS label modules and velocity-compensated ones in the control scans. Sensitivities to the gradient imperfections (e.g., eddy currents) were evaluated through phantom studies. CBV quantification procedures based on simulated labeling efficiencies for arteriolar, capillary, and venular blood as a function of cutoff velocity (Vc) are presented. Experiments were conducted on healthy volunteers at 3T to examine the effects of unbalanced diffusion weighting, cerebrospinal (CSF) contamination and variation of Vc.

RESULTS

Phantom results of the used VS pulse trains demonstrated robustness to eddy currents. The mean CBV values of gray matter and white matter for the experiments using Vc = 3.5 mm/s and velocity-compensated control with CSF-nulling were 5.1 ± 0.6 mL/100 g and 2.4 ± 0.2 mL/100 g, respectively, which were 23% and 32% lower than results from the experiment with velocity-insensitive control, corresponding to 29% and 25% lower in averaged temporal signal-to-noise ratio values.

CONCLUSION

A novel technique using VS pulse trains was demonstrated for CBV mapping. The results were both qualitatively and quantitatively close to those from existing methods. Magn Reson Med 77:92-101, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Computational Investigation of Cerebrospinal Fluid Dynamics in the Posterior Cranial Fossa and Cervical Subarachnoid Space in Patients with Chiari I Malformation

Purpose

Previous computational fluid dynamics (CFD) studies have demonstrated that the Chiari malformation is associated with abnormal cerebrospinal fluid (CSF) flow in the cervical part of the subarachnoid space (SAS), but the flow in the SAS of the posterior cranial fossa has received little attention. This study extends previous modelling efforts by including the cerebellomedullary cistern, pontine cistern, and 4th ventricle in addition to the cervical subarachnoid space.

Methods

The study included one healthy control, Con1, and two patients with Chiari I malformation, P1 and P2. Meshes were constructed by segmenting images obtained from T2-weighted turbo spin-echo sequences. CFD simulations were performed with a previously verified and validated code. Patient-specific flow conditions in the aqueduct and the cervical SAS were used. Two patients with the Chiari malformation and one control were modelled.

Results

The results demonstrated increased maximal flow velocities in the Chiari patients, ranging from factor 5 in P1 to 14.8 in P2, when compared to Con1 at the level of Foramen Magnum (FM). Maximal velocities in the cervical SAS varied by a factor 2.3, while the maximal flow in the aqueduct varied by a factor 3.5. The pressure drop from the pontine cistern to the cervical SAS was similar in Con1 and P1, but a factor two higher in P2. The pressure drop between the aqueduct and the cervical SAS varied by a factor 9.4 where P1 was the one with the lowest pressure jump and P2 and Con1 differed only by a factor 1.6.

Conclusion

This pilot study demonstrates that including the posterior cranial fossa is feasible and suggests that previously found flow differences between Chiari I patients and healthy individuals in the cervical SAS may be present also in the SAS of the posterior cranial fossa.

Multiscale modeling and simulation of brain blood flow

The aim of this work is to present an overview of recent advances in multi-scale modeling of brain blood flow. In particular, we present some approaches that enable the in silico study of multi-scale and multi-physics phenomena in the cerebral vasculature. We discuss the formulation of continuum and atomistic modeling approaches, present a consistent framework for their concurrent coupling, and list some of the challenges that one needs to overcome in achieving a seamless and scalable integration of heterogeneous numerical solvers. The effectiveness of the proposed framework is demonstrated in a realistic case involving modeling the thrombus formation process taking place on the wall of a patient-specific cerebral aneurysm. This highlights the ability of multi-scale algorithms to resolve important biophysical processes that span several spatial and temporal scales, potentially yielding new insight into the key aspects of brain blood flow in health and disease. Finally, we discuss open questions in multi-scale modeling and emerging topics of future research.

MRI-guided epilepsy detection

One of the most common neurological brain disorder is epilepsy that happen as an abrupt seizure. Around 30% of patients with epilepsy resist to all forms of medical treatments and, therefore, the removal of epileptic brain tissue is the only solution to get these patients free from chronical seizures. Discovering the epileptic region is a first key into the treatment. In this paper, we introduced a method for epilepsy detection. In this method superparamagnetic nanoparticle, (SPMN) is used as a sensing material in order to investigate the epileptic area. Based on the magnetic field, first they are crossed through the Blood Brain Barrier (BBB). They can cross the blood-brain barrier into the brain by means of magnetic forces. In this study, the optimal force for crossing to the brain and nanoparticles aggregation by means of MRI magnetic field for crossing and weak magnetic field inside the brain have been considered. Nanoparticles aggregation can be used as a marker to increase the contrast of MRI images in the epileptic brain area.

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.

Cerebrospinal fluid outflow along lumbar nerves and possible relevance for pain research: case report and review

CSF outflow through the cribriform plate near the olfactory nerves and the outflow along brain and spinal nerves are together known as peripheral CSF outflow pathway (PCOP). It is still not clear whether the PCOP has pathogenetic relevance. Our previous clinical observations have indicated that CSF may interact with nerves along the PCOP and in this article we present our finding of CSF outflow demonstrated by myelography in a single patient. We also discuss unexplained experimental pain pathomechanisms against the background of the PCOP hypothesis. We observed that CSF flowed along lumbar nerves in distal direction at a speed of about 10 cm per hour on its way through the tissues, mainly muscles. Total CSF outflow volume at the lumbar site was remarkable. CSF outflow at lumbar nerves was also documented by neuroradiology. It is plausible that CSF signaling serves for interaction with nerves along the PCOP, which could explain previously unknown pathomechanisms in pain generation. Experimental findings of tactile pain hypersensitivity within lumbosacral pain pathways could be explained by releasing of molecules, microparticles, or exosomes into the CSF by mast cells, which then move with CSF outflow along the PCOP and interact with nerves, initiating even retrograde synaptic stripping.

CNS wide simulation of flow resistance and drug transport due to spinal microanatomy

Spinal microstructures are known to substantially affect cerebrospinal fluid patterns, yet their actual impact on flow resistance has not been quantified. Because the length scale of microanatomical aspects is below medical image resolution, their effect on flow is difficult to observe experimentally. Using a computational fluid mechanics approach, we were able to quantify the contribution of micro-anatomical aspects on cerebrospinal fluid (CSF) flow patterns and flow resistance within the entire central nervous system (CNS). Cranial and spinal CSF filled compartments were reconstructed from human imaging data; microscopic trabeculae below the image detection threshold were added artificially. Nerve roots and trabeculae were found to induce regions of microcirculation, whose location, size and vorticity along the spine were characterized. Our CFD simulations based on volumetric flow rates acquired with Cine Phase Contrast MRI in a normal human subject suggest a 2-2.5 fold increase in pressure drop mainly due to arachnoid trabeculae. The timing and phase lag of the CSF pressure and velocity waves along the spinal canal were also computed, and a complete spatio-temporal map encoding CSF volumetric flow rates and pressure was created. Micro-anatomy induced fluid patterns were found responsible for the rapid caudo-cranial spread of an intrathecally administered drug. The speed of rostral drug dispersion is drastically accelerated through pulsatile flow around microanatomy induced vortices. Exploring massive parallelization on a supercomputer, the feasibility of computational drug transport studies was demonstrated. CNS-wide simulations of intrathecal drugs administration can become a practical tool for in silico design, interspecies scaling and optimization of experimental drug trials.

Spinal fluid biomechanics and imaging: an update for neuroradiologists

Flow imaging with cardiac-gated phase-contrast MR has applications in the management of neurologic disorders. Together with computational fluid dynamics, phase-contrast MR has advanced our understanding of spinal CSF flow. Phase-contrast MR is used to evaluate patients with Chiari I malformation who are candidates for surgical treatment. In theory, abnormal CSF flow resulting from the abnormal tonsil position causes syringomyelia and other neurologic signs and symptoms in patients with Chiari I. CSF flow imaging also has research applications in syringomyelia and spinal stenosis. To optimize MR acquisition and interpretation, neuroradiologists must have familiarity with healthy and pathologic patterns of CSF flow. The purpose of this review is to update concepts of CSF flow that are important for the practice of flow imaging in the spine.

Flow induced by ependymal cilia dominates near-wall cerebrospinal fluid dynamics in the lateral ventricles

While there is growing experimental evidence that cerebrospinal fluid (CSF) flow induced by the beating of ependymal cilia is an important factor for neuronal guidance, the respective contribution of vascular pulsation-driven macroscale oscillatory CSF flow remains unclear. This work uses computational fluid dynamics to elucidate the interplay between macroscale and cilia-induced CSF flows and their relative impact on near-wall dynamics. Physiological macroscale CSF dynamics are simulated in the ventricular space using subject-specific anatomy, wall motion and choroid plexus pulsations derived from magnetic resonance imaging. Near-wall flow is quantified in two subdomains selected from the right lateral ventricle, for which dynamic boundary conditions are extracted from the macroscale simulations. When cilia are neglected, CSF pulsation leads to periodic flow reversals along the ventricular surface, resulting in close to zero time-averaged force on the ventricle wall. The cilia promote more aligned wall shear stresses that are on average two orders of magnitude larger compared with those produced by macroscopic pulsatile flow. These findings indicate that CSF flow-mediated neuronal guidance is likely to be dominated by the action of the ependymal cilia in the lateral ventricles, whereas CSF dynamics in the centre regions of the ventricles is driven predominantly by wall motion and choroid plexus pulsation.

Idiopathic intracranial hypertension; research progress and emerging themes

Idiopathic intracranial hypertension (IIH) is a condition characterised by increased intracranial pressure of unknown cause predominantly seen in obese women of childbearing age and associated with a history of recent weight gain. The aetiology is poorly understood and there are no evidence-based guidelines on the management of the disease. We aim to provide a review of the recent literature outlining the latest advances in this field over the past few years. Areas of emerging interest related to the pathophysiology of IIH will be discussed, such as the role of obesity, adipose tissue and 11β-hydroxysteroid dehydrogenase type 1. We consider the latest research on the role of venous sinus stenosis in IIH and ex vivo advances into cerebrospinal fluid drainage via the arachnoid granulation tissue. The latest techniques for optic nerve head evaluation and the role of optical coherence tomography will be summarised. Finally, we will discuss recent advances in the management of IIH, including weight loss, and medical and surgical treatment strategies.

Updated physiology and pathophysiology of CSF circulation--the pulsatile vector theory

INTRODUCTION

Hydrocephalus is still a not well-understood diagnostic and a therapeutic dilemma because of the lack of sufficient and comprehensive model of cerebrospinal fluid circulation and pathological alterations.

CONCLUSIONS

Based on current studies, reviews, and knowledge of cerebrospinal fluid dynamics, brain water dynamics, intracranial pressure, and cerebral perfusion physiology, a new concept is deducted that can describe normal and pathological changes of cerebrospinal fluid circulation and pathophysiology of idiopathic intracranial hypertension.

Updating the mild encephalitis hypothesis of schizophrenia

Schizophrenia seems to be a heterogeneous disorder. Emerging evidence indicates that low level neuroinflammation (LLNI) may not occur infrequently. Many infectious agents with low overall pathogenicity are risk factors for psychoses including schizophrenia and for autoimmune disorders. According to the mild encephalitis (ME) hypothesis, LLNI represents the core pathogenetic mechanism in a schizophrenia subgroup that has syndromal overlap with other psychiatric disorders. ME may be triggered by infections, autoimmunity, toxicity, or trauma. A 'late hit' and gene-environment interaction are required to explain major findings about schizophrenia, and both aspects would be consistent with the ME hypothesis. Schizophrenia risk genes stay rather constant within populations despite a resulting low number of progeny; this may result from advantages associated with risk genes, e.g., an improved immune response, which may act protectively within changing environments, although they are associated with the disadvantage of increased susceptibility to psychotic disorders. Specific schizophrenic symptoms may arise with instances of LLNI when certain brain functional systems are involved, in addition to being shaped by pre-existing liability factors. Prodrome phase and the transition to a diseased status may be related to LLNI processes emerging and varying over time. The variability in the course of schizophrenia resembles the varying courses of autoimmune disorders, which result from three required factors: genes, the environment, and the immune system. Preliminary criteria for subgrouping neurodevelopmental, genetic, ME, and other types of schizophrenias are provided. A rare example of ME schizophrenia may be observed in Borna disease virus infection. Neurodevelopmental schizophrenia due to early infections has been estimated by others to explain approximately 30% of cases, but the underlying pathomechanisms of transition to disease remain in question. LLNI (e.g. from reactivation related to persistent infection) may be involved and other pathomechanisms including dysfunction of the blood-brain barrier or the blood-CSF barrier, CNS-endogenous immunity and the volume transmission mode balancing wiring transmission (the latter represented mainly by synaptic transmission, which is often described as being disturbed in schizophrenia). Volume transmission is linked to CSF signaling; and together could represent a common pathogenetic link for the distributed brain dysfunction, dysconnectivity, and brain structural abnormalities observed in schizophrenia. In addition, CSF signaling may extend into peripheral tissues via the CSF outflow pathway along brain nerves and peripheral nerves, and it may explain the peripheral topology of neuronal dysfunctions found, like in olfactory dysfunction, dysautonomia, and even in peripheral tissues, i.e., the muscle lesions that were found in 50% of cases. Modulating factors in schizophrenia, such as stress, hormones, and diet, are also modulating factors in the immune response. Considering recent investigations of CSF, the ME schizophrenia subgroup may constitute approximately 40% of cases.

Space physiology IV: mathematical modeling of the cardiovascular system in space exploration

Mathematical modeling represents an important tool for analyzing cardiovascular function during spaceflight. This review describes how modeling of the cardiovascular system can contribute to space life science research and illustrates this process via modeling efforts to study postflight orthostatic intolerance (POI), a key issue for spaceflight. Examining this application also provides a context for considering broader applications of modeling techniques to the challenges of bioastronautics. POI, which affects a large fraction of astronauts in stand tests upon return to Earth, presents as dizziness, fainting and other symptoms, which can diminish crew performance and cause safety hazards. POI on the Moon or Mars could be more critical. In the field of bioastronautics, POI has been the dominant application of cardiovascular modeling for more than a decade, and a number of mechanisms for POI have been investigated. Modeling approaches include computational models with a range of incorporated factors and hemodynamic sophistication, and also physical models tested in parabolic and orbital flight. Mathematical methods such as parameter sensitivity analysis can help identify key system mechanisms. In the case of POI, this could lead to more effective countermeasures. Validation is a persistent issue in modeling efforts, and key considerations and needs for experimental data to synergistically improve understanding of cardiovascular responses are outlined. Future directions in cardiovascular modeling include subject-specific assessment of system status, as well as research on integrated physiological responses, leading, for instance, to assessment of subject-specific susceptibility to POI or effects of cardiovascular alterations on muscular, vision and cognitive function.

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.

Venous oxygenation mapping using velocity-selective excitation and arterial nulling

A new MRI technique to map the oxygenation of venous blood is presented. The method uses velocity-selective excitation and arterial nulling pulses, combined with phase sensitive signal detection to isolate the venous blood signal. T₂ of this signal along with a T₂-Y calibration curve yields estimates of venous oxygenation in situ. Results from phantoms and healthy human subjects under normoxic and hypoxic conditions are shown, and venous saturation levels estimated from both sagittal sinus and gray matter-based regions of interest are compared with the related techniques TRUST and QUIXOTIC. In addition, combined with an additional scan without arterial nulling pulses, the oxygen saturation level on arterial side can also be estimated.

Microscopic particles in two fractions of fresh cerebrospinal fluid in twins with schizophrenia or bipolar disorder and in healthy controls

Background

Using scanning electron microscopy, microscopic structures have been identified in fresh cerebrospinal fluid (CSF) in patients with schizophrenia and bipolar disorder, but only rarely in control subjects. However, it has not been determined whether these microscopic particles represent state or trait markers, i.e. if their presence is related to clinical manifestations of the disease or if they also can be found in as yet asymptomatic individuals with a genetic liability. This question can be addressed by studying twins discordant or concordant for schizophrenia or bipolar disorder.

Methodology/Principal Findings

We investigated microscopic structures in CSF in 102 individuals: 21 monozygotic and 16 dizygotic twins affected or not affected with schizophrenia, schizoaffective disorder or bipolar disorder and in 65 healthy singleton controls. A first and a second fraction of CSF was freshly applied on filters and examined by scanning electron microscopy technique. Spherical particles with lipid appearance averaging between 0.1 to 8.0 µm in diameter were detected in the center of the filter as well as located in the margins of larger aggregates binding in a viscous state. Structures were found in 12 of 17 probands, 5 of 12 healthy co-twins and 3 of 73 healthy controls. Thus, a positive microscopic finding significantly increased the likelihood of belonging to the proband group (OR = 48, 95% CL: 8.2–550, p<0.0001) and the co-twin-group (OR = 16, 95% CL: 2.0–218, p = 0.006). Age, sex, history of alcohol abuse or anxiety syndrome, somatic disorder and markers of acute inflammatory activity did not account for group differences; nor did exposure to psychotropic medication.

Conclusion

Presence of microscopic particles in CSF may possibly reflect trait dependent genetic or environmental vulnerability in patients with schizophrenia, schizoaffective disorder or bipolar disorder.

In vitro effects of fetal rat cerebrospinal fluid on viability and neuronal differentiation of PC12 cells

Background

Fetal cerebrospinal fluid (CSF) contains many neurotrophic and growth factors and has been shown to be capable of supporting viability, proliferation and differentiation of primary cortical progenitor cells. Rat pheochromocytoma PC12 cells have been widely used as an in vitro model of neuronal differentiation since they differentiate into sympathetic neuron-like cells in response to growth factors. This study aimed to establish whether PC12 cells were responsive to fetal CSF and therefore whether they might be used to investigate CSF physiology in a stable cell line lacking the time-specific response patterns of primary cells previously described.

Methods

In vitro assays of viability, proliferation and differentiation were carried out after incubation of PC12 cells in media with and without addition of fetal rat CSF. An MTT tetrazolium assay was used to assess cell viability and/or cell proliferation. Expression of neural differentiation markers (MAP-2 and β-III tubulin) was determined by immunocytochemistry. Formation and growth of neurites was measured by image analysis.

Results

PC12 cells differentiate into neuronal cell types when exposed to bFGF. Viability and cell proliferation of PC12 cells cultured in CSF-supplemented medium from E18 rat fetuses were significantly elevated relative to the control group. Neuronal-like outgrowths from cells appeared following the application of bFGF or CSF from E17 and E19 fetuses but not E18 or E20 CSF. Beta-III tubulin was expressed in PC12 cells cultured in any media except that supplemented with E18 CSF. MAP-2 expression was found in control cultures and in those with E17 and E19 CSF. MAP2 was located in neurites except in E17 CSF when the whole cell was positive.

Conclusions

Fetal rat CSF supports viability and stimulates proliferation and neurogenic differentiation of PC12 cells in an age-dependent way, suggesting that CSF composition changes with age. This feature may be important in vivo for the promotion of normal brain development. There were significant differences in the effects on PC12 cells compared to primary cortical cells. This suggests there is an interaction in vivo between developmental stage of cells and the composition of CSF. The data presented here support an important, perhaps driving role for CSF composition, specifically neurotrophic factors, in neuronal survival, proliferation and differentiation. The effects of CSF on PC12 cells can thus be used to further investigate the role of CSF in driving development without the confounding issues of using primary cells.