Photoacoustic-guided focused ultrasound for accurate visualization of brachytherapy seeds with the photoacoustic needle

An important problem in minimally invasive photoacoustic (PA) imaging of brachytherapy seeds is reflection artifacts caused by the high signal from the optical fiber/needle tip reflecting off the seed. The presence of these artifacts confounds interpretation of images. In this letter, we demonstrate a recently developed concept called photoacoustic-guided focused ultrasound (PAFUSion) for the first time in the context of interstitial illumination PA imaging to identify and remove reflection artifacts. In this method, ultrasound (US) from the transducer is focused on the region of the optical fiber/needle tip identified in a first step using PA imaging. The image developed from the US diverging from the focus zone at the tip region visualizes only the reflections from seeds and other acoustic inhomogeneities, allowing identification of the reflection artifacts of the first step. These artifacts can then be removed from the PA image. Using PAFUSion, we demonstrate reduction of reflection artifacts and thereby improved interstitial PA visualization of brachytherapy seeds in phantom and ex vivo measurements on porcine tissue. © The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI. [DOI: 10.1117/1.JBO.21.12.120501]


Photoacoustic-guided focused ultrasound for accurate visualization of brachytherapy seeds with the photoacoustic needle
Prostate cancer is one of the most common forms of cancer among men. 1 Depending on the stage and severity of prostate cancer, several therapeutic options such as radiation therapy, hormone therapy, and radical prostatectomy are generally used.Brachytherapy is a targeted radiation therapy treatment which enables high-precision delivery of radiation doses to the tumor area in the prostate.In brachytherapy, the localized tumor area instead of the whole prostate is targeted with defined radiation dose, using multiple (50 to 100) tiny metallic radioactive sources implanted inside the prostate.In brachytherapy, complications may arise when there is a deviation between the planned and delivered radiation dose. 2 Potential reasons for this deviation are seed migration, prostate motion, edema, or implantation errors. 2 Mismatch between planned and delivered radiation dose can be avoided to a large extent by accurate visualization of seeds during the implantation procedure.
Real-time seed visualization during the implantation procedure is currently performed using transrectal ultrasound imaging.However, seeds are often challenging to visualize using ultrasound (US) due to their small size, their similarity in acoustic impedance with microcalcifications, and due to the presence of characteristic acoustic artifacts such as shadowing and reverberations. 3n photoacoustic (PA) imaging, pulsed light absorbed by the target emits US by the thermoelastic mechanism. 4This US is detected using a US detector array to build up an image.PA imaging offers spectroscopic optical absorption contrast with US resolution. 5Owing to the high optical absorption of metal compared to tissue, PA imaging offers excellent contrast in imaging brachytherapy seeds. 3,6Thus, PA imaging can complement US imaging for better visualization of brachytherapy seeds. 3,6,7onsidering the imaging depth required, the PA needle 8,9 with interstitial light delivery via a cutting percutaneous needle is generating much interest in imaging brachytherapy seeds. 7A challenge in interstitial PA imaging of brachytherapy seeds is the presence of reflection artifacts caused by the high PA signal from the optical fiber/needle tip reflecting off the seeds as shown in Fig. 1(b).These artifacts may appear in the region of interest and can reduce the contrast and causes difficulty in image interpretation as reflection artifacts can be misinterpreted as seeds. 7It is important to identify and reduce reflection artifacts for accurate interstitial PA visualization of brachytherapy seeds.
We recently proposed a method called photoacoustic-guided focused ultrasound (PAFUSion) to identify reflection artifacts in epi-mode PA imaging. 10In this work, we study the use of PAFUSion in identifying and reducing in-plane reflection artifacts with the use of the PA needle in the visualization of brachytherapy seeds.We present the results of the combined approach in visualizing nonradioactive brachytherapy seeds in a tissue mimicking phantom and an ex vivo porcine tissue.In addition, we also show the potential of using multispectral PA for separating seed signals from other optical absorbers like blood in a phantom measurement.
The system consists of a commercial handheld US imager (MyLab_One, Esaote Europe BV, Maastricht, The Netherlands) with a linear array probe SL3323 (128 elements, 7.5 MHz, 100% bandwidth). 11The system is used in the research mode with 64 elements where US transmission, laser pulse transmission, and data acquisition are controlled using custom-made software running on a PC.Illumination uses an Nd:YAG laser (Quanta-Ray Pro 250, Spectra Physics, Mountain View, California) pumping an optical parametric oscillator (VersaScan-L532, GWU, Erftstadt, Germany) at 10 Hz.The energy per pulse is monitored using a photodetector.A multimode optical fiber (NA 0.22, diameter 600 μm) coupled to the laser delivers 10 ns laser light with 1.2 mJ∕pulse at 800 nm.For the multiwavelength studies, the range of 720 to 860 nm was covered with 5 nm steps; single wavelength studies used 800 nm.The fiber is held in a biopsy needle (14 G × 13 cm, Bard Ltd., UK) and slide-fits into the 14 G cutting cannula [Fig.1(a)].The PA frame rate is 10 Hz and images can be displayed online using a GPU-based delay and sum reconstruction algorithm.Images presented here are averaged over 100 frames and reconstructed offline using a frequency domain reconstruction algorithm, 12 assuming a constant speed of sound in the medium for all measurements (phantom measurement: 1480 m∕s, ex vivo measurement: 1540 m∕s).
The PAFUSion approach specifically in the context of the PA needle is shown in Figs.shows the schematic of the experiment to validate the proposed combined approach in a scattering phantom.The phantom (10 cm × 7.5 cm × 5 cm) was prepared with Agar dissolved in an aqueous solution of Intralipid (I141, 20% emulsion, Sigma-Aldrich) to provide a bulk reduced scattering coefficient of μ 0 s ¼ 1 mm −1 .Within the phantom, three nonradioactive brachytherapy seeds (length ¼ 4.5 mm × 0.8 mm, Bard Benelux N.V., The Netherlands) were embedded.The seeds were coated with a single layer of India ink (∼200 μm) by dipping them in the ink and drying in air.India ink coating was used in this proof-of-principle study to improve the contrast of the brachytherapy seeds by around 10 dB in the PA image.For clinical implementation, a biocompatible coating with high optical absorption could be more suitable.A single bead (2.2-mm diameter) of sodium alginate carrying gold nanorods (AuNR) 13 was also embedded in the phantom.The reference absorption spectrum of the bead [Fig.2(d)], measured using a spectrophotometer (UV-2401PC, Shimadzu, Kyoto, Japan), shows a characteristic longitudinal plasmon absorption peak at 795 nm.Water coupling was used in between the US probe and the phantom.The needle was oriented parallel to the US detector array.PA and plane-wave US images are acquired in the first step.Then after locating the fiber/needle tip, PAFUSion (focused US) data was also collected for offline processing.PA images for different excitation wavelengths (720 to 860 nm) were also acquired.
Ex vivo validation of the approach was performed using four ink-coated brachytherapy seeds embedded inside porcine tissue.Using US guidance, the seeds were implanted inside the tissue using a standard 18G brachytherapy needle (Bard, Brachystar, Covington, Georgia) by one of the authors (E.M.).The same data acquisition and processing protocol were followed as for the phantom experiment.For both the phantom and ex vivo experiments, the seeds were aligned inside the imaging plane of the US transducer array.Orientation of the needle and arrangement of seeds inside the porcine tissue are schematically shown in Fig. 3(a).For better acoustic coupling, the sample was covered with transparent kitchen cling film (∼20 μm) after smearing with US gel [Fig.3(d)].PAFUSion image uncovered this feature as a reflection artifact.This would have been caused by the fiber/needle tip signal reflecting off the point-like tissue feature visible just above seed 1 (lateral distance ¼ 7.5 mm, depth ¼ 10 mm) in both the PA and US images.
Our approach has been successfully demonstrated in phantom and in ex vivo settings.The method has been developed in the context of PA needle imaging of brachytherapy seeds, but is eminently suited for all applications of the PA needle where reflection artifacts could result due to the presence of acoustic reflectors.However, the approach needs to be further investigated.One issue is that high absorption of heterogeneous tissue proximal to the needle tip could result in multiple bright PA features near the tip.This may create uncertainties in accurately identifying the fiber/needle tip which is essential for using the approach.Another issue that requires investigation is the approach used for correction of the PA images.Subtraction of the PAFUSion images from the PA image worked well in the situations studied, but will require optimization in the case of real PA features mixed with any reflection artifacts.
In addition, PAFUSion will be able to identify and reduce reflection artifacts generated inside the imaging plane only.If the fiber/needle tip is outside the imaging plane, it will be impossible to reduce reflection artifacts using our approach.However, in-plane needle insertion is commonly used in clinical practice for accurate guidance of the needle tip, making this less a limitation.
In linear US array-based PA imaging, one can expect significant contributions from absorbing objects outside the imaging plane as direct clutter. 14These out-of-plane signals may even reflect off echogenic structures inside the imaging plane to generate indirect echo clutter. 14However, these out-of-plane PA signals will be significant only when the tissue is illuminated from the outside in a noninvasive setting unlike the present interstitial illumination setting.
Angular distribution of ultrasonic wavefield after the focus depth is an important parameter when using PAFUSion in the context of the PA needle, as this angle defines the fieldof-view of the method in identifying and eliminating reflection artifacts.In our case, where the center frequency was 7.5 MHz and the pitch was 0.245 mm, this angle limit was at AE24 deg, and was sufficient to identify the reflection artifacts in this study.
In both the phantom and ex vivo experiments, light energy was 1.2 mJ∕pulse, resulting in a fluence of 425 mJ∕cm 2 at the fiber tip.The maximum permissible exposure (MPE) for NIR light inside tissue is not defined (for skin the MPE is 20 mJ∕cm 2 for the wavelength used).However, the high exposure at the tip raises safety concerns from a clinical standpoint.We envisage addressing this issue in future using fiber tip modifications, e.g., using a diffusing tip. 15oating the needle tip with a dye/paint with a spectral signature and identifying the reflection artifacts caused by it using spectral decomposition and correcting for it will potentially be an alternate solution for the reflection artifact reduction in interstitial PA imaging.However, this requires further investigation.
We conclude that PAFUSion enables identification and substantial reduction of reflection artifacts when using the PA needle, considerably improving the visualization of brachytherapy seeds.Our phantom results also confirm that multispectral PA 6 can be potentially a useful tool to separate signals from the brachytherapy seeds and other optical absorbers in the tissue.The PAFUSion approach reported here is suitable for real-time implementation since it requires only a single focus US measurement and an image subtraction step.In a future clinical scenario, we envisage that the PA needle in combination with realtime PAFUSion will provide accurate visualization of the distribution of brachytherapy seeds, which in turn provide feedback to the radiotherapist regarding adherence to the therapy plan.

