Characterizing cellular morphology by photoacoustic spectrum analysis with an ultra-broadband optical ultrasonic detector

: Photoacoustic spectrum analysis (PASA) has been demonstrated as a new method for quantitative tissue imaging and characterization. The ability of PASA in evaluating micro-size tissue features was limited by the bandwidth of detectors for photoacoustic (PA) signal acquisition. We improve upon such a limit, and report on developments of PASA facilitated by an optical ultrasonic detector based on micro-ring resonator. The detector’s broad and flat frequency response significantly improves the performance of PASA and extents its characterization capability from the tissue level to cellular level. The performance of the system in characterizing cellular level (a few microns) stochastic objects was first shown via a study on size-controlled optically absorbing phantoms. As a further demonstration of PASA’s potential clinical application, it was employed to characterize the morphological changes of red blood cells (RBCs) from a biconcave shape to a spherical shape as a result of aging. This work demonstrates that PASA equipped with the micro-ring ultrasonic detectors is an effective technique in characterizing cellular-level micro-features of biological samples.


Introduction
Photoacoustic Imaging (PAI) is a hybrid medical imaging modality based on laser-induced thermo-acoustic effect in biological tissues. PAI combines the advantages of both optical and ultrasound imaging as it is able to present optical contrast in subsurface tissues with an excellent spatial resolution. A majority of previous studies on PAI were focused on measurements of photoacoustic (PA) signal amplitude which is related to the overall optical absorption of the target tissue. Recent studies have demonstrated that the power spectrum of the broadband radio frequency (RF) PA signal contains the micro-structural information of the tissue [1][2][3][4][5][6][7][8][9]. Aiming at quantitative imaging and tissue characterization, a technique termed 'photoacoustic spectral analysis (PASA)' has been developed. Similar to ultrasound spectral analysis (USSA), PASA uses the slope, mid-band fit, and intercept of the linear regression model to characterize the main features of the signal power spectrum [4]. PASA is based on the fact that the frequency components of the PA signals are closely correlated to the morphological properties of optically absorbing objects in tissues, including their sizes, shapes, and densities [4,10]. For example, smaller objects generate shorter PA signals in the time domain, which leads to broader power spectra containing more high frequency components than the PA signals from larger objects. Quantitative analysis of the PA signals in the frequency domain would enable an objective evaluation of the morphologies of optically absorbing objects stochastically distributed in biological tissues [3,11]. PASA offers solutions to a number of practical issues faced by conventional PAI. First, the averaged power spectrum provides a robust method of quickly characterizing the stochastic nature of tissue microstructures, and can lead to measurements that are more repeatable. Second, PASA also minimizes the contributions of system components, and leads to measurements that are more objective and quantitative.
Previous PASA was based on conventional piezoelectric transducers, which usually had limited detection bandwidth. For example, by using commercially available or homefabricated piezoelectric transducers, PASA has been explored for applications in diagnosis of prostate cancer [6,12,13] and evaluation of liver conditions [7,8,14]. In these studies, PASA was performed on the tissue level, and the targeted biological tissues (e.g., small vessels, cellular groups, and lipid clusters) have features on the scale of over 100 microns. In order to achieve more sensitive and accurate evaluation of cells with sizes on the order of a few microns, the detector(s) for PA signal acquisition should have a bandwidth of a few hundreds of MHz. Photoacoustic microscopy systems with specially-designed high-frequency piezoelectric transducers have been used to study the morphology of single red blood cells [9, [15][16][17]. Detecting high-frequency ultrasound signals (e.g., ~hundreds of MHz) usually requires carefully manufactured transducers with thin and fragile piezoelectric elements. Besides, due to the resonance effect in the piezoelectric part, these transducers usually have maximum responses at a high frequency, and therefore, a reduced response at the low frequency regime [18][19][20].
In this work, we demonstrate, for the first time to the best of our knowledge, an all-optical PASA system using a polystyrene (PS) micro-ring resonator as the ultrasonic detector [21]. Micro-ring resonator has been demonstrated as a high-performance ultrasound detector with both broad bandwidth and great detection sensitivity. Our previous studies have shown that the micro-ring detector has an intrinsic acoustic response from nearly DC up to 350MHz at −3dB [22], and a noise equivalent detectable pressure over this bandwidth as low as 105 Pa. This sensitivity is comparable to or even better than those of conventional piezoelectric transducers. PA signals from biological samples contain information of both micro-and macro-size structures, and are broadband in nature. Micro-ring ultrasonic detectors can largely improve the performance of PASA for tissue characterizations by extending its scope from the tissue level to cellular level. With such a system, we can achieve size characterization of PS microspheres below 20 µm, and can differentiate cells with different morphologies.

Experimental setup and phantoms preparation
The setup for our PASA system for tissue characterization is shown in Fig. 1(a). Laser pulses from an optical parametric oscillator (OPO) pumped by a Nd: YAG laser (Vibrant B, Opotek) were used to illuminate the sample and generate PA signals. The pulses had a wavelength of 532 nm, pulse duration of 5 ns, and repetition rate of 10Hz. To achieve a relative uniform light illumination (laser fluence in the range of 15-20 mJ/cm 2 ), the light beam from the laser was first focused by a lens with shorter focal length. At the focal point, a small aperture (a pinhole) was used to filter the high spectral frequency component. Then the light beam was expanded by the second lens with longer focal length. A second larger aperture was added to further reduce the noise at the edge of the beam. Each sample held in a plastic container was made from 10% porcine gel. The samples contained stochastically distributed polystyrenebased latex microspheres or red blood cells (RBCs). The container had openings at top and bottom. The PA signals propagated downward to the ring and upwards to the hydrophone. The micro-ring placed underneath the sample was covered by a thin light-blocking film (silver Mylar balloon) to protect from being damaged by the scattered laser light. DI water was injected between the micro-ring and the film for better coupling of acoustic signal. To compare the performance of the micro-ring to conventional piezoelectric transducers, a calibrated needle hydrophone (HNC-1500, Onda) with a detection band from 0 to 37 MHz (−10dB) was placed on the top of the sample for measurement of PA signals. The distances from the center of the laser beam to the two ultrasonic detectors (i.e., the micro-ring and the hydrophone) were kept the same (3 mm). The PA signals received by the micro-ring and the hydrophone were recorded at the same time by a digital oscilloscope (TDS 540B, Tektronics, sampling frequency 500MHz) and averaged 300 times to increase the signal to noise ratio (SNR). To explore the ability of PASA equipped with a micro-ring, experiments on phantom and biology samples were conducted. The phantom samples we used were microsphere particles (Fluoresbrite plain YG micro-spheres, Polysciences) with sizes of 3 µm, 6 µm, 10 µm, 20 µm, and 45 µm. With these phantoms, we explored the capability of PASA in characterizing the size information of the optically absorbing objects. The biological samples we used were two groups (fresh and aged) of red blood cell samples, one using fresh human blood specimens (< 2 weeks in storage) and the other using aged human blood specimens (3 months in storage).
To prepare the phantoms with stochastically distributed microspheres or blood cells, we mixed the samples in porcine gel. The densities were controlled so that the generated PA signals had similar signal noise ratios. No aggregation of either microspheres or RBCs happened, as confirmed by checking the samples under a microscopy. The properly controlled density also allowed light beam to pass through the entire phantom without significant attenuation, so that the illumination intensity is relatively uniform over the whole phantom. The microspheres were diluted by the porcine gel directly. The blood samples were diluted 1:100 in phosphate buffered saline (PBS) and then mixed 1:30 with 10% porcine gel. The microspheres were added at 46°C and the blood cells were added a temperature below 37°C and the gel was stirred using a magnetic stirring bar to achieve stochastic distribution of optical absorbers. After the gel was cooled and solidified, a phantom containing randomly distributed objects were ready for use. The densities of microsphere particles in the gel were 1.7 × 10 6 , 2.1 × 10 5 , 4.6 × 10 4 , 1.9 × 10 4 , and 1.7 × 10 3 particles/mL, respectively, for sizes of 3 µm, 6 µm, 10 µm, 20 µm, and 45 µm. The densities of RBCs in the gel were 1.2-1.8 × 10 6 cells/mL.

Principles
The details of the micro-ring based ultrasonic detectors have been discussed in our previous publications [24,25]. A micro-ring resonator consists of a ring-shaped waveguide coupled with a straight bus waveguide and there is light coupling in and out of the ring in the gap region. The output intensity spectrum shows dips occurred at specific wavelengths when the corresponding round trip phase accumulation in the ring waveguide equals to multiples of 2π radians. This round trip phase accumulation depends on many parameters, such as the ring radius, cross-section size, and refractive index of the polymer waveguide core/cladding, all of which can affect the effective index of the propagating mode. For polymer micro-rings in contact with ultrasound waves, the intrinsic refractive index will be modulated through the elastic-optic effect [26]. At the same time, the ring is deformed by the ultrasound waves, affecting the effective index of the propagating waveguide mode. These two effects altogether shift the ring resonant wavelength. By setting the probe laser (TLB 6700, Velocity Tunable Laser) wavelength at the slope of the ring transmission dip, the incident ultrasound waves could be recorded by the varying output intensity detected by a high speed photodetector (125MHz, Model 1801, Newport) from the device. Micro-ring has been utilized in highresolution PA imaging [27,28], as well as real-time THz pulse detection [29]. The device has demonstrated an ultra-broadband intrinsic frequency response from nearly DC to 350 MHz at −3dB [22]. However, the whole system response bandwidth is not only limited by the bandwidth of the photodetector but also by the duration of the laser pulses. To evaluate the impulse response of the PASA system equipped with the micro-ring, the 5-ns laser pulse was illuminated on a super-thin Cr film (200-nm thick), and then the generated PA signals were detected by the micro-ring placed at 3 mm from the Cr film. The PSD of the PA signal from the Cr film is shown in Fig. 1(b) which illustrates a system bandwidth of 123 MHz at −10dB.
The working principle of PASA is based on a linear regression on the power spectrum of PA signal from a target sample. The power spectrum represents the distribution of power along the ultrasound frequency, which is calculated as the square of the fast fourier transform of the time-domain PA signal magnitude. PA signals from spherical objects have been derived analytically and the power spectrum of those time-domain bipolar signals show periodic dip pattern at certain frequencies [30]. A general approach to differentiate two spheres with different sizes is to identify the locations of the minima in power spectrum. This method will become invalid when sample the concentration of optically absorbing objects goes up and then the details in the spectrum such as the minima in the spectrum becomes less obvious to observe, leaving only a contour of the curve. In most cases, the power spectrum over frequency show an overall linear descending trend, and a simple and robust linear regression can be used to evaluate this trend when a broad enough spectrum is fitted. The linear fitting is sensitive to the ratio between the high frequency and the low frequency components which is directly related to the morphologies of the microstructures in the tissue. In order to make this linear fitting accurate, the bandwidth of the ultrasound detector must be broad enough so that the descending contour of the PSD can be captured. In addition, to reduce the impact of individual measurement system to the quantified results, the system impulse response needs to be removed from the PSD before linear regression is conducted. Among the three spectral parameters quantified from the linear fitting, the slope of the fitting line is the one most sensitive to the size and the shape of the optically absorbing objects in the sample, and is less susceptible to their concentration.

Simulation
For individual microspheres and RBCs experiments, we compared the measured spectrum with the spectrum simulated by the k-wave toolbox in Matlab described in [30].
For individual microspheres, we simulated the PA spectrums of the microspheres of sizes 20 μm, 45 μm and 100 μm. For RBCs samples, the biconcave shape used to describe healthy RBCs has a diameter of 7.82 μm, a minimum height of 0.81 μm, a maximum height of 2.58 μm, and a volume of 94 μm 3 [23]. When RBCs are aged, their shapes are changed to spheroid, while their volumes are kept the same [31]. Both samples contained stochastically distributed RBCs with a concentration of 1.7 × 10 6 cells/mL. The concentration of RBC here is consistent with the concentration in the experiment. Due to the gravity and the bi-concave shape of the RBCs, their orientations in the space were not uniform. As observed under the microscope, more RBCs had their long axis oriented along the horizontal plane, and fewer oriented along the vertical direction. Compared to uniform distribution, Gaussian distribution better describes the orientations of the RBCs. This was considered in the simulation. The density, the speed of sound (SOS), and the ultrasound attenuation of the sample were set as 1000 kg/m 3 , 1540 m/s, 0.2 dB/(MHz·cm), respectively. For each group (i.e., biconcave and spheroid), 100 samples (N = 100) were studied for statistical analysis.

PASA of individual microspheres
The performance of the PASA system was first evaluated by recording the PA signals from individual microspheres (Fluoresbrite plain YG micro-spheres, Polysciences) of different sizes (i.e., 20 µm, 45 µm, and 100 µm in diameter, respectively). Example PA signals from the individual microspheres with different sizes are shown in Figs. 2(a)-2(c). Each one shows a characteristic bipolar shape. The normalized PSD for microspheres with three different sizes are shown in Figs. 2(d)-2(f) with the solid curves. In this study, facilitated by the ultra-broad detection band of the micro-ring, the PSD from 20 µm, 45 µm and 100 µm microspheres were in reasonably good agreement with the results from simulation, as shown by the dashed curves in Figs. 2(d)-2(f). Experimental results show the first "dips" in the PSD profiles at 113 MHz, 48 MHz, and 26 MHz for 20-μm, 45-μm, and 100-μm microspheres, respectively, which are relative close to the numbers from the simulation (i.e., 100 MHz, 50 MHz, and 23 MHz, respectively). The mild disagreement between the results from simulation and experiment is at least in part due to the difference between the claimed size and the actual size of the microspheres. For example, according to the information from the manufacture, the microspheres with a mean size of 45 μm have a coefficient of variation of up to 20%, meaning the actual size can be anywhere between 36 and 54 μm.

PASA of stochastically distributed microspheres
The phantoms were prepared using the methods described in section 2.1. For each size, a group of 4 phantoms were measured. PA signals from each phantom were recorded by the micro-ring and the needle hydrophone simultaneously. The linear regression was conducted within the −10 dB spectral range of each system (micro-ring vs. hydrophone). In addition, The PSD below 2MHz was excluded in the linear regression in order to remove the interference from the low frequency noise. This low frequency noise was from the PA signals generated by the light illumination on the the gel. Therefore, the frequency range of linear regression fitting was 2-123MHz for the micro-ring based system, and was 2-37 MHz for the hydrophone based system. Before performing the linear fitting, the impulse response of each system was removed from the PSD. The spectral parameter slope of each fitting line was quantified, and the measures from 4 phantoms in each group were averaged.
Example PA signals from a phantom containing 20 µm microspheres as measured by micro-ring and hydrophone are shown in Figs. 3(a) and 3(d), respectively. The normalized and averaged PSD curves of the PA signals from the phantoms with different sizes of microspheres are plotted in Figs. 3(b) and 3(e). The difference in the descending rate of PSD curve corresponding to the difference in object size from 3 to 45 µm is well demonstrated by the measurements using the micro-ring detector. In contrast, the measurements from the hydrophone which has a narrower detection bandwidth are not able to differentiate the phantoms containing microspheres with such small sizes. The qualified spectral parameter slopes for the results from the micro-ring and the hydrophone are compared in Figs. 3(c) and 3(f). It can be seen that PASA using the hydrophone can hardly distinguish the phantoms containing microspheres smaller than 20 μm. In contrast, PASA using micro-ring shows an excellent capability in differentiating the phantoms with microsphere sizes ranging from 3 to 45 µm. This study on phantoms with well-controlled particle sizes indicates that micro-ring's ultra-broadband response facilitates reliable PA characterization of stochastically distributed micron-size objects.

PASA of ex vivo human blood specimens
In this study, we explored the feasibility of PASA in evaluating the morphological change of RBCs as a function of the blood storage time. The procedure of PASA was the same as that in the experiments on microspheres. Healthy RBCs have biconcave, disc-shape and therefore have larger surface area to volume ratio. Previous studies have shown that the geometry of RBCs becomes more spherical with storage, and therefore have a smaller surface area to volume ratio [32]. This morphology change can be reflected in the power spectrum of the PA signals. Before the experiment, we performed simulation using the toolbox mentioned in section 2.3. Two groups of samples each containing different shapes of RBCs (biconcave and spheroid) reflecting different levels of freshness (fresh and aged) were studied. The detailed geometries of biconcave and spheroid RBCs are shown in Fig. 4(a) and simulation details are discussed in 2.3. Similarly, the impulse response of each system was removed from the measured PSD before doing the linear fitting. The PA spectral parameter slopes for the two groups, as quantified in the simulation, are compared in Fig. 4(b). For the fresh RBCs, most of them sit flat in each sample. For these RBCs, what measured by the micro-ring located at the bottom of the sample were their small vertical thicknesses which produced large high frequency components. For the aged RBCs, even though their diameters appear smaller in Fig. 4(c), their vertical thicknesses were actually larger than those of the fresh ones. This difference in shapes explains why the aged RBCs, although appearing smaller under the microscope, produced weaker high frequency components and, therefore, smaller spectral parameter slope. In the experiment, two groups of samples each containing different shapes of RBCs (biconcave and spheroid for fresh and aged blood) were studied. The RBCs in both specimens were photographed under a microscope (FMA050, Amscope), as the example pictures shown in Fig. 4(c). The sample preparation was described in section 2.1. Figs. 4(d) and 4(e) shows an example PA signal and PSD curve from each group. As expected, aged RBC sample has high frequency component in the PSD. By performing the linear regression fit of the power spectrum in the range of 2-123 MHz, the quantified spectral parameter slopes from the two groups were obtained, as compared in Fig. 4(f), where aged RBC sample showed reduced slope. This experimental result matches with the finding from the simulation.

Discussion
The attenuation of the ultrasound and light could affect this linear fitting model's accuracy by modifying the slope, especially for dense medium. However, in our experiments, we controlled the concentrations of the optical absorbers (either microspheres or RBCs) low. This ensured that all the samples were highly optically transparent. Also, we controlled the uniformity of light illumination on the sample. Therefore, we don't expect the small optical attenuation will affect the results. To achieve similar ultrasound attenuation for all measurements, we measured all the samples at a same distance, and the gel concentration in the samples was the same. Since the microspheres occupied <0.01% of the sample volume, they would not strongly affect the ultrasound attenuation. As a result, all samples had the same attenuation with that of the gel itself.
The findings from this study on well-controlled phantoms and ex vivo biological tissues suggest that PASA based on micro-ring ultrasonic detectors, which have an ultra-broad and flat detection band, has the ability of characterizing morphological features including both size and shape in biological samples on the micron-size level. Besides, compared to extremely high-frequency piezoelectric transducers that usually involve sophisticated fabrication and, hence, are usually expensive, PS micro-ring is fabricated with much lower cost and better repeatability using nano-imprinting technology [25,33]. The micro-ring based PASA technique holds potential to detect morphological abnormalities of cells and relate them with corresponding pathological conditions. Unlike existing technologies such as flow cytometry and peripheral blood smear examination, which are based on the analysis of individual cells, PASA is essentially a statistical method and, hence, may better deal with the stochastic nature of biological samples. Moreover, PASA technique could be further developed into a fastspeed cell morphology analysis system for clinical use, and could provide quantitative diagnostic information that is not accessible by the current qualitative imaging methods. For example, certain diseases, infections, or exposure to toxic chemicals can alter the morphological features of RBCs, impairing their functionalities. Since various morphological abnormalities of RBCs are closely correlated with specific disease states, previous studies have been focused on developing techniques that can characterize the morphological changes in RBCs in vivo or ex vivo [30,[34][35][36]. Therefore, one potential application of PASA, as explored initially in this study, is the in vivo diagnosis of RBC-related diseases or ex vivo characterization of donated blood. Besides RBCs, other cells with strong intrinsic optical absorption contrast, such as melanoma cells, could also be evaluated with PASA working at proper optical wavelengths.
In the future, employment of a PA light source providing a shorter pulse duration and a photodetector with a larger bandwidth could further increase the bandwidth of the micro-ring based PASA system up to 350 MHz. Working at such a high frequency, however, the penetration depth in biological tissues is limited due to the large acoustic attenuation in tissues. Therefore, the applications of cellular-level PASA should be either focused on superficial tissues, such as illness detection of mucous membrane in the gastrointestinal track, or realized via a miniaturized PA device that can be placed in deep tissue in a minimally invasive manner. For example, when a PA probe can be miniaturized and integrated into a biopsy needle, early detection and characterization of cancer cells may be achieved. All of these potential applications could benefit from the all-optical design facilitated by the microring detector.