The cornea, the outermost window of our visual system, is made of five layers (from the anterior to the posterior cornea: the epithelium, the Bowman’s membrane, the stroma, the Descemet’s membrane (DM), and the endothelium). As a result, the cornea is easily accessible for imaging, although highly vulnerable to various types of infections and diseases—keratoconus (KC), Fuchs endothelial corneal dystrophy (FECD), pseudophakic bullous keratopathy (PBK), and lattice corneal dystrophy (LCD) to cite a few. The extremely large number of corneal diseases with unclear etiology and the almost impossibility of performing biopsies (except for the scraping of the epithelium) explain the need for high-performance imaging in order to optimize the clinical diagnosis. In addition, the cornea is the most commonly transplanted tissue worldwide with transplants each year in the United States alone,1 and objective quality control of corneal grafts in eye banks mainly relies on corneal imaging. For instance, endothelial cell density (ECD) is a commonly used metric for evaluating the endothelial reserve, related to the capacity to maintain corneal transparency with time.
Specular microscopy (SM) has been a standard of care in eye banks and clinical practice for evaluating endothelial cell (EC) health and cell attrition, following various types of intraocular surgery or treatment.2,3 However, this imaging modality is limited to the evaluation of corneal epithelium and endothelium and does not allow the evaluation of all layers of the cornea. Slit-lamp biomicroscopy is a complementary imaging modality to SM, and is routinely used in the eye bank and clinical practice to assess the other layers of the cornea, but its limited resolution does not offer the opportunity to identify cellular structural changes of the cornea. In the 1990s, confocal microscopy (CM) was adopted in clinical settings as a noninvasive imaging modality capable of producing high-resolution images of all corneal layers.45.6.7.–8 One of the limitations of CM, however, is the difficulty in localizing the axial position of the corneal layers under investigation, specifically because, due to the en face orientation of the CM images, real-time cross-sectional visualization of the layers, which can help in sample alignment and in evaluation of the extent of sample thickness, is not available.9 Furthermore, despite significant advancements in montaging multiple images, CM and SM have a limited field of view that introduces sampling errors and significant variability in ECD estimates.1011.–12 Optical coherence tomography (OCT) is an optical imaging technique that has led to rapid and pervasive developments during the past decades and still presents a great potential for the future.1314.–15 OCT has the advantage of providing high depth sectioning of the corneal layers, resulting in large cross sectional and en face views of the tissue, as well as a long penetration depth that enables the simultaneous visualization of all the layers.16 Nevertheless, a limitation of conventional OCT instruments remains the insufficient lateral resolution, which is limited to tens of micrometers and prevents the use of OCT in a wide range of applications that require cellular resolution. This has led to a strong need for combining the capabilities of CM and OCT into a single imaging modality. Recent innovations in OCT have been pushing the boundaries of resolution limit to achieve isotropic cellular resolution that is comparable to that of CM and enables the visualization of individual cells in the cornea. A first step toward this goal was the development of ultrahigh axial resolution spectral-domain OCT (SD-OCT), driven by application fields where fine sectioning of the layers is required.1718.19.–20 The second step toward this goal has been to achieve cellular lateral resolution required to visualize individual cell in en face view of the corneal layers. In this perspective, full-field OCT (FF-OCT) was proposed as a solution to provide high lateral resolution imaging of the cornea.2122.–23 This technique is based on time-domain OCT where the reference mirror is moved to match individual layer depths in the tissue. Despite impressive advancements during the past years in improving the field of view, sensitivity, and speed enabling in vivo imaging,2425.–26 further efforts are needed to acquire images of deep layers of the cornea, as the relative position of the layer, available on cross-sectional SD-OCT images and valuable for sample alignment purpose, is not available in FF-OCT until the data collection and restacking are complete, unless combined with another imaging modality (e.g., SD-OCT), as reported by Xiao et al.27 Optical coherence microscopy (OCM) was introduced to achieve cellular resolution and high-contrast imaging of deep layers in tissue by combining the high axial resolution of SD-OCT with the high lateral resolution of CM gating.28 Bizheva et al.29,30 have recently demonstrated the capability of this technique for in vivo imaging of the cellular structures of healthy and KC human corneas. However, the gain in lateral resolution with OCM is achieved at the expense of a limited depth of focus on the order of 100 to , insufficient to image the full depth of the cornea.31 Gabor-domain (GD) OCM was introduced by our group in 2008 to achieve isotropic volumetric resolution and to advance the capability of OCM by extending the depth of imaging to the millimeter range, well beyond the depth of focus of the microscope.32 A liquid lens embedded in a custom optical design is dynamically driven to refocus at different depths; multiple three-dimensional (3-D) images are then acquired and combined into a single volumetric image.3334.–35 Table 1 summarizes the relevant imaging characteristics of GD-OCM and of other technologies used for corneal imaging—CM, FF-OCT, SD-OCT, and OCM.
Comparison of the corneal imaging technologies. CM, Confocal microscopy; FF-OCT, full-field OCT; SD-OCT, spectral domain optical coherence tomography; OCM, optical coherence microscopy; and GD-OCM, Gabor-domain optical coherence microscopy.
|Technology||Axial resolution (μm)||Lateral resolution (μm)||Imaging depth (mm)||Field of view (mm)||Primary image orientation||Contact/noncontact||In vivo/ex vivo|
|CM8||7.6||1 to 2||En face||Contact||Both|
|SD-OCT16||1-10||10 to 20||6||6 to 16||Cross sectional||Noncontact||In vivo|
|GD-OCMa||2.8||2.2||1.2||Cross sectional||Contact||Ex vivo|
Our prior work related to corneal imaging using GD-OCM revealed that stromal keratocytes in FECD corneas were highly reflective, and larger in size and number compared to normal corneas.36 In the current study, we investigated the capability of GD-OCM to assess morphological modifications occurring at different stages of FECD. The capability of GD-OCM was also extended to assess the morphological changes occurring in PBK, LCD, and KC. Note that, independent from the imaging modality aiming to visualize cellular details of the cornea, the main challenge is to retrieve the en face view of an individual layer of the cornea, due to its curvature. En face images always cut through several layers of the cornea. For instance, it is not possible to visualize the corneal nerve plexus, nor the ECs layer, on a single en face image, which would be necessary for the full evaluation of the layer (e.g., to compute the ECD). Corneal curvature can be overcome by designing a conical scan pattern together with a scanning lens in which all the chief rays are perpendicular to the cornea,37 yielding not only en face visualization of individual layers of the cornea but also an optimal signal collection at the periphery of the field. This approach is particularly interesting for wide field-of-view imaging of the cornea but becomes limited in overcoming the high-order corneal curvature and folding induced by pathological conditions, as it is the case in several corneal diseases. To address the challenge of corneal layers visualization on single en face images, we have also developed in this study a framework to numerically flatten the cornea relative to a targeted layer of interest–the corneal endothelium in this specific case study.
Materials and Methods
Corneal Imaging with Gabor-Domain Optical Coherence Microscopy
GD-OCM (see Fig. 1) is based on spectral-domain OCT, with the particularity of providing a volumetric cellular resolution in the cornea as well as an extended depth of imaging beyond the depth of focus of the objective, necessary to image the full depth of the cornea at cellular resolution. For this study, the light source was a superluminescent diode laser centered at 840 nm with 100 nm FWHM (BroadLighter D-840-HP-I, Superlum, Ireland), which enabled a theoretical axial point spread function (PSF) of in the cornea (average refractive index of 1.387 at 840 nm). The interferometer consisted of a fiber optics coupler (FC850-40-50-APC, Thorlabs) with a 50:50 split ratio. The two outputs of the interferometer are spliced to two custom fiber collimators (LPC-04-840-5/125 with a 3-m-long and 3-mm OD cable, Oz Optics, Canada) to serve as reference and sample arms. The reference arm was equipped with a polarization controller (FPC030, Thorlabs), and a dispersion compensation block (a diffractive grating—GR25-0608, a focusing lens—AC254-050-B, and a flat mirror—PFSQ10-03-G01, Thorlabs) to compensate the dispersion from fiber mismatch between the reference and sample arms, as well as the dispersion from the GD-OCM objective.38 The custom spectrometer was based on the Czerny–Turner model39 and was made of an achromatic doublet—L1 (AC127-019-B-ML, Thorlabs), two concave mirrors—CM1 and CM2 (CM508-200-P01 and CM508-150-P01, Thorlabs), a diffractive grating—DG1 (GR25-1208), a cylindrical lens (LJ1105L2-B, Thorlabs), and a high-speed CMOS line camera (spl4096-140 km, Basler Inc., Exton, Pennsylvania) operating at the line rate of 70 kHz.
The microscope objective probe incorporated an electrowetting liquid lens, which used dynamic focusing to image different depths of the cornea while maintaining a lateral resolution of within the imaging depth of up to .33 For this lateral resolution, the depth of focus of the scanning lens was at any set focus. Consequently, to image all the layers of the cornea at each location, at the full lateral resolution capability of the system, a field of view, corresponding to A-scans, was scanned multiple times with different focal lengths of the liquid lens. Only the desired zone (the confocal parameter of the microscope) in each acquired volume contributed to the final 3-D image.34,35 For high speed imaging, the acquisition of each zone was done simultaneously with the processing and fusing of the precedent zones, using a parallelized multigraphics processing unit architecture.40 A negligible 20 ms delay was allowed between the acquisition of the different zones, corresponding to the automatic adjustment and stability of the liquid lens. The total time for data acquisition was thus equal to the number of zones acquired, multiplied by the acquisition time of each zone (15 s).40 In this study, the average number of zones required to image the full depth of pathologic corneas was eight, yielding a total acquisition time of 2 min. It is worth noting that for a targeted layer of the cornea, a single volume with a unique set focus near the layer of interest could be sufficient, yielding 15 s for the acquisition of a -depth-resolved volumetric image of the depth of interest. During the GD-OCM imaging, the optical contact between the microscope objective and the cornea ( gap) was made using saline. The incident power of the cornea was measured to be 2.2 mW, which is at least three times below the threshold for thermal tissue damage and in vivo eye safety by American National Standards Institute (ANSI).41 Furthermore, the beam focusing on the cornea is widely spread out when it reaches the retina, further reducing any risk of retinal phototoxicity. The 3-D images of a ( A-scans) field of view were captured at the center of the corneas using a custom Python program on Linux platform. The free open-source image-processing platform ImageJ42 was used for analyzing all data and for optimizing the contrast of the images. The contrast enhancement was performed using a volumetric histogram equalization that yielded 3-D images with a full dynamic range of 16 bits. In addition, Voxx,43 the volume rendering program for 3-D microscopy developed at the Indiana Center for Biological Microscopy (Indianapolis, Indiana), was used to visualize the datasets as solid 3-D images and as two-dimensional (2-D) images of the different layers.
Corneal Endothelium Flattening and Analysis
As an OCT-based system, GD-OCM not only generated 2-D cross-section images to evaluate the thickness of the layers but also produced 2-D en face (cut-through) images of each layer, which were essential to evaluate the cellular structures of the layers across a large field of view. As mentioned in the introduction, the curvature of the cornea often compromises the visualization of the entire layer on a single 2-D en face image. In addition, the pathological conditions, such as FECD and KC, induce a significant folding of the endothelium in vivo, which makes the visualization of ECs even more challenging. This constraint is more pronounced on excised human corneas, as the endothelium becomes considerably folded ex vivo due to the loss of intraocular pressure as well as stromal swelling. In this study, we have developed an algorithm to flatten the entire 3-D images of the cornea relative to the endothelial layer. The first step of the algorithm consists of layer segmentation and computation of the 3-D surface of the endothelium. Several approaches for corneal layer segmentation have successfully been developed in the context of the thickness evaluation of corneal layers.4445.–46 However, to the best of our knowledge, these implementations have often been limited to 2-D cross-sectional images and not extended to en face layer visualization. In the proposed framework, the layer of the endothelium was detected for each B-scan of the 3-D volume by applying a modified version of a peak detection algorithm47 to individual A-scans. Figure 2(b) presents the depth profile of the A-scan along the blue line in Fig. 2(a), in which the endothelium is marked with the red star. The algorithm was extended to all the frames of the 3-D image to generate the 3-D surface of the endothelium. A polynomial fit was applied to the resulting peaks of the endothelium to smooth the surface and eliminate eventual artifacts from the peak detection. The red curve in Fig. 2(a) shows the overlap between the endothelium and the resulting polynomial fit on a given frame. The second step of the algorithm consisted of pixel shifting. For each A-scan of the volume, the difference in terms of number of pixel between the corresponding position of the endothelium and the apex of the overall 3-D surface from the polynomial fit was evaluated. The entire A-scan was then circularly shifted by the difference in the pixel number to reconstruct the final 3-D image with a flat endothelium. Figure 2(c) shows the result of the pixel shifting on one frame, in which the endothelium now appears flat. Figures 2(a)–2(i) present the case of a cornea, in which the initial endothelial surface was folded [Figs. 2(a), 2(d), and 2(g)]. The endothelial surface was computed from the initial volume [Figs. 2(e) and 2(h)]. Unlike the initial cut-through en face image [Fig. 2(d)] that includes both the endothelium and the DM, all ECs can be visible on a single 2-D en face image of the endothelium after the flattening [Fig. 2(f)]. Figures 2(j)–2(l) illustrate the case of a FECD cornea, for which the flattening algorithm helped to visualize the en face view of ECs.
One control cornea and eight corneas with four well-characterized diseases were analyzed in this study: FECD, PBK, KC, and type IIIA LCD. FECD is the second most prevalent and well-established corneal dystrophy that affects 4% of the American population over 40 years of age. This condition is associated with an unexplained appearance of drops (guttae, corresponding to excrescences of collagen fibers) at the posterior surface of the DM, yielding the increase in DM thickness, the decrease in ECD, and the loss of the endothelium pump function leading to corneal edema. PBK is a decompensation of the endothelium following complication of cataract surgery, which can also lead to corneal edema and epithelial bullae. KC, the most common corneal dystrophy worldwide that mainly affects teenagers and young adults, is characterized by a progressive thinning of the cornea that affects the corneal shape, with subtle structural changes from a microscopic point of view. Type IIIA LCD is a genetic stromal corneal dystrophy that starts as an accumulation of amyloid deposits, which then branches as linear opacities in and under the Bowman’s membrane and slowly spreads deeper into the stroma. Table 2 presents the characteristics of the corneas investigated.
Characteristics of corneas observed with GD-OCM. The stage of FECD has been evaluated according to the Krachmer grading scale.48
|Disease||Patient/donor age||Disease stage|
|FECD||92||Early stage (Grade 0: No apparent disease. Up to 11 central guttae on each cornea)|
|FECD||69||Moderate stage (Grade 3: A zone of confluent central guttae 2- to 5-mm wide)|
|FECD||81||Advanced stage (Grade 4: A zone of confluent central guttae wide)|
|FECD||78||Latest stage (Grade 5: A zone of confluent central guttae wide plus edema of the corneal stroma and/or corneal epithelium)|
|PBK||79||Advanced stage (endothelium decompensation and corneal swelling)|
|PBK||76||Advanced stage (endothelium decompensation after cataract surgery)|
|KC||34||Stage IV (central opacity)|
All procedures to obtain the corneas conformed to the tenets of the Declaration of Helsinki for biomedical research involving human subjects. Central corneal buttons of patients were removed at the time of treatment with penetrating keratoplasty and were collected as per the usual protocol in place at the University Hospital of Saint Etienne, France. This protocol was written by the hospital’s commission for clinical research and innovation and accepted by the local ethics committee (CPP Sud Est I, CHU Saint Etienne, France, IRB00010220). Normal corneas assigned to scientific use were procured from bodies donated to science (Laboratory of Anatomy, Faculty of Medicine, Jean Monnet University, Saint Etienne, France), as permitted by French law. Each donor volunteered their body and gave written consent to the Laboratory of Anatomy; no additional approval by the ethics committee was required. The disease stage of corneas with FECD and PBK was clinically evaluated using SM, and classified based on the Krachmer grading scale as well as the confluence of guttae. Patients with LCD and KC were evaluated using the clinical tests light microscopy and corneal tomography, respectively. Finally, a normal cornea from a donor was also evaluated using light microscopy after dissection. All corneas were fixed in 4% paraformaldehyde (Sigma, Saint-Quentin Fallavier, France) for 12 h, rinsed in balanced salt solution (Alcon, Rueil Malmaison, France), and shipped to Rochester, New York for experimentations using GD-OCM. Except for the normal cornea chosen as the control, which was fixed 9 h postmortem, the fixation of all corneas was performed within 1 h after retrieval.
Gabor-Domain Optical Coherence Microscopy Performance Evaluation
After system optimization for dispersion and polarization [Fig. 3(a)], the axial resolution of GD-OCM was evaluated using the glass window of the probe as the sample. The mirror of the reference arm was then axially moved to capture the PSF at different depth locations. The average axial PSF of in air (equivalent to in the cornea) was measured across the 1.2-mm depth. The sensitivity roll-off was about two times lower at the depth of 1.2 mm compared to 0.2 mm, which corresponds to about 3-dB sensitivity drop [Fig. 3(b)]. Figure 3(c) highlights the PSF at the depth of about 0.3 mm. The lateral resolution of the system was also evaluated by imaging a standard USAF 1951 target (R1L1S1P, Thorlabs) at the same depth of 0.3 mm. Figure 3(d) presents the en face view of the target imaged with GD-OCM, showing that the system can clearly resolve features of group-7 and element-6, which correspond to a resolution of .
All the corneas were then successfully imaged with GD-OCM. The average thickness of the different layers was computed to evaluate the relative change in thickness at different stages of the disease. Except for the display of the DM and endothelial layer that leveraged the flattening algorithm, all the images are presented in their initial shape that included corneal curvature. As can be seen in Fig. 2(a), the folding of the posterior layers due to FECD, for instance, does not translate to the anterior cornea, and as such, compensating the anterior layers from the endothelium folding will introduce undesirable folding in the anterior layers [Fig. 2(c)]. A similar approach can be used to compensate the curvature or eventual folding of the anterior layers using the epithelial–Bowman’s membrane interface as the target for the flattening algorithm in case one is interested in visualizing the full en face view of corneal epithelial nerve plexus, for instance.49,50
Microstructural Organization of the Control Cornea
A control cornea from a 61-year-old donor was first evaluated using GD-OCM that provided the high-resolution 3-D image of a 1-mm2 central cornea, as well as high contrast cross-sectional and en face images of the different layers of the cornea (see Fig. 4). The anterior [Fig. 4(e)], middle [Fig. 4(f)], and posterior [Fig. 4(g)] stroma represent , 50%, and 90% of the stromal thickness, respectively. These images demonstrate the capacities of GD-OCM to assess the microstructural organization of the different layers of the cornea.51 We believe that the nonhexagonal pavement of ECs [Fig. 4(h)] may have been caused by the fixation. Indeed, Fig. 4(i) presents the case of a fresh (unfixed) ex vivo human cornea imaged with GD-OCM, which displays the hexagonal shape of the ECs. For this comparative study, it was necessary to use a fixed cornea as baseline.
Evaluation of Corneas with Fuchs Endothelial Corneal Dystrophy
The disease stage of all corneas was determined clinically. Grades 1 and 2 FECD were not reported in the present study, as these grades were not indications for penetrating keratoplasty, and such cases did not occur in the donor corneas retrieved during this period.
The first case of FECD evaluated in the study was an excised cornea procured from a 92-year-old male donor. The cornea was initially classified as normal during the routine evaluation prior to GD-OCM imaging, and the donor had no previous history of eye disease. First, a area of the central cornea was imaged using GD-OCM (see Fig. 5). Except for the change in reflectivity and size of the keratocytes located in the posterior stroma [Fig. 5(f)], the rest of the layers appeared normal, including the endothelial–Descemet complex. The elongated shape and hyperactivation of keratocytes in the posterior stroma suggested some abnormalities35 that triggered the need for further investigations. After acquisition of multiple areas, guttae were found only in one localized paracentral area [Fig. 5(h)], suggesting a minor form of FECD.
For the second case, a corneal button with FECD at grade 3 (moderate stage) was excised from a patient (69 years old) during penetrating keratoplasty for FECD and was imaged using GD-OCM (see Fig. 6). The cornea was thicker than the normal cornea, with specifically significant increase in the stroma and DM thicknesses [Fig. 6(a)]. The hyperactivation and elongation of keratocytes appeared to be more pronounced in the middle stroma [Fig. 6(d)], compared to the posterior stroma [Fig. 6(e)]. In addition, the elongation of keratocytes in this case was less pronounced in the posterior stroma, as compared to the early stage of FECD (see Fig. 5). Several guttae were present in the central corneal endothelial–Descemet complex [Figs. 6(f) and 6(g)]. We observed that bright spots at the endothelial surface [Fig. 6(g)] corresponded to guttae covered by ECs. These bright spots were detected in the GD-OCM images due to the high gradient of refractive index between the extremely thin layer of the EC cytoplasm located at the top of the guttae and the guttae themselves. In contrast, a cut-through view of the guttae showed hyporeflective rounded formations surrounded by hyper-reflective bands, known to be cellular bodies of ECs.
The third case of FECD was a corneal button at an advanced stage (see Fig. 7). Stroma and DM thickening was also observed in this case. The elongated shape and hyperactivation of keratocytes appeared to have moved further toward the anterior stroma and was more pronounced in the anterior side of the stroma [Figs. 7(c), 7(d)], compared to the middle [Fig. 7(e)] and posterior [Fig. 7(f)] stroma. Guttae were highly confluent across the central cornea [Fig. 7(g)] and were also easily observed from the 2-D view of the endothelium [Fig. 7(h)].
Finally, a corneal button with FECD at the latest stage was evaluated using GD-OCM (see Fig. 8). Small epithelial bullae as a consequence of advanced edema reaching the epithelium were observed [Fig. 8(b)]. The increased thicknesses of the stroma and DM were also confirmed in this case. An apparent loss in the number of keratocytes was observed in the whole stroma and specifically in the posterior stroma, where keratocytes were almost absent [Fig. 8(f)]. Large guttae were visible on the whole endothelial surface with almost no remaining ECs [Fig. 8(g)].
Evaluation of Corneas with Pseudophakic Bullous Keratopathy
To investigate cases of corneal edema not caused by FECD, we evaluated two cases of PBK (see Fig. 9). In both cases, we observed a high increase in corneal thickness. The elongated shape and hyperactivation of keratocytes appeared to be localized in the posterior stroma [Figs. 9(f) and 9(m)], whereas the anterior [Figs. 9(d) and 9(k)], and middle [Figs. 9(e) and 9(l)] stromal keratocytes appeared normal. In the first case, it was difficult to establish a clear boundary between ECs and DM, which appeared to be thinner than the control and presented filamentous-like structures [Fig. 9(g)], presumably similar to a case of a non-guttate endothelial degeneration, as previously described by Abbott et al.52 and Yuen et al.53 No evidence of guttae was found in this case [Fig. 9(g)]. The second case was an unusual situation, likely an early stage of FECD that was decompensated after cataract surgery. This case revealed an advanced stage of corneal edema with the development of epithelial bullae [Fig. 9(i)]. Guttae-like structures were observed in the endothelium [Fig. 9(n)]. However, no increase in DM thickness, a feature associated with FECD, was observed in these two cases.
Evaluation of a Keratoconus Cornea
GD-OCM image of a KC cornea revealed a significant amount of stromal striae54 (see Fig. 10) and other artifacts caused by the low rigidity of the thin corneal button, which easily deformed and retracted in a Petri dish. Keratocytes were less visible than in the control cornea. ECs were morphologically normal.
Evaluation of a Cornea with Lattice Cornea Dystrophy
A corneal button from a patient with LCD was evaluated using GD-OCM (see Fig. 11). Images revealed lattice branching filaments in the anterior stroma. These filaments appeared as reflective lines creating shadows in the posterior cornea. As compared to stromal striae that are localized in depth, the shadow in this case is starting from the location of the lattice fibers in the anterior stroma and going all the way down to the endothelium. We can also see that the shadow maintains the same shape of the lattice fibers as we go deeper into the stroma, which is not typical of stromal striae.
In this study, we have demonstrated the capabilities of a high volumetric resolution and high-speed GD-OCM prototype in evaluating structural changes characteristic of four well-described corneal diseases. Specifically, even in the case of advanced corneal edema, the system has achieved high contrast imaging of all the layers of the cornea, which is necessary to detect and evaluate microstructural changes that occur during the course of the disease. GD-OCM, combined with the flattening algorithm, provides a wide view of the flat endothelial surface as well as the transition between the DM and the endothelium. The flattening algorithm was essential in suppressing the global curvature or folding of the endothelium surface without compromising the visualization of local discontinuity created by guttae. While preserving the characteristic features of the disease, the flattening of the endothelium helped to highlight features that were difficult to visualize and interpret until then. For instance, until the flattening algorithm was applied, it was difficult to capture the co-localization of ECs and guttae, as it was not possible to visualize all the features on a single en face image [Fig. 2(j)]. ECs and guttae were not clearly visible in some areas of Fig. 2(j), due to the endothelium folding. After flattening, we were able to isolate the full view of ECs [Fig. 6(g)] from the view of guttae [Fig. 6(f)], thus uncovering that the bright spots on Fig. 6(g) were actually ECs lying on top of the guttae [Fig. 6(f)]. We can also clearly delimitate from Fig. 6(f) the area of the confluence of guttae, which was not possible before the flattening of the endothelium.
GD-OCM revealed the microstructural organization of the different corneal layers. An increase in the DM thickness was observed in corneas with FECD, a feature that has been well described in several studies of FECD.55,56 We also observed that the elongated shape and hyperactivation of the stromal keratocytes were persistent across the pool of samples evaluated in this study at different stages of FECD. This modification appeared to start in the posterior layer of the corneal stroma very early in the disease process and progressed toward the anterior corneal stroma during the course of the disease. Because of the proximity of the DM and the posterior stroma, the modification of keratocytes in the posterior stroma in the early onset of FECD is consistent with previous findings supporting that the abnormalities of collagen VIII in the DM were involved in the development of FECD.5758.–59 We also observed that guttae at moderate and advanced stages were still covered by ECs, suggesting that the development of guttae was not complete at these stages of the disease. At the latest stage, we observed a qualitative decrease in keratocytes density, most notably in the posterior stroma, which may be indicative of cell death, as previously observed in the anterior stroma by Hecker et al.60 Furthermore, the change in keratocytes does not seem to be a consequence of corneal swelling, as the hyperactivation of keratocytes starts very early in the disease process, when the endothelium is still capable of maintaining the normal deswelling of the stroma. Future studies may evaluate in vivo the impact of hyperactivation of keratocytes on the outcomes of different FECD treatments (e.g., DSAEK, DMEK).61 It would be interesting to evaluate the preoperative and postoperative state of keratocytes to understand whether endothelial keratoplasty may restore the normal state of the stromal keratocytes over time. One of the key observations of this study is that guttae, as the primary indicator of FECD, can be very difficult to locate at the very early onset of the disease. These results suggest that the hyperactivation of keratocytes in the posterior stroma may potentially be an indicator of FECD at a very early stage when the disease is not yet symptomatic. The clinical evaluation of FECD solely based on the confluence of guttae may not be sufficient for the objective classification of the disease stage.62 Quantitative evaluation of keratocytes at different depths of the stroma may be an alternative.
The DM appeared thin in one case of PBK, likely due to an early decompensation of the low-grade FECD after cataract surgery, indicating that GD-OCM could distinguish between different corneal edema etiologies. Lattice branching filaments in the anterior stroma of LCD were consistent between light macroscopy and GD-OCM images. GD-OCM further highlighted the hyper-reflectivity of these fibers, which induces significant attenuation of light to the posterior cornea. The shadow created by lattice fibers in this case is clearly distinct from the stromal striae. Finally, the thinning of the cornea and normal EC morphology were confirmed in KC.
Although this study provides some interesting insights into the disease progression, there are several limitations that should be addressed to provide the etiological relevance of these findings. Our observations are solely based on structural changes and should be further investigated with additional studies or functional imaging modalities. The feature of early hyper-reflectivity of keratocytes observed in the posterior stroma, before corneal edema, remains to be explored on fresh (unfixed) specimens in order to evaluate the potential impact of the fixation on these modifications. Although it has been shown that fixation can induce shrinkage of tissue and increase in light scattering,63 fixed specimens present the advantages of being as close as possible to the living tissue at the time of fixation and remaining stable over time. They were appropriate for this comparative study as our primary goal was to evaluate the performance of the system in detecting relative changes at different stages of the disease.
As penetrating keratoplasty is becoming a less frequent procedure, at this point of the study, we were able to obtain only one corneal sample at each stage of the disease. The observations in FECD have to be further investigated on multiple samples at each stage of the disease in order to evaluate the consistency and variability among different subjects.
We report on the ability of GD-OCM to identify microstructural changes of the cornea in four common corneal diseases. The evaluation of all layers of the cornea offers the possibility to further discriminate the disease stages and advance our understanding of the disease process. These findings open a wide range of applications for GD-OCM. For instance, the objective evaluation of donor corneas received by eye banks is critical in order to detect abnormalities of the corneal endothelium (ECD, DM thickness, and evidence of guttae) and of the stroma (infraclinic opacity, LASIK interface). Combined with the flattening algorithm, GD-OCM can assess all the corneal layers across a large field of view, including ECs that remain the main means to evaluate corneas prior to transplantation. In particular, GD-OCM can assess ex vivo specimens before fixation and could be complementary to conventional evaluation techniques of pathology by allowing volumetric visualization and restacking in any direction. Furthermore, examination can be performed without fixation and without staining and as such does not interfere with subsequent conventional histology.
Cristina Canavesi and Jannick Rolland are the president and CTO of the start-up LighTopTech Corp.; other authors have no conflict of interest.
This research was supported by the Center for Visual Science NIH Training Grant No. EY007125, the Centre Hospitalier Universitaire de Saint-Etienne, and the National Science Foundation STTR Phase I and SBIR Phase II Grant Nos. IIP-1346453 and IIP-1534701 awarded to LighTopTech Corp. We thank NVIDIA for the donation of GeForce GTX Titan GPUs to LighTopTech Corp. and the University of Rochester. Authors would like to thank Adam Hayes, Andrea Cogliati, and Anand Santhanam for their efforts in targeted aspects of system development.
Patrice Tankam is an assistant professor at Indiana University School of Optometry. He received his PhD in optics in 2010. He had been a postdoctoral fellow at the Institut d’Optique Rhone-Alpes, France, from 2011 to 2012; and a research associate at the Institute of Optics and Center for Visual Science, Rochester, from 2012 to 2016. His research interests include anterior segment imaging, optical metrology, and optical coherence tomography. He is a member of OSA, ARVO, and SPIE.
Zhiguo He received his medical university diploma in 2001 (University of Yan’an, China); his PhD in life sciences in 2011 (Jean Monnet University, France). He then worked as a research fellow in the Laboratory of Biology, Engineering and Imaging of Corneal Graft, Saint-Etienne, France. He has published more than 30 original scientific articles and specializes in imaging, bioengineering, and pathology of human corneas, especially in the corneal endothelium.
Gilles Thuret holds his MD degree in ophthalmology and his PhD in biomedical science. He is a professor at the Université Jean Monnet, Saint Etienne. He is the codirector of the Laboratory of Biology, Engineering and Imaging of Corneal Graft. He is also a cofounder of the start-up Keranova and a member of the Institut Universitaire de France.
Holly B. Hindman specializes in the medical and surgical care of corneal diseases. Her research collaborations have focused on using advanced corneal imaging technology to better understand the corneal diseases and their impacts on ocular optics with the ultimate goal of improving visual outcomes and patient experience. Apart from her career, she has also provided local and national service to the eye banking community.
Cristina Canavesi holds a PhD in optics from the Institute of Optics, University of Rochester, a master’s degree in technical entrepreneurship and management from Simon School of Business, University of Rochester, and an MBA from Simon School of Business, University of Rochester. She was awarded the SPIE Optical Design and Engineering Scholarship in 2012. She is the cofounder and president of LighTopTech, an optics start-up developing high-definition imaging solutions for industrial and medical applications.
Johana Coyoc Escudero holds a BS degree in optical engineering ‘17 from the University of Rochester and is a Paul F. Forman graduate fellow in the Technical Entrepreneurship and Management Master’s Program with a concentration in optics ‘18. Her work experience includes areas such as adaptive optics, virtual reality, free-space optical communications, and synthetic optical holography.
Thierry Lépine received his diploma in engineering in 1986, and his PhD in optics from the Institut d’Optique in 1991. He is a professor of optics at the Institut d’Optique, since 1992. His research and teaching interests include the theory of optical aberrations, optical systems design, and radiometry. He is also an associate professor at the National Observatory of Thailand. Finally, he cofounded two start-ups and he is a consultant for several companies.
Philippe Gain holds his MD degree in ophthalmology and his PhD in biomedical science. He is a professor at the Université Jean Monnet, Saint Etienne. He is the director of the Laboratory of Biology, Engineering and Imaging of Corneal Graft. He is also a cofounder of the start-up Keranova.
Jannick P. Rolland is the Brian J. Thompson Professor of Optical Engineering at the Institute of Optics, University of Rochester and the CTO of LighTopTech Corp., a start-up she cofounded in 2013 in biotechnology. She is the director of the NSF Center for Freeform Optics (CeFO) and the R.E. Hopkins Center for Optical Design and Engineering. She is a fellow of the OSA and SPIE and the recipient of the 2014 OSA David Richardson Medal and the 2017 Hajim Outstanding Faculty Award.