1.IntroductionIn recent years, an increasing number of studies have shown that tumors are closely related to tumor microenvironments (TMEs).1–3 The TME is a complex system mainly composed of blood vessels, collagen, lipid, etc.3–5 Lipid can provide energy for tumors and collagen mainly constitutes the scaffold of the tumor.6,7 Collagen and lipid content can affect the proliferation and differentiation of tumor cells.8–11 Alterations in the content of lipid and collagen are hallmarks of many diseases, including breast cancer,12 prostate cancer,13,14 renal cell carcinoma,15 skin cancers,16 and others. Therefore, the collagen and lipid content in the TME can be used for diagnosis and differentiation of tumors. Currently, several methods are available for detecting the lipid and collagen of tumors. Second-harmonic generation microscopy, which relies on the nonlinear interaction of a laser with non-centrosymmetric molecules, has been used to image fibrillar collagen within tissues.17 It can provide submicron resolution but has limited detection depth and does not provide information about lipid. Magnetic resonance imaging using chemoselective fat-suppression pulse sequences enables lipid detection. However, it fails to provide effective contrast when there is less fat within the tissue.18 In recent years, Raman spectroscopy has been successfully used to determine the molecular composition and structure based on the inelastic scattering of different molecules to light.19 However, the penetration depth of Raman spectroscopy is limited, and only superficial skin information can be obtained.19,20 To overcome these issues, new diagnostic methods with high sensitivity that can detect both collagen and lipid with sufficient penetration depths are required. Photoacoustics (PA) is a novel non-invasive detection technique that combines high optical absorption contrast with the high penetration depth of ultrasound.21 PA is a physical process of “light in and sound out.”22 A pulsed laser is used to irradiate biological tissues wherein the energy is wavelength-selectively absorbed by endogenous chromophores within the tissues, generating ultrasonic waves (PA signals) through thermoelastic expansion. Because ultrasonic waves scatter much less than optical waves, PA technology has a greater detection depth than optical detection technology, which shows great promise for clinical applications.23,24 The endogenous chromophores, such as hemoglobin, collagen, lipid, and water, have different optical absorption spectra in the visible and infrared bands.25 By irradiating biological tissues with pulsed lasers of different wavelengths, PA can provide rich endogenous chromophores information about tumors.26–29 Lei et al.30 investigated the feasibility of assessing collagen contents to detect fibrosis in Crohn’s disease using PA imaging. Wilson et al.31 implemented multiparametric spectroscopic PA imaging to assess the lipid content, oxygen saturation, and total hemoglobin to identify the development of four types of breast cancer. Conventional PA imaging is mainly based on the amplitude of the envelope of time-domain PA signals to quantify the endogenous chromophore concentrations in biological tissues,32,33 ignoring the frequency and phase information of PA signals associated with the absorbers. Besides, the envelope of the time-domain PA signal is subjected to the effects of noise and transducer response. It is quite difficult to achieve reliable results in quantifying the size and concentration of absorber with sizes smaller than the system resolution. Considering the different sizes of the absorbers, the ultrasonic spectrum shows significant advantages. PA spectral analysis (PASA) method can remove the low-frequency system noise and high-frequency measurement noise to provide objective results and repeatable measurements.34,35 Further, frequency analysis has proved feasible in detecting absorbers with sizes smaller than the system resolution.36,37 Recently, PASA, which analyzes the frequency domain power distribution of PA signals, has demonstrated the ability to assess the content and corresponding microstructure of endogenous chromophore in biological tissues simultaneously.38–41 Moreover, the spectral parameters extracted from the PA spectrum, e.g., slope, power-weighted mean frequency, can be used to characterize the microstructures of endogenous optical absorbers.38,42–44 Xu et al.39 implemented the PASA to assess the changes of lipid for fatty liver identification. They further quantified the Gleason score of prostate cancer based on the tissue microscopic architecture using PASA.45 Our group combined the PASA with machine leaning to better mine the data information and achieved a high diagnostic accuracy of prostate cancer,46–49 osteoporosis,50 and breast cancer.51,52 PASA has shown considerable potential in evaluating the endogenous chromophore in biological tissues for tumor diagnosis. Therefore, in this study, we investigated the feasibility of non-invasive PASA for characterizing the tumor-related features of lipid and collagen content in the TME to identify different types of tumors. We took skin cancers as the research objects, for which an invasive biopsy is the gold standard for diagnosis. Ex vivo experiments were conducted using three types of human skin tissue: normal, squamous cell carcinoma (SCC), and basal cell carcinoma (BCC). The content of lipid and collagen in skin tissue was calculated semi-quantitatively by PASA at different wavelengths. With the help of machine learning classification methods, tumors were successfully identified and tumor types were automatically classified based on quantified PASA parameters. 2.Materials and Methods2.1.Ethics StatementThe study protocol was approved by the Ethical Committee of the Shanghai Skin Disease Hospital and was performed in accordance with the tenets of the Declaration of Helsinki. All patients were informed of the purpose of the study, and written consent was obtained before recruitment and sampling. 2.2.Sample CollectionA total of 39 patients were enrolled in this study, including 15 with suspected SCC and 12 with suspected BCC. Normal samples were collected from the skin collection areas of 12 patients who received skin grafts. All samples were procured from the Institute of Photomedicine, Shanghai Skin Disease Hospital, School of Medicine, Tongji University, China. After surgical excision of the skin tissue, the residual blood on the tissue surface was cleaned with sterile gauze. Each sample had a diameter of . The skin tissues were placed in sterile sample tubes, stored in a portable medical cryostat (2°C to 8°C), and transported to the Institute of Acoustics of Tongji University laboratory for testing within 1 h. 2.3.Experiment Setup and PA MeasurementsFigure 1(a) shows a schematic of the PA experimental setup. An optical parametric oscillator system pumped by a Nd:YAG laser (Phocus Mobile, OPOTEK, Carlsbad, California, United States) was used to provide laser pulses with wavelengths ranging from 1200 to 1700 nm in 10 nm intervals, covering the strong absorption ranges of lipid and collagen. The laser energy output over the entire wavelength range was controlled to 0.1 to 0.5 mJ per pulse, with a pulse duration of 5 ns and pulse repetition rate of 10 Hz. A laser beam with a diameter of 3 mm illuminated the skin tissue, leading to an optical energy density of 7 to , which was below the safety limit specified by the American National Standards Institute. As shown in Fig. 1(a), skin tissue was placed on the phantom to avoid strong scattering of the sound signal by any hard boundary. The PA signals generated by the entire laser irradiation of the skin tissues were received by a needle hydrophone with a bandwidth of 1 to 20 MHz (HNC1500, ONDA Corp., Sunnyvale, California, United States). The laser energy varies with the wavelength and fluctuates only slightly over time. To determine the laser energy variations during the PA measurements, 10% of the laser energy was projected onto a black body, and the PA signals generated by the black body were received by a 5 MHz focused transducer (V326, Immersion Transducers, Olympus Corp., Tokyo, Japan). After amplification (5072PR, Olympus Corp., Tokyo, Japan) to boost the PA signals of the skin tissue samples by 25 dB, they were averaged 64 times and recorded using a digital oscilloscope (HDO6000, Teledyne Lecroy, New York, United States) at a sampling rate of 2500 MHz. To improve the stability and reduce the measurement error, PA signals from each skin tissue sample were detected at two different positions. We developed an efficient automated experimental program that includes laser wavelength switching, triggering, and data acquisition. With this program, PA data acquisition of each skin tissue sample and the blackbody could be completed in , covering 51 wavelengths (1200 to 1700 nm), with an average of 64 times per wavelength. 2.4.PA Spectral AnalysisThe PA signals were analyzed using MATLAB software (R2019b, MathWorks, Natick, Massachusetts, United States). First, the PA signal generated by the skin tissue was calibrated using the peak-to-peak value of the PA signal generated by the black body at each wavelength. Second, based on the calibrated PA signals for each skin tissue sample, the power spectra of the PA signals [Fig. 2(a)] were calculated using Welch’s method with a moving Hamming window and 60% overlap, as shown in Fig. 2(b). Considering the signal-to-noise ratio, the ultrasound frequency was first analyzed in the 1 to 8 MHz range. The power spectra of the PA signals obtained in the 1200 to 1700 nm wavelength range were combined to form a PA physiochemical spectrogram (PAPCS), as depicted in Fig. 2(c). The horizontal axis of the PAPCS is the optical wavelength, representing the relative optical absorption of different endogenous chromophores in skin tissue, whereas the vertical axis is the ultrasonic frequency, revealing the structural distribution corresponding to different optical absorptions in skin tissue. The color bar represents the amplitude of the power spectrum. Differences in the lipid and collagen content and in the microstructure of the TMEs of different tumors will form unique PAPCSs. Furthermore, the changes of relative lipid and collagen content in skin tissues were quantitatively analyzed. Based on the PA power spectral analysis at a wavelength of nm [Fig. 2(b)], we calculated the relative area of power spectrum density (APSD) at the corresponding wavelength as follows: where is the power spectral density at each frequency, , and . The PA signal is affected by low-frequency system noise and high-frequency measurement noise. Based on the PA power spectrum analysis, setting a high cut-off frequency can help avoid high-frequency noise, whereas setting a low cut-off frequency can help minimize system noise. Considering the signal-to-noise ratio and the cluster size of lipid and collagen obtained from histology images (see Fig. 4), we set the low and high cutoff frequencies to 1 and 8 MHz, respectively. The PA power spectral density in the specified frequency range was then summarized as the PA absorption value at each wavelength. A reference wavelength (690 nm) was used to eliminate systematic errors. We then obtained the relative PA absorption of each skin tissue sample with reduced system noise. To improve the stability and reduce the measurement error, PA signals from each skin tissue sample were detected at two different positions. The relative APSD values from two different directions were then obtained and averaged for further analysis. The relative APSD obtained at wavelength (nm) reflects the relative optical absorption of the corresponding endogenous chromophore; thus, it is related to the relative endogenous chromophore content of the skin tissue. According to the literature, lipid exhibits strong absorption at to 1240 nm.25,53,54 Some studies have also shown that at 1600 to 1700 nm excitation,55–57 the lipid PA signal is enhanced compared to that at 1200 nm excitation. In addition, the selective detection of collagen at to 1340 nm can be achieved.46,58 Thus, to reduce the measurement error caused by a single wavelength, the relative APSD obtained at the absorption wavelength range of lipid (or collagen) were averaged. To examine whether the changes in collagen and lipid content of TME in different types of skin tissues were statistically significant, unpaired -tests were performed using GraphPad Prism 9.0.2.5.Support Vector Machine AnalysisSupport vector machine (SVM) analysis is one of the simplest machine learning classification method that is supported by rigorous mathematical theory, is highly interpretable, and can identify the key factors for classification tasks. In this study, a SVM classifier was applied to perform automatic different types of tumors classification by combining relative APSD values at different wavelengths. The APSD values obtained in the three wave bands (1200 to 1240 nm, 1600 to 1700 nm, and 1300 to 1340 nm) were used as the input characteristic parameters for the classification of different types of skin tissues, thus realizing their automatic discrimination, as shown in Fig. 3. The input data comprised 39 datasets with 21 features. Considering the limited number of clinical skin tissue samples available for categorization, we used the C-type SVM model with a radial basis function (RBF) kernel, which is suitable for classifying small amounts of data. The RBF kernel can be defined as follows: where is the variable matrix, and are the indices of the matrix, and is the kernel parameter. There are two parameters for an RBF kernel: the regularization parameter and kernel parameter . Parameter selection for optimal SVM categorization was achieved by grid searching59 using cross-validation. Various pairs of values were tested, and the pair with the best cross-validation accuracy was selected. To increase the classification reliability, a three-fold cross-validation approach was applied by rotating the datasets used to train and test the SVM model. The samples were randomly divided into two groups: two-thirds of the samples were used for training, and the remaining samples were used for prediction. This process was implemented thrice, with each group being tested once. Finally, the performance of the classifier model was evaluated using confusion matrices and receiver operating characteristics (ROCs) curves.2.6.HistopathologyTo validate the PA measurements, the samples were collected and prepared for histological analysis after performing the PA measurements. The pathological analysis was divided into two parts. Some tissue samples were fixed by immersion in a 10% formalin solution and subsequent sectioning after tissue dehydration and paraffin embedding, completing the production of paraffin sections. Next, the sections were stained with hematoxylin and eosin (H&E) and others stained with Masson’s trichrome. Masson’s trichrome can stain collagen with blue color, which can be used for collagen detection. The other tissue sample was treated with an optical coherence tomography embedding agent (optical cutting temperature compound, Sakura Americas, New Mexico, United States) and flushed with distilled water to obtain frozen sections. Serial longitudinal sections ( thick) were subsequently cut and stained with Nile red for lipid detection. Nile red can stain the lipid with yellow color. The staining results were observed by microscopy, as shown in Fig. 4. The relative lipid and collagen contents of the skin tissues were obtained from histological images. The ratio of the positively stained area to the whole area in the histological image, calculated using ImageJ software for each sample slice, was used as the gold standard for the relative content of collagen or lipid in the skin tissues, presenting a quantitative comparison with the relative APSD results. 3.Results3.1.PAPCSs of Different Types of Skin TissueFigures 5(a)–5(c) show the averaged PAPCSs of the normal, SCC, and BCC tissues, respectively. Compared with the normal tissue, in the 1200 to 1240 nm and 1600 to 1700 nm wave bands, the PAPCSs of tumor (SCC and BCC) tissues both show lower spectral magnitudes, indicating decreases in lipid content. In addition to the changes in lipid content, the PAPCSs of tumors exhibit differences in the collagen absorption wavelength range (1300 to 1340 nm). The PAPCSs of SCC and BCC show lower PA spectral magnitudes, indicating reductions in collagen contents. Interestingly, compared with the PAPCS of the BCC tissues, that of the SCC tissues shows higher spectral magnitudes in the lipid- and collagen-absorption wavelength ranges, implying that SCC tissues contain higher lipid and collagen content than BCC tissues. Specifically, Figs. 5(e)–5(g) show the PA power spectrum of different skin tissues at the absorption peaks of lipid (1210 and 1700 nm) and collagen (1310 nm), respectively. The pathological images of the skin tissues [Figs. 4(d)–4(i)] verify these differences in the lipid and collagen contents of the three types of skin tissues. The differences in PAPCSs between different skin tumors suggest that lipid and collagen in TMEs can be used as biomarkers of tumor diversity. 3.2.PASA Semi-Quantitative ResultsBased on the PASA results for skin tissues obtained at different wavelengths, the relative APSD was calculated to semi-quantify the relative lipid and collagen contents in the skin tissues. We calculated the wavelength dependence of the relative APSD for each skin tissue sample. Figures 6(a)–6(c) show the content characterization fluctuations at different wavelength. Thus, we averaged the relative APSD values in the absorption band of lipid or collagen to eliminate the measurement error caused by a single wavelength. Figures 6(d)–6(f) show the statistical results for the averaged relative APSD values, which can be used to distinguish between different types of skin tissues. As shown in Figs. 6(d) and 6(e), the lower lipid content of tumors leads to a decrease in its spectral amplitude and, consequently, a lower relative APSD for the 1200 to 1240 nm and 1600 to 1700 nm wavebands. The relative APSD values of the tumors are also lower for the 1300 to 1340 nm wavelength band, corresponding to their decreased collagen content, as shown in Fig. 6(f). In addition, different types of tumors exhibit significant differences in their collagen and lipid contents (). For example, as shown in Figs. 6(d)–6(f), compared to those of SCC tissues, the relative APSD values of BCC tissues are lower in the regions corresponding to the optical absorption of lipid (1200 to 1240 nm and 1600 to 1700 nm) and collagen (1300 to 1340 nm), indicating that the lipid and collagen contents in the BCC tissues are lower than those contents in SCC. Figures 4(j) and 4(k) depict the relative lipid and collagen contents of the three different skin tissues obtained from the histological images and reveal notable differences. The PA results were corroborated by the changes in the lipid and collagen contents revealed in the histological images. 3.3.SVM Classification ResultsAn SVM was employed to classify skin cancers based on the lipid and collagen contents in the skin tissues. The relative lipid and collagen contents in the TME were characterized by the relative APSD values at different wavelengths. The classification results of the three-fold cross-validation approach are presented in Table 1. As shown in Table 2, the average classification accuracies of the SVM trained based only on the parameter at a single wavelength are 69.2%, 64%, and 89%. However, combining the parameters of the three bands results in an accuracy of 92.3%, the results show that the combination of multiple biomarkers can achieve more accurate intelligent diagnoses. Table 1SVM categorization. A three-fold cross-validation approach was conducted in which the samples were divided into three groups at random. Each cycle used two groups for training and one group for testing. Cate. = Categorization. Samp. = Sample.
Table 2Mean SVM categorization accuracy in Table 1.
The classification performance of the SVM was evaluated using a confusion matrix, as shown in Fig. 7(b). The confusion matrix provides the classification accuracy of the SVM for individual classes. Each column of the matrix corresponds to a true label, and each row corresponds to a predicted label. The main diagonal shows the classification accuracy of the SVM for the individual classes. The off-diagonal values indicate the misclassification rates of the SVM for classifying the individual classes. The results show that SVM-based multiclass categorization can achieve diagnostic accuracies of 91.7% (normal), 93.3% (SCC), and 91.7% (BCC). The ROC curves yielded by the SVM for the normal, SCC, and BCC tissues are presented in Fig. 7(c), revealing areas under the ROC curves of 0.94, 0.96, and 0.92, respectively. Overall, SVM-based categorization proved effective in diagnosing tumors with a high level of accuracy. 4.DiscussionTME is a complex system in which collagen and lipid are the important components of TME. It has been reported that collagen and lipid can be used as biomarker to assess tumors.61,62 Several techniques have been available to detect lipid and collagen, but they all have limitations, including: invasiveness, time-consumption, and limited detection depths. Based on the contrasting optical absorption spectra of different endogenous chromophores in tissue, PA detection technique can provide rich content and microstructure information of the endogenous chromophores in biological tissues, allowing more accurate tumor diagnoses. In addition, PA has a greater detection depth than optical detection technology, benefiting from the fact that the PA signal originates from the absorption of laser energy by biological tissue, independent of the phase of the light wave. Experiments on different types of skin tissues showed that the PAPCS obtained by PASA contained rich diagnostic information. The loss of collagen and lipid content can be observed in the corresponding PAPCS and can be used to characterize different types of tumors. Owing to the differences in lipid and collagen content of different tumors, each type of tumor has a unique PAPCS. The PASA parameter APSD is correlated with the collagen and lipid contents. The decreases in the lipid and collagen contents in skin cancer tissues cause the PA spectral amplitude of the power spectrum to decrease, corresponding to a lower relative APSD value. Furthermore, SCC originates from higher differentiated epithelial cells than BCC, with higher amounts of lipid and collagen,63 resulting in their APSD values higher than those for BCC in the absorption band of lipid and collagen. Statistical results of the content and microstructure changes of endogenous chromophores related to histology can be obtained according to PASA parameters. Thus, PASA provides an objective method of classifying SCC and BCC independently of the experience of a physician or pathologist. In addition, the PASA eliminates the effects of noise and system errors, providing system-independent semi-quantitative results. Our results suggest that PASA can provide sensitive, semi-quantitative content of endogenous chromophore in tissues to achieve non-invasive diagnosis and identification of different types of tumors. Frequency spectrum of PA signal is related to the size and content of the optical absorber. In addition, our previous work has investigated the frequency anisotropy caused by the direction of the structure35. There appeared well distributed and more parallel collagen in normal skin tissues. While in the tumor tissue (SCC and BCC), the collagen was obviously unorderly and the amount of collagen appears reduced.64 There was no statistical difference in the structure and distribution of collagen cluster between SCC and BCC. In the future, we can use multiple parameters to characterize the distribution and content of collagen and lipid simultaneously to better identify tumor tissues from normal tissues. Machine learning is an effective feature extraction tool. Even in small samples, the characteristics of PA signals can be fully learned to achieve good disease diagnosis, and machine learning has been widely used in the detection of various diseases.65,66 In this study, an SVM, one of the simplest machine learning tools, was utilized to distinguish automatically between different types of tumors based on the extracted PA parameters. The proposed method takes advantage of quantitative parameters of multiple wavelengths, contains rich diagnostic information, and improves the separation of high-dimensional datasets. As the number of clinical skin tissue samples was limited, we used the SVM algorithm with the RBF kernel, which is appropriate when the amount of data to be classified is small. Furthermore, to prevent overfitting based on a small number of training datasets, -fold cross-validation was performed. The choice of the number of -fold depends on the sample size, number of parameters, structure of data, and so on.67 Due to the limitation of the sample size, three-fold cross-validation was used in our study. In future work, more training datasets should be included and a higher and appropriate -fold cross-validation should be selected, which would improve the robustness of the SVM model and make the classification of skin cancer types using the SVM approach more accurate. Considering the accessibility of human tumor samples, the object of this study is human skin tumor tissue in vitro for pathological examination. In vitro experiments significantly improved the operability of the samples and the stability of the experiments. However, the disadvantage is that the optical absorption of hemoglobin cannot be measured. The optical absorption of hemoglobin is dominant in the visible to near infrared wavelength range (i.e., 400 to 900 nm). Therefore, the results acquired at the wavelengths of 1200 to 1700 nm were less affected by the hemoglobin. In the future, we plan to use the PASA method in vivo. Considering the influence of water in vivo experiment, which has high absorption in the 1200 to 1700 nm, a spectral unmixing procedure should be used to enable more precise evaluation of the relative contents of lipid and collagen. In fact, this approach also has considerable potential to be extended to the diagnosis of other types of cancers with altered collagen or lipid contents, e.g., breast cancer,12 prostate cancer.13,14. Furthermore, there were clear differences in the PA power spectrum between skin cancers and normal tissues. This finding not only confirms that the PA signals of tumor and healthy tissue are different but also suggests the possibility of using the PASA to define tumor boundaries in the future. Currently, the primary treatment for skin cancers is surgical excision. However, the existing clinical techniques cannot accurately define the boundaries of tumors, leading to a postoperative residual tumor tissue, which may cause disease recurrence. In future studies, we also intend to apply the PASA in vivo to identify the boundaries of skin cancers more accurately and objectively, to help doctors remove such tumors completely and reduce complications. 5.ConclusionIn summary, we introduced the PASA method to characterize the content of endogenous chromophore in the TME to better understand the interaction and regulation of tumors and TME. The parameters extracted from PASA were used to semi-quantify the content of endogenous chromophore. It was found that the collagen and lipid content can be used as biomarkers of tumor diversity, and we used this to successfully diagnose different types of tumors with improved accuracy. Considering the PA technique is non-invasive and has greater detection depth than pure optical technology, the proposed method shows considerable potential for non-invasive and more accurate diagnosis of tumors in vivo. Code, Data, and Materials AvailabilityThe data that support the findings of this study are available from the corresponding author upon reasonable request. AcknowledgmentsThis project was supported by the National Natural Science Foundation of China (Grant No. 12034015, 62088101, 81872212, and 82073013), National Key Research and Development Program of China (Grant No. 2017YFC0111400), Shanghai Municipal Science and Technology Major Project (Grant No. 2021SHZDZX0100), Shanghai Municipal Commission of Science and Technology Project (Grant No. 19511132101), Program of Shanghai Academic Research Leader (Grant No. 21XD1403600). We would like to thank Editage68 for English language editing. ReferencesZ. Tan et al.,
“The role of tumor inflammatory microenvironment in lung cancer,”
Front. Pharmacol., 12 688625 https://doi.org/10.3389/fphar.2021.688625
(2021).
Google Scholar
X. Huang et al.,
“Regulation of tumor microenvironment for pancreatic cancer therapy,”
Biomaterials, 270 120680 https://doi.org/10.1016/j.biomaterials.2021.120680 BIMADU 0142-9612
(2021).
Google Scholar
N. M. Anderson and M. C. Simon,
“The tumor microenvironment,”
Curr. Biol., 30
(16), R921
–R925 https://doi.org/10.1016/j.cub.2020.06.081 CUBLE2 0960-9822
(2020).
Google Scholar
F. R. Balkwill, M. Capasso and T. Hagemann,
“The tumor microenvironment at a glance,”
J. Cell Sci., 125
(23), 5591
–5596 https://doi.org/10.1242/jcs.116392 JNCSAI 0021-9533
(2012).
Google Scholar
E. Brown, J. Brunker and S. E. Bohndiek,
“Photoacoustic imaging as a tool to probe the tumour microenvironment,”
Dis. Model. Mech., 12
(7), dmm039636 https://doi.org/10.1242/dmm.039636
(2019).
Google Scholar
J. Sapudom and T. Pompe,
“Biomimetic tumor microenvironments based on collagen matrices,”
Biomater. Sci., 6
(8), 2009
–2024 https://doi.org/10.1039/C8BM00303C
(2018).
Google Scholar
M. Fang et al.,
“Collagen as a double-edged sword in tumor progression,”
Tumor Biol., 35
(4), 2871
–2882 https://doi.org/10.1007/s13277-013-1511-7
(2014).
Google Scholar
P. Cirri and P. Chiarugi,
“Cancer associated fibroblasts: the dark side of the coin,”
Am. J. Cancer Res., 1
(4), 482
–497
(2011).
Google Scholar
K. R. Levental et al.,
“Matrix crosslinking forces tumor progression by enhancing integrin signaling,”
Cell, 139
(5), 891
–906 https://doi.org/10.1016/j.cell.2009.10.027 CELLB5 0092-8674
(2009).
Google Scholar
F. Mollinedo and C. Gajate,
“Lipid rafts as signaling hubs in cancer cell survival/death and invasion: implications in tumor progression and therapy,”
J. Lipid Res., 61
(5), 611
–635 https://doi.org/10.1194/jlr.TR119000439 JLPRAW 0022-2275
(2020).
Google Scholar
M. Martin-Perez et al.,
“The role of lipids in cancer progression and metastasis,”
Cell Metab., 34
(11), 1675
–1699 https://doi.org/10.1016/j.cmet.2022.09.023 1550-4131
(2022).
Google Scholar
H. Abramczyk et al.,
“The role of lipid droplets and adipocytes in cancer. Raman imaging of cell cultures: MCF10A, MCF7, and MDA-MB-231 compared to adipocytes in cancerous human breast tissue,”
Analyst, 140
(7), 2224
–2235 https://doi.org/10.1039/C4AN01875C ANLYAG 0365-4885
(2015).
Google Scholar
Y. Pu et al.,
“Changes of collagen and nicotinamide adenine dinucleotide in human cancerous and normal prostate tissues studied using native fluorescence spectroscopy with selective excitation wavelength,”
J. Biomed. Opt., 15
(4), 047008 https://doi.org/10.1117/1.3463479 JBOPFO 1083-3668
(2010).
Google Scholar
N. Burns-Cox et al.,
“Change in collagen metabolism in prostate cancer: a host response that may alter progression,”
J. Urol., 166
(5), 1698
–1701 https://doi.org/10.1016/S0022-5347(05)65656-X JOURAA 0022-5347
(2001).
Google Scholar
S. L. Best et al.,
“Collagen organization of renal cell carcinoma differs between low and high grade tumors,”
BMC Cancer, 19
(1), 490 https://doi.org/10.1186/s12885-019-5708-z BCMACL 1471-2407
(2019).
Google Scholar
L. A. Austin, S. Osseiran and C. L. Evans,
“Raman technologies in cancer diagnostics,”
Analyst, 141
(2), 476
–503 https://doi.org/10.1039/C5AN01786F ANLYAG 0365-4885
(2016).
Google Scholar
A. Zoumi, A. Yeh and B. J. Tromberg,
“Imaging cells and extracellular matrix in vivo by using second-harmonic generation and two-photon excited fluorescence,”
Proc. Natl. Acad. Sci. U.S.A., 99
(17), 11014
–11019 https://doi.org/10.1073/pnas.172368799 PNASA6 0027-8424
(2002).
Google Scholar
S. S. Pokharel et al.,
“Current MR imaging lipid detection techniques for diagnosis of lesions in the abdomen and pelvis,”
Radiographics, 33
(3), 681
–702 https://doi.org/10.1148/rg.333125068
(2013).
Google Scholar
C. A. Lieber et al.,
“In vivo nonmelanoma skin cancer diagnosis using Raman microspectroscopy,”
Lasers Surg. Med., 40
(7), 461
–467 https://doi.org/10.1002/lsm.20653 LSMEDI 0196-8092
(2008).
Google Scholar
A. M. Ferreira Lima et al.,
“Discrimination of non-melanoma skin cancer and keratosis from normal skin tissue in vivo and ex vivo by Raman spectroscopy,”
Vib. Spectrosc., 100 131
–141 https://doi.org/10.1016/j.vibspec.2018.11.009 VISPEK 0924-2031
(2019).
Google Scholar
L. V. Wang and J. Yao,
“A practical guide to photoacoustic tomography in the life sciences,”
Nat. Methods, 13
(8), 627
–638 https://doi.org/10.1038/nmeth.3925
(2016).
Google Scholar
S. Zackrisson, S. M. W. Y. van de Ven and S. S. Gambhir,
“Light in and sound out: emerging translational strategies for photoacoustic imaging,”
Cancer Res., 74
(4), 979
–1004 https://doi.org/10.1158/0008-5472.CAN-13-2387 CNREA8 0008-5472
(2014).
Google Scholar
L. V. Wang and S. Hu,
“Photoacoustic tomography: in vivo imaging from organelles to organs,”
Science, 335
(6075), 1458
–1462 https://doi.org/10.1126/science.1216210 SCIEAS 0036-8075
(2012).
Google Scholar
V. Ntziachristos et al.,
“Looking and listening to light: the evolution of whole-body photonic imaging,”
Nat. Biotechnol., 23
(3), 313
–320 https://doi.org/10.1038/nbt1074 NABIF9 1087-0156
(2005).
Google Scholar
H.-W. Wang et al.,
“Label-free bond-selective imaging by listening to vibrationally excited molecules,”
Phys. Rev. Lett., 106
(23), 238106 https://doi.org/10.1103/PhysRevLett.106.238106 PRLTAO 0031-9007
(2011).
Google Scholar
C. Kim, C. Favazza and L. V. Wang,
“In vivo photoacoustic tomography of chemicals: high-resolution functional and molecular optical imaging at new depths,”
Chem. Rev., 110
(5), 2756
–2782 https://doi.org/10.1021/cr900266s CHREAY 0009-2665
(2010).
Google Scholar
S. Hu and L. V. Wang,
“Optical-resolution photoacoustic microscopy: auscultation of biological systems at the cellular level,”
Biophys. J., 105
(4), 841
–847 https://doi.org/10.1016/j.bpj.2013.07.017 BIOJAU 0006-3495
(2013).
Google Scholar
J. Yao et al.,
“High-speed label-free functional photoacoustic microscopy of mouse brain in action,”
Nat. Methods, 12
(5), 407
–410 https://doi.org/10.1038/nmeth.3336 1548-7091
(2015).
Google Scholar
LED-Based Photoacoustic Imaging: From Bench to Bedside, Springer Singapore, Singapore
(2020). Google Scholar
H. Lei et al.,
“Characterizing intestinal inflammation and fibrosis in Crohn’s disease by photoacoustic imaging: feasibility study,”
Biomed. Opt. Express, 7
(7), 2837 https://doi.org/10.1364/BOE.7.002837 BOEICL 2156-7085
(2016).
Google Scholar
K. E. Wilson et al.,
“Multiparametric spectroscopic photoacoustic imaging of breast cancer development in a transgenic mouse model,”
Theranostics, 4
(11), 1062
–1071 https://doi.org/10.7150/thno.9922
(2014).
Google Scholar
H. F. Zhang et al.,
“Functional photoacoustic microscopy for high-resolution and noninvasive in vivo imaging,”
Nat. Biotechnol., 24
(7), 848
–851 https://doi.org/10.1038/nbt1220 NABIF9 1087-0156
(2006).
Google Scholar
X. Wang et al.,
“Noninvasive laser-induced photoacoustic tomography for structural and functional in vivo imaging of the brain,”
Nat. Biotechnol., 21
(7), 803
–806 https://doi.org/10.1038/nbt839 NABIF9 1087-0156
(2003).
Google Scholar
T. Feng et al.,
“Characterization of bone microstructure using photoacoustic spectrum analysis,”
Opt. Express, 23
(19), 25217 https://doi.org/10.1364/OE.23.025217 OPEXFF 1094-4087
(2015).
Google Scholar
M. Zhang et al.,
“Photoacoustic power Azimuth spectrum for microvascular evaluation,”
Photoacoustics, 22 100260 https://doi.org/10.1016/j.pacs.2021.100260
(2021).
Google Scholar
S. Wang et al.,
“Quantitative detection of stochastic microstructure in turbid media by photoacoustic spectral matching,”
Appl. Phys. Lett., 102
(11), 114102 https://doi.org/10.1063/1.4798272 APPLAB 0003-6951
(2013).
Google Scholar
X. Gao et al.,
“Quantitative imaging of microvasculature in deep tissue with a spectrum-based photo-acoustic microscopy,”
Opt. Lett., 40
(6), 970 https://doi.org/10.1364/OL.40.000970 OPLEDP 0146-9592
(2015).
Google Scholar
G. Xu et al.,
“Photoacoustic spectrum analysis for microstructure characterization in biological tissue: a feasibility study,”
Appl. Phys. Lett., 101
(22), 221102 https://doi.org/10.1063/1.4768703 APPLAB 0003-6951
(2012).
Google Scholar
G. Xu et al.,
“The functional pitch of an organ: quantification of tissue texture with photoacoustic spectrum analysis,”
Radiology, 271
(1), 248
–254 https://doi.org/10.1148/radiol.13130777 RADLAX 0033-8419
(2014).
Google Scholar
T. Feng et al.,
“Detection of collagen by multi-wavelength photoacoustic analysis as a biomarker for bone health assessment,”
Photoacoustics, 24 100296 https://doi.org/10.1016/j.pacs.2021.100296
(2021).
Google Scholar
K. Ma et al.,
“Myocardial infarct border demarcation by dual-wavelength photoacoustic spectral analysis,”
Photoacoustics, 26 100344 https://doi.org/10.1016/j.pacs.2022.100344
(2022).
Google Scholar
W. Xie et al.,
“Wavelet transform-based photoacoustic time-frequency spectral analysis for bone assessment,”
Photoacoustics, 22 100259 https://doi.org/10.1016/j.pacs.2021.100259
(2021).
Google Scholar
W. Xie et al.,
“Photoacoustic characterization of bone physico-chemical information,”
Biomed. Opt. Express, 13
(5), 2668 https://doi.org/10.1364/BOE.457278 BOEICL 2156-7085
(2022).
Google Scholar
Y. Yang et al.,
“Photoacoustic tomography of tissue subwavelength microstructure with a narrowband and low frequency system,”
Appl. Phys. Lett., 101
(3), 034105 https://doi.org/10.1063/1.4736994 APPLAB 0003-6951
(2012).
Google Scholar
G. Xu et al.,
“Quantifying Gleason scores with photoacoustic spectral analysis: feasibility study with human tissues,”
Biomed. Opt. Express, 6
(12), 4781 https://doi.org/10.1364/BOE.6.004781 BOEICL 2156-7085
(2015).
Google Scholar
S. Huang et al.,
“Interstitial assessment of aggressive prostate cancer by physio-chemical photoacoustics: an ex vivo study with intact human prostates,”
Med. Phys., 45
(9), 4125
–4132 https://doi.org/10.1002/mp.13061 MPHYA6 0094-2405
(2018).
Google Scholar
Y. Chen et al.,
“Prostate cancer identification via photoacoustic spectroscopy and machine learning,”
Photoacoustics, 23 100280 https://doi.org/10.1016/j.pacs.2021.100280
(2021).
Google Scholar
S. Wu et al.,
“Photoacoustic spectrum analysis for quick identification and grading of prostate cancer,”
in IEEE Int. Ultrasonics Symp.,
1
–4
(2020). https://doi.org/10.1109/IUS46767.2020.9251610 Google Scholar
S. Wu et al.,
“Quick identification of prostate cancer by wavelet transform-based photoacoustic power spectrum analysis,”
Photoacoustics, 25 100327 https://doi.org/10.1016/j.pacs.2021.100327
(2022).
Google Scholar
T. Feng et al.,
“Characterization of multi-biomarkers for bone health assessment based on photoacoustic physicochemical analysis method,”
Photoacoustics, 25 100320 https://doi.org/10.1016/j.pacs.2021.100320
(2022).
Google Scholar
J. Li et al.,
“Photoacoustics spectral analysis for in vivo detection of collagen contents in cancers,”
in IEEE Int. Ultrasonics Symp.,
1
–4
(2022). https://doi.org/10.1109/IUS54386.2022.9958255 Google Scholar
J. Li et al.,
“Molecular breast cancer subtype identification using photoacoustic spectral analysis and machine learning at the biomacromolecular level,”
Photoacoustics, 30 100483 https://doi.org/10.1016/j.pacs.2023.100483
(2023).
Google Scholar
B. Wang et al.,
“Detection of lipid in atherosclerotic vessels using ultrasound-guided spectroscopic intravascular photoacoustic imaging,”
Opt. Express, 18
(5), 4889 https://doi.org/10.1364/OE.18.004889 OPEXFF 1094-4087
(2010).
Google Scholar
J. Hui et al.,
“Bond-selective photoacoustic imaging by converting molecular vibration into acoustic waves,”
Photoacoustics, 4
(1), 11
–21 https://doi.org/10.1016/j.pacs.2016.01.002
(2016).
Google Scholar
P. Wang et al.,
“Bond-selective imaging of deep tissue through the optical window between 1600 and 1850 nm,”
J. Biophotonics, 5
(1), 25
–32 https://doi.org/10.1002/jbio.201100102
(2012).
Google Scholar
B. Wang et al.,
“Intravascular photoacoustic imaging of lipid in atherosclerotic plaques in the presence of luminal blood,”
Opt. Lett., 37
(7), 1244 https://doi.org/10.1364/OL.37.001244 OPLEDP 0146-9592
(2012).
Google Scholar
C. K. N. Patel and A. C. Tam,
“Optical absorption coefficients of water,”
Nature, 280
(5720), 302
–304 https://doi.org/10.1038/280302a0
(1979).
Google Scholar
Y. Zhu et al.,
“Identifying intestinal fibrosis and inflammation by spectroscopic photoacoustic imaging: an animal study in vivo,”
Biomed. Opt. Express, 9
(4), 1590 https://doi.org/10.1364/BOE.9.001590 BOEICL 2156-7085
(2018).
Google Scholar
C.-W. Hsu, C.-C. Chang and C.-J. Lin,
“A practical guide to support vector classification,”
www.csie.ntu.edu.tw/cjlin/papers/guide/guide
(2003).
Google Scholar
W. F. Cheong, S. A. Prahl and A. J. Welch,
“A review of the optical properties of biological tissues,”
IEEE J. Quantum Electron., 26
(12), 2166
–2185 https://doi.org/10.1109/3.64354 IEJQA7 0018-9197
(1990).
Google Scholar
S. A. Arnold et al.,
“Lack of host SPARC enhances vascular function and tumor spread in an orthotopic murine model of pancreatic carcinoma,”
Dis. Model. Mech., 3
(1–2), 57
–72 https://doi.org/10.1242/dmm.003228
(2010).
Google Scholar
X. Li et al.,
“The immunological and metabolic landscape in primary and metastatic liver cancer,”
Nat. Rev. Cancer, 21
(9), 541
–557 https://doi.org/10.1038/s41568-021-00383-9 NRCAC4 1474-175X
(2021).
Google Scholar
S. Heuke et al.,
“Detection and discrimination of non-melanoma skin cancer by multimodal imaging,”
Healthcare, 1
(1), 64
–83 https://doi.org/10.3390/healthcare1010064
(2013).
Google Scholar
S. Y. Xiong, J. G. Yang and J. Zhuang,
“Nonlinear spectral imaging of human normal skin, basal cell carcinoma and squamous cell carcinoma based on two-photon excited fluorescence and second-harmonic generation,”
Laser Phys., 21
(10), 1844
–1849 https://doi.org/10.1134/S1054660X11170269 LAPHEJ 1054-660X
(2011).
Google Scholar
J. Rodrigues,
“Exploring photoacoustic spectroscopy-based machine learning together with metabolomics to assess breast tumor progression in a xenograft model ex vivo,”
Lab. Invest., 101
(7), 952
–965 https://doi.org/10.1038/s41374-021-00597-3 LAINAW 0023-6837
(2021).
Google Scholar
A. R. Rajanna et al.,
“Prostate cancer detection using photoacoustic imaging and deep learning,”
Electron. Imaging, 2016
(15), 1
–6 https://doi.org/10.2352/ISSN.2470-1173.2016.15.IPAS-189
(2016).
Google Scholar
Y. Jung,
“Multiple predicting -fold cross-validation for model selection,”
J. Nonparametr. Stat., 30
(1), 197
–215 https://doi.org/10.1080/10485252.2017.1404598 NOSTEK 1048-5252
(2018).
Google Scholar
BiographyMengjiao Zhang received her BS degree in marine technology from Ocean University of China, in 2013. She is currently pursuing a PhD in physics at the Institute of Acoustics, School of Physics Science and Engineering, in Tongji University. From 2021 to 2022, she was a visiting student in the School of Biomedical Engineering and Imaging Sciences, King’s College London. Her current research focuses on photoacoustic (PA) measurement of microvascular evaluation and PA imaging. Long Wen received a master’s degree in dermatology and venereology from Anhui Medical University, China, in 2019. Currently, he is a junior physician fellow at Shanghai Skin Disease Hospital affiliated with Tongji University. His research interests include the non-invasive diagnosis of skin diseases, especially skin tumors, and photodynamic therapy of skin diseases. Chu Zhou received his master’s degree in dermatology and venereology from Anhui Medical University, Hefei, China. His research focuses on optical coherence tomography and PA imaging of skin diseases. Jing Pan received her master’s degree in physics from Tongji University, China, in 2019. Her research interest includes the medical imaging and spectrum analysis, especially PA imaging and analysis of skin tumors and prostate cancer. Shiying Wu is a PhD candidate student at Institute of Acoustics, School of Physics Science and Engineering of Tongji University. She received her bachelor’s degree from Tongji University. During her undergraduate study, she won the scholarship for outstanding students and the second prize of mathematical modeling contest. During her postgraduate study, she won the honor of doctoral student scholarship. Her current research focuses on PA spectrum analysis and PA imaging. Peiru Wang is a dermatologist in Shanghai Skin Disease Hospital and an associate professor of Tongji University. She is engaged clinical transformation of photomedical technology and treatment and perform research on skin aging and skin cancer. She has published 32 peer-reviewed articles in English and 8 articles in Chinese about dermatology and ALA-PDT. She is now involved in 8 academic projects including the Natural Science Foundation program (NSFC) and several clinical trials regarding skin diseases and cosmetics. Haonan Zhang received his PhD in physics from Tongji University, Shanghai, China, in 2019. His research focuses on PAs imaging and ultrasonic therapy. Panpan Chen received her BS degree in photoelectric information science and engineering from Anhui Normal University, China, in 2018. She is currently pursuing her MS degree at the Institute of Acoustics, Tongji University, China. Her current research interests include PA image reconstruction, bone imaging and deep learning. Her most recent research focused on the PA imaging of numerical bone model based on the modified deep learning network. Qi Chen received his bachelor’s degree in medical imaging from Wenzhou Medical University, Wenzhou, China. Currently, he is a graduate student at Tongji University School of Medicine, Shanghai, China. His research focuses on the application of PA imaging in dermatology. Xiuli Wang is a doctoral tutor, the director of the Institute of Photomedicine, Tongji University School of Medicine, the chief expert of the Photodynamic Therapy Research Center of the Dermatology Branch of the Chinese Medical Association, and the executive director of the International Photodynamic Medicine Association. She was the first to carry out photodynamic clinical treatment, basic research, and technology promotion in dermatology in China. She took the lead in establishing the China Photodynamic Clinical Treatment Center and formulated a Chinese plan for the photodynamic treatment of skin diseases. Qian Cheng received her BS degree in physics, and MS and PhD degrees in acoustics from Tongji University, China, in 2000, 2003, and 2006, respectively. She is currently a professor at Tongji University. Since 2006, her research interests have been optoacoustic phenomena, near-field acoustic imaging technique, Schlieren imaging technique, and the development of the acoustic detecting instruments. Her most recent research has focused on the clinical translation of PA imaging and quantitative analysis, and in particular for tumor diagnosis and evaluation. |
CITATIONS
Cited by 2 scholarly publications.
Tissues
Tumors
Collagen
Skin
Biological samples
Photoacoustic spectroscopy
Absorption