Multispectral optoacoustic tomography of peripheral arterial disease based on muscle hemoglobin gradients—a pilot clinical study
Introduction
Intermittent claudication, the clinical hallmark of peripheral artery disease (PAD), is the result of impaired tissue perfusion and oxygenation due to partial or complete arterial obstruction in the affected limb (1). PAD is a common disease, associated with serious complications such as impaired walking endurance due to pain, development of chronic wounds with poor wound healing and increased risk for infections, decrease in quality of life or even limb amputation. Furthermore, PAD is known to increase the risk for cardiovascular events and death (2).
Diagnosis and grading of PAD includes clinical assessment and peripheral blood pressure assessments in form of ankle-brachial index (ABI) before and after treadmill exercise (3). Besides clinical assessment, PAD evaluation is based on non-invasive imaging playing a key role for the characterization of PAD and subsequent therapy planning (4). Current diagnostic strategies are based on duplex ultrasonography (DUS), magnetic resonance angiography (MRA) or computed tomography angiography (CTA). DUS provides tomographic imaging of arterial wall morphology and blood flow profiles of the diseased artery. MRA or CTA produce three-dimensional luminographies of the complete arterial tree of the lower limb, after the administration of contrast agents.
However, the above-mentioned methods provide only indirect assessments of disease stage by visualizing only macrovascular anatomical and hemodynamic information, such as the location and degree of arterial stenosis, without providing direct and objective information about the perfusion and oxygenation of the affected muscle. Therefore, novel diagnostic approaches are being explored to objectively assess muscle perfusion and oxygenation: contrast-enhanced ultrasound (CEUS) requires the intravenous injection of microbubbles to visualize the perfusion of skeletal muscle (5). Moreover, several magnetic resonance techniques have been developed to measure muscle perfusion (arterial spin labeling, ASL) and oxygenation (blood oxygenation-level dependent, BOLD) without the need for exogenous contrast (6). Near-infrared spectroscopy (NIRS) employs light in the near-infrared range (NIR, 680–980 nm) to resolve the distributions of oxy- (HbO2) and deoxyhemoglobin (Hb) in tissue or, less often, the distribution of the injected indocyanine green (ICG) (7,8). Nevertheless, NIRS perfusion and oxygenation recordings come with high uncertainty about the clear contributions of muscle or superficial soft tissues (especially subcutaneous fat and skin) in the measured signal. Diffuse optical tomography (DOT) employs near-infrared light to record tomographic images of perfusion and oxygenation for several centimeters through the examined tissue (9). However, both NIRS and DOT provide only low spatial resolution images (e.g., 5–10 mm) due to light scattering in tissue, and thus have not entered routine clinical practice.
Multispectral optoacoustic tomography (MSOT) has been introduced to overcome these limitations by taking advantage of the homonymous effect, wherein the emission and absorption of near infrared light pulses lead to the generation of ultrasound waves within the scanned tissue. This way, MSOT produces high-resolution (50–300 µm depending on the setup) images of blood vessels (10-12), as well as tissue perfusion and oxygenation based purely on hemoglobin gradients in space and time (13,14). MSOT has recently shown its unique ability to characterize vascular malformations (15) and to assess perfusion and oxygenation of soft tissues such as skin in systemic sclerosis (16), brown adipose tissue (17), thyroid gland (18), breast cancer tumors (19) and intestines (20) in patients with Crohn’s disease. Moreover, tissue perfusion is also disturbed under vascular occlusion conditions. First applications in healthy volunteers showed the potential of MSOT to quantitatively image skeletal muscle perfusion in iatrogenic occlusions of arteries and veins (13).
Thus, we hypothesized in this study that MSOT could also be used for the non-invasive and label-free assessment of skeletal muscle perfusion and oxygenation in PAD. By means of a clinical hybrid MSOT-ultrasound (US) system we measured the HbO2 and Hb content within the calf muscles of healthy volunteers under cuff-induced conditions of arterial occlusion, as well as of patients with previously diagnosed PAD before and after revascularization therapy.
Methods
Cuff-induced arterial occlusion measurements
To investigate the capability of MSOT to accurately image calf muscle perfusion and oxygenation over time, we conducted a set of controlled cuff-induced arterial occlusion measurements in healthy volunteers. Four (n1=4, 2 females, 2 males, age 31±3.2 years) non-smokers, without any vascular or muscular disease were included (Table 1). All participants were asked to consume no caffeine for at least 6 hours before the measurements to avoid its influence on vascular reactivity. Subjects lied in the prone position and a cuff of appropriate size was placed just above the knee of their dominant leg. The cuff was inflated to 40 mmHg above the systolic blood pressure of each subject. The MSOT-US probe was placed at the level of the largest calf circumference for each subject and kept stable during scanning in a transversal position for 1 minute before cuff inflation, 3 minutes during the occlusion and 2 minutes after cuff deflation. Each cuff measurement lasted for 15 min.
Full table
PAD study population and measurements
Four patients (n2=4, 1 female, 3 males, age 72±5.9) with PAD (Fontaine Stage IIb; Rutherford Grade I, Category 2–3), previously diagnosed by means of US, MRA or CTA were included (Table 1). Patients were asked to lie in supine position with the knees bent. They were then scanned over the affected limb before and 24 h after the revascularization with balloon and stent angioplasty of the superficial femoral and popliteal artery. No patient suffered from acute limb ischemia at the time point of the preoperative examination. Additionally, four healthy volunteers (n3=4, 1 female, 3 males, age 37.3±13.9) were also examined as controls (Table 1). The hand-held scanning probe was placed for both patients and healthy volunteers at the level of the largest calf circumference in a transversal position. Each scan lasted for approximately 10 seconds. Total duration of MSOT examinations was approximately 15 min, which is comparable to ultrasound.
The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). All participants gave written and informed consent prior to study enrolment. The study was reviewed by the local ethics committee (Ethikkommission der Westfälischen Wilhelms-Universität Münster, Protocol No. 2017-538-f-S).
MSOT setup and image acquisition
Technical details of the clinical hand-held MSOT-US system (MSOT Acuity Echo, iThera Medical GmbH, Munich) have been previously described (20,21). The central frequency of the ultrasound transducer was 4 MHz. The output light pulse energy (~15 mJ illuminating an area of 1 cm × 4 cm) adhered to the American National Standards Institute (ANSI) limits of maximum permissible human exposure (22). The central frequency MSOT achieved an imaging penetration depth of ~3–4 cm with a spatial resolution of ~250 µm. Co-registered US images (8 Hz frame rate) were generated and recorded for all MSOT frames (25 Hz frame rate). Experiments took place at normal room temperature (~22 °C). Both healthy volunteers and patients rested for at least 5 min before examination.
MSOT images were acquired at multiple different wavelengths in the NIR including 700, 730, 760, 780, 800 and 850 nm, where the HbO2 and the Hb demonstrate prominent spectral differences. Optoacoustic data for each pulse (wavelength) were reconstructed into tomographic images by means of a customized model-based algorithm (23). Regions of interest (ROIs) were delineated in the US images over the calf muscle in consensus between two examiners blinded for patient or healthy volunteer data. For each subject one characteristic US image was selected based on image quality and low motion artifacts. ROIs were defined based on the local anatomy and the known texture of muscle tissue in US. Levels of HbO2, Hb, total hemoglobin (THb) and hemoglobin saturation (SO2) were extracted from corresponding MSOT images via a spectral unmixing algorithm used in previous clinical MSOT studies (17). The THb image was calculated as the sum of the spectrally unmixed HbO2 and Hb images. The SO2 image was calculated in a per-pixel manner by dividing the HbO2 image by the resulted THb image. The MSOT principle of operation and the data workflow are graphically presented in Figure 1.
Statistical analysis
Data were analyzed using the R environment for statistical analysis and MATLAB© (Version R2018b, The MathWorks, Inc., Natick, MA, USA). A two-sided student’s t-tests was used for group comparisons when data where normal, and Wilcoxon rank sum tests where used when data deviated from normality. Results are shown in corresponding boxplots (Figure 2 for cuff measurements and Figure 3 for PAD imaging). P values <0.05 were considered to be statistically significant. Reported results are represented as the mean value ± one standard deviation.
Results
Cuff-induced arterial occlusion measurements
First, we scanned the calf muscle of four healthy volunteers during a 3-minute cuff-induced arterial occlusion, to investigate and demonstrate the capability of MSOT in monitoring skeletal muscle perfusion and oxygenation changes over time (Figure 2). Figure 2A illustrates the temporal fluctuations of Hb, HbO2, THb and SO2 for a healthy volunteer, as an example. Corresponding MSOT images of the same healthy volunteer as in Figure 2A, as well as the MSOT readouts across the entire cohort show a gradual increase of intramuscular Hb (Figure 2B,F) during the arterial occlusion period followed by a decrease after the release of the occlusion. Correspondingly, a gradual decrease in HbO2 (Figure 2C,F) is followed by a significant increase after the deflation of the cuff or else during the reactive hyperemia period. Furthermore, the THb (Figure 2D,F) remains relatively stable with increasing trend over the occlusion period with an increase after flow release and SO2 (Figure 2E,F) follows the same fluctuation patters as that of HbO2 described above. Table 2 summarizes our findings regarding the cuff-induced arterial occlusion measurements.
Full table
Measurements of patients with PAD
Next, we used MSOT to examine the differences between healthy volunteers and patients with PAD, as well as the impact of the interventional revascularization therapy on muscle hemodynamics and oxygenation (Figure 3). We extracted and plotted the levels of Hb (Figure 3A), HbO2 (Figure 3B), THb (Figure 3C) and SO2 (Figure 3D) from the calf muscle ROIs of the MSOT images. Our results show statistically significant differences between healthy volunteers and PAD patients before intervention for HbO2 (0.022±0.004 for patients vs. 0.036±0.006 for healthy, P=0.028), THb (0.064±0.008 vs. 0.045±0.006, P=0.009) and SO2 (0.56±0.04 vs. 0.484±0.042, P=0.041). Furthermore, patients with PAD show statistically significant differences for HbO2 before and after treatment (0.022±0.004 before, 0.036±0.006 after, P=0.029), while the change in THb before and after treatment is marginally significant or (0.045±0.006 before, 0.067±0.009 after, P=0.057) and the corresponding change in SO2 is non-significant (P=0.287). Finally, MSOT readouts show no significant difference for all measured parameters HbO2, Hb, THb and SO2 between the groups of the healthy volunteers and the PAD patients after the application of the interventional treatment. Figure 3E shows exemplary images of a patient with PAD preoperatively and postoperatively for all MSOT-extracted parameters.
Discussion
In this pilot study, we employed clinical MSOT hybridized with conventional ultrasound as a novel imaging technique for PAD assessment and therapy evaluation. PAD is a serious condition, that carries a substantial morbidity and is additionally considered to be an indicator for coronary artery disease and stroke (24,25). PAD affects more than 200 million people worldwide (26), with significant restraints in life quality and huge healthcare costs. Currently established techniques for PAD imaging focus on the macrovascular pathology of atherosclerotic lesions, without providing any functional or metabolic assessment (e.g., perfusion or oxygenation) of the affected muscle. A timely assessment of 41,882 hospitalized PAD patients revealed that, especially in advanced stages of disease, a lack of appropriate diagnostics limits the use of revascularization procedure resulting in an unacceptable high rate of amputations (1). In detail, the study reports an amputation rate of 18.7% and an amputation risk of 59.6% within 1 year in patients with Rutherford stage 6 PAD.
Herein, we demonstrate that clinical hand-held MSOT can be used for extracting image-based biomarkers of muscle hemodynamics and oxygenation in patients with PAD. First, we show that MSOT is capable of detecting and quantifying the arterial occlusion-associated changes within muscle over time in real-time. Moreover, by sensing the reactive post-occlusive hyperemic response, MSOT may be developed as a novel tool for assessing microvascular endothelial dysfunction. MSOT can also produce high resolution maps of Hb and HbO2 deep in skeletal muscle over time, surpassing the information accessible using NIRS. Thus, MSOT could enable the definition of new image-based biomarkers as possible indicators of microvascular disease that appear before the manifestation of anatomical macrovascular changes in the affected limb, aiding the early identification of high-risk patients.
Second, we show that MSOT enables the differentiation between healthy volunteers and patients with PAD via direct and non-invasive muscle imaging without the need for exogenous contrast agents by extraction of perfusion- and oxygenation-related functional tissue biomarkers. Thus, MSOT shows potential as a novel clinical tool for the detection of PAD and differentiation of other diseases causing similar clinical symptoms, such as spinal stenosis, peripheral neuropathy or deep vein thrombosis (27). This study focused on assessing muscle perfusion in healthy volunteers and PAD patients. However, skin perfusion might be also altered in PAD patients. Thus, further investigation of skin versus muscle perfusion in PAD might be promising in future studies. In this context, although MSOT/US has been used for skin and subcutaneous tissue assessment in diseases associated with microvascular dysfunction, dedicated optoacoustic techniques for skin imaging like raster-scan optoacoustic mesoscopy (RSOM) might be promising to assess skin perfusion in PAD patients.
Third, our results demonstrate that MSOT is able to visualize and quantify the effect of endovascular treatment on several functional parameters of muscle tissue perfusion. In our study patients with PAD showed, with an increase of about 1.6-fold, a significantly higher MSOT HbO2 value after than compared to before revascularization therapy. Although not fully comparable, the dimension of improvement is in line with other studies using for, e.g., thermal load tests with transcutaneous oxygen tension (tcPO2) measurement observing an improvement of 1.3-fold to up to 11-fold (28,29) or of NIRS measurements showing an HbO2 improvement of about 1.8-fold (30). Thus, MSOT-extracted parameters may well be used for the objective quantification of different therapeutic effects.
Nevertheless, even if MSOT offers rich functional information via hemoglobin-based molecular contrast, this information should be evaluated with conventional imaging techniques. For example, our clinical hybrid MSOT-US scanner facilitates such a combination by overlaying functional MSOT data with anatomical information from US. Such system configurations open up excellent possibilities for revealing correlations between MSOT and other imaging modalities in order to promote its clinical translation as a hand-held molecular imaging technique.
MSOT technology does not come without limitations. Even with unprecedented penetration depths compared to purely optical methods, the depth achieved (3–4 cm) is still poor compared to traditional modalities used for PAD imaging in everyday clinical practice (US, MRA, CTA). The development of new illumination regimes or schemes for light attenuation with depth is expected to address these limitations in the future. Furthermore, the spectral unmixing step is conducted on a per-pixel basis and is thus vulnerable to motion artifacts, which may be further emphasized by noise due to vessel reactivity. For this reason, several motion correction algorithms have been recently introduced (31,32). Further, a recent study showed that MSOT provided consistent and reproducible functional soft tissue characterization, independent on the investigating personnel (33). However, robustness and reproducibility of MSOT should be further investigated in future studies.
Notably, our pilot clinical study is limited by a low patient number. Future studies are expected to test the specificity, sensitivity and accuracy of MSOT technology for diagnosis or therapy monitoring in PAD. Nevertheless, current results demonstrate the capability of MSOT as a portable, non-invasive, label- and radiation-free molecular imaging modality for quantitative imaging of functional biomarkers in PAD with novel implications for both basic research on PAD pathophysiology and clinical molecular imaging.
Conclusions
The potential use of MSOT for PAD opens up new possibilities to investigate muscle oxygenation and metabolism in PAD, and to shift the paradigm from only imaging the cause of the disease (arterial stenosis), towards imaging the target organ (muscle) for more objective PAD diagnostics, treatment selection and therapy assessment.
Acknowledgments
We thank Dr. Jing Claussen (iThera Medical GmbH, Munich, Germany) and Dr. Niklas Markwardt (Institute for Biological and Medical Imaging, Helmholtz Center Munich, Neuherberg, Germany) for the invaluable technical support and fruitful discussions.
Funding: This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 694968 (PREMSOT), by the Helmholtz Zentrum München through Physician Scientists for Groundbreaking Projects, in part by the Helmholtz Association of German Research Center, through the Initiative and Networking Fund, i3 (ExNet-0022-Phase2-3) and was supported by the DZHK (German Centre for Cardiovascular Research).
Footnote
Data Sharing Statement: Available at http://dx.doi.org/10.21037/atm-20-3321
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at http://dx.doi.org/10.21037/atm-20-3321). Prof. VN reports to be a stock owner of iThera Medical GmbH, Munich, Germany. The other authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). All participants gave written and informed consent prior to study enrolment. Healthy volunteers freely consented to participate in full accordance with the work safety regulations of the Helmholtz Center Munich (Neuherberg, Germany). The study was reviewed by the local ethics committee (Ethikkommission der Westfälischen Wilhelms-Universität Münster, Protocol No 2017-538-f-S).
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Reinecke H, Unrath M, Freisinger E, et al. Peripheral arterial disease and critical limb ischaemia: still poor outcomes and lack of guideline adherence. Eur Heart J 2015;36:932-8. [Crossref] [PubMed]
- Eberhardt RT, Coffman JD. Cardiovascular morbidity and mortality in peripheral arterial disease. Curr Drug Targets Cardiovasc Haematol Disord 2004;4:209-17. [Crossref] [PubMed]
- Aday AW, Kinlay S, Gerhard-Herman MD. Comparison of different exercise ankle pressure indices in the diagnosis of peripheral artery disease. Vasc Med 2018;23:541-8. [Crossref] [PubMed]
- Pollak AW, Norton PT, Kramer CM. Multimodality imaging of lower extremity peripheral arterial disease: current role and future directions. Circ Cardiovasc Imaging 2012;5:797-807. [Crossref] [PubMed]
- Harvey C. Ultrasound with microbubbles. Cancer Imaging 2015;15:O19. [Crossref]
- Suo S, Zhang L, Tang H, et al. Evaluation of skeletal muscle microvascular perfusion of lower extremities by cardiovascular magnetic resonance arterial spin labeling, blood oxygenation level-dependent, and intravoxel incoherent motion techniques. J Cardiovasc Magn Reson 2018;20:18. [Crossref] [PubMed]
- Jones S, Chiesa ST, Chaturvedi N, et al. Recent developments in near-infrared spectroscopy (NIRS) for the assessment of local skeletal muscle microvascular function and capacity to utilise oxygen. Artery Res 2016;16:25-33. [Crossref] [PubMed]
- Keller E, Nadler A, Alkadhi H, et al. Noninvasive measurement of regional cerebral blood flow and regional cerebral blood volume by near-infrared spectroscopy and indocyanine green dye dilution. Neuroimage 2003;20:828-39. [Crossref] [PubMed]
- Khalil MA, Kim HK, Hoi JW, et al. Detection of Peripheral Arterial Disease Within the Foot Using Vascular Optical Tomographic Imaging: A Clinical Pilot Study. Eur J Vasc Endovasc Surg 2015;49:83-9. [Crossref] [PubMed]
- Karlas A, Reber J, Diot G, et al. Flow-mediated dilatation test using optoacoustic imaging: a proof-of-concept. Biomedical Optics Express 2017;8:3395-403. [Crossref] [PubMed]
- Yang H, Jüstel D, Prakash J, et al. Soft ultrasound priors in optoacoustic reconstruction: Improving clinical vascular imaging. Photoacoustics 2020;19:100172. [Crossref] [PubMed]
- Karlas A, Fasoula NA, Paul-Yuan K, et al. Cardiovascular optoacoustics: From mice to men - A review. Photoacoustics 2019;14:19-30. [Crossref] [PubMed]
- Karlas A, Kallmayer M, Fasoula NA, et al. Multispectral Optoacoustic Tomography of Muscle Perfusion and Oxygenation under Arterial and Venous Occlusion: A Human Pilot Study. J Biophotonics 2020;13:e201960169. [Crossref] [PubMed]
- Taruttis A, Ntziachristos V. Advances in real-time multispectral optoacoustic imaging and its applications. Nat Photonics 2015;9:219-27. [Crossref]
- Masthoff M, Helfen A, Claussen J, et al. Use of Multispectral Optoacoustic Tomography to Diagnose Vascular Malformations. JAMA Dermatol 2018;154:1457-62. [Crossref] [PubMed]
- Masthoff M, Helfen A, Claussen J, et al. Multispectral optoacoustic tomography of systemic sclerosis. J Biophotonics 2018;11:e201800155. [Crossref] [PubMed]
- Reber J, Willershäuser M, Karlas A, et al. Non-invasive Measurement of Brown Fat Metabolism Based on Optoacoustic Imaging of Hemoglobin Gradients. Cell Metabolism 2018;27:689-701.e4. [Crossref] [PubMed]
- Roll W, Markwardt NA, Masthoff M, et al. Multispectral optoacoustic tomography of benign and malignant thyroid disorders: a pilot study. J Nucl Med 2019;60:1461-6. [Crossref] [PubMed]
- Diot G, Metz S, Noske A, et al. Multispectral Optoacoustic Tomography (MSOT) of Human Breast Cancer. Clin Cancer Res 2017;23:6912-22. [Crossref] [PubMed]
- Knieling F, Neufert C, Hartmann A, et al. Multispectral Optoacoustic Tomography for Assessment of Crohn's Disease Activity. N Engl J Med 2017;376:1292-4. [Crossref] [PubMed]
- Becker A, Masthoff M, Claussen J, et al. Multispectral optoacoustic tomography of the human breast: characterisation of healthy tissue and malignant lesions using a hybrid ultrasound-optoacoustic approach. Eur Radiol 2018;28:602-9. [Crossref] [PubMed]
- Ansi Z. American National Standard for Safe Use of Lasers. Laser Institute of America, Orlando; 2000.136.1.
- Rosenthal A, Ntziachristos V, Razansky D. Model-based optoacoustic inversion with arbitrary-shape detectors. Med Phys 2011;38:4285-95. [Crossref] [PubMed]
- Grenon SM, Vittinghoff E, Owens CD, et al. Peripheral artery disease and risk of cardiovascular events in patients with coronary artery disease: Insights from the Heart and Soul Study. Vasc Med 2013;18:176-84. [Crossref] [PubMed]
- Kolls BJ, Sapp S, Rockhold FW, et al. Stroke in Patients With Peripheral Artery Disease. Stroke 2019;50:1356-63. [Crossref] [PubMed]
- Shu J, Santulli G. Update on peripheral artery disease: Epidemiology and evidence-based facts. Atherosclerosis 2018;275:379-81. [Crossref] [PubMed]
- Weinberg I, Jaff MR. Nonatherosclerotic arterial disorders of the lower extremities. Circulation 2012;126:213-22. [Crossref] [PubMed]
- Katsui S, Inoue Y, Yamamoto Y, et al. In Patients with Severe Peripheral Arterial Disease, Revascularization-Induced Improvement in Lower Extremity Ischemia Can Be Detected by Laser Speckle Contrast Imaging of the Fluctuation in Blood Perfusion after Local Heating. Ann Vasc Surg 2018;48:67-74. [Crossref] [PubMed]
- Yamamoto Y, Inoue Y, Igari K, et al. Assessment of the Severity of Ischaemia and the Outcomes of Revascularisation in Peripheral Arterial Disease Patients Based on the Skin Microcirculatory Response to a Thermal Load Test. EJVES Short Rep 2019;42:21-5. [Crossref] [PubMed]
- Manfredini F, Lamberti N, Ficarra V, et al. Biomarkers of Muscle Metabolism in Peripheral Artery Disease: A Dynamic NIRS-Assisted Study to Detect Adaptations Following Revascularization and Exercise Training. Diagnostics (Basel) 2020;10:312. [Crossref] [PubMed]
- Ulrich L, Held KG, Jaeger M, et al. Reliability assessment for blood oxygen saturation levels measured with optoacoustic imaging. J Biomed Opt 2020;25:1-15. [Crossref] [PubMed]
- Jeng GS, Li ML, Kim M, et al. Real-time spectroscopic photoacoustic/ultrasound (PAUS) scanning with simultaneous fluence compensation and motion correction for quantitative molecular imaging. bioRxiv 2019. doi: 10.1101/2019.12.20.885251. [Crossref]
- Helfen A, Masthoff M, Claussen J, et al. Multispectral Optoacoustic Tomography: Intra- and Interobserver Variability Using a Clinical Hybrid Approach. J Clin Med 2019;8:63. [Crossref] [PubMed]