Relationship Between Cerebrospinal Fluid and Blood Dynamics in Healthy Volunteers and Patients with Communicating Hydrocephalus

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.

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.

Impact of Shunt Placement on CSF Dynamics

BACKGROUND

CSF dynamics are disturbed in chronic hydrocephalus (NPH). We hypothesise that these alterations reflect a disturbance of intracranial compliance. The aim of our study is to investigate the variations in intracranial hydrodynamics in NPH after ventricular shunt surgery.

PATIENTS AND METHOD

We included 14 patients with definite NPH. All patients improved after ventriculoperitoneal shunting. The patients underwent an analysis of intracranial haemodynamics by phase-contrast MRI (pcMRI) preoperatively, at 6 months postoperatively, and at 1 year postoperatively. We analysed the dynamics of intraventricular CSF at the level of the aqueduct of Sylvius (SVAQU) and CSF at the level of the high cervical subarachnoid spaces (SVCERV). We calculated the ratio between SVAQU and SVCERV, called CSFRATIO, which reflects the participation of intraventricular pulsatility in overall intracranial CSF pulsatility.

RESULTS

SVAQU significantly (p = 0.003) decreased from 240 ± 114 μL/cc to 214 ± 157 μL/cc 6 months after shunt placement. Six months after shunt placement, SVCERV significantly (p = 0.007) decreased from 627 ± 229 μL/cc to 557 ± 234 μL/cc. Twelve months after shunt placement, SVCERV continued to significantly (p = 0.001) decrease to 496 ± 234 μL/cc. CSFRATIO was not changed by surgery.

CONCLUSIONS

CSF dynamics are altered by shunt placement and might be a useful marker of the shunt's effectiveness-especially if pressure values start to rise again. The detection of changes in CSF dynamics would require a reference postoperative pcMRI measurement for each patient.

Cerebral microcirculatory pulse wave propagation and pulse wave amplitude mapping in retrospectively gated MRI

To analyze cerebral arteriovenous pulse propagation and to generate phase-resolved pulse amplitude maps from a fast gradient-echo sequence offering flow-related enhancement (FREE). Brain MRI was performed using a balanced steady-state free precession sequence at 3T followed by retrospective k-space gating. The time interval of the pulse wave between anterior-, middle- and posterior cerebral artery territories and the superior sagittal sinus were calculated and compared between and older and younger groups within 24 healthy volunteers. Pulse amplitude maps were generated and compared to pseudo-Continuous Arterial Spin Labeling (pCASL) MRI maps by voxel-wise Pearson correlation, Sørensen-Dice maps and in regards to signal contrast. The arteriovenous delays between all vascular territories and the superior sagittal sinus were significantly shorter in the older age group (11 individuals, ≥ 31 years) ranging between 169 ± 112 and 246 ± 299 ms versus 286 ± 244 to 419 ± 299 ms in the younger age group (13 individuals) (P ≤ 0.04). The voxel-wise pulse wave amplitude values and perfusion-weighted pCASL values correlated significantly (Pearson-r = 0.33, P < 0.01). Mean Dice overlaps of high (gray) and low (white matter) regions were 73 ± 3% and 59 ± 5%. No differences in image contrast were seen in the whole brain and the white matter, but significantly higher mean contrast of 0.73 ± 0.23% in cortical gray matter in FREE-MRI compared to 0.52 ± 0.12% in pCASL-MRI (P = 0.01). The dynamic information of flow-related enhancement allows analysis of the cerebral pulse wave propagation potentially providing information about the (micro)circulation on a regional level. However, the pulse wave amplitude reveals weaknesses in comparison to true perfusion-weighting and could rather be used to calculate a pulsatility index.

Hydrodynamic and Hemodynamic Interactions in Chronic Hydrocephalus

BACKGROUND

During a cardiac cycle, intracranial pressure is related to arterial entry into the cranium and its interaction with intracranial compliance. The arterial inflow is compensated by intracranial compliance and, initially, the flushing of cerebrospinal fluid (CSF) into the cervical subarachnoid spaces. Our objective is to analyze the interactions between intracranial arteriovenous exchange and cerebrospinal fluid oscillations.

METHOD

A total of 23 patients (73 ± 8 years) with suspected chronic hydrocephalus (CH) underwent an infusion test and phase-contrast MRI. Rout is an important factor in the diagnosis of CH. Patients were divided into 2 populations: probableCH (Rout: resistance to CSF outflow) (Rout > 12 mmHg/mL/min, 13 patients) and unlikelyCH (Rout < 12 mmHg/mL/min, 10 patients). We measured the intracranial vascular volume (arteriovenous stroke volume: SVvasc) and CSF (CSF stroke volume at upper cervical level: SVCSF) volume variations during the cardiac cycle.

RESULTS

In the whole population, we observed a significant correlation between SVvasc and SVCSF (R2 = 0.43; p = 0.0007). In the population unlikelyCH, this correlation was significant (R2 = 0.76; p = 0.001). In the population probableCH, this correlation was not significant (R2 = 0.17, p = 0.16).

CONCLUSIONS

These results show that the link between the compliance of the oscillating CSF and the abrupt arterial inflow seems to be altered in CH. CSF oscillations between intracranial and cervical fluid spaces limit the impact of the abrupt arterial inflow.

Managing Idiopathic Normal Pressure Hydrocephalus: Need for a Change of Mindset

Idiopathic normal pressure hydrocephalus (iNPH) refers to a complex brain disorder characterized by ventricular enlargement and the classic Hakim's triad of gait and balance difficulties, urinary incontinence, and cognitive impairment. It predominantly affects older patients in the absence of an identified cause. As the elderly population continues to increase, iNPH becomes a growing concern in the complex spectrum of neuro-geriatric care, with significant socio-economic implications. However, unlike other well-structured management approaches for neurodegenerative disorders, the management of iNPH remains largely uncodified, leading to suboptimal care in many cases. In this article, we highlighted the challenges of current practice and identify key points for an optimal structuration of care for iNPH. Adopting a global approach to iNPH could facilitate a progressive shift in mindset, moving away from solely aiming to cure an isolated neurological disease with uncertain outcomes to providing comprehensive care that focuses on improving the daily life of frail patients with complex neurodegenerative burdens, using tailored goals.

Cerebral arterial flow dynamics during systole and diastole phases in young and older healthy adults

BACKGROUND

Since arterial flow is the leading actor in neuro-fluids flow dynamics, it might be interesting to assess whether it is meaningful to study the arterial flow waveform in more detail and whether this provides new important information. Few studies have focused on determining the influence of heart rate variation over time on the arterial flow curve. Therefore, this study aimed to evaluate cerebral arterial flow waveforms at extracranial and intracranial compartments in young and elderly healthy adults, also considering systole and diastole phases.

METHODS

Cine phase-contrast magnetic resonance imaging (CINE-PC MRI) was performed on twenty-eight healthy young volunteers (HYV) and twenty healthy elderly volunteers (HEV) to measure arterial blood flows at the extracranial and intracranial planes. A semi-automated protocol using MATLAB scripts was implemented to identify the main representative points in the arterial flow waveforms. Representative arterial profiles were estimated for each group. Moreover, the effects of age and sex on flow times, amplitude-related parameters, and parameters related to systole and diastole phases were evaluated at the extracranial and intracranial compartments. Student's t-test or Wilcoxon's test (depending on the normality of the distribution) was used to detect significant differences.

RESULTS

In HYVs, significant differences were observed between extracranial and intracranial levels in parameters related to the AP1 amplitude. Besides the detected differences in pulsatility index (extracranial: 0.92 ± 0.20 vs. 1.28 ± 0.33; intracranial: 0.79 ± 0.15 vs. 1.14 ± 0.18, p < .001) and average flow (715 ± 136 vs. 607 ± 125 ml/min, p = .008) between HYV and HEV, differences in the amplitude value of the arterial flow profile feature points were also noted. Contrary to systole duration (HYV: 360 ± 29 ms; HEV: 364 ± 47 ms), diastole duration presented higher inter-individual variability in both populations (HYV: 472 ± 145 ms; HEV: 456 ± 106 ms). Our results also showed that, with age, it is mainly the diastolic phase that changes. Although no significant differences in duration were observed between the two populations, the mean flow value in the diastolic phase was significantly lower in HEV (extracranial: 628 ± 128 vs. 457 ± 111 ml/min; intracranial: 599 ± 121 vs. 473 ± 100 ml/min, p < .001). No significant differences were observed in the arterial flow parameters evaluated between females and males in either HYV or HEV.

CONCLUSION

Our study provides a novel contribution on the influence of the cardiac cycle phases on cerebral arterial flow. The main contribution in this study concerns the identification of age-related alterations in cerebral blood flow, which occur mainly during the diastolic phase. Specifically, we observed that mean flow significantly decreases with age during diastole, whereas mean flow during systole is consistent.

