Perivascular network segmentations derived from high

Jun 26, 2023

Scientific Reports volume 13, Article number: 9205 (2023) Cite this article

Metrics details

A custom segmentation workflow was applied to ex vivo high-field MR images of rat brains acquired following in vivo intraventricular contrast agent infusion to generate maps of the perivascular spaces (PVS). The resulting perivascular network segmentations enabled analysis of perivascular connections to the ventricles, parenchymal solute clearance, and dispersive solute transport within PVS. Numerous perivascular connections between the brain surface and the ventricles suggest the ventricles integrate into a PVS-mediated clearance system and raise the possibility of cerebrospinal fluid (CSF) return from the subarachnoid space to the ventricles via PVS. Assuming rapid solute exchange between the PVS and CSF spaces primarily by advection, the extensive perivascular network decreased the mean clearance distance from parenchyma to the nearest CSF compartment resulting in an over 21-fold reduction in the estimated diffusive clearance time scale, irrespective of solute diffusivity. This corresponds to an estimated diffusive clearance time scale under 10 min for amyloid-beta which suggests that the widespread distribution of PVS may render diffusion an effective parenchymal clearance mechanism. Additional analysis of oscillatory solute dispersion within PVS indicates that advection rather than dispersion is likely the primary transport mechanism for dissolved compounds greater than 66 kDa in the long (> 2 mm) perivascular segments identified here, although dispersion may be significant for smaller compounds in shorter perivascular segments.

Blood vessels in the brain are surrounded by slender perivascular spaces (PVS) that enable fluid exchange between interstitial fluid and cerebrospinal fluid (CSF) compartments1. These structures have garnered much attention recently because of the role they may play in a brain-wide clearance mechanism for toxic metabolic waste such as amyloid-beta (Aβ), the protein that accumulates in Alzheimer's disease2. While rapid imaging tracer uptake from the CSF3,4 and clearance from parenchyma5 have been observed, there is uncertainty1 regarding the mechanism and direction of transport6,7,8, the anatomy of arterial, capillary, and venous perivascular transport routes5, and the effect of aquaporin water channels2 and sleep9 on transport. Nonetheless, PVS-mediated transport in the brain may have significant implications not only for neurodegenerative diseases, but also for drug delivery to brain tissue10 and the migration of brain cancer11,12 and immune cells13.

While a number of studies have demonstrated uptake of imaging tracers into PVS near the brain surface3,4,14, few have examined deeper PVS and their connections to CSF in the cerebral ventricles and cisterns15,16. Histological sections following tracer uptake suggest an intricate and extensive network of PVS throughout the brain17, but the resolution of whole-brain imaging in vivo has limited analysis of the intact perivascular network to only the largest vessels17,18,19. There is a need for a whole-brain 3D map of major perivascular structures to analyze perivascular network properties relevant to clearance such as connections to internal CSF spaces, the distribution of PVS in parenchyma, and perivascular segment length. A map of PVS would enable potential perivascular and parenchymal transport mechanisms such as diffusion, dispersion, and advection to be evaluated via multi-scale mechanical modeling. This is required to better understand waste clearance and to accurately plan a variety of drug delivery techniques including intravenous, intrathecal, and convection-enhanced delivery to parenchyma.

Although a segmentation of the intact perivascular network in rats has not been published to the authors’ knowledge, several semi-automatic strategies have been developed to segment human PVS in clinical MR images20,21,22,23,24,25. In many of these strategies, the "tubeness" or "vesselness" of the image is determined by applying Frangi filters26 that are based on the spatial curvature of the image intensity. By applying a threshold to these tubeness images, a segmentation of PVS is produced and incorporated into a larger segmentation methodology often relying on deep learning techniques20,21. Despite its prevalence in human perivascular segmentation, tubeness thresholding has not been previously applied to segment PVS in rats or mice mainly because the resolution of in vivo MR images acquired during cerebrospinal contrast agent administration is not sufficiently high to resolve most PVS containing the contrast agent.

Magdoom, et al.16 present a set of ex vivo rat brain images with 40 μm isotropic voxels acquired at 17.6 T following intraventricular contrast agent infusion, revealing numerous connections between the perivascular network and the cerebral ventricles. In the present study, these images were newly processed to produce clearer visualizations of the perivascular network that facilitate the identification of major blood vessels surrounded by patent PVS. A custom segmentation workflow relying on tubeness was developed to create high-resolution 3D maps of the PVS containing contrast agent following intraventricular infusion. The resulting maps were used to isolate those perivascular segments in the vicinity of the ventricles and cisterns which may serve as clearance routes from parenchyma to deep CSF spaces. The extent to which the PVS facilitate parenchymal waste clearance by reducing the transport distance to the nearest CSF sink was quantified. An analytical transport model was also developed to evaluate dispersion caused by oscillatory perivascular flow as a potential solute transport mechanism in the longest PVS segments measured.

Extensive perivascular uptake of Gd-DTPA-albumin throughout the brain was observed following intraventricular infusion. Long perivascular segments, some over 3 mm in length, extend from the surface of the brain to the deep tissues surrounding the ventricular system (Fig. 1). The PVS are visible as bright, slender streaks in the coronal (Fig. 1a–h, Supplementary Video S1 and S2), sagittal (Fig. 1i–n, Supplementary Video S3 and S4), and transverse (Fig. 1o–t, Supplementary Video S5 and S6) planes. PVS surrounding subcortical penetrating arteries (scop) traverse the cortex to the subcortical white matter and periventricular tissue thus providing potential pathways for tracer between the dorsal subarachnoid space and the lateral ventricles (Fig. 1a–c, f, j, m). Perivascular uptake along smaller cortical vessels (cop) likely occurred given the dense array of bright, radial streaks (Fig. 1e–h), though individual vessels are not easily distinguishable.

Partial maximum intensity projections (pMIPs). Each panel is a projection across 30 slices parallel to the coronal (a–h), sagittal (i–n), and transverse (o–t) planes for Rat 3 (a–d, i–k, and o–q) and Rat 4 (e–h, l-n, and r–t). Proceeding alphabetically, each coronal panel is posterior to the preceding panel, each sagittal panel is lateral to the preceding panel, and each transverse panel is ventral to the preceding panel for each animal. Notable vessels are labeled (see Supplementary Table S1).

Perivascular uptake was evident along major ventral surface vessels including the anterior cerebral artery (acer), middle cerebral artery (mcer), posterior cerebral artery (pcer), internal carotids (ictd), and the basilar artery (bas) (Fig. 1b, f, i, q, t). Branches from these main vessels extending dorsally into the caudate and midbrain are surrounded by the brightest perivascular segments in the images and include the anterior (astr), medial (mstr), and posterior (pstr) striate arteries, ventral thalamic arteries (vth), and anterior choroidal arteries (ach) (Fig. 1a–c, e–g, j–k, m–n). A number of these PVS appear continuous with the ventricles (Fig. 1b, c, f–g) although the resolution is not sufficient to determine whether they are anatomically continuous with the ventricular CSF or they reside in periventricular tissue through which solutes could feasibly diffuse or flow to the ventricles. PVS surround numerous dorsally penetrating vessels in the midbrain and brain stem, especially near the longitudinal fissure, including the thalamo-perforating arteries (thp), median mesencephalic arteries (mmes), median pontine arteries (mpn), and the median medullary arteries (mmd) (Fig. 1i, l). The latter of these extend towards the fourth ventricle dorsal to the brainstem. PVS are also arranged radially about the aqueduct surrounding the lateral (lpaq) and dorsal (dpaq) periaqueductal arteries as well as the median pontine arteries (mpn) (Fig. 1d, h). In addition to strong contrast enhancement in the ventricles and PVS, weaker enhancement was present in the quadrigeminal and ambient cisterns (ac) (Fig. 1d, h, i–n, o, r). Contrast agent uptake in these cisterns is notable because they are continuous with external CSF spaces and are lined with numerous blood vessels27. The perivascular spaces surrounding some cisternal blood vessels (cbv) are evident in Fig. 1d, h, k, n.

