Investigation of the pulsatility of cerebrospinal fluid using cardiac-gated Intravoxel Incoherent Motion imaging

Diffusion magnetic resonance imaging of cerebrospinal fluid dynamics: Current techniques and future advancements

Cerebrospinal fluid (CSF) plays a critical role in metabolic waste clearance from the brain, requiring its circulation throughout various brain pathways, including the ventricular system, subarachnoid spaces, para-arterial spaces, interstitial spaces, and para-venous spaces. The complexity of CSF circulation has posed a challenge in obtaining noninvasive measurements of CSF dynamics. The assessment of CSF dynamics throughout its various circulatory pathways is possible using diffusion magnetic resonance imaging (MRI) with optimized sensitivity to incoherent water movement across the brain. This review presents an overview of both established and emerging diffusion MRI techniques designed to measure CSF dynamics and their potential clinical applications. The discussion offers insights into the optimization of diffusion MRI acquisition parameters to enhance the sensitivity and specificity of diffusion metrics on underlying CSF dynamics. Lastly, we emphasize the importance of cautious interpretations of diffusion-based imaging, especially when differentiating between tissue- and fluid-related changes or elucidating structural versus functional alterations.

Exploring neurodegenerative disorders using advanced magnetic resonance imaging of the glymphatic system

The glymphatic system, a macroscopic waste clearance system in the brain, is crucial for maintaining neural health. It facilitates the exchange of cerebrospinal and interstitial fluid, aiding the clearance of soluble proteins and metabolites and distributing essential nutrients and signaling molecules. Emerging evidence suggests a link between glymphatic dysfunction and the pathogenesis of neurodegenerative disorders, including Alzheimer's, Parkinson's, and Huntington's disease. These disorders are characterized by the accumulation and propagation of misfolded or mutant proteins, a process in which the glymphatic system is likely involved. Impaired glymphatic clearance could lead to the buildup of these toxic proteins, contributing to neurodegeneration. Understanding the glymphatic system's role in these disorders could provide insights into their pathophysiology and pave the way for new therapeutic strategies. Pharmacological enhancement of glymphatic clearance could reduce the burden of toxic proteins and slow disease progression. Neuroimaging techniques, particularly MRI-based methods, have emerged as promising tools for studying the glymphatic system in vivo. These techniques allow for the visualization of glymphatic flow, providing insights into its function under healthy and pathological conditions. This narrative review highlights current MRI-based methodologies, such as motion-sensitizing pulsed field gradient (PFG) based methods, as well as dynamic gadolinium-based and glucose-enhanced methodologies currently used in the study of neurodegenerative disorders.

Clinical magnetic resonance imaging evaluation of glymphatic function

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

D* from diffusion MRI reveals a correspondence between ventricular cerebrospinal fluid volume and flow in the ischemic rodent model

Quantitative measurement of cerebrospinal fluid (CSF) flow and volume and longitudinal monitoring of CSF dynamics provide insights into the compensatory characteristics of post-stroke CSF. In this study, we compared the MRI pseudo-diffusion index (D*) of live and sacrificed rat brains to confirm the effect of ventricular CSF flow on diffusion signals. We observed the relationship between the CSF peak velocities and D* through Monte Carlo (MC) simulations to further understand the source of D* contrast. We also determined the dominant CSF flow using D* in three directions. Finally, we investigated the dynamic evolutions of ventricular CSF flow and volume in a stroke rat model (n = 8) from preoperative to up to 45 days after surgery and determined the correlation between ventricular CSF volume and flow. MC simulations showed a strong positive correlation between the CSF peak velocity and D* (r = 0.99). The dominant CSF flow variations in the 3D ventricle could be measured using the maximum D* map. A longitudinal positive correlation between ventricular CSF volume and D* was observed in the lateral (r = 0.74) and ventral-third (r = 0.81) ventricles, respectively. The directional D* measurements provide quantitative CSF volume and flow information, which would provide useful insights into ischemic stroke with diffusion MRI.

Dynamic brain ADC variations over the cardiac cycle and their relation to tissue strain assessed with DENSE at high-field MRI

Purpose

The ADC of brain tissue slightly varies over the cardiac cycle. This variation could reflect physiology, including mixing of the interstitial fluid, relevant for brain waste clearance. However, it is known from cardiac diffusion imaging that tissue deformation by itself affects the magnitude of the MRI signal, leading to artificial ADC variations as well. This study investigates to what extent tissue deformation causes artificial ADC variations in the brain.

Theory and Methods

