Saline-Infused Radiofrequency Ablation: A Review on the Key Factors for a Safe and Reliable Tumour Treatment

Radiofrequency ablation (RFA) combined with saline infusion into tissue is a promising technique to ablate larger tumours. Nevertheless, the application of saline-infused RFA remains at clinical trials due to the contradictory findings as a result of the inconsistencies in experimental procedures. These inconsistencies not only magnify the number of factors to consider during the treatment, but also obscure the understanding of the role of saline in enlarging the coagulation zone. Consequently, this can result in major complications, which includes unwanted thermal damages to adjacent tissues and also incomplete ablation of the tumour. This review aims to identify the key factors of saline responsible for enlarging the coagulation zone during saline-infused RFA, and provide a proper understanding on their effects that is supported with findings from computational studies to ensure a safe and reliable cancer treatment.

Enlarging the thermal coagulation volume during thermochemical ablation with alternating acid-base injection by shortening the injection interval: A computational study

BACKGROUND AND OBJECTIVES

Thermochemical ablation (TCA) is a cancer treatment that utilises the heat released from the neutralisation of acid and base to raise tissue temperature to levels sufficient to induce thermal coagulation. Computational studies have demonstrated that the coagulation volume produced by sequential injection is smaller than that with simultaneous injection. By injecting the reagents in an ensuing manner, the region of contact between acid and base is limited to a thin contact layer sandwiched between the distribution of acid and base. It is hypothesised that increasing the frequency of acid-base injections into the tissue by shortening the injection interval for each reagent can increase the effective area of contact between acid and base, thereby intensifying neutralisation and the exothermic heat released into the tissue.

METHODS

To verify this hypothesis, a computational model was developed to simulate the thermochemical processes involved during TCA with sequential injection. Four major processes that take place during TCA were considered, i.e., the flow of acid and base, their neutralisation, the release of exothermic heat and the formation of thermal damage inside the tissue. Equimolar acid and base at 7.5 M was injected into the tissue intermittently. Six injection intervals, namely 3, 6, 15, 20, 30 and 60 s were investigated.

RESULTS

Shortening of the injection interval led to the enlargement of coagulation volume. If one considers only the coagulation volume as the determining factor, then a 15 s injection interval was found to be optimum. Conversely, if one places priority on safety, then a 3 s injection interval would result in the lowest amount of reagent residue inside the tissue after treatment. With a 3 s injection interval, the coagulation volume was found to be larger than that of simultaneous injection with the same treatment parameters. Not only that, the volume also surpassed that of radiofrequency ablation (RFA); a conventional thermal ablation technique commonly used for liver cancer treatment.

CONCLUSION

The numerical results verified the hypothesis that shortening the injection interval will lead to the formation of larger thermal coagulation zone during TCA with sequential injection. More importantly, a 3 s injection interval was found to be optimum for both efficacy (large coagulation volume) and safety (least amount of reagent residue).

The role of shear viscosity as a biomarker for improving chronic kidney disease detection using shear wave elastography: A computational study using a validated finite element model

The application of ultrasound shear wave elastography for detecting chronic kidney disease, namely renal fibrosis, has been widely studied. A good correlation between tissue Young's modulus and the degree of renal impairment has been established. However, the current limitation of this imaging modality pertains to the linear elastic assumption used in quantifying the stiffness of renal tissue in commercial shear wave elastography systems. As such, when underlying medical conditions such as acquired cystic kidney disease, which may potentially influence the viscous component of renal tissue, is present concurrently with renal fibrosis, the accuracy of the imaging modality in detecting chronic kidney disease may be affected. The findings in this study demonstrate that quantifying the stiffness of linear viscoelastic tissue using an approach similar to those implemented in commercial shear wave elastography systems led to percentage errors as high as 87%. The findings presented indicate that use of shear viscosity to detect changes in renal impairment led to a reduction in percentage error to values as low as 0.3%. For cases in which renal tissue was affected by multiple medical conditions, shear viscosity was found to be a good indicator in gauging the reliability of the Young's modulus (quantified through a shear wave dispersion analysis) in detecting chronic kidney disease. The findings show that percentage error in stiffness quantification can be reduced to as low as 0.6%. The present study demonstrates the potential use of renal shear viscosity as a biomarker to improve the detection of chronic kidney disease.

A computational framework for the multiphysics simulation of microbubble-mediated sonothrombolysis using a forward-viewing intravascular transducer

Sonothrombolysis is a technique that utilises ultrasound waves to excite microbubbles surrounding a clot. Clot lysis is achieved through mechanical damage induced by acoustic cavitation and through local clot displacement induced by acoustic radiation force (ARF). Despite the potential of microbubble-mediated sonothrombolysis, the selection of the optimal ultrasound and microbubble parameters remains a challenge. Existing experimental studies are not able to provide a complete picture of how ultrasound and microbubble characteristics influence the outcome of sonothrombolysis. Likewise, computational studies have not been applied in detail in the context of sonothrombolysis. Hence, the effect of interaction between the bubble dynamics and acoustic propagation on the acoustic streaming and clot deformation remains unclear. In the present study, we report for the first time the computational framework that couples the bubble dynamic phenomena with the acoustic propagation in a bubbly medium to simulate microbubble-mediated sonothrombolysis using a forward-viewing transducer. The computational framework was used to investigate the effects of ultrasound properties (pressure and frequency) and microbubble characteristics (radius and concentration) on the outcome of sonothrombolysis. Four major findings were obtained from the simulation results: (i) ultrasound pressure plays the most dominant role over all the other parameters in affecting the bubble dynamics, acoustic attenuation, ARF, acoustic streaming, and clot displacement, (ii) smaller microbubbles could contribute to a more violent oscillation and improve the ARF simultaneously when they are stimulated at higher ultrasound pressure, (iii) higher microbubbles concentration increases the ARF, and (iv) the effect of ultrasound frequency on acoustic attenuation is dependent on the ultrasound pressure. These results may provide fundamental insight that is crucial in bringing sonothrombolysis closer to clinical implementation.

