ORIGINAL RESEARCH article
Bone Chemical Composition Analysis Using Photoacoustic Technique
- 1School of Electronic and Optical Engineering, Nanjing University of Science and Technology, Nanjing, China
- 2Institute of Acoustics, School of Physics Science and Engineering, Tongji University, Shanghai, China
- 3Department of Electronic Engineering, Fudan University, Shanghai, China
- 4Key Laboratory of Spine and Spinal Cord Injury Repair and Regeneration of Ministry of Education, Orthopaedic Department of Tongji Hospital, Tongji University School of Medicine, Shanghai, China
Photoacoustic (PA) signal analysis based on ultrasonic wave detection can provide both high-sensitivity optical contrast information and micro-architectural information which is highly related with the chemical composition of tissue. In this study, the feasibility assessment of bone composition assessment was investigated using the multi-wavelength PA analysis (MWPA) method which could reflect the molecular information. By illuminating a bone specimen using a laser light with wavelength over an optical spectrum ranging from 680 to 950 nm, the optical absorption spectrum of the bone was acquired. Then, with the optical absorption spectra of all optical absorption chemical components in the known bone, a spectral unmixing procedure was performed to quantitatively assess the relative content of each chemical component. The experimental results from rabbit bones show that MWPA method can be used to assess chemical components related to bone metabolism. Our study confirmed that PA technique can be used as a novel bone diagnostic technique by providing new information about the quantity of bone and identifying biomarkers of bone that can improve the current diagnostic imaging techniques.
Osteoporosis, a serious public health threat with significant physical, psychological and economic impacts, is expected to increase in association with worldwide aging of the population. In osteoporosis, the bone mineral density (BMD) decreases, bone microarchitecture (BMA) deteriorates, and the amount and type of proteins in bone alter . Currently, most clinically used non-invasive assessment methods are based on the use of X-ray or ultrasound [2, 3]. These methods, in spite of the applicability to measure bone mineral density (BMD) as well as some mechanical properties, have limited sensitivity to monitor the chemical or molecular changes in the bone. In addition, X-ray based techniques use ionizing radiation, which is not ideal for pediatric, or long-term repetitive monitoring. Quantitative ultrasound (QUS) technology as a practical, low-cost alternative has already led to clinical instrumentation [4, 5]. The QUS bone assessment method is primarily based on the measurement of sound velocity (SOS) and broadband ultrasound attenuation (BUA) through a given tissue. However, the specificity of QUS is limited when pathogenic bone diseases are determined by microstructure and chemical changes [6–8]. Bone quantity and quality are dependent on not only the mass and structure of non-organic mineral matrix but also the organic matrix which is associated with the bone blood flow and cellular metabolism. Recently, it has been reported that magnetic resonance imaging (MRI) can distinguish changes in bone marrow lipid content and bone microarchitecture between normal bone and osteopenia bone [9–11]. Due to the high cost and complexity, MRI examinations are impossible to replace standard DEXA (Dual-Energy X-ray Absorptiometry) measurements with more advanced MRI analysis. In prior studies, optical spectroscopic techniques have been used to evaluate how the alterations of bone composition contribute to bone quality changes related to aging, disease, or injury [12–18]. However, traditional optical techniques suffer from limited spatial resolution and the overwhelming optical scattering in biological tissues, thus reducing the efficacy of skeletal imaging in vivo.
The emerging biomedical photoacoustic (PA) techniques have a unique ability to probe the highly sensitive optical absorption contrast in deep biological tissues [19–23]. The PA signal generated by the bone contains both the microstructural information and molecular information, which are both highly correlated with bone health. Lashkari et al. evaluated the cortical and trabecular bone structure and density variations by using a dual backscattered ultrasound and PA radar system [24, 25]. Furthermore, our group has studied the feasibility of accessing BMD and BMA of the trabecular bone in rat models through using thermal photoacoustic (TPA) and photoacoustic spectral analysis (PASA). Recently, Idan Steinberg et al. used the dual-modality multispectral photoacoustic system to quantify the blood/fat ratio present in the marrow, which has been correlated with molecular changes in the long bone.
In this study, the feasibility of the multi-wavelength PA analysis (MWPA) technique in quantifying the molecular information of trabecular bone based on rabbit models with different BMDs were studied. The experimental measurements on the rabbit models of control and osteoporosis groups were performed. The PA spectroscopic curves of the bone from different groups were obtained and decomposed. Then, the MWPA parameter “relative composition ratio” of different chemical components including hydroxyapatite, lipid, hemoglobin, oxy-hemoglobin, and whole blood were quantified and compared with the gold-standard X-ray images and the relative optical absorption spectrum obtained by the commercial spectrometer system.
Materials and Methods
The bone composition of the two groups was measured by using MWPA at the range of 690–950 nm at 10 nm intervals. The experimental setup for studying the chemical components in bone is shown in Figure 1A. The light beam generated by an Nd:YAG laser pumped OPO (Vibrant B, Opotek) was divided into two parts. 10% of the laser energy was projected to a black rubber by a beam splitter and recorded by the ultrasound transducer (V310-SU, Olympus) for the calibration of subsequent signal magnitude. The remaining 90% illuminated the bone from one side on the surface of the bone. The diameter of the beam was 4 mm and the light fluency was controlled at 15–20 mJ/cm2. Light passed through the bone and then excited PA signal. The PA signal was received by a needle hydrophone (HNC-1500, Onda Co., Sunnyvale, CA, United States) with a broad bandwidth from 0 to 10 MHz. A pre-amplifier was connected after the hydrophone to improve the signal-to-noise ratio (SNR), and then it was digitized and recorded by a digital oscilloscope (HDO6000, oscilloscope, Teledyne Lecroy, United States). To enhance the signal-to-noise ratio (SNR), the PA signal was averaged over 50 laser pulses. Furthermore, the gold-standard DEXA images was conducted by the commercial DEXA imaging system (InAlyzer) and the relative optical absorption spectrum obtained by the commercial spectrometer system were conducted for each bone samples.
FIGURE 1. (A,B) The schematic diagram of the experiment setup. (C) The photograph of the bone sample. (D) An example of the PA signal generated by the bone specimen. (E) An example of the power spectrum density (PSD) of the PA signal generated by the bone.
In this study, the animal bone models we used in the MWPA measurements were distal end of forelimb, as shown in Figure 1C. 8 five-month-old, skeletally mature, female New Zealand white rabbits were divided randomly into two groups: osteoporosis group and control group. Bilateral ovariectomy was used in the osteoporosis group to simulate the symptoms of osteoporosis in elderly women, and sham surgery was used in the control group to avoid other factors affecting the experimental results. Twenty weeks after surgery, rabbits were euthanized, and the distal end of forelimb were dissected and subject to PA assessment.
First, the PA signal received by the Onda hydrophone was calibrated with the PA signal amplitude generated by the black rubber. The signal was related to the laser power of each wavelength. Secondly, the PA signal was transmitted to the frequency domain via fast Fourier transform, and the intensity of the PA frequency power at 0.2–5 MHz was summarized as the PA absorption value of each wavelength without being affected by low frequency (<0.2 MHz) noise or high frequency (>5 MHz) noise. The PA signal of each wavelength of light was then quantified and a PA absorption curve was obtained, which represented the spectral PA absorption of each bone. Thirdly, each bone sample was tested from three different directions in order to reduce measurement errors. The PA spectrum curves were obtained in at three different locations as shown in Figure 1C and were averaged for further analysis. Finally, spectral unmixing based on the least-square regression method was conducted, the PA absorption spectrum of each group of bones was decomposed to obtain the proportion of corresponding chemical components.
The DEXA imaging results from osteoporosis group and control group are shown in Figure 2. DEXA images showing the regions of interests (ROI) is marked by the yellow circle as shown in Figure 2B. Correlations between the BMD results from DEXA were studied by using unpaired two-tailed independent samples t-tests (with Welch’s correction in cases of unequal variances), which lead to p < 0.01, as shown in Figure 2C. This study based on commercial DEXA technologies confirmed the pathologic conditions of the osteoporosis group as well as the difference between the two groups of rabbit models.
FIGURE 2. DEXA images results. (A) Commercial DEXA imaging system. (B) DEXA images showing the regions of interests (ROI). (C) DEXA imaging at the ROI demonstrates significantly reduced BMD in the osteoporosis group vs. control group. ** stands for p < 0.01 in unpaired t-test comparing the results from the two bone groups (N = 8 for osteoporosis group, N = 8 for control group).
Optical Absorption Measurements
The results of optical absorption measurements are shown in Figure 3. Figures 3AB shows the commercial spectrometer system and the working principle. The mean and standard deviations of the relative optical absorption spectrum of the bone samples from the osteoporosis group and control group are compared in Figure 3C. By comparing with the control group, the bone sample from osteoporosis group has higher optical absorption at the range of 800–950 nm and lower optical absorption at the range of 690–800 nm.
FIGURE 3. The optical absorption system and results. (A) The commercial transmission and reflection spectrometer system used in this study. (B) The schematic diagram of the integrating sphere used in the commercial spectrometer system. (C)The relative optical absorption spectrum of the bone from the control group and osteoporosis group measured by the transmission and reflection spectrometer, each normalized at 800 nm.
Multi-Wavelength Photoacoustic Measurements
In the spectral range of 690–950 nm, the main optical absorption components in the bone are oxygenated hemoglobin, deoxygenated hemoglobin, mineral (mostly hydroxyapatite), lipid. Figure 4 shows the MWPA results of the osteoporosis group and control group. The two solid lines in Figure 4A show the averaged PA spectra measured from the two groups of bone samples, while the standard deviation is shown by the shaded area next to each curve. It is obviously that the absorption reached its peaks at 700, 760, and 930 nm. The difference is that the absorption of the osteoporosis group near 700–760 nm is stronger than that of the control group, while the absorption of the control group near 930 nm is stronger than that of the osteoporosis group. As expected, the corresponding content of the components in the osteoporosis group and the control group is consistent with the optical absorption results obtained by the commercial spectrometer system shown in Figure 3C. The optical absorption spectra of the main chemical components in the bone are shown in Figure 4B. By comparing the wavelength positions of the PA signal peaks and the component absorption peaks, it can be found that the strong absorption peak at 700 nm is mainly caused by the absorption of deoxygenated hemoglobin and hydroxyapatite. The absorption peak at 930 nm is mainly contributed by the lipid. Therefore, it can be preliminarily estimated that the lipid content of the osteoporosis group is higher than that of the control group, which is consistent with the past study [1, 9]. Besides, the deoxygenated hemoglobin content and the hydroxyapatite content are lower than that of the control group in the animal bone models we used.
FIGURE 4. (A) The averaged PA spectra measured from the osteoporosis group and control group. (B) Optical absorption spectra of major chemical components in the bone, including deoxygenated hemoglobin (Hb), lipids, hydroxyapatite, and oxygenated hemoglobin (HbO2).
For the PA absorption spectrum for each bone sample, it is given as 
FIGURE 5. (A) The averaged PA spectrum in comparison with the fitted spectrum by using the least-square method. (B) The unmixed spectrum of the MWPA curves for the sample from osteoporosis group and control group, respectively. (C)The relative contents of the hydroxyapatite, lipids, deoxygenated hemoglobin (Hb), and whole blood (Hb+HbO2). * stands for p < 0.05.
Conclusion and Discussion
The results of this study indicate the feasibility of MWPA method in assessing bone chemical composition. In addition to measuring the mineral content in bones, it can also measure the content of chemical components such as lipid, oxygenated hemoglobin and deoxygenated hemoglobin. The content and changes of these ingredients are inextricably related to bone health, and has the potential to used as a new way for bone health assessment.
This study currently has some limitations. First, the number of bone samples is not large enough. There are only four rabbit samples in each group, and only eight sets of data can be obtained, so data analysis is not universal. In the future, additional data should be collected and analyzed on a large group of samples. Second, the chemical composition has not been verified using the bone histomorphometry. In the future work, we plan to further indicate the differences in lipid and Hb content of those two groups though other pathways. For example, the MRI imaging can be used to get the lipid fraction as the gold-standard for lipid content [10, 36, 37]. The reticular fiber staining for quantification of blood vessels can be used as the gold-standard to reveal the differences in blood content of those two groups. Third, the influence of different types of fatty acids as well as the content of collagen which are also related with bone health are not considered in this study. In the future, with the wider range of the laser wavelength and higher resolution of scanning, it has the potential to distinguish the changes in different types of fatty acids and contents of collagen in the bone tissue. Forth, this study is mainly focusing on the isolated bones. In the future, this technology has the potential to be applied for clinical study in vivo. For in vivo study, the bone assessment technique based on the PA detecting method described in our manuscript may be not available and should be optimized with better design. For example, the transmission PA mode for multi-wavelength PA measurement of bone in vivo in our study published recently . In addition, the effects of light attenuation and ultrasound attenuation need to be considered. Besides, not only the bone but also the overlying soft tissue absorb the light and leads to an increase of light attenuation. Therefore, it is necessary to consider the effect of light attenuation in the bone tissue as well as the overlying soft tissue. At the same time, it has a high attenuation of PA signal for the in vivo PA measurements, especially for the high frequency components. Thus, the optimized center frequency of the transducer at a relative lower frequency should be used.
Despite these limitations, this study has successfully demonstrated the feasibility of the emerging PA technology in assessing the chemical information in bones. Compared with DEXA, MRI, QUS and light based techniques, MWPA has the potential to provides more comprehensive bone assessment results, including not only the chemical compositions of bone but also the micro-structure which is highly correlated with bone health as reported in our previous paper . For the gold standard DEXA imaging, since the X-rays are attenuated by both bone and bone marrow fat, therefore, the BMD results quantified by DEXA is altered by the changes in the lipid content of bones. For MRI, in particular MRI spectroscopy, it allows the quantification of bone marrow fats, however, it cannot provide the information of BMD . In addition, the QUS cannot be used for assessing the components of bone composition. Besides, compared with the light based techniques, the PA sensing is based on the detection of light-induced ultrasonic signals which are much less scattered in biological tissues, it can present more spatial information in deep tissues than pure optical techniques. Furthermore, the PA bone assessment method proposed in this article has many advantages such as target-specific, non-ionizing, low-cost, and patient friendly. With all these unique advantages, this technology is expected to be improved and developed into a more accurate bone evaluation method in the future.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
The animal study was reviewed and approved by Animal Ethics Committee of Tongji University, Shanghai.
TF, YX, and WX conceived the study. TF and YX performed data acquisition. TF and YX analyzed all the PA data. TF and YX developed the data processing algorithms, TF, YX, WX, DT, and QC drafted the manuscript.
This work was supported by the National Key Research and Development Project (nos. 2017YFC0111400 and 2016YFA0100800), National Natural Science Foundation of China (nos. 11704188 and 11827808, 12034015, 11674249, 81702962), the Natural Science Foundation of Jiangsu, China (no. BK 20170826), and the Postdoctoral Science Foundation of China under grant no. 2019M651564.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
1. Di Pietro G, Capuani S, Manenti G, Vinicola V, Fusco A, Baldi J, et al. . Bone marrow lipid profiles from peripheral skeleton as potential biomarkers for osteoporosis: a 1H-MR spectroscopy study. Acad Radiol (2016) 23(3):273–83. doi:10.1016/j.acra.2015.11.009
4. Raum K, Leguerney I, Chandelier F, Bossy E, Talmant M, Saïed A, et al. . Bone microstructure and elastic tissue properties are reflected in QUS axial transmission measurements. Ultrasound Med Biol (2005) 31(9):1225–35. doi:10.1016/j.ultrasmedbio.2005.05.002
5. Laugier P. Quantitative ultrasound instrumentation for bone in vivo characterization. In: P. Laugier, and G. Haïat, editors Bone quantitative ultrasound. Dordrecht, Netherlands: Springer (2011). p 47–71.
6. Töyräs J, Nieminen M, Kröger H, Jurvelin J. Bone mineral density, ultrasound velocity, and broadband attenuation predict mechanical properties of trabecular bone differently. Bone (2002) 31(4):503–7. doi:10.1016/s8756-3282(02)00843-8
7. Liu C, Ta D, Fujita F, Hachiken T, Matsukawa M, Mizuno K, et al. . The relationship between ultrasonic backscatter and trabecular anisotropic microstructure in cancellous bone. J Appl Phys (2014) 115(6):064906. doi:10.1063/1.4865173
9. Griffith JF, Yeung DK, Antonio GE, Lee FK, Hong AW, Wong SY, et al. . Vertebral bone mineral density, marrow perfusion, and fat content in healthy men and men with osteoporosis: dynamic contrast-enhanced MR imaging and MR spectroscopy. Radiology (2005) 236(3):945–51. doi:10.1148/radiol.2363041425
10. Griffith JF, Yeung DK, Antonio GE, Wong SY, Kwok TC, Woo J, et al. . Vertebral marrow fat content and diffusion and perfusion indexes in women with varying bone density: MR evaluation. Radiology (2006) 241(3):831–8. doi:10.1148/radiol.2413051858
12. Pifferi A, Torricelli A, Taroni P, Bassi AL, Chikoidze E, Giambattistelli E, et al. . Optical biopsy of bone tissue: a step toward the diagnosis of bone pathologies. J Biomed Optic (2004) 9(3):474–80. doi:10.1117/1.1691029
13. Schulmerich MV, Dooley K, Morris MD, Vanasse TM, Goldstein SA. Transcutaneous fiber optic Raman spectroscopy of bone using annular illumination and a circular array of collection fibers. J Biomed Optic (2006) 11(6):060502. doi:10.1117/1.2400233
15. Draper ER, Morris MD, Camacho NP, Matousek P, Towrie M, Parker AW, et al. . Novel assessment of bone using time‐resolved transcutaneous Raman spectroscopy. J Bone Miner Res (2005) 20(11):1968–72. doi:10.1359/JBMR.050710
19. Wang X, Pang Y, Ku G, Xie X, Stoica G, Wang LV. Noninvasive laser-induced photoacoustic tomography for structural and functional in vivo imaging of the brain. Nat Biotechnol (2003) 21(7):803–6. doi:10.1038/nbt839
22. Zhang W, Oraiqat I, Lei H, Carson PL, EI Naqa I, Wang X. Dual-modality X-ray-induced radiation acoustic and ultrasound imaging for real-time monitoring of radiotherapy. BME Front (2020) 2020:9853609. doi:10.34133/2020/9853609
23. Zhu Y, Lu X, Dong X, Yuan J, Fabiilli ML, Wang X. LED-based photoacoustic imaging for monitoring angiogenesis in fibrin scaffolds. Tissue Eng C Methods (2019) 25(9):523–31. doi:10.1089/ten.TEC.2019.0151
24. Lashkari B, Yang L, Mandelis A. The application of backscattered ultrasound and photoacoustic signals for assessment of bone collagen and mineral contents. Quant Imaging Med Surg (2015) 5(1):46–56. doi:10.3978/j.issn.2223-4292.2014.11.11
26. Kim MW, Jeng GS, O’Donnell M, Pelivanov I. Correction of wavelength-dependent laser fluence in swept-beam spectroscopic photoacoustic imaging with a hand-held probe. Photoacoustics (2020) 19:100192. doi:10.1016/j.pacs.2020.100192
29. Udelhoven T, Emmerling C, Jarmer T. Quantitative analysis of soil chemical properties with diffuse reflectance spectrometry and partial least-square regression: a feasibility study. Plant Soil (2003) 251(2):319–29. doi:10.1023/A:1023008322682
30. Patsch JM, Li X, Baum T, Yap SP, Karampinos DC, Schwartz AV, et al. . Bone marrow fat composition as a novel imaging biomarker in postmenopausal women with prevalent fragility fractures. J Bone Miner Res (2013) 28(8):1721–8. doi:10.1002/jbmr.1950
33. Shen W, Scherzer R, Gantz M, Chen J, Punyanitya M, Lewis CE, et al. . Relationship between MRI-measured bone marrow adipose tissue and hip and spine bone mineral density in African-American and Caucasian participants: the CARDIA study. J Clin Endocrinol Metab (2012) 97(4):1337–46. doi:10.1210/jc.2011-2605
34. Bredella MA, Daley SM, Kalra MK, Brown JK, Miller KK, Torriani M. Marrow adipose tissue quantification of the lumbar spine by using dual-energy CT and single-Voxel (1)H MR spectroscopy: a feasibility study. Radiology (2015) 277(1):230–5. doi:10.1148/radiol.2015142876
36. Patsch JM, Li X, Baum T, Yap SP, Karampinos DC, Schwartz AV, et al. . Bone marrow fat composition as a novel imaging biomarker in postmenopausal women with prevalent fragility fractures. J Bone Miner Res (2013) 28(8):1721–8. doi:10.1002/jbmr.1950
37. Baofeng L, Zhi Y, Bei C, Guolin M, Qingshui Y, Jian L. Characterization of a rabbit osteoporosis model induced by ovariectomy and glucocorticoid. Acta Orthop (2010) 81(3):396–401. doi:10.3109/17453674.2010.483986
Keywords: photoacoustic, bone assessment, chemical composition, multi-wavelength photoacoustic analysis, osteoporosis
Citation: Feng T, Xie Y, Xie W, Ta D and Cheng Q (2020) Bone Chemical Composition Analysis Using Photoacoustic Technique. Front. Phys. 8:601180. doi: 10.3389/fphy.2020.601180
Received: 01 September 2020; Accepted: 18 November 2020;
Published: 16 December 2020.
Edited by:Jun Xia, University at Buffalo, United States
Reviewed by:Silvia Capuani, National Research Council (CNR), Italy
Shu-Chi A. Yeh, Massachusetts General Hospital, United States
Copyright © 2020 Feng, Xie, Xie, Ta and Cheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qian Cheng, email@example.com