Assessment of Different Sampling Methods for Measuring and Representing Macular Cone Density Using Flood-Illuminated Adaptive Optics

The effect of sampling window size on topographical maps of foveal cone density

Purpose

To characterize the effect of sampling window size on maps of foveal cone density derived from adaptive optics scanning light ophthalmoscope (AOSLO) images of the cone mosaic.

Methods

Forty-four AOSLO-derived montages of the foveal cone mosaic (300 x 300µm) were used for this study (from 44 individuals with normal vision). Cone photoreceptor coordinates were semi-automatically identified by one experienced grader. From these coordinates, cone density matrices across each foveal montage were derived using 10 different sampling window sizes containing 5, 10, 15, 20, 40, 60, 80, 100, 150, or 200 cones. For all 440 density matrices, we extracted the location and value of peak cone density (PCD), the cone density centroid (CDC) location, and cone density at the CDC.

Results

Across all window sizes, PCD values were larger than those extracted at the CDC location, though the difference between these density values decreased as the sampling window size increased (p<0.0001). Overall, both PCD (r=-0.8099, p=0.0045) and density at the CDC (r=-0.7596, p=0.0108) decreased with increasing sampling window size. This reduction was more pronounced for PCD, with a 27.8% lower PCD value on average when using the 200-cone versus the 5-cone window (compared to only a 3.5% reduction for density at the CDC between these same window sizes). While the PCD and CDC locations did not occur at the same location within a given montage, there was no significant relationship between this PCD-CDC offset and sampling window size (p=0.8919). The CDC location was less variable across sampling windows, with an average per-participant 95% confidence ellipse area across the 10 window sizes of 47.56µm² (compared to 844.10µm² for the PCD location, p<0.0001).

Conclusion

CDC metrics appear more stable across varying sampling window sizes than PCD metrics. Understanding how density values change according to the method used to sample the cone mosaic may facilitate comparing cone density data across different studies.

Reduced Cone Density Is Associated with Multiple Sclerosis

Purpose

Multiple sclerosis (MS) is an inflammatory neurodegenerative disease of the central nervous system. Recent evidence suggests that degeneration of the inner layers of the retina occurs in MS. This study aimed to examine whether there are outer retinal changes in patients living with MS.

Design

This was a single center, cross-sectional study.

Participants

Sixteen patients with MS and 25 controls (volunteers without diagnosed MS) were recruited for the study.

Methods

We acquired volumetric spectral domain-OCT scans of the macula and a circular scan around the optic nerve head (ONH). We also captured adaptive optics (AO) images at 0° (centered on the foveola), 2°, 4°, and 6° temporal to the fovea.

Main Outcome Measures

We calculated the thickness of the different retinal layers in the macula and around the ONH using the inbuilt software of the OCT. We evaluated changes in cone photoreceptors by calculating cone density and spacing by the inbuilt AO automatic segmentation algorithm with manual correction. We compared patients with and without optic neuritis and controls.

Results

We found significant thinning of the inner retina and a thickening of the outer retina in the eye with a history of optic neuritis (eyes of patients with MS with a history of optic neuritis; mean difference [MD]: -11.13 ± 3.61 μm, P = 0.002 and MD: 2.86 ± 0.89 μm, P = 0.001; respectively). We did not observe changes in retinal layers without optic neuritis in eyes of patients with MS without a history of optic neuritis. However, regional differences were detected in the peripapillary retinal nerve fiber layer. Analyzing AO images revealed a significantly lower cone outer-segment density at all eccentricities in all patients compared with control eyes (P < 0.05), independent of optic neuritis history.

Conclusions

Our results showed that all MS cases were associated with decreased cone densities. Future longitudinal studies will help to elucidate whether this is a specific and sensitive method to detect and monitor the development and progression of MS.

Financial Disclosures

Proprietary or commercial disclosure may be found after the references.

Change in Cone Structure Over 24 Months in USH2A-Related Retinal Degeneration

PURPOSE

To describe cone structure changes using adaptive optics scanning laser ophthalmoscopy (AOSLO) in the Rate of Progression of USH2A-related Retinal Degeneration (RUSH2A) study.

DESIGN

Multicenter, longitudinal natural history study.

METHODS

AOSLO images were acquired at 4 centers, twice at baseline and annually for 24 months in this natural history study. For each eye, at least 10 regions of interest (ROIs) with ≥50 contiguous cones were analyzed by masked, independent graders. Cone spacing Z-scores, standard deviations from the normal mean at the measured location, were compared between graders and tests at baseline. The association of cone spacing with clinical characteristics was assessed using linear mixed effects regression models weighted by image quality score. Annual rates of change were calculated based on differences between visits.

RESULTS

Fourteen eyes of 14 participants were imaged, with 192 ROIs selected at baseline. There was variability among graders, which was greater in images with lower image quality score (P < .001). Cone spacing was significantly correlated with eccentricity, quality score, and disease duration (P < .02). On average, the cone spacing Z-score increased 0.14 annually (about 9%, P < .001). We observed no significant differences in rate of change between disease type (Usher syndrome or retinitis pigmentosa), imaging site, or grader.

CONCLUSIONS

Using current methods, the analysis of quantitative measures of cone structure showed some challenges, yet showed promise that AOSLO images can be used to characterize progressive change over 24 months. Additional multicenter studies using AOSLO are needed to advance cone mosaic metrics as sensitive outcome measures for clinical trials. NOTE: Publication of this article is sponsored by the American Ophthalmological Society.

Pearls and Pitfalls of Adaptive Optics Ophthalmoscopy in Inherited Retinal Diseases