Gold nanorods assisted photothermal therapy of bladder cancer in mice: A computational study on the effects of gold nanorods distribution at the centre, periphery, and surface of bladder cancer

BACKGROUND AND OBJECTIVES

Gold nanorod-assisted photothermal therapy (GNR-PTT) is a cancer treatment whereby GNRs incorporated into the tumour act as photo-absorbers to elevate the thermal destruction effect. In the case of bladder, there are few possible routes to target the tumour with GNRs, namely peri/intra-tumoural injection and intravesical instillation of GNRs. These two approaches lead to different GNR distribution inside the tumour and can affect the treatment outcome.

METHODOLOGY

The present study investigates the effects of heterogeneous GNR distribution in a typical setup of GNR-PTT. Three cases were considered. Case 1 considered the GNRs at the tumour centre, while Case 2 represents a hypothetical scenario where GNRs are distributed at the tumour periphery; these two cases represent intratumoural accumulation with different degree of GNR spread inside the tumour. Case 3 is achieved when GNRs target the exposed tumoural surface that is invading the bladder wall, when they are delivered by intravesical instillation.

RESULTS

Results indicate that for a laser power of 0.6 W and GNR volume fraction of 0.01%, Case 2 and 3 were successful in achieving complete tumour eradication after 330 and 470 s of laser irradiation, respectively. Case 1 failed to form complete tumour damage when the GNRs are concentrated at the tumour centre but managed to produce complete tumour damage if the spread of GNRs is wider. Results from Case 2 also demonstrated a different heating profile from Case 1, suggesting that thermal ablation during GNR-PTT is dependant on the GNRs distribution inside the tumour. Case 3 shows similar results to Case 2 whereby gradual but uniform heating is observed. Cases 2 and 3 show that uniformly heating the tumour can reduce damage to the surrounding tissues.

CONCLUSIONS

Different GNR distribution associated with the different methods of introducing GNRs to the bladder during GNR-PTT affect the treatment outcome of bladder cancer in mice. Insufficient spreading during intratumoural injection of GNRs can render the treatment ineffective, while administered via intravesical instillation. GNR distribution achieved through intravesical instillation present some advantages over intratumoural injection and is worthy of further exploration.

An in silico assessment on the potential of using saline infusion to overcome non-confluent coagulation zone during two-probe, no-touch bipolar radiofrequency ablation of liver cancer

No-touch bipolar radiofrequency ablation (bRFA) is known to produce incomplete tumour ablation with a 'butterfly-shaped' coagulation zone when the interelectrode distance exceeds a certain threshold. Although non-confluent coagulation zone can be avoided by not implementing the no-touch mode, doing so exposes the patient to the risk of tumour track seeding. The present study investigates if prior infusion of saline into the tissue can overcome the issues of non-confluent or butterfly-shaped coagulation. A computational modelling approach based on the finite element method was carried out. A two-compartment model comprising the tumour that is surrounded by healthy liver tissue was developed. Three cases were considered; i) saline infusion into the tumour centre; ii) one-sided saline infusion outside the tumour; and iii) two-sided saline infusion outside the tumour. For each case, three different saline volumes were considered, i.e. 6, 14 and 22 ml. Saline concentration was set to 15% w/v. Numerical results showed that saline infusion into the tumour centre can overcome the butterfly-shaped coagulation only if the infusion volume is sufficient. On the other hand, one-sided infusion outside the tumour did not overcome this. Two-sided infusion outside the tumour produced confluent coagulation zone with the largest volume. Results obtained from the present study suggest that saline infusion, when carried out correctly, can be used to effectively eradicate liver cancer. This presents a practical solution to address non-confluent coagulation zone typical of that during two-probe bRFA treatment.

An in silico derived dosage and administration guide for effective thermochemical ablation of biological tissues with simultaneous injection of acid and base

BACKGROUND AND OBJECTIVES

Thermochemical ablation (TCA) is a thermal ablation technique involving the injection of acid and base, either sequentially or simultaneously, into the target tissue. TCA remains at the conceptual stage with existing studies unable to provide recommendations on the optimum injection rate, and reagent concentration and volume. Limitations in current experimental methodology have prevented proper elucidation of the thermochemical processes inside the tissue during TCA. Nevertheless, the computational TCA framework developed recently by Mak et al. [Mak et al., Computers in Biology and Medicine, 2022, 145:105494] has opened new avenues in the development of TCA. Specifically, a recommended safe dosage is imperative in driving TCA research beyond the conceptual stage.

METHODS

The aforesaid computational TCA framework for sequential injection was applied and adapted to simulate TCA with simultaneous injection of acid and base at equimolar and equivolume. The developed framework, which describes the flow of acid and base, their neutralisation, the rise in tissue temperature and the formation of thermal damage, was solved numerically using the finite element method. The framework will be used to investigate the effects of injection rate, reagent concentration, volume and type (weak/strong acid-base combination) on temperature rise and thermal coagulation formation.

RESULTS

A higher injection rate resulted in higher temperature rise and larger thermal coagulation. Reagent concentration of 7500 mol/m³ was found to be optimum in producing considerable thermal coagulation without the risk of tissue overheating. Thermal coagulation volume was found to be consistently larger than the total volume of acid and base injected into the tissue, which is beneficial as it reduces the risk of chemical burn injury. Three multivariate second-order polynomials that express the targeted coagulation volume as functions of injection rate and reagent volume, for the weak-weak, weak-strong and strong-strong acid-base combinations were also derived based on the simulated data.

