Umeå University's logo

Background: Infusion testing is an established method for assessing CSF resistance in patients with idiopathic normal pressure hydrocephalus (iNPH). To what extent the increased resistance is related to the glymphatic system is an open question. Here we introduce a computational model that includes the glymphatic system and enables us to determine the importance of (1) brain geometry, (2) intracranial pressure, and (3) physiological parameters on the outcome of and response to an infusion test.

Methods: We implemented a seven-compartment multiple network porous medium model with subject specific geometries from MR images using the finite element library FEniCS. The model consists of the arterial, capillary and venous blood vessels, their corresponding perivascular spaces, and the extracellular space (ECS). Both subject specific brain geometries and subject specific infusion tests were used in the modeling of both healthy adults and iNPH patients. Furthermore, we performed a systematic study of the effect of variations in model parameters.

Results: Both the iNPH group and the control group reached a similar steady state solution when subject specific geometries under identical boundary conditions was used in simulation. The difference in terms of average fluid pressure and velocity between the iNPH and control groups, was found to be less than 6% during all stages of infusion in all compartments. With subject specific boundary conditions, the largest computed difference was a 75% greater fluid speed in the arterial perivascular space (PVS) in the iNPH group compared to the control group. Changes to material parameters changed fluid speeds by several orders of magnitude in some scenarios. A considerable amount of the CSF pass through the glymphatic pathway in our models during infusion, i.e., 28% and 38% in the healthy and iNPH patients, respectively.

Conclusions: Using computational models, we have found the relative importance of subject specific geometries to be less important than individual differences in resistance as measured with infusion tests and model parameters such as permeability, in determining the computed pressure and flow during infusion. Model parameters are uncertain, but certain variations have large impact on the simulation results. The computations resulted in a considerable amount of the infused volume passing through the brain either through the perivascular spaces or the extracellular space.

Background: The pressure difference between the eye and brain in upright postures may be affected by compartmentalization of the optic nerve subarachnoid space (ONSAS). Both pressure and deformation will depend on the microstructures of the ONSAS, and most likely also on ocular glymphatic clearance. Studying these factors could yield important knowledge regarding the translaminar pressure difference, which is suspected to play a role in normal-tension glaucoma.

Methods: A compartment model coupling the ONSAS with the craniospinal CSF system was used to investigate the effects of microstructures on the pressure transfer through the ONSAS during a posture change from supine to upright body postures. ONSAS distensibility was based on MRI measurements. We also included ocular glymphatic flow to investigate how local pressure gradients alter this flow with changes in posture.

Results: A compartmentalization of the ONSAS occurred in the upright posture, with ONSAS porosity (degree of microstructural content) affecting the ONSAS pressure (varying the supine/baseline porosity from 1.0 to 0.75 yielded pressures between − 5.3 and − 2 mmHg). Restricting the minimum computed porosity (occurring in upright postures) to 0.3 prevented compartmentalization, and the ONSAS pressure could equilibrate with subarachnoid space pressure (− 6.5 mmHg) in ≤ 1 h. The ocular glymphatics analysis predicted that substantial intraocular-CSF flows could occur without substantial changes in the ONSAS pressure. The flow entering the ONSAS in supine position (both from the intraocular system and from the cranial subarachnoid space) exited the ONSAS through the optic nerve sheath, while in upright postures the flow through the ONSAS was redirected towards the cranial subarachnoid space.

Conclusions: Microstructures affect pressure transmission along the ONSAS, potentially contributing to ONSAS compartmentalization in upright postures. Different pathways for ocular glymphatic flow were predicted for different postures.

Background: Maintaining cerebral perfusion pressure in the brain when a carotid artery is closed during vascular surgery is critical for avoiding intraoperative hypoperfusion and risk of ischemic stroke. Here we propose and evaluate a method based on computational fluid dynamics for predicting patient-specific cerebral perfusion pressures at carotid clamping during carotid endarterectomy.

