Optical tuning of copolymer-in-oil tissue-mimicking materials for multispectral photoacoustic imaging.

OBJECTIVE: The availability of tissue-mimicking materials (TMMs) for manufacturing high-quality phantoms is crucial for standardization, evaluating novel quantitative approaches, and clinically translating new imaging modalities, such as photoacoustic imaging (PAI). Recently, a gel comprising the copolymer styrene-ethylene/butylene-styrene (SEBS) in mineral oil has shown significant potential as TMM due to its optical and acoustic properties akin to soft tissue. We propose using artists' oil-based inks dissolved and diluted in balsam turpentine to tune the optical properties. APPROACH: A TMM was fabricated by mixing a SEBS copolymer and mineral oil, supplemented with additives to tune its optical absorption and scattering properties independently. A systematic investigation of the tuning accuracies and relationships between concentrations of oil-based pigments and optical absorption properties of the TMM across visible and near-infrared wavelengths using collimated transmission spectroscopy was conducted. The photoacoustic spectrum of various oil-based inks was studied to analyze the effect of increasing concentration and depth. MAIN RESULTS: Artists' Oil-based inks dissolved in turpentine proved effective as additives to tune the optical absorption properties of mineral oil SEBS-gel with high accuracy. The TMMs demonstrated long-term stability and suitability for producing phantoms with desired optical absorption properties for PAI studies. SIGNIFICANCE: The findings, including tuning of optical absorption and spectral shape, suggest that this TMM facilitates the development of more sophisticated phantoms of arbitrary shapes. This approach holds promise for advancing the development of PAI, including investigation of the spectral coloring effect. In addition, it can potentially aid in the development and clinical translation of ultrasound optical tomography.

A Computer Program for Assessing Histoanatomical Morphometrics in Ultra-High-Frequency Ultrasound Images of the Bowel Wall in Children: Development and Inter-Observer Variability.

Ultra-high-frequency ultrasound (UHFUS) has a reported potential to differentiate between aganglionic and ganglionic bowel wall, referred to as histoanatomical differences. A good correlation between histoanatomy and UHFUS of the bowel wall has been proven. In order to perform more precise and objective histoanatomical morphometrics, the main research objective of this study was to develop a computer program for the assessment and automatic calculation of the histoanatomical morphometrics of the bowel wall in UHFUS images. A computer program for UHFUS diagnostics was developed and presented. A user interface was developed in close collaboration between pediatric surgeons and biomedical engineers, to enable interaction with UHFUS images. Images from ex vivo bowel wall samples of 23 children with recto-sigmoid Hirschsprung's disease were inserted. The program calculated both thickness and amplitudes (image whiteness) within different histoanatomical bowel wall layers. Two observers assessed the images using the program and the inter-observer variability was evaluated. There was an excellent agreement between observers, with an intraclass correlation coefficient range of 0.970-0.998. Bland-Altman plots showed flat and narrow distributions. The mean differences ranged from 0.005 to 0.016 mm in thickness and 0 to 0.7 in amplitude units, corresponding to 1.1-3.6% and 0.0-0.8% from the overall mean. The computer program enables and ensures objective, accurate and time-efficient measurements of histoanatomical thicknesses and amplitudes in UHFUS images of the bowel wall. The program can potentially be used for several bowel wall conditions, accelerating research within UHFUS diagnostics.

The impact of geometry, intramural friction, and pressure on the antegrade longitudinal motion of the arterial wall: A phantom and finite element study.

Longitudinal motion of the carotid arterial wall, as measured with ultrasound, has shown promise as an indicator of vascular health. The underlying mechanisms are however not fully understood. We have found, in in vivo studies, that blood pressure has a strong relation to the antegrade longitudinal displacement in early systole. Further, we have identified that a tapered geometry and the intramural friction in-between two parts of a vessel wall influence the longitudinal displacement. We therefore studied the interaction between pressure, vessel geometry and intramural friction, tapered and straight ultrasound phantoms in a paralleled hydraulic bench study and corresponding numerical models. Profound antegrade longitudinal motion was induced in the innermost part of both tapered phantoms and the numerical models, but to a lesser extent when intramural friction was increased in the simulations. Strong correlations (R = 0.82-0.96; p < 1e-3; k = 9.3-14 µm/mmHg) between longitudinal displacement and pulse pressure were found in six of seven regions of interest in tapered phantoms. The motion of the straight phantom and the corresponding numerical model was smaller, on average zero or close to zero. This study demonstrates that tapering of the lumen, low intramural friction, and pressure might be important conducive features to the antegrade longitudinal motion of the arterial wall in vivo.

