Material Properties of Human Ocular Tissue at 7-µm Resolution

Biomechanical changes occur in myopic choroidal stroma and mirror those in the adjacent sclera

Retina-derived growth signals relayed from the choroid to the sclera cause remodeling of the extracellular scleral matrix, resulting in myopic ocular elongation. However, to the best of our knowledge, no studies have assessed changes in choroidal stromal biomechanical properties during myopia progression. Here we utilized 7 µm-resolution scanning acoustic microscopy (SAM) to assess biomechanical properties (bulk modulus (K) and mass density (rho)) of choroidal stroma from guinea pig eyes with form-deprivation (FD) induced myopia. The choroidal stroma had considerable intrinsic strength arising from its biomechanical properties and these were differentially affected by myopia in central and peripheral regions. Choroidal stromal biomechanical values were also highly correlated with those in adjacent scleral regions, and the choroidal stromal-scleral association was stronger in myopic eyes. Biomechanical changes observed in the choroidal stroma of myopic eyes were mirrored to those observed in the adjacent sclera. These findings suggest that choroidal stromal remodeling may accompany myopia and open the door to the source of the signals that cause scleral remodeling in myopia.

Scanning Laser Ophthalmoscopy Demonstrates Disc and Peripapillary Strain During Horizontal Eye Rotation in Adults

PURPOSE

We used automated image analysis of scanning laser ophthalmoscopy (SLO) to investigate mechanical strains imposed on disc, and retinal and choroidal vessels during horizontal duction in adults.

DESIGN

Deep learning analysis of optical images.

METHODS

The peripapillary region was imaged by SLO in central gaze, and 35° abduction and adduction, in younger and older healthy adults. Automated image registration was followed by deep learning-based optical flow analysis to track determine local tissue deformations quantified as horizontal, vertical, and shear strain maps relative to central gaze. Choroidal vessel displacements were observed when fundus pigment was light.

RESULTS

Strains in the retina and disc could be quantified in 22 younger (mean ± SEM, 26 ± 5 years) and 19 older (64 ± 10 years) healthy volunteers. Strains were predominantly horizontal and greater for adduction than for abduction. During adduction, maximum horizontal strain was tensile in the nasal hemi-disc, and declined progressively with distance from it. Strain in the temporal hemi-retina during adduction was minimal, except for compressive strain on the disc of older subjects. In abduction, horizontal strains were less and largely confined to the disc, greater in older subjects, and generally tensile. Vertical and shear strains were small. Nasal to the disc, choroidal vessels shifted nasally relative to overlying peripapillary retinal vessels.

CONCLUSIONS

Strain analysis during horizontal duction suggests that the optic nerve displaces the optic canal, choroid, and peripapillary sclera relative to the overlying disc and retina. This peripapillary shearing of the optic nerve relative to the choroid and sclera may be a driver of disc tilting and peripapillary atrophy.

Ultrasound-Mediated Ocular Drug Delivery: From Physics and Instrumentation to Future Directions

Drug delivery to the anterior and posterior segments of the eye is impeded by anatomical and physiological barriers. Increasingly, the bioeffects produced by ultrasound are being proven effective for mitigating the impact of these barriers on ocular drug delivery, though there does not appear to be a consensus on the most appropriate system configuration and operating parameters for this application. In this review, the fundamental aspects of ultrasound physics most pertinent to drug delivery are presented; the primary phenomena responsible for increased drug delivery efficacy under ultrasound sonication are discussed; an overview of common ocular drug administration routes and the associated ocular barriers is also given before reviewing the current state of the art of ultrasound-mediated ocular drug delivery and its potential future directions.

A systematic review of ultrasound-mediated drug delivery to the eye and critical insights to facilitate a timely path to the clinic

