Modeling cerebrospinal fluid dynamics across the entire intracranial space through integration of four-dimensional flow and intravoxel incoherent motion magnetic resonance imaging.

BACKGROUND: Bidirectional reciprocal motion of cerebrospinal fluid (CSF) was quantified using four-dimensional (4D) flow magnetic resonance imaging (MRI) and intravoxel incoherent motion (IVIM) MRI. To estimate various CSF motions in the entire intracranial region, we attempted to integrate the flow parameters calculated using the two MRI sequences. To elucidate how CSF dynamics deteriorate in Hakim's disease, an age-dependent chronic hydrocephalus, flow parameters were estimated from the two MRI sequences to assess CSF motion in the entire intracranial region. METHODS: This study included 127 healthy volunteers aged ≥ 20 years and 44 patients with Hakim's disease. On 4D flow MRI for measuring CSF motion, velocity encoding was set at 5 cm/s. For the IVIM MRI analysis, the diffusion-weighted sequence was set at six b-values (i.e., 0, 50, 100, 250, 500, and 1000 s/mm2), and the biexponential IVIM fitting method was adapted. The relationships between the fraction of incoherent perfusion (f) on IVIM MRI and 4D flow MRI parameters including velocity amplitude (VA), absolute maximum velocity, stroke volume, net flow volume, and reverse flow rate were comprehensively evaluated in seven locations in the ventricles and subarachnoid spaces. Furthermore, we developed a new parameter for fluid oscillation, the Fluid Oscillation Index (FOI), by integrating these two measurements. In addition, we investigated the relationship between the measurements and indices specific to Hakim's disease and the FOIs in the entire intracranial space. RESULTS: The VA on 4D flow MRI was significantly associated with the mean f-values on IVIM MRI. Therefore, we estimated VA that could not be directly measured on 4D flow MRI from the mean f-values on IVIM MRI in the intracranial CSF space, using the following formula; e0.2(f-85) + 0.25. To quantify fluid oscillation using one integrated parameter with weighting, FOI was calculated as VA × 10 + f × 0.02. In addition, the FOIs at the left foramen of Luschka had the strongest correlations with the Evans index (Pearson's correlation coefficient: 0.78). The other indices related with Hakim's disease were significantly associated with the FOIs at the cerebral aqueduct and bilateral foramina of Luschka. FOI at the cerebral aqueduct was also elevated in healthy controls aged ≥ 60 years. CONCLUSIONS: We estimated pulsatile CSF movements in the entire intracranial CSF space in healthy individuals and patients with Hakim's disease using FOI integrating VA from 4D flow MRI and f-values from IVIM MRI. FOI is useful for quantifying the CSF oscillation.

Automatic assessment of disproportionately enlarged subarachnoid-space hydrocephalus from 3D MRI using two deep learning models.

Background: Disproportionately enlarged subarachnoid-space hydrocephalus (DESH) is a key feature for Hakim disease (idiopathic normal pressure hydrocephalus: iNPH), but subjectively evaluated. To develop automatic quantitative assessment of DESH with automatic segmentation using combined deep learning models. Methods: This study included 180 participants (42 Hakim patients, 138 healthy volunteers; 78 males, 102 females). Overall, 159 three-dimensional (3D) T1-weighted and 180 T2-weighted MRIs were included. As a semantic segmentation, 3D MRIs were automatically segmented in the total ventricles, total subarachnoid space (SAS), high-convexity SAS, and Sylvian fissure and basal cistern on the 3D U-Net model. As an image classification, DESH, ventricular dilatation (VD), tightened sulci in the high convexities (THC), and Sylvian fissure dilatation (SFD) were automatically assessed on the multimodal convolutional neural network (CNN) model. For both deep learning models, 110 T1- and 130 T2-weighted MRIs were used for training, 30 T1- and 30 T2-weighted MRIs for internal validation, and the remaining 19 T1- and 20 T2-weighted MRIs for external validation. Dice score was calculated as (overlapping area) × 2/total area. Results: Automatic region extraction from 3D T1- and T2-weighted MRI was accurate for the total ventricles (mean Dice scores: 0.85 and 0.83), Sylvian fissure and basal cistern (0.70 and 0.69), and high-convexity SAS (0.68 and 0.60), respectively. Automatic determination of DESH, VD, THC, and SFD from the segmented regions on the multimodal CNN model was sufficiently reliable; all of the mean softmax probability scores were exceeded by 0.95. All of the areas under the receiver-operating characteristic curves of the DESH, Venthi, and Sylhi indexes calculated by the segmented regions for detecting DESH were exceeded by 0.97. Conclusion: Using 3D U-Net and a multimodal CNN, DESH was automatically detected with automatically segmented regions from 3D MRIs. Our developed diagnostic support tool can improve the precision of Hakim disease (iNPH) diagnosis.

