Correlation between conductivity and prognostic factors in invasive breast cancer using magnetic resonance electric properties tomography (MREPT)

Bias correction for phase-based cr-MREPT using low resolution B1+ magnitude

ObjectiveFull-form Magnetic Resonance Electrical Properties Tomography (MREPT) requires bothB1+magnitude and phase information. SinceB1+phase can be obtained faster and with higher SNR compared toB1+magnitude, several phase-based methods have been developed for conductivity imaging. However, phase-based methods suffer from a concave bias due to the assumption that∇|B1+|is negligible in the ROI.ApproachIn this paper, we re-derive the central equation of phase-based cr-MREPT without assuming that∇|B1+|is negligible and thus propose a correction method directly integrated into the equation system.Main resultsProposed method successfully corrects the concave bias on both simulated and experimental data and significantly increases image quality.SignificanceThe proposed correction method depends on a very low-resolution|B1+|map, and therefore the imaging time does not increase significantly for obtainingB1+magnitude. Moreover, correction can be achieved using simulatedB1+magnitude, hence completely removing the additional imaging requirement.

Feasibility of undersampled spiral trajectories in MREPT for fast conductivity imaging

PURPOSE

To investigate spiral-based imaging including trajectories with undersampling as a fast and robust alternative for phase-based magnetic resonance electrical properties tomography (MREPT) techniques.

METHODS

Spiral trajectories with various undersampling ratios were prescribed to acquire images from an experimental phantom and a healthy volunteer at 3T. The non-Cartesian acquisitions were reconstructed using SPIRiT, and conductivity maps were derived using phase-based cr-MREPT. The resulting maps were compared between different sampling trajectories. Additionally, a conductivity map was obtained using a Cartesian balanced SSFP acquisition from the volunteer to comparatively demonstrate the robustness of the proposed method.

RESULTS

The phantom and volunteer results illustrate the benefits of the spiral acquisitions. Specifically, undersampled spiral acquisitions display improved robustness against field inhomogeneity artifacts and lowered SD values with shortened readout times. Furthermore, average of conductivity values measured for the cerebrospinal fluid with the spiral acquisitions were 1.703 S/m, indicating a close agreement with the theoretical values of 1.794 S/m.

CONCLUSION

A spiral-based acquisition framework for conductivity imaging with and without undersampling is presented. Overall, spiral-based acquisitions improved robustness against field inhomogeneity artifacts, while achieving whole head coverage with multiple averages in less than a minute.

Editorial for "Problem Solving MRI to Reduce False-Positive Biopsy Related to Breast US: Conductivity vs. DWI vs. Abbreviated Contrast-Enhanced MRI"

Level of Evidence

5

Technical Efficacy Stage

3

Data-driven electrical conductivity brain imaging using 3 T MRI

Magnetic resonance electrical properties tomography (MR-EPT) is a non-invasive measurement technique that derives the electrical properties (EPs, e.g., conductivity or permittivity) of tissues in the radiofrequency range (64 MHz for 1.5 T and 128 MHz for 3 T MR systems). Clinical studies have shown the potential of tissue conductivity as a biomarker. To date, model-based conductivity reconstructions rely on numerical assumptions and approximations, leading to inaccuracies in the reconstructed maps. To address such limitations, we propose an artificial neural network (ANN)-based non-linear conductivity estimator trained on simulated data for conductivity brain imaging. Network training was performed on 201 synthesized T2-weighted spin-echo (SE) data obtained from the finite-difference time-domain (FDTD) electromagnetic (EM) simulation. The dataset was composed of an approximated T2-w SE magnitude and transceive phase information. The proposed method was tested three in-silico and in-vivo on two volunteers and three patients' data. For comparison purposes, various conventional phase-based EPT reconstruction methods were used that ignoreB 1 + magnitude information, such as Savitzky-Golay kernel combined with Gaussian filter (S-G Kernel), phase-based convection-reaction EPT (cr-EPT), magnitude-weighted polynomial-fitting phase-based EPT (Poly-Fit), and integral-based phase-based EPT (Integral-based). From the in-silico experiments, quantitative analysis showed that the proposed method provides more accurate and improved quality (e.g., high structural preservation) conductivity maps compared to conventional reconstruction methods. Representatively, in the healthy brain in-silico phantom experiment, the proposed method yielded mean conductivity values of 1.97 ± 0.20 S/m for CSF, 0.33 ± 0.04 S/m for WM, and 0.52 ± 0.08 S/m for GM, which were closer to the ground-truth conductivity (2.00, 0.30, 0.50 S/m) than the integral-based method (2.56 ± 2.31, 0.39 ± 0.12, 0.68 ± 0.33 S/m). In-vivo ANN-based conductivity reconstructions were also of improved quality compared to conventional reconstructions and demonstrated network generalizability and robustness to in-vivo data and pathologies. The reported in-vivo brain conductivity values were in agreement with literatures. In addition, the proposed method was observed for various SNR levels (SNR levels = 10, 20, 40, and 58) and repeatability conditions (the eight acquisitions with the number of signal averages = 1). The preliminary investigations on brain tumor patient datasets suggest that the network trained on simulated dataset can generalize to unforeseen in-vivo pathologies, thus demonstrating its potential for clinical applications.

Correlation between standardized uptake value of 18F-FDG PET/CT and conductivity with pathologic prognostic factors in breast cancer

