Quantitative assessment of human donor corneal endothelium with Gabor domain optical coherence microscopy.

We report on a pathway for Gabor domain optical coherence microscopy (GD-OCM)-based metrology to assess the donor’s corneal endothelial layers ex vivo. Six corneas from the Lions Eye Bank at Albany and Rochester were imaged with GD-OCM. The raw 3-D images of the curved corneas were flattened using custom software to enhance the 2-D visualization of endothelial cells (ECs); then the ECs within a circle of 500-μm-diameter were analyzed using a custom corner method and a cell counting plugin in ImageJ. The EC number, EC area, endothelial cell density (ECD), and polymegethism (CV) were quantified in five different locations for each cornea. The robustness of the method (defined as the repeatability of measurement together with interoperator variability) was evaluated by independently repeating the entire ECD measurement procedure six times by three different examiners. The results from the six corneas show that the current modality reproduces the ECDs with a standard deviation of 2.3% of the mean ECD in every location, whereas the mean ECD across five locations varies by 5.1%. The resolution and imaging area provided through the use of GD-OCM may help to ultimately better assess the quality of donor corneas in transplantation.

Introduction

The corneal endothelium is the innermost layer of the cornea, which serves as a leaky barrier to aqueous humor flow that supplies the stroma with necessary nutrients. Corneal transparency is maintained by healthy endothelial cells (ECs) that contain fluid pumps that control the level of stromal hydration. A malfunction in the deturgescence state of the cornea can lead to corneal edema and eventually to blindness. From penetrating keratoplasty, first developed in 1905, to endothelial keratoplasty,, corneal transplantation remains the main method of treatment for endothelial dysfunction and failure. Endothelial cell density (ECD) is a key determinant in tissue selection/placement and graft survival.

Noncontact specular microscopy (SM) is a standard imaging apparatus used in eye banks to assess the donor corneal endothelium. Although the imaging modality is user-friendly, cost-effective, and fast, most SMs have a small field of view (e.g., Konan EB-10 of ), which introduces sampling error in ECD evaluation. Furthermore, given the curvature of the cornea, SM images often provide a partial view of ECs, which further limits the area of interest for EC analysis and a mixture of proximate layers near the endothelium. The latter issue has been mitigated using the confocal principle, found in various in vivo confocal microscopy (IVCM) such as tandem scanning confocal microscopy,, slit scanning confocal microscopy,, and laser scanning confocal microscopy. More thorough explanations of the utilized techniques and the limitations of various IVCMs for imaging the cornea were reviewed in previous articles,, and IVCMs showed an ECD analysis comparable to SM., Optical coherence tomography (OCT) has been actively adopted in various ophthalmic applications to visualize cross-sectional views of the sample and is used at eye banks for corneal pachymetry. For the histological study of a cornea (e.g., endothelial cells), an embodiment of OCT with a high numerical aperture (NA) microscope objective called optical coherence microscopy (OCM) is gaining attention from researchers for its improved spatial resolution relative to that of IVCM. OCM has steadily evolved into two platforms, depending on the type of illumination and operating principles: full-field OCM (FF-OCM) and Fourier-domain OCM (FD-OCM). FF-OCM is a time-domain OCT, which provides direct access to en face images of the tissue. The volumetric image is obtained by scanning in depth and stacking all en face images from different depths. Cross-sectional images can be extracted from the volumetric image. On the other hand, FD-OCM is a spectral-domain OCT with high-spatial resolution comparable to that of confocal microscopy, which gives direct access to cross-sectional images. En face images can be extracted by restacking the volumetric image. Because of their high-spatial resolution typically obtained with an increase in NA, OCM methods have limited imaging depth that is set by the depth of focus (DOF) of the microscope objective. The DOF is estimated at tens of micrometers (i.e., 40 to depending on the criterion for contrast considered); the curvature of the cornea creates a depth which the ECs lie within . So, in principle, if focused carefully, one could image all ECs across the curved cornea within two zones of the Gabor domain optical coherence microscopy (GD-OCM) using the strictest criterion for DOF. Two or three zones allow the zones to be separated by less than of the DOF of GD-OCM for achieving the highest quality in-focus imaging throughout the endothelium. All imaging was conducted with at least two and up to three zones around the ECs locations for the full field of view of the GD-OCM.