Adaptive optics (AO) retinal imaging enables individual photoreceptors to be visualized in the clinical setting. AO imaging can be a powerful clinical tool for detecting photoreceptor degeneration at a cellular level that might be overlooked through conventional structural assessments, such as spectral-domain optical coherence tomography (SD-OCT). Therefore, AO imaging has gained significant interest in the study of photoreceptor degeneration, one of the most common causes of inherited blindness. Growing evidence supports that AO imaging may be useful for diagnosing early-stage retinal dystrophy before it becomes apparent on fundus examination or conventional retinal imaging. In addition, serial AO imaging may detect structural disease progression in early-stage disease over a shorter period compared to SD-OCT. Although AO imaging is gaining popularity as a structural endpoint in clinical trials, the results should be interpreted with caution due to several pitfalls, including the lack of standardized imaging and image analysis protocols, frequent ocular comorbidities that affect image quality, and significant interindividual variation of normal values. Herein, we summarize the current state-of-the-art AO imaging and review its potential applications, limitations, and pitfalls in patients with inherited retinal diseases.

Comprehensive automatic processing and analysis of adaptive optics flood illumination retinal images on healthy subjects

This work presents a novel fully automated method for retinal analysis in images acquired with a flood illuminated adaptive optics retinal camera (AO-FIO). The proposed processing pipeline consists of several steps: First, we register single AO-FIO images in a montage image capturing a larger retinal area. The registration is performed by combination of phase correlation and the scale-invariant feature transform method. A set of 200 AO-FIO images from 10 healthy subjects (10 images from left eye and 10 images from right eye) is processed into 20 montage images and mutually aligned according to the automatically detected fovea center. As a second step, the photoreceptors in the montage images are detected using a method based on regional maxima localization, where the detector parameters were determined with Bayesian optimization according to manually labeled photoreceptors by three evaluators. The detection assessment, based on Dice coefficient, ranges from 0.72 to 0.8. In the next step, the corresponding density maps are generated for each of the montage images. As a final step, representative averaged photoreceptor density maps are created for the left and right eye and thus enabling comprehensive analysis across the montage images and a straightforward comparison with available histological data and other published studies. Our proposed method and software thus enable us to generate AO-based photoreceptor density maps for all measured locations fully automatically, and thus it is suitable for large studies, as those are in pressing need for automated approaches. In addition, the application MATADOR (MATlab ADaptive Optics Retinal Image Analysis) that implements the described pipeline and the dataset with photoreceptor labels are made publicly available.

Intergrader agreement of foveal cone topography measured using adaptive optics scanning light ophthalmoscopy

The foveal cone mosaic can be directly visualized using adaptive optics scanning light ophthalmoscopy (AOSLO). Previous studies in individuals with normal vision report wide variability in the topography of the foveal cone mosaic, especially the value of peak cone density (PCD). While these studies often involve a human grader, there have been no studies examining intergrader reproducibility of foveal cone mosaic metrics. Here we re-analyzed published AOSLO foveal cone images from 44 individuals to assess the relationship between the cone density centroid (CDC) location and the location of PCD. Across 5 graders with variable experience, we found a measurement error of 11.7% in PCD estimates and higher intergrader reproducibility of CDC location compared to PCD location (p < 0.0001). These estimates of measurement error can be used in future studies of the foveal cone mosaic, and our results support use of the CDC location as a more reproducible anchor for cross-modality analyses.

Comparison of Cone Mosaic Metrics From Images Acquired With the SPECTRALIS High Magnification Module and Adaptive Optics Scanning Light Ophthalmoscopy

Purpose

To compare cone mosaic metrics derived from adaptive optics scanning light ophthalmoscopy (AOSLO) images with those derived from Heidelberg Engineering SPECTRALIS High Magnification Module (HMM) images.

Methods

Participants with contiguous cone mosaics had HMM imaging performed at locations superior and temporal to the fovea. These images were registered and averaged offline and then aligned to split-detection AOSLO images; 200 × 200-µm regions of interest were extracted from both modalities. Cones were semi-automatically identified by two graders to provide estimates of cone density and spacing.

Results

Thirty participants with contiguous cone mosaics were imaged (10 males, 20 females; age range, 11-67 years). Image quality varied, and 80% of our participants had analyzable HMM images. The intergrader intraclass correlation coefficients for cone metrics were good for both modalities (0.688-0.757 for HMM; 0.805-0.836 for AOSLO). Cone density estimates from HMM images were lower by 2661 cones/mm2 (24.1%) on average compared to AOSLO-derived estimates. Accordingly, HMM estimates of cone spacing were increased on average compared to AOSLO.

Conclusions

The cone mosaic can be visualized in vivo using the SPECTRALIS HMM, although image quality is variable and imaging is not successful in every individual. Metrics extracted from HMM images can differ from those from AOSLO, although excellent agreement is possible in individuals with excellent optical quality and precise co-registration between modalities.

Translational Relevance

Emerging non-adaptive optics-based photoreceptor imaging is more clinically accessible than adaptive optics techniques and has potential to expand high-resolution imaging in a clinical environment.

Application of Adaptive Optics in Ophthalmology