We investigated the correlation between standardized uptake value (SUV) of 18F-fluorodeoxyglucose (FDG) positron emission tomography (PET)/computed tomography (CT) and conductivity parameters in breast cancer and explored the feasibility of conductivity as an imaging biomarker. Both SUV and conductivity have the potential to reflect the tumors' heterogeneous characteristics, but their correlations have not been investigated until now. Forty four women diagnosed with breast cancer who underwent breast MRI and 18F-FDG PET/CT at the time of diagnosis were included. Among them, 17 women received neoadjuvant chemotherapy followed by surgery and 27 women underwent upfront surgery. For conductivity parameters, maximum and mean values of the tumor region-of-interests were examined. For SUV parameters, SUVmax, SUVmean, and SUVpeak of the tumor region-of-interests were examined. Correlations between conductivity and SUV were evaluated, and among them, the highest correlation was observed between mean conductivity and SUVpeak (Spearman's correlation coefficient = 0.381). In a subgroup analysis for 27 women with upfront surgery, tumors with lymphovascular invasion (LVI) showed higher mean conductivity than those without LVI (median: 0.49 S/m vs 0.06 S/m, p < 0.001). In conclusion, our study shows a low positive correlation between SUVpeak and mean conductivity in breast cancer. Furthermore, conductivity showed a potential to noninvasively predict LVI status.

Technical note: Multi-receiver combination method for phase-based electrical property tomography of the breast

BACKGROUND

Phase-based electrical property tomography (EPT) is a technique that allows conductivity reconstruction with only phase of the B₁ field under the assumption that the magnitude of the B₁ fields are homogeneous. The more this assumption is violated, the less accurate the reconstructed conductivity. Thus, a method that ensures homogeneity of | B 1 - | $| {{\rm{B}}_1^ - } |$ field is important for breast image using multi-receiver coil.

PURPOSE

To develop a method for multi-receiver combination for phase-based EPT usable for breast EPT imaging in the clinic.

METHODS

Theory of the proposed method is presented. To validate the proposed method, the phantom and in-vivo experiments were conducted. Conductivity images were reconstructed using the transceive phase of the combined image and results were compared with another combination method.

RESULTS

The proposed method's conductivity results were more stable than those of the previous method when | B 1 + | $| {{\rm{B}}_1^ + } |$ was not homogeneous and when the homogeneous contrast region was small. The phantom and in-vivo results indicate that the proposed method produces improved conductivity images than the previous method. The proposed combination method also increased the conductivity contrast between benign and cancerous tissues.

CONCLUSION

The proposed method produced more stable multi-receiver combination for phase-based EPT of the breast in a clinical environment.

Repeatability and Reproducibility Uncertainty in Magnetic Resonance-Based Electric Properties Tomography of a Homogeneous Phantom

Uncertainty assessment is a fundamental step in quantitative magnetic resonance imaging because it makes comparable, in a strict metrological sense, the results of different scans, for example during a longitudinal study. Magnetic resonance-based electric properties tomography (EPT) is a quantitative imaging technique that retrieves, non-invasively, a map of the electric properties inside a human body. Although EPT has been used in some early clinical studies, a rigorous experimental assessment of the associated uncertainty has not yet been performed. This paper aims at evaluating the repeatability and reproducibility uncertainties in phase-based Helmholtz-EPT applied on homogeneous phantom data acquired with a clinical 3 T scanner. The law of propagation of uncertainty is used to evaluate the uncertainty in the estimated conductivity values starting from the uncertainty in the acquired scans, which is quantified through a robust James-Stein shrinkage estimator to deal with the dimensionality of the problem. Repeatable errors are detected in the estimated conductivity maps and are quantified for various values of the tunable parameters of the EPT implementation. The spatial dispersion of the estimated electric conductivity maps is found to be a good approximation of the reproducibility uncertainty, evaluated by changing the position of the phantom after each scan. The results underpin the use of the average conductivity (calculated by weighting the local conductivity values by their uncertainty and taking into account the spatial correlation) as an estimate of the conductivity of the homogeneous phantom.

Application of High-Frequency Conductivity Map Using MRI to Evaluate It in the Brain of Alzheimer's Disease Patients

Background

The previous studies reported increased concentrations of metallic ions, imbalanced Na+ and K+ ions, and the increased mobility of protons by microstructural disruptions in Alzheimer's disease (AD).

Purpose

(1) to apply a high-frequency conductivity (HFC) mapping technique using a clinical 3T MRI system, (2) compare HFC values in the brains of participants with AD, amnestic mild cognitive impairment (MCI), and cognitively normal (CN) elderly people, (3) evaluate the relationship between HFC values and cognitive decline, and (4) explore usefulness of HFC values as an imaging biomarker to evaluate the differentiation of AD from CN.

Materials and Methods

This prospective study included 74 participants (23 AD patients, 27 amnestic MCI patients, and 24 CN elderly people) to explore the clinical application of HFC mapping in the brain from March 2019 to August 2021. We performed statistical analyses to compare HFC maps between the three participant groups, evaluate the association of HFC maps with Mini-Mental State Examination (MMSE) scores, and to evaluate the differentiation between the participant groups for HFC values for some brain areas.

Results

We obtained a good HFC map non-invasively. The HFC value was higher in the AD group than in the CN and MCI groups. MMSE scores were negatively associated with HFC values. Age was positively associated with HFC values. The HFC value in the insula has a high area under the receiver operating characteristic (ROC) curve (AUC) value to differentiate AD patients from the CN participants (Sensitivity [SE] = 82, Specificity [SP] =97, AUC = 0.902, p < 0.0001), better than gray matter volume (GMV) in hippocampus (SE = 79, SP = 83, AUC = 0.880, p < 0.0001). The classification for differentiating AD from CN was highest by adding the hippocampal GMV to the insular HFC value (SE = 87, SP = 87, AUC = 0.928, p < 0.0001).

Conclusion

High-frequency conductivity values were significantly increased in the AD group compared to the CN group and increased with age and disease severity. HFC values of the insula along with the GMV of the hippocampus can be used as an imaging biomarker to improve the differentiation of AD from CN.

Combining deep learning and 3D contrast source inversion in MR-based electrical properties tomography

