Optoretinogram: optical measurement of human cone and rod photoreceptor responses to light.

Noninvasive, objective measurement of rod function is as significant as that of cone function, and for retinal diseases such as retinitis pigmentosa and age-related macular degeneration, rod function may be a more sensitive biomarker of disease progression and efficacy of treatment than cone function. Functional imaging of single human rod photoreceptors, however, has proven difficult because their small size and rapid functional response pose challenges for the resolution and speed of the imaging system. Here, we describe light-evoked, functional responses of human rods and cones, measured noninvasively using a synchronized adaptive optics optical coherence tomography (OCT) and scanning light ophthalmoscopy (SLO) system. The higher lateral resolution of the SLO images made it possible to confirm the identity of rods in the corresponding OCT volumes.

Coextensive synchronized SLO-OCT with adaptive optics for human retinal imaging.

We describe the details of a multimodal retinal imaging system which combines adaptive optics (AO) with an integrated scanning light ophthalmoscopy (SLO) and optical coherence tomography (OCT) imaging system. The OCT subsystem consisted of a swept-source, Fourier-domain mode-locked (FDML) laser, with a very high A-scan rate (1.6 MHz), whose beam was raster scanned on the retina by two scanners-one resonant scanner and one galvanometer. The high sweep rate of the FDML permitted the SLO and OCT to utilize the same scanners for in vivo retinal imaging and, unlike existing multimodal systems, concurrently acquired SLO frames and OCT volumes with approximate en face correspondence at a rate of 6 Hz. The AO provided diffraction-limited cellular resolution for both imaging channels.

Megahertz-rate optical coherence tomography angiography improves the contrast of the choriocapillaris and choroid in human retinal imaging.

Angiographic imaging of the human eye with optical coherence tomography (OCT) is becoming an increasingly important tool in the scientific investigation and clinical management of several blinding diseases, including age-related macular degeneration and diabetic retinopathy. We have observed that OCT angiography (OCTA) of the human choriocapillaris and choroid with a 1.64 MHz A-scan rate swept-source laser yields higher contrast images as compared to a slower rate system operating at 100 kHz. This result is unexpected because signal sensitivity is reduced when acquisition rates are increased, and the incident illumination power is kept constant. The contrast of angiography images generated by acquiring multiple sequential frames and calculating the variation caused by blood flow, however, appears to be improved significantly when lower-contrast images are taken more rapidly. To demonstrate that the acquisition rate plays a role in the quality improvement, we have imaged five healthy subjects with a narrow field of view (1.2 mm) OCTA imaging system using two separate swept-source lasers of different A-line rates and compared the results quantitatively using the radially-averaged power spectrum. The average improvement in the contrast is 23.0% (+/-7.6%). Although the underlying cause of this enhancement is not explicitly determined here, we speculate that the higher-speed system suppresses the noise contribution from eye motion in subjects and operates with an inter-scan time that better discriminates the flow velocities present in the choroid and choriocapillaris. Our result informs OCT system developers on the merits of ultrahigh-speed acquisition in functional imaging applications.

Functional retinal imaging using adaptive optics swept-source OCT at 1.6 MHz.

Objective optical assessment of photoreceptor function may permit earlier diagnosis of retinal disease than current methods such as perimetry, electrophysiology, and clinical imaging. In this work, we describe an adaptive optics (AO) optical coherence tomography (OCT) system designed to measure functional responses of single cones to visible stimuli. The OCT subsystem consisted of a raster-scanning Fourier-domain mode-locked laser that acquires A scans at 1.64 MHz with a center wavelength of 1063 nm and an AO system operating in closed-loop. Analysis of serial volumetric images revealed phase changes of cone photoreceptors consistent with outer segment elongation and proportional to stimulus intensity, as well as other morphological changes in the outer segment and retinal pigment epithelium.

Intraframe motion correction for raster-scanned adaptive optics images using strip-based cross-correlation lag biases.