The eye, the photoreceptive organ used to perceive the external environment, is of great importance to humans. It has been proven that some diseases in humans are accompanied by fundus changes; therefore, the health status of people may be interpreted from retinal images. However, the human eye is not a perfect refractive system for the existence of ocular aberrations. These aberrations not only affect the ability of human visual discrimination and recognition, but restrict the observation of the fine structures of human eye and reduce the possibility of exploring the mechanisms of eye disease. Adaptive optics (AO) is a technique that corrects optical wavefront aberrations. Once integrated into ophthalmoscopes, AO enables retinal imaging at the cellular level. This paper illustrates the principle of AO in correcting wavefront aberrations in human eyes, and then reviews the applications and advances of AO in ophthalmology, including the adaptive optics fundus camera (AO-FC), the adaptive optics scanning laser ophthalmoscope (AO-SLO), the adaptive optics optical coherence tomography (AO-OCT), and their combined multimodal imaging technologies. The future development trend of AO in ophthalmology is also prospected.

Impact of Reference Center Choice on Adaptive Optics Imaging Cone Mosaic Analysis

Purpose

Foveal center marking is a key step in retinal image analysis. We investigated the discordance between the adaptive optics (AO) montage center (AMC) and the foveal pit center (FPC) and its implications for cone mosaic analysis using a commercial flood-illumination AO camera.

Methods

Thirty eyes of 30 individuals (including 15 healthy and 15 patients with rod–cone dystrophy) were included. Spectral-domain optical coherence tomography was used to determine the FPC, and flood-illumination AO imaging was performed with overlapping image frames to create an AO montage. The AMC was determined by averaging the (0,0) coordinates in the four paracentral overlapping AO image frames. Cone mosaic measurements at various retinal eccentricities were compared between corresponding retinal loci relative to the AMC or FPC.

Results

AMCs were located temporally to the FPCs in 14 of 15 eyes in both groups. The average AMC–FPC discordance was 0.85° among healthy controls and 0.33° among patients with rod-cone dystrophy (P < 0.05). The distance of the AMC from the FPC was a significant determinant of the cone density (β estimate = 218 cells/deg²/deg; 95% confidence interval [CI], 107–330; P < 0.001) and inter-cone distance (β estimate = 0.28 arcmin/deg; 95% CI, 0.15–0.40; P < 0.001), after adjustment for age, sex, axial length, spherical equivalent, eccentricity, and disease status.

Conclusions

There is a marked mismatch between the AMC and FPC in healthy eyes that may be modified by disease process such as rod–cone dystrophy. We recommend users of AO imaging systems carefully align the AO montage with a foveal anatomical landmark, such as the FPC, to ensure precise and reproducible localization of the eccentricities and regions of interest for cone mosaic analysis.

Short-Term Parafoveal Cone Loss Despite Preserved Ellipsoid Zone in Rod Cone Dystrophy

Purpose

Rod–cone dystrophy (RCD) is characterized by centripetal loss of rod followed by cone photoreceptors. In this prospective, observational cohort, we used flood-illumination adaptive optics (AO) imaging to investigate parafoveal cone loss in regions with preserved ellipsoid zone (EZ) in patients with RCD.

Methods

Eight patients with RCD and 10 age-matched healthy controls underwent spectral-domain optical coherence tomography and AO imaging. The RCD cohort underwent a follow-up examination after 6 months. Cone density (CD) and intercone distance (ICD) measurements were performed at 2° temporal from the fovea. Baseline CD and ICD values were compared between the control and patient groups, and longitudinal changes were calculated in the patient group. Residual EZ span in patients was measured in horizontal foveal B-scans.

Results

Between the control and patient groups, there was no significant difference in the baseline CD (2094 vs. 1750 cells/deg², respectively; P = 0.09) and ICD (1.46 vs. 1.62 arcmin, respectively; P = 0.08). Mean CD declined by 198 cells/deg² (−11.3%; P < 0.01), and mean ICD increased by 0.09 arcmin (+5.6%; P = 0.01) at the 6-month follow-up in the patient group. Mean baseline and follow-up residual EZ spans in the six patients with EZ defect were 3189 µm and 3065 µm, respectively (−3.9%; P = 0.08).

Conclusions

AO imaging detected significant parafoveal cone loss over 6-month follow-up even in regions with preserved EZ. Further studies to refine AO imaging protocol and validate cone metrics as a structural endpoint in early RCD are warranted.

Translational Relevance

CD and ICD may change prior to EZ span shortening in RCD.

Cone photoreceptor density in the Copenhagen Child Cohort at age 16-17 years

PURPOSE

To examine cone density in relation to gestational and morphological parameters in the Copenhagen Child Cohort (CCC2000).

METHODS

The macula was imaged using adaptive optics in 1,296 adolescents aged 16-17 years. Axial length and distance visual acuity were determined. Absolute and angular cone photoreceptor density were analysed for an 80 × 80-pixel area, 2 degrees temporal to the fovea. Association with axial length was analysed with linear regression. Correlation with visual acuity was described with a Pearson correlation coefficient. Associations of cone density with gestational parameters, maternal smoking, sex and age were analysed using multiple regression adjusted for axial length.

RESULTS

Mean absolute cone density was 30,007 cones/mm² (SD ± 3,802) and mean angular cone density was 2,383 cones/deg² (SD ± 231). Peri- and postnatal parameters, sex and age had no statistically significant effect on cone density (p > 0.05). Absolute cone density decreased with longer axial length (-2,855 cones/mm² per mm or -9.7% per mm, p < 0.0001). For angular density, which included a correction for the geometrical enlargement of the eye with axial length, a decrease with axial length was detectable, but it was small (-20 cones/deg² per mm or -0.84% per mm, p = 0.009).

CONCLUSIONS