Disclosures
Authors declare that there is no conflict of interest regarding the publication of this paper.
1(b) and 1(c).Once the needle tip has been identified from the PA image (by manually identifying the position with the highest pixel intensity at the tip region in the image) in step 1 [Fig.1(b)], US is focused to this point in a second step.The resultant US data processed using the PAFUSion algorithm 10 generates an image with only reflections from deeper structures (PAFUSion image) hence representing the artifacts in the PA image from step 1.The average image pixel value of the reflection artifact with the highest image value is then scaled to the image value of its counterpart in the PA image, and the same scaling is applied to all other identified artifacts in the PAFUSion image.After this, reconstructed and enveloped PAFUSion data are subtracted from the PA data and the positive parts of the result are plotted to obtain the corrected PA image.The negative parts of the result are assigned an image value of zero.

Fig. 1
Fig. 1 (a) Block diagram of setup.RM, reflecting mirror; BS, glass plate to sample beam; PD, photodetector; NDF, neutral density filter; L, lens; OF, optical fiber; and UST, ultrasound transducer.(b) Concept of the PA needle, where pulsed light is applied in tissue via an optical fiber in a cutting needle with external US detection.By this, the tissue at the optical fiber tip, the tip itself, and the brachytherapy seeds can be visualized.Additionally, however, artifacts are also visualized due to reflection of the high PA signal from the fiber/needle tip on the seeds.(c) Concept of PAFUSion where in a second step; US from the transducer is focused at the optical fiber/needle tip.The US diverging from the focus zone at the tip is in a first approximation a mimic of the PA signal from tip in the first step.The detected signals visualize only the reflections from seed and other inhomogeneities allowing identification of the reflection artifacts of the first step.