We implemented a high‐field MRI sequence with stimulated echo acquisition mode that simultaneously measures brain tissue deformation and ADC. Based on the measured tissue deformation, we simulated the artificial ADC variation by combining established theoretical frameworks and compared the results with the measured ADC variation. We acquired data in 8 healthy volunteers with diffusion weighting b = 300 and b = 1000 s/mm².

Results

Apparent diffusion coefficient variation was largest in the feet‐to‐head direction and showed the largest deviation from the mean ADC at peak systole. Artificial ADC variation estimated from tissue deformation was 1.3 ± 0.37·10⁻⁵ mm²/s in the feet‐to‐head direction for gray matter, and 0.75 ± 0.29·10⁻⁵ mm²/s for white matter. The measured ADC variation in the feet‐to‐head direction was 5.6·10⁻⁵ ± 1.5·10⁻⁵ mm²/s for gray matter and 3.2·10⁻⁵ ± 1.0·10⁻⁵ mm²/s for white matter, which was a factor of 3.5 ± 0.82 and 3.4 ± 0.57 larger than the artificial diffusion variations. The measured diffusion variations in the right‐to‐left/anterior‐to‐posterior direction were a factor of 1.5 ± 1.0/1.7 ± 1.4 and 2.0 ± 0.91/2.5 ± 0.94 larger than the artificial diffusion variations for gray matter and white matter, respectively.

Conclusion

Apparent diffusion coefficient variations in the brain likely largely reflect physiology.

Signal to noise and b-value analysis for optimal intra-voxel incoherent motion imaging in the brain

Intravoxel incoherent motion (IVIM) is a method that can provide quantitative information about perfusion in the human body, in vivo, and without contrast agent. Unfortunately, the IVIM perfusion parameter maps are known to be relatively noisy in the brain, in particular for the pseudo-diffusion coefficient, which might hinder its potential broader use in clinical applications. Therefore, we studied the conditions to produce optimal IVIM perfusion images in the brain. IVIM imaging was performed on a 3-Tesla clinical system in four healthy volunteers, with 16 b values 0, 10, 20, 40, 80, 110, 140, 170, 200, 300, 400, 500, 600, 700, 800, 900 s/mm², repeated 20 times. We analyzed the noise characteristics of the trace images as a function of b-value, and the homogeneity of the IVIM parameter maps across number of averages and sub-sets of the acquired b values. We found two peaks of noise of the trace images as function of b value, one due to thermal noise at high b-value, and one due to physiological noise at low b-value. The selection of b value distribution was found to have higher impact on the homogeneity of the IVIM parameter maps than the number of averages. Based on evaluations, we suggest an optimal b value acquisition scheme for a 12 min scan as 0 (7), 20 (4), 140 (19), 300 (9), 500 (19), 700 (1), 800 (4), 900 (1) s/mm².

Intravoxel Incoherent Motion Magnetic Resonance Imaging in Skeletal Muscle: Review and Future Directions

Throughout the body, muscle structure and function can be interrogated using a variety of noninvasive magnetic resonance imaging (MRI) methods. Recently, intravoxel incoherent motion (IVIM) MRI has gained momentum as a method to evaluate components of blood flow and tissue diffusion simultaneously. Much of the prior research has focused on highly vascularized organs, including the brain, kidney, and liver. Unique aspects of skeletal muscle, including the relatively low perfusion at rest and large dynamic range of perfusion between resting and maximal hyperemic states, may influence the acquisition, postprocessing, and interpretation of IVIM data. Here, we introduce several of those unique features of skeletal muscle; review existing studies of IVIM in skeletal muscle at rest, in response to exercise, and in disease states; and consider possible confounds that should be addressed for muscle-specific evaluations. Most studies used segmented nonlinear least squares fitting with a b-value threshold of 200 sec/mm² to obtain IVIM parameters of perfusion fraction (f), pseudo-diffusion coefficient (D*), and diffusion coefficient (D). In healthy individuals, across all muscles, the average ± standard deviation of D was 1.46 ± 0.30 × 10^-3  mm² /sec, D* was 29.7 ± 38.1 × 10^-3  mm² /sec, and f was 11.1 ± 6.7%. Comparisons of reported IVIM parameters in muscles of the back, thigh, and leg of healthy individuals showed no significant difference between anatomic locations. Throughout the body, exercise elicited a positive change of all IVIM parameters. Future directions including advanced postprocessing models and potential sequence modifications are discussed. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY: Stage 2.

