Quantifying the influence of respiration and cardiac pulsations on cerebrospinal fluid dynamics using real-time phase-contrast MRI

Real-time quantification of in vivo cerebrospinal fluid velocity using the functional magnetic resonance imaging inflow effect

In vivo estimation of cerebrospinal fluid (CSF) velocity is crucial for understanding the glymphatic system and its potential role in neurodegenerative disorders such as Alzheimer's disease and Parkinson's disease. Current cardiac or respiratory-gated approaches, such as 4D flow magnetic resonance imaging (MRI), cannot capture CSF movement in real time because of limited temporal resolution and, in addition, deteriorate in accuracy at low fluid velocities. Other techniques like real-time phase-contrast-MRI or time-spatial labeling inversion pulse are not limited by temporal averaging but have limited availability, even in research settings. This study aims to quantify the inflow effect of dynamic CSF motion on functional MRI (fMRI) for in vivo, real-time measurement of CSF flow velocity. We considered linear and nonlinear models of velocity waveforms and empirically fit them to fMRI data from a controlled flow experiment. To assess the utility of this methodology in human data, CSF flow velocities were computed from fMRI data acquired in eight healthy volunteers. Breath-holding regimens were used to amplify CSF flow oscillations. Our experimental flow study revealed that CSF velocity is nonlinearly related to inflow effect-mediated signal increase and well estimated using an extension of a previous nonlinear framework. Using this relationship, we recovered velocity from in vivo fMRI signal, demonstrating the potential of our approach for estimating CSF flow velocity in the human brain. This novel method could serve as an alternative approach to quantifying slow flow velocities in real time, such as CSF flow in the ventricular system, thereby providing valuable insights into the glymphatic system's function and its implications for neurological disorders.

Numerical analysis of the effect of Syringomyelia on cerebrospinal fluid dynamics

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

Research Trends and Visualization of Cerebrospinal Fluid Dynamics (2013-2023)

OBJECTIVE

This study aims to analyze cerebrospinal fluid (CSF) dynamics using VOSviewer, CiteSpace, and the Bibliometrix R-package software to identify research hotspots and future directions.

METHODS

Search by Web of Science Core Collection Database for related literature on CSF dynamics from 2013 to 2023. Bibliometric and visual analysis of data on number of citations, number of publications, most productive countries and institutions, important authors and journals, time of publication, popular topics, and keywords were performed by CiteSpace and VOSviewer.

RESULTS

In the field of CSF dynamics, there is a clear upward trend in annual publications. The United States, Japan, and China are among the top three countries in publishing output. The University of Copenhagen, the University of Idaho, and the University of Zurich are leading institutions in research publications. The most prolific writers in this field are Bryn A. Martin, and Olivier Baledent. Active authors and institutions in the field form multiple structurally stable research teams with each other, but the collaboration between different authors and institutional teams needs to be further strengthened. The literature with the highest citation rates in the past decade is "Blood-Brain Barrier Breakdown in the Aging Human Hippocampus," "Blood-Brain Barrier Breakdown Is an Early Biomarker of Human Cognitive Dysfunction," "Serum Neurofilament Dynamics Predicts Neurodegeneration and Clinical Progression in Presymptomatic Alzheimer's Disease," and Coupled Electrophysiological, Hemodynamic, and Cerebrospinal Fluid Oscillations in Human Sleep." Key research keywords such as CSF, hydrocephalus, dynamics, brain, blood flow, CSF, pressure, CSF flow, and MRI highlight focal areas for CSF dynamics studies. These keywords represent current research priorities and research frontiers in this field.

CONCLUSIONS

This bibliometric analysis reveals hot and future research issues in the field of CSF fluid dynamics, demonstrating the need for enhanced international collaboration and interdisciplinary research to deepen the field. Keyword analysis further clarified the research focus and provided useful guidance for subsequent studies.

Clinical magnetic resonance imaging evaluation of glymphatic function

The glymphatic system is a system of specialized perivascular spaces in the brain that facilitates removal of toxic waste solutes from the brain. Evaluation of glymphatic system function by means of magnetic resonance imaging (MRI) has thus far been largely focused on rodents because of the limitations of intrathecal delivery of gadolinium-based contrast agents to humans. This review discusses MRI methods that can be employed clinically for glymphatic-related measurements intended for early diagnosis, prevention, and the treatment of various neurological conditions. Although glymphatic system-based MRI research is in its early stages, recent studies have identified promising noninvasive MRI markers associated with glymphatic system alterations in neurological diseases. However, further optimization in data acquisition, validation, and modeling are needed to investigate the glymphatic system within the clinical setting.

Evaluating the effect of injection protocols on intrathecal solute dispersion in non-human primates: an in vitro study using a cynomolgus cerebrospinal fluid system

BACKGROUND

Achieving effective drug delivery to the central nervous system (CNS) remains a challenge for treating neurological disorders. Intrathecal (IT) delivery, which involves direct injection into the cerebrospinal fluid (CSF), presents a promising strategy. Large animal studies are important to assess the safety and efficacy of most drugs and treatments and translate the data to humans. An understanding of the influence of IT injection parameters on solute distribution within the CNS is essential to optimize preclinical research, which would potentially help design human clinical studies.

METHODS

A three-dimensional (3D) in vitro model of a cynomolgus monkey, based on MRI data, was developed to evaluate the impact of lumbar injection parameters on intrathecal solute dispersion. The parameters evaluated were (a) injection location, (b) bolus volume, (c) flush volume, (d) bolus rate, and (e) flush rate. To simulate the CSF flow within the subarachnoid space (SAS), an idealized CSF flow waveform with both cardiac and respiratory-induced components was input into the model. A solution of fluorescein drug surrogate tracer was administered in the lumbar region of the 3D in vitro model filled with deionized water. After injection of the tracer, the CSF system wide-solute dispersion was imaged using high-resolution cameras every thirty seconds for a duration of three hours. To ensure repeatability each injection protocol was repeated three times. For each protocol, the average spatial-temporal distribution over three hours post-injection, the area under the curve (AUC), and the percent injected dose (%ID) to extra-axial CSF (eaCSF) at three hours were determined.

RESULTS

The changes to the lumbar injection parameters led to variations in solute distribution along the neuro-axis. Specifically, injection location showed the most impact, enhancing the delivery to the eaCSF up to + 10.5%ID (p = 0.0282) at three hours post-injection. Adding a post-injection flush of 1.5 ml at 1 ml/min increased the solute delivery to the eaCSF by + 6.5%ID (p = 0.0218), while the larger bolus volume resulted in a + 2.3%ID (p = 0.1910) increase. The bolus and flush rates analyzed had minimal, statistically non-significant effects.

CONCLUSION

These results predict the effects of lumbar injection parameters on solute distribution in the intrathecal space in NHPs. Specifically, the choice of injection location, flush, and bolus volume significantly improved solute delivery to eaCSF. The in vitro NHP CSF model and results offer a system to help predict and optimize IT delivery protocols for pre-clinical NHP studies.

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

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

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.

Glymphatic Imaging in Pediatrics

The glymphatic system, which facilitates cerebrospinal fluid (CSF) flow through the brain parenchyma, is important for brain development and waste clearance. Advances in imaging techniques, particularly magnetic resonance imaging, have make it possible to evaluate glymphatic structures and functions in vivo. Recently, several studies have focused on the development and alterations of the glymphatic system in pediatric disorders. This review discusses the development of the glymphatic system, advances of imaging techniques and their applications in pediatric disorders. First, the results of the reviewed studies indicate that the development of the glymphatic system is a long-lasting process that continues into adulthood. Second, there is a need for improved glymphatic imaging techniques that are non-invasive and fast to improve suitability for pediatric applications, as some of existing methods use contrast injection and are susceptible to motion artifacts from long scanning times. Several novel techniques are potentially feasible for pediatric patients and may be used in the future. Third, the glymphatic dysfunction is associated with a large number of pediatric disorders, although only a few have recently been investigated. In conclusion, research on the pediatric glymphatic system remains an emerging field. The preliminary applications of glymphatic imaging techniques have provided unique insight into the pathological mechanism of pediatric diseases, but mainly limited in visualization of enlarged perivascular spaces and morphological measurements on CSF volumes. More in-depth studies on glymphatic functions are required to improve our understanding of the mechanisms underlying brain development and pediatric diseases. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 3.

Validating the accuracy of real-time phase-contrast MRI and quantifying the effects of free breathing on cerebrospinal fluid dynamics

BACKGROUND

Understanding of the cerebrospinal fluid (CSF) circulation is essential for physiological studies and clinical diagnosis. Real-time phase contrast sequences (RT-PC) can quantify beat-to-beat CSF flow signals. However, the detailed effects of free-breathing on CSF parameters are not fully understood. This study aims to validate RT-PC's accuracy by comparing it with the conventional phase-contrast sequence (CINE-PC) and quantify the effect of free-breathing on CSF parameters at the intracranial and extracranial levels using a time-domain multiparametric analysis method.

