Purpose:

An accurate and noninvasive assessment of tumor response following treatment other than traditional anatomical imaging techniques is essential. Deuterium magnetic resonance spectroscopic (MRS) imaging has been demonstrated as an alternative for cancer metabolic imaging by high-field MRI using deuterium-labeled molecules. The study aim was to use 2H tissue labeling and deuterium MRI at clinical field strength for tumor visualization and assessment of three anticancer therapies in pancreatic cancer model mice.

Experimental Design:

MIA PaCa-2 pancreatic carcinoma and C26 colorectal carcinoma models of BALB/c-nu mice was prepared, and repeated deuterium MRI was performed during the first 10 days of free drinking of 30% D2O to track 2H distribution in tissues. 2H accumulation in the tumor after irradiation, bevacizumab administration, or gemcitabine administration was also measured in MIA PaCa-2–bearing mice. Confirmatory proton MRI, ex vivo metabolic hyperpolarization 13C-MRS, and histopathology were performed.

Results:

The mouse's whole-body distribution of 2H was visible 1 day after drinking, and the signal intensity increased daily. Although the tumor size did not change 1 and 3 days after irradiation, the amount of 2H decreased significantly. The 2H image intensity of the tumor also significantly decreased after the administration of bevacizumab or gemcitabine. Metabolic hyperpolarization 13C-MRS, proton MRI, and 2H-NMR spectroscopy confirmed the efficacy of the anticancer treatments.

Conclusions:

Deuterium MRI at 1.5T proved feasible to track 2H distribution throughout mouse tissues during D2O administration and revealed a higher 2H accumulation in the tumor xenografts. This research demonstrated a promising successful method for preliminary assessment of radiotherapy and chemotherapy of cancer.

Translational Relevance

We introduce an imaging approach that applies 2H-labeling followed by deuterium MRI for the visualization of body tissues including tumor xenografts in mice. These xenografts showed a day-by-day cumulative 2H enrichment. We could also investigate the efficacy of three anticancer treatment strategies. The initial response of the tumor to radiation therapy and chemotherapeutic agents, in the treatment of human pancreatic cancer model, MIA PaCa-2 was assessed. We observed significantly decreased accumulation of 2H in xenografts of treated mice at preliminary stages of the treatments. This research demonstrated a promising clinical use of deuterium MRI for providing unique metabolic information essential for assessment of anticancer treatment protocols and monitoring the treatment of patients with pancreatic cancer. The approach could be used as a supporting or potentially an alternative method to detect treatment responses of patients with cancer to drugs in early-stage clinical trials and to guide clinical treatment of individual patients.

Accurate and noninvasive assessment of tumor response following radiation and/or chemotherapy is essential for managing therapeutic regimens in patients. Traditionally, anatomic imaging techniques, such as CT and MRI have been used to plan radiotherapy and assess tumor response. In addition, morphologic evaluation by MRI is now an indispensable technique for accurately diagnosing diseases, such as cancer and tracking therapeutic effects. However, anatomic imaging techniques are insufficient for the early detection of therapeutic responses after therapy because a tumor's size may not change for weeks or months after treatment.

Molecular imaging visualizes and monitors metabolic processes in the living body at the molecular level. Diagnostic methods, such as PET and MRI that target molecular processes, have been clinically applied in a wide range of conditions, and expectations for use in personalized medicine and effective drug discovery have been increasing. Although PET can track the biochemical changes in tumors after therapy, it is frequently postponed because of the confounding effects of inflammation and occasionally low predictive values. Moreover, PET-CT scans have relatively poor spatial resolution and their frequent use may expose patients to excessive ionizing radiation (1).

MRI can also be used to visualize bodily structures, but it has limited molecular-level imaging capabilities because of its low detection sensitivity for trace amounts of substances. The sensitivity of MRI depends on the polarization rate of the nuclear spin. At normal temperatures, the directions of nuclear spins are complicated, and the signals cancel each other out. Therefore, it is necessary to increase nuclear spin polarization to increase MRI sensitivity. Since its development more than a decade ago (2), hyperpolarization via dynamic nuclear polarization (DNP) has rapidly become a well-established technique to enhance the NMR signals of low-gamma nuclei that have a long liquid-state longitudinal relaxation time (T1), such as 13C, 15N, 6Li, and 129Xe (3–6). Hyperpolarization MRI enables metabolic imaging of administered probes, so it is used to visualize the microenvironment of the living body by observing differences in the metabolism of molecules in diseased areas, such as those in cancer. For example, when cancer is tested by DNP–magnetic resonance spectroscopy (MRS), pyruvic acid, a major metabolic intermediate of glucose, is labeled with 13C, hyperpolarized, and administered intravenously as a metabolic probe. By measuring real-time metabolism by MRI, the difference in metabolism between normal cells and cancerous cells can be evaluated (6). Currently, DNP is in practical use and clinical research, but it is laborious to disseminate the DNP device to general hospitals because the introductory cost is high, and its technical requirements restrict its applications to a limited number of suitable metabolites.

In recent years, deuterium metabolic imaging (DMI) has been demonstrated, and recent MR studies have concluded that deuterated metabolites possess desirable properties for the noninvasive mapping of biochemical processes of healthy and diseased organs (7). A promising metabolic imaging approach is to track the fate of orally administered deuterated glucose because it is taken up and metabolized into different products in tissues, such as the heart, brain, pancreas, liver, and cancerous tissue (8–12).

The stable hydrogen isotope, 2H, is used as a tracer in metabolic imaging for studying the dynamic metabolism of compounds. D2O was used in research as a 2H uptake label by administration of D2O to mice and as a free diffusivity tracer for blood flow and tissue circulation (13). Deuterium MR measurements are of growing interest in metabolic imaging because 2H has an integer spin of one, which causes a quadrupolar magnetic moment, and it is sensitive to molecular orientation and produces faster relaxation than 1H (14, 15). When compared with 1H, 2H has a lower natural abundance (∼0.015%) and a gyromagnetic ratio (6.54 MHz/T). These characteristics cause a reduced MR signal. However, the short longitudinal relaxation time (T1) enables rapid signal averaging with an SNR per unit acquisition time gain that is approximately 3×–10× higher for 2H than for 1H. This has led to the use of 2H-labeled substrates as in vivo metabolic tracers and to an increased interest in deuterium MR measurements (9, 16). Therefore, DMI has become an emerging MR modality for investigating energy metabolism in vivo (9, 12, 17). MR measurements that depend on tracing 2H are featured by its insensitivity to magnetic field inhomogeneity, increased sensitivity through signal averaging, and technical simplicity. Moreover, 2H is a stable nonradioactive isotope that does not require special handling facilities. Accordingly, 2H can be traced at clinical magnetic field strengths without the need for special techniques or facilities. In addition, D2O administration is known to lead to the incorporation of 2H into metabolic pathways and tissue labeling. Tracing the 2H-labeled organs by MRI should also help characterize D2O kinetics in normal and diseased body tissues.

