skip to main content

Title: Mapping electrical properties heterogeneity of tumor using boundary informed electrical properties tomography (BIEPT) at 7T

To develop and evaluate a boundary informed electrical properties tomography (BIEPT) technique for high‐resolution imaging of tumor electrical properties (EPs) heterogeneity on a rodent tumor xenograft model.


Tumor EP distributions were inferred from a reference area external to the tumor, as well as internal EP spatial variations derived from a plurality of relative transmit B1measurements at 7T. Edge sparsity constraint was enforced to enhance numerical stability. Phantom experiments were performed to determine the imaging accuracy and sensitivity for structures of various EP values, as well as geometrical sizes down to 1.5 mm. Numerical simulation of a realistic rodent model was used to quantify the algorithm performance in the presence of noise. Eleven athymic rats with human breast cancer xenograft were imaged in vivo, and representative pathological samples were acquired for comparison.


Reconstructed EPs of the phantoms correspond well to the ground truth acquired from dielectric probe measurements, with the smallest structure reliably detectable being 3 mm. EPs heterogeneity inside a tumor is successfully retrieved in both simulated and experimental cases. In vivo tumor imaging results demonstrate similar local features and spatial patterns to anatomical MRI and pathological slides. The imaged conductivity of necrotic tissue is higher than that of viable tissues, which agrees with our expectation.