Hemodynamic and Hydrodynamic Pathophysiology in Chiari Type 1 Malformations: Towards Understanding the Genesis of Syrinx

BACKGROUND

The pathophysiology of this association of type 1 Chiari malformation (CM1) and syrinxes is still unknown. There is an alteration in the dynamics of neurofluids (cerebrospinal fluid, arterial and venous blood) during the cardiac cycle in CM1. Our objective is to quantify CSF or arterial blood or venous blood flow in patients with Chiari syndrome (CS) with and without syrinxes using phase-contrast MRI (PCMRI).

METHODS

We included 28 patients with CM1 (9 with syrinxes, 19 without). Morphological MRI with complementary PCMRI sequences was performed. We analyzed intraventricular CSF, subarachnoid spaces CSF, blood, and tonsillar pulsatility.

RESULTS

There is a highly significant correlation (p < 0.001) between cerebral blood flow, cerebral vascular expansion volume and venous drainage distribution. Venous drainage distribution is significantly inversely correlated with oscillatory CSF volume at the level of the foramen magnum plane [-0.37 (0.04)] and not significantly correlated at the C2C3 level [-0.37 (0.05)] over our entire population. This correlation maintained the same trend in patients with syrinxes [-0.80 (<0.01)] and disappeared in patients without a syrinx [-0.05 (0.81)].

CONCLUSION

The distribution of venous drainage is an important factor in intracranial homeostasis. Impaired venous drainage would lead to greater involvement of the CSF in compensating for arterial blood influx, thus contributing to syrinx genesis.

Cerebrospinal fluid and blood profiles of transfer RNA fragments show age, sex, and Parkinson's disease-related changes

Transfer RNA fragments (tRFs) have recently been shown to be an important family of small regulatory RNAs with diverse functions. Recent reports have revealed modified tRF blood levels in a number of nervous system conditions including epilepsy, ischemic stroke, and neurodegenerative diseases, but little is known about tRF levels in the cerebrospinal fluid (CSF). To address this issue, we studied age, sex, and Parkinson's disease (PD) effects on the distributions of tRFs in the CSF and blood data of healthy controls and PD patients from the NIH and the Parkinson's Progression Markers Initiative (PPMI) small RNA-seq datasets. We discovered that long tRFs are expressed in higher levels in the CSF than in the blood. Furthermore, the CSF showed a pronounced age-associated decline in the level of tRFs cleaved from the 3'-end and anti-codon loop of the parental tRNA (3'-tRFs, i-tRFs), and more pronounced profile differences than the blood profiles between the sexes. In comparison, we observed moderate age-related elevation of blood 3'-tRF levels. In addition, distinct sets of tRFs in the CSF and in the blood segregated PD patients from controls. Finally, we found enrichment of tRFs predicted to target cholinergic mRNAs (Cholino-tRFs) among mitochondrial-originated tRFs, raising the possibility that the neurodegeneration-related mitochondrial impairment in PD patients may lead to deregulation of their cholinergic tone. Our findings demonstrate that the CSF and blood tRF profiles are distinct and that the CSF tRF profiles are modified in a sex-, age-, and disease-related manner, suggesting that they reflect the inter-individual cerebral differences and calling for incorporating this important subset of small RNA regulators into future studies.

Characterising spinal cerebrospinal fluid flow in the pig with phase-contrast magnetic resonance imaging

BACKGROUND

Detecting changes in pulsatile cerebrospinal fluid (CSF) flow may assist clinical management decisions, but spinal CSF flow is relatively understudied. Traumatic spinal cord injuries (SCI) often cause spinal cord swelling and subarachnoid space (SAS) obstruction, potentially causing pulsatile CSF flow changes. Pigs are emerging as a favoured large animal SCI model; therefore, the aim of this study was to characterise CSF flow along the healthy pig spine.

METHODS

Phase-contrast magnetic resonance images (PC-MRI), retrospectively cardiac gated, were acquired for fourteen laterally recumbent, anaesthetised and ventilated, female domestic pigs (22-29 kg). Axial images were obtained at C2/C3, T8/T9, T11/T12 and L1/L2. Dorsal and ventral SAS regions of interest (ROI) were manually segmented. CSF flow and velocity were determined throughout a cardiac cycle. Linear mixed-effects models, with post-hoc comparisons, were used to identify differences in peak systolic/diastolic flow, and maximum velocity (cranial/caudal), across spinal levels and dorsal/ventral SAS. Velocity wave speed from C2/C3 to L1/L2 was calculated.

RESULTS

PC-MRI data were obtained for 11/14 animals. Pulsatile CSF flow was observed at all spinal levels. Peak systolic flow was greater at C2/C3 (dorsal: - 0.32 ± 0.14 mL/s, ventral: - 0.15 ± 0.13 mL/s) than T8/T9 dorsally (- 0.04 ± 0.03 mL/s; p < 0.001), but not different ventrally (- 0.08 ± 0.08 mL/s; p = 0.275), and no difference between thoracolumbar levels (p > 0.05). Peak diastolic flow was greater at C2/C3 (0.29 ± 0.08 mL/s) compared to T8/T9 (0.03 ± 0.03 mL/s, p < 0.001) dorsally, but not different ventrally (p = 1.000). Cranial and caudal maximum velocity at C2/C3 were greater than thoracolumbar levels dorsally (p < 0.001), and T8/T9 and L1/L2 ventrally (p = 0.022). Diastolic velocity wave speed was 1.41 ± 0.39 m/s dorsally and 1.22 ± 0.21 m/s ventrally, and systolic velocity wave speed was 1.02 ± 0.25 m/s dorsally and 0.91 ± 0.22 m/s ventrally.

CONCLUSIONS

In anaesthetised and ventilated domestic pigs, spinal CSF has lower pulsatile flow and slower velocity wave propagation, compared to humans. This study provides baseline CSF flow at spinal levels relevant for future SCI research in this animal model.

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.

A New Definition for Intracranial Compliance to Evaluate Adult Hydrocephalus After Shunting

The clinical application of intracranial compliance (ICC), ∆V/∆P, as one of the most critical indexes for hydrocephalus evaluation was demonstrated previously. We suggest a new definition for the concept of ICC (long-term ICC) where there is a longer amount of elapsed time (up to 18 months after shunting) between the measurement of two values (V1 and V2 or P1 and P2). The head images of 15 adult patients with communicating hydrocephalus were provided with nine sets of imaging in nine stages: prior to shunting, and 1, 2, 3, 6, 9, 12, 15, and 18 months after shunting. In addition to measuring CSF volume (CSFV) in each stage, intracranial pressure (ICP) was also calculated using fluid–structure interaction simulation for the noninvasive calculation of ICC. Despite small increases in the brain volume (16.9%), there were considerable decreases in the ICP (70.4%) and CSFV (80.0%) of hydrocephalus patients after 18 months of shunting. The changes in CSFV, brain volume, and ICP values reached a stable condition 12, 15, and 6 months after shunting, respectively. The results showed that the brain tissue needs approximately two months to adapt itself to the fast and significant ICP reduction due to shunting. This may be related to the effect of the “viscous” component of brain tissue. The ICC trend between pre-shunting and the first month of shunting was descending for all patients with a “mean value” of 14.75 ± 0.6 ml/cm H2O. ICC changes in the other stages were oscillatory (nonuniform). Our noninvasive long-term ICC calculations showed a nonmonotonic trend in the CSFV–ICP graph, the lack of a linear relationship between ICC and ICP, and an oscillatory increase in ICC values during shunt treatment. The oscillatory changes in long-term ICC may reflect the clinical variations in hydrocephalus patients after shunting.

Intrathecal Contrast-Enhanced Magnetic Resonance Imaging of Cerebrospinal Fluid Dynamics and Glymphatic Enhancement in Idiopathic Normal Pressure Hydrocephalus

