Characterization of Nonlinearity and Dispersion in Tissue Impedance During High-Frequency Electroporation

Spatiotemporal estimations of temperature rise during electroporation treatments using a deep neural network

The nonthermal mechanism for irreversible electroporation has been paramount for treating tumors and cardiac tissue in anatomically sensitive areas, where there is concern about damage to nearby bowels, ducts, blood vessels, or nerves. However, Joule heating still occurs as a secondary effect of applying current through a resistive tissue and must be minimized to maintain the benefits of electroporation at high voltages. Numerous thermal mitigation protocols have been proposed to minimize temperature rise, but intraoperative temperature monitoring is still needed. We show that an accurate and robust temperature prediction AI model can be developed using estimated tissue properties (bulk and dynamic conductivity), known geometric properties (probe spacing), and easily measurable treatment parameters (applied voltage, current, and pulse number). We develop the 2-layer neural network on realistic 2D finite element model simulations with conditions encompassing most electroporation applications. Calculating feature contributions, we found that temperature prediction is mostly dependent on current and pulse number and show that the model remains accurate when incorrect tissue properties are intentionally used as input parameters. Lastly, we show that the model can predict temperature rise within ex vivo perfused porcine livers, with error <0.5 °C. This model, using easily acquired parameters, is shown to predict temperature rise in over 1000 unique test conditions with <1 °C error and no observable outliers. We believe the use of simple, readily available input parameters would allow this model to be incorporated in many already available electroporation systems for real-time temperature estimations.

Characterization of dispersion and anisotropic-conductivity in tissue model during electroporation pulses

Electroporation occurs when biological cells are exposed to intensive, short-duration pulses, which can be used to ablate biological tumor tissues. Based on the traditional numerical models, the isotropic conductivity model with the non-dispersion effect (ICND), the anisotropic conductivity model with the dispersion effect (ACD) is developed in this study. The second-order Debye function is introduced to manifest the dielectric relaxation effect, and the two-dimensional Cartesian conductivity matrix is applied to describe the anisotropy of the tissue conductivity during the electroporation process. A monopolar pulse and a monopolar burst are applied to the breast tumor model through the two-needle electrodes configuration. The results show that taking the dispersion into account can increase the total electroporated area more than 2.31%. Considering the conductivity anisotropy, the total electroporated area increases, but the irreversible electroporation (IRE) area decreases by more than 3.99%. The ACD model can achieve a larger electroporated area but a relatively smaller IRE area than those of the ICND model, and comparably minor maximum thermal damage is evaluated in the ACD model. Our model analyzes the effects of the dielectric dispersion and anisotropic conductivity of tissue, which have strong guiding significance for making the treatment planning before clinical practice.

A Theoretical Argument for Extended Interpulse Delays in Therapeutic High-Frequency Irreversible Electroporation Treatments

High-frequency irreversible electroporation (H-FIRE) is a tissue ablation modality employing bursts of electrical pulses in a positive phase-interphase delay (d1)-negative phase-interpulse delay (d2) pattern. Despite accumulating evidence suggesting the significance of these delays, their effects on therapeutic outcomes from clinically-relevant H-FIRE waveforms have not been studied extensively.

OBJECTIVE

We sought to determine whether modifications to the delays within H-FIRE bursts could yield a more desirable clinical outcome in terms of ablation volume versus extent of tissue excitation.

METHODS

We used a modified spatially extended nonlinear node (SENN) nerve fiber model to evaluate excitation thresholds for H-FIRE bursts with varying delays. We then calculated non-thermal tissue ablation, thermal damage, and excitation in a clinically relevant numerical model.

RESULTS

Excitation thresholds were maximized by shortening d1, and extension of d2 up to 1,000 μs increased excitation thresholds by at least 60% versus symmetric bursts. In the ablation model, long interpulse delays lowered the effective frequency of burst waveforms, modulating field redistribution and reducing heat production. Finally, we demonstrate mathematically that variable delays allow for increased voltages and larger ablations with similar extents of excitation as symmetric waveforms.

CONCLUSION

Interphase and interpulse delays play a significant role in outcomes resulting from H-FIRE treatment.

SIGNIFICANCE

Waveforms with short interphase delays (d1) and extended interpulse delays (d2) may improve therapeutic efficacy of H-FIRE as it emerges as a clinical tissue ablation modality.