Magnetic resonance electrical properties tomography (MR-EPT) is a technique used to estimate the conductivity and permittivity of tissues from MR measurements of the transmit magnetic field. Different reconstruction methods are available; however, all these methods present several limitations, which hamper the clinical applicability. Standard Helmholtz-based MR-EPT methods are severely affected by noise. Iterative reconstruction methods such as contrast source inversion electrical properties tomography (CSI-EPT) are typically time-consuming and are dependent on their initialization. Deep learning (DL) based methods require a large amount of training data before sufficient generalization can be achieved. Here, we investigate the benefits achievable using a hybrid approach, that is, using MR-EPT or DL-EPT as initialization guesses for standard 3D CSI-EPT. Using realistic electromagnetic simulations at 3 and 7 T, the accuracy and precision of hybrid CSI reconstructions are compared with those of standard 3D CSI-EPT reconstructions. Our results indicate that a hybrid method consisting of an initial DL-EPT reconstruction followed by a 3D CSI-EPT reconstruction would be beneficial. DL-EPT combined with standard 3D CSI-EPT exploits the power of data-driven DL-based EPT reconstructions, while the subsequent CSI-EPT facilitates a better generalization by providing data consistency.

In vivo electrical conductivity measurement of muscle, cartilage, and peripheral nerve around knee joint using MR-electrical properties tomography

This study aimed to investigate whether in vivo MR-electrical properties tomography (MR-EPT) is feasible in musculoskeletal tissues by evaluating the conductivity of muscle, cartilage, and peripheral nerve around the knee joint, and to explore whether these measurements change after exercise. This prospective study was approved by the institutional review board. On February 2020, ten healthy volunteers provided written informed consent and underwent MRI of the right knee using a three-dimensional balanced steady-state free precession (bSSFP) sequence. To test the effect of loading, the subjects performed 60 squatting exercises after baseline MRI, immediately followed by post-exercise MRI with the same sequences. After reconstruction of conductivity map based on the bSSFP sequence, conductivity of muscles, cartilages, and nerves were measured. Measurements between the baseline and post-exercise MRI were compared using the paired t-test. Test–retest reliability for baseline conductivity was evaluated using the intraclass correlation coefficient. The baseline and post-exercise conductivity values (mean ± standard deviation) [S/m] of muscles, cartilages, and nerves were 1.73 ± 0.40 and 1.82 ± 0.50 (p = 0.048), 2.29 ± 0.47 and 2.51 ± 0.37 (p = 0.006), and 2.35 ± 0.57 and 2.36 ± 0.57 (p = 0.927), respectively. Intraclass correlation coefficient for the baseline conductivity of muscles, cartilages, and nerves were 0.89, 0.67, and 0.89, respectively. In conclusion, in vivo conductivity measurement of musculoskeletal tissues is feasible using MR-EPT. Conductivity of muscles and cartilages significantly changed with an overall increase after exercise.

Brain Tissue Conductivity in Focal Cerebral Ischemia

BACKGROUND

Cerebral ischemia leads to oxygen depletion with rapid breakdown of transmembrane transporters and subsequent impaired electrolyte haemostasis. Electric properties tomography (EPT) is a new contrast in MRI which delivers information on tissue electrical conductivity. In the clinical realm it has been mostly used for tumour mapping. Ischemic cerebral stroke is another promising but neglected application. It might deliver additional information on tissue viability and possible response to therapy.

AIM

The aim of this study was to demonstrate tissue conductivity in a rodent model of stroke. Further, we aimed to compare electric conductivity in ischemic and non-ischemic cerebral tissue.

MATERIALS AND METHODS

Two male Wistar rats were used in this study and were subjected to permanent MCAO. The animals were scanned in a 3 Tesla system (Philips Achieva/Best, the Netherlands) using a dedicated solenoid animal coil (Philips/Hamburg, Germany). In addition to diffusion weighted imaging (DWI), EPT was performed using a steady-state free-precession (SSFP) sequence (repetition time/echo time = 4.5/2.3 ms, measured voxel size = 0.6 × 0.6 × 1.2 mm3, flip angle = 38°, number of excitations = 4). From the transceive phase ϕ of these SSFP scans, conductivity σ was estimated by the equation σ = Δϕ/(2μ0ω) with Δ the Laplacian operator, μ0 the magnetic permeability, and ω the Larmor frequency. Subsequently, a median filter was applied, which was locally restricted to voxels with comparable signal magnitude.

RESULTS

The animals exhibited an infarct as demonstrated on DWI. Conductivity within the infarcted region was 60-70 % of the conductivity of not affected contralateral tissue (0.39 ± 0.07 S/m and 0.31 ± 0.14 S/m vs. 0.64 ± 0.15 S/m and 0.66 ± 0.16 S/m, respectively).

DISCUSSION

Infarcted tissue exhibited decreased conductivity. Further in-vivo experiments with examination of the influence of reperfusion status and temporal evolution of the infarcted areas should be conducted. Depiction of the ischemic penumbra and possibly subclassification of the DWI lesion still seems to be a fruitful target for further studies.

Mapping electric bulk conductivity in the human heart

PURPOSE

To explore the technical feasibility of mapping the electric bulk conductivity in the human heart, and to determine quantitative conductivity values of myocardium and blood from a small group of volunteers.

METHODS

Using a 3T MR system, 6 healthy male volunteers were measured. For all volunteers, a time-resolved 2D sequence over the cardiac cycle was applied (electrocardiogram [ECG]-triggered SSFP acquired in breath-hold). From these data, a dedicated, so-called "2D conductivity" has been derived in the framework of electrical properties tomography (EPT). To validate the concept of 2D conductivity, a static 3D sequence (ECG-triggered and respiratory-gated SSFP 3D whole heart acquisition, allowing the full 3D reconstruction of conductivity) as well as a Q-flow sequence (for investigating the relation between flow and reconstruction errors of the conductivity) have been applied for one of the volunteers.

RESULTS