METHODS

Thirty-six healthy participants underwent MRI in a 3T scanner for CSF oscillations quantification at the cervical spine (C2-C3) and Sylvian aqueduct, using CINE-PC and RT-PC. CINE-PC uses 32 velocity maps to represent dynamic CSF flow over an average cardiac cycle, while RT-PC continuously quantifies CSF flow over 45-seconds. Free-breathing signals were recorded from 25 participants. RT-PC signal was segmented into independent cardiac cycle flow curves (Qt) and reconstructed into an averaged Qt. To assess RT-PC's accuracy, parameters such as segmented area, flow amplitude, and stroke volume (SV) of the reconstructed Qt from RT-PC were compared with those derived from the averaged Qt generated by CINE-PC. The breathing signal was used to categorize the Qt into expiratory or inspiratory phases, enabling the reconstruction of two Qt for inspiration and expiration. The breathing effects on various CSF parameters can be quantified by comparing these two reconstructed Qt.

RESULTS

RT-PC overestimated CSF area (82.7% at aqueduct, 11.5% at C2-C3) compared to CINE-PC. Stroke volumes for CINE-PC were 615 mm³ (aqueduct) and 43 mm³ (spinal), and 581 mm³ (aqueduct) and 46 mm³ (spinal) for RT-PC. During thoracic pressure increase, spinal CSF net flow, flow amplitude, SV, and cardiac period increased by 6.3%, 6.8%, 14%, and 6%, respectively. Breathing effects on net flow showed a significant phase difference compared to the other parameters. Aqueduct-CSF flows were more affected by breathing than spinal-CSF.

CONCLUSIONS

RT-PC accurately quantifies CSF oscillations in real-time and eliminates the need for cardiac synchronization, enabling the quantification of the cardiac and breathing components of CSF flow. This study quantifies the impact of free-breathing on CSF parameters, offering valuable physiological references for understanding the effects of breathing on CSF dynamics.

CSF-to-blood toxins clearance is modulated by breathing through cranio-spinal CSF oscillation

Clearance of brain toxins occurs during sleep, although the mechanism remains unknown. Previous studies implied that the intracranial aqueductal cerebrospinal fluid (CSF) oscillations are involved, but no mechanism was suggested. The rationale for focusing on the aqueductal CSF oscillations is unclear. This study focuses on the cranio-spinal CSF oscillation and the factors that modulate this flow. We propose a mechanism where increased cranio-spinal CSF movements enhance CSF-to-blood metabolic waste clearance through the spinal CSF re-absorption sites. A recent study demonstrating that disturbed sleep impairs CSF-to-blood but not brain-to-CSF clearance, supports the fundamentals of our proposed mechanism. Eight healthy subjects underwent phase-contrast magnetic resonance imaging to quantify the effect of respiration on the cranio-spinal CSF oscillations. Maximal CSF volume displaced from the cranium to the spinal canal during each respiration and cardiac cycle were derived as measures of cranio-spinal CSF mixing level. Transition from normal to slow and abdominal breathing resulted in a 56% increase in the maximal displaced CSF volume. Maximal change in the arterial-venous blood volume, which is the driving force of the CSF oscillations, was increased by 41% during slow abdominal breathing. Cranio-spinal CSF oscillations are driven by the momentary difference between arterial inflow and venous outflow. Breathing modulates the CSF oscillation through changes in the venous outflow. The amount of toxins being transferred to the spinal canal during each respiratory cycle is significantly increased during slow and deeper abdominal breathing, which explains enhanced CSF-to-blood toxins clearance during slow-wave sleep and poor clearance during disrupted sleep.

An analytic model for the flow induced in syringomyelia cavities

A simple two-dimensional fluid-structure-interaction problem, involving viscous oscillatory flow in a channel separated by an elastic membrane from a fluid-filled slender cavity, is analyzed to shed light on the flow dynamics pertaining to syringomyelia, a neurological disorder characterized by the appearance of a large tubular cavity (syrinx) within the spinal cord. The focus is on configurations in which the velocity induced in the cavity, representing the syrinx, is comparable to that found in the channel, representing the subarachnoid space surrounding the spinal cord, both flows being coupled through a linear elastic equation describing the membrane deformation. An asymptotic analysis for small stroke lengths leads to closed-form expressions for the leading-order oscillatory flow, and also for the stationary flow associated with the first-order corrections, the latter involving a steady distribution of transmembrane pressure. The magnitude of the induced flow is found to depend strongly on the frequency, with the result that for channel flow rates of non-sinusoidal waveform, as those found in the spinal canal, higher harmonics can dominate the sloshing motion in the cavity, in agreement with previous in vivo observations. Under some conditions, the cycle-averaged transmembrane pressure, also showing a marked dependence on the frequency, changes sign on increasing the cavity transverse dimension (i.e. orthogonal to the cord axis), underscoring the importance of cavity size in connection with the underlying hydrodynamics. The analytic results presented here can be instrumental in guiding future numerical investigations, needed to clarify the pathogenesis of syringomyelia cavities.

Concomitant evaluation of cardiovascular and cerebrovascular controls via Geweke spectral causality to assess the propensity to postural syncope

The evaluation of propensity to postural syncope necessitates the concomitant characterization of the cardiovascular and cerebrovascular controls and a method capable of disentangling closed loop relationships and decomposing causal links in the frequency domain. We applied Geweke spectral causality (GSC) to assess cardiovascular control from heart period and systolic arterial pressure variability and cerebrovascular regulation from mean arterial pressure and mean cerebral blood velocity variability in 13 control subjects and 13 individuals prone to develop orthostatic syncope. Analysis was made at rest in supine position and during head-up tilt at 60°, well before observing presyncope signs. Two different linear model structures were compared, namely bivariate autoregressive and bivariate dynamic adjustment classes. We found that (i) GSC markers did not depend on the model structure; (ii) the concomitant assessment of cardiovascular and cerebrovascular controls was useful for a deeper comprehension of postural disturbances; (iii) orthostatic syncope appeared to be favored by the loss of a coordinated behavior between the baroreflex feedback and mechanical feedforward pathway in the frequency band typical of the baroreflex functioning during the postural challenge, and by a weak cerebral autoregulation as revealed by the increased strength of the pressure-to-flow link in the respiratory band. GSC applied to spontaneous cardiovascular and cerebrovascular oscillations is a promising tool for describing and monitoring disturbances associated with posture modification.

Real-time Quantification of in vivo cerebrospinal fluid velocity using fMRI inflow effect

In vivo estimation of cerebrospinal fluid (CSF) velocity is crucial for understanding the glymphatic system and its potential role in neurodegenerative disorders such as Alzheimer's disease and Parkinson's disease. Current cardiac or respiratory gated approaches, such as 4D flow MRI, cannot capture CSF movement in real time due to limited temporal resolution and in addition deteriorate in accuracy at low fluid velocities. Other techniques like real-time PC-MRI or time-spatial labeling inversion pulse are not limited by temporal averaging but have limited availability even in research settings. This study aims to quantify the inflow effect of dynamic CSF motion on functional magnetic resonance imaging (fMRI) for in vivo, real-time measurement of CSF flow velocity. We considered linear and nonlinear models of velocity waveforms and empirically fit them to fMRI data from a controlled flow experiment. To assess the utility of this methodology in human data, CSF flow velocities were computed from fMRI data acquired in eight healthy volunteers. Breath holding regimens were used to amplify CSF flow oscillations. Our experimental flow study revealed that CSF velocity is nonlinearly related to inflow effect-mediated signal increase and well estimated using an extension of a previous nonlinear framework. Using this relationship, we recovered velocity from in vivo fMRI signal, demonstrating the potential of our approach for estimating CSF flow velocity in the human brain. This novel method could serve as an alternative approach to quantifying slow flow velocities in real time, such as CSF flow in the ventricular system, thereby providing valuable insights into the glymphatic system's function and its implications for neurological disorders.

Intrathecal delivery of Macromolecules: Clinical status and emerging technologies

The proximity and association of cerebrospinal fluid (CSF) and the intrathecal (IT) space with deep targets in the central nervous system (CNS) parenchyma makes IT injection an attractive route of administration for brain drug delivery. However, the extent to which intrathecally administered macromolecules are effective in treating neurological diseases is a question of both clinical debate and technological interest. We present the biological, chemical, and physical properties of the intrathecal space that are relevant to drug absorption, distribution, metabolism, and elimination from CSF. We then analyze the evolution of IT drug delivery in clinical trials over the last 20 years. Our analysis revealed that the percentage of clinical trials assessing IT delivery for the delivery of biologics (i.e., macromolecules, cells) for treatment of chronic conditions (e.g., neurodegeneration, cancer, and metabolic diseases) has steadily increased. Clinical trials exploring cell or macromolecular delivery within the IT space have not evaluated engineering technologies, such as depots, particles, or other delivery systems. Recent pre-clinical studies have evaluated IT macromolecule delivery in small animals, postulating that delivery efficacy can be assisted by external medical devices, micro- or nanoparticles, bulk biomaterials, and viral vectors. Further studies are necessary to evaluate the extent to which engineering technologies and IT administration improve CNS targeting and therapeutic outcome.