The decrease in cone density per unit solid angle with increasing axial length was small, less than 1 percent per mm, indicating that expansion of the posterior pole during the development of refraction takes place without a clinically significant loss of cones. Perinatal parameters, within the spectrum presented by the study population, had no detectable effect on cone density.

Evaluation of focus and deep learning methods for automated image grading and factors influencing image quality in adaptive optics ophthalmoscopy

Adaptive optics flood illumination ophthalmoscopy (AO-FIO) is an established imaging tool in the investigation of retinal diseases. However, the clinical interpretation of AO-FIO images can be challenging due to varied image quality. Therefore, image quality assessment is essential before interpretation. An image assessment tool will also assist further work on improving the image quality, either during acquisition or post processing. In this paper, we describe, validate and compare two automated image quality assessment methods; the energy of Laplacian focus operator (LAPE; not commonly used but easily implemented) and convolutional neural network (CNN; effective but more complex approach). We also evaluate the effects of subject age, axial length, refractive error, fixation stability, disease status and retinal location on AO-FIO image quality. Based on analysis of 10,250 images of 50 × 50 μm size, at 41 retinal locations, from 50 subjects we demonstrate that CNN slightly outperforms LAPE in image quality assessment. CNN achieves accuracy of 89%, whereas LAPE metric achieves 73% and 80% (for a linear regression and random forest multiclass classifier methods, respectively) compared to ground truth. Furthermore, the retinal location, age and disease are factors that can influence the likelihood of poor image quality.

Retinal Differential Light Sensitivity Variation Across the Macula in Healthy Subjects: Importance of Cone Separation and Loci Eccentricity

Purpose

Microperimetry measures differential light sensitivity (DLS) at specific retinal locations. The aim of this study is to examine the variation in DLS across the macula and the contribution to this variation of cone distribution metrics and retinal eccentricity.

Methods

Forty healthy eyes of 40 subjects were examined by microperimetry (MAIA) and adaptive optics imaging (rtx1). Retinal DLS was measured using the grid patterns: foveal (2°–3°), macular (3°–7°), and meridional (2°–8° on horizontal and vertical meridians). Cone density (CD), distribution regularity, and intercone distance (ICD) were calculated at the respective test loci coordinates. Linear mixed-effects regression was used to examine the association between cone distribution metrics and loci eccentricity, and retinal DLS.

Results

An eccentricity-dependent reduction in DLS was observed on all MAIA grids, which was greatest at the foveal-parafoveal junction (2°–3°) (−0.58 dB per degree, 95% confidence interval [CI]; −0.91 to −0.24 dB, P < 0.01). Retinal DLS across the meridional grid changed significantly with each 1000 cells/deg² change in CD (0.85 dB, 95% CI; 0.10 to 1.61 dB, P = 0.03), but not with each arcmin change in ICD (1.36 dB, 95% CI; −2.93 to 0.20 dB, P = 0.09).

Conclusions

We demonstrate significant variation in DLS across the macula. Topographical change in cone separation is an important determinant of the variation in DLS at the foveal-parafoveal junction. We caution the extrapolation of changes in DLS measurements to cone distribution because the relationship between these variables is complex.

Translational Relevance

Cone density is an independent determinant of DLS in the foveal-parafoveal junction in healthy eyes.

Multimodal imaging analysis of macular dystrophy in patient with maternally inherited diabetes and deafness (MIDD) with m.3243A>G mutation

Purpose: Maternally inherited diabetes and deafness (MIDD) is caused by a heteroplasmic m.3243A>G mutation in the mitochondrial DNA. The main ocular feature in MIDD is macular dystrophy. The purpose of this study was to identify the phenotypical characteristics of a patient with MIDD by multimodal high-resolution imaging analyses.Methods: A detailed history and ophthalmic examination were performed on a 39-year-old patient with MIDD. Multi-modal imaging included fundus photography, fundus autofluorescence imaging, fluorescein angiography, spectral-domain optical coherence tomography, OCT-angiography, and adaptive optics imaging. The PCR-invader and whole exome sequencing (WES) methods were performed on the DNA of the patient.Results: A 39-year-old woman with sensorineural hearing loss, diabetes mellitus presented with atrophic perifoveal changes and MIDD was suspected. The PCR-invader and WES methods showed that the patient had a m.3243A>G mutation in the mitochondrial DNA with 29% and 16.7% of the heteroplasmy in the peripheral blood, respectively. Morphological analyses revealed that the areas of photoreceptor degeneration and chorioretinal atrophy were present mainly in the perifoveal region. Multifocal ERGs showed that the perifoveal responses were reduced. Goldmann visual field was significant for a cecocentral scotoma in the right eye and an enlarged blind spot in the left eye. The central isopter was constricted bilaterally. The results of high-resolution retinal imaging by AO revealed that the densities of the cone photoreceptor were significantly reduced in the fovea where no obvious atrophy of the RPE and choroid was observed.Conclusions: Our findings indicate that WES analysis can be used to detect the m.3243A>G mutation in the mtDNA. The results of multimodal imaging analyses indicated that the primary dysfunction of the photoreceptors in the fovea might precede the dysfunction of the RPE in patient with MIDD.

Adaptive optics ophthalmoscopy, multifocal ERG and OCTA in unique case of suspected torpedo maculopathy presenting with vitelliform lesion

Purpose

To report the case of a ten-year old girl with torpedo maculopathy with a complete vitelliform lesion and describe associated optical coherence tomography (OCT), OCT angiography (OCTA), multifocal electroretinogram (ERG) and adaptive optics ophthalmoscopy (AOO) imaging of the lesion.

Observations