For both, blood and myocardium, quantitative values of obtained 2D conductivity were approximately two-thirds of the obtained 3D conductivity, as expected from Maxwell's equations. Furthermore, the quantitative conductivity values agreed with corresponding literature values. Conductivity of left-ventricular blood volume showed characteristic over- and under-shooting at specific time points during the cardiac cycle for all volunteers investigated. This over- and under-shooting correlated with the phase pattern caused by blood flow into/out of the ventricle.

CONCLUSION

The study demonstrated the technical feasibility of cardiac conductivity measurements using standard MR systems and standard MR sequences, and therefore, may open new options for MR-based cardiac diagnosis.

High-frequency electrical properties tomography at 9.4T as a novel contrast mechanism for brain tumors

PURPOSE

To establish high-frequency magnetic resonance electrical properties tomography (MREPT) as a novel contrast mechanism for the assessment of glioblastomas using a rat brain tumor model.

METHODS

Six F98 intracranial tumor bearing rats were imaged longitudinally 8, 11 and 14 days after tumor cell inoculation. Conductivity and mean diffusivity maps were generated using MREPT and Diffusion Tensor Imaging. These maps were co-registered with T2 -weighted images and volumes of interests (VOIs) were segmented from the normal brain, ventricles, edema, viable tumor, tumor rim, and tumor core regions. Longitudinal changes in conductivity and mean diffusivity (MD) values were compared in these regions. A correlation analysis was also performed between conductivity and mean diffusivity values.

RESULTS

The conductivity of ventricles, edematous area and tumor regions (tumor rim, viable tumor, tumor core) was significantly higher (P < .01) compared to the contralateral cortex. The conductivity of the tumor increased over time while MD from the tumor did not change. A marginal positive correlation was noted between conductivity and MD values for tumor rim and viable tumor, whereas this correlation was negative for the tumor core.

CONCLUSION

We demonstrate a novel contrast mechanism based on ionic concentration and mobility, which may aid in providing complementary information to water diffusion in probing the microenvironment of brain tumors.

High-frequency conductivity at Larmor-frequency in human brain using moving local window multilayer perceptron neural network

Magnetic resonance electrical properties tomography (MREPT) aims to visualize the internal high-frequency conductivity distribution at Larmor frequency using the B1 transceive phase data. From the magnetic field perturbation by the electrical field associated with the radiofrequency (RF) magnetic field, the high-frequency conductivity and permittivity distributions inside the human brain have been reconstructed based on the Maxwell’s equation. Starting from the Maxwell’s equation, the complex permittivity can be described as a second order elliptic partial differential equation. The established reconstruction algorithms have focused on simplifying and/or regularizing the elliptic partial differential equation to reduce the noise artifact. Using the nonlinear relationship between the Maxwell’s equation, measured magnetic field, and conductivity distribution, we design a deep learning model to visualize the high-frequency conductivity in the brain, directly derived from measured magnetic flux density. The designed moving local window multi-layer perceptron (MLW-MLP) neural network by sliding local window consisting of neighboring voxels around each voxel predicts the high-frequency conductivity distribution in each local window. The designed MLW-MLP uses a family of multiple groups, consisting of the gradients and Laplacian of measured B1 phase data, as the input layer in a local window. The output layer of MLW-MLP returns the conductivity values in each local window. By taking a non-local mean filtering approach in the local window, we reconstruct a noise suppressed conductivity image while maintaining spatial resolution. To verify the proposed method, we used B1 phase datasets acquired from eight human subjects (five subjects for training procedure and three subjects for predicting the conductivity in the brain).

Improving phase-based conductivity reconstruction by means of deep learning-based denoising of B 1 + phase data for 3T MRI

PURPOSE

To denoise B 1 + phase using a deep learning method for phase-based in vivo electrical conductivity reconstruction in a 3T MR system.

METHODS