CONCLUSIONS

A guideline for a safe and effective implementation of TCA with simultaneous injection of acid and base was recommended based on the numerical results of the computational model developed. The guideline correlates the coagulation volume with the reagent volume and injection rate, and may be used by clinicians in determining the safe dosage of reagents and optimum injection rate to achieve a desired thermal coagulation volume during TCA.

A computational framework to simulate the thermochemical process during thermochemical ablation of biological tissues

Thermochemical ablation (TCA) is a thermal ablation therapy that utilises heat released from acid-base neutralisation reaction to destroy tumours. This procedure is a promising low-cost solution to existing thermal ablation treatments such as radiofrequency ablation (RFA) and microwave ablation (MWA). Studies have demonstrated that TCA can produce thermal damage that is on par with RFA and MWA when employed properly. Nevertheless, TCA remains a concept that is tested only in a few animal trials due to the risks involved as the result of uncontrolled infusion and incomplete acid-base reaction. In this study, a computational framework that simulates the thermochemical process of TCA is developed. The proposed framework consists of three physics, namely chemical flow, neutralisation reaction and heat transfer. An important parameter in the TCA framework is the neutralisation reaction rate constant, which has values in the order of 10⁸ m³/(mol⋅s). The present study will demonstrate that since the rate constant impacts only the rate and direction of the reaction but has little influence on the extent of reaction, it is possible to replicate the thermochemical process of TCA by employing significantly smaller values of rate constant that are numerically tractable. Comparisons of the numerical results against experimental studies from the literature supports this. The aim of this framework is for researchers to advance and develop TCA to gain an in-depth understanding of the fundamental mechanisms of TCA and to develop a safe treatment protocol of TCA in the hope of advancing TCA into clinical trials.

How does saline backflow affect the treatment of saline-infused radiofrequency ablation?

BACKGROUND AND OBJECTIVE

Saline infusion is applied together with radiofrequency ablation (RFA) to enlarge the ablation zone. However, one of the issues with saline-infused RFA is backflow, which spreads saline along the insertion track. This raises the concern of not only thermally ablating the tissue within the backflow region, but also the loss of saline from the targeted tissue, which may affect the treatment efficacy.

METHODS

In the present study, 2D axisymmetric models were developed to investigate how saline backflow influence saline-infused RFA and whether the aforementioned concerns are warranted. Saline-infused RFA was described using the dual porosity-Joule heating model. The hydrodynamics of backflow was described using Poiseuille law by assuming the flow to be similar to that in a thin annulus. Backflow lengths of 3, 4.5, 6 and 9 cm were considered.

RESULTS

Results showed that there is no concern of thermally ablating the tissue in the backflow region. This is due to the Joule heating being inversely proportional to distance from the electrode to the fourth power. Results also indicated that larger backflow lengths led to larger growth of thermal damage along the backflow region and greater decrease in coagulation volume. Hence, backflow needs to be controlled to ensure an effective treatment of saline-infused RFA.

CONCLUSIONS

There is no risk of ablating tissues around the needle insertion track due to backflow. Instead, the risk of underablation as a result of the loss of saline due to backflow was found to be of greater concern.

A numerical study to investigate the effects of tumour position on the treatment of bladder cancer in mice using gold nanorods assisted photothermal ablation

Gold nanorods assisted photothermal therapy (GNR-PTT) is a new cancer treatment technique that has shown promising potential for bladder cancer treatment. The position of the bladder cancer at different locations along the bladder wall lining can potentially affect the treatment efficacy since laser is irradiated externally from the skin surface. The present study investigates the efficacy of GNR-PTT in the treatment of bladder cancer in mice for tumours growing at three different locations on the bladder, i.e., Case 1: closest to skin surface, Case 2: at the bottom half of the bladder, and Case 3: at the side of the bladder. Investigations were carried out numerically using an experimentally validated framework for optical-thermal simulations. An in-silico approach was adopted due to the flexibility in placing the tumour at a desired location along the bladder lining. Results indicate that for the treatment parameters considered (laser power 0.3 W, GNR volume fraction 0.01% v/v), only Case 1 can be used for an effective GNR-PTT. No damage to the tumour was observed in Cases 2 and 3. Analysis of the thermo-physiological responses showed that the effectiveness of GNR-PTT in treating bladder cancer depends not only on the depth of the tumour from the skin surface, but also on the type of tissue that the laser must pass through before reaching the tumour. In addition, the results are reliant on GNRs with a diameter of 10 nm and an aspect ratio of 3.8 - tuned to exhibit peak absorption for the chosen laser wavelength. Results from the present study can be used to highlight the potential for using GNR-PTT for treatment of human bladder cancer. It appears that Cases 2 and 3 suggest that GNR-PTT, where the laser passes through the skin to reach the bladder, may be unfeasible in humans. While this study shows the feasibility of using GNRs for photothermal ablation of bladder cancer, it also identifies the current limitations needed to be overcome for an effective clinical application in the bladder cancer patients.

Shear Wave Elastography: A Review on the Confounding Factors and Their Potential Mitigation in Detecting Chronic Kidney Disease

Early detection of chronic kidney disease is important to prevent progression of irreversible kidney damage, reducing the need for renal transplantation. Shear wave elastography is ideal as a quantitative imaging modality to detect chronic kidney disease because of its non-invasive nature, low cost and portability, making it highly accessible. However, the complexity of the kidney architecture and its tissue properties give rise to various confounding factors that affect the reliability of shear wave elastography in detecting chronic kidney disease, thus limiting its application to clinical trials. The objective of this review is to highlight the confounding factors presented by the complex properties of the kidney, in addition to outlining potential mitigation strategies, along with the prospect of increasing the versatility and reliability of shear wave elastography in detecting chronic kidney disease.