Development of a 3D Vascular Network Visualization Platform for One-Dimensional Hemodynamic Simulation.

Recent advancements in computational performance and medical simulation technology have made significant strides, particularly in predictive diagnosis. This study focuses on the blood flow simulation reduced-order models, which provide swift and cost-effective solutions for complex vascular systems, positioning them as practical alternatives to 3D simulations in resource-limited medical settings. The paper introduces a visualization platform for patient-specific and image-based 1D-0D simulations. This platform covers the entire workflow, from modeling to dynamic 3D visualization of simulation results. Two case studies on, respectively, carotid stenosis and arterial remodeling demonstrate its utility in blood flow simulation applications.

Deformed Popliteal Artery Due to Highly Flexed Knee Position Can Cause Kinks, Creating an Unfavorable Hemodynamic State.

BACKGROUND: Endovascular treatment devices of the femoropopliteal artery have evolved, improving clinical results. However, the effects of dynamic changes in the popliteal artery during knee flexion have not been sufficiently investigated. In this study we performed a 3-dimensional analysis to clarify the dynamic changes in the popliteal artery during knee flexion and their effects on hemodynamics.Methods and Results: To analyze dynamic changes in the popliteal artery in the knee flexion position, a computed tomography protocol was developed in the right-angled and maximum flexion knee positions. Thirty patients with lower extremity artery disease were recruited. V-Modeler software was used for anatomical and hemodynamic analyses. Various types of deformations of the popliteal artery were revealed, including hinge points and accessory flexions. Kinks can occur in the maximum flexion position; however, they rarely occur in the right-angled flexion position. In addition, hemodynamic analysis revealed a tendency for lower minimum wall shear stress and a higher maximum oscillatory shear index at the maximum curvature of the popliteal artery. CONCLUSIONS: Kinks in the maximum flexion position suggested that the outcome of endovascular treatment may change in areas such as Japan, where knee flexion is customary. Hemodynamics at the maximum curvature of the popliteal artery indicated that the luminal condition was unfavorable for endovascular treatment.

Development of a Model for Plaque Induction in Rat Carotid Arteries.

Objective Plaque induction through intimal injury using a balloon catheter in small animals and by artificial ligation of the carotid artery in large animals have been reported. However, these reports have not yet succeeded in inducing stable plaques nor creating a high degree of intimal thickening to be used as animal models. We have previously developed a plaque induction model in rats but have failed to obtain a plaque incidence frequency that can be used as a model. Thus, in the current study, we aimed to create a versatile disease model to examine the pharmacokinetics of drug administration, determine the efficacy of treatment, and examine the process of intimal thickening. We also attempted to create an improved model with shorter, more frequent, and more severe intimal thickening. Materials and Methods The common carotid artery of male Wistar rats was surgically exposed and completely ligated with a wire and 6-0 nylon thread. Then, the wire was removed to create a partial ligation. To create a high frequency and high degree of intimal thickening, 72 rats were divided into two groups: a single lesion group with a 0.25-mm wire and a single ligature point, and a tandem lesion group with a 0.3-mm wire and two ligature points. Each group was further divided into normal diet and high cholesterol diet groups. The presence and frequency of intimal thickening were examined for each group after 4, 8, and 16 weeks of growth. Results In the single lesion group, intimal thickening was observed in 42% of the 4-week group and 75% of the 8-week group. In the tandem lesion group, intimal thickening was observed in 75% of the 4-week group and 50% of the 8-week group. In addition, 50% of the individuals reared for 16 weeks developed intimal thickening. Conclusion We successfully induced intimal thickening in the carotid arteries of rats with high frequency in the single lesion and tandem lesion groups. The results also showed that the tandem lesion group tended to induce intimal thickening earlier than the single lesion group.