The study aim was to investigate D2O kinetics in bodily tissues and track 2H-labeled tissues during oral D2O administration using a magnetic field strength of 1.5T, which is widely used in clinical practice. We also hoped to demonstrate the feasibility of tumor imaging and delineation of the efficacy of anticancer radiation and chemotherapeutics in a human pancreatic cancer model based on a 1.5T deuterium MRI.

Chemicals

Heavy water (D2O, 98.8%) and deuterated methanol (Methanol-D4, 99.95%) were obtained from Sigma-Aldrich Corp. Pyruvic acid (1–13C, 99%) was obtained from Cambridge Isotope Laboratories, Inc. Trityl radical OX063 was supplied by GE Healthcare. Bevacizumab (trade name, Avastin) was bought from Chugai Pharmaceutical Co., Ltd. Gemcitabine hydrochloride was bought from Tokyo Chemical Industry Co., Ltd. All other chemicals were of reagent-grade quality and obtained from commercial sources.

Animals and dietary treatments

All animal care and experimental protocols were approved by the Committee for Animal Research and Welfare of Medical School, Gifu University and conducted according to the recommendations of the Committee for Care and Use of Laboratory Animals, Gifu University. Female 4- to 6-week-old BALB/c-nu mice (RRID:IMSR_CRL:194) weighing 15–20 g were purchased from Charles River Laboratories Japan, Inc. Before the experiments, all animals were acclimated to a normal diet for 1 week with free access to distilled water and appropriate food (MF diet, Oriental Yeast Co.). Every effort was made to minimize the number of mice used and their suffering.

Cell culture and tumor implantation