Idiopathic normal pressure hydrocephalus (iNPH) is a neurodegenerative disease, characterized by cerebrospinal fluid (CSF) flow disturbance. Today, the only available treatment is CSF diversion surgery (shunt surgery). While traditional imaging biomarkers typically assess CSF space anatomy, recently introduced imaging biomarkers of CSF dynamics and glymphatic enhancement, provide imaging of CSF dynamics and thereby more specifically reveal elements of the underlying pathophysiology. The biomarkers address CSF ventricular reflux grade as well as glymphatic enhancement and derive from intrathecal contrast-enhanced MRI. However, the contrast agent serving as CSF tracer is administered off-label. In medicine, the introduction of new diagnostic or therapeutic methods must consider the balance between risk and benefit. To this end, we performed a prospective observational study of 95 patients with iNPH, comparing different intrathecal doses of the MRI contrast agent gadobutrol (0.10, 0.25, and 0.50 mmol, respectively), aiming at the lowest reasonable dose needed to retrieve diagnostic information about the novel MRI biomarkers. The present observations disclosed a dose-dependent enrichment of subarachnoid CSF spaces (cisterna magna, vertex, and velum interpositum) with dose-dependent ventricular reflux of tracer in iNPH, as well as dose-dependent glymphatic tracer enrichment. The association between tracer enrichment in CSF and parenchymal compartments were as well dose-related. Intrathecal gadobutrol in a dose of 0.25 mmol, but not 0.10 mmol, was at 1.5T MRI considered sufficient for imaging altered CSF dynamics and glymphatic enhancement in iNPH, even though 3T MRI provided better sensitivity. Tracer enrichment in CSF at the vertex and within the cerebral cortex and subcortical white matter was deemed too low for maintaining diagnostic information from a dose of 0.10 mmol. We conclude that reducing the intrathecal dose of gadobutrol from 0.50 to 0.25 mmol gadobutrol improves the safety margin while maintaining the necessary diagnostic information about disturbed CSF homeostasis and glymphatic failure in iNPH.

Analysis of relative changes in pulse shapes of intracranial pressure and cerebral blood flow velocity

Objective.Analysis of relative changes in the shapes of pulse waveforms of intracranial pressure (ICP) and transcranial Doppler cerebral blood flow velocity (CBFV) may provide information on intracranial compliance. We tested this hypothesis, introducing an index named the ratio of pulse slopes (RPS) that is based on inclinations of the ascending parts of the ICP and CBFV pulse waveforms. It has hypothetically a simple interpretation: a value of 1 indicates good compliance and a value less than 1, reduced compliance. Here, we investigated the usefulness of RPS for assessment of intracranial compliance.Approach.ICP and CBFV signals recorded simultaneously in 30 normal-pressure hydrocephalus patients during infusion tests were retrospectively analysed. CBFV was measured in the middle cerebral artery. Changes in RPS during the test were compared with changes in the height ratio of the first and second peak of the ICP pulse (P1/P2) and the shape of the ICP pulse was classified from normal (1) to pathological (4). Values are medians (lower, upper quartiles).Main results.There was a significant correlation between baseline RPS and intracranial elasticity (R = -0.55,p = 0.0018). During the infusion tests, both RPS and P1/P2 decreased with rising ICP [RPS, 0.80 (0.56, 0.92) versus 0.63 (0.44, 0.80),p = 0.00015; P1/P2, 0.58 (0.50, 0.91) versus 0.52 (0.36, 0.71),p = 0.00009] while the ICP pulses became more pathological in shape [class: 3 (2, 3) versus 3 (3, 4),p = 0.04]. The magnitude of the decrease in RPS during infusion was inversely correlated with baseline P1/P2 (R = -0.40,p < 0.03).Significance.During infusion, the slopes of the ascending parts of ICP and CBFV pulses become increasingly divergent with a shift in opposite directions. RPS seems to be a promising methodological tool for monitoring intracranial compliance with no additional volumetric manipulation required.

White Matter Microstructural Damage Associated With Gait Abnormalities in Idiopathic Normal Pressure Hydrocephalus

Background: Idiopathic normal pressure hydrocephalus (iNPH) is a common disease in elderly adults. Patients with iNPH are generally characterized by progressive gait impairment, cognitive deficits, and urinary urgency and/or incontinence. A number of radiographic studies have shown that iNPH patients have enlarged ventricles and altered brain morphology; however, few studies have focused on the relationships between altered brain structure and gait dysfunction due to iNPH. Thus, this study aimed to evaluate the abnormalities of white matter (WM) correlated with gait impairment in iNPH patients and to gain a better understanding of its underlying pathology.

Methods: Fifteen iNPH patients (five women, 10 men) were enrolled in this study, and each patient’s demographic and gait indices were collected. First, we performed a correlation analysis between the demographic and gait indices. Then, all gait indices were grouped according to the number of WM hyperintensities (WMH) among each WM tract (JHU WM tractography atlas), to perform comparative analysis.

Results: Considering sex and illness duration as covariates, correlation analysis showed a significantly negative correlation between step length (r = −0.80, p = 0.001), pace (r = −0.84, p = 2.96e-4), and age. After removing the effects of age, sex, and illness duration, correlation analysis showed negative correlation between step length (r = −0.73, p = 0.007), pace (r = −0.74, p = 0.005), and clinical-grade score and positive correlation between 3-m round trip time (r = 0.66, p = 0.021), rising time (r = 0.76, p = 0.004), and clinical-grade score. Based on WMH of each white matter tract, gait indices showed significant differences (p < 0.05/48, corrected by Bonferroni) between fewer WMH patients and more WMH in the middle cerebellar peduncle, left medial lemniscus, left posterior limb of the internal capsule (IC), and right posterior limb of the IC.

Conclusions: Our results indicated that iNPH patients exhibited gait-related WM abnormalities located in motor and sensory pathways around the ventricle, which is beneficial to understand the underlying pathology of iNPH.

Transmantle Pressure Computed from MR Imaging Measurements of Aqueduct Flow and Dimensions

BACKGROUND AND PURPOSE

Measuring transmantle pressure, the instantaneous pressure difference between the lateral ventricles and the cranial subarachnoid space, by intracranial pressure sensors has limitations. The aim of this study was to compute transmantle pressure noninvasively with a novel nondimensional fluid mechanics model in volunteers and to identify differences related to age and aqueductal dimensions.

MATERIALS AND METHODS

Brain MR images including cardiac-gated 2D phase-contrast MR imaging and fast-spoiled gradient recalled imaging were obtained in 77 volunteers ranging in age from 25-92 years of age. Transmantle pressure was computed during the cardiac cycle with a fluid mechanics model from the measured aqueductal flow rate, stroke volume, aqueductal length and cross-sectional area, and heart rate. Peak pressures during caudal and rostral aqueductal flow were tabulated. The computed transmantle pressure, aqueductal dimensions, and stroke volume were estimated, and the differences due to sex and age were calculated and tested for significance.

RESULTS

Peak transmantle pressure was calculated with the nondimensional averaged 14.4 (SD, 6.5) Pa during caudal flow and 6.9 (SD, 2.8) Pa during rostral flow. It did not differ significantly between men and women or correlate significantly with heart rate. Peak transmantle pressure increased with age and correlated with aqueductal dimensions and stroke volume.

CONCLUSIONS

The nondimensional fluid mechanics model for computing transmantle pressure detected changes in pressure related to age and aqueductal dimensions. This novel methodology can be easily used to investigate the clinical relevance of the transmantle pressure in normal pressure hydrocephalus, pediatric communicating hydrocephalus, and other CSF disorders.

Using deep learning convolutional neural networks to automatically perform cerebral aqueduct CSF flow analysis

Since the development of phase-contrast magnetic resonance imaging (PC-MRI), quantification of cerebrospinal fluid (CSF) flow across the cerebral aqueduct has been utilized for diagnosis of conditions such as normal pressure hydrocephalus (NPH). This study aims to develop an automated method of aqueduct CSF flow analysis using convolution neural networks (CNNs), which can replace the current standard involving manual segmentation of aqueduct region of interest (ROI). Retrospective analysis was performed on 333 patients who underwent PC-MRI, totaling 353 imaging studies. Aqueduct flow measurements using manual ROI placement was performed independently by two radiologists. Two types of CNNs, MultiResUNet and UNet, were trained using ROI data from the senior radiologist, with PC-MRI studies being randomly divided into training (80%) and validation (20%) datasets. Segmentation performance was assessed using Dice similarity coefficient (DSC), and CSF flow parameters were calculated from both manual and CNN-derived ROIs. MultiResUNet, UNet and second radiologist (Rater 2) had DSCs of 0.933, 0.928, and 0.867, respectively, with p < 0.001 between CNNs and Rater 2. Comparison of CSF flow parameters showed excellent intraclass correlation coefficients (ICCs) for MultiResUNet, with lowest correlation being 0.67. For UNet, lower ICCs of -0.01 to 0.56 were observed. Only 3/353 (0.8%) studies failed to have appropriate ROIs placed by MultiResUNet, compared to 12/353 (3.4%) failed cases for UNet. In conclusion, CNNs were able to measure aqueductal CSF flow with similar performance to a senior neuroradiologist. MultiResUNet demonstrated fewer cases of segmentation failure and more consistent flow measurements compared to the widely adopted UNet.

Cerebrospinal fluid flow in normal beagle dogs analyzed using magnetic resonance imaging

Background