In retinal raster imaging modalities, fixational eye movements manifest as image warp, where the relative positions of the beam and retina change during the acquisition of single frames. To remove warp artifacts, strip-based registration methods-in which fast-axis strips from target images are registered to a reference frame-have been applied in adaptive optics (AO) scanning light ophthalmoscopy (SLO) and optical coherence tomography (OCT). This approach has enabled object tracking and frame averaging, and methods have been described to automatically select reference frames with minimal motion. However, inconspicuous motion artifacts may persist in reference frames and propagate themselves throughout the processes of registration, tracking, and averaging. Here we test a previously proposed method for removing movement artifacts in reference frames, using biases in stripwise cross-correlation statistics. We applied the method to synthetic retinal images with simulated eye motion artifacts as well as real AO-SLO images of the cone mosaic and volumetric AO-OCT images, both affected by eye motion. In the case of synthetic images, the method was validated by direct comparison with motion-free versions of the images. In the case of real AO images, performance was validated by comparing the correlation of uncorrected images with that of corrected images, by quantifying the effect of motion artifacts on the image power spectra, and by qualitative examination of AO-OCT B-scans and en face projections. In all cases, the proposed method reduced motion artifacts and produced more faithful images of the retina.

The Properties of Outer Retinal Band Three Investigated With Adaptive-Optics Optical Coherence Tomography.

Purpose: Optical coherence tomography's (OCT) third outer retinal band has been attributed to the zone of interdigitation between RPE cells and cone outer segments. The purpose of this paper is to investigate the structure of this band with adaptive optics (AO)-OCT. Methods: Using AO-OCT, images were obtained from two subjects. Axial structure was characterized by measuring band 3 thickness and separation between bands 2 and 3 in segmented cones. Lateral structure was characterized by correlation of band 3 with band 2 and comparison of their power spectra. Band thickness and separation were also measured in a clinical OCT image of one subject. Results: Band 3 thickness ranged from 4.3 to 6.4 µm. Band 2 correlations ranged between 0.35 and 0.41 and power spectra of both bands confirmed peak frequencies that agree with histologic density measurements. In clinical images, band 3 thickness was between 14 and 19 µm. Measurements of AO-OCT of interband distance were lower than our corresponding clinical OCT measurements. Conclusions: Band 3 originates from a structure with axial extent similar to a single surface. Correlation with band 2 suggests an origin within the cone photoreceptor. These two observations indicate that band 3 corresponds predominantly to cone outer segment tips (COST). Conventional OCT may overestimate both the thickness of band 3 and outer segment length.

Fluorescence laminar optical tomography for brain imaging: system implementation and performance evaluation.

We present our effort in implementing a fluorescence laminar optical tomography scanner which is specifically designed for noninvasive three-dimensional imaging of fluorescence proteins in the brains of small rodents. A laser beam, after passing through a cylindrical lens, scans the brain tissue from the surface while the emission signal is captured by the epi-fluorescence optics and is recorded using an electron multiplication CCD sensor. Image reconstruction algorithms are developed based on Monte Carlo simulation to model light­tissue interaction and generate the sensitivity matrices. To solve the inverse problem, we used the iterative simultaneous algebraic reconstruction technique. The performance of the developed system was evaluated by imaging microfabricated silicon microchannels embedded inside a substrate with optical properties close to the brain as a tissue phantom and ultimately by scanning brain tissue in vivo. Details of the hardware design and reconstruction algorithms are discussed and several experimental results are presented. The developed system can specifically facilitate neuroscience experiments where fluorescence imaging and molecular genetic methods are used to study the dynamics of the brain circuitries.

Calibration of digital optical phase conjugation setups based on orthonormal rectangular polynomials.