The Relationship Between Imbalance Symptom and Cardiac Pulsation Induced Mechanical Strain in the Brainstem and Cerebellum for Chiari Malformation Type I

Chiari malformation Type I (CMI) is known to have an altered biomechanical environment for the brainstem and cerebellum; however, it is unclear whether these altered biomechanics play a role in the development of CMI symptoms. We hypothesized that CMI subjects have a higher cardiac-induced strain in specific neurological tracts pertaining to balance, and postural control. We measured displacement over the cardiac cycle using displacement encoding with stimulated echoes magnetic resonance imaging in the cerebellum, brainstem, and spinal cord in 37 CMI subjects and 25 controls. Based on these measurements, we computed strain, translation, and rotation in tracts related to balance. The global strain on all tracts was small (<1%) for CMI subject and controls. Strain was found to be nearly doubled in three tracts for CMI subjects compared to controls (p < 0.03). The maximum translation and rotation were ∼150 μm and ∼1 deg, respectively and 1.5-2 times greater in CMI compared to controls in four tracts (p < 0.005). There was no significant difference between strain, translation, and rotation on the analyzed tracts in CMI subjects with imbalance compared to those without imbalance. A moderate correlation was found between cerebellar tonsillar position and strain on three tracts. The lack of statistically significant difference between strain in CMI subjects with and without imbalance could imply that the magnitude of the observed cardiac-induced strain was too small to cause substantial damage to the tissue (<1%). Activities such as coughing, or Valsalva may produce a greater strain.

A mathematical framework for the dynamic interaction of pulsatile blood, brain, and cerebrospinal fluid

BACKGROUND

Shedding light on less-known aspects of intracranial fluid dynamics may be helpful to understand the hydrocephalus mechanism. The present study suggests a mathematical framework based on in vivo inputs to compare the dynamic interaction of pulsatile blood, brain, and cerebrospinal fluid (CSF) between the healthy subject and the hydrocephalus patient.

METHOD

The input data for the mathematical formulations was pulsatile blood velocity, which was measured using cine PC-MRI. Tube law was used to transfer the created deformation by blood pulsation in the vessel circumference to the brain domain. The pulsatile deformation of brain tissue with respect to time was calculated and considered to be inlet velocity in the CSF domain. The governing equations in all three domains were continuity, Navier-Stokes, and concentration. We used Darcy law with defined permeability and diffusivity values to define the material properties in the brain.

RESULTS

We validated the preciseness of the CSF velocity and pressure through the mathematical formulations with cine PC-MRI velocity, experimental ICP, and FSI simulated velocity and pressure. We used the analysis of dimensionless numbers including Reynolds, Womersley, Hartmann, and Peclet to evaluate the characteristics of the intracranial fluid flow. In the mid-systole phase of a cardiac cycle, CSF velocity had the maximum value and CSF pressure had the minimum value. The maximum and amplitude of CSF pressure, as well as CSF stroke volume, were calculated and compared between the healthy subject and the hydrocephalus patient.

CONCLUSION

The present in vivo-based mathematical framework has the potential to gain insight into the less-known points in the physiological function of intracranial fluid dynamics and the hydrocephalus mechanism.

CSF Flow and Spinal Cord Motion in Patients With Spontaneous Intracranial Hypotension: A Phase Contrast MRI Study

BACKGROUND AND OBJECTIVES

Spontaneous intracranial hypotension (SIH) is characterized by loss of CSF volume. We hypothesize that in this situation of low volume, a larger CSF flow and spinal cord motion at the upper spine can be measured by noninvasive phase contrast MRI.

METHODS

A prospective, age-, sex-, and body mass index (BMI)-matched controlled cohort study on patients with SIH presenting with spinal longitudinal extradural fluid collection (SLEC) was conducted from October 2021 to February 2022. Cardiac-gated 2D phase contrast MRI sequences were acquired at segment C2/C3, and C5/C6 for CSF flow, and spinal cord motion analysis. Data processing was fully automated. CSF flow and spinal cord motion were analyzed by peak-to-peak amplitude and total displacement per segment and heartbeat, respectively. Clinical data included age, height, BMI, duration of symptoms, Bern score according to Dobrocky et al., and type of the spinal CSF leak according to Schievink et al. Groups were compared via the Mann-Whitney U test; multiple linear regression analysis was performed to address possible relations.

RESULTS

Twenty patients with SIH and 40 healthy controls were analyzed; each group consisted of 70% women. Eleven patients with SIH presented with type 1 leak, 8 with type 2, and 1 was indeterminate. CSF flow per heartbeat was increased at C2/C3 (peak-to-peak amplitude 65.68 ± 18.3 vs 42.50 ± 9.8 mm/s, total displacement 14.32 ± 3.5 vs 9.75 ± 2.7 mm, p < 0.001, respectively). Craniocaudal spinal cord motion per heartbeat was larger at segment C2/C3 (peak-to-peak amplitude 7.30 ± 2.4 vs 5.82 ± 2.0 mm/s, total displacement 1.01 ± 0.4 vs 0.74 ± 0.4 mm, p = 0.006, respectively) and at segment C5/C6 (total displacement 1.41 ± 0.7 vs 0.97 ± 0.4 mm, p = 0.021).

DISCUSSION

SLEC-positive patients with SIH show higher CSF flow and higher spinal cord motion at the upper cervical spine. This increased craniocaudal motion of the spinal cord per heartbeat might produce increased mechanical strain on neural tissue and adherent structures, which may be a mechanism leading to cranial nerve dysfunction, neck pain, and stiffness in SIH. Noninvasive phase contrast MRI of CSF flow and spinal cord motion is a promising diagnostic tool in SIH.

TRIAL REGISTRATION INFORMATION

German Clinical Trials Register, identification number: DRKS00017351.

CLASSIFICATION OF EVIDENCE

This study provides Class III evidence that noninvasive phase contrast MRI of the upper spine identifies differences in CSF flow and spinal cord motion in patients with SIH compared with healthy controls.

Real-time phase contrast MRI versus conventional phase contrast MRI at different spatial resolutions and velocity encodings

PURPOSES

To compare the accuracy of real-time phase-contrast echo-planar MRI (EPI-PC) and conventional cine phase-contrast MRI (Conv-PC) and to assess the influence of spatial resolutions (pixel size) and velocity encoding on flow measurements obtained with the two sequences.

METHODS

Flow quantification was assessed using a pulsatile flow phantom (diameter: 9.5 mm; mean flow rate: 1150 mm³/s; mean flow velocity: 1.6 cm/s). Firstly, the accuracy of the EPI-PC was checked by comparing it with the flow rate in the calibrated phantom and the pulsation index from Conv-PC. Secondly, flow data from the two sequences were compared quantitatively as a function of the pixel size and the velocity encoding.

RESULTS

The mean percentage difference between the EPI-PC flow rate and calibrated phantom flow rate was -2.9 ± 2.1% (Mean ± SD). The pulsatility indices for EPI-PC and Conv-PC were respectively 0.64 and 0.59. In order to keep the flow rate measurement error within 10%, the ROI in Conv-PC had to contain at least 13 pixels, while the ROI in EPI-PC had to contain at least 9 pixels. Furthermore, Conv-PC had a higher velocity-to-noise ratio and could use a higher velocity encoding than EPI-PC (20 cm/s and 15 cm/s, respectively).

CONCLUSIONS

The result of this in vitro study confirmed the accuracy of EPI-PC, and found that EPI-PC can adapt to lower spatial resolutions, but is more sensitive to velocity encoding than Conv-PC.

Fluid dynamics in aging-related dementias

Recent human and animal model experimental studies revealed novel pathways for fluid movement, immune cell trafficking and metabolic waste clearance in CNS. These studies raise the intriguing possibility that the newly discovered pathways, including the glymphatic system, lymphatic meningeal vessels and skull-brain communication channels, are impaired in aging and neurovascular and neurodegenerative diseases associated with dementia, including Alzheimer's disease (AD) and AD-related dementia. We provide an overview of the glymphatic and dural meningeal lymphatic systems, review current methods and approaches used to study glymphatic flow in humans and animals, and discuss current evidence and controversies related to its role in CNS flow homeostasis under physiological and pathophysiological conditions. Non-invasive imaging approaches are needed to fully understand the mechanisms and pathways driving fluid movement in CNS and their roles across lifespan including healthy aging and aging-related dementia.

Measurements of cerebrospinal fluid production: a review of the limitations and advantages of current methodologies

