Ultrasonic actuation of a fine-needle improves biopsy yield

Despite the ubiquitous use over the past 150 years, the functions of the current medical needle are facilitated only by mechanical shear and cutting by the needle tip, i.e. the lancet. In this study, we demonstrate how nonlinear ultrasonics (NLU) extends the functionality of the medical needle far beyond its present capability. The NLU actions were found to be localized to the proximity of the needle tip, the SonoLancet, but the effects extend to several millimeters from the physical needle boundary. The observed nonlinear phenomena, transient cavitation, fluid streams, translation of micro- and nanoparticles and atomization, were quantitatively characterized. In the fine-needle biopsy application, the SonoLancet contributed to obtaining tissue cores with an increase in tissue yield by 3–6× in different tissue types compared to conventional needle biopsy technique using the same 21G needle. In conclusion, the SonoLancet could be of interest to several other medical applications, including drug or gene delivery, cell modulation, and minimally invasive surgical procedures.


Scientific Reports
| (2021) 11:8234 | https://doi.org/10.1038/s41598-021-87303-x www.nature.com/scientificreports/ extending beyond the needle tip has been overlooked as a way to add value to the function of medical needles. In fact, the NLU can generate acoustic force fields providing a precise manipulation of entities selectively closer or at a farther distance from the sound source. Fine-needle aspiration biopsy (FNAB) is a common biopsy method in which a hypodermic needle and a syringe are employed to aspirate tissue constructs from a target tissue, e.g. suspected tumor. The needle is translated by the operator, while the needle tip is inside the target tissue. Suction by the syringe causes the tissue to bulge towards the needle lumen and adhere to the inner walls of the needle, while the rapid movement of the needle tip will cause the protruded portion of tissue to be sliced off by the cutting edges of the needle and to be aspirated into the lumen 28 . Here the source of the external force is the operator hand, which translates to compressive, shear and tensile forces at the needle tip and explain the detachment of cells and tissue constructs from the target tissue. Contrary to FNAB, we hypothesize that ultrasound could serve as an external source of forces near the needle tip further accelerating the detachment of cells and tissue constructs.
Coupling flexural waves to a needle induces sideways motion of the needle tip. This is anticipated to exert direct shear, compression and tensile stresses to the needle tip environment as well as induce sound emission beyond the forces induced by the operator's hand movement. The sound emission could generate NLU phenomena such as acoustic radiation force 29 , acoustic streaming 30 and cavitation 31 . In the context of biopsy, the direct and indirect forces are expected to contribute to extracting cells and tissue constructs, which could enhance the yield in biopsy.
The aim of this study was to investigate how adding ultrasound to a conventional medical needle could change its conventional functions. We first demonstrated how ultrasound emission could be localized at a tip of a conventional medical needle, as an energy source for NLU around the needle tip. From micrometer to millimeter scale, we quantitatively characterized the NLU effects potentially contributing to tissue actuation, such as cavitation, acoustic radiation force and atomization spatially reaching beyond the physical lancet; we defined the confined volume around the needle tip exhibiting pronounced nonlinear effects as the SonoLancet. Finally, we demonstrated the capability of SonoLancet to increase the sample yield in the FNAB application.