Rapid Impedance Spectroscopy for Monitoring Tissue Impedance, Temperature, and Treatment Outcome During Electroporation-Based Therapies

OBJECTIVE

Electroporation-based therapies (EBTs) employ high voltage pulsed electric fields (PEFs) to permeabilize tumor tissue; this results in changes in electrical properties detectable using electrical impedance spectroscopy (EIS). Currently, commercial potentiostats for EIS are limited by impedance spectrum acquisition time (  ∼ 10 s); this timeframe is much larger than pulse periods used with EBTs (  ∼ 1 s). In this study, we utilize rapid EIS techniques to develop a methodology for characterizing electroporation (EP) and thermal effects associated with high-frequency irreversible EP (H-FIRE) in real-time by monitoring inter-burst impedance changes.

METHODS

A charge-balanced, bipolar rectangular chirp signal is proposed for rapid EIS. Validation of rapid EIS measurements against a commercial potentiostat was conducted in potato tissue using flat-plate electrodes and thereafter for the measurement of impedance changes throughout IRE treatment. Flat-plate electrodes were then utilized to uniformly heat potato tissue; throughout high-voltage H-FIRE treatment, low-voltage inter-burst impedance measurements were used to continually monitor impedance change and to identify a frequency at which thermal effects are delineated from EP effects.

RESULTS

Inter-burst impedance measurements (1.8 kHz - 4.93 MHz) were accomplished at 216 discrete frequencies. Impedance measurements at frequencies above  ∼ 1 MHz served to delineate thermal and EP effects in measured impedance.

CONCLUSION

We demonstrate rapid-capture ( 1 s) EIS which enables monitoring of inter-burst impedance in real-time. For the first time, we show impedance analysis at high frequencies can delineate thermal effects from EP effects in measured impedance.

SIGNIFICANCE

The proposed waveform demonstrates the potential to perform inter-burst EIS using PEFs compatible with existing pulse generator topologies.

Establishing an immunocompromised porcine model of human cancer for novel therapy development with pancreatic adenocarcinoma and irreversible electroporation

New therapies to treat pancreatic cancer are direly needed. However, efficacious interventions lack a strong preclinical model that can recapitulate patients’ anatomy and physiology. Likewise, the availability of human primary malignant tissue for ex vivo studies is limited. These are significant limitations in the biomedical device field. We have developed RAG2/IL2RG deficient pigs using CRISPR/Cas9 as a large animal model with the novel application of cancer xenograft studies of human pancreatic adenocarcinoma. In this proof-of-concept study, these pigs were successfully generated using on-demand genetic modifications in embryos, circumventing the need for breeding and husbandry. Human Panc01 cells injected subcutaneously into the ears of RAG2/IL2RG deficient pigs demonstrated 100% engraftment with growth rates similar to those typically observed in mouse models. Histopathology revealed no immune cell infiltration and tumor morphology was highly consistent with the mouse models. The electrical properties and response to irreversible electroporation of the tumor tissue were found to be similar to excised human pancreatic cancer tumors. The ample tumor tissue produced enabled improved accuracy and modeling of the electrical properties of tumor tissue. Together, this suggests that this model will be useful and capable of bridging the gap of translating therapies from the bench to clinical application.

Treatment Planning Optimization in Irreversible Electroporation for Complete Ablation of Variously Sized Cervical Tumors: A Numerical Study

Irreversible electroporation (IRE), a relatively new energy-based tumor ablation technology, has shown itself in the last decade to be able to safely ablate tumors with favorable clinical outcomes, yet little work has been done on optimizing the IRE protocol to variously sized tumors. Incomplete tumor ablation has been shown to be the main reason leading to the local recurrence and thus treatment failure. The goal of this study was to develop a general optimization approach to optimize the IRE protocol for cervical tumors in different sizes, while minimizing the damage to normal tissues. This kind of approach can lay a foundation for future personalized treatment of IRE. First, a statistical IRE cervical tumor death model was built using previous data in our group. Then, a multi-objective optimization problem model was built, in which the decision variables are five IRE-setting parameters, namely, the pulse strength (U), the length of active tip (H), the number of pulses delivered in one round between a pair of electrodes (A), the distance between electrodes (D), and the number of electrodes (N). The domains of the decision variables were determined based on the clinical experience. Finally, the problem model was solved by using nondominated sorting genetic algorithms II (NSGA-II) algorithm to give respective optimal protocol for three sizes of cervical tumors. Every protocol was assessed by the evaluation criterion established in the study to show the efficacy in a more straightforward way. The results of the study demonstrate this approach can theoretically provide the optimal IRE protocol for different sizes of tumors and may be generalizable to other types, sizes, and locations of tumors.