Bipolar radiofrequency ablation treatment of liver cancer employing monopolar needles: A comprehensive investigation on the efficacy of time-based switching

Radiofrequency ablation (RFA) is a thermal ablative treatment method that is commonly used to treat liver cancer. However, the thermal coagulation zone generated using the conventional RFA system can only successfully treat tumours up to 3 cm in diameter. Switching bipolar RFA has been proposed as a way to increase the thermal coagulation zone. Presently, the understanding of the underlying thermal processes that takes place during switching bipolar RFA remains limited. Hence, the objective of this study is to provide a comprehensive understanding on the thermal ablative effects of time-based switching bipolar RFA on liver tissue. Five switch intervals, namely 50, 100, 150, 200 and 300 s were investigated using a two-compartment 3D finite element model. The study was performed using two pairs of RF electrodes in a four-probe configuration, where the electrodes were alternated based on their respective switch interval. The physics employed in the present study were verified against experimental data from the literature. Results obtained show that using a shorter switch interval can improve the homogeneity of temperature distribution within the tissue and increase the rate of temperature rise by delaying the occurrence of roll-off. The coagulation volume obtained was the largest using switch interval of 50 s, followed by 100, 150, 200 and 300 s. The present study demonstrated that the transient thermal response of switching bipolar RFA can be improved by using shorter switch intervals.

Thermal and thermal damage responses during switching bipolar radiofrequency ablation employing bipolar needles: A computational study on the effects of different electrode configuration, input voltage and ablation duration

Recent studies have demonstrated the effectiveness of switching bipolar radiofrequency ablation (bRFA) in treating liver cancer. Nevertheless, the clinical use of the treatment remains less common than conventional monopolar RFA - likely due to the lack of understanding of how the tissues respond thermally to the switching effect. The problem is exacerbated by the numerous possible switching combinations when bRFA is performed using bipolar needles, thus making theoretical deduction and experimental studies difficult. This article addresses this issue via computational modelling by examining if significant variation in the treatment outcome exists amongst six different electrode configurations defined by the X-, C-, U-, N-, Z- and O-models. Results indicated that the tissue thermal and thermal damage responses varied depending on the electrode configuration and the operating conditions (input voltage and ablation duration). For a spherical tumour, 30 mm in diameter, complete ablation could not be attained in all configurations with 70 V input voltage and 5 minutes ablation duration. Increasing the input voltage to 90 V enlarged the coagulation zone in the X-model only. With the other configurations, extending the ablation duration to 10 minutes was found to be the better at enlarging the coagulation zone.

Rapid sperm capture: high-throughput flagellar waveform analysis

STUDY QUESTION

Can flagellar analyses be scaled up to provide automated tracking of motile sperm, and does knowledge of the flagellar waveform provide new insight not provided by routine head tracking?

SUMMARY ANSWER

High-throughput flagellar waveform tracking and analysis enable measurement of experimentally intractable quantities such as energy dissipation, disturbance of the surrounding medium and viscous stresses, which are not possible by tracking the sperm head alone.

WHAT IS KNOWN ALREADY

The clinical gold standard for sperm motility analysis comprises a manual analysis by a trained professional, with existing automated sperm diagnostics [computer-aided sperm analysis (CASA)] relying on tracking the sperm head and extrapolating measures. It is not currently possible with either of these approaches to track the sperm flagellar waveform for large numbers of cells in order to unlock the potential wealth of information enclosed within.

STUDY DESIGN, SIZE, DURATION

The software tool in this manuscript has been developed to enable high-throughput, repeatable, accurate and verifiable analysis of the sperm flagellar beat.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Using the software tool [Flagellar Analysis and Sperm Tracking (FAST)] described in this manuscript, we have analysed 176 experimental microscopy videos and have tracked the head and flagellum of 205 progressive cells in diluted semen (DSM), 119 progressive cells in a high-viscosity medium (HVM) and 42 stuck cells in a low-viscosity medium. Unscreened donors were recruited at Birmingham Women's and Children's NHS Foundation Trust after giving informed consent.

MAIN RESULTS AND THE ROLE OF CHANCE

We describe fully automated tracking and analysis of flagellar movement for large cell numbers. The analysis is demonstrated on freely motile cells in low- and high-viscosity fluids and validated on published data of tethered cells undergoing pharmacological hyperactivation. Direct analysis of the flagellar beat reveals that the CASA measure 'beat cross frequency' does not measure beat frequency; attempting to fit a straight line between the two measures gives ${\mathrm{R}}^2$ values of 0.042 and 0.00054 for cells in DSM and HVM, respectively. A new measurement, track centroid speed, is validated as an accurate differentiator of progressive motility. Coupled with fluid mechanics codes, waveform data enable extraction of experimentally intractable quantities such as energy dissipation, disturbance of the surrounding medium and viscous stresses. We provide a powerful and accessible research tool, enabling connection of the mechanical activity of the sperm to its motility and effect on its environment.

LARGE SCALE DATA

The FAST software package and all documentation can be downloaded from www.flagellarCapture.com.

LIMITATIONS, REASONS FOR CAUTION

The FAST software package has only been tested for use with negative phase contrast microscopy. Other imaging modalities, with bright cells on a dark background, have not been tested but may work. FAST is not designed to analyse raw semen; it is specifically for precise analysis of flagellar kinematics, as that is the promising area for computer use. Flagellar capture will always require that cells are at a dilution where their paths do not frequently cross.