Results
Localization of ultrasound near the needle tip. To demonstrate that the effects potentially contributing to tissue actuation can be generated near an ultrasonically actuated medical needle ( Fig. 1A(1-3)), a 21 gauge hypodermic needle was selected, because it represents a common medical needle. Moreover, needles of this diameter or less (fine needles) induce less tissue trauma than larger needles. The selected geometry permits access to the needle hub with a variety of adjuvant devices commonly used in medical practice, e.g. tubing, syringe, vacutainer; for this study we selected a 10 mL syringe.
We employed an axially translating Langevin transducer as a sound source, which was coupled to an S-shaped waveguide connected to a medical needle ( Fig. 1A(1)). Using the waveguide, the longitudinal mode waves of the ultrasound transducer are converted to flexural mode waves in the needle. Reflection of the wave from the needle tip allowed generation of a flexural standing wave between the tips of the needle and the waveguide. Moreover, the converging structure of the medical needle bevel gives geometric amplification of the wave towards the needle extremity. As a consequence of the wave amplification, the needle tip is made to act as a dipole-like sound source oscillating at a large displacement amplitude (Fig. 1D), exhibiting considerably greater time-averaged acoustic intensity at the needle tip ( Fig. 1B(1,2)) than elsewhere near the needle lumen, which carries the energy. These characteristics have specific relevance to nonlinear ultrasonics associated with time-averaged intensity (e.g. acoustic radiation pressure) or characterized by threshold behavior (e.g. cavitation), since the NLU phenomena can be limited to the proximity of the needle tip.
As the ultrasound field near the needle is largely confined within a distance of few mm near the needle bevel ( Fig. 1B(1)), the localization of the SonoLancet is spatially precise and microtrauma in tissues farther away from the needle tip is expected to be minimized. In fact, the threshold phenomenon of cavitation and the fast geometric attenuation of the acoustic field near the needle exclude explanations where cavitation could unexpectedly occur very far from the needle. Cells and tissue components influenced by the SonoLancet can be withdrawn through the needle lumen directly from the site of ultrasound action or alternatively, cells and agents can be delivered to a spatially well-defined site in conjunction with ultrasonic actions.
Cavitation. Large displacements in the direction of x-axis (Fig. 1D) detected near the lancet point produce pronounced acoustic intensity radiating outwards from the center axis of the needle. The high magnitude of the peak-negative-pressure, associated with the acoustic intensity, elevates the probability of ultrasound-microbubble interactions, i.e. cavitation ( Fig. 1C(1,2)). Cavitation activity exhibited within a region extending less than 2 mm from the tip along the positive z-direction. Optical high-speed (HS) imaging (300,000 fps) revealed that the temporal probability of observing cavitation across a time window of 100 ms can be up to 50% in this region ( Fig. 1C(2)), while in the proximity of the rest of the needle there were no signs of cavitation. This spatiotemporal behaviour arises from the inertial cavitation being a threshold phenomenon 32 . Cavitation, i.e. growth, oscillation and collapse of gas bubbles, is a consequence of the fluctuating pressure amplitude being momentarily the lowest in this area due to the pronounced amplification of needle tip motion compared to other locations along the needle lumen ( Fig. 1B(1,2)).
While the needle tip displacement was observed to be < 100 µm (Fig. 1D), the bubble-water boundary displacement was more pronounced, extending beyond 300 µm along the positive x-direction. This is due to the different compliances of water and an air bubble. The activity of the primary bubble within the SonoLancet produced fluid movements of several m s −1 (Fig. 1G), accelerations up to five orders of magnitude of gravity (Fig. 1H), thus allowing production of high shear stresses localized near the needle tip 33  www.nature.com/scientificreports/ produce secondary sound formation, a potential mechanism of shock wave formation 34 and strong and transient shear forces have great potential to actuate matter, such as micro-and nanoparticles or tissue such as cells, groups of cells or localized regions of an organ.
Acoustic radiation pressure. Acoustic radiation force can be employed to push or pull medium, objects and interfaces 35 . It is a force exerted on a target arising from acoustic radiation pressure, which is due to a change in momentum of the acoustic wave, when the wave interacts with such target. We here studied the capability of SonoLancet to generate acoustic radiation pressure in a suspension of microparticles. Two main jets ( Fig. 2A) of microparticles (models for micro-vehicles), were emanating from opposite sides of the needle bevel with a maximum velocity of approximately 1 mm s −1 at 0.6 mm from the needle tip, decreasing to 0.5 mm s −1 at a distance of 2 mm away from the maximum. Considering the wavelength, the small object size and similar acoustic impedance of particles compared to that of water, the main mechanism of translation arises from acoustic streaming, i.e. acoustic radiation force exerted on the liquid, rather than acoustic radiation force directly pushing the particles. This is supported by a finding that SonoLancet in water ( 22 • C ) without microparticles induced acoustic streaming; such observation was confirmed under Schlieren imaging (Fig. 2B), which revealed a disturbance of water laminar flow ( 15 • C ) close to the sonicating needle. The results demonstrate that SonoLancet is capable of inducing acoustic radiation pressure leading to mass transfer near the needle tip, which could be relevant to Atomization. Capillary waves propagating at the water-air interface can induce micro-droplet formation at large displacement amplitudes 36 . To investigate the potential of SonoLancet to produce micro-droplets, we first introduced a water droplet near the needle tip, while the needle was operated in air. Capillary waves were observed at low driving power levels (electrical output power < 1 W ) of the transducer. At increased powers (electrical output power ∼ 1 W ) crests of the waves started to extrude, eventually forming droplets ejected from the water-air interface 36 . By adjusting the driving frequency, we modified the size distribution of the deionized water droplets ejected from the surface (Fig. 2C,D) from a median size of 39.5 µm down to 17.5 µm . The trend in median droplet size is in line with the theoretical prediction: where d p is the droplet size, σ is the surface tension, ρ is the density and f is the excitation frequency 37 . These results demonstrate the capability of SonoLancet to produce water micro-droplets with a controlled size distribution, relevant to e.g. delivery of therapeutic agents such as drugs or antiseptics in pulmonary conditions.
Influence of SonoLancet on a tissue phantom and tissue. As demonstrated above, the SonoLancet can induce nonlinear effects in well-defined media such as water and air. In the following, we demonstrate SonoLancet's influence on a tissue phantom and on excised tissue. Ballistic gelatin was selected as a model for visualizing the needle activity inside a soft-tissue-like material. We subjected a hydrogel tissue phantom (10% w/v porcine gelatin/deionized water, 22 • C ) to the SonoLancet in conjunction with fluorescent nanoparticles (diameter 63 nm) introduced to the proximity of the needle tip through the lumen. When a continuous wave of 31.7 kHz and total acoustic power (TAP) ∼ 0.8 W were applied to the needle, the generation of cavitation bubbles and acoustic streaming was optically confirmed. The boundary of the influenced volume extended outward from the needle tip over time (Fig. 3A). After sonication was applied together with fluorescent nanoparticles, a 3D reconstruction of the nanoparticle distributions characterized by optical projection tomography revealed a volume www.nature.com/scientificreports/ containing nanoparticles is localized near the needle tip (Fig. 3B). Dissection of the phantom was used to verify that the gel was liquefied. At the highest TAP levels, ∼ 0.8 W, we recorded a temperature rise of < 2 • C or < 4 • C during a 5 or 20 s sonication, respectively, suggesting that cavitation was the main mechanism of liquefaction 38 , rather than ultrasound-induced temperature rise. Even though the observed acoustic phenomena in gelatin may occur in a different fashion in an ex vivo tissue, this experiment gave us an insight on how the acoustic energy is manifesting and influencing the surroundings of the needle tip, when the needle has been confined by material with similar acoustic impedance to that of soft tissue. To demonstrate the influence of SonoLancet on tissue, we performed FNABs with and without ultrasound at the needle tip (ultrasound-enhanced fine-needle aspiration biopsy, USeFNAB) in different bovine tissues: liver, kidney, spleen and striated muscle. USeFNAB contributed to increased sample yield, whose weight increased with increasing TAP (p < 0.0125), being on average up to 4× the yield obtained with FNAB in liver, 5× in kidney, 3× in spleen and 6× in muscle (Fig. 3C). When ultrasoundenhanced biopsies were performed at the highest TAP (0.8 W), the average sample weight obtained in liver was 88.6 ± 9.3 mg (average ± standard deviation; n = 6) ( Fig. 3C(1)), 53.7 ± 8.6 mg (n = 6) in kidney ( Fig. 3C(2)), 62.2 ± 5.4 mg (n = 6) in spleen ( Fig. 3C(3)) and 11.9 ± 1.9 mg (n = 6) in muscle ( Fig. 3C (4)); the sample masses obtained with FNAB were respectively 22.7 ± 1.8 mg (n = 6), 11.7 ± 2.0 mg (n = 6), 16.6 ± 1.0 mg (n = 6) and 1.8 ± 0.6 mg (n = 6). Based on a histopathological evaluation, performed by a pathologist (K.P.H.P.) with > 40 years of histopathology and cytopathology experience, all tissue types included intact cells and tissue ( Fig. 4) containing histologically relevant structures that could be observed in all samples. A computer-assisted method was adopted to quantify the fragmentation in the liver samples (Fig. 5A). The area of tissue fragments (large regions of adjacent hepatocytes), cell clusters (small aggregation of hepatocytes), single cells and debris pieces (organic waste) detected in the control group slides (FNAB) were measured and compared to the experimental groups (USeFNAB, TAP levels: 0 W, 0.2 W, 0.5 W and 0.8 W). On average the total areas of tissue fragments, cell clusters, single cells and debris pieces were slightly influenced by ultrasound; however, they were not statistically significantly different from the respective areas of FNAB (p > 0.0125). At the highest TAP the relative total area of tissue fragments was lower ( 57.9 ± 9.5 % , average ± standard deviation, n = 6, p = 0.0306) as compared to FNAB ( 69.9 ± 5.8 % , n = 5) (Fig. 5B(1)). The relative total area of cell clusters and individual cells were not different from those of FNAB when the highest TAP was employed (p = 0.0656 and p = 0.8102, respectively) ( Fig. 5B(2,3)). The relative total debris area ( 17.2 ± 4.6 % , n = 6, p = 0.0131) was higher as compared to FNAB ( 11.0 ± 2.3 % , n = 5) (Fig. 5B(4)).
The results show a systematic increase in the sample mass with increasing TAP, while inducing a minimal effect on the morphological structure of the sample as compared to the FNAB approach, as a manifestation of the actuation of tissue induced by SonoLancet.  1)) liver, (C(2)) kidney, (C(3)) spleen and (C(4)) muscle at different power levels compared to the standard FNAB. The needle employed in FNAB and ultrasound-enhanced FNAB (USeFNAB) was the same 21G needle. In the bar charts, the bar height represents the mean of the data set and the error bar indicates the standard deviation. Individual data points are shown as white-filled circles.