Imaging for central nervous system (CNS) interstitial fluidopathy: disorders with impaired interstitial fluid dynamics

After the introduction of the glymphatic system hypothesis, an increasing number of studies on cerebrospinal fluid and interstitial fluid dynamics within the brain have been investigated and reported. A series of diseases are known which develop due to abnormality of the glymphatic system including Alzheimer's disease, traumatic brain injury, stroke, or other disorders. These diseases or disorders share the characteristics of the glymphatic system dysfunction or other mechanisms related to the interstitial fluid dynamics. In this review article, we propose "Central Nervous System (CNS) Interstitial Fluidopathy" as a new concept encompassing diseases whose pathologies are majorly associated with abnormal interstitial fluid dynamics. Categorizing these diseases or disorders as "CNS interstitial fluidopathies," will promote the understanding of their mechanisms and the development of potential imaging methods for the evaluation of the disease as well as clinical methods for disease treatment or prevention. In other words, having a viewpoint of the dynamics of interstitial fluid appears relevant for understanding CNS diseases or disorders, and it would be possible to develop novel common treatment methods or medications for "CNS interstitial fluidopathies."

Correlation of histological marrow characteristics and intravoxel incoherent motion-derived parameters in benign and malignant hematological disorders

PURPOSE

To investigate the relationship between intravoxel incoherent motion (IVIM) parameters and histological parameters of vascularity and cellularity in marrow of hyperplasia disease and compare the difference between benign and malignant marrow disorders.

METHODS

From August 2016 to March 2017, 43 newly diagnosed patients were recruited, which included 15 anemia patients (benign hemopathy) and 28 patients with acute leukemia (AL) (malignant hemopathy). All patients underwent IVIM in the lumbar marrow and the D, D*, f values were measured. The microvessel density (MVD) and cellularity of marrow were calculated from the sample of iliac crest biopsy. Pearson correlation analysis was used to study the relationship between IVIM-derived and histology-derived parameters. We performed unpaired t test to analyze the differences of all parameters between AL and anemia.

RESULT

The MVD was positively correlated with f in patients with AL, anemia and both of them (r = 0.692, P < 0.001; r = 0.595, P = 0.019; r = 0.673, P < 0.001, respectively). But there was no correlation between D* and MVD in three groups. D was not related to bone marrow cellularity (BMC) in all groups. In addition, the f and MVD were higher in AL than anemia (t = 3.546, P = 0.001; t = 6.695, P < 0.001, respectively). The BMC was significantly higher in AL than in anemia (t = 3.330, P = 0.004), but D and D* value had no significant difference between the two groups.

CONCLUSION

The blood volume fraction f was positively correlated with the histological features of marrow in hematological disorders, while f can show the difference of vascularity between benign and malignant marrow disease.

Diffusion tensor imaging of the abdominal organs: Influence of oriented intravoxel flow compartments

Water flow in partially oriented intravoxel compartments mimics an anisotropic fast-diffusion regime, which contributes to the signal attenuation in diffusion-weighted images. In the abdominal organs, this flow may reflect physiological fluid movements (eg, tubular urine flow in kidneys, or bile flow through the liver) and have a clinical relevance. This study investigated the influence of anisotropic intravoxel water flow on diffusion tensor imaging (DTI) of the abdominal organs. Diffusion-weighted images were acquired in five healthy volunteers using an EPI sequence with diffusion preparation (TR/TE: 1000 ms/71 ms; b-values: 0, 10, 20, 40, 70, 120, 250, 450, 700, 1000 s/mm² ; 12 noncollinear diffusion-encoding directions). DTI of liver and kidneys was performed assuming (i) monoexponential decay of the diffusion-weighted signal, and (ii) accounting for potential anisotropy of the fast-diffusion compartments using a tensorial generalization of the IVIM model. Additionally, potential dependency of the metrics of the tensors from the anatomical location was evaluated. Significant differences in the metrics of the diffusion tensor (DT) were found in both liver and kidneys when comparing the two models. In both organs, the trace and the fractional anisotropy of the DT were significantly higher in the monoexponential model than when accounting for perfusion. The comparison of areas of the liver proximal to the hilum with distal regions and of renal cortex with the medulla also proved a location dependency of the size of the fast-diffusion compartments. Pseudo-diffusion correction in DTI enables the assessment of the solid parenchyma regardless of the organ perfusion or other pseudo-diffusive fluid movements. This may have a clinical relevance in the assessment of parenchymal pathologies (eg, liver fibrosis). The fast pseudo-diffusion components present a detectable anisotropy, which may reflect the hepatic microcirculation or other sources of mesoscopic fluid movement in the abdominal organs.