WIDER IMPLICATIONS OF THE FINDINGS

Combining tracked flagella with mathematical modelling has the potential to reveal new mechanistic insight. By providing the capability as a free-to-use software package, we hope that this ability to accurately quantify the flagellar waveform in large populations of motile cells will enable an abundant array of diagnostic, toxicological and therapeutic possibilities, as well as creating new opportunities for assessing and treating male subfertility.

STUDY FUNDING/COMPETING INTEREST(S)

M.T.G., G.C., J.C.K-B. and D.J.S. gratefully acknowledge funding from the Engineering and Physical Sciences Research Council, Healthcare Technologies Challenge Award (Rapid Sperm Capture EP/N021096/1). J.C.K-B. is funded by a National Institute of Health Research (NIHR) and Health Education England, Senior Clinical Lectureship Grant: The role of the human sperm in healthy live birth (NIHRDH-HCS SCL-2014-05-001). This article presents independent research funded in part by the NIHR and Health Education England. The views expressed are those of the authors and not necessarily those of the NHS, the NIHR or the Department of Health. The data for experimental set (2) were funded through a Wellcome Trust-University of Birmingham Value in People Fellowship Bridging Award (E.H.O.).The authors declare no competing interests.

Cross validated serum small extracellular vesicle microRNAs for the detection of oropharyngeal squamous cell carcinoma

BACKGROUND

Oropharyngeal squamous cell carcinoma (OPSCC) is often diagnosed at an advanced stage because the disease often causes minimal symptoms other than metastasis to neck lymph nodes. Better tools are required to assist with the early detection of OPSCC. MicroRNAs (miRNAs, miRs) are potential biomarkers for early head and neck squamous cell cancer diagnosis, prognosis, recurrence, and presence of metastatic disease. However, there is no widespread agreement on a panel of miRNAs with clinically meaningful utility for head and neck squamous cell cancers. This could be due to variations in the collection, storage, pre-processing, and isolation of RNA, but several reports have indicated that the selection and reproducibility of biomarkers has been widely affected by the methods used for data analysis. The primary analysis issues appear to be model overfitting and the incorrect application of statistical techniques. The purpose of this study was to develop a robust statistical approach to identify a miRNA signature that can distinguish controls and patients with inflammatory disease from patients with human papilloma virus positive (HPV +) OPSCC.

METHODS

Small extracellular vesicles were harvested from the serum of 20 control patients, 20 patients with gastroesophageal reflux disease (GORD), and 40 patients with locally advanced HPV + OPSCC. MicroRNAs were purified, and expression profiled on OpenArray™. A novel cross validation method, using lasso regression, was developed to stabilise selection of miRNAs for inclusion in a prediction model. The method, named StaVarSel (for Stable Variable Selection), was used to derive a diagnostic biomarker signature.

RESULTS

A standard cross validation approach was unable to produce a biomarker signature with good cross validated predictive capacity. In contrast, StaVarSel produced a regression model containing 11 miRNA ratios with potential clinical utility. Sample permutations indicated that the estimated cross validated prediction accuracy of the 11-miR-ratio model was not due to chance alone.

CONCLUSIONS

We developed a novel method, StaVarSel, that was able to identify a panel of miRNAs, present in small extracellular vesicles derived from blood serum, that robustly cross validated as a biomarker for the detection of HPV + OPSCC. This approach could be used to derive diagnostic biomarkers of other head and neck cancers.

Influence of natural convection on gold nanorods-assisted photothermal treatment of bladder cancer in mice

Background: The thermally-induced urine flow can generate cooling that may alter the treatment outcome during hyperthermic treatments of bladder cancer. This paper investigates the effects of natural convection inside the bladder and at skin surface during gold nanorods (GNR) - assisted photothermal therapy (PTT) of bladder cancer in mice. Methods: 3D models of mouse bladder at orientations corresponding to the mouse positioned on its back, its side and its abdomen were examined. Numerical simulations were carried out for GNR volume fractions of 0.001, 0.005 and 0.01% and laser power of 0.2 and 0.3 W. Results: The obtained results showed that cooling due to natural convection inside the bladder and above the skin depends on the mouse orientation. For a mouse positioned on its back, on its side or on its abdomen, the maximum temperature achieved inside the tumour at 0.001% GNR volume fraction and 0.2 W laser power was 55.2°C, 50.0°C and 52.2°C, respectively compared to 56.8°C when natural convection was not considered. The average thermal gradients when natural convection was considered were also lower, suggesting a more homogenous temperature distribution. Conclusions: Natural convection inside the bladder can be beneficial but also detrimental to GNR-assisted PTT depending on the level of heating. At low levels of heating due to low GNR volume fraction and/or laser power, flow inside the bladder may dissipate heat from the targeted tissue; making the treatment ineffective. At high levels of heating due to high GNR volume fraction and/or laser power, cooling may prevent excessive thermal damage to surrounding tissues.

A numerical study on the no-touch bipolar radiofrequency ablation

The commonly used radiofrequency ablation (RFA) technique for treating liver cancer is in the monopolar mode. This requires the insertion of the RF electrode directly into the tumor tissue, which increases the risks of tumor track seeding (TTS). One way to overcome TTS is by employing the bipolar RFA, implemented in the no-touch mode. In the no-touch mode, two RF electrodes are inserted into the healthy tissue that surrounds the tumor. The distance between the electrodes and the tumor is defined as the no-touch gap. The ability of the no-touch bipolar RFA to overcome TTS has been demonstrated in laboratory studies; however, little is known about the thermo-physiological responses of the tissue during the ablation process of the no-touch procedure. This will be investigated numerically in the present study. A 3D model of the liver tissue is developed and the no-touch bipolar RFA implemented using a pair of RF electrodes is simulated using the finite element method. In particular, the effects of the no-touch gap on the treatment outcome of the RFA procedure are investigated. Results show that a larger no-touch gap may result incomplete tumor destruction due to the central region of the tumor not being directly affected by the Joule heating phenomenon that is more prominent around the electrodes. This suggests that an improperly selected no-touch gap may result in a reduced efficiency of the no-touch bipolar RFA.