3-D H-Scan Ultrasound Imaging and Use of a Convolutional Neural Network for Scatterer Size Estimation

H-Scan ultrasound (US) is a new imaging technology that estimates the relative size of acoustic scattering objects and structures. The purpose of this study was to introduce a three-dimensional (3-D) H-scan US imaging approach for scatterer size estimation in volume space. Using a programmable research scanner (Vantage 256, Verasonics Inc, Kirkland, WA, USA) equipped with a custom volumetric imaging transducer (4 DL7, Vermon, Tours, France), raw radiofrequency (RF) data was collected for offline processing to generate H-scan US volumes. A deep convolutional neural network (CNN) was modified and used to achieve voxel mapping from the input H-scan US image to underlying scatterer size. Preliminary studies were conducted using homogeneous gelatin-based tissue-mimicking phantom materials embedded with acoustic scatterers of varying size (15 to 250 μm) and concentrations (0.1 to 1%). Two additional phantoms were embedded with 63 or 125 µm-sized microspheres and used to test CNN estimation accuracy. In vitro results indicate that 3-D H-scan US imaging can visualize the spatial distribution of acoustic scatterers of varying size at different concentrations (R² > 0.85, p < 0.03). The result of scatterer size estimation reveals that a CNN can achieve an average mapping accuracy of 93.3%. Overall, our preliminary in vitro findings reveal that 3-D H-scan US imaging allows the visualization of tissue scatterer patterns and incorporation of a CNN can be used to help estimate size of the acoustic scattering objects.

Starting a Fire Without Flame: The Induction of Cell Death and Inflammation in Electroporation-Based Tumor Ablation Strategies

New therapeutic strategies and paradigms are direly needed for the treatment of cancer. While the surgical removal of tumors is favored in most cancer treatment plans, resection options are often limited based on tumor localization. Over the last two decades, multiple tumor ablation strategies have emerged as promising stand-alone or combination therapeutic options for patients. These strategies are often employed to treat tumors in areas where surgical resection is not possible or where chemotherapeutics have proven ineffective. The type of cell death induced by the ablation modality is a critical aspect of therapeutic success that can impact the efficacy of the treatment and systemic anti-tumor immune system responses. Electroporation-based ablation technologies include electrochemotherapy, irreversible electroporation, and other modalities that rely on pulsed electric fields to create pores in cell membranes. These pores can either be reversible or irreversible depending on the electric field parameters and can induce cell death either alone or in combination with a therapeutic agent. However, there have been many controversial findings among these technologies as to the cell death type initiated, from apoptosis to pyroptosis. As cell death mechanisms can impact treatment side effects and efficacy, we review the main types of cell death induced by electroporation-based treatments and summarize the impact of these mechanisms on treatment response. We also discuss potential reasons behind the variability of findings such as the similarities between cell death pathways, differences between cell-types, and the variation in electric field strength across the treatment area.

Real-Time Impedance Monitoring During Electroporation Processes in Vegetal Tissue Using a High-Performance Generator

Classical application of electroporation is carried out by using fixed protocols that do not clearly assure the complete ablation of the desired tissue. Nowadays, new methods that pursue the control of the treatment by studying the change in impedance during the applied pulses as a function of the electric field are being developed. These types of control seek to carry out the treatment in the fastest way, decreasing undesired effects and treatment time while ensuring the proper tumour ablation. The objective of this research is to determine the state of the treatment by continuously monitoring the impedance by using a novel versatile high-voltage generator and sensor system. To study the impedance dynamics in real time, the use of pulses of reduced voltage, below the threshold of reversible electroporation, is tested to characterise the state-of-the-treatment without interfering with it. With this purpose, a generator that provides both low voltage for sense tissue changes and high voltage for irreversible electroporation (IRE) was developed. In conclusion, the characterisation of the effects of electroporation in vegetal tissue, combined with the real-time monitoring of the state-of-the-treatment, will enable the provision of safer and more effective treatments.

Experimental High-Frequency Irreversible Electroporation Using a Single-Needle Delivery Approach for Nonthermal Pancreatic Ablation In Vivo