Ultra-High Frequency Ultrasound Imaging of Bowel Wall in Hirschsprung's Disease-Correlation and Agreement Analyses of Histoanatomy.

Hirschsprung's disease (HD) is characterized by aganglionosis in the bowel wall, requiring resection. Ultra-high frequency ultrasound (UHFUS) imaging of the bowel wall has been suggested to be an instantaneous method of deciding resection length. The aim of this study was to validate UHFUS imaging of the bowel wall in children with HD by exploring the correlation and systematic differences between UHFUS and histopathology. Resected fresh bowel specimens of children 0-1 years old, operated on for rectosigmoid aganglionosis at a national HD center 2018-2021, were examined ex vivo with UHFUS center frequency 50 MHz. Aganglionosis and ganglionosis were confirmed by histopathological staining and immunohistochemistry. Histoanatomical layers of bowel wall in histopathological and UHFUS images, respectively, were outlined using MATLAB programs. Both histopathological and UHFUS images were available for 19 aganglionic and 18 ganglionic specimens. The thickness of muscularis interna correlated positively between histopathology and UHFUS in both aganglionosis (R = 0.651, p = 0.003) and ganglionosis (R = 0.534, p = 0.023). The muscularis interna was systematically thicker in histopathology than in UHFUS images in both aganglionosis (0.499 vs. 0.309 mm; p < 0.001) and ganglionosis (0.644 versus 0.556 mm; p = 0.003). Significant correlations and systematic differences between histopathological and UHFUS images support the hypothesis that UHFUS reproduces the histoanatomy of the bowel wall in HD accurately.

Design of a Pediatric Rectal Ultrasound Probe Intended for Ultra-High Frequency Ultrasound Diagnostics.

It has been shown that ultra-high frequency (UHF) ultrasound applied to the external bowel wall can delineate the histo-anatomic layers in detail and distinguish normal bowel from aganglionosis. This would potentially reduce or lessen the need for biopsies that are currently mandatory for the diagnosis of Hirschsprung's disease. However, to our knowledge, no suitable rectal probes for such a use are on the market. The aim was to define the specifications of an UHF transrectal ultrasound probe (50 MHz center frequency) suitable for use in infants. Probe requirements according to patient anatomy, clinicians' requests, and biomedical engineering UHF prerequisites were collected within an expert group. Suitable probes on the market and in clinical use were reviewed. The requirements were transferred into the sketching of potential UHF ultrasound transrectal probes followed by their 3D prototype printing. Two prototypes were created and tested by five pediatric surgeons. The larger and straight 8 mm head and shaft probe was preferred as it facilitated stability, ease of anal insertion, and possible UHF technique including 128 piezoelectric elements in a linear array. We hereby present the procedure and considerations behind the development of a proposed new UHF transrectal pediatric probe. Such a device can open new possibilities for the diagnostics of pediatric anorectal conditions.

Response of the carotid artery longitudinal motion to submaximal physical activity in healthy humans-Marked changes already at low workload.