Digital optical phase conjugation (DOPC) has proven to be a promising technique in deep tissue fluorescence imaging. Nonetheless, DOPC optical setups require precise alignment of all optical components to accurately read the wavefront of scattered light in a turbid medium and playback the conjugated beam toward the sample. Minor misalignments and possible imperfections in the arrangement or the structure of the optical components significantly reduce the performance of the method. In this paper, a calibration procedure based on orthogonal rectangular polynomials is introduced to compensate major imperfections including the optical aberration in the wavefront of the reference beam and the substrate curvature of the spatial light modulator without adding extra optical components to the original setup. The proposed algorithm also provides a systematic calibration procedure for mechanical fine tuning of DOPC systems. It is shown experimentally that the proposed calibration process improves the peak-to-background ratio when focusing light after passing through a highly scattering medium.

OptogenSIM: a 3D Monte Carlo simulation platform for light delivery design in optogenetics.

Optimizing light delivery for optogenetics is critical in order to accurately stimulate the neurons of interest while reducing nonspecific effects such as tissue heating or photodamage. Light distribution is typically predicted using the assumption of tissue homogeneity, which oversimplifies light transport in heterogeneous brain. Here, we present an open-source 3D simulation platform, OptogenSIM, which eliminates this assumption. This platform integrates a voxel-based 3D Monte Carlo model, generic optical property models of brain tissues, and a well-defined 3D mouse brain tissue atlas. The application of this platform in brain data models demonstrates that brain heterogeneity has moderate to significant impact depending on application conditions. Estimated light density contours can show the region of any specified power density in the 3D brain space and thus can help optimize the light delivery settings, such as the optical fiber position, fiber diameter, fiber numerical aperture, light wavelength and power. OptogenSIM is freely available and can be easily adapted to incorporate additional brain atlases.

Effect of blood vessels on light distribution in optogenetic stimulation of cortex.

In this Letter, the impact of blood vessels on light distribution during photostimulation of cortical tissue in small rodents is investigated. Brain optical properties were extracted using a double-integrating sphere setup, and optical coherence tomography was used to image cortical vessels and capillaries to generate a three-dimensional angiogram of the cortex. By combining these two datasets, a complete volumetric structure of the cortical tissue was developed and linked to a Monte Carlo code which simulates light propagation in this inhomogeneous structure and illustrates the effect of blood vessels on the penetration depth and pattern preservation in optogenetic stimulation.

Closed-Loop Optogenetic Brain Interface.

This paper presents a new approach for implementation of closed-loop brain-machine interface algorithms by combining optogenetic neural stimulation with electrocorticography and fluorescence microscopy. We used a new generation of microfabricated electrocorticography (micro-ECoG) devices in which electrode arrays are embedded within an optically transparent biocompatible substrate that provides optical access to the brain tissue during electrophysiology recording. An optical setup was designed capable of projecting arbitrary patterns of light for optogenetic stimulation and performing fluorescence microscopy through the implant. For realization of a closed-loop system using this platform, the feedback can be taken from electrophysiology data or fluorescence imaging. In the closed-loop systems discussed in this paper, the feedback signal was taken from the micro-ECoG. In these algorithms, the electrophysiology data are continuously transferred to a computer and compared with some predefined spatial-temporal patterns of neural activity. The computer which processes the data also readjusts the duration and distribution of optogenetic stimulating pulses to minimize the difference between the recorded activity and the predefined set points so that after a limited period of transient response the recorded activity follows the set points. Details of the system design and implementation of typical closed-loop paradigms are discussed in this paper.

Monitoring hemodynamic changes in stroke-affected muscles using near-infrared spectroscopy.

The oxygenation level of a tissue is an important marker of the health of the tissue and has a direct effect on performance. It has been shown that the blood flow to the paretic muscles of hemiparetic post-stroke patients is significantly reduced compared to non-paretic muscles. It is hypothesized that hemodynamic activity in paretic muscles is suppressed as compared to non-paretic muscles, and that oximetry can be used to measure this disparity in real-time. In order to test this hypothesis, a custom-made oximetry device was used to measure hemodynamic activity in the forearm extensor muscles in post-stroke patients' paretic and non-paretic sides and in a control population during three exercise levels calibrated to the subject's maximum effort. The change in oxygenation (ΔOxy) and blood volume (ΔBV) were calculated and displayed in real-time. Results show no apparent difference in either ΔOxy or ΔBV between control subjects' dominant and non-dominant muscles. However, the results show a significant difference in ΔOxy between paretic and non-paretic muscles, as well as a significant difference between normalized post-stroke and control data. Further work will be necessary to determine if the observed difference between the paretic and non-paretic muscles changes over the course of physical therapy and can be correlated with functional improvements.

