Assessment of liver function by gadoxetic acid avidity in MRI in a model of rapid liver regeneration in rats

BACKGROUND

This animal study investigates the hypothesis of an immature liver growth following ALPPS (associating liver partition and portal vein ligation for staged hepatectomy) by measuring liver volume and function using gadoxetic acid avidity in magnetic resonance imaging (MRI) in models of ALPPS, major liver resection (LR) and portal vein ligation (PVL).

METHODS

Wistar rats were randomly allocated to ALPPS, LR or PVL. In contrast-enhanced MRI scans with gadoxetic acid (Primovist®), liver volume and function of the right median lobe (=future liver remnant, FLR) and the deportalized lobes (DPL) were assessed until post-operative day (POD) 5. Liver functionFLR/DPL was defined as the inverse value of time from injection of gadoxetic acid to the blood pool-corrected maximum signal intensityFLR/DPL multiplied by the volumeFLR/DPL.

RESULTS

In ALPPS (n = 6), LR (n = 6) and PVL (n = 6), volumeFLR and functionFLR increased proportionally, except on POD 1. Thereafter, functionFLR exceeded volumeFLR increase in LR and ALPPS, but not in PVL. Total liver function was significantly reduced after LR until POD 3, but never undercuts 60% of its pre-operative value following ALPPS and PVL.

DISCUSSION

This study shows for the first time that functional increase is proportional to volume increase in ALPPS using gadoxetic acid avidity in MRI.

Influence of age on the relation between body position and noninvasively acquired intracranial pulse waves

The capacitive measurement of the head's dielectric properties has been recently proposed as a noninvasive method for deriving surrogates of craniospinal compliance (CC), a parameter used in the evaluation of space-occupying neurological disorders. With the higher prevalence of such disorders in the older compared to the younger population, data on the head's dielectric properties of older healthy individuals would be of particularly high value before assessing pathologic changes. However, so far only measurements on young volunteers (< 30 years) were reported. In the present study, we have investigated the capacitively obtained electric signal known as W in older healthy individuals. Thirteen healthy subjects aged > 60 years were included in the study. W was acquired in the resting state (supine horizontal position), and during head-up and head-down tilting. AMP, the peak-to-valley amplitude of W related to cardiac action, was extracted from W. AMP was higher in this older cohort compared to the previously investigated younger one (0°: 5965 ± 1677 arbitrary units (au)). During head-up tilting, AMP decreased (+ 60°: 4446 ± 1620 au, P < 0.001), whereas it increased during head-down tilting (- 30°: 7600 ± 2123 au, P < 0.001), as also observed in the younger cohort. Our observation that AMP, a metric potentially reflective of CC, is higher in the older compared to the younger cohort aligns with the expected decrease of CC with age. Furthermore, the robustness of AMP is reinforced by the consistent relative changes observed during tilt testing in both cohorts.

Current Understanding of the Anatomy, Physiology, and Magnetic Resonance Imaging of Neurofluids: Update From the 2022 "ISMRM Imaging Neurofluids Study group" Workshop in Rome

Neurofluids is a term introduced to define all fluids in the brain and spine such as blood, cerebrospinal fluid, and interstitial fluid. Neuroscientists in the past millennium have steadily identified the several different fluid environments in the brain and spine that interact in a synchronized harmonious manner to assure a healthy microenvironment required for optimal neuroglial function. Neuroanatomists and biochemists have provided an incredible wealth of evidence revealing the anatomy of perivascular spaces, meninges and glia and their role in drainage of neuronal waste products. Human studies have been limited due to the restricted availability of noninvasive imaging modalities that can provide a high spatiotemporal depiction of the brain neurofluids. Therefore, animal studies have been key in advancing our knowledge of the temporal and spatial dynamics of fluids, for example, by injecting tracers with different molecular weights. Such studies have sparked interest to identify possible disruptions to neurofluids dynamics in human diseases such as small vessel disease, cerebral amyloid angiopathy, and dementia. However, key differences between rodent and human physiology should be considered when extrapolating these findings to understand the human brain. An increasing armamentarium of noninvasive MRI techniques is being built to identify markers of altered drainage pathways. During the three-day workshop organized by the International Society of Magnetic Resonance in Medicine that was held in Rome in September 2022, several of these concepts were discussed by a distinguished international faculty to lay the basis of what is known and where we still lack evidence. We envision that in the next decade, MRI will allow imaging of the physiology of neurofluid dynamics and drainage pathways in the human brain to identify true pathological processes underlying disease and to discover new avenues for early diagnoses and treatments including drug delivery. Evidence level: 1 Technical Efficacy: Stage 3.

Incorporating Unresolved Stresses in Blood Damage Modeling: Energy Dissipation More Accurate Than Reynolds Stress Formulation

OBJECTIVE

Reynolds Averaged Navier Stokes (RANS) models are often used as the basis for modeling blood damage in turbulent flows. To predict blood damage by turbulence stresses that are not resolved in RANS, a stress formulation that represents the corresponding scales is required. Here, we compare two commonly employed stress formulations: a scalar stress representation that uses Reynolds stresses as a surrogate for unresolved fluid stresses, and an effective stress formulation based on energy dissipation.

METHODS

We conducted unsteady RANS simulations of the CentriMag blood pump with three different closure models and a Large Eddy Simulation (LES) for reference. We implemented both stress representations in all models and compared the resulting total stress distributions in Eulerian and Lagrangian frameworks.

RESULTS

The Reynolds-stress-based approach overestimated the contribution of unresolved stresses in RANS, with differences between closure models of up to several orders of magnitude. With the dissipation-based approach, the total stresses predicted with RANS deviated by about 50% from the LES reference, which was more accurate than only considering resolved stresses.

CONCLUSION

The Reynolds-stress-based formulation proved unreliable for estimating scalar stresses in our RANS simulations, while the dissipation-based approach provided an accuracy improvement over simply neglecting unresolved stresses.

SIGNIFICANCE

Our results suggest that dissipation-based inclusion of unresolved stresses should be the preferred choice for blood damage modeling in RANS.

Evolution of hypoxia and hypoxia-inducible factor asparaginyl hydroxylase regulation in chronic kidney disease

BACKGROUND

The roles of hypoxia and hypoxia inducible factor (HIF) during chronic kidney disease (CKD) are much debated. Interventional studies with HIF-α activation in rodents have yielded contradictory results. The HIF pathway is regulated by prolyl and asparaginyl hydroxylases. While prolyl hydroxylase inhibition is a well-known method to stabilize HIF-α, little is known about the effect asparaginyl hydroxylase factor inhibiting HIF (FIH).

METHODS

We used a model of progressive proteinuric CKD and a model of obstructive nephropathy with unilateral fibrosis. In these models we assessed hypoxia with pimonidazole and vascularization with three-dimensional micro-computed tomography imaging. We analysed a database of 217 CKD biopsies from stage 1 to 5 and we randomly collected 15 CKD biopsies of various severity degrees to assess FIH expression. Finally, we modulated FIH activity in vitro and in vivo using a pharmacologic approach to assess its relevance in CKD.

RESULTS

In our model of proteinuric CKD, we show that early CKD stages are not characterized by hypoxia or HIF activation. At late CKD stages, some areas of hypoxia are observed, but these are not colocalizing with fibrosis. In mice and in humans, we observed a downregulation of the HIF pathway, together with an increased FIH expression in CKD, according to its severity. Modulating FIH in vitro affects cellular metabolism, as described previously. In vivo, pharmacologic FIH inhibition increases the glomerular filtration rate of control and CKD animals and is associated with decreased development of fibrosis.

CONCLUSIONS

The causative role of hypoxia and HIF activation in CKD progression is questioned. A pharmacological approach of FIH downregulation seems promising in proteinuric kidney disease.

Terabyte-scale supervised 3D training and benchmarking dataset of the mouse kidney