An asymptomatic ten-year old girl with visual acuity of 20/15 OU was referred for evaluation of possible Best's disease of her left eye. The unilaterality, location, and shape of the lesion was consistent with torpedo maculopathy. OCT and autofluorescence (AF) revealed that the entire lesion was composed of subretinal hyperreflective material that was hyperautofluorescent, consistent with vitelliform material. Within the boundary of the lesion, OCTA showed reduced choriocapillaris density while adjacent to the lesion, the choriocapillaris density was slightly increased. Microperimetry demonstrated normal sensitivity in both eyes, electrooculograms (EOG) were normal and multifocal ERG showed symmetrical mildly supernormal amplitudes. Additionally, AOO demonstrated that nasal to the lesion there were clusters of hyper-reflective areas, and immediately adjacent to the lesion cones were poorly resolved. However, there was a return to more normal photoreceptor architecture outside of the lesion.

Conclusions and Importance

Torpedo maculopathy lesions typically present with outer retinal attenuation and retinal pigmented epithelium (RPE) atrophy. Vitelliform material was recently observed for the first time in association with Torpedo maculopathy in a case report that described small vitelliform material at the periphery of the lesion. We report the second case of torpedo maculopathy associated with a vitelliform lesion and the first description of a torpedo lesion composed fully of presumed vitelliform material. We also describe findings of OCTA, multifocal ERG and AOO imaging in torpedo maculopathy with vitelliform lesion.

Three-dimensional assessment of para- and perifoveal photoreceptor densities and the impact of meridians and age in healthy eyes with adaptive-optics optical coherence tomography (AO-OCT)

An adaptive optics optical coherence tomography (AO-OCT) system is used to assess sixty healthy eyes of thirty subjects (age 22 to 75) to evaluate how the outer retinal layers, foveal eccentricity and age effect the mean cone density. The cone mosaics of different retinal planes (the cone inner segment outer segment junction (IS/OS), the cone outer segment combined with the IS/OS (ISOS+), the cone outer segment tips (COST), and the full en-face plane (FEF)) at four main meridians (superior, nasal, inferior, temporal) and para- and perifoveal eccentricities (ecc 2.5° and 6.5°) were analyzed quantitatively. The mean overall cone density was 19,892/mm² at ecc 2.5° and 13,323/mm² at ecc 6.5°. A significant impact on cone density was found for eccentricity (up to 6,700/mm² between ecc 2.5° and 6.5°), meridian (up to 3,700/mm² between nasal and superior meridian) and layer (up to 1,400/mm² between FEF and IS/OS). Age showed only a weak negative effect. These factors as well as inter-individual variability have to be taken into account when comparing cone density measurements between healthy and pathologically changed eyes, as their combined effect on density can easily exceed several thousand cones per mm² even in parafoveal regions.

High-resolution photoreceptor imaging analysis of patients with autosomal dominant retinitis pigmentosa (adRP) caused by HK1 mutation

PURPOSE

The hexokinase 1 (HK1) gene encodes one of the four human hexokinases that play essential roles in glucose metabolism. Recently, several cases of E847K mutation in the HK1 gene were reported to cause inherited retinal dystrophy. The purpose of this study was to identify the phenotypical characteristics of patients with a recurrent E847K mutation in the HK1 gene.

METHODS

Three generations of one family with autosomal dominant retinitis pigmentosa were examined. Whole exome sequencing was performed on the DNA. Fundus imaging by an adaptive optics fundus camera was used to obtain high-resolution photoreceptor images.

RESULTS

Fundus examination of the proband showed degeneration of the mid-peripheral retina, and SD-OCT images showed an absence of the ellipsoid zone (EZ) and interdigitation zone (IZ) in the parafovea and more peripherally. SD-OCT images of the mother of the proband showed an absence of the EZ and IZ, and fundus autofluorescence images showed hypo-autofluorescence surrounding the macular region. One daughter of the proband had only mild night blindness, however, the density of the cone photoreceptors was reduced in the parafoveal region. Whole exome sequencing identified a heterozygous variant, E847K, in the HK1 gene. This variant was found to co-segregate with the disease in three family members.

CONCLUSIONS

Although the systemic phenotypes were found to be associated with the HK1 mutations, only the E847K mutation can cause a non-syndromic photoreceptor degeneration. Our study strengthened the hypothesis that the amino acid E847 might play a critical role in the maintenance of the morphology and function of the photoreceptors.

Cellular imaging of inherited retinal diseases using adaptive optics

Adaptive optics (AO) is an insightful tool that has been increasingly applied to existing imaging systems for viewing the retina at a cellular level. By correcting for individual optical aberrations, AO offers an improvement in transverse resolution from 10-15 μm to ~2 μm, enabling assessment of individual retinal cell types. One of the settings in which its utility has been recognised is that of the inherited retinal diseases (IRDs), the genetic and clinical heterogeneity of which warrants better cellular characterisation. In this review, we provide a summary of the basic principles of AO, its integration into multiple retinal imaging modalities and its clinical applications, focusing primarily on IRDs. Furthermore, we present a comprehensive summary of AO-based cellular findings in IRDs according to their associated disease-causing genes.

Repeatability of Adaptive Optics Automated Cone Measurements in Subjects With Retinitis Pigmentosa and Novel Metrics for Assessment of Image Quality

Purpose

We determine the intersession repeatability of cone measurements via flood-illuminated adaptive optics (AO) imaging in patients with retinitis pigmentosa (RP), to better differentiate variation due to imaging inaccuracies versus pathology-driven change.

Methods