PURPOSE

To investigate the feasibility of single-needle high-frequency irreversible electroporation (SN-HFIRE) to create reproducible tissue ablations in an in vivo pancreatic swine model.

MATERIALS AND METHODS

SN-HFIRE was performed in swine pancreas in vivo in the absence of intraoperative paralytics or cardiac synchronization using 3 different voltage waveforms (1-5-1, 2-5-2, and 5-5-5 [on-off-on times (μs)], n = 6/setting) with a total energized time of 100 μs per burst. At necropsy, ablation size/shape was determined. Immunohistochemistry was performed to quantify apoptosis using an anticleaved caspase-3 antibody. A numerical model was developed to determine lethal thresholds for each waveform in pancreas.

RESULTS

Mean tissue ablation time was 5.0 ± 0.2 minutes, and no cardiac abnormalities or muscle twitch was detected. Mean ablation area significantly increased with increasing pulse width (41.0 ± 5.1 mm² [range 32-66 mm²] vs 44 ± 2.1 mm² [range 38-56 mm²] vs 85.0 ± 7.0 mm² [range 63-155 mm²]; 1-5-1, 2-5-2, 5-5-5, respectively; p < 0.0002 5-5-5 vs 1-5-1 and 2-5-2). The majority of the ablation zone did not stain positive for cleaved caspase-3 (6.1 ± 2.8% [range 1.8-9.1%], 8.8 ± 1.3% [range 5.5-14.0%], and 11.0 ± 1.4% [range 7.1-14.2%] cleaved caspase-3 positive 1-5-1, 2-5-2, 5-5-5, respectively), with significantly more positive staining at the 5-5-5 pulse setting compared with 1-5-1 (p < 0.03). Numerical modeling determined a lethal threshold of 1114 ± 123 V/cm (1-5-1 waveform), 1039 ± 103 V/cm (2-5-2 waveform), and 693 ± 81 V/cm (5-5-5 waveform).

CONCLUSIONS

SN-HFIRE induces rapid, predictable ablations in pancreatic tissue in vivo without the need for intraoperative paralytics or cardiac synchronization.

Recent Advances in Electrochemotherapy

Electrochemotherapy is gaining recognition as an effective local therapy that uses systemically or intratumorally injected bleomycin or cisplatin with electroporation as a delivery system that brings drugs into the cells to exert their cytotoxic effects. Preclinical work is still ongoing, testing new drugs, seeking the best treatment combination with other treatment modalities, and exploring new sets of pulses for effective tissue electroporation. The applications of electrochemotherapy are being fully exploited in veterinary oncology, where electrochemotherapy, because of its simple execution, has a relatively good cost-benefit ratio and is used in the treatment of cutaneous tumors. In human oncology, electrochemotherapy is fully recognized as a local therapy for cutaneous tumors and metastases. Its effectiveness is being explored in combination with immunomodulatory drugs. However, the development of electrochemotherapy is directed into the treatment of deep-seated tumors with a percutaneous approach. Because of the vast number of reports, this review discusses the articles published in the past 5 years.

Temporal Characterization of Blood-Brain Barrier Disruption with High-Frequency Electroporation

Treatment of intracranial disorders suffers from the inability to accumulate therapeutic drug concentrations due to protection from the blood–brain barrier (BBB). Electroporation-based therapies have demonstrated the capability of permeating the BBB, but knowledge of the longevity of BBB disruption (BBBD) is limited. In this study, we quantify the temporal, high-frequency electroporation (HFE)-mediated BBBD in an in vivo healthy rat brain model. 40 male Fisher rats underwent HFE treatment; two blunt tipped monopolar electrodes were advanced into the brain and 200 bursts of HFE were delivered at a voltage-to-distance ratio of 600 V/cm. BBBD was verified with contrast enhanced T1W MRI (gadopentetate dimeglumine) and pathologically (Evans blue dye) at time points of 1, 24, 48, 72, and 96 h after HFE. Contrast enhanced T1W scans demonstrated BBBD for 1 to 72 h after HFE but intact BBB at 96 h. Histologically, tissue damage was restricted to electrode insertion tracks. BBBD was induced with minimal muscle contractions and minimal cell death attributed to HFE. Numerical modeling indicated that brief BBBD was induced with low magnitude electric fields, and BBBD duration increased with field strength. These data suggest the spatiotemporal characteristics of HFE-mediated BBBD may be modulated with the locally applied electric field.