The longitudinal motion of the arterial wall, that is, the displacement of the arterial wall along the artery, parallel to blood flow, is still largely unexplored. The magnitude and nature of putative changes in longitudinal motion of the arterial wall in response to physical activity in humans remain unknown. The aim of this study was therefore to study the longitudinal motion of the carotid artery wall during physical activity in healthy humans. Using in-house developed non-invasive ultrasonic methods, the longitudinal motion of the intima-media complex and the diameter changes of the right common carotid artery (CCA) in 40 healthy volunteers (20 volunteers aged 22-35 years; 20 volunteers aged 55-68 years) were assessed at rest and during submaximal supine bicycle exercise. In a subset of the subjects (n = 18) also intramural shear strain were analyzed. The longitudinal motion of the intima-media complex underwent marked changes in response to physical activity, already at low workload; with most evident a marked increase of the first antegrade displacement (p < 0.001) in early systole. Likewise, the corresponding shear strain also increased significantly (p = 0.004). The increase in longitudinal motion showed significant correlation to increase in blood pressure, but not to blood flow velocity or wall shear stress. In conclusion, physical activity markedly influences the longitudinal motion of the carotid artery wall in healthy humans already at low load. A possible "cushioning" function as well as possible implications for the function of the vasa vasorum, endothelium, and smooth muscle cells and extracellular matrix of the media, are discussed.

Histopathological dimensions differ between aganglionic and ganglionic bowel wall in children with Hirschsprung's disease.

BACKGROUND: In the validation of new imaging technology for children with Hirschsprung's disease (HSCR), basic anatomical parameters of the bowel wall must be established specifically for this patient group. AIM: To explore differences in histoanatomical layers of bowel wall, comparing ganglionic and aganglionic bowel walls, and to examine if the bowel wall thickness is linked to patient weight. METHODS: This was an observational study of bowel specimens from children weighing 0-10 kg, operated on consecutively during 2018-2020. Ganglionic and aganglionic bowel walls were measured in digitalized microscopy images from 10 sites per trans-sectional specimen and compared regarding the thickness of their histoanatomical layers. RESULTS: Bowel walls were measured in 21 children. Full bowel wall thickness did not differ between aganglionic and ganglionic bowel (2.20 vs 2.04; p = 0.802) while weight at surgery correlated positively with both ganglionic and aganglionic bowel wall thickness (r = 0.688 and 0.849, respectively), and age at surgery with ganglionic bowel wall thickness (r = 0.517). In aganglionic segments, the muscularis externa layer was thicker compared to that in ganglionosis (0.45 vs 0.31 mm, p = 0.012) whereas the muscularis interna was thinner (0.45 vs 0.62 mm, p < 0.001). A diagnostic index was identified whereby a lower ratio of muscularis interna/externa thickness followed by a thinner muscularis interna differed between aganglionic and ganglionic bowel in all specimens. CONCLUSION: Thicknesses of the bowel wall's muscle layers differ between aganglionic and ganglionic bowel walls in children with HSCR. These findings support a diagnostic index that could be validated for transfer to instant diagnostic imaging techniques. LEVEL OF EVIDENCE: Diagnostic: 3.

Automatic threshold selection algorithm to distinguish a tissue chromophore from the background in photoacoustic imaging.

The adaptive matched filter (AMF) is a method widely used in spectral unmixing to classify different tissue chromophores in photoacoustic images. However, a threshold needs to be applied to the AMF detection image to distinguish the desired tissue chromophores from the background. In this study, we propose an automatic threshold selection (ATS) algorithm capable of differentiating a target from the background, based on the features of the AMF detection image. The mean difference between the estimated thickness, using the ATS algorithm, and the known values was 0.17 SD (0.24) mm for the phantom inclusions and -0.05 SD (0.21) mm for the tissue samples of malignant melanoma. The evaluation shows that the thickness and the width of the phantom inclusions and the tumors can be estimated using AMF in an automatic way after applying the ATS algorithm.

Photoacoustic imaging for the monitoring of local changes in oxygen saturation following an adrenaline injection in human forearm skin.

Clinical monitoring of blood oxygen saturation (sO2) is traditionally performed using optical techniques, such as pulse oximetry and diffuse reflectance spectroscopy (DRS), which lack spatial resolution. Photoacoustic imaging (PAI) is a rapidly developing biomedical imaging technique that is superior to previous techniques in that it combines optical excitation and acoustic detection, providing a map of chromophore distribution in the tissue. Hitherto, PAI has primarily been used in preclinical studies, and only a few studies have been performed in patients. Its ability to measure sO2 with spatial resolution during local vasoconstriction after adrenaline injection has not yet been investigated. Using PAI and spectral unmixing we characterize the heterogeneous change in sO2 after injecting a local anesthetic containing adrenaline into the dermis on the forearm of seven healthy subjects. In comparison to results obtained using DRS, we highlight contrasting results obtained between the two methods arising due to the so-called 'window effect' caused by a reduced blood flow in the superficial vascular plexus. The results demonstrate the importance of spatially resolving sO2 and the ability of PAI to assess the tissue composition in different layers of the skin.