A total of 25 4° × 4° AO images were acquired three times on the same day in 10 subjects with RP, registered in i2K Retina, and cones were identified using a custom-built MATLAB algorithm. Nine equally spaced regions of interest were selected for each imaging set. A subset of subjectively “poor” and “good” quality images was selected by three independent graders, analyzed using cone density, cone location similarity (CLS) and cone spacing, and compared to age-matched normals.

Results

The coefficient of variation (CoV), repeatability, and percent repeatability of automated cone density were slightly higher in patients with RP compared to age-matched normals, but showed no statistically significant difference. The standard deviation of CLS and cone spacing of nearest-neighbor distance demonstrated a statistically significant difference between good- and poor-quality images.

Conclusions

Repeatability of automated cone density measurements in patients with RP is comparable to normals. Misidentification of cones due to image quality variability is a major limitation of automated cone counting algorithms in patients with RP. Our study suggests that CLS and cone spacing metrics could be used to help define image quality and, thus, increase confidence in automated cone counts in patients with RP.

Translational Relevance

The novel AO image quality assessment metrics described in our study could help to improve patient image interpretation, prognosis, and longitudinal care.

High-Resolution Retinal Imaging Reveals Preserved Cone Photoreceptor Density and Choroidal Thickness in Female Carriers of Choroideremia

BACKGROUND AND OBJECTIVE

To characterize the photoreceptors and choroidal morphology of heterozygous female carriers of choroideremia who typically do not have any visual defects but can have severe funduscopic changes.

PATIENTS AND METHODS

This was a clinical case series study. Detailed ophthalmic examinations were performed on six female carriers from four families with choroideremia. The subfoveal choroidal thickness (SFCT) was determined by spectral-domain optical coherence tomography (SD-OCT) and the cone photoreceptor density by adaptive optics (AO) retinal imaging. SFCT and cone densities of the carriers were compared to that of normal eyes of healthy subjects.

RESULTS

The mean age of the carriers was 42.5 years. Fundus photographs showed diffuse, patchy depigmentation; however, the SFCT was within the normal limits. AO retinal imaging revealed preserved cone densities at temporal eccentricities from 2 to 8 angular degrees.

CONCLUSIONS

The findings indicate that despite the presence of distinctive depigmentation of the retinal pigment epithelium in female carriers of choroideremia, their cone photoreceptor densities and SFCT are well-preserved. These observations may account for the good visual acuity and lack of an awareness of visual disturbances. [Ophthalmic Surg Lasers Imaging Retina. 2019;50:76-85.].

Distribution of mid-peripheral cones in emmetropic and myopic subjects using adaptive optics flood illumination camera

PURPOSE

We measured in vivo cone photoreceptors up to 24° of eccentricity along the horizontal meridian of healthy human retina. We also investigated the impact on cone densities of axial eye length elongation occurring with myopia.

METHODS

Using a flood illumination device coupled with an adaptive optics system, rtx1™, ( www.imagine-eyes.com), 55 right healthy retinas were imaged along the horizontal (i.e. nasal and temporal) meridian over a 48° field (i.e. from 3° to 24° each 3°). Then, cones were manually detected within 80 × 80 pixel regions of interest. Cone density and packing geometry (i.e. number of neighbours) were calculated (AOdetect software™). Subjects were divided into three groups: a group of 36 emmetropic (i.e. refractive error from -0.25D to +0.50D) subjects; a group of 10 low myopic subjects (i.e. refractive error from -0.50D to -2.50D); and a group of nine high myopic subjects (i.e. >-2.50D).

RESULTS

Cone density decreased with eccentricity in both semi-meridians. The decrease in cone photoreceptors occurred mainly in the first 9°. The difference of cone density between the nasal and temporal semi-meridian increased with eccentricity from 0.6% at 3° to 26% at 24°. Average cone density of emmetropes (850 cones deg^-2 or 11 087 cones mm^-2 ), low myopes (830 cones deg^-2 or 9731 cones mm^-2 ), and high myopes (912 cones deg^-2 or 9744 cones mm^-2 ), suggested that the retinas of the high myopic subjects were more stretched than the low myopic subjects retinas and even more stretched than that of the emmetropes. The axial eyeball elongation (square of the ratio of the axial eye length of 9%) seems to explain the cone density (11%) difference between emmetropes and low myopes. However, while the eyeball elongation between low and high myopes is still important (i.e. 11%), cone density difference between both populations was negligible (i.e. 3%). The ratio of cone density varied from -17% to 22% as a function of eccentricity involving that the retinal stretching is not uniform along the horizontal meridian.

CONCLUSION

The difference of cone density (i.e. cone mm^-2 ) between groups supports the hypothesis that the retina is stretched with the eyeball elongation. However, this elongation does not seem to be uniform along the horizontal meridian favouring the hypothesis of a local elongation of the retina.

CEP250 mutations associated with mild cone-rod dystrophy and sensorineural hearing loss in a Japanese family

BACKGROUND

CEP250 encodes the C-Nap1 protein which belongs to the CEP family of proteins. C-Nap1 has been reported to be expressed in the photoreceptor cilia and is known to interact with other ciliary proteins. Mutations of CEP250 cause atypical Usher syndrome which is characterized by early-onset sensorineural hearing loss (SNHL) and a relatively mild retinitis pigmentosa. This study tested the hypothesis that the mild cone-rod dystrophy (CRD) and SNHL in a non-consanguineous Japanese family was caused by CEP250 mutations.

METHODS

Detailed ophthalmic and auditory examinations were performed on the proband and her family members. Whole exome sequencing (WES) was used on the DNA obtained from the proband.