Diseases related to cerebrospinal fluid flow, such as hydrocephalus, syringomyelia, and Chiari malformation, are often found in small dogs. Although studies in human medicine have revealed a correlation with cerebrospinal fluid flow in these diseases by magnetic resonance imaging, there is little information and no standard data for normal dogs.

Objectives

The purpose of this study was to obtain cerebrospinal fluid flow velocity data from the cerebral aqueduct and subarachnoid space at the foramen magnum in healthy beagle dogs.

Methods

Six healthy beagle dogs were used in this experimental study. The dogs underwent phase-contrast and time-spatial labeling inversion pulse magnetic resonance imaging. Flow rate variations in the cerebrospinal fluid were observed using sagittal time-spatial labeling inversion pulse images. The pattern and velocity of cerebrospinal fluid flow were assessed using phase-contrast magnetic resonance imaging within the subarachnoid space at the foramen magnum level and the cerebral aqueduct.

Results

In the ventral aspect of the subarachnoid space and cerebral aqueduct, the cerebrospinal fluid was characterized by a bidirectional flow throughout the cardiac cycle. The mean ± SD peak velocities through the ventral and dorsal aspects of the subarachnoid space and the cerebral aqueduct were 1.39 ± 0.13, 0.32 ± 0.12, and 0.76 ± 0.43 cm/s, respectively.

Conclusions

Noninvasive visualization of cerebrospinal fluid flow movement with magnetic resonance imaging was feasible, and a reference dataset of cerebrospinal fluid flow peak velocities was obtained through the cervical subarachnoid space and cerebral aqueduct in healthy dogs.

Diffusion and Flow MR Imaging to Investigate Hydrocephalus Patients Before and After Endoscopic Third Ventriculostomy

INTRODUCTION

In patients with noncommunicating hydrocephalus, dilation of the ventricles stresses white matter fibers and alters the cerebral blood flow (CBF) and cerebrospinal fluid (CSF) dynamics. The purpose of this work was to investigate, non-invasively, how endoscopic third ventriculostomy (ETV) impacts white matter, CSF oscillations, and CBF.

METHODS

Eleven patients presenting with chronic headaches and noncommunicating hydrocephalus due to aqueductal stenosis were treated by ETV. Phase Contrast-MRI (PCMRI) and Diffusion Tensor Imaging (DTI) were performed before and after surgery to evaluate CSF and CBF as well as white matter stresses in the Corpus Callosum (CC) and Corona Radiata (CR). ETV success was confirmed by quantification of the CSF oscillations through the aperture in the third ventricle.

RESULTS

All patients improved after surgery. CSF stroke volume was five times greater than normal ventricular stroke volume. Decrease in cervical CSF oscillations and increase in CBF were observed after ETV. In CR, fiber anisotropy decreased, while water diffusion increased. In CC, anisotropy did not vary, while water diffusion also increased.

CONCLUSION

Even if static ICP typically do not increase, CSF and blood flow are impacted. PCMRI and DTI can provide useful information to help neurosurgeons select patients with good chance to improve after ETV.

Exploring mechanisms of ventricular enlargement in idiopathic normal pressure hydrocephalus: a role of cerebrospinal fluid dynamics and motile cilia

Idiopathic normal pressure hydrocephalus (iNPH) is considered an age-dependent chronic communicating hydrocephalus associated with cerebrospinal fluid (CSF) malabsorption; however, the aetiology of ventricular enlargement in iNPH has not yet been elucidated. There is accumulating evidence that support the hypothesis that various alterations in CSF dynamics contribute to ventricle dilatation in iNPH. This review focuses on CSF dynamics associated with ventriculomegaly and summarises the current literature based on three potential aetiology factors: genetic, environmental and hydrodynamic. The majority of gene mutations that cause communicating hydrocephalus were associated with an abnormal structure or dysfunction of motile cilia on the ventricular ependymal cells. Aging, alcohol consumption, sleep apnoea, diabetes and hypertension are candidates for the risk of developing iNPH, although there is no prospective cohort study to investigate the risk factors for iNPH. Alcohol intake may be associated with the dysfunction of ependymal cilia and sustained high CSF sugar concentration due to uncontrolled diabetes increases the fluid viscosity which in turn increases the shear stress on the ventricular wall surface. Sleep apnoea, diabetes and hypertension are known to be associated with the impairment of CSF and interstitial fluid exchange. Oscillatory shear stress to the ventricle wall surfaces is considerably increased by reciprocating bidirectional CSF movements in iNPH. Increased oscillatory shear stress impedes normal cilia beating, leading to motile cilia shedding from the ependymal cells. At the lack of ciliary protection, the ventricular wall is directly exposed to increased oscillatory shear stress. Additionally, increased oscillatory shear stress may be involved in activating the flow-mediated dilation signalling of the ventricular wall. In conclusion, as the CSF stroke volume at the cerebral aqueduct increases, the oscillatory shear stress increases, promoting motor cilia shedding and loss of ependymal cell coverage. These are considered to be the leading causes of ventricular enlargement in iNPH.

Strain Tensor Imaging: Cardiac-induced brain tissue deformation in humans quantified with high-field MRI

The cardiac cycle induces blood volume pulsations in the cerebral microvasculature that cause subtle deformation of the surrounding tissue. These tissue deformations are highly relevant as a potential source of information on the brain's microvasculature as well as of tissue condition. Besides, cyclic brain tissue deformations may be a driving force in clearance of brain waste products. We have developed a high-field magnetic resonance imaging (MRI) technique to capture these tissue deformations with full brain coverage and sufficient signal-to-noise to derive the cardiac-induced strain tensor on a voxel by voxel basis, that could not be assessed non-invasively before. We acquired the strain tensor with 3 mm isotropic resolution in 9 subjects with repeated measurements for 8 subjects. The strain tensor yielded both positive and negative eigenvalues (principle strains), reflecting the Poison effect in tissue. The principle strain associated with expansion followed the known funnel shaped brain motion pattern pointing towards the foramen magnum. Furthermore, we evaluate two scalar quantities from the strain tensor: the volumetric strain and octahedral shear strain. These quantities showed consistent patterns between subjects, and yielded repeatable results: the peak systolic volumetric strain (relative to end-diastolic strain) was 4.19⋅10^-4 ± 0.78⋅10^-4 and 3.98⋅10^-4 ± 0.44⋅10^-4 (mean ± standard deviation for first and second measurement, respectively), and the peak octahedral shear strain was 2.16⋅10^-3 ± 0.31⋅10^-3 and 2.31⋅10^-3 ± 0.38⋅10^-3, for the first and second measurement, respectively. The volumetric strain was typically highest in the cortex and lowest in the periventricular white matter, while anisotropy was highest in the subcortical white matter and basal ganglia. This technique thus reveals new, regional information on the brain's cardiac-induced deformation characteristics, and has the potential to advance our understanding of the role of microvascular pulsations in health and disease.

Direction and magnitude of cerebrospinal fluid flow vary substantially across central nervous system diseases

Background

Several central nervous system diseases are associated with disturbed cerebrospinal fluid (CSF) flow patterns and have typically been characterized in vivo by phase-contrast magnetic resonance imaging (MRI). This technique is, however, limited by its applicability in space and time. Phase-contrast MRI has yet to be compared directly with CSF tracer enhanced imaging, which can be considered gold standard for assessing long-term CSF flow dynamics within the intracranial compartment.

Methods

Here, we studied patients with various CSF disorders and compared MRI biomarkers of CSF space anatomy and phase-contrast MRI at level of the aqueduct and cranio-cervical junction with dynamic intrathecal contrast-enhanced MRI using the contrast agent gadobutrol as CSF tracer. Tracer enrichment of cerebral ventricles was graded 0–4 by visual assessment. An intracranial pressure (ICP) score was used as surrogate marker of intracranial compliance.

Results

The study included 94 patients and disclosed marked variation of CSF flow measures across disease categories. The grade of supra-aqueductal reflux of tracer varied, with strong reflux (grades 3–4) in half of patients. Ventricular tracer reflux correlated with stroke volume and aqueductal CSF pressure gradient. CSF flow in the cerebral aqueduct was retrograde (from 4th to 3rd ventricle) in one third of patients, with estimated CSF net flow volume about 1.0 L/24 h. In the cranio-cervical junction, net flow was cranially directed in 78% patients, with estimated CSF net flow volume about 4.7 L/24 h.

Conclusions

The present observations provide in vivo quantitative evidence for substantial variation in direction and magnitude of CSF flow, with re-direction of aqueductal flow in communicating hydrocephalus, and significant extra-cranial CSF production. The grading of ventricular reflux of tracer shows promise as a clinical useful method to assess CSF flow pattern disturbances in patients.

Graphic abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s12987-021-00251-6.

Variations in the cerebrospinal fluid dynamics of the American alligator (Alligator mississippiensis)