Image segmentation revealed a dense network of PVS throughout the brain (Fig. 2, Supplementary Fig. S1) and numerous PVS that extend from near the ventricles to the brain surface and cisterns (Fig. 3, labeled vessels). The segmented PVS in sets of 30 coronal slices (Fig. 2a, c, h, j) include the vast majority of the PVS visible in the pMIPs for the same 30 slices (Fig. 1b, c, f, g). The 3D structure of the perivascular segments (Fig. 2b, d, i, k) depicts both long segments (striate arteries and subcortical penetrating arteries) and shorter segments within the cortex and caudate putamen. These shorter segments appear to be both smaller perivascular spaces near the limit of the imaging resolution and vessels oriented perpendicular to the coronal plane in the caudate putamen. While the segmented PVS in the cortex, brainstem, and surrounding the aqueduct were generally less than 1 mm in length, several perivascular segments originating from surface vessels and terminating near the ventricles were over 2.5 mm in length. The solid renderings indicate that no parenchymal region is conspicuously devoid of PVS, although inter-animal variability in cortical uptake is apparent (Fig. 2e–g, l–m, Supplementary Fig. S1). The perivascular uptake of contrast agent in the cortex, striatum, periaqueductal tissue, and the brain stem was quantified by calculating the percent of voxels segmented as PVS within volumes of interest (VOI) in each anatomical region (Table 1). The highest average percent volume (2.86 ± 0.934) was observed in the striatum while the remaining regions had average percent volumes between 1.02 (periaqueductal tissue) and 1.89 (cortex). The percent segmented perivascular volume was 25% greater in the cortex of Rat 4 than in Rat 3 which is consistent with the 3D renderings in Fig. 2. The MR images for Rat 1 revealed that in this animal contrast agent was partially infused into the periventricular tissue. This likely explains the higher percent segmented perivascular volume in the striatum and lower percent segmented perivascular volume in the periaqueduct and brain stem, far from the infusion site.

Perivascular and ventricle segmentations. The PVS (gold) and the ventricles (blue) segmented in two sets of 30 contiguous coronal slices are overlayed on the corresponding pMIPs for Rat 3 (a, c) and Rat 4 (h, j). These pMIPs are shown with overlayed segmentations for comparison with Fig. 1b, c, f, g, respectively. These 30-slice segmentations are presented alongside their 3D renderings for Rat 3 (b, d) and Rat 4 (i, k). All PVS (gold) and the ventricles (blue) are rendered in 3D for Rat 3 (e–g) and Rat 4 (l–n). Panels (e) and (l) are a rostral view of the coronal plane, (f, m) are a dorsal view of the transverse plane, and (g) and (n) are a left view of the sagittal plane. Notable vessels are labeled (see Supplementary Table S1).

Connected perivascular spaces occupying periventricular regions. The portions of the perivascular network containing segments that traverse a three-voxel margin around the ventricles are 3D-rendered for Rat 3 (a–c) and Rat 4 (d–f). The left column is a rostral view of the coronal plane, the central column is a dorsal view of the transverse plane, and the right column is a left view of the sagittal plane. The needle track in each rat is labeled with a red star. Notable vessels are labeled (see Supplementary Table S1).

Isolating the portion of the perivascular network extending outward from within a narrow region at most three voxels from the ventricular surface revealed connections between the lateral ventricles and major ventral surface vessels via the anterior choroidal arteries and the striate arteries (Fig. 3a, c, d, f). A margin of three voxels was selected because this region is narrow enough to enable relatively rapid transport of perivascular tracers into the ventricles. Perivascular connections from the lateral ventricles to the dorsal surface via subcortical penetrating arteries (Fig. 3) are also apparent. The PVS surrounding the transverse hippocampal arteries visible in the transverse pMIPs (Fig. 1o, r) appear to connect the caudal boundary of the lateral ventricles to the ambient cistern (Fig. 3b, e). The fourth ventricle is near PVS that traverse the brain stem and extend rostrally (Fig. 3, f).

The minimum distance between each parenchyma voxel and the nearest CSF space was computed both considering and not considering the PVS to be a CSF space (Fig. 4). In Rat 4, an animal exhibiting extensive perivascular tracer uptake, the parenchyma was on average 0.784 ± 0.516 mm away from the nearest CSF space in the absence of PVS, with some parenchymal regions located over 2 mm away from the nearest CSF space (Fig. 4g). These regions consisted of portions of the cortex and midbrain (Fig. 4a–c). When PVS were considered a CSF compartment, the mean minimum distance between parenchyma and the CSF was reduced to 0.169 ± 0.100 mm and all parenchyma voxels were within approximately 0.9 mm of a CSF space (Fig. 4g). This reduction in minimum distance represents an over 21-fold reduction in the time scale for diffusive clearance (see Eq. (2a)) from the parenchyma (Fig. 4h). The lowest average minimum clearance distance (MCD) across all animals was observed in the striatum VOI (0.150 ± 0.00438 mm) while the cortex, periaqueductal tissue, and brain stem VOI had values between 0.215 mm (brain stem) and 0.229 mm (periaqueduct) (Table 1). The maximum standard deviation of the MCD across animals was 62.4 μm in the periaqueductal tissue VOI. Supposing an interstitial flow were present in the parenchyma as posited by the glymphatic theory, the presence of PVS was estimated to reduce nearly five-fold the advective clearance time scale (see Eq. (2b)) from parenchyma.

Minimum distance between parenchyma and CSF spaces. The distances between parenchyma voxels and the nearest CSF space without PVS (a–c) and including PVS (d–f) for Rat 4 are presented as planar color images. The left column is a rostral view of the coronal plane, the central column is a dorsal view of the transverse plane, and the right column is a left view of the sagittal plane. (g) Distributions of minimum distance between parenchyma voxels and the CSF spaces with and without PVS. (h) Diffusive transport time scales for a range of solute effective diffusivities with and without PVS.

A transient transport analysis in a one-dimensional semi-infinite domain (Fig. 5a) was performed to ascertain whether oscillatory solute dispersion can explain the perivascular uptake observed in our intraventricular infusion experiment and following infusion into the cisterna magna as previously reported2,14,19. Assuming the CSF spaces were well-mixed and attained a high concentration shortly after the start of infusion, oscillatory perivascular flow resulted in an evolving perivascular concentration gradient with low-concentration leading regions traveling faster than trailing regions of higher concentration. These regions, or advancing solute fronts, are defined by their characteristic fraction \(\alpha\) of the CSF concentration. For albumin, the \(\alpha\) = 0.1 front traveled about 1 mm in 40 min whereas the \(\alpha\) = 0.8 front traveled about 150 μm in the same amount of time (Fig. 5c). By comparison, the brains with contrast agent contained numerous segments well over 2 mm long, some up to 4 mm (Fig. 5b). The front velocity was maximum at \(t\) = 0 because a step change in concentration was present at the boundary with the CSF (the concentration gradient was greatest here) and was just over 10 µm/s for \(\alpha\) = 0.1. Within a minute, the velocity for all fronts with \(\alpha\) > 0.1 dropped below 2 μm/s (Fig. 5d). The distance traveled due to oscillatory dispersion varies with the square root of the diffusion coefficient, meaning a four-fold change in this coefficient results in a two-fold change in the distance traveled. The times required for the \(\alpha\) = 0.5 front to traverse perivascular lengths of 250 μm and 1000 μm for a physiologically relevant range of molecular diffusivities and two literature values of dispersive enhancement, \(k\), are shown in Fig. 5e. The albumin solute front delayed between 8.11 min (\(k\) = 1.7)6 and 13.14 min (\(k\) = 1.05)7 to traverse 250 μm, but much longer (2.16 h at \(k\) = 1.7) to traverse 1000 μm, consistent with the position vs. time plot (\(k\) = 1.05) in Fig. 5c. The solute front for Aβ monomer (\(D\) = 180 μm2/s)28, being smaller than albumin, traversed both distances more quickly, but still delayed 1.00 h (\(k\) = 1.7) and 1.62 h (\(k\) = 1.05) to travel 1000 μm. The solute front for sodium ions traversed 1000 μm in 11.36 min (\(k\) = 1.05).