The performance of machine learning algorithms, when used for segmenting 3D biomedical images, does not reach the level expected based on results achieved with 2D photos. This may be explained by the comparative lack of high-volume, high-quality training datasets, which require state-of-the-art imaging facilities, domain experts for annotation and large computational and personal resources. The HR-Kidney dataset presented in this work bridges this gap by providing 1.7 TB of artefact-corrected synchrotron radiation-based X-ray phase-contrast microtomography images of whole mouse kidneys and validated segmentations of 33 729 glomeruli, which corresponds to a one to two orders of magnitude increase over currently available biomedical datasets. The image sets also contain the underlying raw data, threshold- and morphology-based semi-automatic segmentations of renal vasculature and uriniferous tubules, as well as true 3D manual annotations. We therewith provide a broad basis for the scientific community to build upon and expand in the fields of image processing, data augmentation and machine learning, in particular unsupervised and semi-supervised learning investigations, as well as transfer learning and generative adversarial networks.

Radiological feature heterogeneity supports etiological diversity among patient groups in Meniere's disease

We aimed to determine the prevalence of radiological temporal bone features that in previous studies showed only a weak or an inconsistent association with the clinical diagnosis of Meniere's disease (MD), in two groups of MD patients (n = 71) with previously established distinct endolymphatic sac pathologies; i.e. the group MD-dg (ES degeneration) and the group MD-hp (ES hypoplasia). Delayed gadolinium-enhanced MRI and high-resolution CT data were used to determine and compare between and within (affected vs. non-affected side) groups geometric temporal bone features (lengths, widths, contours), air cell tract volume, height of the jugular bulb, sigmoid sinus width, and MRI signal intensity alterations of the ES. Temporal bone features with significant intergroup differences were the retrolabyrinthine bone thickness (1.04 ± 0.69 mm, MD-hp; 3.1 ± 1.9 mm, MD-dg; p < 0.0001); posterior contour tortuosity (mean arch-to-chord ratio 1.019 ± 0.013, MD-hp; 1.096 ± 0.038, MD-dg; p < 0.0001); and the pneumatized volume (1.37 [0.86] cm³, MD-hp; 5.25 [3.45] cm³, MD-dg; p = 0.03). Features with differences between the affected and non-affected sides within the MD-dg group were the sigmoid sinus width (6.5 ± 1.7 mm, affected; 7.6 ± 2.1 mm, non-affected; p = 0.04) and the MRI signal intensity of the endolymphatic sac (median signal intensity, affected vs. unaffected side, 0.59 [IQR 0.31-0.89]). Radiological temporal bone features known to be only weakly or inconsistently associated with the clinical diagnosis MD, are highly prevalent in either of two MD patient groups. These results support the existence of diverse-developmental and degenerative-disease etiologies manifesting with distinct radiological temporal bone abnormalities.

Assessment of extracranial carotid artery disease using digital twins - A pilot study

To improve risk stratification in extracranial internal carotid artery disease (CAD), patients who would benefit maximally from revascularization must be identified. In cardiology, the fractional flow reserve (FFR) has become a reference standard for evaluating the functional severity of coronary artery stenosis, and noninvasive surrogates thereof relying on computational fluid dynamics (CFD) have been developed. Here, we present a CFD-based workflow using digital twins of patients' carotid bifurcations derived from computed tomography angiography for the noninvasive functional assessment of CAD. We reconstructed patient-specific digital twins of 37 carotid bifurcations. We implemented a CFD model using common carotid artery peak systolic velocity (PSV) acquired with Doppler ultrasound (DUS) as inlet boundary condition and a two-element Windkessel model as oulet boundary condition. The agreement between CFD and DUS on the PSV in the internal carotid artery (ICA) was then compared. The relative error for the agreement between DUS and CFD was 9% ± 20% and the intraclass correlation coefficient was 0.88. Furthermore, hyperemic simulations in a physiological range were feasible and unmasked markedly different pressure drops along two ICA stenoses with similar degree of narrowing under comparable ICA blood flow. Hereby, we lay the foundation for prospective studies on noninvasive CFD-based derivation of metrics similar to the FFR for the assessment of CAD.

Queckenstedt's test repurposed for the quantitative assessment of the cerebrospinal fluid pulsatility curve

PURPOSE