Cerebrospinal fluid (CSF) is an essential and critical component of the central nervous system (CNS). According to the concept of the "third circulation" originally proposed by Cushing, CSF is mainly produced by the choroid plexus and subsequently leaves the cerebral ventricles via the foramen of Magendie and Luschka. CSF then fills the subarachnoid space from whence it disperses to all parts of the CNS, including the forebrain and spinal cord. CSF provides buoyancy to the submerged brain, thus protecting it against mechanical injury. CSF is also transported via the glymphatic pathway to reach deep interstitial brain regions along perivascular channels; this CSF clearance pathway promotes transport of energy metabolites and signaling molecules, and the clearance of metabolic waste. In particular, CSF is now intensively studied as a carrier for the removal of proteins implicated in neurodegeneration, such as amyloid-β and tau. Despite this key function of CSF, there is little information about its production rate, the factors controlling CSF production, and the impact of diseases on CSF flux. Therefore, we consider it to be a matter of paramount importance to quantify better the rate of CSF production, thereby obtaining a better understanding of CSF dynamics. To this end, we now review the existing methods developed to measure CSF production, including invasive, noninvasive, direct, and indirect methods, and MRI-based techniques. Depending on the methodology, estimates of CSF production rates in a given species can extend over a ten-fold range. Throughout this review, we interrogate the technical details of CSF measurement methods and discuss the consequences of minor experimental modifications on estimates of production rate. Our aim is to highlight the gaps in our knowledge and inspire the development of more accurate, reproducible, and less invasive techniques for quantitation of CSF production.

Blood and cerebrospinal fluid flow oscillations measured with real-time phase-contrast MRI: breathing mode matters

BACKGROUND

Cervical blood and cerebrospinal fluid (CSF) flow rates can be quantified with Phase-contrast (PC) MRI, which is routinely used for clinical studies. Previous MRI studies showed that venous and CSF flow alterations are linked to various pathological conditions. Since it is well known that, besides the heart beating, the thoracic pump influences the blood and CSF dynamics, we studied the effect of different respiration modes on blood and CSF flow rates using a real-time (RT)-PC prototype.

METHODS

Thirty healthy volunteers were examined with a 3 T scanner. A RT-PC sequence was acquired at the first cervical level to quantify the flow rates of internal carotid arteries, internal jugular veins (IJVs) and CSF. Each RT-PC acquisition was repeated three times, while the subjects were asked to breathe in three different ways for 60 s each: freely (F), with a constant rate (PN) and with deep and constant respiration rate (PD). The average flow rates were computed, they were removed from the respective signals and integrated in the inspiratory and expiratory phases (differential volumes). Finally, the power spectral density was computed for each detrended flow rate. High- and very-high frequency peaks were identified on the spectra while their frequencies were compared to the respiratory and cardiac frequencies estimated using a thoracic belt and a pulse oximeter. The area under the spectra was computed in four 0.5 Hz-wide ranges, centered on the high-frequency peak, on very-high frequency peak and its 2nd and 3rd harmonics, and then they were normalized by the flow rate variance. The effect of breathing patterns on average flow rates, on systolic and diastolic peaks, and on the normalized power was tested. Finally, the differential volumes of inspiration were compared to those of expiration.

RESULTS

The frequencies of the high- and very-high spectral peaks corresponded to the respiratory and cardiac frequencies. The average flow rate progressively decreased from F to PN to PD breathing, and the cardiac modulations were less predominant especially for the IJVs. The respiratory modulation increased with PD breathing. The average volumes displaced in the inspiratory phases were not significantly different from those of the expiratory one.

CONCLUSIONS

The spectral analyses demonstrated higher respiratory modulations in PD compared to free breathing, even prevailing the cardiac modulation in the IJVs, showing an increment of the thoracic pump affecting the flow rate shape.

Computational fluid dynamics model to predict the dynamical behavior of the cerebrospinal fluid through implementation of physiological boundary conditions

Cerebrospinal fluid (CSF) dynamics play an important role in maintaining a stable central nervous system environment and are influenced by different physiological processes. Multiple studies have investigated these processes but the impact of each of them on CSF flow is not well understood. A deeper insight into the CSF dynamics and the processes impacting them is crucial to better understand neurological disorders such as hydrocephalus, Chiari malformation, and intracranial hypertension. This study presents a 3D computational fluid dynamics (CFD) model which incorporates physiological processes as boundary conditions. CSF production and pulsatile arterial and venous volume changes are implemented as inlet boundary conditions. At the outlets, 2-element windkessel models are imposed to simulate CSF compliance and absorption. The total compliance is first tuned using a 0D model to obtain physiological pressure pulsations. Then, simulation results are compared with in vivo flow measurements in the spinal subarachnoid space (SAS) and cerebral aqueduct, and intracranial pressure values reported in the literature. Finally, the impact of the distribution of and total compliance on CSF pressures and velocities is evaluated. Without respiration effects, compliance of 0.17 ml/mmHg yielded pressure pulsations with an amplitude of 5 mmHg and an average value within the physiological range of 7-15 mmHg. Also, model flow rates were found to be in good agreement with reported values. However, when adding respiration effects, similar pressure amplitudes required an increase of compliance value to 0.51 ml/mmHg, which is within the range of 0.4-1.2 ml/mmHg measured in vivo. Moreover, altering the distribution of compliance over the four different outlets impacted the local flow, including the flow through the foramen magnum. The contribution of compliance to each outlet was directly proportional to the outflow at that outlet. Meanwhile, the value of total compliance impacted intracranial pressure. In conclusion, a computational model of the CSF has been developed that can simulate CSF pressures and velocities by incorporating boundary conditions based on physiological processes. By tuning these boundary conditions, we were able to obtain CSF pressures and flows within the physiological range.

Effect of Normal Breathing on the Movement of CSF in the Spinal Subarachnoid Space

BACKGROUND AND PURPOSE

Forced respirations reportedly have an effect on CSF movement in the spinal canal. We studied respiratory-related CSF motion during normal respiration.

MATERIALS AND METHODS

Six healthy subjects breathed at their normal rate with a visual guide to ensure an unchanging rhythm. Respiratory-gated phase-contrast MR flow images were acquired at 5 selected axial planes along the spine. At each spinal level, we computed the flow rate voxelwise in the spinal canal, together with the associated stroke volume. From these data, we computed the periodic volume changes of spinal segments. A phantom was used to quantify the effect of respiration-related magnetic susceptibility changes on the velocity data measured.

RESULTS

At each level, CSF moved cephalad during inhalation and caudad during expiration. While the general pattern of fluid movement was the same in the 6 subjects, the flow rates, stroke volumes, and spine segment volume changes varied among subjects. Peak flow rates ranged from 0.60 to 1.59 mL/s in the cervical region, 0.46 to 3.17 mL/s in the thoracic region, and 0.75 to 3.64 mL/s in the lumbar region. The differences in flow rates along the canal yielded cyclic volume variations of spine segments that were largest in the lumbar spine, ranging from 0.76 to 3.07 mL among subjects. In the phantom study, flow velocities oscillated periodically during the respiratory cycle by up to 0.02 cm/s or 0.5%.

CONCLUSIONS

Respiratory-gated measurements of the CSF motion in the spinal canal showed cyclic oscillatory movements of spinal fluid correlated to the breathing pattern.

Cerebrovascular activity is a major factor in the cerebrospinal fluid flow dynamics

Cerebrospinal fluid (CSF) provides physical protection to the central nervous system as well as an essential homeostatic environment for the normal functioning of neurons. Additionally, it has been proposed that the pulsatile movement of CSF may assist in glymphatic clearance of brain metabolic waste products implicated in neurodegeneration. In awake humans, CSF flow dynamics are thought to be driven primarily by cerebral blood volume fluctuations resulting from a number of mechanisms, including a passive vascular response to blood pressure variations associated with cardiac and respiratory cycles. Recent research has shown that mechanisms that rely on the action of vascular smooth muscle cells ("cerebrovascular activity") such as neuronal activity, changes in intravascular CO2, and autonomic activation from the brainstem, may lead to CSF pulsations as well. Nevertheless, the relative contribution of these mechanisms to CSF flow remains unclear. To investigate this further, we developed an MRI approach capable of disentangling and quantifying CSF flow components of different time scales associated with these mechanisms. This approach was evaluated on human control subjects (n = 12) performing intermittent voluntary deep inspirations, by determining peak flow velocities and displaced volumes between these mechanisms in the fourth ventricle. We found that peak flow velocities were similar between the different mechanisms, while displaced volumes per cycle were about a magnitude larger for deep inspirations. CSF flow velocity peaked at around 10.4 s (range 7.1-14.8 s, n = 12) following deep inspiration, consistent with known cerebrovascular activation delays for this autonomic challenge. These findings point to an important role of cerebrovascular activity in the genesis of CSF pulsations. Other regulatory triggers for cerebral blood flow such as autonomic arousal and orthostatic challenges may create major CSF pulsatile movement as well. Future quantitative comparison of these and possibly additional types of CSF pulsations with the proposed approach may help clarify the conditions that affect CSF flow dynamics.

Cerebrospinal fluid dynamics along the optic nerve

