What this is
- This research investigates in the rat brain using H217O MRI.
- It challenges previous assumptions about the speed and extent of glymphatic flow.
- The study confirms the critical role of channels in this transport process.
Essence
- in the rat brain is significantly faster and more extensive than previously understood, as shown by direct imaging with H217O MRI. The study emphasizes the importance of channels in facilitating this transport.
Key takeaways
- using H217O is markedly quicker than with traditional Gd-DTPA tracers, indicating a more efficient waste clearance mechanism in the brain.
- The study shows that blocking channels leads to an 80% reduction in , underscoring their critical role in this process.
Caveats
- The study is limited to a rat model, which may not fully represent human glymphatic function. Further research is needed to confirm these findings in humans.
- The use of anesthesia during the experiments may influence rates, potentially affecting the generalizability of the results.
Definitions
- glymphatic transport: A proposed mechanism for solute transport and waste clearance in the brain involving the flow of cerebrospinal fluid through the interstitial space.
- aquaporin-4: A water channel protein in the brain that facilitates the movement of water across cell membranes, crucial for glymphatic transport.
AI simplified
Introduction
The transport of solutes within the brain parenchyma is of fundamental importance to nutrient delivery and the clearance of metabolites, neurotransmitters, and toxic macromolecules (e.g., β-amyloid). Conventionally, interstitial solutes were considered to be transported via diffusion, but recent evidence suggests an additional bulk flow of the interstitial fluid (ISF). The glymphatic (glial-lymphatic) hypothesis proposes that subarachnoid cerebrospinal fluid (CSF) is driven by arterial pulsation along the perivascular space surrounding penetrating arteries, with influx into the brain interstitium mediated by the astroglial water channel aquaporin-4 (AQP4) (). It is proposed that this influx results in the bulk flow of ISF, which then exits along perivenous spaces, providing an efficient clearance mechanism for waste products from the parenchyma. Glymphatic transport has been shown to increase during sleep, and this increases the clearance of β-amyloid in mice, potentially explaining the decades old mystery of why sleep is restorative (). Glymphatic transport has been further shown to increase in ischemic stroke and play an important role in poststroke edema (). Conversely, decreased glymphatic transport has been shown in animal models of Alzheimer’s disease (), vascular dementia (), and traumatic brain injury (). Moreover, evidence that the glymphatic pathway is present in humans has been obtained via the intrathecal injection of MRI tracers, with delayed tracer clearance in a cohort of patients with dementia (). 1 2 3 4 5 6 7
However, the glymphatic hypothesis has proven to be highly controversial (,), particularly regarding the proposed bulk flow of ISF (–) and the mediating role of AQP4 (,,). A key limitation of previous studies has been the use of tracer molecules that are much larger than water (18 Da) — for example, in ex vivo fluorescence imaging (700–3000 Da) () and in vivo MRI (Gd-DTPA, 938 Da) (). The fact that these tracers cannot be transported by AQP4 channels most likely means that they underestimate the true magnitude of ISF flow, which will depend on the rate at which CSF water molecules enter the interstitium. We reasoned that using water molecules as a tracer would yield insights into glymphatic transport and the role of AQP4 channels. 8 9 10 12 8 13 14 1 15
Water molecules have previously been used as tracers by exploiting various isotopes. For example, tritiated water (HO) has been used to investigate CSF transport by measuring the radiation levels of blood samples using a liquid scintillation counter (). HO has been imaged with positron emission tomography (), though its short half-life (122.2 seconds) makes it impractical for studying glymphatic function. The stableO isotope is not radioactive but possesses a quadrupolar nucleus (I = 5/2) that can be imaged using MRI. Although direct imaging of theO signal is difficult due to its low gyromagnetic ratio and low natural abundance (0.037%) (), the presence of HO can be detected indirectly via the effect of the quadrupolarO nucleus on theH MRI signal (), where the quadrupolarO nucleus reduces the T2 relaxation times of bonded hydrogen nuclei. This effect is further magnified by the exchange of hydrogen atoms with other surrounding water molecules. HO produced fromOgas has been used to image metabolism (), and HO has also been administered as a tracer to image cerebral blood flow (). In a series of studies, the Nakada group employed indirectH detection of HO to investigate the interstitial circulation. Using an i.v. bolus injection of HO (20%O enriched), they showed that exchange of blood water with the ISF and CSF was dependent on AQP4 and not AQP1 (,). They also used the same method to demonstrate that water influx into CSF was significantly impaired in surfactant protein–bearing (SP-bearing) transgenic mice (). 3 15 17 17 17 17 1 17 17 17 17 1 17 17 17 2 2 2 2 2 2 2 2 16 17 18 19 20 21 22 23 24
In order to investigate the glymphatic system more directly, we decided to infuse HO into the CSF at the cisterna magna, thus replicating the delivery approach used by Iliff et al. () in their seminal Gd-DTPA MRI experiments. To avoid elevating the intracranial pressure, the HO has to be infused slowly (1.8 μL/min) (), thus limiting the amount that can be delivered. To make possible the detection of smallH MRI signal changes, it was necessary to use the highest available enrichment of HO (90%O enriched) combined with a high signal-to-noise ratio (SNR) T2-weighted imaging sequence (SNR > 200). 2 2 2 17 17 1 17 17 15 25
Results
To directly image glymphatic water transport within the rat brain, the highly enriched HO tracer (90%O enriched) was infused into the CSF of anesthetized rats at the cisterna magna. Serial MRI revealed glymphatic transport that was strikingly more rapid and extensive than previously observed using conventional Gd-DTPA tracers () (). Movement of the Gd-DTPA tracer through the subarachnoid space was slightly slowed compared with that of HO. Moreover, as previously observed (), the Gd-DTPA tracer was slow to penetrate the parenchyma, resulting in a buildup of Gd-DTPA concentration in the subarachnoid space and ventricles. 2 2 17 17 17 Figure 1B 15 15
By contrast to the Gd-DTPA tracer, HO rapidly penetrates the parenchyma in all brain regions (). Thus, unlike Gd-DTPA, the concentration of HO does not build up in the subarachnoid space and ventricles (). This is best appreciated by viewing(supplemental material available online with this article;) that shows, side-by-side, the temporal evolution of both tracers. For example, HO fully penetrates the rostral cortex within 10 minutes (), whereas, as previously observed (), Gd-DTPA does not penetrate even after 85 minutes. There is much debate in the literature regarding the existence (), or lack thereof (), of a bulk convective flow of ISF. From previous measurements, we estimate, using equation 7 (below), the root mean squared (rms) displacement of water molecules in the brain to be 1.05 mm in 15 minutes, whereas — for Gd-DTPA molecules — it is 0.68 mm. This reflects the larger molecular weight of Gd-DTPA (938 Da) compared with HO (19 Da). In addition, Gd-DTPA molecules are restricted to the extracellular compartment, whereas HO is able to more freely diffuse into cells. However, comparing the estimated water displacements due to only diffusion, with the actual HO images (), it appears unlikely that diffusion alone explains the rapid brain-wide distribution of HO after 15 minutes, supporting the argument for a bulk convective flow of the ISF. 2 2 2 2 2 2 2 17 17 17 17 17 17 17 Figure 1, C and D Figure 1C Supplemental Video 1 Figure 1C 15 9 8 Figure 2D https://doi.org/10.1172/jci.insight.141159DS1↗
We next exploited the HO tracer to investigate the importance of AQP4 channels to glymphatic transport. The group pretreated with an AQP4 inhibitor (TGN-020, i.p., IC= 3.1 μM) () experienced an 80% ± 10% reduction in HO transported into the parenchyma compared with the vehicle group (). It is clear that, as the AQP4 inhibitor prevents HO penetrating the parenchyma, the concentration of HO builds up in the subarachnoid space and ventricles. Indeed, the distribution of HO in the brains of the AQP4 inhibited group () resembles the distribution of the Gd-DTPA tracer in the control group (). 2 50 2 2 2 2 17 17 17 17 17 26 Figure 2, A and B Figure 2B Figure 1B
The distribution of both tracers within the cerebellum is particularly interesting. In the control group, the Gd-DTPA tracer was slow to penetrate the cerebellum, but by 85 minutes after injection, it did penetrate the outer regions (). However, HO quickly penetrated the whole cerebellum but was more rapidly removed from the white matter (). In the group pretreated with an AQP4 inhibitor (TGN-020), the penetration of HO into the gray matter of the cerebellum was considerably slowed but was strong by 85 minutes, whereas there was still no penetration of the white matter even after 85 minutes. This is presumably related to the high concentration of AQP4 in the cerebellum (). Figure 1, B and D Figure 1B 27 2 2 17 17
Discussion
A serious limitation of previous in vivo glymphatic experiments has been their use of large tracer molecules (MRI tracer Gd-DTPA, MW 938 Da) that cannot be transported by AQP4 water channels in the brain. We hypothesized that the use of these large tracer molecules would result in the underestimation of the rate that CSF water molecules enter the brain parenchyma, thereby underestimating the true flow of interstitial water. By labeling actual water molecules, via highlyO enriched water, we were able to detect glymphatic water transport in vivo using high SNR T2-weighted MRI. The ability of our HO MRI method to directly image water transport in vivo provides an accurate means of studying glymphatic transport. 17 17 2
Reasons for the discrepancy between the transport of HO and Gd-DTPA tracers include the large difference in molecular weight (19 Da for HO and 938 Da for Gd-DTPA) and the presence of the astrocytic water channel AQP4. Molecular tracers such as Gd-DTPA that lack a specific transport pathway (such as ion transporters or channels) are able to reach the parenchyma only through the ~20 nm clefts between overlapping astrocytic end feet. Water, however, is also able to travel through the ICS mediated by highly selected AQP4 water channels, which are expressed in astrocytic end feet that cover the entire vasculature of the CNS. For this reason, it is likely that previous experimental studies based on large tracer molecules may have systematically underestimated subarachnoid CSF water flow into the brain and, thus, underestimated the convective bulk flow of ISF. Asgari et al. () modeled the astrocytic syncytium between CSF and the compartment of the brain interstitium by including, in their model, AQP4 on the plasma membranes, an abundance of AQP4 on the perivascular surfaces, and 20 nm inter–end feet gaps. They demonstrated that the resistance to water flow through extracellular space (ECS) is 2 orders of magnitude larger than through the intracellular space (ICS) of astrocytes. This appears to be a likely explanation for the rapid penetration of the parenchymal that is seen using the HO tracer compared with Gd-DTPA. 2 2 2 17 17 17 28
An alternative hypothesis of CSF physiology has been proposed Oreskovic and Klarica (), in which CSF production and absorption occurs at the level of the capillaries and depends on hydrostatic and osmotic forces. That hypothesis was not supported by our experiments using an AQP4 inhibitor (TGN-020); our experiments clearly demonstrated that the rapid penetration of HO into the parenchyma () is strongly dependent on AQP4, which is absent from the endothelium of brain capillaries. However, it does support the glymphatic hypothesis of Iliff et al. (). 29 Figure 2 15 2 17
Klarica et al. have also suggested a continuous bidirectional mixing of water between the blood, ISF, and the CSF compartments, with no unidirectional net flow of CSF along the CSF spaces (). Again, this is not supported by our HO experiments (). Infusing HO into the cisterna magna leads to a high uptake in the subarachnoid spaces (), and this implies that HO flows largely through the CSF compartment and subarachnoid space. Furthermore, the HO signal change was greatest along the ventral brain surface and along the margins of the pineal recess and the olfactory bulb. These anatomical structures were previously identified by Iliff et al. as the proximal glymphatic transport pathway for CSF within the wider subarachnoid space (). Again, these observations provide support for the glymphatic hypothesis of Iliff et al. (). 30 Figure 1B Figure 2A 15 15 2 2 2 2 17 17 17 17
Similar to natural sleep, general anesthesia has been shown to enhance the transport of CSF tracers (e.g., Gd-DTPA). Studies by the Nedergaard group reported CSF tracer transport was highest under ketamine/xylazine (K/X) anesthesia and lower with α-chloralose, Avertin, or isoflurane (). At least 1 of the effects of anesthesia appears to be increased extracellular volume fraction (), which may explain the enhanced transport of extracellular tracers. However, as the HO tracer can more freely travel through intracellular and extracellular space, we would hypothesize that sleep or general anesthesia may have less influence on its rate of transport. 31 32 2 17
In summary, we have conclusively demonstrated that glymphatic flow imaged using our HO tracer is much more rapid and extensive than when imaged using the Gd-DTPA tracer (). This is strong evidence that the ISF experiences a substantial bulk flow, which can more rapidly clear waste molecules from the parenchyma compared with diffusion alone. Furthermore, we were able to conclusively demonstrate that these glymphatic flows are strongly mediated by AQP4. We believe these advances will not only answer much of the controversy surrounding the glymphatic hypothesis, but will also provide a valuable tool for future investigations into associated neurological disorders. 2 17 Supplemental Video 1
Methods
Study design.
Male Wistar rats (280–300 g, 20–24 weeks old) were obtained from Charles River Laboratories and were randomly assigned to 5 experimental groups. The first study was designed to test our hypothesis that the use of the HO tracer would provide a more accurate measurement of glymphatic transport than Gd-DTPA. Two groups were used (HO [= 7], Gd-DTPA [= 7]). The second study was designed to test the hypothesis that AQP4 channels are critical to glymphatic transport using our HO method (AQP4 inhibitor [= 6], vehicle [= 6], artificial CSF [aCSF;= 3]). After the surgical implantation of an intracisternal cannula, each rat was placed inside the MRI scanner, and MRI was continued for a total of 85 minutes (). In the drug-treated group, TGN-020 was administered i.p. 15 minutes prior to the MRI study. Respiration, blood pressure (BP), and heart rate (HR) were continuously monitored during MRI measurements, and body temperature was maintained at 37.0°C ± 0.5°C (). 2 2 2 17 17 17 n n n n n Figure 1A Supplemental Table 1
Supplemental video.
The video displays differences in the brain-wide distribution of HO (90%O enriched water) and Gd-DTPA (Magnevist, 21 mM) over 85 minutes of recording. The influx of the HO tracer was faster and much more extensive across the whole brain. 2 2 17 17 17
Surgery and physiological monitoring.
All animals were initially anesthetized (5% isoflurane in 30:70 O/NOmixture) in an induction chamber, intubated, and artificially ventilated (with 2%–3% isoflurane in 30:70 O/NOmixture). Body temperature was monitored throughout the experiment, with a rectal thermocouple and maintained at 37°C ± 0.5°C. The femoral artery was cannulated with PE-50 tubing for continuous monitoring of the mean arterial BP and HR (Biopac Systems, MP100) and for the measurement of arterial blood gases (Bayer, Rapidlab 248). The animal was then transferred to a stereotaxic frame, where the head was secured with ear and tooth bars. The head was tilted (45°, snout down), and a midline skin incision was made to expose the dura mater overlying the cisterna magna space. A CM cannula (22-gauge PEEK, SAI, RCMC-01) connected to a closed-end PE10 tube loaded with aCSF (NaCl 140 mmol/L, KCl 3 mmol/L, NaH2PO4 12 mmol/L, NaHCO3 18 mmol/L, CaCl2 2.5 mmol/L, pH 7.4) was advanced 2 mm into the intrathecal space and secured with cyanoacrylate glue to avoid leakage. In the drug treatment group, 200 mg/kg TGN-020 (MilliporeSigma), dissolved in DMSO (10 mM) solution, was administered i.p. (200 mg/kg in 5 mL) 15 minutes before starting the MRI study. Animals were placed prone in a cradle, transferred to the MRI scanner, and monitored and maintained under anesthesia within the physiological range. The CM cannula was connected to a length of PE10 tubing (2.5 cm long, 0.28 mm ID × 0.61 mm OD, Braintree Scientific) filled with the desired tracer (13.5 mM Gd-DOTA, 90% HO or aCSF) and attached to an infusion pump (Graseby 3150 Syringe Pump) for tracer delivery. The infusion pump was evaluated for accuracy by measuring the mass of a target volume of a tracer delivered at the determined flow rate to account for losses and/or bubbles in the infusion system. Respiration, BP, and HR were continuously monitored during MRI measurements, and body temperature was maintained at 37.0°C ± 0.8°C (). At the end of the experiment, the animal was euthanized. 2 2 2 2 2 17 Supplemental Table 1
MRI.
MRI data were acquired using a Bruker PharmaScan 7T/16 cm system controlled by Paravision 5.1 software (Bruker BioSpin) with a gradient coil insert (internal diameter = 90 mm, 300 mT/m) and a 4-channel phased-array surface receive coil used for rat brain imaging. Two different tracers were used for this experiment: Gd-DTPA (Magnevist, 21 mM; MW, 938 Da; Bayer HealthCare Pharmaceuticals Inc.) and HO (90%O-enriched water, NUKEM Isotopes). Furthermore, aCSF was used as a negative control. The scanning protocol for all studies consisted of 3 baseline scans, followed by the intrathecal infusion of tracer via the CM catheter (50 μL at 1.8 μL/min; total time, 28 min). MRI data were continually acquired over a period of 85 minutes. 2 17 17
Gd-DTPA tracer imaged with T1-weighted imaging.
To visualize the glymphatic pathways using gadolinium (T1W shortening effects), 3D T1-weighted FLASH images were acquired in the sagittal plane (repetition time/echo time [TR/TE] = 15 ms/3 ms; flip angle = 15°; number of averages [NA] = 1; field of view [FOV] = 3.0 × 3.0 × 3.0 cm; total scan time = 3 minutes, 5 seconds; acquisition matrix size = 128 × 128, yielding an original image resolution of 0.234 × 0.234 × 0.234 mm). A baseline T1-weighted image was acquired before infusion, which was subtracted from the corresponding images after the infusion of Gd-DTPA. The resulting difference image reflected the distribution of Gd-DTPA in the brain. aCSF was used as a negative control.
HO tracer imaged with high-SNR T2-weighted imaging. 2 17
The quadrupolar moment of theO nucleus reduced the transverse relaxation time, T2, of bondedH nuclei within a water molecule via scalar coupling. This effect was further magnified by the rapid exchange of surroundingHions with HO molecules, 17 1 1 + 17 2
(Equation 1)
According to Meiboom (), for low concentrations of HO in HO, the resultingH transverse relaxation time of water is given by, 19 2 2 17 16 1
(Equation 2)
Where
(Equation 3)
Where P is the fraction of water molecules containingO, τ is the proton exchange lifetime in HO, J is the scalar coupling constant forH-O bond, Tis theH longitudinal relaxation time, and τ is theH exchange lifetime. Hence, the presence of HO water can be detected by a reduction in theH signal on T2-weighted MRI. A key aspect of the MRI sequence employed here was the use of high SNR (>200) T2-weighted MRI scans to enable the detection of small signal changes. A baseline T2-weighted image was acquired before infusion. This was subtracted from corresponding images after the infusion of HO or aCSF (negative control), and the resulting difference image reflected the distribution of HO in the observed brain. T2-weighted images were acquired with a fast spin-echo sequence (rapid acquisition with relaxation enhancement [RARE]) (8 slices, slice thickness = 1.5mm, TR/TE = 3,000 ms/61 ms; NA = 4; FOV = 6.0 × 6.0 × 1.2 cm; total scan time = 3 minutes, 40 seconds; acquisition matrix size = 200 × 200, yielding an original image resolution of 0.3 × 0.3 × 1.5 mm). 17 17 1 17 1 1 17 1 17 17 2 1 2 2 2
MRI data analysis.
MATLAB R2018b (MathWorks) code was developed in-house for postprocessing MRI images. The general postprocessing procedure consisted of brain extraction, head motion correction, and voxel-by-voxel conversion to a percentage signal change. Briefly, a brain mask was created for the removal of nonbrain tissue, improving the performance of the following: using rigid body alignment of each scan to the mean precontrast image, scan-to-scan misregistration caused by head movement was corrected. The resulting registrations were visually inspected to ensure adequate alignment. To ensure that voxel intensity represented a percentage change relative to the average baseline images, all time-series images were subtracted and divided by the baseline average image using the following expression:
For the Gd T1W image study:
(Equation 4)
For the HO/aCSF T2W image study: 2 17
(Equation 5)
Where ΔSis the percent signal change from the baseline (B), I is the time-series image after tracer infusion, and (x,y,z) is the voxel position. The T1-weighted and T2-weighted averaged baseline images were used to anatomically guide the placement of regions of interest (ROIs). In order to compare the T2-weighted signal changes produced by HO with the T1-weighted signal changes produced by Gd-DTPA, both were normalized to the maximum signal change observed over the entire measurement period. Maps of the tracer arrival time, which are defined as the time that the tracer arrives at each voxel after injection, were created on a voxel-by-voxel basis. % 2 17
The 1-dimensional Einstein diffusion equation can be solved to give the probability density function,, for particles undergoing Brownian motion: P(x,t)
(Equation 6)
Where D is the diffusion coefficient andis the particle displacement from its initial position, during an observation time, t. For this Gaussian function the root mean square displacement,, of an ensemble of particles is given by: x x 0 x rms
(Equation 7)
From previous diffusion MRI experiments using a long observation time (600 ms) (), we measured the plateau diffusion coefficient of water within the rodent cortex as 0.00062 mm/s. From equation 7, the root mean square displacement, x, of water for experiment times of 3.85 minutes (231 sec), 7.65 minutes (459 sec), 11.08 minutes (665 sec), and 14.83 minutes (893 sec) was calculated to be 0.53 mm, 0.75 mm, 0.91 mm, and 1.05 mm, respectively. Previous phantom studies of Gd-DTPA diffusivity in a tissue mimic material (polyvinyl alcohol-cryogel) have been performed using MRI, giving a diffusion coefficient of 0.00026 mm/s (). From equation 2, the root mean square displacement, x, of Gd-DTPA for experiment times of 7.65 minutes (459 sec), 11.08 minutes (665 sec), and 14.83 minutes (893 sec) were calculated to be 0.49 mm,0.59 mm, and 0.68 mm, respectively. 33 34 2 2 rms rms
Statistics.
No statistical methods were used to predetermine sample size, with sample sizes being similar to those reported in a previous study (). Statistical tests performed using GraphPad Prism (GraphPad Prism Software). Statistical comparisons between groups were performed by repeated measures 2-way ANOVA, followed by Sidak’s test to correct for multiple comparisons. All tests were considered statistically significant for< 0.05. All data are presented as mean ± SD unless otherwise stated. 15 P
Study approval.
Experiments were carried out under license from the UK Home Office in accordance with the Animals (Scientific Procedures) Act, 1986, incorporating European Directive 2010/63/EU and approved by the University of Glasgow Ethical Review Panel.
Author contributions.
MSA and WMH developed the MRI method, contributed to conceiving and designing experiments, analyzed all data, produced the figures, and contributed to discussion and interpretation of results and the writing of the final manuscript. LG contributed to the design of the experiments, performing surgical procedures on rats and acquiring the MRI data. LG also contributed to writing the manuscript. LMW was the secondary Ph.D. supervisor of MSA and contributed to the design of the experiments, to discussion and interpretation of results, and to the writing of the final manuscript.