Comparison of photoacoustic imaging and histopathological examination in determining the dimensions of 52 human melanomas and nevi ex vivo.

Surgical excision followed by histopathological examination is the gold standard for the diagnosis and staging of melanoma. Reoperations and unnecessary removal of healthy tissue could be reduced if non-invasive imaging techniques were available for presurgical tumor delineation. However, no technique has gained widespread clinical use to date due to shallow imaging depth or the absence of functional imaging capability. Photoacoustic (PA) imaging is a novel technology that combines the strengths of optical and ultrasound imaging to reveal the molecular composition of tissue at high resolution. Encouraging results have been obtained from previous animal and human studies on melanoma, but there is still a lack of clinical data. This is the largest study of its kind to date, including 52 melanomas and nevi. 3D multiwavelength PA scanning was performed ex vivo, using 59 excitation wavelengths from 680 nm to 970 nm. Spectral unmixing over this broad wavelength range, accounting for the absorption of several tissue chromophores, provided excellent contrast between healthy tissue and tumor. Combining the results of spectral analysis with spatially resolved information provided a map of the tumor borders in greater detail than previously reported. The tumor dimensions determined with PA imaging were strongly correlated with those determined by histopathological examination for both melanomas and nevi.

Acousto-optic interaction strengths in optically scattering media using high pressure acoustic pulses.

Ultrasound optical tomography (UOT) is a developing medical imaging technique with the potential to noninvasively image tissue oxygenation at depths of several centimeters in human tissue. To accurately model the UOT imaging, it is necessary the calculate the signal produced by the interaction between ultrasound and light in the scattering medium. In this paper we present a rigorous description for modeling this process for ultrasound pulses in the non-linear regime with peak pressures ranging up to the medical safety limit. Simulation results based on the presented model agree well with measurements performed with fully characterized ultrasound pulses. Our results also indicate that the UOT modeling process can be accurately simplified by disregarding the acoustically induced movement of scatterers. Our results suggest that the explored model and its software implementation can be used as a virtual lab to aid future development of pulses and UOT imaging algorithms.

Regional motion correction for in vivo photoacoustic imaging in humans using interleaved ultrasound images.

In translation from preclinical to clinical studies using photoacoustic imaging, motion artifacts represent a major issue. In this study the feasibility of an in-house algorithm, referred to as intensity phase tracking (IPT), for regional motion correction of in vivo human photoacoustic (PA) images was demonstrated. The algorithm converts intensity to phase-information and performs 2D phase-tracking on interleaved ultrasound images. The radial artery in eight healthy volunteers was imaged using an ultra-high frequency photoacoustic system. PA images were motion corrected and evaluated based on PA image similarities. Both controlled measurements using a computerized stepping motor and free-hand measurements were evaluated. The results of the controlled measurements show that the tracking corresponded to 97 ± 6% of the actual movement. Overall, the mean square error between PA images decreased by 52 ± 15% and by 43 ± 19% when correcting for controlled- and free-hand induced motions, respectively. The results show that the proposed algorithm could be used for motion correction in photoacoustic imaging in humans.

Photoacoustic imaging of the spatial distribution of oxygen saturation in an ischemia-reperfusion model in humans.