Correction for fast pseudo-diffusive fluid motion contaminations in diffusion tensor imaging

In this prospective study, we quantified the fast pseudo-diffusion contamination by blood perfusion or cerebrospinal fluid (CSF) intravoxel incoherent movements on the measurement of the diffusion tensor metrics in healthy brain tissue. Diffusion-weighted imaging (TR/TE = 4100 ms/90 ms; b-values: 0, 5, 10, 20, 35, 55, 80, 110, 150, 200, 300, 500, 750, 1000, 1300 s/mm², 20 diffusion-encoding directions) was performed on a cohort of five healthy volunteers at 3 Tesla. The projections of the diffusion tensor along each diffusion-encoding direction were computed using a two b-value approach (2b), by fitting the signal to a monoexponential curve (mono), and by correcting for fast pseudo-diffusion compartments using the biexponential intravoxel incoherent motion model (IVIM) (bi). Fractional anisotropy (FA) and mean diffusivity (MD) of the diffusion tensor were quantified in regions of interest drawn over white matter areas, gray matter areas, and the ventricles. A significant dependence of the MD from the evaluation method was found in all selected regions. A lower MD was computed when accounting for the fast-diffusion compartments. A larger dependence was found in the nucleus caudatus (bi: median 0.86 10^-3 mm²/s, Δ2b: -11.2%, Δmono: -14.4%; p = 0.007), in the anterior horn (bi: median 2.04 10^-3 mm²/s, Δ2b: -9.4%, Δmono: -11.5%, p = 0.007) and in the posterior horn of the lateral ventricles (bi: median 2.47 10^-3 mm²/s, Δ2b: -5.5%, Δmono: -11.7%; p = 0.007). Also for the FA, the signal modeling affected the computation of the anisotropy metrics. The deviation depended on the evaluated region with significant differences mainly in the nucleus caudatus (bi: median 0.15, Δ2b: +39.3%, Δmono: +14.7%; p = 0.022) and putamen (bi: median 0.19, Δ2b: +3.1%, Δmono: +17.3%; p = 0.015). Fast pseudo-diffusive regimes locally affect diffusion tensor imaging (DTI) metrics in the brain. Here, we propose the use of an IVIM-based method for correction of signal contaminations through CSF or perfusion.

Cardiac-gated intravoxel incoherent motion diffusion-weighted magnetic resonance imaging for the investigation of intracranial cerebrospinal fluid dynamics in the lateral ventricle: a feasibility study

PURPOSE

Intravoxel incoherent motion (IVIM) in diffusion-weighted magnetic resonance imaging (DW-MRI) attributes the signal attenuation to the molecular diffusion and to a faster pseudo-diffusion. Purpose of the study was to demonstrate the feasibility of IVIM for the investigation of intracranial cerebrospinal fluid (CSF) dynamics.

METHODS

Cardiac-gated DW-MRI images with fifteen b-values (0-1300s/mm²) along three orthogonal directions (mediolateral (ML), anteroposterior (AP), and craniocaudal (CC)) were acquired during maximum systole and diastole in 10 healthy volunteers (6 males, mean age 36 ± 15 years). A pixel-wise bi-exponential fitting with an iterative nonparametric algorithm was carried out to calculate the following parameters: diffusion coefficient (D), fast diffusion coefficient (D*), and fraction of fast diffusion (f). Region of interest measurements were performed in both lateral ventricles. Comparison of IVIM parameters was performed among two cardiac cycle acquisitions and among the diffusion-encoding directions using a paired Student's t test.

RESULTS

f significantly (p < 0.05) depended on the diffusion-encoding direction and on the cardiac cycle (diastole AP 0.30 ± 0.13, ML 0.22 ± 0.12, CC 0.26 ± 0.17; systole AP 0.45 ± 0.17, ML 0.34 ± 0.15, CC 0.40 ± 0.21). Neither a cardiac cycle nor a direction dependency was found among mean D values (which is in line with the expected intraventricular isotropic diffusion) and D* values (p > 0.05 each).

CONCLUSION

The fraction of fast diffusion from IVIM is feasible to detect a direction-dependent and cardiac-dependent pulsatile CSF flow within the lateral ventricles allowing for quantitative monitoring of CSF dynamics. This technique might provide opportunities to further investigate the pathophysiology of various neurological disorders involving altered CSF dynamics.