Theranostics 2021; 11(13):6240-6250. doi:10.7150/thno.55953 This issue Cite
Research Paper
1. Department of Electronic Science, Fujian Provincial Key Laboratory of Plasma and Magnetic Resonance, School of Electronic Science and Engineering, National Model Microelectronics College, Xiamen University, Xiamen, Fujian, China.
2. Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins University School of Medicine, Baltimore, MD, USA.
3. The First Affiliated Hospital of Jinan University, Guangzhou, Guangdong, China.
4. Department of Diagnostic Radiology and Nuclear Medicine, University of Maryland, Baltimore, MD, USA.
5. F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, MD, USA.
6. Department of Neurosurgery, Johns Hopkins University, Baltimore, MD, USA.
7. Department of Neurology, Johns Hopkins University, Baltimore, MD, USA.
* These authors contributed equally.
Rationale: Endovascular intervention plays an important role in the treatment of various diseases, in which MRI-guidance can potentially improve precision. However, the clinical applications of currently available contrast media, including Gadolinium-based contrast agents and superparamagnetic iron oxide particles (SPIO), are hindered by safety concerns. In the present study, we sought to develop D2O as a novel contrast agent for guiding endovascular neurointervention.
Methods: Animal studies were approved by institutional ACUC and conducted using an 11.7 T Bruker Biospec system and a 3T Siemens Trio clinical scanner for rodent and canine imaging, respectively. The locally selective blood brain barrier opening (BBBO) in rat brains was obtained by intraarterial (IA) injection of mannitol. The dynamic T2w* EPI MRI sequence was used to study the trans-catheter perfusion territory by IA administered SPIO before mannitol administration, whereas a dynamic T1w FLASH sequence was used to acquire Gd contrast-enhanced MRI for assessing BBBO after injection of mannitol. The contrast generated by D2O assessed by either EPI or FLASH methods was compared with the corresponding results assessed by SPIO or Gd. The utility of D2O MRI was also demonstrated to guide drug delivery to glioma in a mouse model. Finally, the clinical utility of D2O-MRI was demonstrated in a canine model.
Results: Our study has shown that the contrast generated by D2O can be used to precisely delineate trans-catheter perfusion territory in both small and large animals. The perfusion territories determined by D2O-MRI show moderate correlation with those by SPIO-MRI (Spearman coefficient r = 0.5234, P < 0.001). Moreover, our results show that the perfusion territory determined by D2O-MRI can successfully predict the areas with BBBO after mannitol treatment similar to that assessed by Gd-MRI (Spearman coefficient r = 0.6923, P < 0.001). Using D2O-MRI as imaging guidance, the optimal infusion rate in the mouse brain was determined to be 150 µL/min to maximize the delivery efficacy to the tumor without serious off-target delivery to the brain parenchyma. The enhanced drug delivery of antibodies to the brain tumor was confirmed by fluorescence imaging.
Conclusion: Our study demonstrated that D2O can be used as a negative MRI contrast medium to guide endovascular neurointervention. The established D2O -MRI method is safe and quantitative, without the concern of contrast accumulation. These qualities make it an attempting approach for a variety of endovascular procedures.
Keywords: deuterium oxide, MRI contrast medium, intra-arterial hyperosmotic blood brain barrier (BBB) opening, endovascular neurointervention, MRI guidance.
Intra-arterial (IA) drug delivery has recently become an important branch of endovascular neurointervention for treating a variety of central nervous system (CNS) diseases [1]. Direct endovascular infusion through a catheter placed in the carotid or vertebral arteries can deliver drugs at high, virtually undiluted, local concentration to the targeted sites [2]. Moreover, IA infusion can be combined with blood-brain barrier (BBB) opening strategies to further improve the drug delivery efficiency in the regions with intact BBB. For example, Charest et al. reported that, compared with IA infusion alone, transiently opening the BBB with mannitol could result in an overall 2-fold higher carboplatin accumulation in experimental brain tumors (an 80-time increase in the nucleus of tumor cells) [3].
Real-time imaging guidance is crucial for accurately visualizing target vessels during catheter advancement and prediction of the affected tissue region, hence warranting the success of endovascular neurointervention. While fluoroscopy is still the standard of practice, increasing interests have been focused on MRI-based guidance because MRI has superb soft-tissue contrast and no ionizing radiation, and can provide morphological and functional assessments that are critical for treatment planning and monitoring. Recent developments in hardware and software have provided modern MRI with a satisfactory spatial and temporal resolution, thereby making MRI a better interactive tool for guiding endovascular intervention in real-time. Very often, contrast media, including Gadolinium-based contrast agents (GBCA) [4], superparamagnetic iron oxide particles (SPIO) [5] and carbon dioxide [6], are used to improve the conspicuity and thereby ensure the accuracy of intervention. However, the widespread clinical use of these agents is hampered by safety concerns.
In the present study, we sought to develop deuterium oxide (D2O) as a novel contrast agent to guide endovascular neurointervention. D2O is a stable isotopic form of water with similar physical and chemical properties as regular water. Previous studies demonstrated that D2O can be used as a diffusible tracer to measure blood flow and tissue perfusion [7]. While traditionally D2O is only visible using 2H-tuned MRI hardware, the negative contrast of D2O in the 1H MRI caused by the H2O-D2O replacement effect (also known as proton replace effect) has been used to noninvasively study water dynamics in Wood Xylem [8] and precisely measure the relaxation times of fat in excised breast tissue samples[9]. Recently, Wang et al. also showed the utility of D2O as a negative contrast medium in 1H-MRI for assessing cerebral perfusion in rat brains [10]. Inspired by these studies, we rationalize that the D2O-based contrast can be utilized to track the spatial distribution of endovascularly infused drug solutions. Transcatheter infusion of D2O, even in a small quantity, is expected to generate sufficiently high contrast as a high local concentration of D2O is present in the arteries distal to the catheter.
All chemicals were purchased from Sigma Aldrich (Saint Louis, MO, USA) unless otherwise stated. For the phantom experiment, two sets of samples were prepared in either aqueous solutions or 2% agarose gel. The D2O (99.9%) was mixed with PBS (pH = 7.4) at the final concentrations (v/v) of 5%, 10%, 20%, 40%, 60% and 80%. All samples were prepared freshly and transferred to 5 mm NMR tubes prior to MRI measurement.
All procedures were approved by our institutional Animal Care and Use Committee and in accordance with guidelines for the care and use of laboratory animals. Sprague-Dawley rats (male, 200-250 g, Charles River, n = 6) were used to demonstrate the ability D2O-MRI to monitor the cerebral IA injection by a means of carotid artery catheterization. The surgical procedures were performed as previously described [11]. Briefly, the animals were anesthetized under 2% isoflurane and positioned supine. The common carotid artery (CCA) bifurcation was dissected, and the occipital artery was coagulated. The external carotid artery (ECA) and pterygopalatine artery (PPA) were temporarily ligated with 4-0 silk sutures. A small arteriotomy in the CCA was made and a catheter (VAH-PU-C20, Instech Solomon Inc.) connected to #30 PTFE tubing was introduced into the internal carotid artery (ICA) via the arteriotomy. The intra-arterial catheter was secured with the animal during MRI scans.
The same catheterization procedures were performed for cerebral IA injection in brain tumor-bearing mice. In brief, six C57BL6 mice (female, 5-6 weeks) were stereotactically injected with 5 x 104 murine glioma CT2A cells. The injection coordinates were -2.0 mm for the anterior-posterior (from Bregma), 1.5 mm for medial-lateral (from midline), and -3.0 mm for dorsal-ventral (from the surface of skull) axes, respectively. Orthotopic brain tumors were allowed to grow ~ 7 days to reach the desirable size (~2-3 mm in diameter). Prior to MRI, three mice were randomly selected to, under anesthesia, put a microcatheter (PE-8-100, SAI Infusion Technologies) in proximal ICA via a small arteriotomy from CCA. After the intra-arterial catheter was secured using two purse-string suture ties applied around CCA, mice were transferred to the MRI scanner. The rest of the mice (n = 3) were injected with 0.6 mL IgG-FITC (1.5 mg/mL) via tail vein and euthanized two hours later to collect brains for ex vivo fluorescence imaging.
For cerebral IA injection into the anterior cerebral artery in a canine (male greyhound dog, weight ~ 25 kg, n = 1), a 5-French femoral arterial sheath was surgically introduced. A 5-French Glide catheter was then advanced over a 0.035-inch guide wire and the right anterior cerebral artery was selectively catheterized using a 1.7-French microcatheter was then advanced over a 0.014-inch microwire under X-ray angiography roadmap guidance.
In vitro MRI was performed on an 11.7 T Bruker Avance system (Bruker, Ettlingen, Germany) equipped with a 15 mm sawtooth RF coil. All measurements were conducted at room temperature (20 oC). T1 and T2 maps were acquired according to our previously published procedures [12]. In brief, longitudinal (T1) relaxation times of the samples were assessed using a Rapid Acquisition with Relaxation Enhancement (RARE) -based saturation recovery sequence with eight TR times ranging between 200 ms to 15,000 ms (TE = 4.3 ms and RARE factor = 4, central encoding). T2 maps were acquired using a modified RARE pulse sequence (TR/TE = 25000/4.3 ms and RARE factor = 16) with a Carr-Purcell-Meiboom-Gill (CPMG) T2 preparation module. The T2 preparation period consisted of an element of CPMG pulse train with tCPMG = 10 ms, and a total of 16 CPMG loops (2 to 1024) were used, corresponding to echo times = 20 ms to 10.24 sec. Other imaging parameters were: 256x256 acquisition matrix with a spatial resolution of c.a. 50x50 µm2, and slice thickness of 1 mm.
Rat and mouse MRI was performed on a horizontal bore 11.7 T Bruker Biospec system using a 72 mm quadrature volume resonator as a transmitter, and a four-element (2×2) phased array coil as a receiver. The timeline of in vivo MRI study is illustrated in Figure 1. Briefly, the perfusion territory was assessed by both SPIO- and D2O- MRI. To precisely determine the perfusion territory, SPIO and D2O were infused at the same flow rate as mannitol. Then, BBB opening (BBBO) was induced by infusion of 1 mL mannitol (25%, Hospira) via the IA route at the speed of 0.6 mL/min. Five minutes after mannitol injection, BBBO was assessed by Gd-based DCE MRI.
For SPIO-MRI, 300 µL ferumoxytol (0.3 mg/mL) was injected through an IA catheter, and dynamic T2w* images were acquired using a multi-slice echo planar imaging (EPI) sequence with TR/TE = 2000/9.7ms, segment factor = 2, bandwidth = 455kHz, slice number = 15, slice thickness = 1mm, a matrix size of 128×128 within a FOV of 3×3 cm2 [13].
For Gd-MRI, a bolus of Prohance (0.1 mmol/kg) was injected intraperitoneally, and dynamic T1w images were acquired using a FLASH sequence with TR/TE = 100 ms/3 ms, flip angle = 25°, average number = 1, slice number = 11, slice thickness = 1 mm, a matrix size of 128×128 within a FOV of 3×3 cm2. The acquisition time for each scan is 19 sec.
D2O was IA injected at the rate of 0.6 mL/min and the total injection volume was 0.4 mL and D2O-MRI was acquired using the same EPI (n = 3) or FLASH sequences (n = 3), which would permit the direct comparison of D2O MRI measures with those of SPIO- or Gd-CE MRI, respectively.
Schematic of MRI studies of intracarotid infusion of mannitol in the rat brain. A) A catheter was inserted into the ICA for IA injection of SPIO, D2O and mannitol. B) Timeline of MRI study. The rat was sequentially scanned after IA injection of SPIO and then D2O. After the mannitol injection, the animal was removed from the scanner to receive intra-peritoneal (i.p.) injection of ProHance (Gadoteridol, 0.1 mmol/kg) and T1w images were acquired.
The same MRI methods were used for mouse imaging (n = 3), except a series of bolus injections (duration = 20 s) of D2O at different infusion rates of 50, 100, 150, or 200 µL/min (total injection volume = 16.7- 66.7 µL) was investigated for the optimal perfusion to the tumor region. At each infusion rate, seven D2O MR images were acquired starting at 40 s before injection till 74 s after the stop of injection (total acquisition time= 2 min 14 s). Using the optimal rate determined by D2O MRI, a mixture (0.6 mL) of 25% mannitol and IgG-FITC (1.5 mg/mL) was then administered. Two hours later, mice were sacrificed, and brains were harvested and fixed in 4% paraformaldehyde solution for ex vivo fluorescence imaging.
Dog MR imaging was performed on a clinical 3T Siemens Trio using a quadrature head coil. The dynamic T2*w images during IA D2O and SPIO (2mL/min) were acquired using a single shot GE-EPI (TE = 36 ms, TR = 3000 ms, FOV = 180x180 mm2, matrix size = 192x192, slice thickness= 2.5 mm, and acquisition time = 3 s and 50-100 repetitions). Two mL of 25% mannitol over 1 min was administrated thereafter at the same rate as D2O injection. Five minutes after IA mannitol, a bolus of Magnevist (gadolinium, 0.125 mmol/kg) was injected intravenously. Standard Gd-CE (Gd-contrast enhancement) images were acquired then.
The uptake of IgG-FITC in the fixed mouse brains was measured using a Spectrum/CT IVIS® in vivo imaging system (Ex/Em = 500/540 nm) and analyzed using the Living Image® processing software from the manufacturer (PerkinElmer, Waltham, MA).
All experiment data were processed by either custom-written scripts in MATLAB (R2020a, Mathworks, Natick, MA, USA) or ImageJ (version 1.51, NIH, https://imagej.nih.gov/ij/, Bethesda, MD, USA). For in vitro data, using ROI masks that were manually drawn on a T2w image, mean ROI values were calculated and plotted. The proton density (PD) of each sample was estimated using the mean ROI value of the final image of the saturation recovery MRI, where TR = 15 sec and effective TE = 4.3 ms.
T1 relaxation times were estimated by fitting the ROI values to Eq. 1,
Eq. 1
where S(TR) is the MRI signal at each TR time, and the theoretical maximal MRI signal S0, T1 times, and constant C are the parameters to be estimated.
The R1 relaxivity (r1) of each compound was subsequently estimated by fitting water proton relaxation rates (R1 = 1/T1) as a linear function of D2O concentration (Eq. 2).
Eq. 2
where R10 is the water intrinsic R1 relaxation rate, a fitting parameter in our study.
T2 relaxation times were estimated by fitting ROI values to Eq. 3,
Eq. 3
where S(TE) are the MRI signal at each CPMG TE time, and the theoretical maximal MRI signal S0, T2 time, and constant C are the parameters to be estimated.
The R2 relaxivity (r2) of each compound was subsequently estimated by fitting R2 relaxation rates (R2 = 1/T2) as a linear function of D2O concentration (Eq. 2).
Eq. 4
where R20 is the water intrinsic R2 relaxation rate, a fitting parameter in our study.
For in vivo data, we calculated the pixel-wise maximum signal change map (∆SImax) after injection of SPIO or D2O using the Matlab function max ([SIpost(t)- SIpre]/SIpre). The Gd-contrast-enhanced (CE) map was calculated by comparing the post-contrast image at 9 min with that of before injection.
The area affected by the first pass of transcatheter injection, namely perfusion territory, was delineated using a histogram-based analysis as described previously for SPIO [14] and Gd agents[13]. In brief, for each animal, a histogram was plotted to represent the frequency distribution (percentage of voxels, y-axis) of the level of signal change (x-axis) within the entire brain, which was then fitted into a 2-Gaussian distribution using custom-written Matlab scripts. A cut-off of signal change was chosen at where the overlap between two Gaussian functions was minimum. Pixels with signal changes larger than the cut-off were assigned as the perfusion territory.
All statistical analyses were conducted using Matlab and GraphPad Prism (version 7, GraphPad Inc, San Diego, USA). The contrast maps generated by D2O, SPIO, and Gd-DOTA were compared using the non-parametric Spearman's rank-order correlation analysis. Non-parametric tests were used because the normality of data distribution was not passed the D'agostino-Pearson test. The Student's t-test was used to compare the D2O contrast changes in the tumor at different infusion rates and the difference between the mean fluorescence intensities in the mouse brains administered using MRI-guided IA injection and those by IA injection. In all these analyses, p-value < 0.05 was regarded as statistically significant.
We first characterized the MRI contrast of D2O in aqueous solution and 2% agarose gel. Figure 2A shows a representative MR image of D2O agarose gel samples at concentrations from 0 to 0.8 (D2O : H2O volume fraction) that were acquired using an 11.7T preclinical scanner. The negative effect of D2O on proton density (PD) contrast is linearly proportional to its concentration (R2 = 0.9952 and 0.9963, in solution and gel respectively, Figure 2B), attributed to the proton replacement effect. D2O also has relatively small but observable effects on water T1 and T2 contrast. At the field strength used in our study (i.e., 11.7T), the r1 relaxivities per D2O fraction (v/v) were estimated to be -0.230 s-1 (R2 = 0.9998) and -0.236 s-1 (R2 = 0.9993) in solution and gel, respectively (Figure 2C). The r2 relaxivities per D2O fraction (v/v) were estimated to be -0.211 s-1 (R2 = 0.8406) and -9.28 s-1 (R2 = 0.9956) in solution and gel, respectively (Figure 2D). The negative sign of relaxivity denotes the tendency of D2O to decrease the R1 and R2 relaxation rates (increase T1 and T2 relaxation times) of water, as opposed to commonly used MRI relaxation contrast agents. It should be noted that, as indicated by the simulations (Figure S1, Supporting Information), the T1/T2 effects are not negligible and, depending TR/TE parameters, may cause errors if simply interpreting the signal decrease (∆SI) as the concentration of D2O.
MRI of deuterated aqueous PBS solutions and agarose (2%). A) Representative proton densities (PD)-weighted image (TR/TE = 15,000/4.3 ms) of D2O-H2O samples in 2% agarose gel, where a small capillary (OD = 1 mm) containing PBS was used as a marker for sample position. B) Proton densities, C) R1 relaxation rates, and D) R2 relaxation rates of aqueous (blue) and gel (red) samples prepared at different D2O: H2O ratios.
As shown in Fig. 3 and the supp. video S1, IA infusing of D2O (injection rate = 0.6 mL/min, total volume = 0.4 mL) to the rat brain resulted in a rapid signal decrease in the ipsilateral brain hemisphere (red ROI), attributed to the proton replacement effect of locally high concentration of D2O. As soon as the infusion was stopped, MRI signal recovered exponentially, reflecting quick washout of D2O from the tissue [15]. At the end of the MRI scan, ~ 4 minutes after the stop of infusion, the MRI signal recovered back to 95.3% of its baseline level. In comparison, a much slower signal reduction was observed in the contralateral hemisphere (Figures 3A, B, blue ROI), i.e., ∆SI ~ 4.0% by the end of MRI scan, attributed to circulating D2O in the bloodstream. The signal changes of the two ROIs reached equilibrium approximately at 2 minutes after the infusion was stopped. It should be noted that in certain regions (Figure 3A, yellow arrow), D2O clearance was quite slow, i.e., remaining ∆SI = 14.5% at the end of the study, indicating that a significant amount of D2O entered and remained in the interstitial and intracellular spaces in these regions. Using the ∆SI% at the time of infusion stop, we can determine the perfusion territory that the infused solution can reach (Figure 3C). It should be noted that due to the T1/T2 relaxation effects as shown in Figure S1, ∆SI is not necessarily the same as D2O concentration. Indeed, assuming that majority of D2O still remained in the blood (where T10/T20 = 2800/46 ms [16]), with the acquisition parameters used to acquire Figure 3 (TR/TE = 100/3 ms), the concentrations of D2O to generate 4.0% (contralateral) and 4.7% (ipsilateral) signal decrease are estimated to be 2.5% and 2.9%, respectively.
Figure 4 shows the comparison between the perfusion territories determined by SPIO- and D2O-MRI in a representative rat (Figures 4B,D), revealing similar contrast enhanced regions. This is further confirmed by comparing the maximal contrast-enhancement maps (Figure 4C) acquired by the two methods, with a good pixel-wise correlation (Spearman coefficient = 0.5234, P < 0.0001, Figure 4F) and good overlap between the contrast-enhanced regions. As shown in Figure S2, all animals showed good agreement between the territories determined by the two methods (r = 0.5692 and 0.5585, for the other two rats), which implies that, even there are fundamental differences between D2O and SPIO nanoparticles, D2O can also be used as an alternative to SPIO for delineating the perfusion territory of endovascular injection.
We then investigated the ability of D2O-MRI to predict the effective regions in response to a transcatheter infusion procedure, in particular, in the present study, the regions where BBB will be opened by mannitol. As depicted in Figure 1, D2O-MRI was performed before mannitol infusion, and Gd-enhanced MRI was carried out after mannitol treatment to assess the resulted BBBO. The comparison of Gd-MRI and D2O-MRI is shown in Figures 5B-E, revealing the contrast-enhanced areas by the two methods were in good agreement, i.e., Spearman coefficient = 0.6923 (P < 0.0001, Figure 5F). As shown in Figure S2, a good correlation was found between the territories determined by the two methods in all animals (r = 0.6558 and 0.5817, for the other two rats). Given Gd-CE MRI is a commonly used method to assess BBBO after intervention, our results imply that, even there are fundamental differences between D2O and Gd agents, D2O-MRI may also be used as the alternative to Gd-CE MRI by predicting the BBBO reliably before the actual injection of mannitol.
Dynamic D2O-enhanced MRI in a representative rat brain after intraarterial injection of 0.6 mL D2O. A) MR images before (t = 0), during (t = 60) and post-injection (t = 240 sec). B) Dynamic signal changes in two representative ROIs that were manually drawn in the ipsilateral (red) and contralateral (blue) hemispheres. C) 3D visualization of the perfusion territory (gold).
Comparison of the perfusion regions measured by SPIO-MRI and D2O-MRI. A) T2w reference image. B) T2*w image and C) Maximum contrast-enhancement map by SPIO infusion. D) T2*w image and E) Maximum contrast-enhancement map by D2O infusion. F) Correlation between SPIO- and D2O- enhanced MRI. The Spearman correlation coefficient (r) was determined to be 0.5234.
Comparison of the BBB opening regions predicted by D2O-MRI and the BBB-disrupted region measured by Gd-CE MRI. A) T2w reference image. B) T1w image and c) Maximum contrast-enhancement map by Gd-CE MRI. D) T1w image and E) Maximum contrast-enhancement map by D2O. F) Correlation between the two measures. The Spearman correlation coefficient (r) was determined to be 0.6923.
D2O-MRI guided neurointervention drug delivery to a brain tumor in a representative mouse. A) T2w image showing the location of the tumor (yellow arrow). Pre-D2O images, post-D2O images, contrast maps (∆SI = Spost- Spre), and maps of relative signal change (%) with D2O injected at 50 µL/min (B), 100 µL/min (C), 150 µL/min (D), and 200 µL/min (E), respectively. F) Bar plots showing the comparison between the mean MRI signal changes in tumor regions of three mice using two different flow rates (P = 0.0382, paired Student's t-test). G) ex vivo fluorescence imaging of IgG-FITC accumulated in the mouse brains with the IgG injected either by D2O-guided IA injection or by IV injection (P = 0.0023, unpaired Student's t-test).
We then demonstrated the ability of D2O MRI to guide interventional drug delivery in an orthotopic mouse brain tumor model. For transcatheter injection, besides the location of the catheter, injection rate is another important parameter determining the area that the injected agent can reach. Our results showed that, when injected at slow rates (i.e., 50 and 100 µL/min), only a small portion of the tumor could be perfused as revealed by the D2O MRI (Figures 6B,C), indicating a low delivery efficiency. Increasing infusion rate significantly increased the delivery efficiency to the tumor (Figures 6D,E). However, when the rate was too high (200 µL/min, Figure 6E), overwhelming perfusion in the brain parenchyma was observed, which may cause unwanted adverse effects. Using D2O-MRI as imaging guidance, we determined that 150 µL/min to be the optimal rate for drug delivery to the tumors without serious off-target delivery. D2O MRI revealed that all animals (n = 3) showed a significantly augmented perfusion in the tumor region as compared to that by a low rate (P = 0.0382, paired Student's t-test, Figure 6F) when this optimal rate was used. The optimized drug delivery was confirmed by fluorescence imaging using fluorescent IgG-FITC (resembling antibody drugs, which cannot permeate the intact BBB). The quantitative analysis of ex vivo fluorescence imaging of the excised mouse brains showed a two-fold higher drug delivery efficiency (as quantified as total radiant efficiency (photons/s)/(μW/cm2)) was achieved to deliver antibody IgG to the tumor-bearing brain compared to the traditional intravenous (IV) injection route (P = 0.0023, unpaired Student's t-test).
Finally, we performed a 3T D2O-MRI in a canine model for predicting interventional BBBO areas (Figures 7 and S3). As can be seen in Figures 7B-C, both SPIO and D2O MRI methods were able to assess the perfusion territories, but the regions revealed by D2O MRI appeared larger than SPIO MRI. A moderate correlation (Spearman r = 0.5350, P < 0.0001) was found between the hypointense regions caused by D2O and SPIO (Figure 6E). Following the infusion of mannitol, areas of BBBO were visualized by Gd-CE MRI (Figure 6D). The hypointense area of D2O correlated well with the hyperintense area shown in Gd-CE MRI (Spearman r = 0.6009, P < 0.0001, Figure 6F), indicating D2O MRI can be used to guide interventional BBBO conducted on 3T MRI scanners. The predictive power of D2O MRI was even better than SPIO MRI, which only showed a relatively poor correlation with Gd-CE MRI (Spearman r = 0.4365, P < 0.0001). This pilot 3T study demonstrates the feasibility of D2O-MRI on clinical MRI scanners and future studies are warranted to fully test and validate the clinical utility of D2O MRI.
Our study demonstrated that the negative contrast generated by D2O can be used as a useful contrast medium to assess the spatial and temporal distribution of the injected solution in the brain, providing a new means to guide endovascular neurointervention. D2O-MRI has the capacity for accurately determining the perfusion territory and predicting the affected areas and treatment responses, allowing for real-time guiding endovascular interventions such as IA drug delivery targeted to a specific region in the brain.
3T D2O MRI guidance for interventionally hyperosmotic BBBO in a dog. A) T2w image pre-injection. B) D2O MRI contrast-enhancement map. C) SPIO MRI contrast-enhancement map. D) Gd-contrast-enhancement map after the administration of mannitol. E) Correlation between SPIO- and D2O- MRI. The Spearman coefficient (r) = 05350. F) Correlation between Gd- and D2O- MRI. The Spearman coefficient (r) = 0.6009.
Deuterium is a non-radioactive isotope of hydrogen with a natural abundance of 0.0156%. As a stable and enriched isotopic form of 1H water, D2O has been used in clinical tests and treatments with few adverse effects [17]. For example, deuterium-dilution technique is the gold standard to measure total body water space [18] and body composition [19], which is even safe for children. Numerous studies showed the low toxicity of D2O even at the level of 20% of total body fluids over a short period [2, 20]. Unlike previous studies required the i.v. injection of D2O at a relatively high dose (e.g., 20 mL/kg [21]), our study employed a relatively low dose (i.e., 1.6 mL/kg). Indeed, this dose level only accounts for ~3.2% plasma concentration (total injection volume of D2O = 0.4 mL, assuming the total blood volume of rats is ~12.5 mL). Considering the high contrast level achieved in the present study, the injection volume can be further decreased substantially. Hence, D2O-MRI can be considered as a highly safe approach for accomplishing image-guided endovascular neurointervention.
D2O can be directly detected by 2H NMR or MRI, an approach explored for measuring blood flow and perfusion in normal organs [21-24] and tumors [15, 23, 25, 26]. However, its wide application is hindered by low sensitivity (i.e., ~6.5 times lower than 1H) and the need for special hardware. In 1H NMR or MRI, D2O is invisible and will reduce the detectable MR signal when mixing with water in a concentration-dependent manner. While D2O produces 1H MRI contrast mainly through changing the proton density of water by the so-called proton replacement effect, it also has observable effects on the T1 and T2 relaxation times of water. Of note, the r2 of D2O in agarose gel (I.e., -9.28 s-1) was found much larger than that in aqueous solution (i.e., -0.21 s-1), whereas the r1 of D2O was largely unaffected (-0.230 and -0.236 s-1 in solution and gel, respectively). The much stronger T2 effect of D2O in the agarose gel can be explained by the fact that its effect to weaken the dipole-dipole interaction between protons would become more pronounced when the mobility of water molecules becomes slower due to either a higher amount of bound water on the surface of macromolecules [27] or a slow diffusion rate [28]. Our observation is consistent with previous studies [28, 29]. For example, Zhong et al. [29] reported the r2 relaxivities of D2O (at 7T) of -0.2 s-1 and -5.3 s-1 in pure solution and rat liver, respectively.
The negative MRI contrast of D2O was used to measure water permeability across cell membranes in vitro [30] and cerebral blood flow in vivo [10]. For the later application, D2O was injected through the tail-vein at a dose of 2 mL/100g body weight, and approximately 5% MRI signal decrease was generated in the brain [10]. In comparison, we explored the application of D2O as a contrast medium for guiding IA drug delivery, in which D2O was injected directly into the blood vessels of interest. As a result, the contrast was remarkably high (> 60% in the areas close to the catheter), which is sufficient for rendering the perfusion territory reliably. Indeed, in the present study, we have successfully demonstrated the utility of D2O MRI in improving the delivery of antibodies to the targeted brain tumor using an image-guided neurointervention approach. Based on our findings, we foresee the great potential of D2O-MRI in endovascular neurointervention.
During an endovascular intervention, the affected areas are largely determined by the transcatheter perfusion territory, which is determined by the location of catheter and infusion rate. The assessment of perfusion territory is crucial for the precision and effectiveness of an intravascular intervention. SPIO has been previously used to detect the perfusion territory in the brain and to guide IA injection [31]. Our present study revealed a good correlation between the perfusion territory determined by SPIO-MRI and that by D2O-MRI. However, it should be noted that perfusion territories by the two methods are only moderately correlated. The discrepancy is not a surprise as D2O is a small molecule and to which BBB is permeable (via an active water transport mechanism [32]), whereas SPIO is a nanoparticle and to which intact BBB is completely impermeable. Moreover, SPIO particles have a well-known spatial “blooming effect” in the T2*w images, which may lead to overestimation of perfusion area. This effect is caused by the quick diphase of diffusing water molecules by the SPIO-generated local B0 disturbances/inhomogeneities [33]. As D2O behaves just like a water molecule and causes no B0 distortion, it would not cause a blooming effect. The perfusion region measured by D2O is hence more accurate. Indeed, our 3T study on a canine model revealed that D2O-MRI provided a better prediction of the BBBO region than SPIO-MRI, suggesting D2O is an excellent contrast medium for rendering the perfusion territory and guiding IA drug delivery.
One potential limitation of using D2O MRI to guide endovascular intervention is that D2O is a diffusible tracer that won't restrictively stay in the bloodstream after injection. Hence, the accurate determination of its concentration in the intravascular (or extravascular) compartment will require compartmental pharmacokinetic analysis [15, 34]. Another potential limitation is the slow washout kinetics of D2O in some regions. This can prevent D2O from being used repeatedly in a short time interval. Reducing the injection volume may partially overcome this challenge, and further studies to optimize the injection volume are warranted. Similar to SPIO, D2O generates a negative contrast on MR images. Generally speaking, positive contrast (e.g., Gd agents) is more preferable because the extent of contrast enhancement can be higher than 100% by a positive contrast agent, whereas maximum change by a negative contrast agent (e.g., SPIO and D2O) is 100%. Moreover, the darkened regions by a negative contrast agent may not be distinguishable from voids in the image caused by other effects, such as image artifacts or intrinsic hypointense tissues. Fortunately, as demonstrated in the present study, when a dynamic imaging scheme is used, the post-contrast images are compared with the baseline image, allowing easy discrimination between the darkened area generated by D2O easily and intrinsic hypointense regions. Finally, it should be noted that the quantification of D2O could be confounded by its effects on water T1 and T2 relaxation times. Depending on the TR and TE times and intrinsic tissue T1 and T2 times, the MRI signal change (contrast) caused by D2O may deviate significantly from the pure proton replacement effect. The T2 effect is even more pronounced when D2O enters tissues than in the blood. Therefore, the accurate quantification of D2O requires taking the relaxation effects into account. If possible, long TR and short TE should be used to reduce the errors.
In summary, our study showed that D2O is a negative MRI contrast medium suitable for guiding endovascular neurointervention. Using a high field 11.7T and low field 3T 1H MRI scanners, we demonstrated the ability of D2O MRI to accurately visualize perfusion territory and predict BBBO regions in different animal models. The D2O-MRI is a safe approach with a good quantitative ability, and hence may be useful for a variety of endovascular procedures.
BBB: blood brain barrier; D2O: deuterium oxide; IA: intra-arterial; BBBO: BBB opening; ROI: region of interest; SPIO: superparamagnetic iron oxide particles.
Supplementary figures.
Supplementary video/movie.
This work was supported by NIH (R01CA211087, R01NS091110, and R21NS106436), Maryland Stem Cell Research Fund (2019-MSCRFD-5031), and National Multiple Sclerosis Society (PP-1908-34973).
The authors have declared that no competing interest exists.
1. Klopp CT, Alford TC, Bateman J, Berry GN, Winship T. Fractionated intra-arterial cancer; chemotherapy with methyl bis amine hydrochloride; a preliminary report. Ann Surg. 1950;132:811-32
2. Joshi S, Meyers PM, Ornstein E. Intracarotid delivery of drugs: the potential and the pitfalls. Anesthesiology. 2008;109:543-64
3. Charest G, Sanche L, Fortin D, Mathieu D, Paquette B. Optimization of the route of platinum drugs administration to optimize the concomitant treatment with radiotherapy for glioblastoma implanted in the Fischer rat brain. J Neurooncol. 2013;115:365-73
4. Frayne R, Omary RA, Unal O, Strother CM. Determination of optimal injection parameters for intraarterial gadolinium-enhanced MR angiography. J Vasc Interv Radiol. 2000;11:1277-84
5. Wacker FK, Reither K, Ebert W, Wendt M, Lewin JS, Wolf K-J. MR Image-guided Endovascular Procedures with the Ultrasmall Superparamagnetic Iron Oxide SH U 555 C as an Intravascular Contrast Agent: Study in Pigs. Radiology. 2003;226:459-64
6. Wacker FK, Maes RM, Jesberger JA, Nour SG, Duerk JL, Lewin JS. MR imaging-guided vascular procedures using CO2 as a contrast agent. AJR Am J Roentgenol. 2003;181:485-9
7. Kim SG, Ackerman JJ. Multicompartment analysis of blood flow and tissue perfusion employing D2O as a freely diffusible tracer: a novel deuterium NMR technique demonstrated via application with murine RIF-1 tumors. Magnetic resonance in medicine. 1988;8:410-26
8. Ilvonen K, Palva L, Peramaki M, Joensuu R, Sepponen R. MRI-based D2O/H2O-contrast method to study water flow and distribution in heterogeneous systems: Demonstration in wood xylem. Journal of Magnetic Resonance. 2001;149:36-44
9. Ali TS, Tourell MC, Hugo HJ, Pyke C, Yang S, Lloyd T. et al. Transverse relaxation-based assessment of mammographic density and breast tissue composition by single-sided portable NMR. Magn Reson Med. 2019;82:1199-213
10. Wang FN, Peng SL, Lu CT, Peng HH, Yeh TC. Water signal attenuation by D2O infusion as a novel contrast mechanism for 1H perfusion MRI. NMR Biomed. 2013;26:692-8
11. Walczak P, Wojtkiewicz J, Nowakowski A, Habich A, Holak P, Xu J. et al. Real-time MRI for precise and predictable intra-arterial stem cell delivery to the central nervous system. J Cereb Blood Flow Metab. 2017;37:2346-58
12. Zhang J, Li Y, Slania S, Yadav NN, Liu J, Wang R. et al. Phenols as Diamagnetic T2 -Exchange Magnetic Resonance Imaging Contrast Agents. Chemistry. 2018;24:1259-63
13. Chu C, Jablonska A, Lesniak WG, Thomas AM, Lan X, Linville RM. et al. Optimization of osmotic blood-brain barrier opening to enable intravital microscopy studies on drug delivery in mouse cortex. J Control Release. 2020;317:312-21
14. Chu C, Liu G, Janowski M, Bulte JWM, Li S, Pearl M. et al. Real-Time MRI Guidance for Reproducible Hyperosmolar Opening of the Blood-Brain Barrier in Mice. Front Neurol. 2018;9:921
15. Kim SG, Ackerman JJ. Multicompartment analysis of blood flow and tissue perfusion employing D2O as a freely diffusible tracer: a novel deuterium NMR technique demonstrated via application with murine RIF-1 tumors. Magn Reson Med. 1988;8:410-26
16. Lin AL, Qin Q, Zhao X, Duong TQ. Blood longitudinal (T1) and transverse (T2) relaxation time constants at 11.7 Tesla. MAGMA. 2012;25:245-9
17. Katz JJ, Crespi HL, Hasterlik RJ, Thomson JF, Finkel AJ. Some observations on biological effects of deuterium, with special reference to effects on neoplastic processes. J Natl Cancer Inst. 1957;18:641-59
18. Lichtenbelt WDV, Westerterp KR, Wouters L. Deuterium Dilution as a Method for Determining Total-Body Water - Effect of Test Protocol and Sampling Time. Brit J Nutr. 1994;72:491-7
19. Pullicino E, Coward WA, Stubbs RJ, Elia M. Bedside and field methods for assessing body composition: comparison with the deuterium dilution technique. Eur J Clin Nutr. 1990;44:753-62
20. Wallace SA, Mathur JN, Allen BJ. The influence of heavy water on boron requirements for neutron capture therapy. Med Phys. 1995;22:585-90
21. Obata T, Ikehira H, Shishido F, Fukuda N, Ueshima Y, Koga M. et al. Deuterium MR in vivo imaging of the rat eye using 2H2O. Acta Radiologica. 1995;36:552-5
22. Detre JA, Subramanian V, Mitchell M, Smith D, Kobayashi A, Zaman A. et al. Measurement of regional cerebral blood flow in cat brain using intracarotid 2H2O and 2H NMR imaging. Magnetic resonance in medicine. 1990;14:389-95
23. Kim SG, Ackerman JJ. Quantification of regional blood flow by monitoring of exogenous tracer via nuclear magnetic resonance spectroscopy. Magn Reson Med. 1990;14:266-82
24. Ackerman JJ, Ewy CS, Becker NN, Shalwitz RA. Deuterium nuclear magnetic resonance measurements of blood flow and tissue perfusion employing 2H2O as a freely diffusible tracer. Proc Natl Acad Sci U S A. 1987;84:4099-102
25. Eskey CJ, Koretsky AP, Domach MM, Jain RK. 2H-nuclear magnetic resonance imaging of tumor blood flow: spatial and temporal heterogeneity in a tissue-isolated mammary adenocarcinoma. Cancer Res. 1992;52:6010-9
26. Assmann JC, Brender JR, Farthing DE, Saito K, Warrick KA, Maglakelidze N. et al. In vivo deuterated water labeling followed by deuterium magnetic resonance imaging of xenografted tumors. bioRxiv. 2019: 809327.
27. Zhong JH, Gore JC, Armitage IM. Quantitative studies of hydrodynamic effects and cross-relaxation in protein solutions and tissues with proton and deuteron longitudinal relaxation times. Magn Reson Med. 1990;13:192-203
28. Tadimalla S, Momot KI. Effect of partial H2O-D2O replacement on the anisotropy of transverse proton spin relaxation in bovine articular cartilage. PLoS One. 2014;9:e115288
29. Zhong JH, Gore JC, Armitage IM. Relative contributions of chemical exchange and other relaxation mechanisms in protein solutions and tissues. Magn Reson Med. 1989;11:295-308
30. Sehy JV, Banks AA, Ackerman JJ, Neil JJ. Importance of intracellular water apparent diffusion to the measurement of membrane permeability. Biophys J. 2002;83:2856-63
31. Zawadzki M, Walecki J, Kostkiewicz B, Kostyra K, Pearl MS, Solaiyappan M. et al. Real-time MRI guidance for intra-arterial drug delivery in a patient with a brain tumor: technical note. BMJ Case Rep. 2019 12
32. Xu X, Xu J, Chan KWY, Liu J, Liu H, Li Y. et al. GlucoCEST imaging with on-resonance variable delay multiple pulse (onVDMP) MRI. Magn Reson Med. 2019;81:47-56
33. Bulte JW, Douglas T, Witwer B, Zhang SC, Strable E, Lewis BK. et al. Magnetodendrimers allow endosomal magnetic labeling and in vivo tracking of stem cells. Nat Biotechnol. 2001;19:1141-7
34. Lorthois S, Duru P, Billanou I, Quintard M, Celsis P. Kinetic modeling in the context of cerebral blood flow quantification by H2(15)O positron emission tomography: the meaning of the permeability coefficient in Renkin-Crones model revisited at capillary scale. J Theor Biol. 2014;353:157-69
Corresponding author: Guanshu Liu, Ph.D., 707 N. Broadway, Baltimore, MD 21205; Phone (office): 443-923-9500, Fax: 410-614-3147; Email: guanshujhu.edu
Received 2020-11-15
Accepted 2021-3-20
Published 2021-4-15