MIA PaCa-2 human pancreatic carcinoma cells (catalog no. CRL-1420TM, RRID: CVCL_0428) were purchased from ATCC and C26 murine colorectal carcinoma cells (Cell# RCB2657, RRID:CVCL_0240) were obtained from RIKEN Cell Bank. Cells were grown in DMEM (Wako Pure Chemical Industries, Ltd.), supplemented with 10% heat-inactivated FBS (MP Bio LLC.) and penicillin–streptomycin (100 U/mL penicillin and 100 mg/mL streptomycin, Gibco Life Technologies) at 37°C in a humidified incubator containing 5% CO2. A pathologic tumor model was developed by subcutaneous administration of 0.8–1.2 × 106 mouse cells into the right leg of each mouse. The tumor diameters were measured with a Vernier caliper. When the tumor volume reached approximately 500 mm3, the experiments were started. The following approximate formula was used to determine the tumor volume: tumor volume = length/2 × width × height. Before the experiments, the body weight was between 18g and 20g.

Phantom study

To evaluate the feasibility of visualizing 2H at 1.5T MRI, a phantom of nine wells was designed and prepared by a 3D printer for use. Different concentrations of D2O in H2O [10%–90% (v/v)] were added to each well. 1H-specific and 2H-specific coils at resonance frequencies of 64 and 9.8 MHz, respectively, were used to scan the phantom. The phantom was imaged by proton and deuterium MRI methods at different repetition times (1H: 500–3,000 ms, 2H: 100–500 ms).

In vivo2H-labeling

For the in vivo imaging studies, mice were provided free access to drinking water containing 30% v/v D2O during imaging and treatment periods. The control mice were allowed free access to regular tap water (0.015% D2O) instead of deuterated water.

Irradiation and drug administration

The effects of X-ray irradiation, bevacizumab, and gemcitabine on tumor progression and the 2H MR signals in tumor tissues were determined in mice randomly divided into four groups (n = 4): the control group (intraperitoneal injection of 100 μL physiologic saline, at 1-day intervals), bevacizumab group (intraperitoneal injection of bevacizumab 10 mg/kg, at 1-day intervals), gemcitabine group (intraperitoneal injection of gemcitabine 120 mg/kg at 1-day interval), and the irradiation group (single dose of 20 Gy X-ray irradiation). For irradiation, the mice were restrained without anesthesia in a custom jig. The xenograft-bearing legs were then irradiated in an X-ray medical linear accelerator (LINAC), Primus (Siemens Healthcare). During all the treatment experiments, mice were allowed free access to deuterated water simultaneously with treatment. Mice were sacrificed by cervical dislocation and xenografts were removed, fixed in 10% formalin for histopathologic examination, or kept at −80°C for the NMR and hyperpolarized 13C-MRS studies.

Proton and deuterium MRI

Xenograft-bearing leg images were obtained on a 1.5T permanent magnet MRI system (Medalist) obtained from Japan Redox Ltd. Mice were anesthetized with isoflurane (3% for induction and 1.0%–2.0% for maintenance in medical air, 350 mL/minute). Using a PicoM nonmagnetic rectal fiber optic temperature sensor (Opsens Inc.), the mouse body temperature was monitored and kept at 36°C ± 1°C using a warm air flow set at 38°C. The mouse body was then fixed in a prone position on a customized holder using adhesive skin tape. The holder was set into a 1H or 2H volume coil of the MR scanner.

Deuterium MR scanning was performed using gradient echo coronal or transverse 2D deuterium imaging. After a quick assessment of the sample position by a fast low-angle shot pilot sequence, 2H images were obtained for the head and forelimbs, abdomen, and hind limbs. During the treatment, images of xenograft-bearing legs were also acquired. One 50-mm–thick coronal or 15-mm–thick transverse slice was acquired at a 100 ms repetition time (TR), echo time (TE) of 10 ms, 64 accumulations, and a 64 × 64 matrix. During each scan, a marker phantom, which was a 5-mm diameter vial containing 5% D2O in H2O (v/v), was placed next to the target ROI as a 2H reference signal. Coronal and transverse T1- and T2-weighted proton images were also acquired by a spin echo sequence with TR of 500 ms or 3,000 ms, TE of 10 ms or 40 ms, two or four accumulations, five 2-mm slices, and a 64 × 64 matrix.

Hyperpolarized 13C-MRS

Hyperpolarization was achieved by using 22 μL of 11 M 13C-pyruvate with 15 mmol/L OX063. A HyperSense DNP polarizer (Oxford Instruments) was set at a microwave irradiation frequency of 93.42 GHz at 100 mW and 1.4K for 45 minutes. The hyperpolarized probe solution was mixed with 3 mL dissolution buffer containing 50 mmol/L trisaminomethane, 0.3 mmol/L EDTA, and 60 mmol/L NaOH. Immediately after dissolution, 250 μL of hyperpolarized probe solution was mixed in a prewarmed NMR tube containing 350 μL of 1:1 tumor homogenate in PBS and 50 μL of 500 mmol/L nicotinamide adenine dinucleotide. 13C spectral acquisition was then done for 260 seconds using a 1.4T Spinsolve 60 Carbon High-Performance Benchtop NMR instrument (Magritek) at a flip angle = 10° and TR = 4 seconds.

2H-NMR spectroscopy

Tumor tissues were collected and homogenized in 100 mmol/L phosphate buffer at a 1:1 rate and centrifuged at 5,000 rpm for 20 minutes at 4°C. Supernatants were collected, and methanol-d4 was added at a final concentration of 100 mmol/L in a total of 600 μL. Samples containing the reference standard were placed into a 5 mm NMR tube and examined for 2H content by 2H-NMR. For that purpose, an Avance 800 NMR spectrometer (Bruker BioSpin) was used. The spectrometer was operated at a 1H frequency of 800.15 MHz and a 5-mm 1H inverse detection probe with triple-axis gradient coils was used for all measurements. The 2H spectra were successively acquired at 25°C, with 4,096 complex points covering 1,476 Hz for 2H. NMR data processing and analysis were performed using the TopSpin NMR software package v3.6.2 (Bruker BioSpin).

Histopathology and IHC

After collection, xenografts were fixed in 10% formalin and embedded in paraffin. Tissues were cut into 3 μm–thick sections, deparaffinized, and stained with hematoxylin and eosin (H&E) or used for IHC analysis. For antigen retrieval, the sections were placed in citrate buffer (pH 6.0) and autoclaved at 110°C for 10 minutes. The sections were then washed three times in PBS for 5 minutes each and then washed and blocked with 3% hydrogen peroxide in methanol for 10 minutes to remove endogenous peroxidase. Next, the sections were rinsed three times in PBS for 5 minutes again. Nonspecific binding sites were blocked by adding PBS (pH 7.4) containing 2% BSA (Fujifilm Holdings Corporation) for 40 minutes at room temperature and then incubated overnight at 4°C with primary antibodies against the cellular proliferation marker, anti-human Ki-67 (dilution, 1:100; M7240; Dako, RRID:AB_2142367), tumor angiogenesis marker, and anti-mouse CD31 (dilution, 1:50; DIA 310; Dinova Gmbh, RRID:AB_2631039). Primary antibodies were detected by using Simple Stain Mouse MAX PO (rat; Nichirei Bioscience) for 60 minutes at room temperature. The bound complex was visualized using SIGMAFAST 3,30-diaminobenzidine tablets (Sigma-Aldrich Corporation) and counterstained with hematoxylin. The sections were mounted and images were obtained by using a Keyence all-in-one fluorescence BZ-X800 microscope (KEYENCE) and analyzed with imaging cellSens software (RRID:SCR_014551, Olympus).

Statistical analysis

All experiments were performed at least three times. All results shown in the graphs are expressed as the mean ± SD. Differences in the group means were determined by a two-tailed Student t test. The minimum level of significance was set at P < 0.05.

Data availability

Representative data of MRI, histopathology, and hyperpolarized 13C-MRS experiments are presented in this article and the Supplementary Data. Raw NMR data can be accessed at https://zenodo.org/record/8333098. All other materials, data, and protocols generated in this study are available upon reasonable request from the corresponding author.

The D2O phantom is visually measurable at 1.5T by deuterium MRI

We first used a phantom containing various concentrations of D2O in H2O (v/v) to evaluate the sensitivity of 2H detection by a 1.5T MR scanner and 2H-specific coil. The D2O was clearly visible at concentrations as low as 10% (Fig. 1A and B). The signal intensity at 9.8 MHz was directly related to the D2O content. Acceptable image quality of both H2O (by proton MRI) and D2O by deuterium MRI were obtained at a repetition time of 500 ms in both cases (Fig. 1B).

Figure 1.

The D2O phantom is visually measurable at 1.5T by deuterium MRI. A, A 3D phantom with nine wells designed to evaluate the feasibility of visualization of 2H at 1.5T MRI. Each well contained a different concentration of D2O. B, Representative coronal MRI images of the deuterium and proton coils (top) of the D2O phantom and graphical representation of the signal intensities obtained (bottom).

Figure 1.

The D2O phantom is visually measurable at 1.5T by deuterium MRI. A, A 3D phantom with nine wells designed to evaluate the feasibility of visualization of 2H at 1.5T MRI. Each well contained a different concentration of D2O. B, Representative coronal MRI images of the deuterium and proton coils (top) of the D2O phantom and graphical representation of the signal intensities obtained (bottom).

Close modal

In vivo2H-labeling by D2O followed by deuterium MRI enables visualization of 2H MR signal in tissues

To evaluate in vivo deuterium MRI at 1.5T, mice with MIA PaCa-2 or C26 xenografts were freely allowed to drink 30% D2O and imaged for D2O kinetics. During the 10-day administration, the distribution and accumulation of 2H in different mouse tissues were studied. Imaging of each of the three body sections was performed. The distribution of 2H in the tissues was observable as early as 1 day after D2O administration (Fig. 2A; Supplementary Fig. S1A). Moreover, upon continuous drinking of D2O, the 2H atoms of D2O gradually replaced the protons in the tissues to various degrees, and the signal strength was gradually increasing on the 3rd, 7th, and 10th days (Fig. 2A; Supplementary Fig. S1A). Proton images of the whole mouse body were also acquired, and each image was superimposed to obtain a clear image of the whole-body 2H distribution (Fig. 2B and C; Supplementary Fig. S1B and S1C). Interestingly, the highest 2H tissue accumulation was observed in the tumor tissues (Fig. 2D; Supplementary Fig. S1D).

Figure 2.

Sensitivity of in vivo deuterium MRI and kinetics of 2H enrichments in BALB/c-nu mice allowed to drink D2O freely. A, Representative coronal deuterium MRI images of the distribution of 2H in the body during D2O administration in mice. Images were collected on days 1, 3, 7, and 10 of the first exposure to D2O. B, Representative T1-weighted MRI images of mouse body parts. C, Full-body MRI images after placing the MRI images of three body parts together show the distribution of 2H and 1H signals collected throughout the body. D, Quantification of 2H MRI signals of body tissues. Analysis was performed using ImageJ software. The data show the means ± SD (n = 6).

Figure 2.

Sensitivity of in vivo deuterium MRI and kinetics of 2H enrichments in BALB/c-nu mice allowed to drink D2O freely. A, Representative coronal deuterium MRI images of the distribution of 2H in the body during D2O administration in mice. Images were collected on days 1, 3, 7, and 10 of the first exposure to D2O. B, Representative T1-weighted MRI images of mouse body parts. C, Full-body MRI images after placing the MRI images of three body parts together show the distribution of 2H and 1H signals collected throughout the body. D, Quantification of 2H MRI signals of body tissues. Analysis was performed using ImageJ software. The data show the means ± SD (n = 6).

Close modal

Deuterium MRI enabled tracking of the efficacy of the radiotherapy

We studied the feasibility of applying deuterium MRI for monitoring the efficacy of the different anticancer interferences. The MIA PaCa-2 xenograft-bearing legs of mice were subjected to a single dose of radiation, and the mice were freely allowed to drink 30% D2O. Deuterium MRI of the tumor was then performed on days 1, 3, and 10. The changes in the tumor size were also measured. The accumulation of 2H in the tumor tissue was significantly reduced in the treated mice at days 1, 3, and 10 posttreatments (Fig. 3A and B). On the other hand, there were no significant differences in the tumor volumes during this early stage of tumor treatment. The effect of treatment on tumor volumes was later confirmed by continuous monitoring of tumor volume and by proton imaging up to day 27. The MRI data also showed a significant decrease in the tumor volume in the treated mice xenografts (Fig. 3C and D).

Figure 3.

Monitoring of the radiotherapy efficacy by deuterium MRI. A, Representative transverse deuterium MRI images of the MIA PaCa-2 model BALB/c-nu mice legs after irradiation and during D2O administration. Images were collected on days 1, 3, and 10 of irradiation. The images shown next to the deuterium images represent transverse proton MRI images of the tumors on day 10 after treatment started. B, Graphical representation of the quantified 2H MRI signals collected from tumor xenografts. Analysis was performed using ImageJ software. The data show the means ± SD. *, P < 0.05; **, P < 0.01 (n = 4, each group). C, Graphical representation of the average tumor volumes as a function of time after irradiation of mice. The data show the means ± SD (n = 4, each group). D, Representative transverse T1- and T2-weighted MR images of mouse hindlimbs with tumor xenografts on day 27 of irradiation obtained by 1.5T animal MRI.

Figure 3.

Monitoring of the radiotherapy efficacy by deuterium MRI. A, Representative transverse deuterium MRI images of the MIA PaCa-2 model BALB/c-nu mice legs after irradiation and during D2O administration. Images were collected on days 1, 3, and 10 of irradiation. The images shown next to the deuterium images represent transverse proton MRI images of the tumors on day 10 after treatment started. B, Graphical representation of the quantified 2H MRI signals collected from tumor xenografts. Analysis was performed using ImageJ software. The data show the means ± SD. *, P < 0.05; **, P < 0.01 (n = 4, each group). C, Graphical representation of the average tumor volumes as a function of time after irradiation of mice. The data show the means ± SD (n = 4, each group). D, Representative transverse T1- and T2-weighted MR images of mouse hindlimbs with tumor xenografts on day 27 of irradiation obtained by 1.5T animal MRI.

Close modal

Deuterium MRI enabled tracking of the efficacy of the chemotherapeutics

On a similar basis, the effect of anticancer drugs on 2H MR signal was assessed in the xenograft-bearing mice after treatment with bevacizumab or gemcitabine. The mice were intraperitoneally injected while allowed to drink D2O for 1 week, and deuterium MRI was performed on days 1, 3, and 7 of treatment. Similar to the pattern for radiotherapy, despite no obvious change in tumor volume, the 2H MR signals in the tumor xenografts were significantly lower during the treatment period in the treated mice than in the control mice (Fig. 4AD). As measured by total ingested volumes per day, there was no significant difference in the amount of D2O consumed by the treated mice groups.

Figure 4.

Monitoring of the efficacy of chemotherapeutics by deuterium MRI. A, Representative transverse deuterium MRI images of the MIA PaCa-2 model BALB/c-nu mice legs during gemcitabine treatment and D2O administration. Images were collected on days 1, 3, and 7 of administration. B, Graphical representation of the quantified 2H-MRI signals collected from tumor xenografts. Analysis was performed using ImageJ software. The data show the means ± SD. *, P < 0.05; **, P < 0.01 (n = 4, each group). C, Representative transverse deuterium MRI images of the MIA PaCa-2 model mice legs during bevacizumab treatment and D2O administration. Images were collected on days 1, 3, and 7 of administration. D, Graphical representation of the quantified 2H-MRI signals collected from tumor xenografts. Analysis was performed using ImageJ software. The data show the means ± SD. *, P < 0.05 (n = 4, each group).

Figure 4.

Monitoring of the efficacy of chemotherapeutics by deuterium MRI. A, Representative transverse deuterium MRI images of the MIA PaCa-2 model BALB/c-nu mice legs during gemcitabine treatment and D2O administration. Images were collected on days 1, 3, and 7 of administration. B, Graphical representation of the quantified 2H-MRI signals collected from tumor xenografts. Analysis was performed using ImageJ software. The data show the means ± SD. *, P < 0.05; **, P < 0.01 (n = 4, each group). C, Representative transverse deuterium MRI images of the MIA PaCa-2 model mice legs during bevacizumab treatment and D2O administration. Images were collected on days 1, 3, and 7 of administration. D, Graphical representation of the quantified 2H-MRI signals collected from tumor xenografts. Analysis was performed using ImageJ software. The data show the means ± SD. *, P < 0.05 (n = 4, each group).

Close modal

Investigation by 800-MHz NMR confirmed the decreased 2H accumulation during treatment

To quantitatively confirm the observed reductions in the 2H content during treatment, we used the 2H-NMR technique to measure the relative 2H content in the extracts of equal volumes of tumor tissues from the treated and control mice. The NMR spectrometer was equipped with a specific probe for 2H acquisition that provided signals relevant to the concentrations of 2H in the samples. As shown in Supplementary Fig. S2, the 2H content was significantly lower in the gemcitabine-treated mice than in the controls.

Hyperpolarized 13C-MRS confirmed the effects of treatments on cancer metabolism

To observe the effect of anticancer treatments, an ex vivo study of hyperpolarized 13C-MRS was conducted to investigate the changes in the Warburg effect in the control and treated mice. Tumor homogenates were mixed with hyperpolarized 13C-pyruvate, and the 13C-metabolite signals were measured by NMR in real time immediately after mixing. In a pattern matching that of the deuterium MR data, the calculated lactate-to-pyruvate ratios were significantly lower in all treated groups than in the controls (Fig. 5A).

Figure 5.

Investigation of the early effects of anticancer treatments by hyperpolarization 13C-MRS and histopathology in MIA PaCa-2 tumor-bearing BALB/c-nu mice. A, Representative time-resolved 13C-NMR spectra after mixing 250 μL of 80 mmol/L hyperpolarized 13C pyruvate with 300 μL of tumor homogenates. Real-time dynamic changes of 13C-metabolites are graphically shown under each spectrum. Spectra were acquired using a 10° pulse with 4-sec temporal resolution. A quantitative representation of 13C-lactate to 13C-pyruvate ratios calculated by the areas under the curves of MR signal amplitudes is shown on the right of the spectra. The data show the mean ± SD. *, P < 0.05; **, P < 0.01 (n = 4, each group). B, Representative H&E images of the MIA PaCa-2 tumor sections of treated and control mice after 7 days of treatments. A graphical representation of the quantified apoptosis areas is shown on the right. Data show mean ± SD. *, P < 0.05; **, P < 0.01 (n = 3, each group). C, Representative IHC images of the MIA PaCa-2 tumor sections 7 days posttreatment. Slides were incubated with primary antibodies against the cellular proliferation marker, Ki67. A graphical representation of the quantified ratio of positive cells in a microscopic view is shown to the right. The data show the means ± SD (n = 3, each group). D, Representative IHC images of the MIA PaCa-2 tumor sections 7 days posttreatment. Slides were incubated with primary antibodies against the tumor angiogenesis marker, CD31. A graphical representation of the quantified positive microvascular density areas is shown to the right. The data show the means ± SD (n = 3, each group).

Figure 5.

Investigation of the early effects of anticancer treatments by hyperpolarization 13C-MRS and histopathology in MIA PaCa-2 tumor-bearing BALB/c-nu mice. A, Representative time-resolved 13C-NMR spectra after mixing 250 μL of 80 mmol/L hyperpolarized 13C pyruvate with 300 μL of tumor homogenates. Real-time dynamic changes of 13C-metabolites are graphically shown under each spectrum. Spectra were acquired using a 10° pulse with 4-sec temporal resolution. A quantitative representation of 13C-lactate to 13C-pyruvate ratios calculated by the areas under the curves of MR signal amplitudes is shown on the right of the spectra. The data show the mean ± SD. *, P < 0.05; **, P < 0.01 (n = 4, each group). B, Representative H&E images of the MIA PaCa-2 tumor sections of treated and control mice after 7 days of treatments. A graphical representation of the quantified apoptosis areas is shown on the right. Data show mean ± SD. *, P < 0.05; **, P < 0.01 (n = 3, each group). C, Representative IHC images of the MIA PaCa-2 tumor sections 7 days posttreatment. Slides were incubated with primary antibodies against the cellular proliferation marker, Ki67. A graphical representation of the quantified ratio of positive cells in a microscopic view is shown to the right. The data show the means ± SD (n = 3, each group). D, Representative IHC images of the MIA PaCa-2 tumor sections 7 days posttreatment. Slides were incubated with primary antibodies against the tumor angiogenesis marker, CD31. A graphical representation of the quantified positive microvascular density areas is shown to the right. The data show the means ± SD (n = 3, each group).

Close modal

Histopathologic examination revealed variable effects of treatments at the early stage of treatments

When the tumor samples were examined histopathologically, the H&E-stained sections revealed significantly bigger necrotic areas in the bevacizumab- and gemcitabine-treated mice but not the irradiated mice. In addition, the IHC analysis revealed no significant differences in the cell proliferation or microvascular density markers (Fig. 5BD).

The growing demand for metabolic-specific imaging strategies has rekindled interest in deuterium MRI, a reliable procedure that involves administering a 2H-enriched substrate (glucose, acetate, fumarate, etc.) followed by imaging to examine its distribution and metabolism. This technique enables the investigation of metabolic processes under healthy and diseased conditions. Despite its low natural abundance, the short relaxation time of 2H enables rapid radiofrequency pulses without saturation and proper image acquisition (18). Furthermore, deuterium MRI can provide insights into the spatial biochemistry of living tissues in conjunction with routine proton imaging for anatomic localization (19–21). Therefore, deuterium MRI is analogous to other techniques that depend on the administration of metabolites as probes to track metabolic processes. Moreover, deuterium MRI may be superior to these techniques. It allows tracking of the uptake and distribution of administered D2O in target tissues. D2O has been used to study cell kinetics, protein synthesis, and metabolism (22–24). The low natural abundance of 2H (0.015%) is a major advantage for the detection of an exogenous deuterated substrate and its products above the background and serves as a quantitative internal reference. For example, when D-[6,6′-2H2]-glucose was administered, its metabolites of lactate, glutamine, and glutamate were visually measurable through 2H-MRS after approximately 1 hour. That made three-dimensional DMI maps of deuterated glucose metabolites using chemical-shift imaging also feasible (9). However, instead of the costly high doses of deuterated glucose, the administration of D2O water in this study showed that it is still possible to clearly visualize the deuterated molecules’ fate in tissues. Moreover, the variations in uptake, accumulation, and distribution among these tissues are a further advantage of D2O for studying the accumulation of newly formed molecules in highly proliferative tissues, such as tumors (21).

Previous reports have shown that ingestion of water with 2H content >40% affects the homeostasis of mouse bodies and reduces daily intake (25). Therefore, in this experiment, we allowed the mice to freely drink water containing 30% D2O for a maximum of 10 days during which image acquisition was conducted. First, imaging was performed separately of the mouse's head and fore limbs, abdomen, and lower limbs. The distribution and accumulation of 2H in these anatomical areas were visually measurable as early as 1 day after D2O administration. Interestingly, the deuterium MR imaging data revealed that accumulation of 2H was highest in tumor xenografts than in other tissues from the first day of administration (Fig. 2D). That accumulation further increased daily during the 10 days of the study. Similar results were also observed in the C26 xenograft-bearing mice (Supplementary Fig. S1D) which confirm that highly proliferative cancer tissues uptake 2H more than other organs. These findings could be explained by cell proliferation in the presence of 2H atoms inducing 2H-labeling of the newly formed cell molecules by replacing the carbon bonds with protium atoms, especially given that cancerous cells are characterized by rapid cell division and growth. These growing cells have a high demand for carbon and hydrogen to produce new cellular biomass and maintain high rates of proliferation (26, 27). During the synthesis of this cellular biomass, enzymatic processes can use 2H instead of 1H as a substrate, leading to the formation of new carbon–deuterium bonds that are not exchangeable with protons. After D2O administration, 2H is also known to be integrated into amino acids through transamination and other intermediary pathways, and new proteins are then synthesized from these amino acids (24). In addition, 2H is incorporated into the DNA of dividing cells during nucleotide synthesis through the constitutive de novo nucleotide synthesis pathway (21–23). Thus, simple oral administration of D2O may enable in vivo labeling of highly proliferating cells, including cancerous cells (23, 28).

The 2H-labeling of cell components after administration of D2O has long been used through the use of deuterated water or glucose for measuring cell proliferation in different tissues, including slow-turnover cells such as mammary epithelial cells, colon epithelial cells, and vascular smooth muscle cells (29). Administration of D2O also has been used for in vivo follow-up of the proliferation rates of pancreatic islets (22) as a key factor in the pathogenesis of diabetes mellitus. The observed results of higher 2H-labeling rates in tumor tissues are in accordance with previous reports in which D2O enrichment was used as a probe and showed increased DNA synthesis rates after the administration of a carcinogen (30). These data led to our hypothesis that hydrogen will usually accumulate in rapidly proliferating tumor cells during progression because of the high metabolic activity and involvement in the build-up of the newly formed cell biomass and DNA. The data also prompted us to hypothesize that tracing 2H signals in tumors may help monitor tumor progression and assess the irradiation efficacy early.

During treatments, tumor morphology may remain unchanged for weeks or months, rendering anatomic imaging methods inadequate for early detection of therapeutic response. In this study, we evaluated a comprehensive application of deuterium MRI for monitoring tumor development and early treatment efficacy using radiotherapy or chemotherapy, with a special focus on its potential clinical applications in tumor theranostics. Systemic D2O labeling in four groups of mice bearing MIA PaCa-2 tumor xenografts was induced. The mice were then treated with radiation, bevacizumab, gemcitabine, or saline and examined by deuterium MRI. Bevacizumab is a recombinant mAb used as a first-line therapy for treating several types of cancers; it blocks angiogenesis by inhibiting VEGF A, a protein that stimulates angiogenesis in cancer (31–33). Gemcitabine is a pyrimidine antagonist that is widely used in chemotherapy regimens for a number of carcinomas. The drug is used as a first-line treatment for pancreatic cancer (34) by interfering with the DNA of cancer cells, inhibiting the synthesis of DNA necessary for cell division and suppressing cancer division and growth. Despite no morphologic changes, significant reductions in 2H accumulation were observed in MIA PaCa-2 tumor xenografts following ionizing radiation or the intraperitoneal injections of bevacizumab or gemcitabine anticancer drugs. This finding suggested the ability of anticancer therapy to induce cell apoptosis, disruption of blood vessels, and loss of the newly formed carbon–deuterium bonds. This ability was partially confirmed by the higher cellular apoptosis observed in the treated tumors according to histopathologic data (Fig. 5B). In addition, the ex vivo hyperpolarization MRS study of the changes in the metabolism of the hyperpolarized 13C-pyruvate probe in the tumor homogenate of the control and treated mice also confirmed the in vivo results (Fig. 5A). Following injection of prepolarized 13C-pyruvate, changes in the relative lactate observed has been shown to be a marker for tumor progression or early treatment response (35). In addition, Kaggie and colleagues recently compared 13C-HP-MRS using DMI as a noninvasive method for the first time and showed how they simultaneously probed metabolism (36). Their data demonstrated that the two methods provided complementary data. Although it may still be challenging to obtain high-quality metabolite images by deuterium MRI as those of 13C-HP-MRS, this issue might be overcome by improving the tools used. In addition, using D2O as a target probe instead of deuterated glucose may have solved this issue partially. Our images showed acceptable to good signal-to-noise ratios obtained in a few minutes at 1.5T in mice administered D2O orally.

The D2O-labeling/deuterium MRI technique is characterized by its operational simplicity, and lower cost of D2O than other deuterated molecules, such as glucose. Moreover, the fact that D2O enrichment of proliferative tissues is measurable by deuterium MRI provides an additional advantage. The ease of visualization of D2O enrichment of body tissues using deuterium MRI is also a great advantage for outpatient or field studies. This technique may be considered suitable for routine applications in monitoring proliferating tissues because of its simplicity, adaptability, reproducibility, and predictability.

In this study, deuterium MRI at 1.5T was successfully used to track the 2H signal distributed throughout mouse tissues during D2O administration with higher 2H observed in both the MIA PaCa-2 and C26 tumor tissues. The feasibility of using deuterium MRI after D2O administration to detect anticancer treatment responses in a pancreatic cancer model was demonstrated in vivo. The tested radiotherapy or chemotherapeutics decreased the 2H in the tumor tissues before macroscopic changes occur. The method is efficient, safe, economical and promising.

Conclusions

Deuterium MRI at 1.5T was successfully used to track the 2H signal distributed throughout mouse tissues during D2O administration. We found that the oral uptake of D2O in mice at 30% followed by deuterium MRI revealed a higher accumulation of 2H in the tumor tissues than in controls. In addition, the administration of radiation treatment or chemotherapeutic agents decreased the observable 2H in the tumor tissue before the occurrence of macroscopic changes in tumor morphology. Our findings suggest that assessment of tumor progression and treatment based on deuterium MRI after ingestion of D2O, as a potential imaging biomarker, could be useful for assessing responses to radiotherapy and chemotherapy. The method should be an efficient, safe, and noninvasive approach using D2O labeling to monitor tumor progression and treatment in vivo and is simple, economical, and applicable in animals and humans. This research indicated the potential of deuterium MRI as a clinical method for understanding and monitoring anticancer agents, particularly for the treatment of pancreatic cancer. Overall, the results indicate that deuterium MRI, a minimally invasive approach, may offer new ways of studying and characterizing tumor progression and treatment.

No disclosures were reported.

H. Asano: Conceptualization, resources, data curation, formal analysis, visualization, methodology, writing–original draft, writing–review and editing. A.E. Elhelaly: Conceptualization, resources, data curation, software, formal analysis, visualization, methodology, writing–original draft, writing–review and editing. F. Hyodo: Conceptualization, resources, supervision, funding acquisition, visualization. R. Iwasaki: Conceptualization, resources, data curation, software, formal analysis, methodology. Y. Noda: Data curation, software, methodology. H. Kato: Conceptualization, data curation, software. K. Ichihashi: Data curation, investigation, methodology. H. Tomita: Resources, data curation, software, formal analysis, supervision. M. Murata: Resources. T. Mori: Conceptualization, resources, data curation, supervision, visualization. M. Matsuo: Conceptualization, resources, supervision, funding acquisition, methodology.

This research was supported by AMED under Grant Number JP19cm0106435h0002 and JP21fk0108555h0001. This work was also supported by JSPS KAKENHI (grant numbers 22K07768 and 20KK0253), MEXT Quantum Leap Flagship Program (MEXT Q-LEAP) grant number JPMXS0120330644, and JST Fusion Oriented Research for disruptive Science and Technology (grant number JPMJFR2168).

The publication costs of this article were defrayed in part by the payment of publication fees. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.

Note: Supplementary data for this article are available at Clinical Cancer Research Online (http://clincancerres.aacrjournals.org/).

1.
Council
NR
.
Health risks from exposure to low levels of ionizing radiation: BEIR VII phase 2
.
Washington, DC
:
The National Academies Press
;
2006
. p.
422
2.
Ardenkjaer-Larsen
JH
,
Fridlund
B
,
Gram
A
,
Hansson
G
,
Hansson
L
,
Lerche
MH
, et al
.
Increase in signal-to-noise ratio of >10,000 times in liquid-state NMR
.
Proc Natl Acad Sci U S A
2003
;
100
:
10158
63
.
3.
Cudalbu
C
,
Comment
A
,
Kurdzesau
F
,
van Heeswijk
RB
,
Uffmann
K
,
Jannin
S
, et al
.
Feasibility of in vivo 15N MRS detection of hyperpolarized 15N labeled choline in rats
.
Phys Chem Chem Phys
2010
;
12
:
5818
23
.
4.
van Heeswijk
RB
,
Uffmann
K
,
Comment
A
,
Kurdzesau
F
,
Perazzolo
C
,
Cudalbu
C
, et al
.
Hyperpolarized lithium-6 as a sensor of nanomolar contrast agents
.
Magn Reson Med
2009
;
61
:
1489
93
.
5.
Capozzi
A
,
Roussel
C
,
Comment
A
,
Hyacinthe
J-N
.
Optimal glass-forming solvent brings sublimation dynamic nuclear polarization to 129Xe hyperpolarization biomedical imaging standards
.
J Phys Chem C
2015
;
119
:
5020
5
.
6.
Karlsson
M
,
Jensen
PR
,
Duus
,
Meier
S
,
Lerche
MH
.
Development of dissolution DNP-MR substrates for metabolic research
.
Appl Magn Reson
2012
;
43
:
223
36
.
7.
Markovic
S
,
Roussel
T
,
Neeman
M
,
Frydman
L
.
Deuterium magnetic resonance imaging and the discrimination of fetoplacental metabolism in normal and L-NAME-induced preeclamptic mice
.
Metabolites
2021
;
11
:
376
.
8.
Lu
M
,
Zhu
XH
,
Zhang
Y
,
Mateescu
G
,
Chen
W
.
Quantitative assessment of brain glucose metabolic rates using in vivo deuterium magnetic resonance spectroscopy.
J Cereb Blood Flow Metab
2017
;
37
:
3518
30
.
9.
De Feyter
HM
,
Behar
KL
,
Corbin
ZA
,
Fulbright
RK
,
Brown
PB
,
McIntyre
S
, et al
.
Deuterium metabolic imaging (DMI) for MRI-based 3D mapping of metabolism in vivo
.
Sci Adv
2018
;
4
:
eaat7314
.
10.
Kreis
F
,
Wright
AJ
,
Hesse
F
,
Fala
M
,
Hu
D-e
,
Brindle
KM
.
Measuring tumor glycolytic flux in vivo by using fast deuterium MRI
.
Radiology
2020
;
294
:
289
96
.
11.
Li
H
,
Zhu
X-H
,
Zhu
W
,
Lee
B-Y
,
Wiesner
HM
,
Zhang
Y
, et al
.
Dynamic deuterium MRS imaging for studying rat heart energy metabolism in vivo—initial experience
.
ISMRM 2019, Montréal, Canada
.
May 11–16, 2019
.
12.
De Feyter
HM
,
Thomas
MA
,
Behar
KL
,
de Graaf
RA
.
NMR visibility of deuterium-labeled liver glycogen in vivo
.
Magn Reson Med
2021
;
86
:
62
8
.
13.
Larcombe McDouall
JB
,
Evelhoch
JL
.
Deuterium nuclear magnetic resonance imaging of tracer distribution in D2O clearance measurements of tumor blood flow in mice
.
Cancer Res
1990
;
50
:
363
9
.
14.
Gursan
A
,
Froeling
M
,
Hendriks
AD
,
Welting
D
,
Kentgens
APM
,
Klomp
DWJ
, et al
.
Residual quadrupolar couplings observed in 7 Tesla deuterium MR spectra of skeletal muscle
.
Magn Reson Med
2022
;
87
:
1165
73
.
15.
Navon
G
,
Shinar
H
,
Eliav
U
,
Seo
Y
.
Multiquantum filters and order in tissues
.
NMR Biomed
2001
;
14
:
112
32
.
16.
De Feyter
HM
,
de Graaf
RA
.
Deuterium metabolic imaging - back to the future
.
J Magn Reson
2021
;
326
:
106932
.
17.
Ruhm
L
,
Avdievich
N
,
Ziegs
T
,
Nagel
AM
,
De Feyter
HM
,
de Graaf
RA
, et al
.
Deuterium metabolic imaging in the human brain at 9.4 Tesla with high spatial and temporal resolution
.
Neuroimage
2021
;
244
:
118639
.
18.
Polvoy
I
,
Qin
H
,
Flavell
RR
,
Gordon
J
,
Viswanath
P
,
Sriram
R
, et al
.
Deuterium metabolic imaging-rediscovery of a spectroscopic tool
.
Metabolites
2021
;
11
:
570
.
19.
de Graaf
RA
,
Hendriks
AD
,
Klomp
DWJ
,
Kumaragamage
C
,
Welting
D
,
Arteaga de Castro
CS
, et al
.
On the magnetic field dependence of deuterium metabolic imaging
.
NMR Biomed
2020
;
33
:
e4235
.
20.
Mahar
R
,
Donabedian
PL
,
Merritt
ME
.
HDO production from [(2)H7]glucose quantitatively identifies warburg metabolism
.
Sci Rep
2020
;
10
:
8885
.
21.
Buxbaum
NP
,
Farthing
DE
,
Maglakelidze
N
,
Lizak
M
,
Merkle
H
,
Carpenter
AC
, et al
.
In vivo kinetics and nonradioactive imaging of rapidly proliferating cells in graft-versus-host disease
.
JCI Insight
2017
;
2
:
e92851
.
22.
Chen
S
,
Turner
S
,
Tsang
E
,
Stark
J
,
Turner
H
,
Mahsut
A
, et al
.
Measurement of pancreatic islet cell proliferation by heavy water labeling
.
Am J Physiol Endocrinol Metab
2007
;
293
:
E1459
64
.
23.
Busch
R
,
Neese
RA
,
Awada
M
,
Hayes
GM
,
Hellerstein
MK
.
Measurement of cell proliferation by heavy water labeling
.
Nat Protoc
2007
;
2
:
3045
57
.
24.
Miller
BF
,
Reid
JJ
,
Price
JC
,
Lin
H-JL
,
Atherton
PJ
,
Smith
K
.
CORP: The use of deuterated water for the measurement of protein synthesis
.
J Appl Physiol
2020
;
128
:
1163
76
.
25.
Katz
JJ
,
Crespi
HL
,
Czajka
DM
,
Finkel
AJ
.
Course of deuteriation and some physiological effects of deuterium in mice
.
Am J Physiol Legacy Content
1962
;
203
:
907
13
.
26.
Lunt
SY
,
Vander Heiden
MG
.
Aerobic glycolysis: meeting the metabolic requirements of cell proliferation
.
Annu Rev Cell Dev Biol
2011
;
27
:
441
64
.
27.
Brender
JR
,
Assmann
JC
,
Farthing
DE
,
Saito
K
,
Kishimoto
S
,
Warrick
KA
, et al
.
In vivo deuterium magnetic resonance imaging of xenografted tumors following systemic administration of deuterated water
.
Scientific Reports
2023
;
13
:
14699
.
28.
Strawford
A
,
Antelo
F
,
Christiansen
M
,
Hellerstein
MK
.
Adipose tissue triglyceride turnover, de novo lipogenesis, and cell proliferation in humans measured with 2H2O
.
Am J Physiol Endocrinol Metab
2004
;
286
:
E577
88
.
29.
Neese
RA
,
Misell
LM
,
Turner
S
,
Chu
A
,
Kim
J
,
Cesar
D
, et al
.
Measurement in vivo of proliferation rates of slow turnover cells by 2H2O labeling of the deoxyribose moiety of DNA
.
Proc Natl Acad Sci USA
2002
;
99
:
15345
50
.
30.
Kim
SJ
,
Turner
S
,
Killion
S
,
Hellerstein
MK
.
In vivo measurement of DNA synthesis rates of colon epithelial cells in carcinogenesis
.
Biochem Biophys Res Commun
2005
;
331
:
203
9
.
31.
Le Rhun
E
,
Preusser
M
,
Roth
P
,
Reardon
DA
,
van den Bent
M
,
Wen
P
, et al
.
Molecular targeted therapy of glioblastoma
.
Cancer Treat Rev
2019
;
80
:
101896
.
32.
Cohen
MH
,
Shen
YL
,
Keegan
P
,
Pazdur
R
.
FDA drug approval summary: bevacizumab (Avastin®) as treatment of recurrent glioblastoma multiforme
.
Oncologist
2009
;
14
:
1131
8
.
33.
Rini
BI
,
Halabi
S
,
Rosenberg
JE
,
Stadler
WM
,
Vaena
DA
,
Archer
L
, et al
.
Phase III trial of bevacizumab plus interferon alfa versus interferon alfa monotherapy in patients with metastatic renal cell carcinoma: final results of CALGB 90206
.
J Clin Oncol
2010
;
28
:
2137
43
.
34.
Zhang
XW
,
Ma
YX
,
Sun
Y
,
Cao
YB
,
Li
Q
,
Xu
CA
.
Gemcitabine in combination with a second cytotoxic agent in the first-line treatment of locally advanced or metastatic pancreatic cancer: a systematic review and meta-analysis
.
Targeted Oncology
2017
;
12
:
309
21
.
35.
Woitek
R
,
Gallagher
FA
.
The use of hyperpolarised 13C-MRI in clinical body imaging to probe cancer metabolism
.
Br J Cancer
2021
;
124
:
1187
98
.
36.
Kaggie
JD
,
Khan
AS
,
Matys
T
,
Schulte
RF
,
Locke
MJ
,
Grimmer
A
, et al
.
Deuterium metabolic imaging and hyperpolarized (13)C-MRI of the normal human brain at clinical field strength reveals differential cerebral metabolism
.
Neuroimage
2022
;
257
:
119284
.
This open access article is distributed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0) license.