Oscillatory solute dispersion in PVS. (a) Dispersion model geometry consisting of a well-mixed CSF region with steady concentration \(C_{1}\) and an adjacent semi-infinite straight perivascular segment initially at concentration \(C\) = 0. (b) Measured lengths for a sample (n = 10) of the longest visually identifiable perivascular segments in Rat 5. The data points have been horizontally scattered to improve visibility. (c) Position and (d) velocity of advancing tracer fronts with concentrations \(\alpha C_{1}\) over time (k = 1.05). (e) Time to traverse 250 μm and 1000 μm for a range of solute molecular diffusivities. Solid and dashed lines denote a 5%7 (k = 1.05) and 70%6 (k = 1.7) increase in the diffusion coefficient (dispersive enhancement) caused by oscillatory flow.

A rapidly growing body of literature suggests the PVS communicate with the CSF spaces to form a brain-wide fluid exchange and waste clearance system1,29. Several pre-clinical studies have focused on dynamic tracer uptake into cortical PVS from the subarachnoid space using either in vivo fluorescence microscopy or post-mortem histology2,4,5. While cranial window microscopy enables high-resolution, real-time tracking of perivascular transport, a fixed field of view and shallow penetration depth limit analysis to a small, superficial portion of the perivascular network. As a result, the extent and connectivity of the brain-wide perivascular network and its functional relationship with the cerebral ventricles and cisterns are poorly understood. Bedussi, et al.15 observed PVS-mediated clearance to the ventricular system following infusion into the striatum, suggesting the ventricles function as a clearance sink for deep parenchymal tissue. Strong tracer signal in periventricular and cisternal tissues, as well as in a dense network of perivascular structures near the ventricles has also been observed following infusion into the cisterna magna17. These results imply that tracer arrives at the ventricles via deep parenchymal PVS since upstream intra-ventricular transport against ventricular flow seems less plausible. Indeed, dynamic contrast-enhanced MRI (DCE-MRI) data acquired during intracisternal infusion show that contrast agent moves through superficial perivascular spaces throughout the brain prior to arriving at the cerebral aqueduct indicating that minimal tracer moves upstream through the ventricular system during the infusion18.

In this study, we address the need for a high-resolution 3D visualization of the rat perivascular network and its connections to internal CSF spaces using high-field ex vivo MR images. While in vivo perivascular uptake17,18,19 of contrast agent following infusion into the cisterna magna has been observed with DCE-MRI, the imaging resolution in these studies limited analysis to large surface vessels. Maximum intensity projections through sub-volumes of the high-resolution (40 μm isotropic voxels) ex vivo MR images from Magdoom, et al.16 revealed numerous perivascular segments appearing to provide a direct pathway between the exterior surfaces of the brain and the deep tissues surrounding the cerebral ventricles (Fig. 1a–c, e–g). A number of these segments were between 2 and 4 mm in length (Fig. 5b). On the dorsal surface these include the subcortical penetrating arteries (scop), some of which, after traversing the cortex, run parallel to white matter tracks. Perivascular branches from large ventral surface vessels including the anterior choroidal arteries (ach), the striate arteries (astr, mstr, pstr), and the median medullary arteries (mmd) appeared to pass closely to or terminate near the ventricular surfaces, although the resolution is not sufficient to determine whether they are anatomically continuous with the ventricular CSF or they reside in periventricular tissue through which solutes could feasibly diffuse or flow to the ventricles. Numerous perivascular connections between the ambient cistern and the ventricles along trans-hippocampal arteries (trhi) and periaqueductal arteries (lpaq, dpaq) were evident. The clearance function of the trhi PVS may be especially relevant to Alzheimer's disease given that plaque formation and atrophy in the hippocampus is a feature of the disease30.

The geometries of the PVS penetrating the brain from major ventral and dorsal surface vessels are consistent with published rat and mouse cerebrovascular atlases. Extensive images of the rat vasculature presented in Scremin31 were referenced when identifying the blood vessel groups that may provide perivascular solute pathways from the brain surface to the ventricles. Subcortical penetrating arteries (scop) traverse the cortex with minimal tapering to the subcortical white matter31,32,33 and to periventricular parenchyma in some cases31. Long branches of relatively constant diameter originating from major ventral surface vessels traverse the caudate putamen and other subcortical structures, many extending to the periventricular parenchyma31,32.

The presence of perivascular pathways from the exterior surface of the brain to tissue surrounding the ventricles along with prior repeated observations of rapid tracer uptake along penetrating arterioles raises the possibility of CSF transport from the subarachnoid space to the ventricles via the PVS. In this case, parenchymal waste would be at least partially cleared to the ventricles by CSF return flow, forming a brain-wide clearance circulation. Although slightly different from the glymphatic circulation proposed in Iliff, et al.3 in which inflowing periarterial CSF flows through the parenchyma and eventually clears via perivenular spaces, these circulations are not mutually exclusive. Unfortunately, our MRI data do not indicate the direction of contrast agent transport to directly support perivascular CSF return to the ventricles. In our experiment, it is possible that perivascular transport during the infusion proceeded from the ventricles outward due to elevated ventricular pressure as opposed to traveling along the ventricular network to the subarachnoid space before uptake into the PVS. However, these experiments were conducted with the goal of outlining the perivascular network, not to track the natural evolution of tracer over time. The former required the infusion of a large amount of contrast agent to develop sufficient contrast in the extensive perivascular network for visualization with MRI. It is nonetheless worth noting that a previous infusion study using a higher flow rate (1.6 μL/min) for a longer period (60 min) produced only subtle changes in intracranial pressure34. The combination of in vivo MRI during contrast agent infusion into CSF and ex vivo whole-brain imaging at higher resolution is needed to observe brain-wide perivascular circulation and evaluate the hypothesized flow patterns described above.

The glymphatic theory posits that waste is cleared from the parenchyma by bulk flow of interstitial fluid from PVS surrounding arteries to PVS surrounding veins2. This assumption of flow is seemingly at odds with evidence from real-time iontophoresis (RTI) studies indicating that diffusion is the primary mode of solute transport in parenchyma35 and was soon challenged by computational studies suggesting unphysical pressure gradients would be required to drive advective parenchymal transport6,36,37,38. A later computational modeling study demonstrated how slow interstitial flow could explain errors observed in earlier RTI studies and argued that parenchymal transport is likely a combination of advection and diffusion39.

In this study, we consider the matter of parenchymal clearance from a geometric standpoint made possible by our perivascular space segmentation. By calculating the shortest distance from every parenchyma voxel to either the brain surface or the ventricles, the spatial distribution of minimum clearance distance (MCD) from parenchyma to the nearest CSF space was determined. For Rat 4, the MCD was on average 0.784 mm corresponding to a diffusion time scale for Aβ (\(D^{*}\) = 62.3 μm2/s28) of 2.74 h. By considering the PVS an extension of the CSF compartment, the average MCD was reduced to 0.169 mm and the diffusion time scale for Aβ was reduced over 21-fold to 7.67 min. The pervasive nature of the perivascular network and its rapid, primarily advective, solute exchange with the CSF appears to render diffusion a viable mechanism of Aβ clearance from parenchyma. This possibility is analogous to the systemic diffusive exchange of dissolved gases between the bloodstream and tissue and has been raised in the context of PVS-mediated clearance in previous reviews40,41. Moreover, the average MCD was likely overestimated due to smaller perivascular segments requiring higher resolution to distinguish. Previously reported spatial distributions of Peclet number in the rat brain for intrathecally infused gadoteric acid (Gd-DOTA) derived from DCE-MRI images suggest a transition from mainly advective transport in CSF spaces and larger PVS to transport that is mainly diffusive in the parenchyma42. This is consistent with hydraulic network models of the PVS and parenchyma38,43 which predict mostly moderate to high Peclet numbers in PVS (\(Pe\) > 0.1) but mostly low Peclet numbers (\(Pe\) < 0.1) in parenchyma. Of course, the reduction in clearance distance caused by the perivascular network also decreased the clearance time scale assuming purely advective parenchymal transport, although less dramatically, by a factor of 4.63. If interstitial bulk flow were indeed present, transport would proceed as a combination of advection and diffusion as dictated by the Peclet number \(Pe\) for a particular substance. For instance, an interstitial velocity of 0.367 µm/s would result in equal advective and diffusive time scales (\(Pe\) = 1) for Aβ over the average MCD when considering PVS to be a CSF compartment.