Before the era of spinal imaging, presence of a spinal canal block was tested through gross changes in cerebrospinal fluid pressure (CSFP) provoked by manual compression of the jugular veins (referred to as Queckenstedt's test; QT). Beyond these provoked gross changes, cardiac-driven CSFP peak-to-valley amplitudes (CSFPp) can be recorded during CSFP registration. This is the first study to assess whether the QT can be repurposed to derive descriptors of the CSF pulsatility curve, focusing on feasibility and repeatability.

METHOD

Lumbar puncture was performed in lateral recumbent position in fourteen elderly patients (59.7±9.3 years, 6F) (NCT02170155) without stenosis of the spinal canal. CSFP was recorded during resting state and QT. A surrogate for the relative pulse pressure coefficient was computed from repeated QTs (i.e., RPPC-Q).

RESULTS

Resting state mean CSFP was 12.3 mmHg (IQR 3.2) and CSFPp was 1.0 mmHg (0.5). Mean CSFP rise during QT was 12.5 mmHg (7.3). CSFPp showed an average 3-fold increase at peak QT compared to the resting state. Median RPPC-Q was 0.18 (0.04). There was no systematic error in the computed metrics between the first and second QT.

CONCLUSION

This technical note describes a method to reliably derive, beyond gross CSFP increments, metrics related to cardiac-driven amplitudes during QT (i.e., RPPC-Q). A study comparing these metrics as obtained by established procedures (i.e., infusion testing) and by QT is warranted.

Large-scale morphometry of the subarachnoid space of the optic nerve

BACKGROUND

The meninges, formed by dura, arachnoid and pia mater, cover the central nervous system and provide important barrier functions. Located between arachnoid and pia mater, the cerebrospinal fluid (CSF)-filled subarachnoid space (SAS) features a variety of trabeculae, septae and pillars. Like the arachnoid and the pia mater, these structures are covered with leptomeningeal or meningothelial cells (MECs) that form a barrier between CSF and the parenchyma of the optic nerve (ON). MECs contribute to the CSF proteome through extensive protein secretion. In vitro, they were shown to phagocytose potentially toxic proteins, such as α-synuclein and amyloid beta, as well as apoptotic cell bodies. They therefore may contribute to CSF homeostasis in the SAS as a functional exchange surface. Determining the total area of the SAS covered by these cells that are in direct contact with CSF is thus important for estimating their potential contribution to CSF homeostasis.

METHODS

Using synchrotron radiation-based micro-computed tomography (SRµCT), two 0.75 mm-thick sections of a human optic nerve were acquired at a resolution of 0.325 µm/pixel, producing images of multiple terabytes capturing the geometrical details of the CSF space. Special-purpose supercomputing techniques were employed to obtain a pixel-accurate morphometric description of the trabeculae and estimate internal volume and surface area of the ON SAS.

RESULTS

In the bulbar segment, the ON SAS microstructure is shown to amplify the MECs surface area up to 4.85-fold compared to an "empty" ON SAS, while just occupying 35% of the volume. In the intraorbital segment, the microstructure occupies 35% of the volume and amplifies the ON SAS area 3.24-fold.

CONCLUSIONS

We provided for the first time an estimation of the interface area between CSF and MECs. This area is of importance for estimating a potential contribution of MECs on CSF homeostasis.

Cerebrospinal Fluid Pressure Dynamics as a Bedside Test in Traumatic Spinal Cord Injury to Assess Surgical Spinal Cord Decompression: Safety, Feasibility, and Proof-of-Concept

BACKGROUND

Sufficient and timely spinal cord decompression is a critical surgical objective for neurological recovery in spinal cord injury (SCI). Residual cord compression may be associated with disturbed cerebrospinal fluid pressure (CSFP) dynamics.

OBJECTIVES

This study aims to assess whether intrathecal CSFP dynamics in SCI following surgical decompression are feasible and safe, and to explore the diagnostic utility.

METHODS

Prospective cohort study. Bedside lumbar CSFP dynamics and cervical MRI were obtained following surgical decompression in N = 9 with mostly cervical acute-subacute SCI and N = 2 patients with non-traumatic SCI. CSFP measurements included mean CSFP, cardiac-driven CSFP peak-to-valley amplitudes (CSFPp), Valsalva maneuver, and Queckenstedt's test (firm pressure on jugular veins, QT). From QT, proxies for cerebrospinal fluid pulsatility curve were calculated (ie, relative pulse pressure coefficient; RPPC-Q). CSFP metrics were compared to spine-healthy patients. computer tomography (CT)-myelography was done in 3/8 simultaneous to CSFP measurements.

RESULTS

Mean age was 45 ± 9 years (range 17-67; 3F), SCI was complete (AIS A, N = 5) or incomplete (AIS B-D, N = 6). No adverse events related to CSFP assessments. CSFP rise during QT was induced in all patients [range 9.6-26.6 mmHg]. However, CSFPp was reduced in 3/11 (0.1-0.3 mmHg), and in 3/11 RPPC-Q was abnormal (0.01-0.05). Valsalva response was reduced in 8/11 (2.6-23.4 mmHg). CSFP dynamics corresponded to CT-myelography.

CONCLUSIONS

Comprehensive bedside lumbar CSFP dynamics in SCI following decompression are safe, feasible, and can reveal distinct patterns of residual spinal cord compression. Longitudinal studies are required to define critical thresholds of impaired CSFP dynamics that may impact neurological recovery and requiring surgical revisions.

The effect of body position change on noninvasively acquired intracranial pulse waves

Objective. Craniospinal compliance (CC) is an important metric for the characterization of space-occupying neurological pathologies. CC is obtained using invasive procedures that carry risks for the patients. Therefore, noninvasive methods for acquiring surrogates of CC have been proposed, most recently based on changes in the head's dielectric properties during the cardiac cycle. Here, we have tested whether changes in body position, which are known to influence CC, are reflected in a capacitively acquired signal (hereinafter referred to as W) originating from dynamic changes of the head's dielectric properties. &#xD;Approach. Eighteen young healthy volunteers were included in the study. After 10 minutes in supine position, subjects were tilted head-up (HUT), back to 0° (horizontal, control), and then head-down (HDT). Metrics related to cardiovascular action were extracted from W, including AMP, the peak-to-valley amplitude of the cardiac modulation of W. Computational electromagnetic simulations were performed to probe the association between intracranial volume change and W. &#xD;Main results. AMP decreased during HUT (0°: 2869±597 arbitrary units (au); +75°: 2307±490 au, P=0.002) and increased during HDT (-30°: 4403±1428 au, P<0.0001). The same behavior was predicted by the electromagnetic model. &#xD;Significance. Tilting affects the distribution of CC between cranial and spinal compartments. Cardiovascular action induces compliance-dependent oscillatory changes in the intracranial fluid composition, which causes corresponding variations in the head's dielectric properties. These manifest as increasing AMP with decreasing intracranial compliance, which suggests that W may contain information related to CC, and that it might be possible to derive CC surrogates therefrom.

Are standing osmotic gradients the main driver of cerebrospinal fluid production? A computational analysis

BACKGROUND

The mechanisms of cerebrospinal fluid (CSF) production by the ventricular choroid plexus (ChP) have not been fully deciphered. One prominent hypothesized mechanism is trans-epithelial water transport mediated by accumulation of solutes at the luminal ChP membrane that produces local osmotic gradients. However, this standing osmotic gradient hypothesis has not been systematically tested.

METHODS

To assess the plausibility of the standing gradient mechanism serving as the main driver of CSF production by the ChP, we developed a three-dimensional (3D) and a one-dimensional (1D) computational model to quantitatively describe the associated processes in the rat ChP inter-microvillar spaces and in CSF pools between macroscopic ChP folds (1D only). The computationally expensive 3D model was used to examine the applicability of the 1D model for hypothesis testing. The 1D model was employed to predict the rate of CSF produced by the standing gradient mechanism for 200,000 parameter permutations. Model parameter values for each permutation were chosen by random sampling from distributions derived from published experimental data.

RESULTS

Both models predict that the CSF production rate by the standing osmotic gradient mechanism is below 10% of experimentally measured values that reflect the contribution of all actual production mechanisms. The 1D model indicates that increasing the size of CSF pools between ChP folds, where diffusion dominates solute transport, would increase the contribution of the standing gradient mechanism to CSF production.

CONCLUSIONS

The models suggest that the effect of standing osmotic gradients is too small to contribute substantially to CSF production. ChP motion and movement of CSF in the ventricles, which are not accounted for in the models, would further reduce this effect, making it unlikely that standing osmotic gradients are the main drivers of CSF production.

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.

Theory for a non-invasive diagnostic biomarker for craniospinal diseases

Monitoring intracranial pressure (ICP) and craniospinal compliance (CC) is frequently required in the treatment of patients suffering from craniospinal diseases. However, current approaches are invasive and cannot provide continuous monitoring of CC. Dynamic exchange of blood and cerebrospinal fluid (CSF) between cranial and spinal compartments due to cardiac action transiently modulates the geometry and dielectric properties of the brain. The resulting impedance changes can be measured and might be usable as a non-invasive CC surrogate. A numerically robust and computationally efficient approach based on the reciprocity theorem was developed to compute dynamic impedance changes resulting from small geometry and material property changes. The approach was successfully verified against semi-analytical benchmarks, before being combined with experimental brain pulsation data to study the information content of the impedance variation. The results indicate that the measurable signal is dominated by the pulsatile displacement of the cortical brain surface, with minor contributions from the ventricular surfaces and from changes in brain perfusion. Different electrode setups result in complementary information. The information content from the investigated three electrode pairs was employed to successfully infer subject-specific brain pulsation and motion features. This suggests that non-invasive CC surrogates based on impedance monitoring could be established.

Membrane transporters control cerebrospinal fluid formation independently of conventional osmosis to modulate intracranial pressure

Background

Disturbances in the brain fluid balance can lead to life-threatening elevation in the intracranial pressure (ICP), which represents a vast clinical challenge. Nevertheless, the details underlying the molecular mechanisms governing cerebrospinal fluid (CSF) secretion are largely unresolved, thus preventing targeted and efficient pharmaceutical therapy of cerebral pathologies involving elevated ICP.

Methods

Experimental rats were employed for in vivo determinations of CSF secretion rates, ICP, blood pressure and ex vivo excised choroid plexus for morphological analysis and quantification of expression and activity of various transport proteins. CSF and blood extractions from rats, pigs, and humans were employed for osmolality determinations and a mathematical model employed to determine a contribution from potential local gradients at the surface of choroid plexus.

Results

We demonstrate that CSF secretion can occur independently of conventional osmosis and that local osmotic gradients do not suffice to support CSF secretion. Instead, the CSF secretion across the luminal membrane of choroid plexus relies approximately equally on the Na⁺/K⁺/2Cl⁻ cotransporter NKCC1, the Na⁺/HCO₃⁻ cotransporter NBCe2, and the Na⁺/K⁺-ATPase, but not on the Na⁺/H⁺ exchanger NHE1. We demonstrate that pharmacological modulation of CSF secretion directly affects the ICP.

Conclusions

CSF secretion appears to not rely on conventional osmosis, but rather occur by a concerted effort of different choroidal transporters, possibly via a molecular mode of water transport inherent in the proteins themselves. Therapeutic modulation of the rate of CSF secretion may be employed as a strategy to modulate ICP. These insights identify new promising therapeutic targets against brain pathologies associated with elevated ICP.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12987-022-00358-4.

Cerebrospinal fluid pressure dynamics reveal signs of effective spinal canal narrowing in ambiguous spine conditions

Spinal canal narrowing with consecutive spinal cord compression is considered a key mechanism in degenerative cervical myelopathy (DCM). DCM is a common spine condition associated with progressive neurological disability, and timely decompressive surgery is recommended. However, the clinical and radiological diagnostic workup is often ambiguous, challenging confident proactive treatment recommendations. Cerebrospinal fluid pressure dynamics (CSFP) are altered by spinal canal narrowing. Therefore, we aim to explore the potential value of bedside CSFP assessments for qualitative and quantitative assessment of spinal canal narrowing in DCM. In this prospective case series, seven patients with DCM underwent bedside lumbar puncture with measurement of CSFP dynamics and routine CSF analysis (NCT02170155). The patients were enrolled when standard diagnostic algorithms did not permit a clear treatment decision. Measurements include baseline CSFP, cardiac-driven CSFP peak-to-trough amplitude (CSFPp), and the Queckenstedt's test (firm pressure on jugular veins) in neutral and reclined head position. From the Queckenstedt's test, proxies for craniospinal elastance (i.e., relative pulse pressure coefficient; RPPC-Q) were calculated analogously to infusion testing. CSFP metrics were deemed suspicious of canal narrowing when numbers were lower than the minimum value from a previously tested elderly spine-healthy cohort (N = 14). Mean age was 56 ± 13 years (range, 38–75; 2F); symptom severity was mostly mild to moderate (mean mJOA, 13.5 ± 2.6; range, 9–17). All the patients showed some extent of cervical stenosis in the MRI of unclear significance (5/7 following decompressive cervical spine surgery with an adjacent level or residual stenosis). Baseline CSFP was normal except for one patient (range, 4.7–17.4 mmHg). Normal values were found for CSFPp (0.4–1.3 mmHg) and the Queckenstedt's test in normal head positioning (9.-25.3 mmHg). During reclination, the Queckenstedt's test significantly decreased in one, and CSFPp in another case (>50% compared to normal position). RPPC-Q (0.07–0.19) aligned with lower values from spine-healthy (0.10–0.44). Routine CSF examinations showed mild total protein elevation (mean, 522 ± 108 mg/ml) without further evidence for the disturbed blood brain barrier. Intrathecal CSFP measurements allow discerning disturbed from normal CSFP dynamics in this population. Prospective longitudinal studies should further evaluate the diagnostic utility of CSFP assessments in DCM.

Noninvasive Monitoring of Intracranial Pulse Waves

OBJECTIVE

The clinical management of several neurological disorders benefits from the assessment of intracranial pressure and craniospinal compliance. However, the associated procedures are invasive in nature. Here, we aimed to assess whether naturally occurring periodic changes in the dielectric properties of the head could serve as the basis for deriving surrogates of craniospinal compliance noninvasively.

METHODS

We designed a device and electrodes for noninvasive measurement of periodic changes of the dielectric properties of the human head. We characterized the properties of the device-electrode-head system by measurements on healthy volunteers, by computational modeling, and by electromechanical modeling. We then performed hyperventilation testing to assess whether the measured signal is of intracranial origin.

RESULTS

Signals obtained with the device on volunteers showed characteristic cardiac and respiratory modulations. Signal oscillations can be attributed primarily to changes in resistive properties of the head during cardiac and respiratory cycles. Reduction of end-tidal CO₂, through hyperventilation, resulted in a decrease in the signal amplitude associated with cardiovascular action.

CONCLUSION

Given the higher CO₂ reactivity of intracranial vessels compared to extracranial ones, the results of hyperventilation testing suggest that the acquired signal is, in part, of intracranial origin.

SIGNIFICANCE

If confirmed in larger cohorts, our observations suggest that noninvasive capacitive acquisition of changes in the dielectric properties of the head could be used to derive surrogates of craniospinal compliance.

Shape Trumps Size: Image-Based Morphological Analysis Reveals That the 3D Shape Discriminates Intracranial Aneurysm Disease Status Better Than Aneurysm Size

Background

To date, it remains difficult for clinicians to reliably assess the disease status of intracranial aneurysms. As an aneurysm's 3D shape is strongly dependent on the underlying formation processes, it is believed that the presence of certain shape features mirrors the disease status of the aneurysm wall. Currently, clinicians associate irregular shape with wall instability. However, no consensus exists about which shape features reliably predict instability. In this study, we present a benchmark to identify shape features providing the highest predictive power for aneurysm rupture status.

Methods

3D models of aneurysms were extracted from medical imaging data (3D rotational angiographies) using a standardized protocol. For these aneurysm models, we calculated a set of metrics characterizing the 3D shape: Geometry indices (such as undulation, ellipticity and non-sphericity); writhe- and curvature-based metrics; as well as indices based on Zernike moments. Using statistical learning methods, we investigated the association between shape features and aneurysm disease status. This processing was applied to a clinical dataset of 750 aneurysms (261 ruptured, 474 unruptured) registered in the AneuX morphology database. We report here statistical performance metrics [including the area under curve (AUC)] for morphometric models to discriminate between ruptured and unruptured aneurysms.

Results

The non-sphericity index NSI (AUC = 0.80), normalized Zernike energies ZNsurf (AUC = 0.80) and the modified writhe-index W¯meanL1 (AUC = 0.78) exhibited the strongest association with rupture status. The combination of predictors further improved the predictive performance (without location: AUC = 0.82, with location AUC = 0.87). The anatomical location was a good predictor for rupture status on its own (AUC = 0.78). Different protocols to isolate the aneurysm dome did not affect the prediction performance. We identified problems regarding generalizability if trained models are applied to datasets with different selection biases.

Conclusions

Morphology provided a clear indication of the aneurysm disease status, with parameters measuring shape (especially irregularity) being better predictors than size. Quantitative measurement of shape, alone or in conjunction with information about aneurysm location, has the potential to improve the clinical assessment of intracranial aneurysms.

Renal blood flow and oxygenation

Our kidneys receive about one-fifth of the cardiac output at rest and have a low oxygen extraction ratio, but may sustain, under some conditions, hypoxic injuries that might lead to chronic kidney disease. This is due to large regional variations in renal blood flow and oxygenation, which are the prerequisite for some and the consequence of other kidney functions. The concurrent operation of these functions is reliant on a multitude of neuro-hormonal signaling cascades and feedback loops that also include the regulation of renal blood flow and tissue oxygenation. Starting with open questions on regulatory processes and disease mechanisms, we review herein the literature on renal blood flow and oxygenation. We assess the current understanding of renal blood flow regulation, reasons for disparities in oxygen delivery and consumption, and the consequences of disbalance between O₂ delivery, consumption, and removal. We further consider methods for measuring and computing blood velocity, flow rate, oxygen partial pressure, and related parameters and point out how limitations of these methods constitute important hurdles in this area of research. We conclude that to obtain an integrated understanding of the relation between renal function and renal blood flow and oxygenation, combined experimental and computational modeling studies will be needed.

Insights Into the Low Rate of In-Pump Thrombosis With the HeartMate 3: Does the Artificial Pulse Improve Washout?

While earlier studies reported no relevant effect of the HeartMate 3 (HM3) artificial pulse (AP) on bulk pump washout, its effect on regions with prolonged residence times remains unexplored. Using numerical simulations, we compared pump washout in the HM3 with and without AP with a focus on the clearance of the last 5% of the pump volume. Results were examined in terms of flush-volume (V_f, number of times the pump was flushed with new blood) to probe the effect of the AP independent of changing flow rate. Irrespective of the flow condition, the HM3 washout scaled linearly with flush volume up to 70% washout and slowed down for the last 30%. Flush volumes needed to washout 95% of the pump were comparable with and without the AP (1.3–1.4 V_f), while 99% washout required 2.1–2.2 V_f with the AP vs. 2.5 V_f without the AP. The AP enhanced washout of the bend relief and near-wall regions. It also transiently shifted or eliminated stagnation regions and led to rapid wall shear stress fluctuations below the rotor and in the secondary flow path. Our results suggest potential benefits of the AP for clearance of fluid regions that might elicit in-pump thrombosis and provide possible mechanistic rationale behind clinical data showing very low rate of in-pump thrombosis with the HM3. Further optimization of the AP sequence is warranted to balance washout efficacy while limiting blood damage.

Analysis of L-leucine amino acid transporter species activity and gene expression by human blood brain barrier hCMEC/D3 model reveal potential LAT1, LAT4, B0AT2 and y+LAT1 functional cooperation

In the CNS, amino acid (AA) neurotransmitters and neurotransmitter precursors are subject to tight homeostatic control mediated by blood-brain barrier (BBB) solute carrier amino acid transporters (AATs). Since the BBB is composed of multiple closely apposed cell types and opportunities for human in vivo studies are limited, we used in vitro and computational approaches to investigate human BBB AAT activity and regulation. Quantitative real-time PCR (qPCR) of the human BBB endothelial cell model hCMEC/D3 (D3) was used to determine expression of selected AAT, tight junction (TJ), and signal transduction (ST) genes under various culture conditions. L-leucine uptake data were interrogated with a computational model developed by our group for calculating AAT activity in complex cell cultures. This approach is potentially applicable to in vitro cell culture drug studies where multiple "receptors" may mediate observed responses. Of 7 Leu AAT genes expressed by D3 only the activity of SLC7A5-SLC3A2/LAT1-4F2HC (LAT1), SLC43A2/LAT4 (LAT4) and sodium-dependent AATs, SLC6A15/B⁰AT2 (B⁰AT2), and SLC7A7/y⁺LAT1 (y⁺LAT1) were calculated to be required for Leu uptake. Therefore, D3 Leu transport may be mediated by a potentially physiologically relevant functional cooperation between the known BBB AAT, LAT1 and obligatory exchange (y⁺LAT1), facilitative diffusion (LAT4), and sodium symporter (B⁰AT2) transporters.

Virtual histology of an entire mouse brain from formalin fixation to paraffin embedding. Part 2: Volumetric strain fields and local contrast changes

BACKGROUND

Fixation and embedding of post mortem brain tissue is a pre-requisite for both gold-standard conventional histology and X-ray virtual histology. This process alters the morphology and density of the brain microanatomy.

NEW METHOD

To quantify these changes, we employed synchrotron radiation-based hard X-ray tomography with 3 μm voxel length to visualize the same mouse brain after fixation in 4% formalin, immersion in ethanol solutions (50%, 70%, 80%, 90%, and 100%), xylene, and finally after embedding in a paraffin block. The volumetric data were non-rigidly registered to the initial formalin-fixed state to align the microanatomy within the entire mouse brain.

RESULTS

Volumetric strain fields were used to characterize local shrinkage, which was found to depend on the anatomical region and distance to external surface. X-ray contrast was altered and enhanced by preparation-induced inter-tissue density changes. The preparation step can be selected to highlight specific anatomical features. For example, fiber tract contrast is amplified in 100% ethanol.

COMPARISON WITH EXISTING METHODS

Our method provides volumetric strain fields, unlike approaches based on feature-to-feature or volume measurements. Volumetric strain fields are produced by non-rigid registration, which is less labor-intensive and observer-dependent than volume change measurements based on manual segmentations. X-ray microtomography provides spatial resolution at least an order of magnitude higher than magnetic resonance microscopy, allowing for analysis of morphology and density changes within the brain's microanatomy.

CONCLUSION

Our approach belongs to three-dimensional virtual histology with isotropic micrometer spatial resolution and therefore complements atlases based on a combination of magnetic resonance microscopy and optical micrographs of serial histological sections.

Intraoperative monitoring of CSF pressure in patients with degenerative cervical myelopathy (COMP-CORD Study): a prospective cohort study

Degenerative cervical myelopathy (DCM) is hallmarked by spinal canal narrowing and related cord compression and myelopathy. CSF pressure dynamics are likely disturbed due to spinal canal stenosis. The study aims to investigate the diagnostic value of continuous intraoperative CSF pressure monitoring during surgical decompression. Prospective single center study (NCT02170155) with enrolment of DCM patients that underwent surgical decompression between December 2019 and May 2021. Data from N=17 patients were analyzed, symptom severity graded with the modified Japanese Orthopedic Score (mJOA). CSF pulsations were continuously monitored with a lumbar intrathecal catheter during surgical decompression. Mean patient age was 62±9 years (range 38-73; 8F), symptoms were mild-moderate in most patients (mean mJOA 14±2, range 10-18). Measurements were well tolerated without safety concerns. In 15/16 (94%) CSF pulsations increased at the time of surgical decompression. In one case, responsiveness could not be evaluated for technical reasons. Unexpected CSF pulsation decrease was related to adverse events (i.e., CSF leakage). Median CSF pulsation amplitudes increased from pre-decompression (0.52 mmHg [IQR 0.71]) to post-decompression (0.72 mmHg [IQR 0.96]) (P=0.001). Mean baseline CSF pressure increased with lower magnitude than pulsations, from 9.5±3.5 to 10.3±3.8 mmHg (P=0.003). Systematic relations of CSF pulsations were confined to surgical decompression, independent of arterial blood pressure (P=0.927) or heart rate (P=0.102). Intraoperative CSF pulsation monitoring was sensitive, timely, and specifically related to surgical decompression while in addition adverse events could be discerned. Further investigation of the clinical value of intraoperative guidance for decompression in complex DCM surgery is promising.

Newborn Incubators Do Not Protect from High Noise Levels in the Neonatal Intensive Care Unit and Are Relevant Noise Sources by Themselves

BACKGROUND

While meaningful sound exposure has been shown to be important for newborn development, an excess of noise can delay the proper development of the auditory cortex.

AIM

The aim of this study was to assess the acoustic environment of a preterm baby in an incubator on a newborn intensive care unit (NICU).

METHODS

An empty but running incubator (Giraffe Omnibed, GE Healthcare) was used to evaluate the incubator frequency response with 60 measurements. In addition, a full day and night period outside and inside the incubator at the NICU of the University Hospital Zurich was acoustically analyzed.

RESULTS

The fan construction inside the incubator generates noise in the frequency range of 1.3-1.5 kHz with a weighted sound pressure level (SPL) of 40.5 dB(A). The construction of the incubator narrows the transmitted frequency spectrum of sound entering the incubator to lower frequencies, but it does not attenuate transient noises such as alarms or opening and closing of cabinet doors substantially. Alarms, as generated by the monitors, the incubator, and additional devices, still pass to the newborn.

CONCLUSIONS

The incubator does protect only insufficiently from noise coming from the NICUThe transmitted frequency spectrum is changed, limiting the impact of NICU noise on the neonate, but also limiting the neonate's perception of voices. The incubator, in particular its fan, as well as alarms from patient monitors are major sources of noise. Further optimizations with regard to the sound exposure in the NICU, as well as studies on the role of the incubator as a source and modulator, are needed to meet the preterm infants' multi-sensory needs.

In Vitro Testing and Comparison of Additively Manufactured Polymer Impellers for the CentriMag Blood Pump

Additive manufacturing (AM) is an effective tool for accelerating knowledge gain in development processes, as it enables the production of complex prototypes at low cost and with short lead times. In the development of mechanical circulatory support, the use of cheap polymer-based AM techniques for prototype manufacturing allows more design variations to be tested, promoting a better understanding of the respective system and its optimization parameters. Here, we compare four commonly used AM processes for polymers with respect to manufacturing accuracy, surface roughness, and shape fidelity in an aqueous environment. Impeller replicas of the CentriMag blood pump were manufactured with each process and integrated into original pump housings. The assemblies were tested for hydraulic properties and hemolysis in reference to the commercially available pump. Computational fluid dynamic simulations were carried out to support the transfer of the results to other applications. In hydraulic testing, the deviation in pressure head and motor current of all additively manufactured replicas from the reference pump remained below 2% over the entire operating range of the pump. In contrast, significant deviations of up to 620% were observed in hemolysis testing. Only the replicas produced by stereolithography showed a nonsignificant deviation from the reference pump, which we attribute to the low surface roughness of parts manufactured thereby. The results suggest that there is a flow-dependent threshold of roughness above which a surface strongly contributes to cell lysis by promoting a hydraulically rough boundary flow.

High-Throughput Lossy-to-Lossless 3D Image Compression

The rapid increase in medical and biomedical image acquisition rates has opened up new avenues for image analysis, but has also introduced formidable challenges. This is evident, for example, in selective plane illumination microscopy where acquisition rates of about 1-4 GB/s sustained over several days have redefined the scale of I/O bandwidth required by image analysis tools. Although the effective bandwidth could, principally, be increased by lossy-to-lossless data compression, this is of limited value in practice due to the high computational demand of current schemes such as JPEG2000 that reach compression throughput of one order of magnitude below that of image acquisition. Here we present a novel lossy-to-lossless data compression scheme with a compression throughput well above 4 GB/s and compression rates and rate-distortion curves competitive with those achieved by JPEG2000 and JP3D.

RAQ: a novel surrogate for the craniospinal pressure-volume relationship

OBJECTIVE

The intracranial pressure-volume relation contains information relevant for diagnostics of hydrocephalus and other space-occupying pathologies. We aimed to design a noise-resilient surrogate for this relationship that can be calculated from intracranial pressure (ICP) signals.

APPROACH

The new surrogate, termed respiratory amplitude quotient (RAQ), characterizes the modulation of the cardiac pulse wave amplitude by the respiratory wave in the ICP time course. RAQ is defined as the ratio of the amplitude of the respiratory wave in the ICP signal to the amplitude of the respiration-induced variation in the course of the cardiac pulse wave amplitude. We validated the calculation of RAQ on synthetically generated ICP waveforms. We further extracted RAQ retrospectively from overnight ICP recordings in a cohort of hydrocephalus patients with aqueductal stenosis, age 55.8 ± 18.0 years, and a comparison group with hydrocephalus diagnosed by morphology in MRI, but not responsive to either external lumbar drainage or ventriculo-peritoneal shunting, age 72.5 ± 6.1 years. RAQ was determined for the full recordings, and separately for periods containing B-waves.

MAIN RESULTS

We found a mean difference of less than 2% between the calculated values of RAQ and the theoretically determined equivalent descriptors of the synthetic ICP waveforms. In the overnight recordings, we found significantly different RAQ values during B-waves in the aqueductal stenosis (0.86 ± 0.11) and non-responsive hydrocephalus patient groups (1.07 ± 0.20), p = 0.027. In contrast, there was no significant difference in other tested parameters, namely pressure-volume index, elastance coefficient, and resistance to outflow. Neither did we find significant difference when considering RAQ over the full recordings.

SIGNIFICANCE

Our results indicate that RAQ may function as a potential surrogate for the intracranial pressure-volume relation.

Shape irregularity of the intracranial aneurysm lumen exhibits diagnostic value

BACKGROUND

Morphological irregularity is linked to intracranial aneurysm wall instability and manifests in the lumen shape. Yet there is currently no consent on how to assess shape irregularity. The aims of this work are to quantify irregularity as perceived by clinicians, to break down irregularity into morphological attributes, and to relate these to clinically relevant factors such as rupture status, aneurysm location, and patient age or sex.

METHODS

Thirteen clinicians and 26 laypersons assessed 134 aneurysm lumen segmentations in terms of overall perceived irregularity and five different morphological attributes (presence/absence of a rough surface, blebs, lobules, asymmetry, complex geometry of the parent vasculature). We examined rater agreement and compared the ratings with clinical factors by means of regression analysis or binary classification.

RESULTS

Using rank-based aggregation, the irregularity ratings of clinicians and laypersons did not differ statistically. Perceived irregularity showed good agreement with curvature (coefficient of determination R² = 0.68 ± 0.08) and was modeled very accurately using the five morphological rating attributes plus shape elongation (R² = 0.95 ± 0.02). In agreement with previous studies, irregularity was associated with aneurysm rupture status (AUC = 0.81 ± 0.08); adding aneurysm location as an explanatory variable increased the AUC to 0.87 ± 0.09. Besides irregularity, perceived asymmetry, presence of blebs or lobules, aneurysm size, non-sphericity, and curvature were linked to rupture. No association was found between morphology and any of patient sex, age, and history of smoking or hypertension. Aneurysm size was linked to morphology.

CONCLUSIONS

Irregular lumen shape carries significant information on the aneurysm's disease status. Irregularity constitutes a continuous parameter that shows a strong association with the rupture status. To improve the objectivity of morphological assessment, we suggest examining shape through six different morphological attributes, which can characterize irregularity accurately.

Extending statistical learning for aneurysm rupture assessment to Finnish and Japanese populations using morphology, hemodynamics, and patient characteristics

OBJECTIVE

Incidental aneurysms pose a challenge for physicians, who need to weigh the rupture risk against the risks associated with treatment and its complications. A statistical model could potentially support such treatment decisions. A recently developed aneurysm rupture probability model performed well in the US data used for model training and in data from two European cohorts for external validation. Because Japanese and Finnish patients are known to have a higher aneurysm rupture risk, the authors' goals in the present study were to evaluate this model using data from Japanese and Finnish patients and to compare it with new models trained with Finnish and Japanese data.

METHODS

Patient and image data on 2129 aneurysms in 1472 patients were used. Of these aneurysm cases, 1631 had been collected mainly from US hospitals, 249 from European (other than Finnish) hospitals, 147 from Japanese hospitals, and 102 from Finnish hospitals. Computational fluid dynamics simulations and shape analyses were conducted to quantitatively characterize each aneurysm's shape and hemodynamics. Next, the previously developed model's discrimination was evaluated using the Finnish and Japanese data in terms of the area under the receiver operating characteristic curve (AUC). Models with and without interaction terms between patient population and aneurysm characteristics were trained and evaluated including data from all four cohorts obtained by repeatedly randomly splitting the data into training and test data.

RESULTS

The US model's AUC was reduced to 0.70 and 0.72, respectively, in the Finnish and Japanese data compared to 0.82 and 0.86 in the European and US data. When training the model with Japanese and Finnish data, the average AUC increased only slightly for the Finnish sample (to 0.76 ± 0.16) and Finnish and Japanese cases combined (from 0.74 to 0.75 ± 0.14) and decreased for the Japanese data (to 0.66 ± 0.33). In models including interaction terms, the AUC in the Finnish and Japanese data combined increased significantly to 0.83 ± 0.10.

CONCLUSIONS

Developing an aneurysm rupture prediction model that applies to Japanese and Finnish aneurysms requires including data from these two cohorts for model training, as well as interaction terms between patient population and the other variables in the model. When including this information, the performance of such a model with Japanese and Finnish data is close to its performance with US or European data. These results suggest that population-specific differences determine how hemodynamics and shape associate with rupture risk in intracranial aneurysms.

Simultaneous Three-Dimensional Vascular and Tubular Imaging of Whole Mouse Kidneys With X-ray μCT

Concurrent three-dimensional imaging of the renal vascular and tubular systems on the whole-kidney scale with capillary level resolution is labor-intensive and technically difficult. Approaches based on vascular corrosion casting and X-ray micro computed tomography (μCT), for example, suffer from vascular filling artifacts and necessitate imaging with an additional modality to acquire tubules. In this work, we report on a new sample preparation, image acquisition, and quantification protocol for simultaneous vascular and tubular μCT imaging of whole, uncorroded mouse kidneys. The protocol consists of vascular perfusion with the water-soluble, aldehyde-fixable, polymeric X-ray contrast agent XlinCA, followed by laboratory-source μCT imaging and structural analysis using the freely available Fiji/ImageJ software. We achieved consistent filling of the entire capillary bed and staining of the tubules in the cortex and outer medulla. After imaging at isotropic voxel sizes of 3.3 and 4.4 μm, we segmented vascular and tubular systems and quantified luminal volumes, surface areas, diffusion distances, and vessel path lengths. This protocol permits the analysis of vascular and tubular parameters with higher reliability than vascular corrosion casting, less labor than serial sectioning and leaves tissue intact for subsequent histological examination with light and electron microscopy.

Crosslinkable polymeric contrast agent for high-resolution X-ray imaging of the vascular system

A contrast agent for X-ray micro computed tomography (μCT), called XlinCA, that combines reliable perfusion and permanent retention and contrast properties, was developed for ex vivo imaging.

Integrated Flow Chamber System for Live Cell Microscopy

In vitro quantification of the effect of mechanical loads on cells by live microscopy requires precise control of load and culture environment. Corresponding systems are often bulky, their setup and maintenance are time consuming, or the cell yield is low. Here, we show the design and initial testing of a new cell culture system that fits on standard light microscope stages. Based on the parallel plate principle, the system allows for live microscopy of cells exposed to flow-induced shear stress, features short setup time and requires little user interaction. An integrated feedback-controlled heater and a bubble trap enable long observation times. The key design feature is the possibility for quick exchange of the cultured cells. We present first test results that focus on verifying the robustness, biocompatibility, and ease of use of the device.

Propagation of Plasma L-Phenylalanine Concentration Fluctuations to the Neurovascular Unit in Phenylketonuria: An in silico Study

Phenylketonuria (PKU) is an inherited metabolic disease characterized by abnormally high concentrations of the essential amino acid L-phenylalanine (Phe) in blood plasma caused by reduced activity of phenylalanine hydroxylase (PAH). While numerous studies have shown association between high plasma Phe concentration and intellectual impairment, it is not clear whether increased Phe fluctuations also observed in PKU affect the brain as well. To investigate this, time-resolved in vivo data on Phe and competing large neutral amino acid (LNAA) concentrations in neurons are needed, but cannot be acquired readily with current methods. We have used in silico modeling as an alternative approach to characterize the interactive dynamics of Phe and competing LNAAs (CL) in the neurovascular unit (NVU). Our results suggest that plasma Phe fluctuations can propagate into the NVU cells and change there the concentration of LNAAs, with the highest magnitude of this effect observed at low frequency and high amplitude-to-mean ratio of the plasma Phe concentration fluctuations. Our model further elucidates the effect of therapeutic LNAA supplementation in PKU, showing how abnormal concentrations of Phe and CL in the NVU move thereby toward normal physiologic levels.

Fluid Dynamics in the HeartMate 3: Influence of the Artificial Pulse Feature and Residual Cardiac Pulsation

Ventricular assist devices (VADs), among which the HeartMate 3 (HM3) is the latest clinically approved representative, are often the therapy of choice for patients with end-stage heart failure. Despite advances in the prevention of pump thrombosis, rates of stroke and bleeding remain high. These complications are attributed to the flow field within the VAD, among other factors. One of the HM3's characteristic features is an artificial pulse that changes the rotor speed periodically by 4000 rpm, which is meant to reduce zones of recirculation and stasis. In this study, we investigated the effect of this speed modulation on the flow fields and stresses using high-resolution computational fluid dynamics. To this end, we compared Eulerian and Lagrangian features of the flow fields during constant pump operation, during operation with the artificial pulse feature, and with the effect of the residual native cardiac cycle. We observed good washout in all investigated situations, which may explain the low incidence rates of pump thrombosis. The artificial pulse had no additional benefit on scalar washout performance, but it induced rapid variations in the flow velocity and its gradients. This may be relevant for the removal of deposits in the pump. Overall, we found that viscous stresses in the HM3 were lower than in other current VADs. However, the artificial pulse substantially increased turbulence, and thereby also total stresses, which may contribute to clinically observed issues related to hemocompatibility.

Functional Polarity of Microvascular Brain Endothelial Cells Supported by Neurovascular Unit Computational Model of Large Neutral Amino Acid Homeostasis

The homeostatic regulation of large neutral amino acid (LNAA) concentration in the brain interstitial fluid (ISF) is essential for proper brain function. LNAA passage into the brain is primarily mediated by the complex and dynamic interactions between various solute carrier (SLC) transporters expressed in the neurovascular unit (NVU), among which SLC7A5/LAT1 is considered to be the major contributor in microvascular brain endothelial cells (MBEC). The LAT1-mediated trans-endothelial transport of LNAAs, however, could not be characterized precisely by available in vitro and in vivo standard methods so far. To circumvent these limitations, we have incorporated published in vivo data of rat brain into a robust computational model of NVU-LNAA homeostasis, allowing us to evaluate hypotheses concerning LAT1-mediated trans-endothelial transport of LNAAs across the blood brain barrier (BBB). We show that accounting for functional polarity of MBECs with either asymmetric LAT1 distribution between membranes and/or intrinsic LAT1 asymmetry with low intraendothelial binding affinity is required to reproduce the experimentally measured brain ISF response to intraperitoneal (IP) L-tyrosine and L-phenylalanine injection. On the basis of these findings, we have also investigated the effect of IP administrated L-tyrosine and L-phenylalanine on the dynamics of LNAAs in MBECs, astrocytes and neurons. Finally, the computational model was shown to explain the trans-stimulation of LNAA uptake across the BBB observed upon ISF perfusion with a competitive LAT1 inhibitor.

Significant Association of Slow Vasogenic ICP Waves with Normal Pressure Hydrocephalus Diagnosis

OBJECTIVE

We aimed to test whether there is an association of slow vasogenic wave (SVW) occurrence with positive response to external lumbar drainage (ELD) and ventriculoperitoneal shunting and to design a method for the recognition and quantification of SVWs in the intracranial pressure (ICP) signal.

MATERIALS AND METHODS

We constructed SVW templates using normalized sine waves. We calculated the cross-correlation between the respective SVW template and the ICP signal. This was followed by shifting the templates forward and performing the cross-correlation analysis again until the end of the recording. Cross-correlation values above a threshold were considered to be indicative of SVWs. This threshold was previously determined and validated on a sample of ICP records of six patients. We calculated the root mean square of the recognized SVW periods as a measure of signal strength. Time-averaged signal strength was calculated over the full recording time (ICP_Smean) and over the wave periods (ICP_S).

RESULTS

We determined ICP_S and ICP_Smean in recordings of 2 groups of patients presenting with Hakim's triad: 26 normal pressure hydrocephalus (NPH) patients and 20 non-NPH patients. We then tested whether there was an association between ICP_S or ICP_Smean and the respective diagnosis using a Mann-Whitney test. We found significant association between ICP_S (p = 0.014) and ICP_Smean (p = 0.022) and the diagnoses.

CONCLUSIONS

The described method based on pattern recognition in the time domain is suitable for the detection and quantification of SVWs in ICP signals. We found a significant association between the occurrence of SVWs and independent NPH diagnosis.

Virtual surgical planning, flow simulation, and 3-dimensional electrospinning of patient-specific grafts to optimize Fontan hemodynamics

BACKGROUND

Despite advances in the Fontan procedure, there is an unmet clinical need for patient-specific graft designs that are optimized for variations in patient anatomy. The objective of this study is to design and produce patient-specific Fontan geometries, with the goal of improving hepatic flow distribution (HFD) and reducing power loss (P_loss), and manufacturing these designs by electrospinning.

METHODS

Cardiac magnetic resonance imaging data from patients who previously underwent a Fontan procedure (n = 2) was used to create 3-dimensional models of their native Fontan geometry using standard image segmentation and geometry reconstruction software. For each patient, alternative designs were explored in silico, including tube-shaped and bifurcated conduits, and their performance in terms of P_loss and HFD probed by computational fluid dynamic (CFD) simulations. The best-performing options were then fabricated using electrospinning.

RESULTS

CFD simulations showed that the bifurcated conduit improved HFD between the left and right pulmonary arteries, whereas both types of conduits reduced P_loss. In vitro testing with a flow-loop chamber supported the CFD results. The proposed designs were then successfully electrospun into tissue-engineered vascular grafts.

CONCLUSIONS

Our unique virtual cardiac surgery approach has the potential to improve the quality of surgery by manufacturing patient-specific designs before surgery, that are also optimized with balanced HFD and minimal P_loss, based on refinement of commercially available options for image segmentation, computer-aided design, and flow simulations.

Influence of Standard Laboratory Procedures on Measures of Erythrocyte Damage

The ability to characterize the mechanical properties of erythrocytes is important in clinical and research contexts: to diagnose and monitor hematologic disorders, as well as to optimize the design of cardiovascular implants and blood circulating devices with respect to blood damage. However, investigation of red blood cell (RBC) properties generally involves preparatory and processing steps. Even though these impose mechanical stresses on cells, little is known about their impact on the final measurement results. In this study, we investigated the effect of centrifuging, vortexing, pipetting, and high pressures on several markers of mechanical blood damage and RBC membrane properties. Using human venous blood, we analyzed erythrocyte damage by measuring free hemoglobin, phosphatidylserine exposure by flow cytometry, RBC deformability by ektacytometry and the parameters of a complete blood count. We observed increased levels of free hemoglobin for all tested procedures. The release of hemoglobin into plasma depended significantly on the level of stress. Elevated pressures and centrifuging also altered mean cell volume (MCV) and mean corpuscular hemoglobin (MCH), suggesting changes in erythrocyte population, and membrane properties. Our results show that the effects of blood handling can significantly influence erythrocyte damage metrics. Careful quantification of this influence as well as other unwanted secondary effects should thus be included in experimental protocols and accounted for in clinical laboratories.

Comparison of anti-siphon devices-how do they affect CSF dynamics in supine and upright posture?

BACKGROUND

Three different types of anti-siphon devices (ASDs) have been developed to counteract siphoning-induced overdrainage in upright posture. However, it is not known how the different ASDs affect CSF dynamics under the complex pressure environment seen in clinic due to postural changes. We investigated which ASDs can avoid overdrainage in upright posture best without leading to CSF accumulation.

METHODS

Three shunts each of the types Codman Hakim with SiphonGuard (flow-regulated), Miethke miniNAV with proSA (gravitational), and Medtronic Delta (membrane controlled) were tested. The shunts were compared on a novel in vitro setup that actively emulates the physiology of a shunted patient. This testing method allows determining the CSF drainage rates, resulting CSF volume, and intracranial pressure in the supine, sitting, and standing posture.

RESULTS

The flow-regulated ASDs avoided increased drainage by closing their primary flow path when drainage exceeded 1.39 ± 0.42 mL/min. However, with intraperitoneal pressure increased in standing posture, we observed reopening of the ASD in 3 out of 18 experiment repetitions. The adjustable gravitational ASDs allow independent opening pressures in horizontal and vertical orientation, but they did not provide constant drainage in upright posture (0.37 ± 0.03 mL/min and 0.26 ± 0.03 mL/min in sitting and standing posture, respectively). Consequently, adaptation to the individual patient is critical. The membrane-controlled ASDs stopped drainage in upright posture. This eliminates the risk of overdrainage, but leads to CSF accumulation up to the volume observed without shunting when the patient is upright.

CONCLUSIONS

While all tested ASDs reduced overdrainage, their actual performance will depend on a patient's specific needs because of the large variation in the way the ASDs influence CSF dynamics: while the flow-regulated shunts provide continuous drainage in upright posture, the gravitational ASDs allow and require additional adaptation, and the membrane-controlled ASDs show robust siphon prevention by a total stop of drainage.

Barrier dysfunction or drainage reduction: differentiating causes of CSF protein increase

BACKGROUND

Cerebrospinal fluid (CSF) protein analysis is an important element in the diagnostic chain for various central nervous system (CNS) pathologies. Among multiple existing approaches to interpreting measured protein levels, the Reiber diagram is particularly robust with respect to physiologic inter-individual variability, as it uses multiple subject-specific anchoring values. Beyond reliable identification of abnormal protein levels, the Reiber diagram has the potential to elucidate their pathophysiologic origin. In particular, both reduction of CSF drainage from the cranio-spinal space as well as blood-CNS barrier dysfunction have been suggested ρas possible causes of increased concentration of blood-derived proteins. However, there is disagreement on which of the two is the true cause.

METHODS

We designed two computational models to investigate the mechanisms governing protein distribution in the spinal CSF. With a one-dimensional model, we evaluated the distribution of albumin and immunoglobulin G (IgG), accounting for protein transport rates across blood-CNS barriers, CSF dynamics (including both dispersion induced by CSF pulsations and advection by mean CSF flow) and CSF drainage. Dispersion coefficients were determined a priori by computing the axisymmetric three-dimensional CSF dynamics and solute transport in a representative segment of the spinal canal.

RESULTS

Our models reproduce the empirically determined hyperbolic relation between albumin and IgG quotients. They indicate that variation in CSF drainage would yield a linear rather than the expected hyperbolic profile. In contrast, modelled barrier dysfunction reproduces the experimentally observed relation.

CONCLUSIONS

High levels of albumin identified in the Reiber diagram are more likely to originate from a barrier dysfunction than from a reduction in CSF drainage. Our in silico experiments further support the hypothesis of decreasing spinal CSF drainage in rostro-caudal direction and emphasize the physiological importance of pulsation-driven dispersion for the transport of large molecules in the CSF.

Is posture-related craniospinal compliance shift caused by jugular vein collapse? A theoretical analysis

Background

Postural changes are related to changes in cerebrospinal fluid (CSF) dynamics. While sitting up leads to a decrease in cranial CSF pressure, it also causes shifts in the craniospinal CSF volume and compliance distribution. We hypothesized that jugular vein collapse in upright posture is a major contributor to these shifts in CSF volume and compliance.

Methods

To test this hypothesis, we implemented a mathematical lumped-parameter model of the CSF system and the relevant parts of the cardiovascular system. In this model, the CSF and the venous system are each divided into a cranial and a spinal part. The pressures in these cranial and spinal portions differ by the posture-dependent hydrostatic pressure columns in the connecting vessels. Jugular collapse is represented by a reduction of the hydrostatic pressure difference between cranial and spinal veins. The CSF pressure–volume relationship is implemented as a function of the local CSF to venous pressure gradient. This implies that an increase in CSF volume leads to a simultaneous displacement of blood from adjacent veins. CSF pulsations driven by the cardiovascular system are introduced through a pulsating cranial arterial volume.

Results

In upright posture, the implemented CSF pressure–volume relationship shifts to lower cranial CSF pressures compared to the horizontal position, leading to a decrease in cranial CSF pressure when sitting up. Concurrently, the compliance of the spinal compartment decreases while the one of the cranial compartment increases. With this, in upright posture only 10% of the CSF system’s compliance is provided by the spinal compartment compared to 35% in horizontal posture. This reduction in spinal compliance is accompanied by a caudal shift of CSF volume. Also, the ability of the spinal CSF compartment to compensate for cerebral arterial volume pulsations reduces in upright posture, which in turn reduces the calculated craniospinal CSF flow pulsations.

Conclusion

The mathematical model enabled us to isolate the effect of jugular collapse and quantify the induced shifts of compliance and CSF volume. The good concordance of the modelled changes with clinically observed values indicates that jugular collapse can be considered a major contributor to CSF dynamics in upright posture.

Quantifying the relative contributions of different solute carriers to aggregate substrate transport

Determining the contributions of different transporter species to overall cellular transport is fundamental for understanding the physiological regulation of solutes. We calculated the relative activities of Solute Carrier (SLC) transporters using the Michaelis-Menten equation and global fitting to estimate the normalized maximum transport rate for each transporter (V_max). Data input were the normalized measured uptake of the essential neutral amino acid (AA) L-leucine (Leu) from concentration-dependence assays performed using Xenopus laevis oocytes. Our methodology was verified by calculating Leu and L-phenylalanine (Phe) data in the presence of competitive substrates and/or inhibitors. Among 9 potentially expressed endogenous X. laevis oocyte Leu transporter species, activities of only the uniporters SLC43A2/LAT4 (and/or SLC43A1/LAT3) and the sodium symporter SLC6A19/B⁰AT1 were required to account for total uptake. Furthermore, Leu and Phe uptake by heterologously expressed human SLC6A14/ATB^0,+ and SLC43A2/LAT4 was accurately calculated. This versatile systems biology approach is useful for analyses where the kinetics of each active protein species can be represented by the Hill equation. Furthermore, its applicable even in the absence of protein expression data. It could potentially be applied, for example, to quantify drug transporter activities in target cells to improve specificity.

The Bohr Effect Is Not a Likely Promoter of Renal Preglomerular Oxygen Shunting

The aim of this study was to evaluate whether possible preglomerular arterial-to-venous oxygen shunting is affected by the interaction between renal preglomerular carbon dioxide and oxygen transport. We hypothesized that a reverse (venous-to-arterial) shunting of carbon dioxide will increase partial pressure of carbon dioxide and decrease pH in the arteries and thereby lead to increased oxygen offloading and consequent oxygen shunting. To test this hypothesis, we employed a segment-wise three-dimensional computational model of coupled renal oxygen and carbon dioxide transport, wherein coupling is achieved by shifting the oxygen-hemoglobin dissociation curve in dependence of local changes in partial pressure of carbon dioxide and pH. The model suggests that primarily due to the high buffering capacity of blood, there is only marginally increased acidity in the preglomerular vasculature compared to systemic arterial blood caused by carbon dioxide shunting. Furthermore, effects of carbon dioxide transport do not promote but rather impair preglomerular oxygen shunting, as the increase in acidity is higher in the veins compared to that in the arteries. We conclude that while substantial arterial-to-venous oxygen shunting might take place in the postglomerular vasculature, the net amount of oxygen shunted at the preglomerular vasculature appears to be marginal.