BACKGROUND

Studies of mammalian CSF dynamics have been focused on three things: paravascular flow, pressure and pulsatility, and "bulk" flow; and three (respective) potential motive forces have been identified: vasomotor, cardiac, and ventilatory. There are unresolved questions in each area, and few links between the different areas. The American alligator (Alligator mississippiensis) has pronounced plasticity in its ventilatory and cardiovascular systems. This study was designed to test the hypothesis that the greater cardiovascular and ventilatory plasticity of A. mississippiensis would result in more variation within the CSF dynamics of this species.

METHODS

Pressure transducers were surgically implanted into the cranial subarachnoid space of 12 sub-adult alligators; CSF pressure and pulsatility were monitored along with EKG and the exhalatory gases. In four of the alligators a second pressure transducer was implanted into the spinal subarachnoid space. In five of the alligators the CSF was labeled with artificial microspheres and Doppler ultrasonography used to quantify aspects of the spinal CSF flow.

RESULTS

Both temporal and frequency analyses of the CSF pulsations showed highly variable contributions of both the cardiac and ventilatory cycles. Unlike the mammalian condition, the CSF pressure pulsations in the alligator are often of long (~ 3 s) duration, and similar duration CSF unidirectional flow pulses were recorded along the spinal cord. Reduction of the duration of the CSF pulsations, as during tachycardia, can lead to a "summation" of the pulsations. There appears to be a minimum duration (~ 1 s) of isolated CSF pulsations. Simultaneous recordings of cranial and spinal CSF pressures reveal a 200 ms delay in the propagation of the pressure pulse from the cranium to the vertebral canal.

CONCLUSIONS

Most of the CSF flow dynamics recorded from the alligators, are similar to what has been reported from studies of the human CSF. It is hypothesized that the link between ventilatory mechanics and CSF pulsations in the alligator is mediated by displacement of the spinal dura. The results of the study suggest that understanding the CSF dynamics of Alligator may provide unique insights into the evolutionary origins and functional regulation of the human CSF dynamics.

Amplified Flow Imaging (aFlow): A Novel MRI-Based Tool to Unravel the Coupled Dynamics Between the Human Brain and Cerebrovasculature

With each heartbeat, periodic variations in arterial blood pressure are transmitted along the vasculature, resulting in localized deformations of the arterial wall and its surrounding tissue. Quantification of such motions may help understand various cerebrovascular conditions, yet it has proven technically challenging thus far. We introduce a new image processing algorithm called amplified Flow (aFlow) which allows to study the coupled brain-blood flow motion by combining the amplification of cine and 4D flow MRI. By incorporating a modal analysis technique known as dynamic mode decomposition into the algorithm, aFlow is able to capture the characteristics of transient events present in the brain and arterial wall deformation. Validating aFlow, we tested it on phantom simulations mimicking arterial walls motion and observed that aFlow displays almost twice higher SNR than its predecessor amplified MRI (aMRI). We then applied aFlow to 4D flow and cine MRI datasets of 5 healthy subjects, finding high correlations between blood flow velocity and tissue deformation in selected brain regions, with correlation values r = 0.61 , 0.59, 0.52 for the pons, frontal and occipital lobe ( ). Finally, we explored the potential diagnostic applicability of aFlow by studying intracranial aneurysm dynamics, which seems to be indicative of rupture risk. In two patients, aFlow successfully visualized the imperceptible aneurysm wall motion, additionally quantifying the increase in the high frequency wall displacement after a one-year follow-up period (20%, 76%). These preliminary data suggest that aFlow may provide a novel imaging biomarker for the assessment of aneurysms evolution, with important potential diagnostic implications.

A novel non-invasive method for estimating the local wave speed at a single site in the internal carotid artery

Objectives Local wave speed is a biomarker which provides an objective analysis of the cardiovascular function. The aim of this study was to determine the local wave speed in the internal carotid artery by a new non-invasive method that measures blood velocity waveform at only one site. Methods For this purpose, the cepstral analysis was employed to determine the arrival time of the reflection wave and the wave speed in the carotid artery. To validate our model, we applied it experimentally in vivo on young and old healthy subjects. The blood velocity waveform was measured by using phase-contrast magnetic resonance for 22 subjects. Results Our experimental results correlated with reference values reported in previous studies conducted on the internal arterial carotid usually adopting the invasive method. They also correlated with those obtained by using the foot-to-foot method (R2=0.72). The wave speed obtained by the method developed in this study and that of the foot-to-foot method increased with age (p<0.001). Conclusions The method developed in this study can be applied in the other arteries and it can also be used with other techniques such as ultrasound imaging.

Cerebrospinal fluid dynamics in idiopathic normal pressure hydrocephalus on four-dimensional flow imaging

OBJECTIVES

To evaluate complex CSF movements and shear stress in patients with idiopathic normal pressure hydrocephalus (iNPH) on four-dimensional (4D) flow MRI.

METHODS

Three-dimensional velocities and volumes of the reciprocating CSF movements through 12 ROIs from the foramen of Monro to the upper cervical spine were measured in 41 patients with iNPH, 23 patients with co-occurrence of iNPH and Alzheimer's disease (AD), and 9 age-matched controls, using 4D flow imaging and application. Stroke volume, reversed-flow rate, and shear stress were automatically calculated. Relationships between flow-related parameters and morphological measurements were also assessed.

RESULTS

Stroke volumes, reversed-flow rates, and shear stress at the cerebral aqueduct were significantly higher in patients with iNPH than in controls. Patients with pure iNPH had significantly higher shear stress at the ventral aspect of the cerebral aqueduct than those with co-occurrence of iNPH and AD. The stroke volume at the upper end of the cerebral aqueduct had the strongest association with the anteroposterior diameter of the lower end of the cerebral aqueduct (r = 0.52). The stroke volume at the foramen of Monro had significant associations with the indices specific to iNPH. The shear stress at the dorsal aspect of the cerebral aqueduct had the strongest association with the diameter of the foramen of Magendie (r = 0.52).

CONCLUSIONS

Stroke volumes, reversed-flow rates, and shear stress through the cerebral aqueduct on 4D flow MRI are useful parameters for iNPH diagnosis. These findings can aid in elucidating the mechanism of ventricular enlargement in iNPH.

KEY POINTS

• The CSF stroke volume and bimodal shear stress at the cerebral aqueduct were considerably higher in patients with iNPH. • The patients with pure iNPH had significantly higher shear stress at the ventral aspect of the cerebral aqueduct than those with co-occurrence of iNPH and AD. • The shear stress at the cerebral aqueduct was significantly associated with the diameter of the foramen of Magendie.

In vitro and numerical simulation of blood removal from cerebrospinal fluid: comparison of lumbar drain to Neurapheresis therapy

BACKGROUND

Blood removal from cerebrospinal fluid (CSF) in post-subarachnoid hemorrhage patients may reduce the risk of related secondary brain injury. We formulated a computational fluid dynamics (CFD) model to investigate the impact of a dual-lumen catheter-based CSF filtration system, called Neurapheresis™ therapy, on blood removal from CSF compared to lumbar drain.

METHODS

A subject-specific multiphase CFD model of CSF system-wide solute transport was constructed based on MRI measurements. The Neurapheresis catheter geometry was added to the model within the spinal subarachnoid space (SAS). Neurapheresis flow aspiration and return rate was 2.0 and 1.8 mL/min, versus 0.2 mL/min drainage for lumbar drain. Blood was modeled as a bulk fluid phase within CSF with a 10% initial tracer concentration and identical viscosity and density as CSF. Subject-specific oscillatory CSF flow was applied at the model inlet. The dura and spinal cord geometry were considered to be stationary. Spatial-temporal tracer concentration was quantified based on time-average steady-streaming velocities throughout the domain under Neurapheresis therapy and lumbar drain. To help verify CFD results, an optically clear in vitro CSF model was constructed with fluorescein used as a blood surrogate. Quantitative comparison of numerical and in vitro results was performed by linear regression of spatial-temporal tracer concentration over 24-h.

RESULTS

After 24-h, tracer concentration was reduced to 4.9% under Neurapheresis therapy compared to 6.5% under lumbar drain. Tracer clearance was most rapid between the catheter aspiration and return ports. Neurapheresis therapy was found to have a greater impact on steady-streaming compared to lumbar drain. Steady-streaming in the cranial SAS was ~ 50× smaller than in the spinal SAS for both cases. CFD results were strongly correlated with the in vitro spatial-temporal tracer concentration under Neurapheresis therapy (R² = 0.89 with + 2.13% and - 1.93% tracer concentration confidence interval).

CONCLUSION