Discussion
We successfully coupled flexural standing waves to a conventional medical needle. The converging structure of the needle bevel seemed acoustically quite optimal at the employed frequency, because it geometrically amplified the wave allowing to generate a highly localized and large displacements of the needle tip. The acoustic intensity of the longitudinal waves emitted from the needle extended beyond the linear regime to nonlinear regime, suggested by large displacements < 100 µm and the observed nonlinear effects. The maximum time-averaged acoustic intensity next to the needle bevel was more than twice compared to that along the lumen. Since the acoustic radiation pressure associated phenomena, such as acoustic radiation force or acoustic streaming, are proportional to intensity, the mass transfer mechanisms are pronounced at the needle tip. The long needle shaft has multitude of resonant eigen-frequencies, which allowed changing the operation frequency. This permitted controlling the size of the atomized droplets contrary to common ultrasound atomizers, which typically operate on a single frequency [39][40][41] . Cavitation being a threshold phenomenon allowed generation of cavitation at will. Importantly, this was spatially restricted to the proximity of the needle tip without observations of cavitation along the needle lumen, suggesting that these mechanisms can be made to localize to the proximity of the needle tip. The high accelerations > 20000× the gravitational acceleration of the bubble-water boundary suggests that extreme shear forces can be induced by this boundary to the surroundings, when cavitation is intended. In addition, the direct effect of the rapidly moving needle tip with sharp edges provide points to exert high stress fields on matter, such as tissue.
Minimally invasive high-frequency mechanical actuation at a known position inside the body has specific relevance to various medical applications. We exemplified the influence of the SonoLancet on tissue in the FNAB application. Depending on the tissue type, the tissue yield was increased 3-6× on average and with increasing Based on this finding, we believe that the ultrasonic actuation of the needle introduces an additional degree of freedom in the needle movement (flexural displacement of the needle tip) that introduces additional forces facilitating the tissue extraction in an analogous way to what the fanning technique would contribute to, but at a smaller length-scale and higher frequency. In addition, the studied NLU phenomena are expected to generate shear forces that contribute to the extraction of cells and tissue constructs from the target, which are then drawn into the needle by the low pressure applied with the syringe. We foresee this to have importance especially in cancer diagnostics, which increasingly relies on molecular assessments for therapy decisions. The failure rate of molecular assessment is up to 50% 14 arising mainly from inadequacy of the nucleic acids required for determining the cancer type. This could be resolved by obtaining a greater sample volume using SonoLancet, permitting www.nature.com/scientificreports/ access to greater nucleic acid quantity. The SonoLancet could contribute also to histopathological assessments, which suffer from failure rates up to 34% due to inadequacy of tissue yield 42,43 . Future studies should address the clinically relevant adequacy of the sample, its micro-level integrity and the influence of ultrasound on the needle tract.
SonoLancet could be also beneficial for other medical applications beyond cancer diagnostics. Controlled cavitation, used in a variety of existing and emerging ultrasonic medical technologies 44 , could be exploited in personalized therapeutic and diagnostic applications, by applying cavitation locally at the needle tip precisely at the target site. Such applications include tumor ablation 45 , permeation of tissue matrix for enhanced stem cell delivery or migration, bubble-enhanced poration of cells, i.e. sonoporation 46 , transfection of genetic material 47 or even modifying gene expression by stimulating ultrasound 48 . Moreover, cavitation-induced shock waves and associated shear forces, arising from imploding cavitation bubbles 49 , could be employed for lithotripsy 50 or controlled softening of calcified vessels 51 with spatially precise and localized mechanical stress gradients 52 . The SonoLancet could potentially serve as a way to induce streaming of biological fluids and translation of entities inside the body. Unidirectional motion of drugs or drug nanocarriers, induced by the transfer of momentum from the wave to the liquid, could enhance pharmaceutical effectiveness 53 and facilitate applications such as targeted and localized drug, stem cell or gene delivery. Atomization provides the ability to transform liquid into microdroplets with control over droplet size distribution and direction of the droplet jets allowing new platforms for delivery of drugs or cells e.g. within bodily cavities such as those in the respiratory system.
Considering the safety, the diameter of the volume of tissue influenced by SonoLancet is comparable to that of a 11G needle (outer diameter = 3.05 mm ), while the energy is delivered through a small tract generated by a 21G needle (outer diameter = 0.82 mm ). This allows influencing a large volume of tissue in respect to the small size of the needle, while minimizing trauma, when accessing target site. Therefore, the SonoLancet could contribute to modifying the tissue permeability for drug delivery applications at precise tissue locations and with extreme ultrasound exposure, or allow minimally invasive surgical interventions via spatially precise histotripsy. Importantly, in all of these experimentations the employed ultrasound frequencies were relatively low, < 100 kHz . Moreover, the risk of cavitation to occur and induce mechanical damage in unwanted regions inside tissue (e.g. along the needle shaft) is minimized by the fact that the acoustic energy is highly concentrated at the tip of the needle, as supported by numerical simulations (Fig. 1B(1,2)). For this reason, cavitation is more likely to appear at the needle tip than elsewhere, and, if such event were to occur along the needle shaft, its effects would be limited in a confined region from the needle boundary, as the intensity attenuates strongly in the direction radial to the needle center axis due to geometric factors. Since at these frequencies the attenuation coefficients in tissue and water are small, the heat deposition to tissue is minimized. This is supported by low temperature rise recordings < 2 • C and in extreme conditions < 4 • C . Therefore, these findings suggest that the SonoLancet provides a minimally invasive way to influence matter via nonlinear ultrasonic mechanisms within a confined volume at the proximity of the needle without major thermal impact. To our knowledge, this study is the first to demonstrate the selected NLU effects near the ultrasonically actuated medical needle. While this study exemplified the use of ultrasonically actuated needle in the biopsy application, its use could extend to other medical applications.
To conclude, we have functionalized the medical needle tip by nonlinear ultrasonics and, at a micrometer to millimeter scale, and quantitatively characterized its physical properties such as cavitation, acoustic radiation pressure and atomization within a confined volume defined as the SonoLancet. In this approach, the conventional hypodermic needle structure employs the lumen as a conduit to bring sound energy to the needle tip. The 3D shape of the conventional needle bevel amplifies the ultrasound action resulting in localized nonlinear effects at the proximity of the needle tip. Limiting the ultrasound action to the SonoLancet provides spatiotemporal control of sound-tissue interaction and with respect to minimizing microtrauma of adjacent tissues. By bringing precisely controlled ultrasound energy to a specific location in the body, the SonoLancet has the potential to give conventional medical needles a new role in many applications, such as tissue biopsy, molecular diagnostics, drug or gene delivery, cell modulation, or minimally invasive surgical procedures. Numerical simulation. In this study, the software COMSOL Multiphysics v5.5 54 was used to simulate the acoustic emission of the needle in a volume of water. The 3D geometry consists of an ultrasound transducer coupled to a 21G hypodermic needle via an S-shaped aluminum waveguide (Fig. 1A(1)). The needle shaft is surrounded by a cylindrical volume of water ( external dimensions = R × H = 5 mm × 50 mm ). The ultrasonic actuation of the needle is obtained by applying a potential difference of 100 V at a frequency of 31.32 kHz across the piezoelectric stack. The pressure field in water, generated by the vibrating needle, was calculated by performing a frequency domain analysis. The x-component of the time-average acoustic intensity, defined as I ta,x = 1 2 Re{pv * x } with p being the complex pressure and v * x the complex conjugate of the acoustic velocity along the x-direction, was then evaluated on the xy-plane coincident with the needle center axis (Fig. 1B(1)) and on a line adjacent to the outer needle surface, starting from the needle tip and ending at the waveguide-needle attachment ( Fig. 1B(2)).  (Fig. 6B). The needle was immersed to a depth of 45 www.nature.com/scientificreports/ mm from the water surface and subsequently actuated by a continuous sinusoidal signal in the frequency range from 31 to 32 kHz and emitted TAP of ∼ 0.2 W. HS recordings were taken in order to visually capture the cavitation activity around the needle at different positions along the needle length. More specifically 45 videos (sample rate = 300,000 frames per second (fps), exposure = 2.8 µs , resolution = 128 pixels × 112 pixels, lens aperture = 3.3; lens model: Canon MP-E 65 mm) were acquired, covering a total length of 45 mm from the needle tip to the waveguide after vertical concatenation. These settings permitted to achieve an image resolution of about 9 µm/pixel . The HS videos were finally analyzed by an algorithm implemented in MATLAB (R2020b) 55 capable of tracking the needle position over time and quantify the cavitation activity. The needle displacement along the x-axis from its rest position was measured using sub-pixel image registration by cross-correlation in the frequency domain 56 between the i th frame I i (x, y) and the first frame I 1 (x, y) of the high speed footage as follows:

Cavitation. Cavitation
where x i is the measured needle x-displacement from the rest position and FFT is the two-dimensional fast Fourier transform. In order to quantify the projected area of cavitation a new set of frames, where only the cavitation activity is present, is created. This is done by subtracting the reference image I 1 (rigidally translated by x i in order to match the needle position in the i th frame) from each frame of the original footage: where iFFT is the two-dimensional discrete inverse Fourier transform. A probability map P cavitation (x, y) , in which the value of each pixel represents the probability to observe a cavitation event at a specific coordinate, is then measured as follows: where N is the total number of frames and I cavitation,bw (x, y) is a frame-set generated by thresholding I cavitation (x, y) with the Otsu method 57 . The projected area of cavitation as a function of time is then obtained: Finally, a total of 45 probability maps were vertically concatenated to reproduce the entire length from the needle tip to the waveguide.
Atomization. The purpose of this experiment was to demonstrate that with an ultrasonically actuated needle one can eject micro-droplets from the needle tip surface and control the droplet size distribution. A droplet of deionized water was introduced near the needle tip, while the needle was operated in air (Fig. 6C). HS recordings of the atomization process were acquired with the following settings: sample rate = 16653 fps, exposure = 10 µs , resolution = 768 pixels × 768 pixels, lens aperture = 16 (lens model: Canon MP-E 65 mm). The video frames were then analyzed in MATLAB to determine the droplet size distribution for different excitation frequencies and to compare it with theoretical predictions. Each frame I i (x, y) was first segmented with an adaptive thresholding 58 . Then, the built-in MATLAB function 'regionprops' was applied to measure the area and the equivalent diameter for each 8-connected component in the binary image I bw,i (x, y) . Finally, the median droplet diameter distribution was calculated over N = 1000 randomly selected frames and compared with theoretical predictions defined by Eq. (1). Acoustic streaming. Polystyrene based microparticles with a diameter of 30 µm (catalogue number: 84135-5ML-F, Micro particles based on polystyrene analytical standard, size: 30 µm , Sigma-Aldrich, St. Louis, MO, United States) were used to visualize their translations near the needle tip. The test was conducted in a plastic cuvette ( external dimensions = L × W × H = 12 mm × 12 mm × 45 mm , wall thickness = 1 mm) filled with deionized water (3 mL) with a particle concentration of 0.05% v/v (Fig. 6D). A low power continuous sinusoidal signal (TAP = 0.1 W) was applied to the needle, whose action was recorded with the following care settings: sample rate = 10,000 fps, exposure = 2 µs , resolution = 1024 pixels × 800 pixels, lens aperture = 2.8 (lens model: Canon MP-E 65 mm). An averaged velocity field V (x, y) =v xx +v yŷ , being x and ŷ unit vectors, was then calculated by employing a cross-correlation technique 59 . Each frame I i (x, y) is divided into small interrogation boxes of the size of 64 pixels × 64 pixels. Considering I i (x, y) and I i+1 (x, y) to be the frames recorded at time t and t + dt , let f(x, y) be an interrogation box centered at specific location (x, y) in I i (x, y) and g(x, y) an interrogation box centered at specific location (x, y) in I i+1 (x, y) . The displacement field u(x, y) of two consecutive frames is calculated by measuring the offset (�x, �y) of the peak of the correlation map between f(x, y) and g(x, y) as follows: The averaged velocity map is then calculated as: The needle was inserted to a depth of 7 mm from the gelatin surface, penetrating three times the sample in order to let the contrast agent wet the walls of the created cavity. Continuous ultrasound wave at 31.4 kHz and TAP ∼ 0.8 W was applied to the needle for 10 s. HS videos of projections of the volume of interaction between the gelatin and the SonoLancet were taken using the following settings: sample rate = 1000 fps, exposure = 10 µs , resolution = 768 pixels × 768 pixels , lens aperture = 16 (lens model: Canon MP-E 65 mm). Optical projection tomography scanning of the sample was performed using an OPT scanner (Bioptonics OPT 3001M Scanner). Reconstructed images (resolution = 512 pixels × 512 pixels), which represent the intersection between planes perpendicular to the needle center-axis and the volume of actuation in the tissue phantom, were generated using NRecon version 1.6.1.0 (Skyscan) software and further analyzed in MATLAB. Each reconstructed image I z (x, y) was first thresholded with the Otsu method. A sequence of morphological erosion and dilation using a diskshaped structuring element of radius r = 10 pixels was then applied to fill eventual gaps in the image and remove small unwanted objects arising from noise. Contours of the sections, defined as sequences of n contour points {P 0 , P 1 , . . . , P n−1 } , are extrapolated from the images and filtered with a Savitzky-Golay finite impulse response (FIR) smoothing filter (polynomial order = 2, window size = 11). Finally, the 3D surface of actuation around the needle tip is constructed with triangular elements and further adjusted in Fusion 360.