The cerebrospinal fluid (CSF) plays an important role in delivering nutrients and eliminating the metabolic wastes of the central nervous system. An interrupted CSF flow could cause disorders of the brain and eyes such as Alzheimer's disease and glaucoma. This review provides an overview of the anatomy and flow pathways of the CSF system with an emphasis on the optic nerve. Imaging technologies used for visualizing the CSF dynamics and the anatomic structures associated with CSF circulation have been highlighted. Recent advances in the use of computational models to predict CSF flow patterns have been introduced. Open questions and potential mechanisms underlying CSF circulation at the optic nerves have also been discussed.

Cardiac and Respiratory Influences on Intracranial and Neck Venous Flow, Estimated Using Real-Time Phase-Contrast MRI

The study of brain venous drainage has gained attention due to its hypothesized link with various neurological conditions. Intracranial and neck venous flow rate may be estimated using cardiac-gated cine phase-contrast (PC)-MRI. Although previous studies showed that breathing influences the neck's venous flow, this aspect could not be studied using the conventional segmented PC-MRI since it reconstructs a single cardiac cycle. The advent of real-time PC-MRI has overcome these limitations. Using this technique, we measured the internal jugular veins and superior sagittal sinus flow rates in a group of 16 healthy subjects (12 females, median age of 23 years). Comparing forced-breathing and free-breathing, the average flow rate decreased and the respiratory modulation increased. The flow rate decrement may be due to a vasoreactive response to deep breathing. The respiratory modulation increment is due to the thoracic pump's greater effect during forced breathing compared to free breathing. These results showed that the breathing mode influences the average blood flow and its pulsations. Since effective drainage is fundamental for brain health, rehabilitative studies might use the current setup to investigate if respiratory exercises positively affect clinical variables and venous drainage.

Treadmill locomotion in the American alligator (Alligator mississippiensis) produces dynamic changes in intracranial cerebrospinal fluid pressure

To examine the influence of movement on cerebrospinal fluid (CSF) dynamics, intracranial subdural pressure recordings were taken from sub-adult alligators (Alligator mississippiensis) locomoting on a treadmill. Pressure recordings documenting the cardiac, ventilatory, and barostatic influences on the CSF were in good agreement with previous studies. During locomotion the CSF exhibits sinusoidal patterns of pressure change that spanned a mean amplitude of 56 mm Hg, some 16 × the amplitude of the cardiac-linked pulsations. These sinusoidal CSF pulsations were closely linked to the locomotor kinematics, particularly the lateral oscillations of the alligator's head. Data recorded from the freely moving alligators suggest that fluid inertia, body cavity pressures, and likely other factors all influence the CSF pressure. The clear relationship between movement and CSF pressure described in this study suggests that the paucity of studies examining human CSF dynamics during movement should be addressed.

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.

Investigation of Human Intrathecal Solute Transport Dynamics Using a Novel in vitro Cerebrospinal Fluid System Analog

BACKGROUND

Understanding the relationship between cerebrospinal fluid (CSF) dynamics and intrathecal drug delivery (ITDD) injection parameters is essential to improve treatment of central nervous system (CNS) disorders.

METHODS

An anatomically detailed in vitro model of the complete CSF system was constructed. Patient-specific cardiac- and respiratory-induced CSF oscillations were input to the model in the subarachnoid space and within the ventricles. CSF production was input at the lateral ventricles and CSF absorption at the superior sagittal sinus. A model small molecule simulated drug product containing fluorescein was imaged within the system over a period of 3-h post-lumbar ITDD injections and used to quantify the impact of (a) bolus injection volume and rate, (b) post-injection flush volume, rate, and timing, (c) injection location, and (d) type of injection device. For each experiment, neuraxial distribution of fluorescein in terms of spatial temporal concentration, area-under-the-curve (AUC), and percent of injected dose (%ID) to the brain was quantified at a time point 3-h post-injection.

RESULTS

For all experiments conducted with ITDD administration in the lumbar spine, %ID to the brain did not exceed 11.6% at a time point 3-h post-injection. Addition of a 12 mL flush slightly increased solute transport to the brain up to +3.9%ID compared to without a flush (p < 0.01). Implantation of a lumbar catheter with the tip at an equivalent location to the lumbar placed needle, but with rostral tip orientation, resulted in a small improvement of 1.5%ID to the brain (p < 0.05). An increase of bolus volume from 5 to 20 mL improved solute transport to the brain from 5.0 to 6.3%ID, but this improvement was not statistically significant. Increasing bolus injection rate from 5 to 13.3 mL/min lacked improvement of solute transport to the brain, with a value of 6.3 compared to 5.7%ID.

CONCLUSION

The in vitro modeling approach allowed precisely controlled and repeatable parametric investigation of ITDD injection protocols and devices. In combination, the results predict that parametric changes in lumbar spine ITDD-injection related parameters and devices can alter %ID to the brain and be tuned to optimize therapeutic benefit to CNS targets.

Real-Time Phase-Contrast MRI to Monitor Cervical Blood and Cerebrospinal Fluid Flow Beat-by-Beat Variability

Beat-by-beat variability (BBV) rhythms are observed in both cardiovascular (CV) and intracranial (IC) compartments, yet interactions between the two are not fully understood. Real-Time Phase-Contrast (RT-PC) MRI sequence was acquired for 30 healthy volunteers at 1st cervical level on a 3T scanner. The arterial (AF), venous (VF), and cerebrospinal fluid (CSF) flow (CSFF) were computed as velocity integrals over the internal carotid artery, internal jugular vein, and CSF. AF, VF, and CSFF signals were segmented in inspiration and expiration beats, to assess the respiration influence. Systolic and diastolic BBV, and heart period series underwent autoregressive power spectral density analysis, to evaluate the low-frequency (LF, Mayer waves) and high frequency (HF, respiratory waves) components. The diastolic VF had the largest BBV. LF power was high in the diastolic AF series, poor in all CSFF series. The pulse wave analyses revealed higher mean amplitude during inspiration. Findings suggests a possible role of LF modulation of IC resistances and propagation of HF waves from VF to AF and CCSF. PC-RT-MRI could provide new insight into the interaction between CV and IC regulation and pave the way for a detailed analysis of the cerebrovascular effects of varied respiration patterns due to exercise and rehabilitation.

Use of real-time phase-contrast MRI to quantify the effect of spontaneous breathing on the cerebral arteries

Quantification of the effect of breathing on the cerebral circulation provides a better mechanistic understanding of the brain's circulatory system and is important in the early diagnosis of certain neurological diseases. However, conventional cine phase-contrast (CINE-PC) MRI cannot be used in this field of study because it only provides an average cardiac cycle flow curve reconstructed from multiple cardiac cycles. Unlike CINE-PC, phase-contrast echo-planar imaging (EPI-PC) can be used to quantify the blood flow rate in "real-time" and thus assess the effect of breathing on blood flow. Here, we first used post-processing software (developed in-house) to determine the feasibility of quantifying cerebral arterial blood flow with EPI-PC (relative to CINE-PC) in 16 participants. In a second step, we developed a new time-domain method for quantifying the intensity and the phase shift of the effects of breathing on the mean flow rate, stroke volume, cardiac period and amplitude of cerebral blood flow (in 10 participants). Our results showed that EPI-PC can quantify cerebral arterial blood flow rate with much the same degree of accuracy as CINE-PC but is more strongly influenced by differences in magnetic susceptibility. We found that breathing affected the mean flow rate, stroke volume and cardiac period of cerebral arterial blood flow.

Real‐time imaging of respiratory effects on cerebrospinal fluid flow in small diameter passageways

Purpose

Respiration‐related CSF flow through the cerebral aqueduct may be useful for elucidating physiology and pathophysiology of the glymphatic system, which has been proposed as a mechanism of brain waste clearance. Therefore, we aimed to (1) develop a real‐time (CSF) flow imaging method with high spatial and sufficient temporal resolution to capture respiratory effects, (2) validate the method in a phantom setup and numerical simulations, and (3) apply the method in vivo and quantify its repeatability and correlation with different respiratory conditions.

Methods

A golden‐angle radial flow sequence (reconstructed temporal resolution 168 ms, spatial resolution 0.6 mm) was implemented on a 7T MRI scanner and reconstructed using compressed sensing. A phantom setup mimicked simultaneous cardiac and respiratory flow oscillations. The effect of temporal resolution and vessel diameter was investigated numerically. Healthy volunteers (n = 10) were scanned at four different respiratory conditions, including repeat scans.

Results

Phantom data show that the developed sequence accurately quantifies respiratory oscillations (ratio real‐time/reference Q_R = 0.96 ± 0.02), but underestimates the rapid cardiac oscillations (ratio Q_C = 0.46 ± 0.14). Simulations suggest that Q_C can be improved by increasing temporal resolution. In vivo repeatability was moderate to very strong for cranial and caudal flow (intraclass correlation coefficient range: 0.55–0.99) and weak to strong for net flow (intraclass correlation coefficient range: 0.48–0.90). Net flow was influenced by respiratory condition (p < 0.01).

Conclusions