A subject-specific CFD model of CSF system-wide solute transport was used to investigate the impact of Neurapheresis therapy on tracer removal from CSF compared to lumbar drain over a 24-h period. Neurapheresis therapy was found to substantially increase tracer clearance compared to lumbar drain. The multiphase CFD results were verified by in vitro fluorescein tracer experiments.

Can pulsatile CSF flow across the cerebral aqueduct cause ventriculomegaly? A prospective study of patients with communicating hydrocephalus

Background

Communicating hydrocephalus is a disease where the cerebral ventricles are enlarged. It is characterized by the absence of detectable cerebrospinal fluid (CSF) outflow obstructions and often with increased CSF pulsatility measured in the cerebral aqueduct (CA). We hypothesize that the cardiac-related pulsatile flow over the CA, with fast systolic outflow and slow diastolic inflow, can generate net pressure effects that could source the ventriculomegaly in these patients. This would require a non-zero cardiac cycle averaged net pressure difference (ΔP_net) over the CA, with higher average pressure in the lateral and third ventricles.

Methods

We tested the hypothesis by calculating ΔP_net across the CA using computational fluid dynamics based on prospectively collected high-resolution structural (FIESTA-C, resolution 0.39 × 0.39 × 0.3 mm³) and velocimetric (2D-PCMRI, in-plane resolution 0.35 × 0.35 mm²) MRI-data from 30 patients investigated for communicating hydrocephalus.

Results

The ΔP_net due to CSF pulsations was non-zero for the study group (p = 0.03) with a magnitude of 0.2 ± 0.4 Pa (0.001 ± 0.003 mmHg), with higher pressure in the third ventricle. The maximum pressure difference over the cardiac cycle ΔP_max was 20.3 ± 11.8 Pa and occurred during systole. A generalized linear model verified an association between ΔP_net and CA cross-sectional area (p = 0.01) and flow asymmetry, described by the ratio of maximum inflow/outflow (p = 0.04), but not for aqueductal stroke volume (p = 0.35).

Conclusions

The results supported the hypothesis with respect to the direction of ΔP_net, although the magnitude was low. Thus, although the pulsations may generate a pressure difference across the CA it is likely too small to explain the ventriculomegaly in communicating hydrocephalus.

Uncertainty quantification of parenchymal tracer distribution using random diffusion and convective velocity fields

Background

Influx and clearance of substances in the brain parenchyma occur by a combination of diffusion and convection, but the relative importance of these mechanisms is unclear. Accurate modeling of tracer distributions in the brain relies on parameters that are partially unknown and with literature values varying by several orders of magnitude. In this work, we rigorously quantified the variability of tracer distribution in the brain resulting from uncertainty in diffusion and convection model parameters.

Methods

Using the convection–diffusion–reaction equation, we simulated tracer distribution in the brain parenchyma after intrathecal injection. Several models were tested to assess the uncertainty both in type of diffusion and velocity fields and also the importance of their magnitude. Our results were compared with experimental MRI results of tracer enhancement.

Results

In models of pure diffusion, the expected amount of tracer in the gray matter reached peak value after 15 h, while the white matter did not reach peak within 24 h with high likelihood. Models of the glymphatic system were similar qualitatively to the models of pure diffusion with respect to expected time to peak but displayed less variability. However, the expected time to peak was reduced to 11 h when an additional directionality was prescribed for the glymphatic circulation. In a model including drainage directly from the brain parenchyma, time to peak occured after 6–8 h for the gray matter.

Conclusion

Even when uncertainties are taken into account, we find that diffusion alone is not sufficient to explain transport of tracer deep into the white matter as seen in experimental data. A glymphatic velocity field may increase transport if a large-scale directional structure is included in the glymphatic circulation.

Sex and Age Dependencies of Aqueductal Cerebrospinal Fluid Dynamics Parameters in Healthy Subjects

Objectives

To assess the influence of age and sex on 10 cerebrospinal fluid (CSF) flow dynamics parameters measured with an MR phase contrast (PC) sequence within the cerebral aqueduct at the level of the intercollicular sulcus.

Materials and Methods

128 healthy subjects (66 female subjects with a mean age of 52.9 years and 62 male subjects with a mean age of 51.8 years) with a normal Evans index, normal medial temporal atrophy (MTA) score, and without known disorders of the CSF circulation were included in the study. A PC MR sequence on a 3T MR scanner was used. Ten different flow parameters were analyzed using postprocessing software. Ordinal and linear regression models were calculated.

Results

The parameters stroke volume (sex: p < 0.001, age: p = 0.003), forward flow volume (sex: p < 0.001, age: p = 0.002), backward flow volume (sex: p < 0.001, age: p = 0.018), absolute stroke volume (sex: p < 0.001, age: p = 0.005), mean flux (sex: p < 0.001, age: p = 0.001), peak velocity (sex: p = 0.009, age: p = 0.0016), and peak pressure gradient (sex: p = 0.029, age: p = 0.028) are significantly influenced by sex and age. The parameters regurgitant fraction, stroke distance, and mean velocity are not significantly influenced by sex and age.

Conclusion

CSF flow dynamics parameters measured in the cerebral aqueduct are partly age and sex dependent. For establishment of reliable reference values for clinical use in future studies, the impact of sex and age should be considered and incorporated.

Predicting the Aqueductal Cerebrospinal Fluid Pulse: A Statistical Approach

The cerebrospinal fluid (CSF) pulse in the Aqueduct of Sylvius (aCSF pulse) is often used to evaluate structural changes in the brain. Here we present a novel application of the general linear model (GLM) to predict the motion of the aCSF pulse. MR venography was performed on 13 healthy adults (9 female and 4 males—mean age = 33.2 years). Flow data was acquired from the arterial, venous and CSF vessels in the neck (C2/C3 level) and from the AoS. Regression analysis was undertaken to predict the motion of the aCSF pulse using the cervical flow rates as predictor variables. The relative contribution of these variables to predicting aCSF flow rate was assessed using a relative weights method, coupled with an ANOVA. Analysis revealed that the aCSF pulse could be accurately predicted (mean (SD) adjusted r2 = 0.794 (0.184)) using the GLM (p < 0.01). Venous flow rate in the neck was the strongest predictor of aCSF pulse (p = 0.001). In healthy individuals, the motion of the aCSF pulse can be predicted using the GLM. This indicates that the intracranial fluidic system has broadly linear characteristics. Venous flow in the neck is the strongest predictor of the aCSF pulse.

Decreased Cerebrospinal Fluid Flow Is Associated With Cognitive Deficit in Elderly Patients

Background: Disruptions in cerebrospinal fluid (CSF) flow during aging could compromise protein clearance from the brain and contribute to the etiology of Alzheimer’s Disease (AD).

Objective: To determine whether CSF flow is associated with cognitive deficit in elderly patients (>70 years).

Methods: We studied 92 patients admitted to our geriatric unit for non-acute reasons using phase-contrast magnetic resonance imaging (PC-MRI) to calculate their ventricular and spinal CSF flow, and assessed their global cognitive status, memory, executive functions, and praxis. Multivariable regressions with backward selection (criterion p < 0.15) were performed to determine associations between cognitive tests and ventricular and spinal CSF flow, adjusting for depression, anxiety, and cardiovascular risk factors.

Results: The cohort comprised 71 women (77%) and 21 (33%) men, aged 84.1 ± 5.2 years (range, 73–96). Net ventricular CSF flow was 52 ± 40 μL/cc (range, 0–210), and net spinal CSF flow was 500 ± 295 μL/cc (range, 0–1420). Ventricular CSF flow was associated with the number of BEC96 figures recognized (β = 0.18, CI, 0.02–0.33; p = 0.025). Spinal CSF flow was associated with the WAIS Digit Span Backward test (β = 0.06, CI, 0.01–0.12; p = 0.034), and categoric verbal fluency (β = 0.53, CI, 0.07–0.98; p = 0.024) and semantic verbal fluency (β = 0.55, CI, 0.07–1.02; p = 0.024).

Conclusion: Patients with lower CSF flow had significantly worse memory, visuo-constructive capacities, and verbal fluency. Alterations in CSF flow could contribute to some of the cognitive deficit observed in patients with AD. Diagnosis and treatment of CSF flow alterations in geriatric patients with neurocognitive disorders could contribute to the prevention of their cognitive decline.

Enhanced in vitro model of the CSF dynamics

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Assessing test-retest reliability of phase contrast MRI for measuring cerebrospinal fluid and cerebral blood flow dynamics

PURPOSE

Pathological states occur when cerebrospinal fluid (CSF) and cerebral blood flow (CBF) dynamics become dysregulated in the brain. Phase-contrast MRI (PC-MRI) is a noninvasive imaging technique that enables quantitative measurements of CSF and CBF flow. While studies have validated PC-MRI as an imaging technique for flow, few studies have evaluated its reliability for CSF and CBF flow parameters commonly associated with neurological disease. The purpose of this study was to evaluate test-retest reliability at the cerebral aqueduct (CA) and C2-C3 area using PC-MRI to assess the feasibility of investigating CSF and CBF flow dynamics.