Photoacoustic imaging (PAI) is a novel hybrid imaging technique that combines the advantages of optical and ultrasound imaging to produce hyperspectral images of the tissue. The feasibility of measuring oxygen saturation (sO2) with PAI has been demonstrated pre-clinically, but has limited use in humans under conditions of ischemia and reperfusion. As an important step towards making PAI clinically available, we present a study in which PAI was used to estimate the spatial distribution of sO2 in vivo during and after occlusion of the finger of eight healthy volunteers. The results were compared with a commercial oxygen saturation monitor based on diffuse reflectance spectroscopy. We here describe the capability of PAI to provide spatially resolved picture of the evolution of sO2 during ischemia following vascular occlusion of a finger, demonstrating the clinical viability of PAI as a non-invasive diagnostic tool for diseases indicated by impaired microvascularization.

Comparison of the multi-phasic longitudinal displacement of the left and right common carotid artery in healthy humans.

BACKGROUND: During the cardiac cycle, there is a multi-phasic bidirectional longitudinal movement (LMov) of the intima-media complex of large arteries, i.e. along the arteries. On the left side the common carotid artery (CCA) arises directly from the aortic arc, whereas on the right side the CCA originate from the innominate artery. AIM: The aim of this study was to compare LMov of the left and right CCA of healthy subjects to investigate whether the difference in anatomy is of importance for LMov. MATERIAL AND METHODS: The CCA's of 93 healthy subjects were investigated using in-house developed ultrasound methods. RESULTS: Although the basic pattern were the same in the majority of subjects, several phases of LMov were significantly larger on the left side (the first retrograde phase, p = 0.0006; the second antegrade, "returning" phase, p < 0.00001; and the rapid retrograde phase of movement at the end of the cardiac cycle, p < 0.000001). In contrast, no significant side-difference in the amplitude of the first antegrade movement was seen. The maximal (peak-to-peak) LMov was significantly larger on the left side (p = 0.002). DISCUSSION AND CONCLUSION: The side-differences found in LMov may be related to the difference in anatomy, including possible difference in distance to the heart and especially the presence of an extra bifurcation on the right side. Our data provide an important base for the further study of the relation between LMov and cardiovascular risk factors and atherosclerosis.

Ultra high frequency ultrasonography to distinguish ganglionic from aganglionic bowel wall in Hirschsprung disease: A first report.

BACKGROUND/PURPOSE: In Hirschsprung disease (HD) surgery, confirming ganglionic bowel is essential. A faster diagnostic method than the current frozen biopsy is desirable. This study investigated whether aganglionic and ganglionic intestinal wall can be distinguished from each other by ultra high frequency ultrasound (UHF ultrasound). METHODS: In an HD center during 2019, intestinal walls of recto-sigmoid specimens from HD patients were examined ex vivo with a 70 MHz UHF ultrasound transducer. Data from four sites were described. Histopathologic analysis was compared to the ultrasonography outcome at each site. Each patient's specimen served as its own control. RESULTS: 11 resected recto-sigmoid specimens (median 20 cm long [range 6.5-33]) with transition zones of 5 cm (2-11 cm) were taken from children aged 22 days (13-48) weighing 3668 g (3500-5508); 44 key sites were analyzed. There was full concordance for 42/44 (95%) key sites and 10 of 11 (91%) specimens. The specimen with discordance of two key sites contained a segment of aganglionosis (3 cm) and a transition zone (1 cm): the site discordance was limited to the transition zone ends. CONCLUSIONS: This first report on UHF ultrasound in recto-sigmoid HD shows promising results in identifying aganglionosis, transition zones and ganglionic bowel. Further in vivo studies are required.

Scaled reassigned spectrograms applied to linear transducer signals.

This study evaluates the applicability of scaled reassigned spectrograms (ReSTS) on ultrasound radio frequency data obtained with a clinical linear array ultrasound transducer. The ReSTS's ability to resolve axially closely spaced objects in a phantom is compared to the classical cross-correlation method with respect to the ability to resolve closely spaced objects as individual reflectors using ultrasound pulses with different lengths. The results show that the axial resolution achieved with the ReSTS was superior to the cross-correlation method when the reflected pulses from two objects overlap. A novel B-mode imaging method, facilitating higher image resolution for distinct reflectors, is proposed.

Photoacoustic imaging for non-invasive examination of the healthy temporal artery - systematic evaluation of visual function in healthy subjects.