The presented real‐time flow MRI method can quantify respiratory‐related variations of CSF flow in the cerebral aqueduct, but it underestimates rapid cardiac oscillations. In vivo, the method showed good repeatability and a relationship between flow and respiration.

Evaluation of Cardiac- and Respiratory-driven Cerebrospinal Fluid Motions by Applying the S-transform to Steady-state Free Precession Phase Contrast Imaging

PURPOSE

To extract the status of hydrocephalus and other cerebrospinal fluid (CSF)-related diseases, a technique to characterize the cardiac- and respiratory-driven CSF motions separately under free breathing was developed. This technique is based on steady-state free precession phase contrast (SSFP-PC) imaging in combination with a Stockwell transform (S-transform).

METHODS

2D SSFP-PC at 3 T was applied to measure the CSF velocity in the caudal-cranial direction within a sagittal slice at the midline (N = 3) under 6-, 10-, and 16-s respiratory cycles and free breathing. The frequency-dependent window width of the S-transform was controlled by a particular scaling factor, which then converted the CSF velocity waveform into a spectrogram. Based on the frequency bands of the cardiac pulsation and respiration, as determined by the electrocardiogram (ECG) and respirator pressure sensors, Gaussian bandpass filters were applied to the CSF spectrogram to extract the time-domain cardiac- and respiratory-driven waveforms.

RESULTS

The cardiac-driven CSF velocity component appeared in the spectrogram clearly under all respiratory conditions. The respiratory-driven velocity under the controlled respiratory cycles was observed as constant frequency signals, compared to a time-varying frequency signal under free breathing. When the widow width was optimized using the scale factor, the temporal change in the respiratory-driven CSF component was even more apparent under free breathing.

CONCLUSION

Velocity amplitude variations and transient frequency changes of both cardiac- and respiratory-driven components were successfully characterized. These findings indicated that the proposed technique is useful for evaluating CSF motions driven by different cyclic forces.

Deep breathing couples CSF and venous flow dynamics

Venous system pathologies have increasingly been linked to clinically relevant disorders of CSF circulation whereas the exact coupling mechanisms still remain unknown. In this work, flow dynamics of both systems were studied using real-time phase-contrast flow MRI in 16 healthy subjects during normal and forced breathing. Flow evaluations in the aqueduct, at cervical level C3 and lumbar level L3 for both the CSF and venous fluid systems reveal temporal modulations by forced respiration. During normal breathing cardiac-related flow modulations prevailed, while forced breathing shifted the dominant frequency of both CSF and venous flow spectra towards the respiratory component and prompted a correlation between CSF and venous flow in the large vessels. The average of flow magnitude of CSF was increased during forced breathing at all spinal and intracranial positions. Venous flow in the large vessels of the upper body decreased and in the lower body increased during forced breathing. Deep respiration couples interdependent venous and brain fluid flow—most likely mediated by intrathoracic and intraabdominal pressure changes. Further insights into the driving forces of CSF and venous circulation and their correlation will facilitate our understanding how the venous system links to intracranial pressure regulation and of related forms of hydrocephalus.

The glymphatic hypothesis: the theory and the evidence

The glymphatic hypothesis proposes a mechanism for extravascular transport into and out of the brain of hydrophilic solutes unable to cross the blood-brain barrier. It suggests that there is a circulation of fluid carrying solutes inwards via periarterial routes, through the interstitium and outwards via perivenous routes. This review critically analyses the evidence surrounding the mechanisms involved in each of these stages. There is good evidence that both influx and efflux of solutes occur along periarterial routes but no evidence that the principal route of outflow is perivenous. Furthermore, periarterial inflow of fluid is unlikely to be adequate to provide the outflow that would be needed to account for solute efflux. A tenet of the hypothesis is that flow sweeps solutes through the parenchyma. However, the velocity of any possible circulatory flow within the interstitium is too small compared to diffusion to provide effective solute movement. By comparison the earlier classical hypothesis describing extravascular transport proposed fluid entry into the parenchyma across the blood-brain barrier, solute movements within the parenchyma by diffusion, and solute efflux partly by diffusion near brain surfaces and partly carried by flow along "preferred routes" including perivascular spaces, white matter tracts and subependymal spaces. It did not suggest fluid entry via periarterial routes. Evidence is still incomplete concerning the routes and fate of solutes leaving the brain. A large proportion of the solutes eliminated from the parenchyma go to lymph nodes before reaching blood but the proportions delivered directly to lymph or indirectly via CSF which then enters lymph are as yet unclear. In addition, still not understood is why and how the absence of AQP4 which is normally highly expressed on glial endfeet lining periarterial and perivenous routes reduces rates of solute elimination from the parenchyma and of solute delivery to it from remote sites of injection. Neither the glymphatic hypothesis nor the earlier classical hypothesis adequately explain how solutes and fluid move into, through and out of the brain parenchyma. Features of a more complete description are discussed. All aspects of extravascular transport require further study.

Human in silico trials for parametric computational fluid dynamics investigation of cerebrospinal fluid drug delivery: impact of injection location, injection protocol, and physiology

Background

Intrathecal drug delivery has a significant role in pain management and central nervous system (CNS) disease therapeutics. A fluid-physics based tool to assist clinicians in choosing specific drug doses to the spine or brain may help improve treatment schedules.

Methods

This study applied computational fluid dynamics (CFD) and in vitro model verification to assess intrathecal drug delivery in an anatomically idealized model of the human CSF system with key anatomic features of the CNS. Key parameters analyzed included the role of (a) injection location including lumbar puncture (LP), cisterna magna (CM) and intracerebroventricular (ICV), (b) LP injection rate, injection volume, and flush volume, (c) physiologic factors including cardiac-induced and deep respiration-induced CSF stroke volume increase. Simulations were conducted for 3-h post-injection and used to quantify spatial–temporal tracer concentration, regional area under the curve (AUC), time to maximum concentration (T_max), and maximum concentration (C_max), for each case.

Results

CM and ICV increased AUC to brain regions by ~ 2 logs compared to all other simulations. A 3X increase in bolus volume and addition of a 5 mL flush both increased intracranial AUC to the brain up to 2X compared to a baseline 5 mL LP injection. In contrast, a 5X increase in bolus rate (25 mL/min) did not improve tracer exposure to the brain. An increase in cardiac and respiratory CSF movement improved tracer spread to the brain, basal cistern, and cerebellum up to ~ 2 logs compared to the baseline LP injection.

Conclusion

The computational modeling approach provides ability to conduct in silico trials representative of CSF injection protocols. Taken together, the findings indicate a strong potential for delivery protocols to be optimized to reach a target region(s) of the spine and/or brain with a needed therapeutic dose. Parametric modification of bolus rate/volume and flush volume was found to have impact on tracer distribution; albeit to a smaller degree than injection location, with CM and ICV injections resulting in greater therapeutic dose to brain regions compared to LP. CSF stroke volume and frequency both played an important role and may potentially have a greater impact than the modest changes in LP injection protocols analyzed such as bolus rate, volume, and flush.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12987-022-00304-4.

A Brief Overview of the Cerebrospinal Fluid System and Its Implications for Brain and Spinal Cord Diseases

A comprehensive understanding of the cerebrospinal fluid (CSF) system is essential for our understanding of health and disease within the central nervous system (CNS). The system of CSF refers to all components involved in CSF production, movement, and absorption. In recent years, extensive research has resulted in vastly improved understanding of the CSF system in health and disease. Yet, several aspects remain to be fully clarified, notably along the spinal cord as the preponderance of research has focused on the brain. This review briefly summarizes the CSF system and its implications for CNS diseases and highlights the knowledge gaps that require further research.

Measuring Arterial Pulsatility With Dynamic Inflow Magnitude Contrast

The pulsatility of blood flow through cerebral arteries is clinically important, as it is intrinsically associated with cerebrovascular health. In this study we outline a new MRI approach to measuring the real-time pulsatile flow in cerebral arteries, which is based on the inflow phenomenon associated with fast gradient-recalled-echo acquisitions. Unlike traditional phase-contrast techniques, this new method, which we dub dynamic inflow magnitude contrast (DIMAC), does not require velocity-encoding gradients as sensitivity to flow velocity is derived purely from the inflow effect. We achieved this using a highly accelerated single slice EPI acquisition with a very short TR (15 ms) and a 90° flip angle, thus maximizing inflow contrast. We simulate the spoiled GRE signal in the presence of large arteries and perform a sensitivity analysis. The sensitivity analysis demonstrates that in the regime of high inflow contrast, DIMAC shows much greater sensitivity to flow velocity over blood volume changes. We support this theoretical prediction with in-vivo data collected in two separate experiments designed to demonstrate the utility of the DIMAC signal contrast. We perform a hypercapnia challenge experiment in order to experimentally modulate arterial tone within subjects, and thus modulate the arterial pulsatile flow waveform. We also perform a thigh-cuff release challenge, designed to induce a transient drop in blood pressure, and demonstrate that the continuous DIMAC signal captures the complex transient change in the pulsatile and non-pulsatile components of flow. In summary, this study proposes a new role for a well-established source of MR image contrast and demonstrates its potential for measuring both steady-state and dynamic changes in arterial tone.