Rapid perivascular uptake of CSF tracers appears to rely on arterial pulsations3,44,45, but the mechanism by which pulsations drive perivascular transport is not clear41. By periodically varying the perivascular volume, arterial pulsations are expected to produce oscillatory fluid motion in the PVS, a prediction which has been indirectly confirmed by the oscillatory motion of perivascular tracers46,47. Vascular pulsations may generate net fluid motion through a peristaltic mechanism48 although this notion has been challenged by computational modeling results6,49 and theoretically by noting the pulse wavelength is much longer than the typical perivascular segment length38. Oscillatory fluid motion has the potential to enhance perivascular tracer transport without net fluid flow through Taylor dispersion50 which can be approximated as an increase in the effective diffusivity of the tracer. Dispersive transport in PVS has been addressed in recent computational studies6,7,51. Asgari, et al.6 predicted an increase in effective diffusivity greater than 27% while Troyetsky, et al.7 predicted a more modest increase of about 5%.

Here, we investigated dispersive perivascular uptake of albumin by considering the perivascular space a straight, semi-infinite channel adjacent to a well-mixed CSF compartment and tracking the motion of advancing tracer fronts with fractions \(\alpha\) of the CSF albumin concentration. Advancing fronts with lower concentration moved more quickly than fronts with higher concentration (see Eq. (6)). For instance, the \(\alpha\) = 0.1 front advanced at slightly over 10 μm/s initially and covered about 1 mm in 40 min whereas the \(\alpha\) = 0.8 front advanced at just under 2 μm/s initially and covered only about 150 μm in the same amount of time (Fig. 5c). This predicted penetration distance is shorter than the lengths of many perivascular segments (2–4 mm) showing strong, relatively uniform contrast agent enhancement after the 40 min infusion (Fig. 5b). The signal in these segments did not exhibit axial variations in concentration as would be expected from a dispersive transport process. Previously observed uptake into periarterial spaces appears to proceed faster than even the lowest concentration (\(\alpha\) = 0.1) front considered in the dispersive transport analysis. Microspheres move in the perivascular spaces surrounding pial arteries at 18.7 μm/s4 on average whereas the \(\alpha\) = 0.1 front drops to 3.43 μm/s after 10 s and to 1.09 μm/s after 100 s. The front velocities in dispersive transport also depend on molecular weight, contrary to evidence of molecular weight independent transport in periarterial spaces18. Estimates of solute speed in superficial CSF spaces including PVS derived from DCE-MRI by optimal mass transport theory (rOMT) are slower (~ 0.2 μm/s)42, but these estimates are derived from voxels larger than individual PVS and likely do not reflect perivascular velocities due to volume averaging. The average Peclet number in the same CSF region was nonetheless about 80, indicating the predominance of advective transport42. Together these results suggest that advection, rather than dispersion, was the primary perivascular uptake mechanism for the tracer in our experiment and in previous CSF infusion studies with comparable molecular weight tracers. Careful experiments by other groups involving a double-syringe system to maintain constant CSF volume also indicate that uptake is not driven by elevated infusion pressure52.

Dispersion will naturally make a greater transport contribution for higher diffusivity tracers in shorter perivascular segments. For instance, the \(\alpha\) = 0.5 front for Aβ monomer (\(D\) = 180 μm2/s)28, being smaller than albumin, traversed 250 μm in 6.06 min (\(k\) = 1.05), but required 1.62 h (\(k\) = 1.05) to travel 1000 μm (Fig. 5e). The \(\alpha\) = 0.5 solute front for sodium ions traversed 1000 μm in 11.36 min (\(k\) = 1.05) due to oscillatory dispersion. These results suggest that oscillatory dispersion may contribute substantially to transport in short perivascular segments or the start of long perivascular segments, but it loses relevance with increasing perivascular segment length and solute molecular weight. Even for the higher dispersive enhancement (\(k\) = 1.7) reported for \(D\) = 10 μm2/s in Asgari, et al.6 (enhancement tends to become smaller with increasing diffusivity), the Aβ monomer required 1.00 h to travel 1000 μm (Fig. 5e). Assuming the flow velocities observed in Mestre, et al.4 (18.7 μm/s) are present in the penetrating PVS observed in this study, solutes would traverse 4 mm (the length of the longest perivascular segments identified) by advection alone in 3.57 min suggesting that effective CSF exchange is more feasible by solute advection in PVS.

The perivascular network visualizations and segmentations presented here provide insight into the structure and function of the glymphatic system. There are numerous long perivascular segments that extend from the ventral and dorsal surfaces of the brain to the vicinity of the ventricles and cisterns that may serve as pathways for waste clearance to these internal CSF spaces and may integrate the glymphatic circulation into the larger CSF circulation. Assuming rapid fluid exchange between the segmented PVS and CSF spaces, the PVS can effectively reduce the average minimum distance from parenchyma to the nearest CSF sink thus reducing the time scale over which parenchymal waste may diffuse that distance over 21-fold, from an estimated time scale over 2 h to one under 10 min. Furthermore, a transient solute transport analysis within the PVS suggests that oscillatory dispersion is not the predominant transport mechanism for large tracers in the longest perivascular segments observed. The assumption of rapid clearance from the PVS in the parenchymal transport analysis and the inability of dispersion to rapidly transport tracer in the long perivascular segments observed imply that flow in the PVS is required to sustain effective clearance throughout the parenchyma by diffusion. The inflowing CSF may filter into the parenchyma as has been previously posited and/or flow towards the ventricles and cisterns as suggested by the geometries of PVS observed in this study. The fluid entering the parenchyma from the PVS is expected to flow at velocities much slower than velocities observed within the pial PVS given the geometric increase in flow cross-section with distance from the penetrating PVS and the low hydraulic conductivity in parenchyma. In this scenario, parenchyma transport would proceed as a combination of advection and diffusion as dictated by the Peclet number \(Pe\) for a particular substance which does not necessarily rule out diffusion as a significant, if not the predominant, transport mechanism in parenchyma.

The experimental design did not include a means to directly distinguish arteries and veins, making it difficult to determine whether tracer was in the PVS of veins. However, in Magdoom, et al.16 Evans blue dye was evident in the tissue surrounding the superior sagittal sinus suggesting the presence of tracer in venous PVS. Only a few infusion studies have shown tracer surrounding veins1,2,17. An in vivo MR angiography experiment prior to ventricular infusion could distinguish arteries and veins, albeit at lower resolution, and guide the interpretation of ex vivo high-resolution images of the PVS.