BIEPT enables robust detection of tumor EPs heterogeneity with high accuracy and sensitivity to small structures. The retrieved quantitative EPs reflect tumor pathological features (e.g., necrosis). These results provide strong rationale to further expand BIEPT studies toward pathological conditions where EPs may yield valuable, non‐invasive biomarkers.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Magnetic Resonance in Medicine
Page Range / eLocation ID:
p. 393-409
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Background

    Current preclinical models of metastatic prostate cancer (PCa) require sophisticated technologies and/or genetically engineered cells for the noninvasive monitoring of tumors in remote sites, such as bone. Recent developments in circulating tumor DNA (ctDNA) analysis provide an alternative method for noninvasive tumor monitoring at a low cost. Here, we sought to evaluate human Alu and LINE‐1 ctDNA for the longitudinal measurement of subcutaneous and intratibial human PCa xenograft growth and response to ionizing radiation (IR) through comparison with standard slide caliper and bioluminescence measurements.

    Material and Methods

    Intratibial and subcutaneous xenografts were established in male athymic nude mice using LNCaP cells that stably express firefly luciferase. A subset of tumors was treated with a single dose of IR (CT‐guided focal IR, 6 Gy). Tumor measurements were simultaneously taken by slide caliper (subcutaneous only), in vivo bioluminescence imaging, and quantitative real‐time PCR (qPCR) of human‐specific Alu and LINE‐1 ctDNA for several weeks.


    Levels of ctDNA and bioluminescence increased concordantly with subcutaneous and intratibial tumor growth. A statistically significant correlation (Spearman) was observed between ctDNA and subcutaneous tumor volume (LINE‐1,r = .94 and Alu,r = .95,p < .0001), ctDNA and bioluminescence (LINE‐1,r = .66 and Alu,r = .60,p < .002), and bioluminescence and tumor volume (r = .66,p = .0003). Bioluminescence and ctDNA were also significantly correlated in intratibial tumors (LINE‐1,r = .82 and Alu,r = .81,p < .0001). Following external beam IR, the tumor responses varied briefly by method of measurement, but followed a similar trend. Statistically significant correlations were maintained between ctDNA and slide caliper measurement in irradiated subcutaneous tumors (LINE‐1,r = .64 and Alu,r = .44,p < .02), and ctDNA and bioluminescence in intratibial tumors (LINE‐1,r = .55,p = .018).


    Real‐time qPCR of circulating human Alu and LINE‐1 DNA provides an accurate measurement of subcutaneous and intratibial xenograft burden that is comparable with conventional bioluminescence imaging and slide caliper measurement. Transient differences in measurements were observed following tumor‐targeted IR, but overall all measurements mirrored tumor growth and response.

    more » « less
  2. Abstract Purpose

    One standard method, proton resonance frequency shift, for measuring temperature using magnetic resonance imaging (MRI), in MRI‐guided surgeries, fails completely below the freezing point of water. Because of this, we have developed a new methodology for monitoring temperature with MRI below freezing. The purpose of this paper is to show that a strong temperature dependence of the nuclear relaxation timeT1in soft silicone polymers can lead to temperature‐dependent changes of MRI intensity acquired withT1weighting. We propose the use of silicone filaments inserted in tissue for measuring temperature during MRI‐guided cryoablations.


    The temperature dependence ofT1in bio‐compatible soft silicone polymers was measured using nuclear magnetic resonance spectroscopy and MRI. Phantoms, made of bulk silicone materials and put in an MRI‐compatible thermal container with dry ice, allowed temperature measurements ranging from –60°C to + 20°C.T1‐weighted gradient echo images of the phantoms were acquired at spatially uniform temperatures and with a gradient in temperature to determine the efficacy of using these materials as temperature indicators in MRI. Ex vivo experiments on silicone rods, 4 mm in diameter, inserted in animal tissue were conducted to assess the practical feasibility of the method.


    Measurements of nuclear relaxation times of protons in soft silicone polymers show a monotonic, nearly linear, change with temperature (R2 > 0.98) and have a significant correlation with temperature (Pearson'sr > 0.99,p < 0.01). Similarly, the intensity of the MR images in these materials, taken with a gradient echo sequence, are also temperature dependent. There is again a monotonic change in MRI intensity that correlates well with the measured temperature (Pearson'sr < ‐0.98 andp < 0.01). The MRI experiments show that a temperature change of 3°C can be resolved in a distance of about 2.5 mm. Based on MRI images and external sensor calibrations for a sample with a gradient in temperature, temperature maps with 3°C isotherms are created for a bulk phantom. Experiments demonstrate that these changes in MRI intensity with temperature can also be seen in 4 mm silicone rods embedded in ex vivo animal tissue.


    We have developed a new method for measuring temperature in MRI that potentially could be used during MRI‐guided cryoablation operations, reducing both procedure time and cost, and making these surgeries safer.

    more » « less
  3. Purpose

    The purpose of this study is to develop double diffusion encoding (DDE) MRI methods for clinical use. Microscopic diffusion anisotropy measurements from DDE promise greater specificity to changes in tissue microstructure compared with conventional diffusion tensor imaging, but implementation of DDE sequences on whole‐body MRI scanners is challenging because of the limited gradient strengths and lengthy acquisition times.


    A custom single‐refocused DDE sequence was implemented on a 3T whole‐body scanner. The DDE gradient orientation scheme and sequence parameters were optimized based on a Gaussian diffusion assumption. Using an optimized 5‐min DDE acquisition, microscopic fractional anisotropy (μFA) maps were acquired for the first time in multiple sclerosis patients.


    Based on simulations and in vivo human measurements, six parallel and six orthogonal diffusion gradient pairs were found to be the minimum number of diffusion gradient pairs necessary to produce a rotationally invariant measurement of μFA. Simulations showed that optimal precision and accuracy of μFA measurements were obtained using b‐values between 1500 and 3000 s/mm2. The μFA maps showed improved delineation of multiple sclerosis lesions compared with conventional fractional anisotropy and distinct contrast from T2‐weighted fluid attenuated inversion recovery and T1‐weighted imaging.


    The μFA maps can be measured using DDE in a clinical setting and may provide new opportunities for characterizing multiple sclerosis lesions and other types of tissue degeneration. Magn Reson Med 80:507–520, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

    more » « less
  4. Abstract

    In the field of optical imaging, the ability to image tumors at depth with high selectivity and specificity remains a challenge. Surface enhanced resonance Raman scattering (SERRS) nanoparticles (NPs) can be employed as image contrast agents to specifically target cells in vivo; however, this technique typically requires time-intensive point-by-point acquisition of Raman spectra. Here, we combine the use of “spatially offset Raman spectroscopy” (SORS) with that of SERRS in a technique known as “surface enhanced spatially offset resonance Raman spectroscopy” (SESORRS) to image deep-seated tumors in vivo. Additionally, by accounting for the laser spot size, we report an experimental approach for detecting both the bulk tumor, subsequent delineation of tumor margins at high speed, and the identification of a deeper secondary region of interest with fewer measurements than are typically applied. To enhance light collection efficiency, four modifications were made to a previously described custom-built SORS system. Specifically, the following parameters were increased: (i) the numerical aperture (NA) of the lens, from 0.2 to 0.34; (ii) the working distance of the probe, from 9 mm to 40 mm; (iii) the NA of the fiber, from 0.2 to 0.34; and (iv) the fiber diameter, from 100 µm to 400 µm. To calculate the sampling frequency, which refers to the number of data point spectra obtained for each image, we considered the laser spot size of the elliptical beam (6 × 4 mm). Using SERRS contrast agents, we performed in vivo SESORRS imaging on a GL261-Luc mouse model of glioblastoma at four distinct sampling frequencies: par-sampling frequency (12 data points collected), and over-frequency sampling by factors of 2 (35 data points collected), 5 (176 data points collected), and 10 (651 data points collected). In comparison to the previously reported SORS system, the modified SORS instrument showed a 300% improvement in signal-to-noise ratios (SNR). The results demonstrate the ability to acquire distinct Raman spectra from deep-seated glioblastomas in mice through the skull using a low power density (6.5 mW/mm2) and 30-times shorter integration times than a previous report (0.5 s versus 15 s). The ability to map the whole head of the mouse and determine a specific region of interest using as few as 12 spectra (6 s total acquisition time) is achieved. Subsequent use of a higher sampling frequency demonstrates it is possible to delineate the tumor margins in the region of interest with greater certainty. In addition, SESORRS images indicate the emergence of a secondary tumor region deeper within the brain in agreement with MRI and H&E staining. In comparison to traditional Raman imaging approaches, this approach enables improvements in the detection of deep-seated tumors in vivo through depths of several millimeters due to improvements in SNR, spectral resolution, and depth acquisition. This approach offers an opportunity to navigate larger areas of tissues in shorter time frames than previously reported, identify regions of interest, and then image the same area with greater resolution using a higher sampling frequency. Moreover, using a SESORRS approach, we demonstrate that it is possible to detect secondary, deeper-seated lesions through the intact skull.

    more » « less
  5. Background

    There is growing interest in detecting cerebro‐cerebellar circuits, which requires adequate blood oxygenation level dependent contrast and signal‐to‐noise ratio (SNR) throughout the brain. Although 7T scanners offer increased SNR, coverage of commercial head coils is currently limited to the cerebrum.


    To improve cerebellar functional MRI (fMRI) at 7T with high permittivity material (HPM) pads extending the sensitivity of a commercial coil.

    Study Type

    Simulations were used to determine HPM pad configuration and assess radiofrequency (RF) safety. In vivo experiments were performed to evaluate RF field distributions and SNR and assess improvements of cerebellar fMRI.


    Eight healthy volunteers enrolled in a prospective motor fMRI study with and without HPM.

    Field Strength/Sequence

    Gradient echo (GRE) echo planar imaging for fMRI, turbo FLASH for flip angle mapping, GRE sequence for SNR maps, and T1‐weighted MPRAGE were acquired with and without HPM pads at 7T.


    Field maps, SNR maps, and anatomical images were evaluated for coverage. Simulation results were used to assess SAR levels of the experiment. Activation data from fMRI experiments were compared with and without HPM pads.

    Statistical Tests

    fMRI data were analyzed using FEAT FSL for each subject followed by group level analysis using paired t‐test of acquisitions with and without HPM.


    Simulations showed 52% improvement in transmit efficiency in cerebellum with HPM and SAR levels well below recommended limits. Experiments showed 27% improvement in SNR in cerebellum and improvement in coverage on T1‐weighted images. fMRI showed greater cerebellar activation in individual subjects with the HPM pad present (Z > = 4), especially in inferior slices of cerebellum, with 59% average increase in number of activated voxels in the cerebellum. Group‐level analysis showed improved functional activation (Z > = 2.3) in cerebellar regions with HPM pads without loss of measured activation elsewhere.

    Data Conclusion

    HPM pads can improve cerebellar fMRI at 7T with a commonly‐used head coil without compromising RF safety.

    Level of Evidence: 2

    Technical Efficacy: Stage 1

    J. MAGN. RESON. IMAGING 2018;48:431–440.

    more » « less