GD-OCM, which is a variation of FD-OCM, was introduced to overcome the DOF limitation. The microscope objective in GD-OCM is custom designed and integrates a liquid lens to effectively extend the depth of imaging beyond the instantaneous DOF of the microscope by scanning the focus of light through depth with an invariant lateral resolution, and acquiring and fusing multiple DOF-limited volumes at different depths for creating an image with high contrast across the extended depth. Fusion of multiple imaging volumes at different depths and high-speed parallel processing provides a 3-D volumetric GD-OCM image with an uncompromised resolution. GD-OCM has successfully revealed the morphologies of human corneas with a isotropic resolution over a millimeter range imaging depth.

In this pilot study, a pathway to assess corneal endothelial cells from ex vivo donor corneas using GD-OCM was demonstrated on six corneas. To assess the variability in ECD across the location of the cornea, the endothelial layer was evaluated at five locations, namely, nasal, temporal, central, superior, and inferior with respect to the apex of the mounted cornea. Four statistical metrics of the corneal endothelium—ECN, ECA, ECD, and CV—were quantified at the five locations. From the image acquisition to cell analysis, the entire process was repeated six times at each location to evaluate the variability in the estimation of ECD.

Methods

Protocol and Sample Preparation

The donor corneas were received from Lions Eye Bank at Albany and Rochester, New York. Corneas were preserved upon recovery in Optisol-GS medium (Bausch and Lomb, Rochester) inside a storage container for corneas and were delivered to the University of Rochester inside a polystyrene box filled with packs of ice. The donor age ranged from 32 to 66 years. Time from death to tissue preservation was about 5 to 22 h. Donor tissue details are summarized in Table 1. Prior to imaging, corneas were stored at room temperature for 1 h to improve the image quality. As a reference, a conventional noncontact specular microscope, Kerato Analyzer (EKA-98, Konan Medical Inc., Japan) was first employed to image the central location of the corneal endothelium through the storage container. A USAF calibration target (R1L1S1P, THORLABS) was used to convert the number of pixels into the physical dimension of endothelial cells area. Corneas were then removed from the storage container and mounted onto a corneal artificial chamber (Moria, Inc., France) equipped with a microfluidic system that was used to control the physiological pressure inside the anterior chamber. The latter was perfused using Optisol to avoid air bubbles inside the anterior chamber.

Table 1

Information on the donor corneas.

Cornea ID	Age	Gender	Cause of death	Death to preservation time (h)	Death to imaging time (h)
1	32	Male	Motocycle accident	22	28
2	59	Female	Chest bleed	9	95
3	59	″	″	″	119
4	60	Male	Nontraumatic intracranial hemorrhage	14	20
5	66	Female	Acute respiratory failure	5	12
6	66	″	″	″	54

Imaging with GD-OCM

The light source of GD-OCM had a central wavelength of 840 nm and a bandwidth of 100 nm. The FWHM of the axial PSF in air was up to 1.2-mm-depth, equivalent to in the cornea. At least lateral separation could be resolved. The NA of the GD-OCM microscope objective was 0.2 and the theoretical DOF (for ) was . The MTF-driven experimental DOF assuming 20% contrast at was in air.

The assembly was secured on the deck of a tip, tilt and rotation stage (TTR001, THORLABS) mounted on two horizontal - and -linear stages (DTS25, THORLABS) and a vertical -linear stage (MVN80, Newport) as shown in Fig. 1(a). The positioning of the sample relative to the GD-OCM probe was monitored using 2-D scanning in and directions. The air gap of between the GD-OCM probe’s window and the cornea was filled with Optisol for index matching, thus reducing high-reflection imaging artifacts resulting from the corneal anterior surface. For imaging the four near-peripheral endothelial layer (superior, inferior, nasal, and temporal) locations, the cornea platform was tilted by 5 deg (up, down, left, and right, respectively) using the tip, tilt and rotation stage as shown in Fig. 1(b). In all cases, the centering of the cornea on the probe was performed using the 2-D scanning in and directions.