Methods: The study consisted of 22 patients with symptomatic carotid stenosis who underwent carotid endarterectomy (73 ± 5 years, 59–80 years, 17 men). The geometry of the circle of Willis was obtained preoperatively from computed tomography angiography and corresponding flow rates from four-dimensional flow magnetic resonance imaging. The patients were also classified as having a present or absent ipsilateral posterior communicating artery based on computed tomography angiography. The predicted mean stump pressures from computational fluid dynamics were compared with intraoperatively measured stump pressures from carotid endarterectomy.

Findings: On group level, there was no difference between the predicted and measured stump pressures (−0.5 ± 13 mmHg, P = 0.86) and the pressures were correlated (r = 0.44, P = 0.039). Omitting two outliers, the correlation increased to r = 0.78 (P < 0.001) (−1.4 ± 8.0 mmHg, P = 0.45). Patients with a present ipsilateral posterior communicating artery (n = 8) had a higher measured stump pressure than those with an absent artery (n = 12) (P < 0.001).

Interpretation: The stump pressure agreement indicates that the computational fluid dynamics approach was promising in predicting cerebral perfusion pressures during carotid clamping, which may prove useful in the preoperative planning of vascular interventions.

An internal carotid artery (ICA) stenosis can potentially decrease the perfusion pressure to the brain. In this study, computational fluid dynamics (CFD) was used to study if there was a hemispheric pressure laterality between the contra- and ipsilateral middle cerebral artery (MCA) in patients with a symptomatic ICA stenosis. We further investigated if this MCA pressure laterality (ΔP_MCA) was related to the hemispheric flow laterality (ΔQ) in the anterior circulation, i.e., ICA, proximal MCA and the proximal anterior cerebral artery (ACA). Twenty-eight patients (73±6 years, range 59–80 years, 21 men) with symptomatic ICA stenosis were included. Flow rates were measured using 4D flow MRI data (PC-VIPR) and vessel geometries were obtained from computed tomography angiography. The ΔP_MCA was calculated from CFD, where patient-specific flow rates were applied at all input- and output boundaries. The ΔP_MCA between the contra- and ipsilateral side was 6.4±8.3 mmHg (p<0.001) (median 3.9 mmHg, range -1.3 to 31.9 mmHg). There was a linear correlation between the ΔP_MCA and ΔQ_ICA (r = 0.85, p<0.001) and ΔQ_ACA (r = 0.71, p<0.001), respectively. The correlation to ΔQ_MCA was weaker (r = 0.47, p = 0.011). In conclusion, the MCA pressure laterality obtained with CFD, is a promising physiological biomarker that can grade the hemodynamic disturbance in patients with a symptomatic ICA stenosis.

PURPOSE: We hypothesize that a collapse of the optic nerve subarachnoid space (ONSAS) in the upright posture may protect the eyes from large translamina cribrosa pressure differences (TLCPD) believed to play a role in various optic nerve diseases (e.g., glaucoma). In this study, we combined magnetic resonance imaging (MRI) and mathematical modeling to investigate this potential ONSAS collapse and its effects on the TLCPD.

METHODS: First, we performed MRI on six healthy volunteers in 6° head-down tilt (HDT) and 13° head-up tilt (HUT) to assess changes in ONSAS volume (measured from the eye to the optic canal) with changes in posture. The volume change reflects optic nerve sheath (ONS) distensibility. Second, we used the MRI data and mathematical modeling to simulate ONSAS pressure and the potential ONSAS collapse in a 90° upright posture.

RESULTS: The MRI showed a 33% decrease in ONSAS volume from the HDT to HUT (P < 0.001). In the upright posture, the simulations predicted an ONSAS collapse 25 mm behind lamina cribrosa, disrupting the pressure communication between the ONSAS and the intracranial subarachnoid space. The collapse reduced the simulated postural increase in TLCPD by roughly 1 mm Hg, although this reduction was highly sensitive to ONS distensibility, varying between 0 and 4.8 mm Hg when varying the distensibility by ± 1 SD.

CONCLUSIONS: The ONSAS volume along the optic nerve is posture dependent. The simulations supported the hypothesized ONSAS collapse in the upright posture and showed that even small changes in ONS stiffness/distensibility may affect the TLCPD.