Ultrasound has long been identified as a promising, non-invasive modality for improving ocular drug delivery across a range of indications. Yet, with 20 years of learnings behind us, clinical translation remains limited. To help address this, and in accordance with PRISMA guidelines, the various mechanisms of ultrasound-mediated ocular drug delivery have been appraised, ranging from first principles to emergent applications spanning both ex vivo and in vivo models. The heterogeneity of study methods precluded meta-analysis, however an extensive characterisation of the included studies allowed for semi-quantitative and qualitative assessments. Methods: In this review, we reflected on study quality of reporting, and risk of bias (RoB) using the latest Animal Research: Reporting of In Vivo Experiments (ARRIVE 2.0) guidelines, alongside the Systematic Review Centre for Laboratory animal Experimentation (SYRCLE) RoB tools. Literature studies from 2002 to 2022 were initially characterised according to methods of ultrasound application, ultrasound parameters applied, animal models employed, as well as safety and efficacy assessments. This exercise contributed to developing a comprehensive understanding of the current state of play within ultrasound-mediated ocular drug delivery. The results were then synthesised and processed into a guide to aid future study design, with the goal of improving the reliability of data, and to support efficient and timely translation to the clinic. Results: Key attributes identified as hindering translation included: poor reporting quality and high RoB, skewed use of animals unrepresentative of the human eye, and the over reliance of reductionist safety assessments. Ex vivo modelling studies were often unable to have comprehensive safety assessments performed on them, which are imperative to determining treatment safety, and represent a pre-requisite for clinical translation. Conclusion: With the use of our synthesised guide, and a thorough understanding of the underlying physicochemical interactions between ultrasound and ocular biology provided herein, this review offers a firm foundation on which future studies should ideally be built, such that ultrasound-mediated ocular drug delivery can be translated from concept to the coalface where it can provide immense clinical benefit.

Advanced Topics in Quantitative Acoustic Microscopy

Quantitative acoustic microscopy (QAM) reconstructs two-dimensional (2D) maps of the acoustic properties of thin tissue sections. Using ultrahigh frequency transducers (≥ 100 MHz), unstained, micron-thick tissue sections affixed to glass are raster scanned to collect radiofrequency (RF) echo data and generate parametric maps with resolution approximately equal to the ultrasound wavelength. 2D maps of speed of sound, mass density, acoustic impedance, bulk modulus, and acoustic attenuation provide unique and quantitative information that is complementary to typical optical microscopy modalities. Consequently, many biomedical researchers have great interest in utilizing QAM instruments to investigate the acoustic and biomechanical properties of tissues at the micron scale. Unfortunately, current state-of-the-art QAM technology is costly, requires operation by a trained user, and is accompanied by substantial experimental challenges, many of which become more onerous as the transducer frequency is increased. In this chapter, typical QAM technology and standard image formation methods are reviewed. Then, novel experimental and signal processing approaches are presented with the specific goal of reducing QAM instrument costs and improving ease of use. These methods rely on modern techniques based on compressed sensing and sparsity-based deconvolution methods. Together, these approaches could serve as the basis of the next generation of QAM instruments that are affordable and provide high-resolution QAM images with turnkey solutions requiring nearly no training to operate.

Biomechanical changes in myopic sclera correlate with underlying changes in microstructure

Myopia alters the microstructural and biomechanical properties of the posterior sclera, which is characterized as a layered structure with potentially different inter-layer collagen fibril characteristics. Scanning acoustic microscopy (SAM) has been used to investigate how the micron-scale bulk mechanical properties of the posterior sclera are affected by myopia. Other investigators have employed second harmonic generation (SHG) imaging to characterize the collagen microstructure of tissues. In the present study, SAM and SHG imaging were used to investigate the existence of biomechanically-distinct scleral layers and identify relationships between mechanical properties and tissue microstructure in myopic guinea pig (GP) eyes. Diffusers were worn over the right eyes of six, 1-week-old GPs for one week to induce unilateral form-deprivation myopia. GPs were euthanized, enucleated, and eyes were cryosectioned. Twelve-micron-thick adjacent vertical cryosections were scanned with SAM or SHG. SAM maps of bulk modulus, mass density, and acoustic attenuation were estimated. A fiber-extraction algorithm applied to SHG images estimated collagen fiber length, width, straightness, alignment, and number density. Results revealed that the posterior sclera may exhibit biomechanically distinct layers that are affected differently in myopia. Specifically, a layered structure was observed in the mechanical-parameter maps of control eyes that was less apparent in myopic eyes. Collagen fibers in myopic eyes had smaller diameters and were more aligned. Myopia-associated biomechanical changes were most significant in the outermost and innermost scleral layers. SAM-measured mechanical parameters were correlated with collagen fiber microstructure, particularly fiber length, alignment, and number density, which may imply the biomechanical parameters estimated from SAM measurements are related to tissue microstructure. Interestingly, some changes were greatest in more-peripheral regions, suggesting interventions to strengthen the sclera may be effective away from the optic nerve and efficacy may be achieved best when intervention is applied to the outermost layer.

High-Frequency Ultrasound Activation of Perfluorocarbon Nanodroplets for Treatment of Glaucoma