Shape-shifting thermal coagulation zone during saline-infused radiofrequency ablation: A computational study on the effects of different infusion location

BACKGROUND AND OBJECTIVE

The majority of the studies on radiofrequency ablation (RFA) have focused on enlarging the size of the coagulation zone. An aspect that is crucial but often overlooked is the shape of the coagulation zone. The shape is crucial because the majority of tumours are irregularly-shaped. In this paper, the ability to manipulate the shape of the coagulation zone following saline-infused RFA by altering the location of saline infusion is explored.

METHODS

A 3D model of the liver tissue was developed. Saline infusion was described using the dual porosity model, while RFA was described using the electrostatic and bioheat transfer equations. Three infusion locations were investigated, namely at the proximal end, the middle and the distal end of the electrode. Investigations were carried out numerically using the finite element method.

RESULTS

Results indicated that greater thermal coagulation was found in the region of tissue occupied by the saline bolus. Infusion at the middle of the electrode led to the largest coagulation volume followed by infusion at the proximal and distal ends. It was also found that the ability to delay roll-off, as commonly associated with saline-infused RFA, was true only for the case when infusion is carried out at the middle. When infused at the proximal and distal ends, the occurrence of roll-off was advanced. This may be due to the rapid and more intense heating experienced by the tissue when infusion is carried out at the electrode ends where Joule heating is dominant.

CONCLUSION

Altering the location of saline infusion can influence the shape of the coagulation zone following saline-infused RFA. The ability to 'shift' the coagulation zone to a desired location opens up great opportunities for the development of more precise saline-infused RFA treatment that targets specific regions within the tissue.

A computational model to investigate the influence of electrode lengths on the single probe bipolar radiofrequency ablation of the liver

BACKGROUND AND OBJECTIVES

Recently, there have been calls for RFA to be implemented in the bipolar mode for cancer treatment due to the benefits it offers over the monopolar mode. These include the ability to prevent skin burns at the grounding pad and to avoid tumour track seeding. The usage of bipolar RFA in clinical practice remains uncommon however, as not many research studies have been carried out on bipolar RFA. As such, there is still uncertainty in understanding the effects of the different RF probe configurations on the treatment outcome of RFA. This paper demonstrates that the electrode lengths have a strong influence on the mechanics of bipolar RFA. The information obtained here may lead to further optimization of the system for subsequent uses in the hospitals.

METHODS

A 2D model in the axisymmetric coordinates was developed to simulate the electro-thermophysiological responses of the tissue during a single probe bipolar RFA. Two different probe configurations were considered, namely the configuration where the active electrode is longer than the ground and the configuration where the ground electrode is longer than the active. The mathematical model was first verified with an existing experimental study found in the literature.

RESULTS

Results from the simulations showed that heating is confined only to the region around the shorter electrode, regardless of whether the shorter electrode is the active or the ground. Consequently, thermal coagulation also occurs in the region surrounding the shorter electrode. This opened up the possibility for a better customized treatment through the development of RF probes with adjustable electrode lengths.

CONCLUSIONS

The electrode length was found to play a significant role on the outcome of single probe bipolar RFA. In particular, the length of the shorter electrode becomes the limiting factor that influences the mechanics of single probe bipolar RFA. Results from this study can be used to further develop and optimize bipolar RFA as an effective and reliable cancer treatment technique.

Pathophysiology of swallowing following oropharyngeal surgery for obstructive sleep apnea syndrome

BACKGROUND

Uvulopalatopharyngoplasty (UPPP) and coblation channeling of the tongue (CCT) are oropharyngeal surgeries used to treat obstructive sleep apnea syndrome. The extent to which UPPP and CCT affect pharyngeal swallow has not been determined. We therefore conducted a novel case series study employing high-resolution impedance manometry (HRIM) to quantify the swallowing-related biomechanics following UPPP and/or CCT surgery.

METHODS

Twelve patients who underwent UPPP+CCT or CCT only were assessed an average 2.5 years postsurgery. Swallow function data were compared with ten healthy controls. All patients completed the Sydney swallow questionnaire (SSQ). Pharyngeal pressure-flow analysis of HRIM recordings captured key distension, contractility and pressure-flow timing swallow parameters testing 5, 10, and 20 mL volumes of thin and thick fluid consistencies.

KEY RESULTS

Postoperative patients had more dysphagia symptoms with five returning abnormal SSQ scores. Swallowing was biomechanically altered compared to controls, consistent with diminished swallowing reserve, largely driven by elevated hypopharyngeal intrabolus pressure due to a reduced capacity to open the upper esophageal sphincter to accommodate larger volumes.

CONCLUSIONS & INFERENCES

Patients who have undergone UPPP and/or CCT surgery appear to have a deficiency in normal modulation of the swallowing mechanism and a reduced swallowing functional reserve. We speculate that these changes may become relevant in later life with the onset of age-related stressors to the swallowing mechanism. This case series strikes a note of caution that further studies are needed to determine the role of preoperative swallow assessment in patients undergoing UPPP and/or CCT surgery.

A systematic review of validated sinus surgery simulators

BACKGROUND