Fig. 1

Illustration of the setup for GD-OCM imaging of the corneal endothelium layer. The configurations for (a) the central location on the cornea configuration and (b) one of the four near-peripheral corneal locations (nasal, temporal, superior, and inferior).

The effective number of focal planes allocated for imaging the curved endothelial layer was two to three depending on the cornea, accounting for the corneal curvature. For estimating the robustness of ECD measurement, each time the cornea platform was realigned mostly along the vertical direction, the GD-OCM images were acquired. Using the GD-OCM 4D™ software (LighTopTech Corp., Rochester, New York), a set of liquid lens voltages corresponding to two or three focal planes around the endothelial layer were configured to partition the effective imaging volume as shown in Fig. 2. The imaging procedure from the sample positioning to the reconfiguration of focal planes was repeated six times at each location by three examiners. We obtained thirty images per cornea (six images per location × five—central, superior, inferior, nasal, and temporal—locations) and a total of one hundred eighty volumetric images for six corneas to analyze.

Fig. 2

(a) A cross-sectional view (B-scan) around the corneal endothelial layer, which is a 3-D fusion of two images by depth from zone 1 to 2 and (b) the corresponding intensity profiles (A-scans) of each zone. Zone 1 (yellow) covers the central region of the endothelial layer and zone 2 (cyan) covers the peripheral region of the endothelial layer.

Flattening the GD-OCM Image

The three-dimensional visualization of the endothelial layer was a curved surface. To easily assess the endothelial cells on the curved surface, the raw GD-OCM images of the endothelial layer were projected onto a single plane using a custom MATLAB® code described in Fig. 3. The imaging flattening code runs on every B-scan image and first finds the depth coordinate of the endothelial layer using the peak intensity of the layer at each A-scan. Using the consecutive coordinates for one B-scan, the endothelium layer is curve-fitted by polynomials as a trial. Then any possible erroneous depth coordinate is automatically found based on the smoothness of the layer and is excluded in the secondary curve-fitting step. After curve-fitting, every A-scan is shifted using the second coordinate to create an image of a flat endothelium layer. This process is simultaneously applied to all B-scans of the volumetric image using the parallel computing toolbox in MATLAB®. The resulting volumetric image is then restacked to extract the wide en face view of the endothelial cells layer [Fig. 4(c)].

Fig. 3

Flowchart of the image flattening code.

Fig. 4

Process of the image flattening for cornea 4 at the center: (a) the raw en face images of the endothelial layer (EN) at four different depths of intervals before applying the flattening algorithm, (b) B-scan images, overlaid with the best-fit line of EN (dotted line), and (c) the en-face image of EN after the flattening process.

Counting the Endothelial Cells

The 2-D images of the flattened endothelial layers from all corneas were circularly cropped around the central region for the cell analysis as shown in Fig. 5(a). The diameter of the analysis area, where the cell borders were reasonably visible, ranged from 500 to depending on the location on the cornea. The smaller diameter of was then chosen for the analysis of all of the corneas to compare the results [Fig. 5(b)]. Thirty en face images of the flattened endothelial layers per cornea (total of 180 images for 6 corneas) were distributed to three blinded examiners. Every corner of the endothelial cells was manually identified by the examiners; then a custom MATLAB® algorithm was used to delineate the cell borders with a rule of connecting each corner to three other adjacent corners, as shown in Fig. 5(c). The corner algorithm often misinterpreted the cell borders in a way that caused the border to cross the center of a cell. This occurred both at all outer cells near the window of analysis and at some inner cells of irregular hexagonality. Three blinded examiners edited the borders using Adobe Illustrator® and generated an image of the cell borders, as shown in Fig. 5(d), after comparing the cropped image of the endothelial layer with the cell borders found with the custom MATLAB® algorithm. For the three examiners, it took about one hour from cropping the image to creating the refined cell borders per image. Finally, the images of the cell borders were analyzed using the “Analyze Particles…” plugin in Fiji (an advanced version of ImageJ). Every single counted cell was labeled with the cell index and the pixelated area was obtained, coupled to the corresponding cell index as shown in Fig. 5(e).

Fig. 5