The central nervous system coordinates movement through forces generated by motor units (MUs) in skeletal muscles. To analyze MUs function is essential in sports, rehabilitation medicine applications, and neuromuscular diagnostics. The MUs and their function are studied using electromyography. Typically, these methods study only a small muscle volume (1 mm³) or only a superficial (< 1 cm) volume of the muscle. Here we introduce a method to identify so-called mechanical units, i.e., the mechanical response of electrically active MUs, in the whole muscle (4x4 cm, cross-sectional) under voluntary contractions by ultrafast ultrasound imaging and spatiotemporal decomposition. We evaluate the performance of the method by simulation of active MUs' mechanical response under weak contractions. We further test the experimental feasibility on eight healthy subjects. We show the existence of mechanical units that contribute to the tissue dynamics in the biceps brachii at low force levels and that these units are similar to MUs described by electromyography with respect to the number of units, territory sizes, and firing rates. This study introduces a new potential neuromuscular functional imaging method, which could be used to study a variety of questions on muscle physiology that previously were difficult or not possible to address.

Background: Infusion testing is a common procedure to determine whether shunting will be beneficial in patients with normal pressure hydrocephalus. The method has a well-developed theoretical foundation and corresponding mathematical models that describe the CSF circulation from the choroid plexus to the arachnoid granulations. Here, we investigate to what extent the proposed glymphatic or paravascular pathway (or similar pathways) modifies the results of the traditional mathematical models.

Methods: We used a compartment model to estimate pressure in the subarachnoid space and the paravascular spaces. For the arachnoid granulations, the cribriform plate and the glymphatic circulation, resistances were calculated and used to estimate pressure and flow before and during an infusion test. Finally, different variations to the model were tested to evaluate the sensitivity of selected parameters.

Results: At baseline intracranial pressure (ICP), we found a very small paravascular flow directed into the subarachnoid space, while 60% of the fluid left through the arachnoid granulations and 40% left through the cribriform plate. However, during the infusion, 80% of the fluid left through the arachnoid granulations, 20% through the cribriform plate and flow in the PVS was stagnant. Resistance through the glymphatic system was computed to be 2.73 mmHg/(mL/min), considerably lower than other fluid pathways, giving non-realistic ICP during infusion if combined with a lymphatic drainage route.

Conclusions: The relative distribution of CSF flow to different clearance pathways depends on ICP, with the arachnoid granulations as the main contributor to outflow. As such, ICP increase is an important factor that should be addressed when determining the pathways of injected substances in the subarachnoid space. Our results suggest that the glymphatic resistance is too high to allow for pressure driven flow by arterial pulsations and at the same time too small to allow for a direct drainage route from PVS to cervical lymphatics.

Neuromuscular diseases hinder muscle function and may be the outcome of damage and dysfunction of the smallest voluntarily activatable units in skeletal muscle tissue, the so-called motor units (MUs). MUs generate electrical signals and analyzing these signals gives a basis to assess and diagnose MUs. The signals are captured using needle electromyography, which is an invasive and non-imaging method. Here, we showultrasound-basedimaging of MUs, via an ultrasound-based spatiotemporal decomposition framework.

Ultrasound medical imaging can be used to visualize and quantify anatomical and functional aspects of internal tissues and organs of the human body. Skeletal muscle tissue is functionally composed by motor units, which are the smallest voluntarily activatable units. In order to capture a transient phenomenon, such as the contraction mechanism, a high sample rate is required. There has been a lot of research on whole-muscle aspects in terms of skeletal muscle contraction characteristics, neuromuscular disorders, and inter-muscle segmentation. Previous studies have shown that small-scale muscle twitches can be detected using ultrasound and there are several reports on ultrasound-based detection of electro-stimulated motor unit activity. However, methods for intra-muscular ultrasound-based analysis of muscle tissue are largely underdeveloped, in particular regarding the level of motor units.Diagnostics of skeletal muscle tissue is based on analyzing features of these units by invasive, non-imaging electrophysiological methods. Here,we make progress by using non-invasive ultrasound imaging to segment motor units, which have the potential to be a non-invasive substitute and where the imaging provides an important contribution.