1
|
Gonzalez EA, Bell MAL. Photoacoustic Imaging and Characterization of Bone in Medicine: Overview, Applications, and Outlook. Annu Rev Biomed Eng 2023; 25:207-232. [PMID: 37000966 DOI: 10.1146/annurev-bioeng-081622-025405] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2023]
Abstract
Photoacoustic techniques have shown promise in identifying molecular changes in bone tissue and visualizing tissue microstructure. This capability represents significant advantages over gold standards (i.e., dual-energy X-ray absorptiometry) for bone evaluation without requiring ionizing radiation. Instead, photoacoustic imaging uses light to penetrate through bone, followed by acoustic pressure generation, resulting in highly sensitive optical absorption contrast in deep biological tissues. This review covers multiple bone-related photoacoustic imaging contributions to clinical applications, spanning bone cancer, joint pathologies, spinal disorders, osteoporosis, bone-related surgical guidance, consolidation monitoring, and transsphenoidal and transcranial imaging. We also present a summary of photoacoustic-based techniques for characterizing biomechanical properties of bone, including temperature, guided waves, spectral parameters, and spectroscopy. We conclude with a future outlook based on the current state of technological developments, recent achievements, and possible new directions.
Collapse
Affiliation(s)
- Eduardo A Gonzalez
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland, USA
| | - Muyinatu A Lediju Bell
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland, USA
- Department of Electrical and Computer Engineering and Department of Computer Science, Johns Hopkins University, Baltimore, Maryland, USA;
| |
Collapse
|
2
|
Porrelli D, Abrami M, Pelizzo P, Formentin C, Ratti C, Turco G, Grassi M, Canton G, Grassi G, Murena L. Trabecular bone porosity and pore size distribution in osteoporotic patients - A low field nuclear magnetic resonance and microcomputed tomography investigation. J Mech Behav Biomed Mater 2021; 125:104933. [PMID: 34837800 DOI: 10.1016/j.jmbbm.2021.104933] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2021] [Revised: 10/10/2021] [Accepted: 10/25/2021] [Indexed: 01/19/2023]
Abstract
The study of bone morphology is of great importance as bone morphology is influenced by factors such as age and underlying comorbidities and is associated with bone mechanical properties and fracture risk. Standard diagnostic techniques used in bone disease, such as Dual-Energy X-ray absorptiometry and ultrasonography do not provide qualitative and quantitative morphological information. In recent years, techniques such as High Resolution Computed Tomography (HR-CT), micro- CT, Magnetic Resonance Imaging (MRI), and Low Field Nuclear Magnetic Resonance (LF-NMR) have been developed for the study of bone structure and porosity. Data obtained from these techniques have been used to construct models to predict bone mechanical properties thanks to finite element analysis. Cortical porosity has been extensively studied and successfully correlated with disease progression and mechanical properties. Trabecular porosity and pore size distribution, however, have increasingly been taken into consideration to obtain a comprehensive analysis of bone pathology and mechanic. Therefore, we have decided to evaluate the ability of micro- CT (chosen for its high spatial resolving power) and LF-NMR (chosen to analyze the behavior of water molecules within trabecular bone pores) to characterize the morphology of trabecular bone in osteoporosis. Trabecular bone samples from human femoral heads collected during hip replacement surgery were from osteoporosis (test group) and osteoarthritis (control group) patients. Our data show that both micro- CT and LF-NMR can detect qualitative changes in trabecular bone (i.e., transition from plate-like to rod-like morphology). Micro- CT failed to detect significant differences in trabecular bone morphology parameters between osteoporotic and osteoarthritic specimens, with the exception of Trabecular Number and Connectivity Density, which are markers of osteoporosis progression. In contrast, LF-NMR was able to detect significant differences in porosity and pore size of trabecular bone from osteoporotic versus osteoarthritic (control) samples. However, only the combination of these two techniques allowed the detection of structural morphometric changes (increase in the larger pore fraction and enlargement of the larger pores) in the trabecular bone of osteoporotic specimens compared to osteoarthritic ones. In conclusion, the combined use of LF-NMR and micro- CT provides a valuable tool for characterizing the morphology of trabecular bone and may offer the possibility for a new approach to the study and modeling of bone mechanics in the context of aging and disease.
Collapse
Affiliation(s)
- Davide Porrelli
- Department of Medicine, Surgery and Health Sciences, Maggiore Ospital, Trieste University, Piazza dell'Ospitale 1, I-34125, Trieste, Italy
| | - Michela Abrami
- Department of Engineering and Architecture, University of Trieste, Via Valerio 6/A, I 34127, Trieste, Italy
| | - Patrizia Pelizzo
- Department of Medicine, Surgery and Health Sciences, University of Trieste, Strada di Fiume 447, I-34149, Trieste, Italy
| | - Cristina Formentin
- Department of Medicine, Surgery and Health Sciences, University of Trieste, Strada di Fiume 447, I-34149, Trieste, Italy
| | - Chiara Ratti
- Department of Medicine, Surgery and Health Sciences, University of Trieste, Strada di Fiume 447, I-34149, Trieste, Italy
| | - Gianluca Turco
- Department of Medicine, Surgery and Health Sciences, Maggiore Ospital, Trieste University, Piazza dell'Ospitale 1, I-34125, Trieste, Italy
| | - Mario Grassi
- Department of Engineering and Architecture, University of Trieste, Via Valerio 6/A, I 34127, Trieste, Italy.
| | - Gianluca Canton
- Department of Medicine, Surgery and Health Sciences, University of Trieste, Strada di Fiume 447, I-34149, Trieste, Italy
| | - Gabriele Grassi
- Department of Life Sciences, Cattinara University Hospital, Trieste University, Strada di Fiume 447, I-34149, Trieste, Italy
| | - Luigi Murena
- Department of Medicine, Surgery and Health Sciences, University of Trieste, Strada di Fiume 447, I-34149, Trieste, Italy
| |
Collapse
|
3
|
Ni Q, Hua R, Holland D, Tinajero A, Han Y, Jiang JX, Wang X. Characterization of Microstructural Changes on Biglycan Induced Mice Bone by Low-Field Nuclear Magnetic Resonance. APPLIED PHYSICS (KOWLOON, CHINA) 2021; 4:58-67. [PMID: 35479609 PMCID: PMC9040680 DOI: 10.31058/j.ap.2021.42004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
A NMR spin-spin (T2) relaxation technique has been described for determining the porosity, and the bound water distribution in biglycan induced mouse bone and correlate to their mechanical properties. The technique of low-field proton NMR involves spin-spin relaxation and free induction decay (FID) measurements, and the computational inversion methods for decay data analysis. The CPMG T2 relaxation data can be inverted to T2 relaxation distribution and this distribution then can be transformed to a pore size distribution with the longer relaxation times corresponding to larger pores. The FID T2 relaxation data of dried bone (mobile water removed) can be inverted to T2 relaxation distribution and this distribution then can be transformed to bound and solid-like water distribution with the longest relaxation time corresponding to bound water component. These techniques are applied to quantify apparent changes in porosity, and bound water in controlled and biglycan knockout mouse bone. Overall bone porosity from CPMG T2 relaxation is determined using the calibrated NMR fluid volume from the proton relaxation data divided by overall bone volume. Ignore the physical sample differences, from the inversion FID T2 relaxation spectrum, the ratio of the bound to solid-like water components is used to calibrate the bound water inside bone, and the results can be used to correlated bone mechanical properties. Hydration status significantly affects the toughness of bone, and bound water has been considered as a biomarker for prediction of bone fragility fractures. In addition to the collagen phase, recent evidence shows that glycosaminoglycans (GAGs) of proteoglycans (PGs) in the extracellular matrix also play a pivotal role in regulating the tissue-level hydration status of bone, there by affecting the tissue-level toughness of bone. Furthermore, biglycan and decorin are two major types of PGs in bone reports. Biglycan knockout induced changes in GAGs, bound water, as well as bone tissue toughness. Among all subtypes of PGs, biglycan is identified as a major subtype in the bone mineral matrix. In this study, we used a biglycan mouse model and the obtained bone samples were measured by low-field NMR to determine the bone porosity and bound water changes, and used to predict if knockout of biglycan may affect the amount of bound water and subsequently lead to reduce toughness of bone.
Collapse
Affiliation(s)
- Qingwen Ni
- Department of Mathematics and Physics,Texas A&M International University, Laredo, TX, USA
| | - Rui Hua
- Department of Mechanical Engineering, The University of Texas at San Antonio, San Antonio, TX, USA
| | - Douglas Holland
- Department of Mathematics and Physics,Texas A&M International University, Laredo, TX, USA
| | - Anahi Tinajero
- Department of Mathematics and Physics,Texas A&M International University, Laredo, TX, USA
| | - Yan Han
- Department of Mechanical Engineering, The University of Texas at San Antonio, San Antonio, TX, USA
| | - Jean X. Jiang
- Department of Biochemistry and Structural Biology, The University of Texas Health Science Center, San Antonio, TX, USA
| | - Xiaodu Wang
- Department of Mechanical Engineering, The University of Texas at San Antonio, San Antonio, TX, USA
| |
Collapse
|
4
|
Nikel O, Poundarik AA, Bailey S, Vashishth D. Structural role of osteocalcin and osteopontin in energy dissipation in bone. J Biomech 2018; 80:45-52. [PMID: 30205977 DOI: 10.1016/j.jbiomech.2018.08.014] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2017] [Revised: 08/13/2018] [Accepted: 08/14/2018] [Indexed: 11/29/2022]
Abstract
Non-collagenous proteins are a vital component of bone matrix. Amongst them, osteocalcin (OC) and osteopontin (OPN) hold special significance due to their intimate interaction with the mineral and collagenous matrix in bone. Both proteins have been associated with microdamage and fracture, but their structural role in energy dissipation is unclear. This study used bone tissue from genetic deficient mice lacking OC and/or OPN and subjected them to a series of creep-fatigue-creep tests. To this end, whole tibiae were loaded in four-point bending to 70% stiffness loss which captured the three characteristic phases of fatigue associated with initiation, propagation, and coalescence of microdamage. Fatigue loading preceded and followed creep tests to determine creep and dampening parameters. Microdamage in the form of linear microcracks and diffuse damage were analyzed by histology. It was shown that OC and OPN were 'activated' following stiffness loss associated with fatigue damage where they facilitated creep and dampening parameters (i.e. increased energy dissipation). More specifically, post-fatigue creep rate and dampening were significantly greater in wild-types (WTs) than genetic deficient mice (p < 0.05). These results were supported by microdamage analysis which showed significant increase in creep-associated diffuse damage formation in WTs compared to genetic deficient groups (p < 0.05). Based on these findings, we propose that during local yield events, OC and OPN rely on ionic interactions of their charged side chains and on hydrogen bonding to dissipate energy in bone.
Collapse
Affiliation(s)
- Ondřej Nikel
- Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute, Troy, NY, USA; Department of Biomedical Engineering, Rensselaer Polytechnic Institute, Troy, NY, USA
| | - Atharva A Poundarik
- Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute, Troy, NY, USA; Department of Biomedical Engineering, Rensselaer Polytechnic Institute, Troy, NY, USA
| | - Stacyann Bailey
- Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute, Troy, NY, USA; Department of Biomedical Engineering, Rensselaer Polytechnic Institute, Troy, NY, USA
| | - Deepak Vashishth
- Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute, Troy, NY, USA; Department of Biomedical Engineering, Rensselaer Polytechnic Institute, Troy, NY, USA.
| |
Collapse
|
5
|
Abstract
The mechanical properties of bone are fundamental to the ability of our skeletons to support movement and to provide protection to our vital organs. As such, deterioration in mechanical behavior with aging and/or diseases such as osteoporosis and diabetes can have profound consequences for individuals' quality of life. This article reviews current knowledge of the basic mechanical behavior of bone at length scales ranging from hundreds of nanometers to tens of centimeters. We present the basic tenets of bone mechanics and connect them to some of the arcs of research that have brought the field to recent advances. We also discuss cortical bone, trabecular bone, and whole bones, as well as multiple aspects of material behavior, including elasticity, yield, fracture, fatigue, and damage. We describe the roles of bone quantity (e.g., density, porosity) and bone quality (e.g., cross-linking, protein composition), along with several avenues of future research.
Collapse
Affiliation(s)
- Elise F Morgan
- Orthopaedic and Developmental Biomechanics Laboratory, Department of Mechanical Engineering, Boston University, Boston, Massachusetts 02215, USA;
| | - Ginu U Unnikrisnan
- Orthopaedic and Developmental Biomechanics Laboratory, Department of Mechanical Engineering, Boston University, Boston, Massachusetts 02215, USA;
| | - Amira I Hussein
- Orthopaedic and Developmental Biomechanics Laboratory, Department of Mechanical Engineering, Boston University, Boston, Massachusetts 02215, USA;
| |
Collapse
|
6
|
Luczynski KW, Steiger-Thirsfeld A, Bernardi J, Eberhardsteiner J, Hellmich C. Extracellular bone matrix exhibits hardening elastoplasticity and more than double cortical strength: Evidence from homogeneous compression of non-tapered single micron-sized pillars welded to a rigid substrate. J Mech Behav Biomed Mater 2015; 52:51-62. [DOI: 10.1016/j.jmbbm.2015.03.001] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2014] [Revised: 02/13/2015] [Accepted: 03/01/2015] [Indexed: 11/17/2022]
|
7
|
Granke M, Does MD, Nyman JS. The Role of Water Compartments in the Material Properties of Cortical Bone. Calcif Tissue Int 2015; 97:292-307. [PMID: 25783011 PMCID: PMC4526331 DOI: 10.1007/s00223-015-9977-5] [Citation(s) in RCA: 152] [Impact Index Per Article: 16.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/16/2014] [Accepted: 02/27/2015] [Indexed: 12/22/2022]
Abstract
Comprising ~20% of the volume, water is a key determinant of the mechanical behavior of cortical bone. It essentially exists in two general compartments: within pores and bound to the matrix. The amount of pore water-residing in the vascular-lacunar-canalicular space-primarily reflects intracortical porosity (i.e., open spaces within the matrix largely due to Haversian canals and resorption sites) and as such is inversely proportional to most mechanical properties of bone. Movement of water according to pressure gradients generated during dynamic loading likely confers hydraulic stiffening to the bone as well. Nonetheless, bound water is a primary contributor to the mechanical behavior of bone in that it is responsible for giving collagen the ability to confer ductility or plasticity to bone (i.e., allows deformation to continue once permanent damage begins to form in the matrix) and decreases with age along with fracture resistance. Thus, dehydration by air-drying or by solvents with less hydrogen bonding capacity causes bone to become brittle, but interestingly, it also increases stiffness and strength across the hierarchical levels of organization. Despite the importance of matrix hydration to fracture resistance, little is known about why bound water decreases with age in hydrated human bone. Using (1)H nuclear magnetic resonance (NMR), both bound and pore water concentrations in bone can be measured ex vivo because the proton relaxation times differ between the two water compartments, giving rise to two distinct signals. There are also emerging techniques to measure bound and pore water in vivo with magnetic resonance imaging (MRI). The NMR/MRI-derived bound water concentration is positively correlated with both the strength and toughness of hydrated bone and may become a useful clinical marker of fracture risk.
Collapse
Affiliation(s)
- Mathilde Granke
- Department of Orthopaedic Surgery & Rehabilitation, Vanderbilt University Medical Center, Nashville, TN 37232
- Center for Bone Biology, Vanderbilt University Medical Center, Nashville, TN 37232
| | - Mark D. Does
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN 37232
- Institute of Imaging Science, Vanderbilt University, Nashville, TN 37232
- Department of Radiology and Radiological Sciences, Vanderbilt University, Nashville, TN 37232
- Department of Electrical Engineering, Vanderbilt University, Nashville, TN 37232
| | - Jeffry S. Nyman
- Department of Orthopaedic Surgery & Rehabilitation, Vanderbilt University Medical Center, Nashville, TN 37232
- Center for Bone Biology, Vanderbilt University Medical Center, Nashville, TN 37232
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN 37232
- Department of Veterans Affairs, Tennessee Valley Healthcare System, Nashville, TN 37212
| |
Collapse
|
8
|
Şipos RS, Fechete R, Moldovan D, Şuş I, Szasz S, Pávai Z. Assessment of femoral bone osteoporosis in rats treated with simvastatin or fenofibrate. Open Life Sci 2015. [DOI: 10.1515/biol-2015-0039] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
AbstractBackground: The effects of two lipidlowering drugs, simvastatin and fenofibrate on osteoporosis in the femurs of healthy and ovariectomized female rats were investigated quantitatively by histological images and
Collapse
|
9
|
Haupert S, Guérard S, Peyrin F, Mitton D, Laugier P. Non destructive characterization of cortical bone micro-damage by nonlinear resonant ultrasound spectroscopy. PLoS One 2014; 9:e83599. [PMID: 24392089 PMCID: PMC3879251 DOI: 10.1371/journal.pone.0083599] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2013] [Accepted: 11/05/2013] [Indexed: 01/22/2023] Open
Abstract
The objective of the study was to evaluate the ability of a nonlinear ultrasound technique, the so-called nonlinear resonant ultrasound spectroscopy (NRUS) technique, for detecting early microdamage accumulation in cortical bone induced by four-point bending fatigue. Small parallelepiped beam-shaped human cortical bone specimens were subjected to cyclic four-point bending fatigue in several steps. The specimens were prepared to control damage localization during four-point bending fatigue cycling and to unambiguously identify resonant modes for NRUS measurements. NRUS measurements were achieved to follow the evolution of the nonlinear hysteretic elastic behavior during fatigue-induced damage. After each fatigue step, a small number of specimens was removed from the protocol and set apart to quantitatively assess the microcrack number density and length using synchrotron radiation micro-computed tomography (SR-µCT). The results showed a significant effect of damage steps on the nonlinear hysteretic elastic behavior. No significant change in the overall length of microcracks was observed in damaged regions compared to the load-free control regions. Only an increased number of shortest microcracks, those in the lowest quartile, was noticed. This was suggestive of newly formed microcracks during the early phases of damage accumulation. The variation of nonlinear hysteretic elastic behavior was significantly correlated to the variation of the density of short microcracks. Our results suggest that the nonlinear hysteretic elastic behavior is sensitive to early bone microdamage. Therefore NRUS technique can be used to monitor fatigue microdamage progression in in vitro experiments.
Collapse
Affiliation(s)
- Sylvain Haupert
- UPMC Univ Paris 06, CNRS UMR7623, Laboratoire d’Imagerie Paramétrique, Paris, France
- * E-mail:
| | | | - Françoise Peyrin
- CREATIS, INSERM U1044, CNRS 5220, INSA Lyon, Université Lyon 1, Lyon, France
- European Synchrotron Radiation Facility, Grenoble, France
| | - David Mitton
- Université de Lyon, IFSTTAR, LBMC, UMR_T 9406, Université Lyon 1, Lyon, France
| | - Pascal Laugier
- UPMC Univ Paris 06, CNRS UMR7623, Laboratoire d’Imagerie Paramétrique, Paris, France
| |
Collapse
|
10
|
Chan KS, Nicolella DP. Micromechanical modeling of R-curve behaviors in human cortical bone. J Mech Behav Biomed Mater 2012; 16:136-52. [PMID: 23158217 DOI: 10.1016/j.jmbbm.2012.09.009] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2012] [Revised: 09/19/2012] [Accepted: 09/22/2012] [Indexed: 12/01/2022]
Abstract
The risk of bone fracture increases with age because of a variety of factors that include, among others, decreasing bone quantity and quality due to increasing porosity and crack density with age. Experimental evidence has indicated that changes in bone microstructure and trace mineralization with age can result in different crack-tip strain field and fracture response, leading to different fracture mechanisms and R-curve behaviors. In this paper, a micromechanical modeling approach is developed to predict the R-curve response of bone tissue by delineating fracture mechanisms that lead to microdamage and ligament bridging by incorporating the influence of increasing porosity and crack density with age. The effects of age on fracture of human femur cortical bone due to porosity (bone quantity) and bone quality (crack density) with age are then examined via the micromechanical model.
Collapse
Affiliation(s)
- Kwai S Chan
- Southwest Research Institute, 6220 Culebra Road, San Antonio, TX 78238, USA.
| | | |
Collapse
|
11
|
Li H, Liang CZ, Shen CC, Chen QX. Decreases in fluid shear stress due to microcracks: a possible primary pathogenesis of Kümmell's disease. Med Hypotheses 2011; 77:897-9. [PMID: 21873000 DOI: 10.1016/j.mehy.2011.08.005] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2011] [Revised: 07/16/2011] [Accepted: 08/02/2011] [Indexed: 12/18/2022]
Abstract
The German doctor Hermann Kümmell described Kümmell's disease as the clinical scenario in which patients suffer a trivial spinal trauma, but develop a symptomatic, progressive, angular kyphosis after a symptom-free period of months to years. Since an intravertebral vacuum phenomenon, which is considered indicative of ischemic osteonecrosis, is often seen in the radiographs of patients with Kümmell's disease, most authors regard ischemic necrosis of the vertebral body as the primary pathogenesis of Kümmell's disease. However, we argue that Kümmell's disease is not commonly associated with typical avascular osteonecrosis of the femoral head and the intravertebral vacuum phenomenon is also present in other diseases. We postulated that even if ischemia plays a role in the pathogenesis of Kümmell's disease, it would not be the proximal cause of Kümmell's disease. In this article, we review the role of fluid shear stress in bone remolding and propose a microcosmic hypothesis in which microcracks lead to decreased fluid shear stress, which acts as the primary cause of Kümmell's disease. This was supported by conclusions drawn from a literature review: (1) fluid shear stress plays a crucial role in bone remodeling, and the osteocyte network is the main sensor of this mechanical stimulus; (2) decreased fluid shear stress will cause disequilibration of bone homeostasis, increasing bone resorption and reducing bone formation; and (3) the fluid flow of lacunar-canalicular porosity (PLC) and fluid shear stress near the microcracks decreases.
Collapse
Affiliation(s)
- Hao Li
- Department of Orthopaedics, 2nd Affiliated Hospital, School of Medicine, Zhejiang University, #88 Jiefang Road, Hangzhou 310009, PR China
| | | | | | | |
Collapse
|