Upright versus supine MRI: effects of body position on craniocervical CSF flow

Background

Cerebrospinal fluid (CSF) circulation between the brain and spinal canal, as part of the glymphatic system, provides homeostatic support to brain functions and waste clearance. Recently, it has been observed that CSF flow is strongly driven by cardiovascular brain pulsation, and affected by body orientation. The advancement of MRI has allowed for non-invasive examination of the CSF hydrodynamic properties. However, very few studies have addressed their relationship with body position (e.g., upright versus supine). It is important to understand how CSF hydrodynamics are altered by body position change in a single cardiac phase and how cumulative long hours staying in either upright or supine position can affect craniocervical CSF flow.

Methods

In this study, we investigate the changes in CSF flow at the craniocervical region with flow-sensitive MRI when subjects are moved from upright to supine position. 30 healthy volunteers were imaged in upright and supine positions using an upright MRI. The cranio-caudal and caudo-cranial CSF flow, velocity and stroke volume were measured at the C2 spinal level over one cardiac cycle using phase contrast MRI. Statistical analysis was performed to identify differences in CSF flow properties between the two positions.

Results

CSF stroke volume per cardiac cycle, representing CSF volume oscillating in and out of the cranium, was ~ 57.6% greater in supine (p < 0.0001), due to a ~ 83.8% increase in caudo-cranial CSF peak velocity during diastole (p < 0.0001) and extended systolic phase duration when moving from upright (0.25 ± 0.05 s) to supine (0.34 ± 0.08 s; p < 0.0001). Extrapolation to a 24 h timeframe showed significantly larger total CSF volume exchanged at C2 with 10 h spent supine versus only 5 h (p < 0.0001).

Conclusions

In summary, body position has significant effects on CSF flow in and out of the cranium, with more CSF oscillating in supine compared to upright position. Such difference was driven by an increased caudo-cranial diastolic CSF velocity and an increased systolic phase duration when moving from upright to supine position. Extrapolation to a 24 h timeframe suggests that more time spent in supine position increases total amount of CSF exchange, which may play a beneficial role in waste clearance in the brain.

Temperature measurement of intracranial cerebrospinal fluid using second-order motion compensation diffusion tensor imaging

To reduce the determination errors of CSF pulsation in diffusion-weighted image (DWI) thermometry, we investigated whether applying second-order motion compensation diffusion tensor imaging (2nd-MC DTI) and fractional anisotropy (FA) processing improves the measurement of intracranial cerebrospinal fluid (CSF) temperature. In a phantom study, we investigated the relationship between temperature and FA in artificial CSF (ACSF) to determine the threshold for FA processing. The calculated temperatures of ACSF were compared with those of water. In a human study, 18 healthy volunteers were scanned using conventional DTI (c-DTI) and 2nd-MC DTI on a 3.0 T magnetic resonance imaging (MRI) system. A temperature map was created using diffusion coefficients from each DWI with/without FA processing. The temperatures of intracranial CSF were compared between each DTI image using Welch's analysis of variance and Games-Howell's multiple comparisons. In the phantom study, FA did not exceed 0.1 at any temperature. Consequently, pixels exceeding the threshold of 0.1 were removed from the temperature map. Intracranial CSF temperatures significantly differed between the four methods (p < 0.0001). The lowest temperature was 2nd-MC DTI with FA processing (mean, 35.62 °C), followed in order by c-DTI with FA processing (mean, 36.16 °C), 2nd-MC DTI (mean, 37.08 °C), and c-DTI (mean, 39.08 °C;p < 0.01 for each). Because the calculated temperature of ACSF was estimated to be lower than that of water, the temperature of 2nd-DTI with FA processing was considered reasonable. The method of 2nd-MC DTI with FA processing enabled determining intracranial CSF temperature with a reduction in CSF pulsation.

Characterization of Cardiac- and Respiratory-driven Cerebrospinal Fluid Motions Using a Correlation Mapping Technique Based on Asynchronous Two-dimensional Phase Contrast MR Imaging

PURPOSE

The cardiac- and respiratory-driven components of cerebrospinal fluid (CSF) motion characteristics and bulk flow are not yet completely understood. Therefore, the present study aimed to characterize cardiac- and respiratory-driven CSF motions in the intracranial space using delay time, CSF velocity waveform correlation, and displacement.

METHODS

Asynchronous two-dimensional phase-contrast at 3T was applied to measure the CSF velocity in the inferior-superior direction in a sagittal slice at the midline (N = 12) and an axial slice at the foramen magnum (N = 8). Volunteers were instructed to engage in six-second respiratory cycles. The calculated delay time and correlation coefficients of the cardiac- and respiratory-driven velocity waveforms, separated in the frequency domain, were applied to evaluate the propagation of the CSF motion. The cardiac- and respiratory-driven components of the CSF displacement and motion volume were calculated during diastole and systole, and during inhalation and exhalation, respectively. The cardiac- and respiratory-driven components of the velocity, correlation, displacement, and motion volume were compared using an independent two-sample t-test.

RESULTS

The ratio of the cardiac-driven CSF velocity to the sum of the cardiac- and respiratory-driven CSF velocities was higher than the equivalent respiratory-driven ratio for all cases (P < 0.01). Delay time and correlation maps demonstrated that the cardiac-driven CSF motion propagated more extensively than the respiratory-driven CSF motion. The correlation coefficient of the cardiac-driven motion was significantly higher in the prepontine (P < 0.01), the aqueduct, and the fourth ventricle (P < 0.05). The respiratory-driven displacement and motion volume were significantly greater than the cardiac-driven equivalents for all observations (P < 0.01).

CONCLUSION

The correlation mapping technique characterized the cardiac- and respiratory-driven CSF velocities and their propagation properties in the intracranial space. Based on these findings, cardiac-driven CSF velocity is greater than respiratory-induced velocity, but the respiratory-driven velocity might displace farther.

Evaluating the Effect of Arterial Pulsation on Cerebrospinal Fluid Motion in the Sylvian Fissure of Patients with Middle Cerebral Artery Occlusion Using Low b-value Diffusion-weighted Imaging

PURPOSE

Decrease in signal of the cerebrospinal fluid (CSF) on low b-value diffusion weighted image (DWI) due to non-uniform flow can provide additional information regarding CSF motion. The purpose of the current study was to evaluate whether arterial pulsations constitute the driving force of CSF motion.

METHODS

We evaluated the CSF signals within the Sylvian fissure on low b-value DWI in 19 patients with unilateral middle cerebral artery (MCA) occlusion. DWI with b-value of 500 s/mm² was evaluated for a decrease in CSF signal within the Sylvian fissure including the Sylvian vallecula and lower, middle, and higher Sylvian fissures and graded as follows: the same as contralateral side; smaller signal decrease than that on contralateral side; and no signal decrease. MR angiography (MRA) findings of MCA were graded as follows: the same as contralateral, lower signal than contralateral signal, and no signal. In 15 patients, regional cerebral blood flow (rCBF) was evaluated using single-photon emission computed tomography (SPECT) studies and graded as >90%, 90%-70%, and <70% rCBF compared to contralateral. The correlations between the gradings were evaluated using G likelihood-ratio test.

RESULTS

There was no statistically significant correlation between the MRA and low b-value DWI gradings of CSF in all areas. There were statistically significant correlations between the decreases in CBF on SPECT and CSF signals in the middle Sylvian fissure.

CONCLUSION

The driving force of CSF pulsation in the Sylvian sinus may be related to the pulsations of the cerebral hemisphere rather than direct arterial pulsations.

Velocity vector reconstruction for real-time phase-contrast MRI with radial Maxwell correction

PURPOSE

To develop an auto-calibrated image reconstruction for highly accelerated multi-directional phase-contrast (PC) MRI that compensates for (1) reconstruction instabilities occurring for phase differences near ± π and (2) phase errors by concomitant magnetic fields that differ for individual radial spokes.

THEORY AND METHODS

A model-based image reconstruction for real-time PC MRI based on nonlinear inversion is extended to multi-directional flow by exploiting multiple flow-encodings for the estimation of velocity vectors. An initial smoothing constraint during iterative optimization is introduced to resolve the ambiguity of the solution space by penalizing phase wraps. Maxwell terms are considered as part of the signal model on a line-by-line basis to address phase errors by concomitant magnetic fields. The reconstruction methods are evaluated using simulated data and cross-sectional imaging of a rotating-disc, as well as in vivo for the aortic arch and cervical spinal canal at 3T.

RESULTS

Real-time three-directional velocity mapping in the aortic arch is achieved at 1.8 × 1.8 × 6 mm³ spatial and 60 ms temporal resolution. Artificial phase wraps are avoided in all cases using the smoothness constraint. Inter-spoke differences of concomitant magnetic fields are effectively compensated for by the model-based image reconstruction with integrated radial Maxwell correction.