METHODS

This study was performed on 27 cognitively normal young adults (ages 20-35 years). Flow data was acquired on a 3T Siemens Prisma using a 2D cine-PC pulse sequence. Three consecutive flow measurements were acquired at the CA and C2-C3 area. Intraclass correlation coefficient (ICC) and coefficient of variance (CV) were used to evaluate intrarater, inter-rater, and test-retest reliability.

RESULTS

Among the 26 flow parameters analyzed, 22 had excellent reliability (ICC > 0.80), including measurements of CSF stroke volume, flush peak, and fill peak, and 4 parameters had good reliability (ICC 0.60-0.79). 16 flow parameters had a mean CV ≤ 10%, 7 had a CV ≤ 15%, and 3 had a CV ≤ 30%. All CSF and CBF flow measurements had excellent inter-rater and intrarater reliability (ICC > 0.80).

CONCLUSION

This study shows that CSF and CBF flow can be reliably measured at the CA and C2-C3 area using PC-MRI, making it a promising tool for studying flow dynamics in the central nervous system.

Recent Advances in Rational Diagnosis and Treatment of Normal Pressure Hydrocephalus: A Critical Appraisal on Novel Diagnostic, Therapy Monitoring and Treatment Modalities

BACKGROUND

Normal pressure hydrocephalus (NPH) is a critical brain disorder in which excess Cerebrospinal Fluid (CSF) is accumulated in the brain's ventricles causing damage or disruption of the brain tissues. Amongst various signs and symptoms, difficulty in walking, slurred speech, impaired decision making and critical thinking, and loss of bladder and bowl control are considered the hallmark features of NPH.

OBJECTIVE

The current review was aimed to present a comprehensive overview and critical appraisal of majorly employed neuroimaging techniques for rational diagnosis and effective monitoring of the effectiveness of the employed therapeutic intervention for NPH. Moreover, a critical overview of recent developments and utilization of pharmacological agents for the treatment of hydrocephalus has also been appraised.

RESULTS

Considering the complications associated with the shunt-based surgical operations, consistent monitoring of shunting via neuroimaging techniques hold greater clinical significance. Despite having extensive applicability of MRI and CT scan, these conventional neuroimaging techniques are associated with misdiagnosis or several health risks to patients. Recent advances in MRI (i.e., Sagittal-MRI, coronal-MRI, Time-SLIP (time-spatial-labeling-inversion-pulse), PC-MRI and diffusion-tensor-imaging (DTI)) have shown promising applicability in the diagnosis of NPH. Having associated with several adverse effects with surgical interventions, non-invasive approaches (pharmacological agents) have earned greater interest of scientists, medical professional, and healthcare providers. Amongst pharmacological agents, diuretics, isosorbide, osmotic agents, carbonic anhydrase inhibitors, glucocorticoids, NSAIDs, digoxin, and gold-198 have been employed for the management of NPH and prevention of secondary sensory/intellectual complications.

CONCLUSION

Employment of rational diagnostic tool and therapeutic modalities avoids misleading diagnosis and sophisticated management of hydrocephalus by efficient reduction of Cerebrospinal Fluid (CSF) production, reduction of fibrotic and inflammatory cascades secondary to meningitis and hemorrhage, and protection of brain from further deterioration.

Measurement of microvascular cerebral blood volume changes over the cardiac cycle with ferumoxytol-enhanced T2 * MRI

PURPOSE

This feasibility study investigates the non-invasive measurement of microvascular cerebral blood volume (BV) changes over the cardiac cycle using cardiac-gated, ferumoxytol-enhanced T 2 ∗ MRI.

METHODS

Institutional review board approval was obtained and all subjects provided written informed consent. Cardiac gated MR scans were prospectively acquired on a 3.0T scanner in 22 healthy subjects using T 2 ∗ -weighted sequences with 2D-EPI and 3D spiral trajectories. Images were collected before and after the intravenous administration of 2 doses of ferumoxytol (1 mg FE/kg and 4 mg FE/kg). Cardiac cycle-induced R 2 ∗ (1/ T 2 ∗ ) changes (Δ R 2 ∗ ) and BV changes (ΔBV) throughout the cardiac cycle in gray matter (GM) and white matter (WM) were quantified and differences assessed using ANOVA followed by post hoc analysis.

RESULTS

Δ R 2 ∗ was found to increase in a dose-dependent fashion. A significantly larger increase was observed in GM compared to WM in both 2D and 3D acquisitions (P < 0.050). In addition, Δ R 2 ∗ increased significantly (P < 0.001) post versus pre-contrast injection in GM in both T 2 ∗ MRI acquisitions. Mean GM Δ R 2 ∗ derived from 2D-EPI images was 0.14 ± 0.06 s^-1 pre-contrast and 0.33 ± 0.13 s^-1 after 5 mg FE/kg. In WM, Δ R 2 ∗ was 0.19 ± 0.06 s^-1 pre-contrast, and 0.23 ± 0.06 s^-1 after 5 mg FE/kg. The fractional changes in BV throughout the cardiac cycle were 0.031 ± 0.019% in GM and 0.011 ± 0.008% in WM (P < 0.001) after 5 mg FE/kg.

CONCLUSION

Cardiac-gated, ferumoxytol-enhanced T 2 ∗ MRI enables characterization of microvascular BV changes throughout the cardiac cycle in GM and WM tissue of healthy subjects.

Opposing CSF hydrodynamic trends found in the cerebral aqueduct and prepontine cistern following shunt treatment in patients with normal pressure hydrocephalus

Background

This study investigated cerebrospinal fluid (CSF) hydrodynamics using cine phase-contrast MRI in the cerebral aqueduct and the prepontine cistern between three distinct groups: pre-shunt normal pressure hydrocephalus (NPH) patients, post-shunt NPH patients, and controls. We hypothesized that the hyperdynamic flow of CSF through the cerebral aqueduct seen in NPH patients was due to a reduction in cisternal CSF volume buffering. Both hydrodynamic (velocity, flow, stroke volume) and peak flow latency (PFL) parameters were investigated.

Methods

Scans were conducted on 30 pre-treatment patients ranging in age from 58 to 88 years along with an additional 12 controls. Twelve patients also received scans following either ventriculoatrial (VA) or ventriculoperitoneal (VP) shunt treatment (9 VP, 3 VA), ranging in age from 74 to 89 years with a mean follow up time of 6 months.

Results

Significant differences in area, velocity, flow, and stroke volume for the cerebral aqueduct were found between the pre-treatment NPH group and the healthy controls. Shunting caused a significant decrease in both caudal and cranial mean flow and stroke volume in the cerebral aqueduct. No significant changes were found in the prepontine cistern between the pre-treatment group and healthy controls. For the PFL, no significant differences were seen in the cerebral aqueduct between any of the three groups; however, the prepontine cistern PFL was significantly decreased in the pre-treatment NPH group when compared to the control group.

Conclusions

Although several studies have quantified the changes in aqueductal flow between hydrocephalic groups and controls, few studies have investigated prepontine cistern flow. Our study was the first to investigate both regions in the same patients for NPH pre- and post- treatment. Following shunt treatment, the aqueductal CSF metrics decreased toward control values, while the prepontine cistern metrics trended up (not significantly) from the normal values established in this study. The opposing trend of the two locations suggests a redistribution of CSF pulsatility in NPH patients. Furthermore, the significantly decreased latency of the prepontine cisternal CSF flow suggests additional evidence for CSF pulsatility dysfunction.

DTI Profiles for Rapid Description of Cohorts at the Clinical-Research Interface

