Predicting the Aqueductal Cerebrospinal Fluid Pulse: A Statistical Approach

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.

Results of Numerical Modeling of Blood Flow in the Internal Jugular Vein Exhibiting Different Types of Strictures

The clinical relevance of nozzle-like strictures in upper parts of the internal jugular veins remains unclear. This study was aimed at understanding flow disturbances caused by such stenoses. Computational fluid dynamics software, COMSOL Multiphysics, was used. Two-dimensional computational domain involved stenosis at the beginning of modeled veins, and a flexible valve downstream. The material of the venous valve was considered to be hyperelastic. In the vein models with symmetric 2-leaflets valve without upstream stenosis or with minor 30% stenosis, the flow was undisturbed. In the case of major 60% and 75% upstream stenosis, centerline velocity was positioned asymmetrically, and areas of reverse flow and flow separation developed. In the 2-leaflet models with major stenosis, vortices evoking flow asymmetry were present for the entire course of the model, while the valve leaflets were distorted by asymmetric flow. Our computational fluid dynamics modeling suggests that an impaired outflow from the brain through the internal jugular veins is likely to be primarily caused by pathological strictures in their upper parts. In addition, the jugular valve pathology can be exacerbated by strictures located in the upper segments of these veins.

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.

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.

Measuring Aqueduct of Sylvius Cerebrospinal Fluid Flow in Multiple Sclerosis Using Different Software

Aqueduct of Sylvius (AoS) cerebrospinal fluid flow can be quantified using phase-contrast (PC) Magnetic Resonance Imaging. The software used for AoS segmentation might affect the PC-derived measures. We analyzed AoS PC data of 30 people with multiple sclerosis and 19 normal controls using three software packages, and estimated cross-sectional area (CSA), average and highest AoS velocity (Vmean and Vmax), flow rate and volume. Our aims were to assess the repeatability and reproducibility of each PC-derived measure obtained with the various software packages, including in terms of group differentiation. All the variables had good repeatability, except the average Vmean, flow rate and volume obtained with one software package. Substantial to perfect agreement was seen when evaluating the overlap between the AoS segmentations obtained with different software packages. No variable was significantly different between software packages, with the exception of Vmean diastolic peak and CSA. Vmax diastolic peak differentiated groups, regardless of the software package. In conclusion, a clinical study should preliminarily evaluate the repeatability in order to interpret its findings. Vmax seemed to be a repeatable and reproducible measure, since the pixel with its value is usually located in the center of the AoS, and is thus unlikely be affected by ROI size.