The maps of PVS are sufficiently resolved to analyze network connectivity, define segment centerlines, measure segment lengths, and calculate changes in MCD. However, the segmentations cannot provide an accurate readout of perivascular volume because of volume averaging within the 40-μm voxels. Contrast agent within perivascular structures smaller than 40 μm has the potential to affect the signal from multiple 40-μm voxels. The true size of the PVS is likely smaller than their apparent size in the ex vivo images because of volume averaging and contrast agent diffusion from the PVS into the parenchyma over the course of the infusion. However, albumin's relatively large size (66 kDa) limits its diffusion in brain tissue (16.3 μm2/s53) and may hinder its passage between astrocytic endfeet lining the exterior perivascular boundary2. The larger apparent size of the PVS produces an underestimation of the MCD. If a 17.85 μm diameter vessel lied at a shared vertex among four adjacent voxels which are segmented as perivascular voxels because of signal enhancement due to volume averaging and radial contrast agent diffusion, the minimum in-plane clearance distance would be underestimated by at most 19.36 μm. Assuming the majority of segmented PVS were for vessels of roughly 17.85 μm diameter54, taking this additional clearance distance into account would increase the average diffusive transport time scale for Aβ from 7.67 min to 9.53 min. It is reasonable to expect that such a vessel would enhance these four voxels because the albumin concentration in these voxels would be at minimum 45% of the perivascular albumin concentration due to radial albumin diffusion assuming a constant perivascular albumin concentration over 40 min (the duration of the intraventricular infusion).

The brain tissue was fixed after the intraventricular infusion procedure to prevent further albumin diffusion prior to imaging. However, fixation tends to shrink brain tissue slightly, leading to parenchyma clearance distances that are somewhat biased small. Ma, et al.55 report an average brain shrinkage due to formalin fixation in mice of 10.6% using volumetric measurements derived from MRI images. This is equivalent to a volume increase of 11.9% from the shrunken geometry to the in vivo geometry. Because volume increases as the cube of length, the corresponding increase in MCD is 3.81%. The diffusion time scale is proportional to length squared and would increase 7.76%, from 7.67 min to 8.26 min. The error in MCD associated with the apparent size of PVS and fixative-induced brain shrinkage produces modest changes in the diffusion time scale insufficient to undermine diffusion as an effective clearance mechanism from parenchyma given the abundance of perivascular spaces.

Extracellular spaces, including PVS, collapse following animal sacrifice presumably because of extracellular fluid influx into cells and consequent cell swelling56,57. It is unlikely, however, that the observed PVS are an artifact of their post-mortem collapse. The collapse of the PVS could potentially squeeze the contrast agent out of the spaces, detracting from their MRI visibility, making such spaces more difficult to see, not more apparent. Alternatively, fluid influx into cells following death might draw CSF and contrast agent into the PVS as in Du, et al.57, improving their visibility, but this should not outline spaces that were absent in vivo.

The ex vivo approach also prevented direct observation of perivascular transport direction and rate. In vivo imaging during the infusion would enable dynamic tracer observation at the expense of image resolution. As an alternative or addition, ex vivo images could be acquired following infusions of various durations. This approach would allow indirect observation of perivascular transport direction and rate without sacrificing imaging resolution. In the oscillatory transport analysis, perivascular inflow from the SAS driven by vessel pulsations or other sources of pressure gradients were not explicitly modeled. These factors could be incorporated into a model of perivascular flow and solute transport that includes the segmented geometry of PVS to further probe the mechanics of the glymphatic system by comparing simulated solute distributions with data from the in vivo and ex vivo experiments described above.

All animal procedures were approved by the University of Florida Institutional Animal Care and Use Committee and comply with relevant guidelines and regulations. The methods reported herein adhere to the recommendations in the ARRIVE guidelines. In Magdoom, et al.16, high-resolution magnetic resonance images (40 \({\upmu }\)m isotropic voxels) of PVS in five excised whole rat brains were acquired. Gadolinium-labeled human serum albumin, Gd-DTPA-albumin58, was infused into the lateral ventricle of each anesthetized rat (4% isoflurane in 1 L/min oxygen) at a rate of 1.5 μL/min for 40 min to allow distribution into the PVS. The animal was immediately exsanguinated then sacrificed via vascular perfusion with 0.9% sodium chloride solution followed by 4% formaldehyde fixative solution to crosslink the tracer to the brain tissue. After 2.5–3 days storage at 4 °C, the rat brain was excised and placed in Fluorinert oil (FC-43, 3 M Corp., St. Paul, MN, USA) in preparation for high-resolution imaging at 17.6 T (Bruker Biospin, Billerica, MA, USA). A \(T_{1}\)-weighted 3D gradient echo image with a FOV of 20 mm × 16 mm × 12 mm, a matrix size of 500 × 400 × 300, TR = 100 ms, TE = 0.3 ms, flip angle = 50°, and 7 averages was acquired in approximately 24 h. Two naïve control animals underwent the same experimental and imaging protocols in all ways except the stereotaxic surgery and intraventricular infusion. All animals (n = 7) were 2-month old male Sprague–Dawley rats weighing 280 – 300 g.

The rat brain was isolated in the MR images with rodent brain extraction tool, rBET59, and registered to Swanson's atlas60 with linear image registration tool, FLIRT61, before maximum intensity projections (MIPs) were produced in ImageJ62,63 as in Magdoom, et al.16. While this demonstrated the extent of perivascular uptake throughout the brain, further analysis of the perivascular network was made possible by the custom visualization and segmentation workflow presented in this study (Fig. 6). Additional details regarding the animal experiments and imaging protocol can be found in Magdoom, et al.16.

Image analysis workflow. Following the 24 h MRI session at 17.6 T, the brain was extracted and visualized with a custom maximum intensity projection technique (blue box). Then, brain regions and the PVS were segmented and rendered in 3D (red box). Lastly, PVS-ventricle connectivity and transport were analyzed (green box). The software applications used in each step, as indicated by the superscript letters, are as follows: (a) rodent Brain Extraction Tool (rBET); (b) FMRIB's Linear Registration Tool (FLIRT); (c) MATLAB r2018a; (d) ImageJ; (e) itk-SNAP.

The ventricles and PVS appear bright in the \(T_{1}\)-weigted 3D images because Gd-DTPA-albumin reduces \(T_{1}\)16. Despite the resulting contrast with surrounding tissue, the PVS are not readily apparent due to their small size, sometimes spanning a single voxel, and varied orientations. Computing a maximum intensity projection (MIP) of the data, as reported previously by Magdoom, et al.16, results in a set of projection images in which the brightest PVS and the ventricles are more apparent. However, smaller, dimmer PVS may not be projected and PVS may be occluded by the bright ventricles in the center of the brain. To overcome these limitations, partial maximum intensity projections (pMIPs) were created by projecting through 30-voxel thick rectangular volumes parallel to the coronal, sagittal, and transverse planes. Because continuous perivascular segments were often present in consecutive projections, shifting pMIPs (spMIPs) were used to visualize the full length of these structures throughout the brain (see Supplemental Video S1 to S6). In an spMIP, the projection span is shifted by a single voxel for each 30-voxel projection image instead of by 30 voxels as in a pMIP.

The PVS pose a segmentation challenge because of their small size and the surrounding tissue's non-uniform signal strength. Because the PVS surround blood vessels, they appear like slender, tubular, bright structures in many brain regions with a range of mean intensity values. To better distinguish them from these varied surroundings, a measure of the tubular quality or tubeness of the intensity field surrounding each voxel was computed in ImageJ64,65. The tubeness of a given voxel was defined as the geometric mean of the two most negative eigenvalues of the Hessian matrix for that voxel. The Hessian matrix elements are the second spatial derivatives of the image intensity field65