Process of the cell counting of the central endothelium for cornea 6: (a) 2-D en face image of the flattened endothelial layer, (b) the cropped image of (a) within a circle of -in diameter, (c) the endothelial cell borders found using a custom corner method with the manual corner identifications, (d) the corrected borders in (c), and (e) the cell counting using the “analyze particles…” plugin in Fiji.

Results

The cell counting procedure was applied to all the processed images of the endothelial layers as shown in Fig. 6.

Fig. 6

Counted endothelial cells via the corner method, overlaid with the 2-D image of the endothelial layer (N, nasal; T, temporal; C, center; S, superior; and I, inferior, with respect to the apex of the mounted cornea).

For the ’th sample of the counted cells as shown in Fig. 5(e), , , , and were first evaluated using Eq. (1)–(4), respectively. Then ECN, ECA, ECD, and CV were calculated with six samples as , respectively. The variability in the ECD measurements with six samples was quantified as a ratio of the standard deviation to the mean of the ECD in percent as shown in Eq. (5). The statistics of the donor endothelial cells are summarized in Table 2 and Fig. 7.

Table 2

Statistics () of the endothelial cells with the variability in the ECD measurement.

Cornea ID	Age	Area	ECN (cell)	ECA (×10^−4  mm^2/cell)	ECD (cell/mm^2)	CV	Variability in ECD (%)
1	32	N	424.7±117.3	2.496±0.033	4007±53.17	0.361±0.023	1.327
		T	359.5±21.69	2.435±0.033	4108±55.22	0.357±0.018	1.344
		C	465.5±133.7	2.412±0.033	4147±57.59	0.368±0.016	1.389
		S	439.0±98.88	2.496±0.035	4007±56.78	0.365±0.035	1.417
		I	355.3±56.66	2.485±0.032	4025±51.71	0.383±0.031	1.285
		All	408.8±98.87	2.465±0.047	4059±77.84	0.367±0.025	1.918
2	59	N	379.3±56.59	2.857±0.025	3500±30.30	0.263±0.018	0.866
		T	365.3±24.05	2.825±0.021	3540±26.17	0.266±0.006	0.739
		C	383.0±22.11	2.975±0.022	3362±24.39	0.277±0.009	0.725
		S	398.2±16.14	2.797±0.022	3576±28.38	0.274±0.019	0.793
		I	359.8±49.98	2.783±0.022	3594±28.10	0.291±0.018	0.782
		All	377.1±37.45	2.847±0.073	3514±87.99	0.274±0.017	2.504
3	59	N	291.7±131.9	2.821±0.042	3546±53.32	0.313±0.043	1.504
		T	394.7±67.01	2.622±0.022	3815±31.35	0.307±0.012	0.822
		C	279.5±50.86	2.824±0.035	3541±43.94	0.289±0.034	1.241
		S	298.2±56.72	2.870±0.041	3485±49.83	0.285±0.031	1.430
		I	294.0±65.14	2.835±0.029	3527±36.06	0.293±0.030	1.022
		All	311.6±85.61	2.794±0.095	3583±126.6	0.297±0.031	3.534
4	60	N	310.0±32.23	3.549±0.037	2818±29.44	0.306±0.010	1.045
		T	391.5±30.51	3.597±0.044	2781±34.12	0.320±0.005	1.227
		C	337.8±62.01	3.665±0.048	2729±35.97	0.333±0.025	1.318
		S	356.8±42.02	3.650±0.048	2740±36.20	0.293±0.022	1.321
		I	414.7±35.92	3.523±0.041	2839±33.05	0.312±0.009	1.164
		All	362.2±54.49	3.597±0.069	2781±53.64	0.313±0.020	1.929
5	66	N	244.5±50.16	3.668±0.063	2727±46.68	0.311±0.024	1.712
		T	273.0±75.55	3.943±0.087	2537±56.17	0.298±0.026	2.214
		C	247.8±37.61	3.982±0.080	2512±50.74	0.298±0.038	2.020
		S	232.0±30.65	4.071±0.087	2457±52.60	0.285±0.023	2.141
		I	284.8±54.20	4.039±0.069	2476±42.37	0.293±0.027	1.711
		All	256.4±52.17	3.941±0.163	2542±108.7	0.297±0.028	4.275