Simplified Non-Thermal Tissue Ablation With a Single Insertion Device Enabled by Bipolar High-Frequency Pulses

OBJECTIVE

To demonstrate the feasibility of a single electrode and grounding pad approach for delivering high frequency irreversible electroporation treatments (H-FIRE) in in-vivo hepatic tissue.

METHODS

Ablations were created in porcine liver under surgical anesthesia by adminstereing high frequency bursts of 0.5-5.0 μs pulses with amplitudes between 1.1-1.7 kV in the absence of cardiac synchronization or intraoperative paralytics. Finite element simulations were used to determine the electric field strength associated with the ablation margins (E_Lethal) and predict the ablations feasible with next generation electronics.

RESULTS

All animals survived the procedures for the protocol duration without adverse events. E_Lethal of 2550, 1650, and 875 V/cm were found for treatments consisting of 100x bursts containing 0.5 μs pulses and 25, 50, and 75 μs of energized-time per burst, respectively. Treatments with 1 μs pulses consisting of 100 bursts with 100 μs energized-time per burst resulted in E_Lethal of 650 V/cm.

CONCLUSION

A single electrode and grounding pad approach was successfully used to create ablations in hepatic tissue. This technique has the potential to reduce challenges associated with placing multiple electrodes in anatomically challenging environments.

SIGNIFICANCE

H-FIRE is an in situ tumor ablation approach in which electrodes are placed within or around a targeted region to deliver high voltage electrical pulses. Electric fields generated around the electrodes induce irrecoverable cell membrane damage leading to predictable cell death in the relative absence of thermal damage. The sparing of architectural integrity means H-FIRE offers potential advantages compared to thermal ablation modalities for ablating tumors near critical structures.

High-Frequency Irreversible Electroporation Using 5,000-V Waveforms to Create Reproducible 2- and 4-cm Ablation Zones-A Laboratory Investigation Using Mechanically Perfused Liver

PURPOSE

To investigate if high-frequency irreversible electroporation (H-FIRE) treatments can be delivered at higher voltages and with greater energy delivery rates than currently implemented in clinical irreversible electroporation protocols.

MATERIALS AND METHODS

Treatments using 3,000 V and 5,000 V were administered to mechanically perfused ex vivo porcine liver via a single applicator and grounding pad (A+GP) as well as a 4-applicator array (4AA). Integrated energized times (IET) 0.01-0.08 seconds and energy delivery rates 25-300 μs/s were investigated. Organs were preserved at 4°C for 10-15 hours before sectioning and gross analysis using a metabolic stain to identify the size and shape of ablation zones.

RESULTS

A+GP ablations measured between 1.6 cm and 2.2 cm, which did not increase when IET was increased from 0.02 seconds to 0.08 seconds (P > .055; range, 1.9-2.1 cm). Changes in tissue color and texture consistent with thermal damage were observed for treatments with energy delivery rates 50-300 μs/s, but not for treatments delivered at 25 μs/s. Use of the 4AA with a 3-cm applicator spacing resulted in ablations measuring 4.4-4.9 cm with energy delivery times of 7-80 minutes.

CONCLUSIONS

H-FIRE treatments can rapidly and reproducibly create 2-cm ablations using an A+GP configuration. Treatments without thermal injury were produced at the expense of extended treatment times. More rapid treatments resulted in ablations with varying degrees of thermal injury within the H-FIRE ablation zone. Production of 4-cm ablations is possible using a 4AA.

Successful Tumor Electrochemotherapy Using Sine Waves

OBJECTIVE

The purpose of this work is to assess the ability of sine waves to perform electrochemotherapy (ECT) and to study the dependence of the frequency of the applied sine wave on the treatment efficacy.

METHODS

A subcutaneous tumor model in mice was used, and the electric field was delivered in combination with bleomycin. Sinusoidal electric fields of different frequencies, amplitudes, and durations were compared to square waves. Computer simulations were additionally performed.

RESULTS

The results confirmed the ability of a sinusoidal electric field to obtain successful ECT responses. A strong dependence on frequency was obtained. The efficacy of the treatment decreased when the frequency of the sine waves was increased. At low sinusoidal frequency, the efficacy of the treatment is very similar to that obtained with a square wave. The collateral effects such as skin burns and muscle contractions decreased for the highest frequency assayed.