RESULTS

Electrophysiological analysis revealed a mild CRD in two family members. Adaptive optics (AO) imaging showed reduced cone density around the fovea. Auditory examinations showed a slight SNHL in both patients. WES of the proband identified compound heterozygous variants c.361C>T, p.R121*, and c.562C>T, p.R188* in CEP250. The variants were found to co-segregate with the disease in five members of the family.

CONCLUSIONS

The variants of CEP250 are both null variants and according to American College of Medical Genetics and Genomics (ACMG) standards and guideline, these variants are classified into the very strong category (PVS1). The criteria for both alleles will be pathogenic. Our data indicate that mutations of CEP250 can cause mild CRD and SNHL in Japanese patients. Because the ophthalmological phenotypes were very mild, high-resolution retinal imaging analysis, such as AO, will be helpful in diagnosing CEP250-associated disease.

Visual resolution and cone spacing in the nasal and inferior retina

PURPOSE

To determine the retinal eccentricity at which cones are no longer an observable substitute for ganglion cells on nasal and inferior parafoveal visual acuity.

METHOD

Visual acuities were measured on 12 healthy volunteers, under dynamic adaptive optic aberrations correction (crx1™) in white light, from 0° to 6°, every two degrees, along the nasal and inferior retinal meridians. Cone spacing was measured on images of the retina acquired using an adaptive optic flood illumination retina camera (rtx1™) at the same eccentricity, except at 0°.

RESULTS

Cone spacing increased by around 0.13 min of arc per degree of eccentricity and a difference of 7% between both meridians was observed (higher cone spacing in the inferior retinal meridian). Visual resolution was higher in the nasal retinal meridian (difference of around 28% or 0.15 logMAR at 6°). Cone spacing can predict minimum angle of resolution (MAR) at 2° in both semi retinal meridians. In the inferior retinal meridian, MAR measurements are fairly well predicted by Watson's 50% mathematical model based on the midget retinal ganglion cell density. Along the nasal retinal meridian, the measured MAR lies between Watson's 50% and 100% models.

CONCLUSIONS

At 2° of eccentricity, cone density accurately predicts visual resolution in both the nasal and inferior retina, supporting the idea that only 50% of the foveal midget retinal ganglion cells determine VA. The 50% model can also predict VA in the inferior retinal meridian at 4° and 6° of eccentricity. However, the 50% model underestimated visual acuity in the nasal retinal meridian at 4° and 6° of eccentricity consistent with the partially overlapping ON and OFF midget retinal ganglion cell receptive fields.

Distribution of cone density, spacing and arrangement in adult healthy retinas with adaptive optics flood illumination

The aim of this article is to analyse cone density, spacing and arrangement using an adaptive optics flood illumination retina camera (rtx1^™) on a healthy population. Cone density, cone spacing and packing arrangements were measured on the right retinas of 109 subjects at 2°, 3°, 4°, 5° and 6° of eccentricity along 4 meridians. The effects of eccentricity, meridian, axial length, spherical equivalent, gender and age were evaluated. Cone density decreased on average from 28 884 ± 3 692 cones/mm², at 2° of eccentricity, to 15 843 ± 1 598 cones/mm² at 6°. A strong inter-individual variation, especially at 2°, was observed. No important difference of cone density was observed between the nasal and temporal meridians or between the superior and inferior meridians. However, the horizontal and vertical meridians differed by around 14% (T-test, p<0.0001). Cone density, expressed in units of area, decreased as a function of axial length (r² = 0.60), but remained constant (r² = 0.05) when cone density is expressed in terms of visual angle supporting the hypothesis that the retina is stretched during the elongation of the eyeball. Gender did not modify the cone distribution. Cone density was slightly modified by age but only at 2°. The older group showed a smaller density (7%). Cone spacing increased from 6,49 ± 0,42 μm to 8,72 ± 0,45 μm respectively between 2° and 6° of eccentricity. The mosaic of the retina is mainly triangularly arranged (i.e. cells with 5 to 7 neighbors) from 2° to 6°. Around half of the cells had 6 neighbors.

Agreement in Cone Density Derived from Gaze-Directed Single Images Versus Wide-Field Montage Using Adaptive Optics Flood Illumination Ophthalmoscopy

Purpose

We compared cone density measurements derived from the center of gaze-directed single images with reconstructed wide-field montages using the rtx1 adaptive optics (AO) retinal camera.

Methods

A total of 29 eyes from 29 healthy subjects were imaged with the rtx1 camera. Of 20 overlapping AO images acquired, 12 (at 3.2°, 5°, and 7°) were used for calculating gaze-directed cone densities. Wide-field AO montages were reconstructed and cone densities were measured at the corresponding 12 loci as determined by field projection relative to the foveal center aligned to the foveal dip on optical coherence tomography. Limits of agreement in cone density measurement between single AO images and wide-field AO montages were calculated.

Results

Cone density measurements failed in 1 or more gaze directions or retinal loci in up to 58% and 33% of the subjects using single AO images or wide-field AO montage, respectively. Although there were no significant overall differences between cone densities derived from single AO images and wide-field AO montages at any of the 12 gazes and locations (P = 0.01-0.65), the limits of agreement between the two methods ranged from as narrow as -2200 to +2600, to as wide as -4200 to +3800 cones/mm².

Conclusions

Cone density measurement using the rtx1 AO camera is feasible using both methods. Local variation in image quality and altered visibility of cones after generating montages may contribute to the discrepancies.

Translational Relevance