Simulation provides a safe and effective opportunity to develop surgical skills. A variety of endoscopic sinus surgery (ESS) simulators has been described in the literature. Validation of these simulators allows for effective utilisation in training.

OBJECTIVE OF REVIEW

To conduct a systematic review of the published literature to analyse the evidence for validated ESS simulation.

SEARCH STRATEGY

Pubmed, Embase, Cochrane and Cinahl were searched from inception of the databases to 11 January 2017.

EVALUATION METHOD

Twelve thousand five hundred and sixteen articles were retrieved of which 10 112 were screened following the removal of duplicates. Thirty-eight full-text articles were reviewed after meeting search criteria. Evidence of face, content, construct, discriminant and predictive validity was extracted.

RESULTS

Twenty articles were included in the analysis describing 12 ESS simulators. Eleven of these simulators had undergone validation: 3 virtual reality, 7 physical bench models and 1 cadaveric simulator. Seven of the simulators were shown to have face validity, 7 had construct validity and 1 had predictive validity. None of the simulators demonstrated discriminate validity.

CONCLUSION

This systematic review demonstrates that a number of ESS simulators have been comprehensively validated. Many of the validation processes, however, lack standardisation in outcome reporting, thus limiting a meta-analysis comparison between simulators.

Validation of virtual-reality-based simulations for endoscopic sinus surgery

BACKGROUND

Virtual reality (VR) simulators provide an alternative to real patients for practicing surgical skills but require validation to ensure accuracy. Here, we validate the use of a virtual reality sinus surgery simulator with haptic feedback for training in Otorhinolaryngology - Head & Neck Surgery (OHNS).

METHODS

Participants were recruited from final-year medical students, interns, resident medical officers (RMOs), OHNS registrars and consultants. All participants completed an online questionnaire after performing four separate simulation tasks. These were then used to assess face, content and construct validity. anova with post hoc correlation was used for statistical analysis.

RESULTS

The following groups were compared: (i) medical students/interns, (ii) RMOs, (iii) registrars and (iv) consultants. Face validity results had a statistically significant (P < 0.05) difference between the consultant group and others, while there was no significant difference between medical student/intern and RMOs. Variability within groups was not significant. Content validity results based on consultant scoring and comments indicated that the simulations need further development in several areas to be effective for registrar-level teaching. However, students, interns and RMOs indicated that the simulations provide a useful tool for learning OHNS-related anatomy and as an introduction to ENT-specific procedures.

CONCLUSIONS

The VR simulations have been validated for teaching sinus anatomy and nasendoscopy to medical students, interns and RMOs. However, they require further development before they can be regarded as a valid tool for more advanced surgical training.

Effects of aqueous humor hydrodynamics on human eye heat transfer under external heat sources

The majority of the eye models developed in the late 90s and early 00s considers only heat conduction inside the eye. This assumption is not entirely correct, since the anterior and posterior chambers are filled aqueous humor (AH) that is constantly in motion due to thermally-induced buoyancy. In this paper, a three-dimensional model of the human eye is developed to investigate the effects AH hydrodynamics have on the human eye temperature under exposure to external heat sources. If the effects of AH flow are negligible, then future models can be developed without taking them into account, thus simplifying the modeling process. Two types of external thermal loads are considered; volumetric and surface irradiation. Results showed that heat convection due to AH flow contributes to nearly 95% of the total heat flow inside the anterior chamber. Moreover, the circulation inside the anterior chamber can cause an upward shift of the location of hotspot. This can have significant consequences to our understanding of heat-induced cataractogenesis.

Cell death, perfusion and electrical parameters are critical in models of hepatic radiofrequency ablation

Purpose: A sensitivity analysis has been performed on a mathematical model of radiofrequency ablation (RFA) in the liver. The purpose of this is to identify the most important parameters in the model, defined as those that produce the largest changes in the prediction. This is important in understanding the role of uncertainty and when comparing the model predictions to experimental data. Materials and methods: The Morris method was chosen to perform the sensitivity analysis because it is ideal for models with many parameters or that take a significant length of time to obtain solutions. A comprehensive literature review was performed to obtain ranges over which the model parameters are expected to vary, crucial input information. Results: The most important parameters in predicting the ablation zone size in our model of RFA are those representing the blood perfusion, electrical conductivity and the cell death model. The size of the 50 °C isotherm is sensitive to the electrical properties of tissue while the heat source is active, and to the thermal parameters during cooling. Conclusions: The parameter ranges chosen for the sensitivity analysis are believed to represent all that is currently known about their values in combination. The Morris method is able to compute global parameter sensitivities taking into account the interaction of all parameters, something that has not been done before. Research is needed to better understand the uncertainties in the cell death, electrical conductivity and perfusion models, but the other parameters are only of second order, providing a significant simplification.

Lysozyme expression is increased in the sinus mucosa of patients with chronic rhinosinusitis

BACKGROUND

The presence of fungi and bacteria in the paranasal sinuses may contribute to ongoing inflammation. Lysozyme is an innate immune peptide with bactericidal and fungicidal activity. The expression of lysozyme in chronic rhinosinusitis (CRS) is poorly understood and deficiencies in lysozyme expression may contribute to the ongoing inflammation in CRS patients.

OBJECTIVE

Determine lysozyme expression in sinus mucosa of normal and CRS patients with (CRSwNP) and without (CRSsNP) nasal polyps.

METHODOLOGY

Sinus mucosa specimens (n = 82) were processed for standard histology, immunohistochemical localisation of lysozyme, immunofluorescent localisation of fungi, and qPCR analysis of lysozyme expression.

RESULTS