The matrix element subscripts \(i\) and \(j\) indicate the axes along which the differentiations are performed (see Supplementary Table S2 for all variables and parameters). The eigenvectors and eigenvalues of the Hessian matrix are thus the directions and magnitudes of the principal curvatures for the intensity field. For a region of high intensity pixels with lower intensity surroundings, the eigenvalues are negative, \(0 > \lambda_{1} > \lambda_{2} > \lambda_{3}\), and each eigenvalue is the image intensity curvature along its associated eigenvector. For a tubular intensity field, \(\lambda_{1} = 0\) and \(\lambda_{2} = \lambda_{3} \ll 0\)65 because the intensity does not vary along the tube axis, but has negative curvature perpendicular to this axis. Accordingly, in the ImageJ implementation, tubeness is defined as \(\sqrt {\lambda_{2} \lambda_{3} }\) and increases in value for more tubular intensity fields64. In this implementation, the tubeness field can be made sensitive to tubular structures of a certain size by applying a Gaussian filter prior to computing tubeness. Because many PVS only span a single voxel, the standard deviation for the Gaussian filter was set to 40 μm. The tubeness calculation is a variety of Frangi filtering26, a technique often applied to clinical MR images for perivascular space segmentation in humans20,21,22,25.

The perivascular space segmentations were generated by computing the tubeness of the 3D images and selecting voxels with tubeness above a threshold. Because the contrast agent greatly increased signal strength above levels in the naïve brains, the peaks of the intensity distribution for each brain were aligned by adjusting the data range before determining the tubeness threshold (Fig. 7a, b).

Signal intensity and tubeness distributions. (a) Signal intensity from the brains with contrast agent (Rat 1–5) and without contrast agent (Naïve 1–2). The signal intensity for each animal was mapped linearly between 0 and 1. (b) Signal intensity following intensity range remapping to align signal distribution peaks and create better perivascular contrast with the surrounding tissue. (c–d) Tubeness distribution for the brain (c) interior and (d) surface. The gray vertical line is the average of the tubeness values corresponding to the 0.95 cumulative distribution value in each naïve brain. This tubeness value is greater than the tubeness in 95% of the naïve voxels on average.

The threshold was set to exceed the tubeness values of 95% of voxels in the naïve brains (tubeness = 3.75; Fig. 7c). This was done to exclude background tubeness not associated with the uptake of contrast agent in the PVS. The brain's exterior surfaces had higher background tubeness levels because the interface between the parenchyma and the dark surroundings produced high spatial intensity gradients in the naïve brains. Therefore, a higher tubeness threshold (tubeness = 8.08) was used to segment the PVS of surface blood vessels (Fig. 7d).

Several segmentations of anatomical structures including the ventricles and cerebellum were created in itk-SNAP66 using both manual outlining and the "snake" semi-autosegmentation tool. The ventricle network appeared bright because of the infused contrast agent which facilitated segmentation by comparison with the Paxinos rat brain atlas27. The tubeness-based segmentation of the PVS was improved by removing the cerebellum because its internal boundaries made PVS difficult to distinguish by means of tubeness alone. Volumes of interest (VOI) in the cortex, striatum, periaqueductal tissue, and the brain stem spanning 30 coronal slices each were manually outlined to compare perivascular voxel volume percent and mean MCD in these anatomical regions across animals (Fig. 8). The cortex and striatum VOI occupied the same set of 30 slices, and the striatum was defined as the tissue ventral and medial to the cortex within these slices.

Brain volumes of interest (VOI) for perivascular space segmentation analysis. A coronal slice from the VOI in the (a) cortex, (b) striatum, (c) periaqueductal tissue, and (d) brain stem. Each VOI spans 30 coronal slices.

The connections between the PVS and ventricle network were examined with a custom region-grow algorithm that iteratively expanded the ventricle segmentation into the perivascular segmentation. In each iteration, perivascular voxels immediately adjacent to a ventricle voxel were added to the growing ventricle segmentation for a total of 100 iterations. Prior to beginning this procedure, the ventricle segmentation was enlarged by three voxels to capture PVS within 208 µm of the ventricular surface, the diagonal distance across three 40-μm cubical voxels aligned diagonally. A margin of three voxels was selected because this length is short enough to enable relatively rapid transport of perivascular tracers into the ventricles. This produced a segmentation that included the ventricles and portions of the perivascular network continuous with perivascular segments near the ventricles.

PVS length estimates were made for a selection of long, continuous segments penetrating from the ventral and dorsal brain surfaces using a modified region-grow algorithm. A seed voxel was placed at the surface end of each perivascular segment and the centroid of voxels added after each growth iteration was recorded. The distances between these centroids were computed and summed to yield an estimate of the total perivascular segment length.

The extent to which the perivascular network facilitates exchange between CSF and interstitial fluid was quantified by calculating the distance between parenchyma voxels and the nearest CSF compartment, \(L\), both considering and not considering the PVS to be a CSF space. In this analysis, the PVS were considered CSF compartments given numerous observations of rapid exchange between PVS and larger CSF spaces. The minimum distances were determined by computing a Euclidean distance transform67 of the CSF segmentation in MATLAB (MATLAB v. 9.4.0.813654 (R2018a), The MathWorks, Inc., Natick, MA). Clearance time scales \(\tau_{d}\) and \(\tau_{a}\) for a range of relevant solute diffusivities were computed assuming either purely diffusive (Eq. (2a)) or advective (Eq. (2b)) transport through the parenchyma.

In Eq. (2a), \(D^{*}\) is the effective diffusivity of the solute in parenchyma and \(u\) is the interstitial velocity magnitude. \(D^{*}\) is lower than the solute free diffusivity \(D\) because of parenchymal structures that hinder random solute motion. The reduction in parenchymal clearance time due to the perivascular network was estimated for a physiologically relevant range of \(D^{*}\) values (101–103 µm2/s37) including that of Aβ (62.3 μm2/s28).

A perivascular transport model was developed to investigate oscillatory solute dispersion in the longest perivascular segments observed in this study. This model assumes the CSF in the subarachnoid space is well mixed and maintains a constant solute concentration \(C_{1}\) while the perivascular segments contain free fluid initially at concentration \(C_{0}\) = 0. Because the perivascular segments are long and slender, a one-dimensional semi-infinite geometry (Fig. 5a) was considered for which the origin is located at the interface between the CSF reservoir and the perivascular segment. Dynamic changes in perivascular geometry due to vessel wall motion were not modeled since these are small relative to the perivascular gap width3,68. Dispersion due to oscillatory flow in an annulus has been previously modeled as an enhanced diffusion process where the solute free diffusivity \(D\) is increased by a factor \(k\)6,7,50. Perivascular inflow from the SAS driven by vessel pulsations or other sources of pressure gradients were not explicitly modeled; only the dispersive effect of oscillatory flow on the evolution of solute concentration was considered. Therefore, the governing equation for this scenario is

Given the absence of a geometric length scale, a self-similar solution was sought for non-dimensional concentration \(\alpha\) in terms of variable \(\eta\) where

The self-similar solution is defined in terms of the error function as

The position \(x\) and velocity \(v\) of the material point with a fraction \(\alpha\) of the CSF concentration \(C_{1}\) in the perivascular segment over time is given by

Equation (6) represents the kinematics of dispersive tracer uptake into the PVS. Bovine serum albumin (\(D\) = 83 μm2/s53) with maximum enhancement \(k\) = 1.057 and \(k\) = 1.76 was considered a representative large molecular weight tracer.

The data is available upon request from the corresponding author.

Wardlaw, J. M. et al. Perivascular spaces in the brain: anatomy, physiology and pathology. Nat. Rev. Neurol. 16, 137–153 (2020).

PubMed Google Scholar

Iliff, J. J. et al. A paravascular pathway facilitates CSF flow through the brain parenchyma and the clearance of interstitial solutes, including amyloid beta. Sci. Transl. Med. 4, 11 (2012).

Google Scholar

Iliff, J. J. et al. Cerebral arterial pulsation drives paravascular CSF-interstitial fluid exchange in the murine brain. J. Neurosci. 33, 18190–18199 (2013).

CAS PubMed Google Scholar

Mestre, H. et al. Flow of cerebrospinal fluid is driven by arterial pulsations and is reduced in hypertension. Nat. Commun. 9 (2018).

Morris, A. W. J. et al. Vascular basement membranes as pathways for the passage of fluid into and out of the brain. Acta Neuropathol. 131, 725–736 (2016).

CAS PubMed Google Scholar