USeFNAB.
We quantified the ultrasound-tissue interaction by measuring the tissue yield of ultrasoundenhanced fine-needle aspiration biopsies obtained from bovine liver, kidney, spleen and muscle at 4 different TAP levels (0 W control; sonication groups: 0.2 W, 0.5 W, 0.8 W) with the yield obtained with the conventional fine-needle aspiration biopsy approach (FNAB control) having the same syringe (10 mL) and needle (21G × 80 mm) as USeFNAB. Liver (from 33 months old female bovine), kidney (from 33 months old male bovine), spleen (from 55 months old female bovine) and striated muscle (diaphragm, from 19 months old female bovine) specimens were obtained from the slaughterhouse (Vainion Teurastamo Oy, Orimattila, Finland) within 1 h post mortem and experiments were performed 3-5 h post mortem at room temperature (22-24 • C ). Only peripheral parts of the same liver were used, to avoid large blood vessels and acquisition of bloody samples. In liver and in spleen, samples were taken at a depth of 2 cm from the organ surface by keeping the needle parallel to the table, where the target tissue is fixed throughout the procedure. Sample acquisition in kidney was limited to the cortical area and each biopsy was performed in different renal lobes, while, in muscle, biopsies were taken along the fibre orientation. The used needles and syringes were weighed before and after the biopsies using a balance with resolution of 0.1 mg (model: Quintix Precision Balance 5100g × 10 mg, Sartorius AG, Germany) to measure the mass of the sample extracted from the tissue. The test was divided into 6 sets of 5 repetition. In each set, the order of the groups to be tested (FNAB (n = 6); USeFNAB: 0 W (n = 6), 0.2 W (n = 6), 0.5 W (n = 6) and 0.8 W (n = 6)) was randomized with a random permutation within a set of five samples (one sample from each group). The FNAB procedure consisted of inserting the tip of the needle into the target tissue using a 'fanning technique' , pulling the plunger back to the 2 mL mark and moving the needle back and forth 7 times with a frequency of 1 pass per second. The penetration depth was approximately 15 mm at a targeted angular offset of approximately 5 • between the passes. The suction was maintained manually throughout the biopsy sampling by keeping the plunger of the syringe at the the 2 mL mark and only released after the needle was withdrawn 60 . The same procedure was adopted for USeFNABs, where ultrasound continuous waves (f = 31.7 kHz, pulse repetition frequency (PRF) = 55 Hz, duty cycle (DC) = 52.05%) were applied throughout the biopsy after the needle was inserted into the tissue. The selected parameters of the protocol were considered appropriate to obtain a sufficient amount of sample for yield analysis or to be processed into agarose cell blocks. The operator skills in performing the biopsy and the protocol were confirmed by a specialist in radiology (J.K.) with considerable experience in FNAB. Assuming that normal distribution does not apply, the FNAB mass was compared to the masses of USeFNAB groups using Mann-Whitney U test. The Bonferroni-corrected level of statistical significance was 0.0125. and subjected to centrifugation at 500 G for 5 min. After a second PBS washing, pellets were resuspended in 500 µL of 10% neutral buffered formalin (15817114, HistoTainer, Simport, Saint-Mathieu-de-Beloeil, Quebec, Canada). After 10 min of incubation time, tubes were centrifuged at 500 G for 5 min, then, the supernatant was carefully discarded, and pellets were washed in PBS twice. Excess PBS was carefully removed and formalin fixed pellets were resuspended in 100 µL of molten agarose (10377033, Agarose Low-Melting, Nucleic Acid Recovery/Molecular Biology Grade, Thermo Fisher Scientific, Waltham, MA, United States) kept at 42 • C and gently centrifuged at 180 G for 1 min. The agarose cell blocks were let solidify at 4 • C and then transferred into 70% v/v ethanol (172381, 70% v/v Ethanol Solution Molecular Biology Grade, Fisher Scientific, Hampton, NH, United States) pre-filled tubes. Several sections of the paraffin embedded cell blocks (4 to 5 for each sample) were cut at 5 µm thickness, mounted on 25 mm × 75 mm slides (J1800AMNZ, SuperFrost Plus Adhesion slides, Thermo Fisher Scientific, Waltham, MA, United States) and stained with hematoxylin and eosin (H&E). Finally, microscopy pictures (bright field, 20× magnification, extended focus) were generated using 3DHISTECH Panoramic 250 FLASH II digital slide scanner. Four biopsy samples (liver, n = 1; spleen, n = 1; muscle, n = 2) failed to be processed for histological analysis due to breakage of the agarose cell block.
Quantification of sample fragmentation. A supervised classification method has been implemented in MATLAB to quantity the fragmentation of the liver samples by means of feature extraction of objects identified in the histological images 62 (Fig. 7). The classification algorithm is able to detect four different histological structures, namely tissue fragments (large regions of adjacent hepatocytes), cell clusters (small aggregation of hepatocytes), single cells and debris (organic waste). The dataset is composed of 29 histological slides belonging to the following sample groups: FNAB (n = 5); USeFNAB: 0 W (n = 6), 0.2 W (n = 6), 0.5 W (n = 6) and 0.8 W (n = 6). Each image I h (x, y) was first converted into gray-scale image and then subjected to thresholding using the Otsu method. The the Moore-Neighbor tracing algorithm modified by Jacob's stopping criteria 63 was applied to trace the boundaries of all 8-connected objects detected in the binary image I bw,h (x, y) . The objects were assigned to a certain class based on a similarity score method defined as follows: www.nature.com/scientificreports/ In the above, C is the the output class which each object is assigned to by the classifier, x i is a vector containing a number M of morphological features (i.e. area, eccentricity, circularity, solidity, etc.) extracted from the i th object and k is an integer number corresponding to a certain class (1 = tissue fragments; 2 = cell clusters; 3 = single cells; 4 = debris). S(x i , k) is the score function whose output value denotes how well an object is represented by a class k (the lower the score, the higher the similarity). Y k is a M × N matrix containing feature vectors belonging to N pre-classified objects that were priorly assigned to a class k in the supervised learning phase. Finally, several indicators of sample fragmentation (i.e. total number of tissue fragments, cell clusters, single cells, debris pieces and their respective average area) were extracted and used for statistical comparison between the FNAB and USeFNAB groups. Assuming that normal distribution does not apply, the Mann-Whitney U test was chosen, with a Bonferroni-corrected level of statistical significance of 0.0125.

Data availability
The datasets are available upon request.

Code availability
The codes are available upon request.