CRS specimens displayed high-levels of lysozyme immunoreactivity in many of the abundant serous cells. Moderate levels were detected in some epithelial cells and inflammatory cells. Low levels were detected in some subepithelial glands of control specimens. No difference in immunoreactivity was detected between CRSwNP and CRSsNP specimens. Fungal elements were not visualised in any sinus specimen. qPCR analysis demonstrated variable lysozyme expression between individuals.

CONCLUSIONS

Lysozyme protein expression is increased in patients with CRS, suggesting a defect in lysozyme expression is not responsible for the microbial colonisation often associated with CRS. The functional activity of lysozyme in CRS patients needs to be further investigated.

Prediction and parametric analysis of thermal profiles within heated human skin using the boundary element method

In this paper, an axisymmetric model of the human skin is developed to simulate the steady-state temperature distribution during contact with a hot solid. Simulations are carried out using the boundary element method. This study seeks to investigate the feasibility of using the boundary element method in the studies of burn. A sensitivity analysis is carried out to examine the effects of various parameters on the temperature distribution inside the skin during burn. Furthermore, a statistical analysis based on the Taguchi method is performed to determine the combination of factors that produce the desired outcome (least increase in temperature). In order to validate the accuracy of the numerical scheme, results obtained using the boundary element method are compared with the solutions obtained using the more established finite-element method.

Boundary element method with bioheat equation for skin burn injury

Burns are second to vehicle crashes as the leading cause of non-intentional injury deaths in the United States. The survival of a burn patient actually depends on the seriousness of the burn. It is important to understand the physiology of burns for a successful treatment of a burn patient. This has prompted researchers to conduct investigations both numerically and experimentally to understand the thermal behaviour of the human skin when subjected to heat injury. In this study, a model of the human skin is developed where the steady state temperature during burns is simulated using the boundary element method (BEM). The BEM is used since it requires boundary only discretion and thus, reduces the requirement of high computer memory. The skin is modeled as three layered in axisymmetric coordinates. The three layers are the epidermis (uppermost), dermis (middle) and subcutaneous fat. Burning is applied via a heating disk which is assumed to be at constant temperature. The results predicted by the BEM model showed very good agreement with the results obtained using the finite element method (FEM). The good agreement despite using only linear elements as compared to quadratic elements in the FEM model shows the versatility of the BEM. A sensitivity analysis was conducted to investigate how changes in the values of certain skin variables such as the thermal conductivity and environmental conditions like the ambient convection coefficient affect the temperature distribution inside the skin. The Taguchi method was also applied to identify the combination of parameters which produces the largest increase in skin temperature during burns.

Variations in the corneal surface temperature with contact lens wear

This paper presents the two-dimensional simulation of heat propagation in the human eye model during contact lens wear with finite element analysis. Three types of contact lens are studied: Lotrafilcon A, Balafilcon A, and Etafilcon A. The models are solved for both steady and transient solutions. The corneal surface temperature during contact lens wear is found to decrease (average, 0.52 ± 0.05°C compared with a bare cornea for all lens types). A contact lens with a higher water content has a lower steady state temperature than a contact lens with a lower water content does. Various initial temperatures for the contact lens are found to affect the first 400 s of the temperature variation. When the initial temperature is lower than the corneal temperature, a reduction in temperature is observed during contact lens insertion while the opposite is observed when the initial temperature is higher than the corneal temperature. The increase in evaporation rate when a contact lens is worn increases the cooling effect on the ocular surface. This is suggested to be the cause of lower corneal surface temperature when wearing a contact lens.

Ocular surface temperature: a 3D FEM prediction using bioheat equation

Computational and mathematical human eye models from previous studies which were constructed in two-dimensions (2D) did not give a precise representation of the actual human eye. This work is an extension from an earlier published work on the 2D model. In this paper, a 3D FEM model of the human eye is simulated for the steady state temperature distribution during normal condition and during electromagnetic (EM) wave radiation. Results show a discrepancy of 0.49% for a normal condition as opposed to 1.9% of a 2D model when compared to experimental results from open literatures. Investigations on the EM wave radiations found an average power absorption density of 15,151 and 22,145 Wm(-3) for the 750 and 1500 MHz radiation, respectively. A peak temperature of 38.18( composite function)C was predicted for the 750 MHz radiation while 41.19( composite function)C was computed for the 1500 MHz radiation. These temperatures are in reasonable agreement with the simulated results computed by another report in the past.

FEM simulation of the eye structure with bio-heat analysis

Computer simulation on medical sciences has gain increasing popularity as computational technology advances. Successful thermal modeling of the human eye will assist in enabling early detections of eye abnormalities such as inflammatory. However, validity of every computer simulated results must be benchmarked with experimental measurement and this can be a daunting task especially in biomedical fields where experimental data is not in abundance. This paper presents a 2D finite element (FE) human eye model developed to simulate its thermal steady state conditions based on the properties and parameters reported in the open literatures. The results are verified with experimental and computational results obtained by previous studies on human as well as animal eyes. Results show discrepancy of only 0.33% when compared to images from infrared (IR) screening and a difference of only 0.127% compared to another finite element model. The sensitivity analysis also provides good agreement with results by previous studies. This promising simulation allows new possibility in computational methods for eye health care.

Persistent buccopharyngeal membrane in an adult

Persistence of the buccopharyngeal membrane (BPM) also called the oropharyngeal membrane is a rare congenital oropharyngeal anomaly. We report a case of an adult aboriginal male patient with a membrane that closed his oropharyngeal isthmus except for a 2 cm diameter central perforation. The patient had no symptoms related to this membrane and no other congenital anomalies were found. This finding has not previously been reported in an adult. The embryology and management of this rare condition is discussed.