Asgari, M., de Zelicourt, D. & Kurtcuoglu, V. Glymphatic solute transport does not require bulk flow. Sci. Rep. 6, 11 (2016).

Google Scholar

Troyetsky, D. E., Tithof, J., Thomas, J. H. & Kelley, D. H. Dispersion as a waste-clearance mechanism in flow through penetrating perivascular spaces in the brain. Sci. Rep. 11, 12 (2021).

Google Scholar

Martinac, A. D. & Bilston, L. E. Computational modelling of fluid and solute transport in the brain. Biomech. Model. Mechanobiol. 19, 781–800 (2020).

PubMed Google Scholar

Xie, L. L. et al. Sleep drives metabolite clearance from the adult brain. Science 342, 373–377 (2013).

ADS CAS PubMed Google Scholar

Foley, C. P., Nishimura, N., Neeves, K. B., Schaffer, C. B. & Olbricht, W. L. Real-time imaging of perivascular transport of nanoparticles during convection-enhanced delivery in the rat cortex. Ann. Biomed. Eng. 40, 292–303 (2012).

PubMed Google Scholar

Pacioni, S. et al. Brain invasion along perivascular spaces by glioma cells: relationship with blood-brain barrier. Cancers 12, 17 (2020).

Google Scholar

Watkins, S. et al. Disruption of astrocyte-vascular coupling and the blood-brain barrier by invading glioma cells. Nat. Commun. 5, 15 (2014).

Google Scholar

Troili, F. et al. Perivascular unit: This must be the place. The anatomical crossroad between the immune, vascular and nervous system. Front. Neuroanat. 14, 17 (2020).

CAS PubMed Google Scholar

Smith, A. J., Yao, X. M., Dix, J. A., Jin, B. J. & Verkman, A. S. Test of the "glymphatic" hypothesis demonstrates diffusive and aquaporin-4-independent solute transport in rodent brain parenchyma. Elife 6, 16 (2017).

Google Scholar

Bedussi, B. et al. Clearance from the mouse brain by convection of interstitial fluid towards the ventricular system. Fluids Barriers CNS 12, 13 (2015).

Google Scholar

Magdoom, K. N. et al. MRI of whole rat brain perivascular network reveals role for ventricles in brain waste clearance. Sci. Rep. 9, 11 (2019).

Google Scholar

Pizzo, M. E. et al. Intrathecal antibody distribution in the rat brain: surface diffusion, perivascular transport and osmotic enhancement of delivery. J. Physiol. Lond. 596, 445–475 (2018).

CAS PubMed Google Scholar

Iliff, J. J. et al. Brain-wide pathway for waste clearance captured by contrast-enhanced MRI. J. Clin. Investig. 123, 1299–1309 (2013).

CAS PubMed Google Scholar

Koundal, S. et al. Optimal mass transport with lagrangian workflow reveals advective and diffusion driven solute transport in the glymphatic system (vol 10, 1990, 2020). Sci. Rep. 10, 1 (2020).

Google Scholar

Park, S. H., Zong, X. P., Gao, Y. Z., Lin, W. L. & Shen, D. G. Segmentation of perivascular spaces in 7 T MR image using auto-context model with orientation-normalized features. Neuroimage 134, 223–235 (2016).

PubMed Google Scholar

Zhang, J. et al. Structured learning for 3-D perivascular space segmentation using vascular features. IEEE Trans. Biomed. Eng. 64, 2803–2812 (2017).

PubMed Google Scholar

Ballerini, L. et al. Perivascular spaces segmentation in brain MRI using optimal 3D filtering. Sci. Rep. 8, 11 (2018).

Google Scholar

Lian, C. F. et al. Multi-channel multi-scale fully convolutional network for 3D perivascular spaces segmentation in 7T MR images. Med. Image Anal. 46, 106–117 (2018).

PubMed Google Scholar

Boutinaud, P. et al. 3D segmentation of perivascular spaces on T1-weighted 3 tesla MR images with a convolutional autoencoder and a U-shaped neural network. Front. Neuroinform. 15, 21 (2021).

Google Scholar

Langan, M. T. et al. Semi-automated segmentation and quantification of perivascular spaces at 7 tesla in COVID-19. Front. Neurol. 13, 11 (2022).

Google Scholar

Frangi, A. F., Niessen, W. J., Vincken, K. L. & Viergever, M. A. Multiscale vessel enhancement filtering. In Medical Image Computing and Computer—Assisted Intervention—Miccai’98 Vol. 1496 (eds Wells, W. M. et al.) 130–137 (Springer, 1998).

Google Scholar

Paxinos, G. & Watson, C. The Rat Brain in Stereotaxic Coordinates Vol. 6 (Academic Press, New York, 2007).

Google Scholar

Waters, J. The concentration of soluble extracellular amyloid-beta protein in acute brain slices from CRND8 mice. PLoS ONE 5, 16 (2010).

Google Scholar

Benveniste, H. et al. The glymphatic system and waste clearance with brain aging: A review. Gerontology 65, 106–119 (2019).

PubMed Google Scholar

DeTure, M. A. & Dickson, D. W. The neuropathological diagnosis of Alzheimer's disease. Mol. Neurodegener. 14, 18 (2019).

Google Scholar

Scremin, O. U. Cerebral vascular system. In The Rat Nervous System 4th edn (ed. Paxinos, G.) 985–1036 (Elsevier, 2015).

Google Scholar

Kirst, C. et al. Mapping the fine-scale organization and plasticity of the brain vasculature. Cell 180, 780-+ (2020).

CAS PubMed Google Scholar

Bonney, S. K. et al. Public volume electron microscopy data: An essential resource to study the brain microvasculature. Front. Cell Dev. Biol. 10 (2022).

Yang, L. J. et al. Evaluating glymphatic pathway function utilizing clinically relevant intrathecal infusion of CSF tracer. J. Transl. Med. 11, 9 (2013).

Google Scholar

Sykova, E. & Nicholson, C. Diffusion in brain extracellular space. Physiol. Rev. 88, 1277–1340 (2008).

CAS PubMed Google Scholar

Jin, B. J., Smith, A. J. & Verkman, A. S. Spatial model of convective solute transport in brain extracellular space does not support a "glymphatic" mechanism. J. Gen. Physiol. 148, 489–501 (2016).

CAS PubMed Google Scholar

Holter, K. E. et al. Interstitial solute transport in 3D reconstructed neuropil occurs by diffusion rather than bulk flow. Proc. Natl. Acad. Sci. USA 114, 9894–9899 (2017).

ADS CAS PubMed Google Scholar

Rey, J. & Sarntinoranont, M. Pulsatile flow drivers in brain parenchyma and perivascular spaces: A resistance network model study. Fluids Barriers CNS 15, 11 (2018).

Google Scholar

Ray, L., Iliff, J. J. & Heys, J. J. Analysis of convective and diffusive transport in the brain interstitium. Fluids Barriers CNS 16, 18 (2019).

Google Scholar

Abbott, N. J., Pizzo, M. E., Preston, J. E., Janigro, D. & Thorne, R. G. The role of brain barriers in fluid movement in the CNS: Is there a "glymphatic" system?. Acta Neuropathol. 135, 387–407 (2018).

CAS PubMed Google Scholar

Thomas, J. H. Fluid dynamics of cerebrospinal fluid flow in perivascular spaces. J. R. Soc. Interface 16, 20190572 (2019).

CAS PubMed Google Scholar

Chen, X. et al. Cerebral amyloid angiopathy is associated with glymphatic transport reduction and time-delayed solute drainage along the neck arteries. Nat. Aging 2, 214–223 (2022).

CAS PubMed Google Scholar

Boster, K. A. S., Tithof, J., Cook, D. D., Thomas, J. H. & Kelley, D. H. Sensitivity analysis on a network model of glymphatic flow. J. R. Soc. Interface 19, 14 (2022).

Google Scholar