CONCLUSION

The use of sine wave burst represents a feasible option for the treatment of cancer by ECT.

SIGNIFICANCE

These results could have important implications for the treatment of cancer in the clinical world where ECT is performed with dc square pulses.

Histopathological and Ultrastructural Changes after Electroporation in Pig Liver Using Parallel-Plate Electrodes and High-Performance Generator

Irreversible electroporation (IRE) has gained attention as a new non-thermal therapy for ablation with important benefits in terms of homogeneous treatment and fast recovery. In this study, a new concept of high voltage generator is used, enabling irreversible electroporation treatment in large tissue volume using parallel plates. Unlike currently available generators, the proposed versatile structure enables delivering high-voltage high-current pulses. To obtain homogeneous results, 3-cm parallel-plates electrodes have also been designed and implemented. IRE ablation was performed on six female pigs at 2000 V/cm electric field, and the results were analysed after sacrifice three hours, three days and seven days after ablation. Histopathological and ultrastructural studies, including transmission and scanning electron microscopy, were carried out. The developed high-voltage generator has proved to be effective for homogeneous IRE treatment using parallel plates. The destruction of the membrane of the hepatocytes and the alterations of the membranes of the cellular organelles seem incompatible with cell death by apoptosis. Although endothelial cells also die with electroporation, the maintenance of vascular scaffold allows repairing processes to begin from the third day after IRE as long as the blood flow has not been interrupted. This study has opened new direction for IRE using high performance generators and highlighted the importance of taking into account ultrastructural changes after IRE by using electron microscopy analysis.

Tumor Ablation Enhancement by Combining Radiofrequency Ablation and Irreversible Electroporation: An In Vitro 3D Tumor Study

We hypothesized and demonstrated for the first time that significant tumor ablation enhancement can be achieved by combining radiofrequency ablation (RFA) and irreversible electroporation (IRE) using a 3D cervical cancer cell model. Three RFA (43, 50, and 60 °C for 2 min) and IRE protocols (350, 700, and 1050 V/cm) were used to study the combining effect in the 3D tumor cell model. The in vitro experiment showed that both RFA enhanced IRE and IRE enhanced RFA can lead to a significant increase in the size of the ablation zone compared to IRE and RFA alone. It was also noted that the sequence of applying ablation energy (RFA → RE or IRE → RFA) affected the efficacy of tumor ablation enhancement. The electrical conductivity of 3D tumor was found to be increased after preliminary RFA or IRE treatment. This increase in tumor conductivity may explain the enhancement of tumor ablation. Another explanation might be that there is repeat injury to the transitional zone of the first treatment by the second one. The promising results achieved in the study can provide us useful clues about the treatment of large tumors abutting large vessels or bile ducts.

Ablation outcome of irreversible electroporation on potato monitored by impedance spectrum under multi-electrode system

BACKGROUND

Irreversible electroporation (IRE) therapy relies on pulsed electric fields to non-thermally ablate cancerous tissue. Methods for evaluating IRE ablation in situ are critical to assessing treatment outcome. Analyzing changes in tissue impedance caused by electroporation has been proposed as a method for quantifying IRE ablation. In this paper, we assess the hypothesis that irreversible electroporation ablation outcome can be monitored using the impedance change measured by the electrode pairs not in use, getting more information about the ablation size in different directions.

METHODS

Using a square four-electrode configuration, the two diagonal electrodes were used to electroporate potato tissue. Next, the impedance changes, before and after treatment, were measured from different electrode pairs and the impedance information was extracted by fitting the data to an equivalent circuit model. Finally, we correlated the change of impedance from various electrode pairs to the ablation geometry through the use of fitted functions; then these functions were used to predict the ablation size and compared to the numerical simulation results.

RESULTS

The change in impedance from the electrodes used to apply pulses is larger and has higher deviation than the other electrode pairs. The ablation size and the change in resistance in the circuit model correlate with various linear functions. The coefficients of determination for the three functions are 0.8121, 0.8188 and 0.8691, respectively, showing satisfactory agreement. The functions can well predict the ablation size under different pulse numbers, and in some directions it did even better than the numerical simulation method, which used different electric field thresholds for different pulse numbers.

CONCLUSIONS

The relative change in tissue impedance measured from the non-energized electrodes can be used to assess ablation size during treatment with IRE according to linear functions.