Aging and Sex Differences in Brain Volume and Cerebral Blood Flow.

How do regional brain volume ratios and cerebral blood flow (CBF, mL/min) change with aging, and are there sex differences? This study aimed to comprehensively evaluate the relationships between regional brain volume ratios and CBF in healthy brains. The study participants were healthy volunteers who underwent three-dimensional T1-weighted MRI, time-of-flight MR angiography, and four-dimensional (4D) flow MRI between 2020 and 2022. The brain was automatically segmented into 21 brain subregions from 3D T1-weighted MRI, and CBF in 16 major intracranial arteries were measured by 4D flow MRI. The relationships between segmented brain volume ratios and CBFs around the circle of Willis were comprehensively investigated in each decade and sex. This study included 129 healthy volunteers (mean age ± SD, 48.2 ± 16.8; range, 22-92 years; 43 males and 86 females). The association was strongest between the cortical gray matter volume ratio and total outflow of the intracranial major arteries distal to the circle of Willis (Pearson's correlation coefficient, r: 0.425). In addition, the mean flow of the total inflow and outflow around the circle of Willis were significantly greater in women than men, and significant left-right differences were observed in CBFs even on the peripheral side of the circle of Willis. Moreover, the correlation was strongest between the left cortical gray matter volume ratio and the combined flows of the left anterior and posterior cerebral arteries distal to the circle of Willis (r: 0.486). There was a trend toward greater total intracranial CBF, especially among women in their 40s and younger, who had a larger cortical gray matter volume. This finding may be one of the reasons for the approximately twofold higher incidence of cerebral aneurysms and subarachnoid hemorrhage, and a threefold higher incidence of migraine headaches.

Computational modeling of multiscale collateral blood supply in a whole-brain-scale arterial network.

The cerebral arterial network covering the brain cortex has multiscale anastomosis structures with sparse intermediate anastomoses (O[102] µm in diameter) and dense pial networks (O[101] µm in diameter). Recent studies indicate that collateral blood supply by cerebral arterial anastomoses has an essential role in the prognosis of acute ischemic stroke caused by large vessel occlusion. However, the physiological importance of these multiscale morphological properties-and especially of intermediate anastomoses-is poorly understood because of innate structural complexities. In this study, a computational model of multiscale anastomoses in whole-brain-scale cerebral arterial networks was developed and used to evaluate collateral blood supply by anastomoses during middle cerebral artery occlusion. Morphologically validated cerebral arterial networks were constructed by combining medical imaging data and mathematical modeling. Sparse intermediate anastomoses were assigned between adjacent main arterial branches; the pial arterial network was modeled as a dense network structure. Blood flow distributions in the arterial network during middle cerebral artery occlusion simulations were computed. Collateral blood supply by intermediate anastomoses increased sharply with increasing numbers of anastomoses and provided one-order-higher flow recoveries to the occluded region (15%-30%) compared with simulations using a pial network only, even with a small number of intermediate anastomoses (≤10). These findings demonstrate the importance of sparse intermediate anastomoses, which are generally considered redundant structures in cerebral infarction, and provide insights into the physiological significance of the multiscale properties of arterial anastomoses.

Age-Related Changes in Cerebrospinal Fluid Dynamics in the Pathogenesis of Chronic Hydrocephalus in Adults.

Cerebrospinal fluid (CSF) dynamics has dramatically changed in this century. In the latest concept of CSF dynamics, CSF is thought to be produced mainly from interstitial fluid excreted from the brain parenchyma and is absorbed in the meningeal lymphatics. Moreover, CSF does not always flow from the ventricles to the subarachnoid space unidirectionally through the foramina of Magendie and Luschka. In an environment of increased intracranial CSF in idiopathic normal pressure hydrocephalus, CSF freely moves through the inferior choroidal point of the choroidal fissure, which interfaces between the inferior horn of the lateral ventricles and the ambient cistern and through the velum interpositum between the third ventricle and the quadrigeminal cistern. The structure of the hippocampus adjacent to the inferior part of the choroidal fissure may be important in preventing the accumulation of waste products in the hippocampus. A recent imaging technology for CSF dynamics, such as four-dimensional flow and intravoxel incoherent motion magnetic resonance imaging, can visualize and quantify the pulsatile complex CSF motion in clinical usage. We present the current concepts of CSF dynamics with advanced magnetic resonance imaging techniques, which will be helpful in the management and understanding of the pathogenesis of chronic hydrocephalus in adults.

Tightened Sulci in the High Convexities as a Noteworthy Feature of Idiopathic Normal Pressure Hydrocephalus.

OBJECTIVE: The presence of tightened sulci in the high-convexities (THC) is a key morphological feature for the diagnosis of idiopathic normal pressure hydrocephalus (iNPH), but the exact localization of THC has yet to be defined. The purpose of this study was to define THC and compare its volume, percentage, and index between iNPH patients and healthy controls. METHODS: According to the THC definition, the high-convexity part of the subarachnoid space was segmented and measured the volume and percentage from the 3D T1-weighted and T2-weighted magnetic resonance images in 43 patients with iNPH and 138 healthy controls. RESULTS: THC was defined as a decrease in the high-convexity part of the subarachnoid space located above the body of the lateral ventricles, with anterior end on the coronal plane perpendicular to the anterior commissure-posterior commissure (AC-PC) line passing through the front edge of the genu of corpus callosum, the posterior end in the bilateral posterior parts of the callosomarginal sulci, and the lateral end at 3 cm from the midline on the coronal plane perpendicular to the AC-PC line passing through the midpoint between AC and PC. Compared to the volume and volume percentage, the high-convexity part of the subarachnoid space volume per ventricular volume ratio < 0.6 was the most detectable index of THC on both 3D T1-weighted and T2-weighted magnetic resonance images. CONCLUSIONS: To improve the diagnostic accuracy of iNPH, the definition of THC was clarified, and high-convexity part of the subarachnoid space volume per ventricular volume ratio <0.6 proposed as the best index for THC detection in this study.

Aging-related volume changes in the brain and cerebrospinal fluid using artificial intelligence-automated segmentation.

OBJECTIVES: To verify the reliability of the volumes automatically segmented using a new artificial intelligence (AI)-based application and evaluate changes in the brain and CSF volume with healthy aging. METHODS: The intracranial spaces were automatically segmented in the 21 brain subregions and 5 CSF subregions using the AI-based application on the 3D T1-weighted images in healthy volunteers aged > 20 years. Additionally, the automatically segmented volumes of the total ventricles and subarachnoid spaces were compared with the manually segmented volumes of those extracted from 3D T2-weighted images using the intra-class correlation and Bland-Altman analysis. RESULTS: In this study, 133 healthy volunteers aged 21-92 years were included. The mean intra-class correlations between the automatically and manually segmented volumes of the total ventricles and subarachnoid spaces were 0.986 and 0.882, respectively. The increase in the CSF volume was estimated to be approximately 30 mL (2%) per decade from 265 mL (18.7%) in the 20s to 488 mL (33.7%) in ages above 80 years; however, the increase in the volume of total ventricles was approximately 20 mL (< 2%) until the 60s and increased in ages above 60 years. CONCLUSIONS: This study confirmed the reliability of the CSF volumes using the AI-based auto-segmentation application. The intracranial CSF volume increased linearly because of the brain volume reduction with aging; however, the ventricular volume did not change until the age of 60 years and above and then gradually increased. This finding could help elucidate the pathogenesis of chronic hydrocephalus in adults. KEY POINTS: • The brain and CSF spaces were automatically segmented using an artificial intelligence-based application. • The total subarachnoid spaces increased linearly with aging, whereas the total ventricle volume was around 20 mL (< 2%) until the 60s and increased in ages above 60 years. • The cortical gray matter gradually decreases with aging, whereas the subcortical gray matter maintains its volume, and the cerebral white matter increases slightly until the 40s and begins to decrease from the 50s.

Usefulness of intravoxel incoherent motion MRI for visualizing slow cerebrospinal fluid motion.

BACKGROUND: In the cerebrospinal fluid (CSF) dynamics, the pulsations of cerebral arteries and brain is considered the main driving force for the reciprocating bidirectional CSF movements. However, measuring these complex CSF movements on conventional flow-related MRI methods is difficult. We tried to visualize and quantify the CSF motion by using intravoxel incoherent motion (IVIM) MRI with low multi-b diffusion-weighted imaging. METHODS: Diffusion-weighted sequence with six b values (0, 50, 100, 250, 500, and 1000 s/mm2) was performed on 132 healthy volunteers aged ≥ 20 years and 36 patients with idiopathic normal pressure hydrocephalus (iNPH). The healthy volunteers were divided into three age groups (< 40, 40 to < 60, and ≥ 60 years). In the IVIM analysis, the bi-exponential IVIM fitting method using the Levenberg-Marquardt algorithm was adapted. The average, maximum, and minimum values of ADC, D, D*, and fraction of incoherent perfusion (f) calculated by IVIM were quantitatively measured in 45 regions of interests in the whole ventricles and subarachnoid spaces. RESULTS: Compared with healthy controls aged ≥ 60 years, the iNPH group had significantly lower mean f values in all the parts of the lateral and 3rd ventricles, whereas significantly higher mean f value in the bilateral foramina of Luschka. In the bilateral Sylvian fossa, which contain the middle cerebral bifurcation, the mean f values increased gradually with increasing age, whereas those were significantly lower in the iNPH group. In the 45 regions of interests, the f values in the bilateral foramina of Luschka were the most positively correlated with the ventricular size and indices specific to iNPH, whereas that in the anterior part of the 3rd ventricle was the most negatively correlated with the ventricular size and indices specific to iNPH. Other parameters of ADC, D, and D* were not significantly different between the two groups in any locations. CONCLUSIONS: The f value on IVIM MRI is useful for evaluating small pulsatile complex motion of CSF throughout the intracranial CSF spaces. Patients with iNPH had significantly lower mean f values in the whole lateral ventricles and 3rd ventricles and significantly higher mean f value in the bilateral foramina of Luschka, compared with healthy controls aged ≥ 60 years.

Uncertainty quantification in cerebral circulation simulations focusing on the collateral flow: Surrogate model approach with machine learning.

Collateral circulation in the circle of Willis (CoW), closely associated with disease mechanisms and treatment outcomes, can be effectively investigated using one-dimensional-zero-dimensional hemodynamic simulations. As the entire cardiovascular system is considered in the simulation, it captures the systemic effects of local arterial changes, thus reproducing collateral circulation that reflects biological phenomena. The simulation facilitates rapid assessment of clinically relevant hemodynamic quantities under patient-specific conditions by incorporating clinical data. During patient-specific simulations, the impact of clinical data uncertainty on the simulated quantities should be quantified to obtain reliable results. However, as uncertainty quantification (UQ) is time-consuming and computationally expensive, its implementation in time-sensitive clinical applications is considered impractical. Therefore, we constructed a surrogate model based on machine learning using simulation data. The model accurately predicts the flow rate and pressure in the CoW in a few milliseconds. This reduced computation time enables the UQ execution with 100 000 predictions in a few minutes on a single CPU core and in less than a minute on a GPU. We performed UQ to predict the risk of cerebral hyperperfusion (CH), a life-threatening condition that can occur after carotid artery stenosis surgery if collateral circulation fails to function appropriately. We predicted the statistics of the postoperative flow rate increase in the CoW, which is a measure of CH, considering the uncertainties of arterial diameters, stenosis parameters, and flow rates measured using the patients' clinical data. A sensitivity analysis was performed to clarify the impact of each uncertain parameter on the flow rate increase. Results indicated that CH occurred when two conditions were satisfied simultaneously: severe stenosis and when arteries of small diameter serve as the collateral pathway to the cerebral artery on the stenosis side. These findings elucidate the biological aspects of cerebral circulation in terms of the relationship between collateral flow and CH.

Computational fluid dynamic analysis of the nasal respiratory function before and after postero-superior repositioning of the maxilla.

Morphological changes in the upper airway and the resulting alteration in the nasal respiratory function after jawbone repositioning during orthognathic surgery have garnered attention recently. In particular, nasopharyngeal stenosis, because of the complex influence of both jaws, the effects of which have not yet been clarified owing to postero-superior repositioning of the maxilla, may significantly impact sleep and respiratory function, necessitating further functional evaluation. This study aimed to perform a functional evaluation of the effects of surgery involving maxillary repositioning, which may result in a larger airway resistance if the stenosis worsens the respiratory function, using CFD for treatment planning. A model was developed from CT images obtained preoperatively (PRE) and postoperatively (POST) in females (n = 3) who underwent maxillary postero-superior repositioning using Mimics and ICEM CFD. Simultaneously, a model of stenosis (STENOSIS) was developed by adjusting the severity of stenosis around the PNS to simulate greater repositioning than that in the POST. Inhalation at rest and atmospheric pressure were simulated in each model using Fluent, whereas pressure drop (ΔP) was evaluated using CFD Post. In this study, ΔP was proportional to airway resistance because the flow rate was constant. Therefore, the magnitude of ΔP was evaluated as the level of airway resistance. The ΔP in the airway was lower in the POST compared to the PRE, indicating that the analysis of the effects of repositioning on nasal ventilation showed that current surgery is appropriate with respect to functionality, as it does not compromise respiratory function. The rate of change in the cross-sectional area of the mass extending pharynx (α) was calculated as the ratio of each neighboring section. The closer the α-value is to 1, the smaller the ΔP, so ideally the airway should be constant. This study identified airway shapes that are favorable from the perspective of fluid dynamics.

Numerical Investigation of T-Shaped Microfluidic Oscillator with Viscoelastic Fluid.

Oscillatory flow has many applications in micro-scaled devices. The methods of realizing microfluidic oscillators reported so far are typically based on the impinging-jet and Coanda effect, which usually require the flow Reynolds number to be at least at the order of unity. Another approach is to introduce elastomeric membrane into the microfluidic units; however, the manufacturing process is relatively complex, and the membrane will become soft after long-time operation, which leads to deviation from the design condition. From the perspective of the core requirement of a microfluidic circuit, i.e., nonlinearity, the oscillatory microfluidic flow can be realized via the nonlinear characteristics of viscoelastic fluid flow. In this paper, the flow characteristics of viscoelastic fluid (Boger-type) in a T-shaped channel and its modified structures are studied by two-dimensional direct numerical simulation (DNS). The main results obtained from the DNS study are as follows: (1) Both Weissenberg (Wi) number and viscosity ratio need to be within a certain range to achieve a periodic oscillating performance; (2) With the presence of the dynamic evolution of the pair of vortices in the upstream near the intersection, the oscillation intensity increases as the elasticity-dominated area in the junction enlarges; (3) Considering the simplicity of the T-type channel as a potential oscillator, the improved structure should have a groove carved toward the entrance near the upper wall. The maximum oscillation intensity measured by the standard deviation of flow rate at outlet is increased by 129% compared with that of the original standard T-shaped channel under the same condition. To sum up, with Wi number and viscosity ratio within a certain range, the regular periodic oscillation characteristics of Oldroyd-B type viscoelastic fluid flow in standard T-shaped and its modified channels can be obtained. This structure can serve as a passive microfluidic oscillator with great potential value at an extremely low Reynolds number, which has the advantages of simplicity, no moving parts and fan-out of two.

Carotid Stump Pressure and Contralateral Internal Carotid Stenosis Ratio During Carotid Endarterectomies: 1D-0D Hemodynamic Simulation of Cerebral Perfusion.

Objective: We selectively place carotid shunting when ipsilateral mean stump pressure is less than 40 mmHg during carotid endarterectomy (CEA). This study aimed to assess the validity of our selective shunting criterion by 1D-0D hemodynamic simulation technology. Materials and Methods: We retrospectively reviewed 88 patients (95 cases) of CEA and divided them into two groups based on the degree of contralateral internal carotid artery (ICA) stenosis ratio, which was determined as severe when the peak systolic velocity ratio of the ICA to the common carotid artery was ≥4 by carotid duplex ultrasonography. Patients with severe stenosis or occlusion in contralateral ICA were classified as hypoperfusion group, and those without such contralateral ICA obstruction were classified as control group. Results: Perioperatively, the mean carotid stump pressures were 33 mmHg in hypoperfusion group and 46 mmHg in the control group (P=0.006). We simulated changes in carotid stump pressure according to the changes in the contralateral ICA stenosis ratio. 1D-0D simulation indicated a sharp decline in carotid stump pressure when the contralateral stenosis ratio was >50%, while peripheral pressure of the middle cerebral arteries declined sharply at a ≥70% contralateral stenosis ratio. At this ratio, the direction of the ipsilateral cerebral arterial flow became inverted, the carotid stump pressure became dependent on the basilar artery circulation, and the ipsilateral middle cerebral artery became hypoperfused. Conclusion: Our clinical and computer-simulated results confirmed the validation of our carotid shunting criterion and suggested that contralateral ICA stenosis ratio over 70% is a safe indication of selective shunting during CEA.

A penalized spline fitting method to optimize geometric parameters of arterial centerlines extracted from medical images.

In order to grasp the spatial and temporal evolution of vascular geometry, three-dimensional (3D) arterial bending structure and geometrical changes of arteries and stent grafts (SG) must be quantified using geometrical parameters such as curvature and torsion along the vasculature centerlines extracted from medical images. Here, we develop a robust method for constructing smooth centerlines based on a spline fitting method (SFM) such that the optimized geometric parameters of curvature and torsion can be obtained independently of digitization noise in the images. Conventional SFM consists of the 3rd degree spline basis function and 2nd derivative penalty term. In contrast, the present SFM uses the 5th degree spline basis function and 3rd and 4th derivative penalty terms, the coefficients of which are derived by the Akaike information criterion. The results show that the developed SFM can reduce the errors of curvature and torsion compared to conventional SFM. We then apply the present SFM to the centerline of the SG in an abdominal aortic aneurysm (AAA), and those of bilateral internal carotid arteries (ICA) in 6 cases: 3 cases with aneurysms and 3 cases without any aneurysm. The SG centerlines were obtained from temporal medical images at three scan times. The strong peak of the curvature could be clearly observed in the distal area of the SG, the inversion of the torsion at 0 months in the middle area of SG disappeared over time, and the torsions around the SG bifurcation at the three time periods were inverted. The curvature-torsion graphs along the ICA centerlines superimposing five aneurysmal positions were useful for investigating the relationship between arterial bending structure and aneurysmal positions. Both ICAs had curvature peak values higher than 0.4 within the ICA syphons. The ICA torsion graphs indicated that left and right ICA tended to be a right- and left-handed helix, respectively. In the left ICA syphon, the biggest aneurysm could be observed downstream of the salient torsion inversion. All aneurysms for 3 cases were positioned at the downstream of the inverted torsion.

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.

Reliability and Interobserver Variability of Evans' Index and Disproportionately Enlarged Subarachnoid Space Hydrocephalus as Diagnostic Criteria for Idiopathic Normal Pressure Hydrocephalus.

BACKGROUND: The image diagnosis of idiopathic normal-pressure hydrocephalus (iNPH) is based on the ventriculomegaly, whose criterion is an Evans' Index (EI) >0.3. Recently, disproportionately enlarged subarachnoid space hydrocephalus (DESH) has been proposed as a morphological characteristic to iNPH. Several studies cast doubt on the reliability of these criteria in the diagnosis of iNPH. Furthermore, interobserver differences of these criteria have not yet been investigated. The objective of this study was to assess the diagnostic reliability and interobserver variability of EI and DESH. MATERIALS AND METHODS: The preoperative magnetic resonance (MR) images of 84 definite iNPH patients were retrospectively evaluated by a neuroradiologist (NR) and physical therapist (PT). They independently assessed the EI and DESH. The MR images were evaluated preoperatively by a neurosurgeon (NS). The results were showed in mean (standard deviation). RESULTS: The mean age was 78.4 (6.3) years (male:female = 49:35). The mean EI was 0.33 (0.04), 0.32 (0.04), and 0.31 (0.03) for NS, NR, and PT, respectively (P < 0.0001). The rate of accurate diagnosis of iNPH with EI >0.3 was 74%, 66%, and 61% for NS, NR, and PT, respectively, and there was a moderate level of agreement. By contrast, there was a substantial lower level of accuracy in assessment with DESH for all three evaluators as 50%, 44%, and 27% for NS, NR, and PT, respectively, again with a moderate level of agreement. However, the rates of patients fulfilling both EI >0.3 and DESH were remarkably lower than either of the two parameters individually at a mere 37%, 30%, and 16% for NS, NR, and PT, respectively, with a low level of agreement between the rates. CONCLUSION: This study suggests that DESH cannot be a diagnostic criterion for iNPH. If EI >0.3 and DESH were both necessary to diagnose iNPH, then more than 70% of patients would have been misdiagnosed and would have been deprived of the chance of treatment and its benefits. These results request a paradigm shift in the concepts of iNPH.

Computational simulation of flow-induced arterial remodeling of the pancreaticoduodenal arcade associated with celiac artery stenosis.

Arterial remodeling of the pancreaticoduodenal arcade, which enables collateral flow to the liver, spleen, and stomach, is a well-recognized clinical sign of celiac artery (CA) stenosis. However, the hemodynamic changes due to remodeling are poorly understood, despite their importance in surgical procedures such as pancreaticoduodenectomy. In this study, a framework to simulate remodeling of the arterial network following pathological flow alterations was developed and applied to investigate the hemodynamic characteristics of patients with CA stenosis. A one-dimensional-zero-dimensional cardiovascular model was used for blood flow simulation. After introducing CA stenosis into the normal network, arterial remodeling was simulated by iteratively changing the diameter of each artery until time-averaged wall shear stress reached its value under normal conditions. A representative case was simulated to validate the present framework, followed by simulation cases to investigate the impact of stenosis severity on remodeling outcome. A markedly dilated arcade was observed whose diameter agreed well with the corresponding values measured in subjects with CA stenosis, confirming the ability of the framework to predict arterial remodeling. A series of simulations clarified how the geometry and hemodynamics after remodeling change with stenosis severity. In particular, the arterial remodeling and resulting blood flow redistribution were found to maintain adequate organ blood supply regardless of stenosis severity. Furthermore, it was suggested that flow conditions in patients with CA stenosis could be estimated from geometric factors, namely, stenosis severity and arcade diameter, which can be preoperatively and non-invasively measured using diagnostic medical images.