Rennels, M. L., Gregory, T. F., Blaumanis, O. R., Fujimoto, K. & Grady, P. A. Evidence for a paravascular fluid circulation in the mammalian central nervous-system, provided by the rapid distribution of tracer protein throughout the brain from the subarachnoid space. Brain Res. 326, 47–63 (1985).

CAS PubMed Google Scholar

Hadaczek, P. et al. The "perivascular pump" driven by arterial pulsation is a powerful mechanism for the distribution of therapeutic molecules within the brain. Mol. Ther. 14, 69–78 (2006).

CAS PubMed Google Scholar

Ichimura, T., Fraser, P. A. & Cserr, H. F. Distribution of extracellular tracers in perivascular spaces of the rat-brain. Brain Res. 545, 103–113 (1991).

CAS PubMed Google Scholar

Bedussi, B., Almasian, M., de Vos, J., VanBavel, E. & Bakker, E. Paravascular spaces at the brain surface: Low resistance pathways for cerebrospinal fluid flow. J. Cereb. Blood Flow Metab. 38, 719–726 (2018).

PubMed Google Scholar

Wang, P. & Olbricht, W. L. Fluid mechanics in the perivascular space. J. Theor. Biol. 274, 52–57 (2011).

ADS MathSciNet PubMed MATH Google Scholar

Daversin-Catty, C., Vinje, V., Mardal, K. A. & Rognes, M. E. The mechanisms behind perivascular fluid flow. PLoS ONE 15, 20 (2020).

Google Scholar

Aris, R. On the dispersion of a solute in pulsating flow through a tube. Proc. R. Soc. Lond. Ser. A Math. Phys. Sci. 259, 370–376 (1960).

ADS MathSciNet CAS MATH Google Scholar

Sharp, M. K., Carare, R. O. & Martin, B. A. Dispersion in porous media in oscillatory flow between flat plates: Applications to intrathecal, periarterial and paraarterial solute transport in the central nervous system. Fluids Barriers CNS 16, 1–17 (2019).

Google Scholar

Raghunandan, A. et al. Bulk flow of cerebrospinal fluid observed in periarterial spaces is not an artifact of injection. Elife 10, 15 (2021).

Google Scholar

Tao, L. & Nicholson, C. Diffusion of albumins in rat cortical slices and relevance to volume transmission. Neuroscience 75, 839–847 (1996).

CAS PubMed Google Scholar

Zhang, M. Q. et al. Ultra-high-resolution 3D digitalized imaging of the cerebral angioarchitecture in rats using synchrotron radiation. Sci. Rep. 5, 12 (2015).

Google Scholar

Ma, Y. et al. In vivo 3D digital atlas database of the adult C57BL/6J mouse brain by magnetic resonance microscopy. Front. Neuroanat. 2, 10 (2008).

Google Scholar

Mestre, H., Mori, Y. & Nedergaard, M. The brain's glymphatic system: current controversies. Trends Neurosci. 43, 458–466 (2020).

CAS PubMed Google Scholar

Du, T. et al. Cerebrospinal fluid is a significant fluid source for anoxic cerebral oedema. Brain 145, 787–797 (2022).

PubMed Google Scholar

Ogan, M. D. et al. Albumin labeled with GD-DTPA—an intravascular contrast-enhancing agent for magnetic-resonance blood pool imaging—preparation and characterization. Investig. Radiol. 22, 665–671 (1987).

CAS Google Scholar

Wood, T.C., Lythgoe, D.J. & Williams, S.C.R. rBET: Making BET work for Rodent Brains. Proc. Intl. Soc. Mag. Reson. Med. 21 (2013).

Swanson, L. W. Brain Maps III: Structure of the Rat Brain (Elsevier, 2004).

Google Scholar

Jenkinson, M., Bannister, P., Brady, M. & Smith, S. Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage 17, 825–841 (2002).

PubMed Google Scholar

Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH image to ImageJ: 25 years of image analysis. Nat. Methods 9, 671–675 (2012).

CAS PubMed Google Scholar

Schindelin, J. et al. Fiji: An open-source platform for biological-image analysis. Nat. Methods 9, 676–682 (2012).

CAS PubMed Google Scholar

Longair, M., Preibisch, S. & Schindelin, J. "Tubeness": PlugIn for ImageJ. https://www.longair.net/edinburgh/imagej/tubeness/

Sato, Y. et al. 3D multi-scale line filter for segmentation and visualization of curvilinear structures in medical images. In Virtual Reality and Robotics in Medicine and Medical Robotics and Computer-Assisted Surgery Vol. 1205 (eds Troccaz, J. et al.) 213–222 (Springer, 1997).

Google Scholar

Yushkevich, P. A. et al. User-guided 3D active contour segmentation of anatomical structures: Significantly improved efficiency and reliability. Neuroimage 31, 1116–1128 (2006).

PubMed Google Scholar

Maurer, C. R., Qi, R. S. & Raghavan, V. A linear time algorithm for computing exact Euclidean distance transforms of binary images in arbitrary dimensions. IEEE Trans. Pattern Anal. Mach. Intell. 25, 265–270 (2003).

Google Scholar

Kedarasetti, R. T., Drew, P. J. & Costanzo, F. Arterial pulsations drive oscillatory flow of CSF but not directional pumping. Sci. Rep. 10, 12 (2020).

Google Scholar

Download references

We thank Dr. Robert Brasch (University of California, San Francisco) for generously supplying the gadolinium contrast agent used in this study. The magnetic resonance images presented were acquired in the Advanced MRI/S (AMRIS) Facility at the McKnight Brain Institute at the University of Florida, a part of the National High Magnetic Field Laboratory (NHMFL) supported by the National Science Foundation Cooperative Agreement DMR-1157490, the State of Florida, and the U.S. Department of Energy. This study was partly funded by the Bobby Jones Chiari & Syringomyelia Foundation (bobbyjonescsf.org) and an NHMFL user grant. We also thank Dr. Tavarekere N. Nagaraja from Henry Ford Hospital in Detroit, MI, USA for sharing his experience with CSF tracer distributions following intraventricular infusion. We also thank Rasheed M. Anjum for assisting with the analysis of radial contrast agent diffusion from the perivascular space.

Department of Mechanical and Aerospace Engineering, University of Florida, PO BOX 116250, Gainesville, FL, 32611, USA

Julian A. Rey, Uzair M. Farid, Christopher M. Najjoum, Kulam Najmudeen Magdoom & Malisa Sarntinoranont

Department of Biochemistry and Molecular Biology, University of Florida, Gainesville, FL, USA

Alec Brown & Thomas H. Mareci

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

J.A.R. assisted with the animal intraventricular infusion experiments, designed and executed the perivascular space segmentation and transport analysis workflow, and wrote the manuscript. U.M.F. assisted in the development and execution of custom perivascular space visualization and segmentation code. C.M.N. was responsible for the manual segmentation of the ventricles and cerebellum. A.B. performed the preliminary brain segmentation and anatomical registration. K.N.M. conducted the intraventricular infusion experiments, prepared and imaged the samples, and lead the preliminary image processing. T.H.M. oversaw the magnetic resonance imaging and contributed to the image post-processing and perivascular space segmentation strategies. M.S. designed and obtained funding for the intraventricular infusion experiments. M.S. also oversaw the perivascular space visualization and segmentation workflow development, the parenchymal and perivascular transport analyses, and the manuscript composition. All authors reviewed the manuscript.

Correspondence to Malisa Sarntinoranont.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Rey, J.A., Farid, U.M., Najjoum, C.M. et al. Perivascular network segmentations derived from high-field MRI and their implications for perivascular and parenchymal mass transport in the rat brain. Sci Rep 13, 9205 (2023). https://doi.org/10.1038/s41598-023-34850-0

Download citation

Received: 31 August 2022

Accepted: 09 May 2023

Published: 06 June 2023

DOI: https://doi.org/10.1038/s41598-023-34850-0

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.