Elevated intraocular pressure (IOP) is the most prevalent risk factor for initiation and progression of neurodegeneration in glaucoma. Ocular hypertension results from increased resistance to aqueous fluid outflow caused by reduced porosity and increased stiffness of tissues of the outflow pathway. Acoustic activation and resulting bioeffects of the perfluorocarbon (PFC) nanodroplets (NDs) introduced into the anterior chamber (AC) of the eye could potentially represent a treatment for glaucoma by increasing permeability in the aqueous outflow track. To evaluate the potential of NDs to enter the outflow track, 100-nm diameter perfluoropentane (PFP) NDs with a lipid shell were injected into the AC of ex vivo pig eyes and in vivo rat eyes. The NDs were activated and imaged with 18- and 28-MHz linear arrays to assess their location and diffusion. NDs in the AC could also be visualized using optical coherence tomography (OCT). Because of their higher density with respect to aqueous humor, some NDs settled into the iridocorneal angle where they entered the outflow pathway. After acoustic activation of the NDs at the highest acoustic pressure, small gas bubbles were observed in the AC. After two days, no acoustic activation events were visible in the AC of the rats and their eyes showed no evidence of inflammation.

Quantifying scattering from dense media using two-dimensional impedance maps

A better understanding of ultrasound scattering in a three-dimensional (3D) medium can provide more accurate methods for ultrasound tissue characterization. The possibility of using two-dimensional impedance maps (2DZMs) based on correlation coefficients has shown promise in the case of isotropic and sparse medium [Luchies and Oelze, J. Acoust. Soc. Am. 139, 1557-1564 (2016)]. The present study investigates the use of 2DZMs in order to quantify 3D scatterer properties of dense media from two-dimensional (2D) histological slices. Two 2DZM approaches were studied: one based on the correlation coefficient and the other based on the 2D Fourier transform of 2DZMs. Both 2DZM approaches consist in estimating the backscatter coefficient (BSC) from several 2DZMs, and then the resulting BSC was fit to the theoretical polydisperse structure factor model to yield 3D scatterer properties. Simulation studies were performed to evaluate the ability of both 2DZM approaches to quantify scattering of a 3D medium containing randomly distributed polydisperse spheres or monodisperse ellipsoids. Experimental studies were also performed using the histology photomicrographs obtained from HT29 cell pellet phantoms. Results demonstrate that the 2DZM Fourier transform-based approach was more suitable than the correlation coefficient-based approach for estimating scatterer properties when using a small number of 2DZMs.

Correction of phase rotation in pulse spectrum method for scanning acoustic microscopy and its application to measurements of cells

Scanning acoustic microscopy (SAM) can measure the mechanical properties, such as sound speed, thickness, and density, of biological tissues, by using the pulse spectrum method. However, the estimation method needs to be modified because of increases in the center frequency of acoustic transducers. In this paper, we proposed a new estimation method combining a time-of-flight method by Wiener filtering with the pulse spectrum method. First, an optimal control parameter β for Wiener filter was chosen based on a simulation by k-wave MATLAB toolbox. Setting the thickness of a layer to be 1.95 μm, a bias error between the estimated and true thickness was 0.0016% and the control parameter β was chosen to be 0.01 based on the simulated result and previous research. Next, the thickness of a film sample was measured by the time-of-flight method with Wiener filtering and was compared with an optically-measured thickness to confirm the estimation accuracy. Thickness was estimated to be 18.3 ± 0.025 μm at a center frequency of 120 MHz and agreed with the optically-measured thickness. Finally, the parameter n, the number of phase rotation in Gaussian plane, is calculated from the thickness and sound speed, and the pulse spectrum method with the correction of the parameter n is applied to the cellular measurements. Also, the mechanical properties estimated by the proposed method was compared with these by the conventional method.

Regional changes in the elastic properties of myopic Guinea pig sclera