PURPOSE: Photoacoustic (PA) imaging has the potential to become a non-invasive diagnostic tool for giant cell arteritis, as shown in pilot experiments on seven patients undergoing surgery. Here, we present a detailed evaluation of the safety regarding visual function and patient tolerability in healthy subjects, and define the spectral signature in the healthy temporal artery. METHODS: Photoacoustic scanning of the temporal artery was performed in 12 healthy subjects using 59 wavelengths (from 680 nm to 970 nm). Visual function was tested before and after the examination. The subjects' experience of the examination was rated on a 0-100 VAS scale. Two- and three-dimensional PA images were generated from the spectra obtained from the artery. RESULTS: Photoacoustic imaging did not affect the best corrected visual acuity, colour vision (tested with Sahlgren's Saturation Test or the Ishihara colour vision test) or the visual field. The level of discomfort was low, and only little heat and light sensation were reported. The spectral signature of the artery wall could be clearly differentiated from those of the subcutaneous tissue and skin. Spectral unmixing provided visualization of the chromophore distribution and overall architecture of the artery. CONCLUSIONS: Photoacoustic imaging of the temporal artery is well tolerated and can be performed without any risk to visual function, including the function of the retina and the optic nerve. The spectral signature of the temporal artery is specific, which is promising for future method development.

Photoacoustic imaging for three-dimensional visualization and delineation of basal cell carcinoma in patients.

BACKGROUND: Photoacoustic (PA) imaging is an emerging non-invasive biomedical imaging modality that could potentially be used to determine the borders of basal cell carcinomas (BCC) preoperatively in order to reduce the need for repeated surgery. METHODS: Two- and three-dimensional PA images were obtained by scanning BCCs using 59 wavelengths in the range 680-970 nm. Spectral unmixing was performed to visualize the tumor tissue distribution. Spectral signatures from 38 BCCs and healthy tissue were compared ex vivo. RESULTS AND DISCUSSION: The PA spectra could be used to differentiate between BCC and healthy tissue ex vivo (p < 0.05). Spectral unmixing provided visualization of the overall architecture of the lesion and its border. CONCLUSION: PA imaging can be used to differentiate between BCC and healthy tissue and can potentially be used to delineate tumors prior to surgical excision.

Unique spectral signature of human cutaneous squamous cell carcinoma by photoacoustic imaging.

Cutaneous squamous cell carcinoma (cSCC) is a common skin cancer with metastatic potential. To reduce reoperations due to nonradical excision, there is a need to develop a technique for identification of tumor margins preoperatively. Photoacoustic (PA) imaging is a novel imaging technology that combines the strengths of laser optics and ultrasound. Our aim was to determine the spectral signature of cSCC using PA imaging and to use this signature to visualize tumor architecture and borders. Two-dimensional PA images of 33 cSCCs and surrounding healthy skin were acquired ex vivo, using 59 excitation wavelengths from 680 to 970 nm. The spectral response of the cSCCs was compared to healthy tissue, and the difference was found to be greatest at wavelengths in the range 765 to 960 nm (P < .05). Three-dimensional PA images were constructed from spectra obtained in the y-z plane using a linear stepper motor moving along the x-plane. Spectral unmixing was then performed which provided a clear three-dimensional view of the distribution of tumor masses and their borders.

Characterization and modeling of acousto-optic signal strengths in highly scattering media.

Ultrasound optical tomography (UOT) is an imaging technique based on the acousto-optic effect that can perform optical imaging with ultrasound resolution inside turbid media, and is thus interesting for biomedical applications, e.g. for assessing tissue blood oxygenation. In this paper, we present near background free measurements of UOT signal strengths using slow light filter signal detection. We carefully analyze each part of our experimental setup and match measured signal strengths with calculations based on diffusion theory. This agreement between experiment and theory allows us to assert the deep tissue imaging potential of ∼ 5 cm for UOT of real human tissues predicted by previous theoretical studies [Biomed. Opt. Express8, 4523 (2017)] with greater confidence, and indicate that future theoretical analysis of optimized UOT systems can be expected to be reliable.