CONCLUSION

Velocity vector reconstructions based on nonlinear inversion allow for high degrees of radial data undersampling paving the way for multi-directional PC MRI in real time. Whether a spoke-wise treatment of Maxwell terms is required or a computationally cheaper frame-wise approach depends on the individual application.

Non-invasive assessment of cerebrospinal fluid flow dynamics using phase-contrast magnetic resonance imaging in communicating hydrocephalus

This work aims to evaluate the changes in cerebrospinal fluid (CSF) hydrodynamics in patients diagnosed with communicating hydrocephalus. Besides, we establish the relationship between CSF flow dynamic parameters on the midbrain aqueduct and intracranial pressure (ICP). CSF hydrodynamics analysis was performed using Phase-Contrast Magnetic Resonance Imaging (PC-MRI) techniques on the midbrain aqueduct of 41 patients diagnosed with communicating hydrocephalus and 22 healthy volunteers. The correlation between CSF average flow in the midbrain aqueduct and intracranial pressure measured by Lumbar Puncture (LP) was assessed in patients with hydrocephalus. Pearson correlation coefficient was used to establish the correction between the average CSF flow of midbrain aqueduct and ICP. CSF dynamic parameters of the midbrain aqueduct in hydrocephalus patients, including peak positive velocity (7.348 cm/s), average velocity (0.623 cm/s), average flow (50.799 mm³/s), and regions of interest (ROI) area (9.978 mm²) were significantly higher than in the healthy controls (p < 0.05). This was after adjusting the age, gender, heart rate, systolic blood pressure, diastolic blood pressure, and body mass index. However, only the peak negative velocity of the midbrain aqueduct did not significantly differ between the groups (p = 0.209). A positive correlation was noted between the average flow (AF) of the midbrain aqueducts and ICP in hydrocephalus patients (y (AF) = 0.386× (ICP)-33.738, r = 0.787, p < 0.05). Reference data of CSF flow dynamic parameters was obtained through the PC-MRI in middle-aged healthy volunteers and communicating hydrocephalus patients. Although the sample size was constrained, this study has significant contributions. For instance, a significant correlation was noted between the average CSF flow of the aqueduct and ICP. This therefore provides a reference for clinicians to monitor ICP in patients with hydrocephalus.

Respiratory-driven Cyclic Cerebrospinal Fluid Motion in the Intracranial Cavity on Magnetic Resonance Imaging: Insights into the Pathophysiology of Neurofluid Dysfunction

Neurofluids, a recently developed term that refers to interstitial fluids in the parenchyma and cerebrospinal fluid (CSF) in the ventricle and subarachnoid space, play a role in draining waste products from the brain. Neurofluids have been implicated in pathological conditions such as Alzheimer’s disease and normal pressure hydrocephalus. Given that CSF moves faster in the CSF cavity than in the brain parenchyma, CSF motion can be detected by magnetic resonance imaging. CSF motion is synchronized to the heartbeat and respiratory cycle, but respiratory cycle-induced CSF motion has yet to be investigated in detail. Therefore, we analyzed CSF motion using dynamic improved motion-sensitized driven-equilibrium steady-state free precession-based analysis. We analyzed CSF motion linked to the respiratory cycle in four women and six men volunteers aged 23 to 38 years. We identified differences between free respiration and tasked respiratory cycle-associated CSF motion in the ventricles and subarachnoid space. Our results indicate that semi-quantitative analysis can be performed using the cranial site at which CSF motion is most prominent as a standard. Our findings may serve as a reference for elucidating the pathophysiology of diseases caused by abnormalities in neurofluids.

Achieving brain clearance and preventing neurodegenerative diseases-A glymphatic perspective

Age-related neurodegenerative diseases are a growing burden to society, and many are sporadic, meaning that the environment, diet and lifestyle play significant roles. Cerebrospinal fluid (CSF)-mediated clearing of brain waste products via perivascular pathways, named the glymphatic system, is receiving increasing interest, as it offers unexplored perspectives on understanding neurodegenerative diseases. The glymphatic system is involved in clearance of metabolic by-products such as amyloid-β from the brain, and its function is believed to lower the risk of developing some of the most common neurodegenerative diseases. Here, we present magnetic resonance imaging (MRI) data on the heart cycle's control of CSF flow in humans which corroborates findings from animal studies. We also review the importance of sleep, diet, vascular health for glymphatic clearance and find that these factors are also known players in brain longevity.

Hydrocephalus Revisited: New Insights into Dynamics of Neurofluids on Macro- and Microscales

New experimental and clinical findings question the historic view of hydrocephalus and its 100-year-old classification. In particular, real-time magnetic resonance imaging (MRI) evaluation of cerebrospinal fluid (CSF) flow and detailed insights into brain water regulation on the molecular scale indicate the existence of at least three main mechanisms that determine the dynamics of neurofluids: (1) inspiration is a major driving force; (2) adequate filling of brain ventricles by balanced CSF upsurge is sensed by cilia; and (3) the perivascular glial network connects the ependymal surface to the pericapillary Virchow-Robin spaces. Hitherto, these aspects have not been considered a common physiologic framework, improving knowledge and therapy for severe disorders of normal-pressure and posthemorrhagic hydrocephalus, spontaneous intracranial hypotension, and spaceflight disease.

Slithering CSF: Cerebrospinal Fluid Dynamics in the Stationary and Moving Viper Boa, Candoia aspera

Simple Summary

The cerebrospinal fluid (CSF) flows through and around the central nervous system to nourish, cleanse, and support the brain and spinal cord. Though abnormalities of this CSF flow have been linked to multiple human neural diseases, little is known about the underlying mechanics of CSF flow. This study was designed to test the hypothesis that movement of the body’s trunk could cause CSF flow; hence, the study was conducted on a snake, an animal with prominent trunk movement. The results demonstrate that the resting snake has a CSF pressure profile that is very similar to what is seen in humans and other mammals, and that the CSF dynamics are changed during either artificial (manual) or natural (locomotor) movement of the snake’s body

Abstract

In the viper boa (Candoia aspera), the cerebrospinal fluid (CSF) shows two stable overlapping patterns of pulsations: low-frequency (0.08 Hz) pulses with a mean amplitude of 4.1 mmHg that correspond to the ventilatory cycle, and higher-frequency (0.66 Hz) pulses with a mean amplitude of 1.2 mmHg that correspond to the cardiac cycle. Manual oscillations of anesthetized C. aspera induced propagating sinusoidal body waves. These waves resulted in a different pattern of CSF pulsations with frequencies corresponding to the displacement frequency of the body and with amplitudes greater than those of the cardiac or ventilatory cycles. After recovery from anesthesia, the snakes moved independently using lateral undulation and concertina locomotion. The episodes of lateral undulation produced similar influences on the CSF pressure as were observed during the manual oscillations, though the induced CSF pulsations were of lower amplitude during lateral undulation. No impact on the CSF was found while C. aspera was performing concertina locomotion. The relationship between the propagation of the body and the CSF pulsations suggests that the body movements produce an impulse on the spinal CSF.

Fluid metabolic pathways after subarachnoid hemorrhage

Aneurysmal subarachnoid hemorrhage (aSAH) is a devastating cerebrovascular disease with high mortality and morbidity. In recent years, a large number of studies have focused on the mechanism of early brain injury (EBI) and delayed cerebral ischemia (DCI), including vasospasm, neurotoxicity of hematoma and neuroinflammatory storm, after aSAH. Despite considerable efforts, no novel drugs have significantly improved the prognosis of patients in phase III clinical trials, indicating the need to further re-examine the multifactorial pathophysiological process that occurs after aSAH. The complex pathogenesis is reflected by the destruction of the dynamic balance of the energy metabolism in the nervous system after aSAH, which prevents the maintenance of normal neural function. This review focuses on the fluid metabolic pathways of the central nervous system (CNS), starting with ruptured aneurysms, and discusses the dysfunction of blood circulation, cerebrospinal fluid (CSF) circulation and the glymphatic system during disease progression. It also proposes a hypothesis on the metabolic disorder mechanism and potential therapeutic targets for aSAH patients.

Harnessing cerebrospinal fluid circulation for drug delivery to brain tissues

Initially thought to be useful only to reach tissues in the immediate vicinity of the CSF circulatory system, CSF circulation is now increasingly viewed as a viable pathway to deliver certain therapeutics deeper into brain tissues. There is emerging evidence that this goal is achievable in the case of large therapeutic proteins, provided conditions are met that are described herein. We show how fluid dynamic modeling helps predict infusion rate and duration to overcome high CSF turnover. We posit that despite model limitations and controversies, fluid dynamic models, pharmacokinetic models, preclinical testing, and a qualitative understanding of the glymphatic system circulation can be used to estimate drug penetration in brain tissues. Lastly, in addition to highlighting landmark scientific and medical literature, we provide practical advice on formulation development, device selection, and pharmacokinetic modeling. Our review of clinical studies suggests a growing interest for intra-CSF delivery, particularly for targeted proteins.