Biomechanical changes in the sclera likely underlie the excessive eye elongation of axial myopia. We studied the biomechanical characteristics of myopic sclera at the microscopic level using scanning acoustic microscopy (SAM) with 7-μm in-plane resolution. Guinea pigs underwent form-deprivation (FD) in one eye from 4 to 12 days of age to induce myopia, and 12-μm-thick scleral cryosections were scanned using a custom-made SAM. Two-dimensional maps of the bulk modulus (K) and mass density (ρ) were derived from the SAM data using a frequency-domain approach. We assessed the effect on K and ρ exerted by: 1) level of induced myopia, 2) region (superior, inferior, nasal or temporal) and 3) eccentricity from the nerve using univariate and multivariate regression analyses. Induced myopia ranged between -3D and -9.3D (Mean intraocular difference of -6.2 ± 1.7D, N = 11). K decreased by 0.036 GPa for every 1.0 D increase in induced myopia across vertical sections (p < 0.001). Among induced myopia right eyes, K values in the inherently more myopic superior region were 0.088 GPa less than the inferior region (p = 0.002) and K in the proximal nasal region containing the central axis were 0.10 GPa less than temporal K (p = 0.036). K also increased 0.12 GPa for every 1 mm increase in superior vertical distance (p < 0.001), an effect that was blunted after 1 week of FD. Overall, trends for ρ were less apparent than for K. ρ values increased by 20.7 mg/cm3 for every 1.00 D increase in induced myopia across horizontal sections (p < 0.001), and were greatest in the region containing the central posterior pole. ρ values in the inherently more myopic superior region were 13.1 mg/cm3 greater than that found in inferior regions among control eyes (p = 0.002), and increased by 11.2 mg/cm3 for every 1 mm increase in vertical distance (p = 0.001). This peripheral increase in ρ was blunted after 1 week of FD. Scleral material properties vary depending on the location in the sclera and the level of induced myopia. Bulk modulus was most reduced in the most myopic regions (both induced myopia and inherent regional myopia), and suggests that FD causes microscopic local decreases in sclera stiffness, while scleral mass density was most increased in the most myopic regions.

Autoregressive Signal Processing Applied to High-Frequency Acoustic Microscopy of Soft Tissues

Quantitative acoustic microscopy (QAM) at frequencies exceeding 100 MHz has become an established imaging tool to depict acoustical and mechanical properties of soft biological tissues at microscopic resolutions. In this study, we investigate a novel autoregressive (AR) model to improve signal processing and parameter estimation and to test its applicability to QAM. The performance of the AR model for estimating acoustical parameters of soft tissues (i.e., acoustic impedance, speed of sound, and attenuation) was compared to the performance of the Hozumi model using simulated ultrasonic QAM signals and using experimentally measured signals from thin (i.e., 12 and ) sections of human lymph-node and pig-cornea tissue specimens. Results showed that the AR and Hozumi methods performed equally well (i.e., produced an estimation error of 0) in signals with low, linear attenuation in the tissue and high impedance contrast between the tissue and the coupling medium. However, the AR model outperformed the Hozumi model in estimation accuracy and stability (i.e., parameter error variation and number of outliers) in cases of 1) thin tissue-sample thickness and high tissue-sample speed of sound, 2) small impedance contrast between the tissue sample and the coupling medium, 3) high attenuation in the tissue sample, and 4) nonlinear attenuation in the tissue sample. Furthermore, the AR model allows estimating the exponent of nonlinear attenuation. The results of this study suggest that the AR model approach can improve current QAM by providing more reliable, quantitative, tissue-property estimates and also provides additional values of parameters related to nonlinear attenuation.

Improved High-Frequency Ultrasound Corneal Biometric Accuracy by Micrometer-Resolution Acoustic-Property Maps of the Cornea

Purpose

Mapping of epithelial thickness (ET) is useful for detection of keratoconus, a disease characterized by corneal thinning and bulging in which epithelial thinning occurs over the apex. In prior clinical studies, optical coherence tomography (OCT) measurements of ET were systematically thinner than those obtained by 40-MHz high-frequency ultrasound (HFU) where a constant speed of sound (c) of 1636 m/s was used for all corneal layers. The purpose of this work was to study the acoustic properties, that is, c, acoustic impedance (Z), and attenuation (α) of the corneal epithelium and stroma independently using a scanning acoustic microscope (SAM) to investigate the discrepancy between OCT and HFU estimates of ET.

Methods

Twelve unfixed pig corneas were snap-frozen and 6-μm sections were scanned using a custom-built SAM with an F-1.08, 500-MHz transducer and a 264-MHz bandwidth. Two-dimensional maps of c, Z, and α with a spatial resolution of 4 μm were derived.

Results

SAM showed that the value of c in the epithelium (i.e., 1548 ± 18 m/s) is substantially lower than the value of c in the stroma (i.e., 1686 ± 33 m/s).

Conclusion

SAM results demonstrated that the assumption of a constant value of c for all corneal layers is incorrect and explains the prior discrepancy between OCT and HFU ET determinations.

Translational Relevance

The findings of this study have important implications for HFU-based ET measurements and will improve future keratoconus diagnosis by providing more-accurate ET estimates.