Patterned optogenetic modulation of neurovascular and metabolic signals.

The hemodynamic and metabolic response of the cortex depends spatially and temporally on the activity of multiple cell types. Optogenetics enables specific cell types to be modulated with high temporal precision and is therefore an emerging method for studying neurovascular and neurometabolic coupling. Going beyond temporal investigations, we developed a microprojection system to apply spatial photostimulus patterns in vivo. We monitored vascular and metabolic fluorescence signals after photostimulation in Thy1-channelrhodopsin-2 mice. Cerebral arteries increased in diameter rapidly after photostimulation, while nearby veins showed a slower smaller response. The amplitude of the arterial response was depended on the area of cortex stimulated. The fluorescence signal emitted at 450/100 nm and excited with ultraviolet is indicative of reduced nicotinamide adenine dinucleotide, an endogenous fluorescent enzyme involved in glycolysis and the citric acid cycle. This fluorescence signal decreased quickly and transiently after optogenetic stimulation, suggesting that glucose metabolism is tightly locked to optogenetic stimulation. To verify optogenetic stimulation of the cortex, we used a transparent substrate microelectrode array to map cortical potentials resulting from optogenetic stimulation. Spatial optogenetic stimulation is a new tool for studying neurovascular and neurometabolic coupling.

Extraction of optical properties and prediction of light distribution in rat brain tissue.

Predicting the distribution of light inside any turbid media, such as biological tissue, requires detailed information about the optical properties of the medium, including the absorption and scattering coefficients and the anisotropy factor. Particularly, in biophotonic applications where photons directly interact with the tissue, this information translates to system design optimization, precision in light delivery, and minimization of unintended consequences, such as phototoxicity or photobleaching. In recent years, optogenetics has opened up a new area in deep brain stimulation with light and the method is widely adapted by researchers for the study of the brain circuitries and the dynamics of neurological disorders. A key factor for a successful optogenetic stimulation is delivering an adequate amount of light to the targeted brain objects. The adequate amount of light needed to stimulate each brain object is identified by the tissue optical properties as well as the type of opsin expressed in the tissue, wavelength of the light, and the physical dimensions of the targeted area. Therefore, to implement a precise light delivery system for optogenetics, detailed information about the optical properties of the brain tissue and a mathematical model that incorporates all determining factors is needed to find a good estimation of light distribution in the brain. In general, three measurements are required to obtain the optical properties of any tissue, namely diffuse transmitted light, diffuse reflected light, and transmitted ballistic beam. In this report, these parameters were measured in vitro using intact rat brain slices of 500 µm thickness via a two-integrating spheres optical setup. Then, an inverse adding doubling method was used to extract the optical properties of the tissue from the collected data. These experiments were repeated to cover the whole brain tissue with high spatial resolution for the three different cuts (transverse, sagittal, and coronal) and three different wavelengths (405, 532, and 635 nm) in the visible range of the spectrum. A three-dimensional atlas of the rat brain optical properties was constructed based on the experimental measurements. This database was linked to a Monte Carlo toolbox to simulate light distribution in the tissue for different light source configurations.

Nonlinear optical signal processing on multiwavelength sensitive materials.

Exploiting salient features in the photodynamics of specific types of light sensitive materials, a new approach is presented for realization of parallel nonlinear operations with optics. We briefly review the quantum structure and mathematical models offered for the photodynamics of two multiwavelength sensitive materials, doped crystals of lithium niobate and thick layers of bacteriorhodopsin. Next, a special mode of these dynamics in each material is investigated and a graphical design procedure is offered to produce highly nonlinear optical responses that can be dynamically reshaped via applying minimum changes in the optical setup.