6	66	N	235.7±28.32	4.068±0.039	2458±23.52	0.270±0.009	0.957
		T	209.0±32.14	4.560±0.070	2194±33.99	0.275±0.018	1.550
		C	256.8±33.25	4.170±0.072	2399±41.02	0.302±0.032	1.710
		S	272.2±50.46	4.389±0.070	2279±36.63	0.325±0.055	1.607
		I	215.5±36.02	3.968±0.060	2521±38.47	0.310±0.038	1.526
		All	237.8±125.3	4.231±0.227	2370±125.3	0.297±0.038	5.288

Note: N, nasal; T, temporal; C, center; S, superior; and I, inferior with respect to the apex of the mounted cornea. ECDs of corneas 1, 2, and 3 were measured as higher than typical ECDs measured with SM. We present the likely cause for this discrepancy in the discussion.

Fig. 7

Statistics of the donor endothelial cells at five different areas of interest for six corneas: (a) ECN (cell), (b) ECA (), (c) ECD (), and (d) CV. The error bars are the SD of each parameter.

For the purpose of validation, six SM images in total were taken near the central location of corneas 5 and 6 (not exactly at the same location as where the GD-OCM image was taken) then an image with the best quality of contrast was used per cornea to analyze the endothelial cells. The identical procedures of the cell count and parameterization were applied to the SM images. The measured , , and were 51 cells, , and , respectively, for cornea 5; and 56 cells, , and , respectively, for cornea 6. For cornea 6, the fell within the range of the average of all areas. However, for cornea 5, the showed a smaller value than the average of all areas by . The results from six corneas showed that (1) the interoperator variability in the estimation of the ECD via our current modality was for each location and for all locations and (2) the ECD value varied with the location of the assessment, which points to the need for evaluating ECD at multiple locations.

Discussion

We presented a pathway to measure the clinical parameters of cornea endothelial cells, i.e., ECN, ECA, ECD, and CV, with GD-OCM. This metrology showed promise in the assessment with the capability of counting more than 400 cells from a curved cornea (but, depending on the image quality, the number of counted cells can be lower). It is worth noting that the standard deviation for the ECD measurement is small and of the mean ECD value at each location. Furthermore, 3-D visualization of corneal structures with cellular resolution and wide field of view () allows for integrative study of corneal structure across the layers. However, the current metrology has four technical limitations for its practical use in eye banks. First, the current flattening algorithm projects the curved endothelial cells to a plane; therefore, the cellular area measured is smaller than its actual area at the curve. This feature will locally result in increased ECD value depending on the curvature. The future flattening algorithm must account for the curvature dependence in the measurement of the cellular area. Second, a donor cornea had to be mounted onto a corneal artificial chamber for the GD-OCM imaging. GD-OCM is a noninvasive imaging technology but the manner of imaging required a mounted cornea. This limitation originated from the short working distance of the GD-OCM probe, which hindered the ability to image a donor cornea through its cornea storage container. However, this inability has been recently surmounted using the GD-OCM probe with a working distance of 15 mm, developed by LighTopTech Co., Rochester, New York. Third, the area of the cellular analysis did not cover the entire field of view of provided by the GD-OCM. This third issue arose from a lower scattering signal at the periphery of the imaging field of view even though the imaging conditions were optimized using multiple focuses across the curved layer. The low scattering signal from the inclined surface is the common phenomena for all intensity-based imaging techniques like SM and IVCM, but the issue might not have been resolved with the interference imaging technique that has an advantage of weak signal detection, herein, GD-OCM. In general, the scattering intensity highly depends on the incident angle of light. If the incident angles at the periphery could be reduced without flattening a donor cornea by force, e.g., adjusting the chief ray angles at the corneal surface via the nontelecentric illumination, this issue could be mitigated. Finally, the cell counting method took about 1 h for one en-face image regardless of the examiner. We attribute this issue to the image quality that influences the correct identification of a single endothelial cell. A robust cell counting algorithm that is less sensitive to image quality is needed to speed up this modality.