Cone densities from single AO images are not interchangeable with wide-field montage derived-measurements.

Perifoveal interdigitation zone loss in hydroxychloroquine toxicity leads to subclinical bull's eye lesion appearance on near-infrared reflectance imaging

Purpose

To characterize the ultrastructural and functional correlates of hydroxychloroquine (HCQ)-induced subclinical bull’s eye lesion seen on near-infrared reflectance (NIR) imaging.

Methods

An asymptomatic 54-year-old male taking HCQ presented with paracentral ring-like scotoma, abnormal multifocal electroretinography (mfERG) and preserved ellipsoid zone on optical coherence tomography (OCT). Dense raster OCT was performed to create en face reflectivity maps of the interdigitation zone. Macular Integrity Assessment (MAIA) microperimetry and mfERG findings were compared with NIR imaging, en face OCT, retinal thickness profiles and wave-guiding cone density maps derived from flood-illumination adaptive optics (AO) retinal photography.

Results

The bull’s eye lesion is an oval annular zone of increased reflectivity on NIR with an outer diameter of 1450 µm. This region corresponds exactly to an area of preserved interdigitation zone reflectivity in en face OCT images and of normal cone density on AO imaging. Immediately surrounding the bull’s eye lesion is an annular zone (3°–12° eccentricity) of depressed retinal sensitivity on MAIA and reduced amplitude density on mfERG. Wave-guiding cone density at 2° temporal was 25,400 per mm². This declined rapidly to 12,900 and 1200 per mm² at 3° and 4°.

Conclusion

Multimodal imaging illustrated pathology in the area surrounding the NIR bull’s eye, characterized by reduced reflectance, wave-guiding cone density and retinal function. Further studies are required to investigate whether the bull’s eye on NIR imaging and en face OCT is prominent or consistent enough for diagnostic use.

Electronic supplementary material

The online version of this article (doi:10.1007/s10633-017-9615-9) contains supplementary material, which is available to authorized users.

Photoreceptor-Based Biomarkers in AOSLO Retinal Imaging

Improved understanding of the mechanisms underlying inherited retinal degenerations has created the possibility of developing much needed treatments for these relentless, blinding diseases. However, standard clinical indicators of retinal health (such as visual acuity and visual field sensitivity) are insensitive measures of photoreceptor survival. In many retinal degenerations, significant photoreceptor loss must occur before measurable differences in visual function are observed. Thus, there is a recognized need for more sensitive outcome measures to assess therapeutic efficacy as numerous clinical trials are getting underway. Adaptive optics (AO) retinal imaging techniques correct for the monochromatic aberrations of the eye and can be used to provide nearly diffraction-limited images of the retina. Many groups routinely are using AO imaging tools to obtain in vivo images of the rod and cone photoreceptor mosaic, and it now is possible to monitor photoreceptor structure over time with single cell resolution. Highlighting recent work using AO scanning light ophthalmoscopy (AOSLO) across a range of patient populations, we review the development of photoreceptor-based metrics (e.g., density/geometry, reflectivity, and size) as candidate biomarkers. Going forward, there is a need for further development of automated tools and normative databases, with the latter facilitating the comparison of data sets across research groups and devices. Ongoing and future clinical trials for inherited retinal diseases will benefit from the improved resolution and sensitivity that multimodal AO retinal imaging affords to evaluate safety and efficacy of emerging therapies.

Parafoveal cone abnormalities and recovery on adaptive optics in posterior uveitis

Purpose

To determine if adaptive optics (AO) flood illumination imaging can detect subclinical changes in 4 cases of posterior uveitis affecting the outer retina.

Observations

In all 4 cases, the affected eye showed altered areas in the photoreceptor mosaic on AO that corresponded to changes on other imaging modalities. Abnormalities not apparent on other imaging modalities were also noted. In one case of multifocal choroiditis with acute outer retinal atrophy, AO revealed decreased visualization of photoreceptors in the unaffected eye that was not noted on spectral domain-optical coherence tomography. In the patient with multiple evanescent white dot syndrome, focal photoreceptor abnormalities were more apparent on AO compared to other imaging modalities, and these areas normalized on AO during follow-up. Five weeks after initiation of high dose prednisone and azathioprine in a patient with serpiginous choroidopathy, AO images showed recovery in apparent parafoveal cone density.

Conclusions and importance

AO detects subclinical changes in the photoreceptor layer in posterior uveitis that can recover over time. AO may be useful in following outer retinal inflammatory conditions.

Hessian-LoG filtering for enhancement and detection of photoreceptor cells in adaptive optics retinal images

Automated analysis of retinal images plays a vital role in the examination, diagnosis, and prognosis of healthy and pathological retinas. Retinal disorders and the associated visual loss can be interpreted via quantitative correlations, based on measurements of photoreceptor loss. Therefore, it is important to develop reliable tools for identification of photoreceptor cells. In this paper, an automated algorithm is proposed, based on the use of the Hessian-Laplacian of Gaussian filter, which allows enhancement and detection of photoreceptor cells. The performance of the proposed technique is evaluated on both synthetic and high-resolution retinal images, in terms of packing density. The results on the synthetic data were compared against ground truth as well as cone counts obtained by the Li and Roorda algorithm. For the synthetic datasets, our method showed an average detection accuracy of 98.8%, compared to 93.9% for the Li and Roorda approach. The packing density estimates calculated on the retinal datasets were validated against manual counts and the results obtained by a proprietary software from Imagine Eyes and the Li and Roorda algorithm. Among the tested methods, the proposed approach showed the closest agreement with manual counting.