Fig. 2
Fig. 2 (a) Schematic of the experiment showing orientation of the PA needle, arrangement of seeds and AuNR bead in the Agar phantom with Intralipid, (b) PA image at a wavelength of 800 nm showing the expected features and with an additional feature at 9.5 × 24 mm 2 (lateral distance × depth), (c) planewave US image, (d) measured PA spectrum of seed 1, bead with error margins (standard deviation of the averaged mean) compared with reference spectrum of the bead, (e) PAFUSion image showing that the unknown feature is a reflection artifact, and (f) corrected PA image obtained by subtracting scaled PAFUSion image from the PA image.

Figure 2 (
Figure2(a) shows the schematic of the experiment to validate the proposed combined approach in a scattering phantom.The phantom (10 cm × 7.5 cm × 5 cm) was prepared with Agar dissolved in an aqueous solution of Intralipid (I141, 20% emulsion, Sigma-Aldrich) to provide a bulk reduced scattering coefficient of μ 0 s ¼ 1 mm −1 .Within the phantom, three nonradioactive brachytherapy seeds (length ¼ 4.5 mm × 0.8 mm, Bard Benelux N.V., The Netherlands) were embedded.The seeds were coated with a single layer of India ink (∼200 μm) by dipping them in the ink and drying in air.India ink coating was used in this proof-of-principle study to improve the contrast of the brachytherapy seeds by around 10 dB in the PA image.For clinical implementation, a biocompatible coating with high optical absorption could be more suitable.A single bead (2.2-mm diameter) of sodium alginate carrying gold nanorods (AuNR)13 was also embedded in the phantom.The reference absorption spectrum of the bead [Fig.2(d)],measured using a spectrophotometer (UV-2401PC, Shimadzu, Kyoto, Japan), shows a characteristic longitudinal plasmon absorption peak at 795 nm.Water coupling was used in between the US probe and the phantom.The needle was oriented parallel to the US detector array.PA and plane-wave US images are acquired in the first step.Then after locating the fiber/needle tip, PAFUSion (focused US) data was also collected for offline processing.PA images for different excitation wavelengths (720 to 860 nm) were also acquired.Ex vivo validation of the approach was performed using four ink-coated brachytherapy seeds embedded inside porcine tissue.Using US guidance, the seeds were implanted inside the tissue using a standard 18G brachytherapy needle (Bard, Brachystar, Covington, Georgia) by one of the authors (E.M.).The same data acquisition and processing protocol were followed as for the phantom experiment.For both the phantom and ex vivo experiments, the seeds were aligned inside the imaging plane of the US transducer array.Figure 2(b) shows the reconstructed PA image (phantom experiment) in which the shaft and the fiber/needle tip and all three seeds are visualized.The AuNR bead is seen only with top and bottom surfaces, as expected in a limited-view situation.The US image [Fig.2(c)] shows the needle shaft

Figure 2 (
Figure2(a) shows the schematic of the experiment to validate the proposed combined approach in a scattering phantom.The phantom (10 cm × 7.5 cm × 5 cm) was prepared with Agar dissolved in an aqueous solution of Intralipid (I141, 20% emulsion, Sigma-Aldrich) to provide a bulk reduced scattering coefficient of μ 0 s ¼ 1 mm −1 .Within the phantom, three nonradioactive brachytherapy seeds (length ¼ 4.5 mm × 0.8 mm, Bard Benelux N.V., The Netherlands) were embedded.The seeds were coated with a single layer of India ink (∼200 μm) by dipping them in the ink and drying in air.India ink coating was used in this proof-of-principle study to improve the contrast of the brachytherapy seeds by around 10 dB in the PA image.For clinical implementation, a biocompatible coating with high optical absorption could be more suitable.A single bead (2.2-mm diameter) of sodium alginate carrying gold nanorods (AuNR)13 was also embedded in the phantom.The reference absorption spectrum of the bead [Fig.2(d)],measured using a spectrophotometer (UV-2401PC, Shimadzu, Kyoto, Japan), shows a characteristic longitudinal plasmon absorption peak at 795 nm.Water coupling was used in between the US probe and the phantom.The needle was oriented parallel to the US detector array.PA and plane-wave US images are acquired in the first step.Then after locating the fiber/needle tip, PAFUSion (focused US) data was also collected for offline processing.PA images for different excitation wavelengths (720 to 860 nm) were also acquired.Ex vivo validation of the approach was performed using four ink-coated brachytherapy seeds embedded inside porcine tissue.Using US guidance, the seeds were implanted inside the tissue using a standard 18G brachytherapy needle (Bard, Brachystar, Covington, Georgia) by one of the authors (E.M.).The same data acquisition and processing protocol were followed as for the phantom experiment.For both the phantom and ex vivo experiments, the seeds were aligned inside the imaging plane of the US transducer array.Figure 2(b) shows the reconstructed PA image (phantom experiment) in which the shaft and the fiber/needle tip and all three seeds are visualized.The AuNR bead is seen only with top and bottom surfaces, as expected in a limited-view situation.The US image [Fig.2(c)] shows the needle shaft

Figure 3 (
b) shows the reconstructed PA image in which PA needle/ fiber tip, all four seeds, and several other bright features are visible.It is difficult to interpret this image since several bright features with shapes similar to the seeds are present.The shaft of the PA needle, all the four seeds, and several other echogenic tissue structures are visible in the US image [Fig.3(c)].

Figure 3 (
Fig. 3 (a) Schematic showing orientation of the PA needle and arrangement of seeds inside the porcine tissue, (b) PA image at a wavelength of 800 nm, (c) plane-wave US image, (d) photograph of the imaging situation with needle inserted in the ex vivo sample.Inset is a photograph of the inner and outer needles alongside the seed.(e) PAFUSion image with the identified reflection artifacts, and (f) corrected PA image obtained by subtracting the scaled PAFUSion image from the PA image.