Normal pressure hydrocephalus (NPH) is a syndrome comprising gait disturbance, cognitive decline and urinary incontinence that is an unique model of reversible brain injury, but it presents as a challenging spectrum of disease cohorts. Diffusion Tensor Imaging (DTI), with its ability to interrogate structural white matter patterns at a microarchitectural level, is a potentially useful tool for the confirmation and characterization of disease cohorts at the clinical-research interface. However, obstacles to its widespread use involve the need for consistent DTI analysis and interpretation tools across collaborator sites. We present the use of DTI profiles, a simplistic methodology to interpret white matter injury patterns based on the morphology of diffusivity parameters. We examined 13 patients with complex NPH, i.e., patients with NPH and overlay from multiple comorbidities, including vascular risk burden and neurodegenerative disease, undergoing extended CSF drainage, clinical assessments, and multi-modal MR imaging. Following appropriate exclusions, we compared the morphology of DTI profiles in such complex NPH patients (n = 12, comprising 4 responders and 8 non-responders) to exemplar DTI profiles from a cohort of classic NPH patients (n = 16) demonstrating responsiveness of white matter injury to ventriculo-peritoneal shunting. In the cohort of complex NPH patients, mean age was 71.3 ± 7.6 years (10 males, 2 females) with a mean MMSE score of 21.1. There were 5 age-matched healthy controls, mean age was 73.4 ± 7.2 years (1 male, 4 females) and mean MMSE score was 26.8. In the exemplar cohort of classic NPH patients, mean age was 74.7 ± 5.9 years (10 males, 6 females) and mean MMSE score was 24.1. There were 9 age-matched healthy controls, mean age was 69.4 ± 9.7 years (4 males, 5 females) and mean MMSE score was 28.6. We found that, despite the challenges of acquiring DTI metrics from differing scanners across collaborator sites and NPH patients presenting as differing cohorts along the spectrum of disease, DTI profiles for responsiveness to interventions were comparable. Distinct DTI characteristics were demonstrated for complex NPH responders vs. non-responders. The morphology of DTI profiles for complex NPH responders mimicked DTI patterns found in predominantly shunt-responsive patients undergoing intervention for classic NPH. However, DTI profiles for complex NPH non-responders was suggestive of atrophy. Our findings suggest that it is possible to use DTI profiles to provide a methodology for rapid description of differing cohorts of disease at the clinical-research interface. By describing DTI measures morphologically, it was possible to consistently compare white matter injury patterns across international collaborator datasets.

Differences in the Calculated Transvenous Pressure Drop between Chronic Hydrocephalus and Idiopathic Intracranial Hypertension

BACKGROUND AND PURPOSE

Chronic hydrocephalus is associated with dilated ventricles despite a normal intracranial pressure. In idiopathic intracranial hypertension, the ventricles are normal despite an elevated intracranial pressure. This apparent paradox has largely remained unexplained. It is suggested that a pressure difference between the superficial and deep venous territories of the brain could account for the variation between the 2 diseases. The purpose of this paper is to investigate the cause of this pressure difference.

MATERIALS AND METHODS

Using MR phase-contrast imaging, we calculated the hydraulic diameters of the sagittal and straight sinuses in 21 patients with hydrocephalus, 20 patients with idiopathic intracranial hypertension, and 20 age-matched controls. The outflow resistance of each sinus was estimated using the Poiseuille equation. The outflow pressure was estimated using the flow data. A smaller subset of the patients with hydrocephalus had these studies repeated after successful shunt insertion.

RESULTS

In hydrocephalus, the sagittal sinuses were 21% smaller than those in controls (P < .001); the straight sinuses were not significantly different. In idiopathic intracranial hypertension, both sinuses were not significantly different from those of controls. The pressure drop from the sagittal sinus to the end of the straight sinus was elevated by 1.2 mm Hg in hydrocephalus (P = .001) but not significantly different from that in controls in idiopathic intracranial hypertension. Shunt insertion dilated the sagittal sinuses in hydrocephalus, leaving them 18% larger than normal and eliminating the transvenous pressure change.

CONCLUSIONS

There is a transvenous pressure difference in hydrocephalus that is absent in idiopathic intracranial hypertension. This difference is eliminated by shunt insertion. The findings may have a bearing on ventricular dilation.

Cerebrospinal fluid volumetric net flow rate and direction in idiopathic normal pressure hydrocephalus

The aim of the present study was to examine cerebrospinal fluid (CSF) volumetric net flow rate and direction at the cranio-cervical junction (CCJ) and cerebral aqueduct in individuals with idiopathic normal pressure hydrocephalus (iNPH) using cardiac-gated phase-contrast magnetic resonance imaging (PC-MRI). An in-depth, pixel-by-pixel analysis of regions of interest from the CCJ and cerebral aqueduct, respectively, was done in 26 iNPH individuals, and in 4 healthy subjects for validation purposes. Results from patients were compared with over-night measurements of static and pulsatile intracranial pressure (ICP). In iNPH, CSF net flow at CCJ was cranially directed in 17/22 as well as in 4/4 healthy subjects. Estimated daily CSF volumetric net flow rate at CCJ was 6.9 ± 9.9 L/24 h in iNPH patients and 4.5 ± 5.0 L/24 h in healthy individuals. Within the cerebral aqueduct, the CSF net flow was antegrade in 7/21 iNPH patients and in 4/4 healthy subjects, while it was retrograde (i.e. towards ventricles) in 14/21 iNPH patients. Estimated daily CSF volumetric net flow rate in cerebral aqueduct was 1.1 ± 2.2 L/24 h in iNPH while 295 ± 53 mL/24 h in healthy individuals. Magnitude of cranially directed CSF net flow in cerebral aqueduct was highest in iNPH individuals with signs of impaired intracranial compliance. The study results indicate CSF flow volumes and direction that are profoundly different from previously assumed. We hypothesize that spinal CSF formation may serve to buffer increased demand for CSF flow through the glymphatic system during sleep and during deep inspiration to compensate for venous outflow.

Phase contrast MRI measurements of net cerebrospinal fluid flow through the cerebral aqueduct are confounded by respiration

Background

Net cerebrospinal fluid (CSF) flow through the cerebral aqueduct may serve as a marker of CSF production in the lateral ventricles, and changes that occur with aging and in disease.

Purpose

To investigate the confounding influence of the respiratory cycle on net CSF flow and stroke volume measurements.

Study Type

Cross‐sectional study.

Subjects

Twelve young, healthy subjects (seven male, age range 19–39 years, average age 28.3 years).

Field Strength/Sequence

Phase contrast MRI (PC‐MRI) measurements were performed at 7T, with and without respiratory gating on expiration and on inspiration. All measurements were repeated.

Assessment

Net CSF flow and stroke volume in the aqueduct, over the cardiac cycle, was determined.

Statistical Tests

Repeatability was determined using the intraclass correlation coefficient (ICC) and linear regression analysis between the repeated measurements. Repeated measures analysis of variance (ANOVA) was performed to compare the measurements during inspiration/expiration/no gating. Linear regression analysis was performed between the net CSF flow difference (inspiration minus expiration) and stroke volume.

Results

Net CSF flow (average ± standard deviation) was 0.64 ± 0.32 mL/min (caudal) during expiration, 0.12 ± 0.49 mL/min (cranial) during inspiration, and 0.31 ± 0.18 mL/min (caudal) without respiratory gating. Respiratory gating did not affect stroke volume measurements (41 ± 18, 42 ± 19, 42 ± 19 μL/cycle for expiration, no respiratory gating and inspiration, respectively). Repeatability was best during inspiration (ICC = 0.88/0.56/–0.31 for gating on inspiration/expiration/no gating). A positive association was found between average stroke volume and net flow difference between inspiration and expiration (R = 0.678/0.605, P = 0.015/0.037 for the first/second repeated measurement).

Data Conclusion

Measured net CSF flow is confounded by respiration effects. Therefore, net CSF flow measurements with PC‐MRI cannot in isolation be directly linked to CSF production.

Level of Evidence: 1

Technical Efficacy: Stage 2

J. Magn. Reson. Imaging 2019;49:433–444.

Numerical Cerebrospinal System Modeling in Fluid-Structure Interaction

OBJECTIVE

Cerebrospinal fluid (CSF) stroke volume in the aqueduct is widely used to evaluate CSF dynamics disorders. In a healthy population, aqueduct stroke volume represents around 10% of the spinal stroke volume while intracranial subarachnoid space stroke volume represents 90%. The amplitude of the CSF oscillations through the different compartments of the cerebrospinal system is a function of the geometry and the compliances of each compartment, but we suspect that it could also be impacted be the cardiac cycle frequency. To study this CSF distribution, we have developed a numerical model of the cerebrospinal system taking into account cerebral ventricles, intracranial subarachnoid spaces, spinal canal and brain tissue in fluid-structure interactions.

MATERIALS AND METHODS

A numerical fluid-structure interaction model is implemented using a finite-element method library to model the cerebrospinal system and its interaction with the brain based on fluid mechanics equations and linear elasticity equations coupled in a monolithic formulation. The model geometry, simplified in a first approach, is designed in accordance with realistic volume ratios of the different compartments: a thin tube is used to mimic the high flow resistance of the aqueduct. CSF velocity and pressure and brain displacements are obtained as simulation results, and CSF flow and stroke volume are calculated from these results.

RESULTS

Simulation results show a significant variability of aqueduct stroke volume and intracranial subarachnoid space stroke volume in the physiological range of cardiac frequencies.

CONCLUSIONS

Fluid-structure interactions are numerous in the cerebrospinal system and difficult to understand in the rigid skull. The presented model highlights significant variations of stroke volumes under cardiac frequency variations only.