For B 1 + phase deep-learning denoising, a convolutional neural network (U-net) was chosen. Training was performed on data sets from 10 healthy volunteers. Input data were the real and imaginary components of single averaged spin-echo data (SNR = 45), which was used to approximate the B 1 + phase. For label data, multiple signal-averaged spin-echo data (SNR = 128) were used. Testing was performed on in silico and in vivo data. Reconstructed conductivity maps were derived using phase-based conductivity reconstructions. Additionally, we investigated the usability of the network to various SNR levels, imaging contrasts, and anatomical sites (ie, T₁ , T₂ , and proton density-weighted brain images and proton density-weighted breast images. In addition, conductivity reconstructions from deep learning-based denoised data were compared with conventional image filters, which were used for data denoising in electrical properties tomography (ie, the Gaussian filtering and the Savitzky-Golay filtering).

RESULTS

The proposed deep learning-based denoising approach showed improvement for B 1 + phase for both in silico and in vivo experiments with reduced quantitative error measures compared with other methods. Subsequently, this resulted in an improvement of reconstructed conductivity maps from the denoised B 1 + phase with deep learning.

CONCLUSION

The results suggest that the proposed approach can be used as an alternative preprocessing method to denoise B 1 + maps for phase-based conductivity reconstruction without relying on image filters or signal averaging.

Transceive phase corrected 2D contrast source inversion-electrical properties tomography

PURPOSE

To remove the necessity of the tranceive phase assumption for CSI-EPT and show electrical properties maps reconstructed from measured data obtained using a standard 3T birdcage body coil setup.

METHODS

The existing CSI-EPT algorithm is reformulated to use the transceive phase rather than relying on the transceive phase assumption. Furthermore, the radio frequency (RF)-shield is numerically implemented to accurately model the RF fields inside the MRI scanner. We verify that the reformulated two-dimensional (2D) CSI-EPT algorithm can reconstruct electrical properties maps given 2D electromagnetic simulations. Afterward, the algorithm is tested with three-dimensional (3D) FDTD simulations to investigate if the 2D CSI-EPT can retrieve the electrical properties for 3D RF fields. Finally, an MR experiment at 3T with a phantom is performed.

RESULTS

From the results of the 2D simulations, it is seen that CSI-EPT can reconstruct the electrical properties using MRI accessible quantities. For 3D simulations, it is observed that the electrical properties are underestimated, nonetheless, CSI-EPT has a lower standard deviation than the standard Helmholtz based methods. Finally, the first CSI-EPT reconstructions based on measured data are presented showing comparable accuracy and precision to reconstructions based on simulated data, and demonstrating the feasibility of CSI-EPT.

CONCLUSIONS

The CSI-EPT algorithm was rewritten to use MRI accessible quantities. This allows for CSI-EPT to fully exploit the benefits of the higher static magnetic field strengths with a standard quadrature birdcage coil setup.

Noncontrast-Enhanced MR-Based Conductivity Imaging for Breast Cancer Detection and Lesion Differentiation

BACKGROUND

There is increasing interest in noncontrast-enhanced MRI due to safety concerns for gadolinium contrast agents.

PURPOSE

To investigate the clinical feasibility of MR-based conductivity imaging for breast cancer detection and lesion differentiation.

STUDY TYPE

Prospective.

SUBJECTS

One hundred and ten women, with 112 known cancers and 17 benign lesions (biopsy-proven), scheduled for preoperative MRI.

FIELD STRENGTH/SEQUENCE

Non-fat-suppressed T2-weighted turbo spin-echo sequence (T2WI), dynamic contrast-enhanced MRI and diffusion-weighted imaging (DWI) at 3T.

ASSESSMENT

Cancer detectability on each imaging modality was qualitatively evaluated on a per-breast basis: the conductivity maps derived from T2WI were independently reviewed by three radiologists (R1-R3). T2WI, DWI, and pre-operative digital mammography were independently reviewed by three other radiologists (R4-R6). Conductivity and apparent diffusion coefficient (ADC) measurements (mean, minimum, and maximum) were performed for 112 cancers and 17 benign lesions independently by two radiologists (R1 and R2). Tumor size was measured from surgical specimens.

STATISTICAL TESTS

Cancer detection rates were compared using generalized estimating equations. Multivariable logistic regression analysis was performed to identify factors associated with cancer detectability. Discriminating ability of conductivity and ADC was evaluated by using the areas under the receiver operating characteristic curve (AUC).

RESULTS

Conductivity imaging showed lower cancer detection rates (20%-32%) compared to T2WI (62%-71%), DWI (85%-90%), and mammography (79%-88%) (all P < 0.05). Fatty breast on MRI (odds ratio = 11.8, P < 0.05) and invasive tumor size (odds ratio = 1.7, P < 0.05) were associated with cancer detectability of conductivity imaging. The maximum conductivity showed comparable ability to the mean ADC in discriminating between cancers and benign lesions (AUC = 0.67 [95% CI: 0.59, 0.75] vs. 0.84 [0.76, 0.90], P = 0.06 (R1); 0.65 [0.56, 0.73] vs. 0.82 [0.74, 0.88], P = 0.07 (R2)).

DATA CONCLUSION

Although conductivity imaging showed suboptimal performance in breast cancer detection, the quantitative measurement of conductivity showed the potential for lesion differentiation.

EVIDENCE LEVEL

1 TECHNICAL EFFICACY: Stage 2.

Low conductivity on electrical properties tomography demonstrates unique tumor habitats indicating progression in glioblastoma

OBJECTIVES

Tissue conductivity measurements made with electrical properties tomography (EPT) can be used to define temporal changes in tissue habitats on longitudinal multiparametric MRI. We aimed to demonstrate the added insights for identifying tumor habitats obtained by including EPT with diffusion- and perfusion-weighted MRI, and to evaluate the use of these tumor habitats for determining tumor treatment response in post-treatment glioblastoma.

METHODS

Tumor habitats were developed from EPT, diffusion-weighted, and perfusion-weighted MRI in 60 patients with glioblastoma who underwent concurrent chemoradiotherapy. Voxels from EPT, apparent diffusion coefficient (ADC), and cerebral blood volume (CBV) maps were clustered into habitats, and each habitat was serially examined to assess its temporal change. The usefulness of temporal changes in tumor habitats for diagnosing tumor progression and treatment-related change was investigated using logistic regression. The performance of significant predictors was measured using the area under the curve (AUC) from receiver-operating-characteristics analysis with 1000-fold bootstrapping.

RESULTS

Five tumor habitats were identified, and of these, the hypervascular cellular habitat (odds ratio [OR] 5.45; 95% CI, 1.75-31.42; p = .02), hypovascular low conductivity habitat (OR 2.00; 95% CI, 1.45-3.05; p < .001), and hypovascular intermediate habitat (OR 1.57; 95% CI, 1.18-2.30; p = .006) were predictive of tumor progression. Low EPT and low CBV reflected a unique hypovascular low conductivity habitat that showed the highest diagnostic performance (AUC 0.86; 95% CI, 0.76-0.96). The combined habitats showed high performance (AUC 0.90; 95% CI, 0.82-0.98) in the differentiation of tumor progression from treatment-related change.

CONCLUSION

EPT reveals low conductivity habitats that can improve the diagnosis of tumor progression in post-treatment glioblastoma.

KEY POINTS

• Electrical properties tomography (EPT) demonstrated lower conductivity in tumor progression than in treatment-related change. • EPT allowed identification of a unique hypovascular low conductivity habitat when combined with cerebral blood volume mapping. • Tumor habitats with a hypovascular low conductivity habitat, hypervascular cellular habitat, and hypovascular intermediate habitat yielded high diagnostic performance for diagnosing tumor progression.

Decomposition of high-frequency electrical conductivity into extracellular and intracellular compartments based on two-compartment model using low-to-high multi-b diffusion MRI

Background

As an object’s electrical passive property, the electrical conductivity is proportional to the mobility and concentration of charged carriers that reflect the brain micro-structures. The measured multi-b diffusion-weighted imaging (Mb-DWI) data by controlling the degree of applied diffusion weights can quantify the apparent mobility of water molecules within biological tissues. Without any external electrical stimulation, magnetic resonance electrical properties tomography (MREPT) techniques have successfully recovered the conductivity distribution at a Larmor-frequency.

Methods

This work provides a non-invasive method to decompose the high-frequency conductivity into the extracellular medium conductivity based on a two-compartment model using Mb-DWI. To separate the intra- and extracellular micro-structures from the recovered high-frequency conductivity, we include higher b-values DWI and apply the random decision forests to stably determine the micro-structural diffusion parameters.

Results

To demonstrate the proposed method, we conducted phantom and human experiments by comparing the results of reconstructed conductivity of extracellular medium and the conductivity in the intra-neurite and intra-cell body. The phantom and human experiments verify that the proposed method can recover the extracellular electrical properties from the high-frequency conductivity using a routine protocol sequence of MRI scan.

Conclusion

We have proposed a method to decompose the electrical properties in the extracellular, intra-neurite, and soma compartments from the high-frequency conductivity map, reconstructed by solving the electro-magnetic equation with measured B1 phase signals.

A robust electrical conductivity imaging method with total variation and wavelet regularization

PURPOSE

This study aims to develop and evaluate a robust conductivity imaging method that combines total variation and wavelet regularization to enhance the accuracy of conductivity maps.

THEORY AND METHODS

The proposed approach is based on a gradient-based method. The central equation is derived from Maxwell's equation and describes the relationship between conductivity and the transceive phase. A linear system equation is obtained via a finite-difference method and solved using a least-squares method. Total variation and wavelet transform regularization terms are added to the minimization problem and solved using the Split Bregman method to improve reconstruction stability. The proposed approach is compared with conventional and gradient-based methods. Numerical simulations are performed to validate the accuracy of the developed method, and the effects of noise are determined. Phantom and in vivo experiments are conducted at 3 T to verify the clinical applicability of the proposed method.

RESULTS

Numerical simulations show that the proposed method is more robust than other methods and can suppress the effects of noise. The quantitative conductivity value of the phantom experiment agrees with the measured value. The in vivo experiment results present a clear structure, and the conductivity value of the tumor region is significantly higher than that around healthy tissues.

CONCLUSION

The proposed electrical conductivity imaging method can improve the quality of conductivity reconstruction, and thus, has future clinical applications.

Magnetic resonance electrical property mapping at 21.1 T: a study of conductivity and permittivity in phantoms, ex vivo tissue and in vivo ischemia

Electrical properties (EP), namely conductivity and permittivity, can provide endogenous contrast for tissue characterization. Using electrical property tomography (EPT), maps of EP can be generated from conventional MRI data. This report investigates the feasibility and accuracy of EPT at 21.1 T for multiple RF coils and modes of operation using phantoms. Additionally, it demonstrates the EP of the in vivo rat brain with and without ischemia. Helmholtz-based EPT was implemented in its Full-form, which demands the complex [Formula: see text] field, and a simplified form requiring either just the [Formula: see text] field phase for conductivity or the [Formula: see text] field magnitude for permittivity. Experiments were conducted at 21.1 T using birdcage and saddle coils operated in linear or quadrature transceive mode, respectively. EPT approaches were evaluated using a phantom, ex and in vivo Sprague-Dawley rats under naïve conditions and ischemic stroke via transient middle cerebral artery occlusion. Different conductivity reconstruction approaches applied to the phantom displayed average errors of 12%-73% to the target acquired from dielectric probe measurements. Permittivity reconstructions showed higher agreement and an average 3%-8% error to the target depending on reconstruction approach. Conductivity and permittivity of ex and in vivo rodent brain were measured. Elevated EP in the ischemia region correlated with the increased sodium content and the influx of water intracellularly following ischemia in the lesion were detected. The Full-form technique generated from the linear birdcage provided the best accuracy for EP of the phantom. Phase-based conductivity and magnitude-based permittivity mapping provided reasonable estimates but also demonstrated the limitations of Helmholtz-based EPT at 21.1 T. Permittivity reconstruction was improved significantly over lower fields, suggesting a novel metric for in vivo brain studies. EPT applied to ischemic rat brain proved sensitivity to physiological changes, motivating the future application of more advanced reconstruction approaches.

Numerical Experiments on the Contrast Capability of Magnetic Resonance Electrical Property Tomography

Purpose:

Magnetic resonance electrical property tomography (MR EPT) is a technique for non-invasively obtaining the electric property (EP) distribution of biological tissues, with a promising potential for application in the early detection of tumors. However, the contrast capability (CC) of this technique has not been fully studied. This work aims to theoretically explore the CC for detecting the variation of EP values and the size of the imaging region.

Methods:

A simulation scheme was specifically designed to evaluate the CC of MR EPT. The simulation study has the advantage that the magnetic field can be accurately obtained. EP maps of the designed phantom embedded with target regions of designated various sizes and EPs were reconstructed using the homogeneous Helmholtz equation based on B1+ with different signal-to-noise ratios (SNRs). The CC was estimated by determining the smallest detectable EP contrast when the target region size was as large as the Laplacian kernel and the smallest detectable target region size when the EP contrast was the same as the difference between healthy and malignant tissues in the brain, based on the reconstructed EP maps.

Results:

Using noise free B1+, the smallest detectable contrast_σ and contrast_εr were 1% and 3%, respectively, and the smallest detectable target region size was 1 mesh unit (the base unit size used in the simulation) for conductivity and relative permittivity. The smallest detectable EP contrast and target region size were decreased as the B1+ SNR increased.

Conclusion:

The CC of MR EPT was related with the SNR of the magnetic field. A small EP contrast and size were necessary for detection at a high-SNR magnetic field. Obtaining a high-SNR magnetic field is important for improving the CC of MR EPT.

Integral MR-EPT With the Calculation of Coil Current Distributions

Many integral equation (IE)-based magnetic resonance electrical property tomography (MR-EPT) methods use unloaded incident radio-frequency (RF) fields from simulations that may not fully reflect the real situation and thus lead to reconstruction errors. To improve the accuracy of IE-based MR-EPT methods, a novel approach that enables the calculation of loaded coil current distributions and avoids the explicit use of incident RF fields is presented in this paper. In the proposed method, a hybrid source composed of the current source from the coil and the contrast source from the subject are introduced in the integral equations. Because the loaded coil current distributions can be extracted from the reconstructed hybrid source, the simulated incident RF fields are eliminated from the problem formulations. To improve the convergence performance, a modified conjugate gradient (CG) scheme was used where the gradients of the current source and contrast source were balanced through using different weighting parameters. The proposed method was verified through full-wave simulations at 9.4 and 7 T involving a heterogeneous ball and an anatomical head phantom. The numerical results indicated that by using the proposed method, an accurate coil current distributions and EPs profiles can be reconstructed and the desirable robustness against noise can also be achieved.

Opening a new window on MR-based Electrical Properties Tomography with deep learning

In the radiofrequency (RF) range, the electrical properties of tissues (EPs: conductivity and permittivity) are modulated by the ionic and water content, which change for pathological conditions. Information on tissues EPs can be used e.g. in oncology as a biomarker. The inability of MR-Electrical Properties Tomography techniques (MR-EPT) to accurately reconstruct tissue EPs by relating MR measurements of the transmit RF field to the EPs limits their clinical applicability. Instead of employing electromagnetic models posing strict requirements on the measured MRI quantities, we propose a data driven approach where the electrical properties reconstruction problem can be casted as a supervised deep learning task (DL-EPT). DL-EPT reconstructions for simulations and MR measurements at 3 Tesla on phantoms and human brains using a conditional generative adversarial network demonstrate high quality EPs reconstructions and greatly improved precision compared to conventional MR-EPT. The supervised learning approach leverages the strength of electromagnetic simulations, allowing circumvention of inaccessible MR electromagnetic quantities. Since DL-EPT is more noise-robust than MR-EPT, the requirements for MR acquisitions can be relaxed. This could be a major step forward to turn electrical properties tomography into a reliable biomarker where pathological conditions can be revealed and characterized by abnormalities in tissue electrical properties.

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

PURPOSES

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.

METHODS

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 B1 measurements 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.

RESULTS

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.

CONCLUSION

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.

Electrical properties tomography: Available contrast and reconstruction capabilities

MR-based electrical properties tomography converts the MRI transmit/receive RF field measurements to tissue electrical property maps through dedicated reconstruction algorithms. Recent reports showed that despite limitations, electrical properties tomography holds promise for generating additional contrast for tumor detection and patient-specific modeling of tissue-RF field interactions. This review summarizes the available tissue electrical property contrasts and compares them with the capabilities of the most commonly used electrical properties tomography reconstruction method. Future directions and prospects of clinical translation are discussed.

Redesign of the Laplacian kernel for improvements in conductivity imaging using MRI

PURPOSE

To develop an electrical property tomography reconstruction method that achieves improvements over standard method by redesigning the Laplacian kernel.

THEORY AND METHODS

A decomposition property of the governing PET equation shows the possibility of redesigning the Laplacian kernel for conductivity reconstruction. Hence, the discrete Laplacian operator used for electrical property tomography reconstruction is redesigned to have a Gaussian-like envelope, which enables manipulation of the spatial and spectral response. The characteristics of the proposed kernel are investigated through numerical simulations and in vivo brain experiments.

RESULTS

The proposed method reduces textured noise, which hampers observing features of the conductivity image. Furthermore, the proposed scheme can mitigate the propagation of local phase error such as flow-induced phase. By doing so, the proposed method can recover feature information in conductivity (or resistivity) images. Lastly, the proposed kernel can be extended to other electrical property tomography reconstructions, improving the quality of images.

CONCLUSION

An alternative design of the Laplacian kernel for conductivity imaging has been developed to mitigate the textured noise and the propagation of local phase artifact.

Dictionary-based electric properties tomography

PURPOSE

To develop and validate a new algorithm called "dictionary-based electric properties tomography" (dbEPT) for deriving tissue electric properties from measured B₁ maps.

METHODS

Inspired by Magnetic Resonance fingerprinting, dbEPT uses a dictionary of local patterns ("atoms") of B₁ maps and corresponding electric properties distributions, derived from electromagnetic field simulations. For reconstruction, a pattern from a measured B₁ map is compared with the B₁ atoms of the dictionary. The B₁ atom showing the best match with the measured B₁ pattern yields the optimum electric properties pattern that is chosen for reconstruction. Matching was performed through machine learning algorithms. Two dictionaries, using transmit and transceive phases, were evaluated. The spatial distribution of local matching distance between optimal atom and measured pattern yielded a reconstruction reliability map. The method was applied to reconstruct conductivity of 4 volunteers' brains. A conventional, Helmholtz-based Electric properties tomography (EPT) reconstruction was performed for reference. Noise performance was studied through phantom simulations.

RESULTS

Quantitative values of conductivity agree with literature values. Results of the 2 dictionaries exhibit only minor differences. Somewhat larger differences are visible between dbEPT and Helmholtz-based EPT. Quantified by the correlation between conductivity and anatomic images, dbEPT depicts brain details more clearly than Helmholtz-based EPT. Matching distance is minimal in homogeneous brain ventricles and increases with tissue heterogeneity. Central processing unit time was approximately 2 minutes per dictionary training and 3 minutes per brain conductivity reconstruction using standard hardware equipment.

CONCLUSION

A new, dictionary-based approach for reconstructing electric properties is presented. Its conductivity reconstruction is able to overcome the EPT transceive-phase problem.

3-D Contrast Source Inversion-Electrical Properties Tomography

Contrast source inversion-electrical properties tomography (CSI-EPT) is an iterative reconstruction method to retrieve the electrical properties (EPs) of tissues from magnetic resonance data. The method is based on integral representations of the electromagnetic field and has been shown to allow EP reconstructions of small structures as well as tissue boundaries with compelling accuracy. However, to date, the CSI-EPT has been implemented for 2-D configurations only, which limits its applicability. In this paper, a full 3-D extension of the CSI-EPT method is presented, to enable CSI-EPT to be applied to realistic 3-D scenarios. Here, we demonstrate a proof-of-principle of 3-D CSI-EPT and present the reconstructions of a 3-D abdominal body section and a 3-D head model using different settings of the transmit coil. Numerical results show that the full 3-D approach yields accurate reconstructions of the EPs, even at tissue boundaries and is most accurate in regions where the absolute value of the electric field is highest.

Electrical Properties Tomography Based on $B_{{1}}$ Maps in MRI: Principles, Applications, and Challenges

OBJECTIVE

The purpose is to provide a comprehensive review of the electrical properties tomography (EPT) technique, which was introduced to image the electrical properties (EPs) of tissue noninvasively by exploiting the measured field data of MRI.

METHODS

We reviewed the principle of EPT, reconstruction methods, biomedical applications such as tumor imaging, and existing challenges. As a key application of EPT, the estimation of specific absorption rate (SAR) due to MRI was discussed in the background of elevated risk of tissue heating at high field.

RESULTS AND CONCLUSION

Since the originally proposed local, homogeneous Helmholtz equation-based reconstruction algorithm, advanced EPT algorithms have emerged to address the challenges of EPT, including reconstruction error near tissue boundaries, noise sensitivity, inaccurate phase estimation, and elimination of the unmeasurable component, along with demonstrations of in vivo experiments. EPT techniques have been applied to investigate EPs of both healthy and pathological tissues in vivo and factors contributing to various EP value, including sodium, water content, etc. More studies are anticipated to consolidate the current findings. EPT-based subject-specific SAR estimation has led to in vivo demonstration of its feasibility and prediction of temperature increase of phantom during MRI scans merely using measured data.

SIGNIFICANCE

EPT has the advantage of high resolution and practical feasibility in a clinical setup for imaging the biomedically interesting EPs of tissue in the radiofrequency range. EPT-based SAR estimation is another promising topic for predicting tissue heating of individual subjects during a specific MRI scan.

Electric properties tomography: Biochemical, physical and technical background, evaluation and clinical applications

Electric properties tomography (EPT) derives the patient's electric properties, i.e. conductivity and permittivity, using standard magnetic resonance (MR) systems and standard MR sequences. Thus, EPT does not apply externally mounted electrodes, currents or radiofrequency (RF) probes, as is the case in competing techniques. EPT is quantitative MR, i.e. it yields absolute values of conductivity and permittivity. This review summarizes the physical equations underlying EPT, the corresponding basic and advanced reconstruction techniques and practical numerical aspects to realize these reconstruction techniques. MR sequences which map the field information required for EPT are outlined, and experiments to validate EPT in phantom and in vivo studies are described. Furthermore, the review describes the clinical findings which have been obtained with EPT so far, and attempts to understand the physiologic background of these findings.

Fast Spin Echo Imaging-Based Electric Property Tomography With K-Space Weighting via ${T}_{2}$ Relaxation (rEPT)

Magnetic resonance electrical property tomography (MREPT) is a technique used to extract the electrical properties of tissues (conductivity in particular) using a magnetic resonance imaging system. In this paper, we propose an improved data acquisition scheme for the electrical property tomography technique by utilizing T ₂ modulation in fast spin echo (FSE) imaging. This technique was motivated by a numerical analysis of conductivity reconstruction in the frequency domain; results reveal the spatial frequency-dependent noise texture of conventional methods. A data-acquisition scheme using the FSE sequence was formulated to concentrate the signal within a specific frequency range where notable noise amplification is observed in the conventional method. Through numerical studies, the performance of the proposed acquisition was investigated. Furthermore, a compensation scheme was applied to reduce quantification errors due to tissue-specific T ₂ modulation, which is inherent in FSE imaging. The technique was applied to phantom and in vivo experiments. Results showed improved conductivity contrasts in both experiments, as compared with conventional MREPT methods.

Noninvasive electrical conductivity measurement by MRI: a test of its validity and the electrical conductivity characteristics of glioma

OBJECTIVES

This study noninvasively examined the electrical conductivity (σ) characteristics of diffuse gliomas using MRI and tested its validity.

METHODS

MRI including a 3D steady-state free precession (3D SSFP) sequence was performed on 30 glioma patients. The σ maps were reconstructed from the phase images of the 3D SSFP sequence. The σ histogram metrics were extracted and compared among the contrast-enhanced (CET) and noncontrast-enhanced tumour components (NCET) and normal brain parenchyma (NP). Difference in tumour σ histogram metrics among tumour grades and correlation of σ metrics with tumour grades were tested. Validity of σ measurement using this technique was tested by correlating the mean tumour σ values measured using MRI with those measured ex vivo using a dielectric probe.

RESULTS

Several σ histogram metrics of CET and NCET of diffuse gliomas were significantly higher than NP (Bonferroni-corrected p ≤ .045). The maximum σ of NCET showed a moderate positive correlation with tumour grade (r = .571, Bonferroni-corrected p = .018). The mean tumour σ measured using MRI showed a moderate positive correlation with the σ measured ex vivo (r = .518, p = .040).

CONCLUSIONS

Tissue σ can be evaluated using MRI, incorporation of which may better characterise diffuse gliomas.

KEY POINTS

• This study tested the validity of noninvasive electrical conductivity measurements by MRI. • This study also evaluated the electrical conductivity characteristics of diffuse glioma. • Gliomas have higher electrical conductivity values than the normal brain parenchyma. • Noninvasive electrical conductivity measurement can be helpful for better characterisation of glioma.