Abstract
purpose. To investigate and map the polarizing properties of keratoconus corneas in vitro and to compare the results with those obtained in normal corneas.
methods. Corneal buttons of five keratoconus corneas were investigated by polarization-sensitive optical coherence tomography (PS-OCT). The instrument measures backscattered intensity (conventional OCT), retardation, and (cumulative) slow axis distribution simultaneously. Three-dimensional (3-D) data sets of the polarizing parameters are recorded, and two-dimensional (2-D) cross-sectional images as well as en face images of the distribution of these parameters at the posterior corneal surface are derived. The results are compared to similar maps obtained from normal corneas.
results. Compared with normal corneas, the retardation and slow axis orientation patterns are heavily distorted in keratoconus corneas. Larger areas of increased and decreased retardation can be found in keratoconus corneas, markedly increased retardation (up to >50°) can especially be found near the rim of corneal thinning. Contrary to normal corneas, regions where the slow axis markedly changes with depth (by up to 50°–90°) are observed in keratoconus.
conclusions. The observed changes in the cornea’s birefringence properties indicate a change in the arrangement of collagen fibrils in the corneal stroma associated with keratoconus. PS-OCT may be a useful tool for the study and diagnosis of corneal disease.
The tissue of the human cornea is arranged in several layers. More than 90% of its thickness is made up by the stroma, which consists of approximately 200 lamellae with a thickness of 1.5 to 2.5 μm. Each of these lamellae is composed of parallel collagen fibrils embedded in an optically homogenous ground substance.
1 2 The fibrils in successive lamellae are usually oriented at large, approximately orthogonal angles to each other. X-ray diffraction studies have revealed two preferred directions of fibril orientation in the cornea. Approximately 60% of the fibrils are oriented within the 45° sectors around the inferior–superior and nasal–temporal directions, whereas ∼40% are oriented in the oblique areas between those sectors.
3 4 The regular arrangement of fibrils within each lamella is responsible for the transparency of the tissue,
1 5 whereas the orientation of successive fibril layers is important for establishing the mechanical stability of the cornea.
6 7
Because of the critical dependence of the corneal shape and mechanical stability on the orientation of fibrils in corneal lamellae, diseases like keratoconus, in which corneal shape changes, are associated with changes in the regular fibril arrangement.
4 8 Keratoconus is a common dystrophy of the cornea and is characterized by a progressive thinning and ectasia of the central cornea that causes myopia and irregular astigmatism. With an incidence of approximately 50 to 230 per 100,000 of the population,
9 10 it is a frequent indication for keratoplasty.
Although the knowledge of details of local fibril orientation is based on electron microscopy studies,
1 11 information on structural characteristics of the entire cornea is essentially based on x-ray diffraction methods.
3 4 8 12 13 With this technique, an intense beam of x-rays transmits the cornea and is scattered to produce a diffraction pattern characteristic of the ordered structure of the tissue’s constituents. Details on diameter and orientation of the collagen fibrils, integrated over the transmitted volume, can be obtained. With this method, it has been shown that the angle between the preferred collagen orientations in keratoconus corneas were changed from 90° and 180° to 60° and 120°, indicating structural changes in these corneas.
4 Furthermore, a recent study using high-intensity synchrotron radiation enabling a transverse resolution of 0.25 mm demonstrated a dramatic change of the gross organization of stromal lamellae and an uneven distribution of collagen fibrillar mass in keratoconus corneas.
8 Despite the detailed knowledge on collagen structure that has been provided by x-ray diffraction techniques, this technology has some severe drawbacks: Measurement times are very long (several hours), which, in connection with the ionizing nature of x-rays, prevents the use of this technique for in vivo diagnostic applications. Furthermore, the integral nature of the measurement (being performed in transmission) prevents depth-resolved information.
Another method that can be used to obtain information on the corneal structure is probing the cornea with optical radiation. Although scattering methods can also be used,
14 another effect that can be exploited is the birefringence of the cornea. It has long been known that the corneal stroma is optically birefringent, and several papers on corneal birefringence have been published (for an overview, see Ref.
15 ). Different models on the nature of corneal birefringence have been reported that are partly contradictory (e.g., the cornea was modeled as a uniaxial
16 and a biaxial crystal).
17
In a recent study, we suggested a new model of corneal birefringence.
18 This model is based on detailed studies of corneal polarization properties with polarization-sensitive optical coherence tomography (PS-OCT),
19 and on the knowledge of the collagen fibril orientation derived from x-ray diffraction studies. In modeling corneal birefringence, each lamella of the corneal stroma is regarded as a birefringent plate with its slow axis lying along the direction of the collagen fibrils.
20 The fibrils in successive lamellae are usually oriented at large, approximately orthogonal angles to each other, so that the birefringence properties of successive lamellae cancel each other largely. A slight prevalence of one lamella orientation causes a net retardation and a net optic axis orientation.
This model describes the observation of corneal birefringence with low retardation and a preferred axis near the corneal apex under the assumption that the probing light beam impinges perpendicularly on the corneal surface. Toward the periphery, the beam illuminates the cornea with increasing inclination to the corneal surface’s normal vector (the beam is parallel to the optic axis of the eye). For a nonperpendicular probing beam, the birefringences of successive lamellae no longer compensate each other, and retardation increases toward the periphery. Also the slow axis shows a characteristic pattern with axis orientation varying linearly with azimuth angle.
18
PS-OCT is a variant of OCT
21 that measures, in addition to intensity, the polarization state of the light backscattered by the sample, allowing three-dimensional (3-D) information to be obtained on various polarizing properties of tissue (e.g., birefringence,
22 23 24 diattenuation,
25 26 and depolarization, or polarization scrambling).
27 We developed a PS-OCT method that allows the measurement of three parameters simultaneously: reflectivity, retardation, and birefringent axis orientation.
24 In this study, we used our PS-OCT technique to measure and image the polarizing properties of keratoconus corneas in vitro. We obtained 3-D data sets of reflectivity, retardation, and slow axis orientation and derived cross-sectional images of the birefringence parameters, as well as images showing the distribution of these parameters across the posterior surface of the cornea. We compared these results with those previously obtained in normal corneas and, in this article, discuss implications for diagnostic use of PS-OCT.
PS-OCT imaging was performed with a system previously described.
18 24 The system is based on a classic bulk optics Michelson interferometer and works in time domain. A superluminescent diode emits a low-coherence beam (center wavelength, 828 nm; full-width-at-half-maximum bandwidth, 22.4 nm; coherence length, 13.5 μm; depth resolution in tissue, ∼10 μm) that is vertically polarized before entering the interferometer. A nonpolarizing beam splitter cube splits the beam into sample and reference beams. The polarization plane of the reference beam is rotated by 45° on double passing a quarter-wave plate, providing equal power to both channels of a polarization-sensitive detection unit. The sample beam transmits a quarter-wave plate oriented to have a circularly polarized beam incident on the sample. After recombination of reference and sample beams, the interfering beams are split by a polarizing beam splitter into the horizontal and vertical detection channel, where they are phase sensitively recorded by separate photodetectors.
The cornea was illuminated with a power of ∼600 μW. Between the individual A-scans, the sample was moved by a motorized x–y translation stage to obtain A-scans distributed over the cornea along an evenly spaced x–y grid. In this way, 3-D data sets consisting of 80 × 80 A-scans (with a spacing of 75 μm in x and y directions) were recorded, covering a total volume of 6 mm (x) × 6 mm (y) × 2.5 mm (z, optical depth). A-scan rate was 2 Hz.
To obtain sample reflectivity, retardation δ, and (cumulative) slow axis orientation θ, an algorithm based on a Hilbert transform and on a Jones vector analysis was used. While reflectivity and retardation are obtained from the amplitudes of the two polarization channels, the axis orientation is derived from their phase difference.
18 24 On OCT tomograms, the backscattered signal intensity is displayed on a logarithmic false color scale, δ and θ are displayed on linear false color scales. The unambiguous ranges of δ and θ determination are 0° to 90° and 0° to 180°, respectively.
It should be pointed out that the interpretation of the slow axis data requires some care: (1) The algorithm used to calculate θ causes a 90° change in θ at depth position z, where δ crosses 90° (or multiples of 90°). (2) If several layers with different slow axes are stacked on top of each other, the measured value θ (z), at a given depth position z, is a cumulative, or effective, axis: It is the slow axis of a single layered birefringent plate that would give rise to the same change in polarization state (on double pass of the sample) as the multilayered sample does, down to depth z. To differentiate the axis orientations of individual layers would require that (1) these individual layers be resolved (i.e., are thicker than the coherence length) and (2) an algorithm be used that is based on the information obtained on the overlying layers and propagates the polarization state from layer to layer.
For comparison purposes, corneal thickness maps were recorded before tissue explantation in the in vivo cornea by the Orbscan system (Bausch & Lomb, Rochester, NY). The Orbscan is based on the principle of slit projection. Two scanning slit lamps project 40 calibrated beams onto the cornea, and the light diffusely backscattered from the anterior and posterior corneal surfaces is measured to obtain topographic maps of the two surfaces. A corneal thickness map is obtained by calculating the difference between these two topographic maps.
We used PS-OCT to measure retardation and axis orientations of keratoconus corneas in vitro. In a previous paper,
24 we had shown that our instrument provides results with high precision (SD for δ, ∼0.5°, for θ, ∼1.8°) and good accuracy for δ (maximum deviation, ∼3°−4° at δ = 0° and 90°), although there was a systematic deviation of ∼12° for axis orientation, caused by imperfect polarizing elements in the instrument. Because only relative axis orientations are measured in this work (the absolute orientation of the explant with respect to the in vivo eye was not known) we decided not to correct for the systematic offset. Because low signal intensity can degrade the results of PS-OCT measurements,
29 areas with poor signal quality are displayed in gray in our images. An additional factor that degrades image quality in coherent images is speckle noise, which causes the dotted appearance of the PS-OCT images. Averaging techniques can improve the image smoothness, but with the drawback of reduced resolution. Therefore, we decided to present the raw, unsmoothed images.
The results of our measurements show that the normal cornea has a very characteristic birefringence pattern,
18 showing low retardation near the corneal apex (at perpendicular beam incidence) that increases approximately radially and symmetrically toward the periphery. At peripheral locations, a pronounced increase of retardation with depth was observed. The slow axis also showed a characteristic pattern: Near the corneal apex, a predominant axis can be observed, whereas off the apex the slow axis varied approximately linearly with the azimuth angle (in case of a sampling beam parallel to the vision axis). There was little change of axis orientation with depth. This pattern is explained by a birefringence model based on a stack of thin birefringent lamellae with two preferential, nearly orthogonal fibril orientations, superimposed on a background of lamellae with generally randomly oriented fibrils that had a slightly preferential orientation.
The results of this study have shown that this normal pattern is heavily distorted in keratoconus corneas. Typical findings were asymmetric or totally irregular retardation patterns and locally increased retardation near the rim of corneal thinning. The slow axis pattern typically showed larger areas with a predominating range of axis orientations, areas where the axis changed rapidly in a transverse direction, and areas where the axis changed markedly with depth.
It should be mentioned that the observed corneal birefringence depends on the local inclination of the cornea with respect to the sampling beam. The normal pattern is a result of the cornea’s being composed of individual birefringent lamellae stacked in a deliberate way and of the local inclination of these lamellae. It is not easy to separate these effects. Our goal is to obtain information on how far a change in the structure and/or the arrangement of the lamellae occurs in keratoconus. The change in corneal shape is probably a secondary effect of minor interest for the pathogenesis of keratoconus. Since the corneal explants consist of partly very soft and unstable tissue whose shape can be changed during handling and mounting, we have to take care to compare only those parts of the tissue whose inclination is approximately regular (i.e., we have to avoid areas with indentations or folds; this information is provided by the surface elevation data). Taking this into account, deviations of the keratoconus birefringence patterns from normal subjects are still clearly visible.
It is, therefore, likely that the observed birefringence changes are caused by a change in the lamellar structure in keratoconus corneas. Such a change was reported almost a decade ago by Daxer and Fratzl,
4 who used x-ray scattering methods. Meek et al.
8 recently published results obtained by a similar method with improved lateral resolution. They reported a dramatic change in the organization of stromal lamellae and an uneven distribution of collagen fibrillar mass distribution, especially around the apex of the cone.
Care should be used in direct comparison of x-ray scattering data and birefringence measurements by PS-OCT. Although both x-ray scattering and birefringence depend on the orientation of collagen fibrils, this dependence is different. X-ray scattering shows directly the distribution of collagen fibril orientations within the tissue volume transmitted by the x-ray beam. This is an advantage; however, a drawback is that a depth-resolved measurement is not possible. PS-OCT has the advantage that depth resolution is possible (i.e., changes of the parameters δ and θ with depth can be observed). However, since the individual lamellae cannot be resolved (the depth resolution of our technique is ∼ 10 μm, the lamella thickness ∼ 2 μm), only the birefringence integrated over several lamellae is observed. Because the birefringence orientations of lamellae with orthogonal fibril orientation are also orthogonal, their birefringences cancel each other (at perpendicular beam incidence), and so only an excess of one fibril orientation can be detected, the amount of this excess is proportional to the retardation, and its orientation is equal to the birefringent slow axis. No information is obtained on the orientation of those fibers that cancel each other. This indirect and incomplete information provided by birefringence measurements is certainly a disadvantage of our method (in addition to the sensitivity to local surface inclination).
Keeping this difference in information provided by both methods in mind, comparisons of the results show quite good agreement. In both cases, heavy distortions of the patterns—scattering and birefringence—are observed in keratoconus corneas. Considerable changes in the local orientation of x-ray scattering patterns correspond to changes in the slow axis orientations observed by PS-OCT. X-ray diffraction demonstrates that lamellar axis orientation tends to curve around the cones. This means a predominant axis orientation in these regions, and a predominant axis orientation leads to increased birefringence which is observed near the rim of thinned areas. In addition, an associated axis orientation change can be partly observed
(Figs. 3C 3D) . Finally, keratoconus corneas with locally increased and decreased aligned scattering are observed by the x-ray technique. These areas have a larger and a lower amount of fibrils of a predominant orientation, respectively. In a birefringence pattern, increased and decreased retardation would be expected, respectively, in these areas, and, indeed, such areas can be observed
(Fig. 2D) . To summarize, the results of our study support the hypothesis that considerable changes in collagen fibril arrangement are associated with keratoconus.
Our study has also implications for the use of ophthalmic diagnostic instruments operating with polarized light. In particular, scanning laser polarimetry (SLP) of the retinal nerve fiber layer requires that the corneal birefringence be compensated for.
31 The newest generation of SLP has a variable corneal compensator capable of compensating the individual cornea.
32 However, if the instrument is to be used to image through a keratoconus cornea, the corneal birefringence will vary with scan angle, and the corneal compensation will not work properly. Future studies should determine whether this problem also arises in moderate or early-stage keratoconus.
The PS-OCT instrument used in this study is a slow, time-domain system, therefore it could only be used for in vitro measurements. However, OCT technology has undergone considerable improvement recently. Spectral domain OCT has been shown to have huge sensitivity and speed advantages compared with time domain OCT.
33 34 This enables 3-D imaging of the human retina with a data-acquisition time on the order of a few seconds.
35 36 37 Recently, we developed a spectral domain PS-OCT system for retinal imaging.
38 This instrument acquires a 3-D PS-OCT data set of the human retina in 3 seconds. The modification of this instrument for the purposes of anterior segment imaging will allow in vivo applications of this technology in cornea diagnostics, as well. Advanced evaluation algorithms quantifying the local birefringence changes compared with the normal pattern and taking into account the information on corneal surface shape and inclination have to be developed to simplify the interpretation of the PS-OCT retardation and axis orientation maps.
Because PS-OCT is sensitive to changes in the lamellar structure of the cornea, a potential clinical relevance of the technique might be the detection of abnormal corneal biomechanical properties in vivo. These properties can deviate from the normal state as a result of surgery or disease. Until now, there has been no method of determining the biomechanical properties of the cornea in vivo. Potential uses of such measurements might be a better understanding of laser refractive surgery outcomes and the results of nonablative refractive correction such as conductive keratoplasty (CK), and determining other potential factors in ectasia. Furthermore, it may help to understand the impact of corneal physical properties on intraocular pressure measurement and glaucoma risk.
Supported by the Austrian Science Fund FWF Grant P 16776-N02.
Submitted for publication June 29, 2006; revised January 16 and February 13, 2007; accepted May 25, 2007.
Disclosure:
E. Götzinger, None;
M. Pircher, None;
I. Dejaco-Ruhswurm, None;
S. Kaminski, None;
C. Skorpik, None;
C.K. Hitzenberger, None
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked “
advertisement” in accordance with 18 U.S.C. §1734 solely to indicate this fact.
Corresponding author: Christoph K. Hitzenberger, Center for Biomedical Engineering and Physics, Medical University of Vienna, Währinger Str. 13, A-1090 Vienna, Austria;
christoph.hitzenberger@meduniwien.ac.at.
The authors thank Bernhard Baumann for performing the monitoring measurements on corneal tissue.
MauriceDM. The structure and transparency of the cornea. J Physiol (Lond). 1957;136:263–286.
[CrossRef] [PubMed]DonohueDJ, StoyanovBJ, McCallyRL, FarrellRA. Numerical modeling of the cornea’s lamellar structure and birefringence properties. J Opt Soc Am A. 1995;12:1425–1438.
[CrossRef] NewtonRH, MeekKM. The integration of the corneal limbal fibrils in the human eye. Biophys J. 1998;75:2508–2512.
[CrossRef] [PubMed]DaxerA, FratzlP. Collagen fibril orientation in the human corneal stroma and its implication in keratoconus. Invest Ophthalmol Vis Sci. 1997;38:121–129.
[PubMed]BenedekGB. Theory of transparency of the eye. Appl Opt. 1971;10:459–473.
[CrossRef] [PubMed]NyquistGW. Rheology of the cornea: experimental techniques and results. Exp Eye Res. 1968;7:183–188.
[CrossRef] [PubMed]NashI, GreeneP, FosterS. Comparison of mechanical properties of keratoconus and normal corneas. Exp Eye Res. 1982;35:413–423.
[CrossRef] [PubMed]MeekKM, TuftSJ, HuangY, et al. Changes in collagen orientation and distribution in keratoconus corneas. Invest Ophthalmol Vis Sci. 2005;46:1948–1956.
[CrossRef] [PubMed]KrachmerJH, FederRS, BelinMW. Keratoconus and related non-inflammatory corneal thinning disorders. Surv Ophthalmol. 1984;28:293–322.
[CrossRef] [PubMed]KomaiY, UshikiT. The three-dimensional organization of collagen fibrils in the human cornea and sclera. Invest Ophthalmol Vis Sci. 1991;32:2244–2258.
[PubMed]MeekKM, BlamiresT, ElliotGF, GyiTJ, NaveC. The organization of collagen fibrils in the human corneal stroma: a synchrotron x-ray diffraction study. Curr Eye Res. 1987;6:841–846.
[CrossRef] [PubMed]FratzlP, DaxerA. Structural transformation of collagen fibrils in corneal stroma during drying: an x-ray scattering study. Biophys J. 1993;64:1210–1214.
[CrossRef] [PubMed]McCallyRL, FarrellRA. Structural implications of small-angle light scattering from cornea. Exp Eye Res. 1982;34:99–113.
[CrossRef] [PubMed]BourLJ. Polarized light and the eye.CharmanWN eds. Visual Optics and Instrumentation. 1991;310–325.CRC Press Boca Raton, FL.
StanworthA, NaylorEJ. The polarization optics of the isolated cornea. Br J Ophthalmol. 1950;34:201–211.
[CrossRef] [PubMed]Van BloklandGJ, VerhelstSC. Corneal polarization in the living human eye explained with a biaxial model. J Opt Soc Am A. 1987;4:82–90.
[CrossRef] [PubMed]GötzingerE, PircherM, StickerM, FercherAF, HitzenbergerCK. Measurement and imaging of birefringent properties of the human cornea with phase-resolved, polarization-sensitive optical coherence tomography. J Biomed Opt. 2004;9:94–102.
[CrossRef] [PubMed]De BoerJF, MilnerTE, Van GemertMJC, NelsonJS. Two-dimensional birefringence imaging in biological tissue by polarization-sensitive optical coherence tomography. Opt Lett. 1997;22:934–936.
[CrossRef] [PubMed]DucrosMG, de BoerJF, HuangHE, et al. Polarization sensitive optical coherence tomography of the rabbit eye. IEEE J Sel Top Quantum Electron. 1999;5:1159–1167.
[CrossRef] HuangD, SwansonEA, LinCP, et al. Optical coherence tomography. Science. 1991;254:1178–1181.
[CrossRef] [PubMed]De BoerJF, MilnerTE, NelsonJS. Determination of the depth resolved Stokes parameters of light backscattered from turbid media using polarization sensitive optical coherence tomography. Opt Lett. 1999;24:300–302.
[CrossRef] [PubMed]YaoG, WangLV. Two-dimensional depth-resolved Mueller matrix characterization of biological tissue by optical coherence tomography. Opt Lett. 1999;24:537–539.
[CrossRef] [PubMed]HitzenbergerCK, GoetzingerE, StickerM, PircherM, FercherAF. Measurement and imaging of birefringence and optic axis orientation by phase resolved polarization sensitive optical coherence tomography. Opt Express. 2001;9:780–790.
[CrossRef] [PubMed]TodorovicM, JiaoS, WangLV, StoicaG. Determination of local polarization properties of biological samples in the presence of diattenuation by use of Mueller optical coherence tomography. Opt Lett. 2004;29:2402–2404.
[CrossRef] [PubMed]KempNJ, ZaatariHN, ParkJ, RylanderHG, III, MilnerTE. Form-biattenuance in fibrous tissues measured with polarization-sensitive optical coherence tomography (PS-OCT). Opt Express. 2005;13:4611–4628.
[CrossRef] [PubMed]PircherM, GötzingerE, LeitgebR, SattmannH, FindlO, HitzenbergerCK. Imaging of polarization properties of human retina in vivo with phase resolved transversal PS-OCT. Opt Express. 2004;12:5940–5951.
[CrossRef] [PubMed]ClaerhoutI, BeeleH, Van den AbeeleK, KestelynP. Therapeutic penetrating keratoplasty: clinical outcome and evolution of endothelial cell density. Cornea. 2002;21:637–642.
[CrossRef] [PubMed]EverettMJ, SchoenenbergerK, ColstonBW, Jr, Da SilvaLB. Birefringence characterization of biological tissue by use of optical coherence tomography. Opt Lett. 1998;23:228–230.
[CrossRef] [PubMed]DrexlerW, HitzenbergerCK, BaumgartnerA, FindlO, SattmannH, FercherAF. Investigation of dispersion effects in ocular media by multiple wavelength partial coherence interferometry. Exp Eye Res. 1998;66:25–33.
[CrossRef] [PubMed]GreenfieldDS, KnightonRW, HuangX-R. Effect of corneal polarization axis on assessment of retinal nerve fiber layer thickness by scanning laser polarimetry. Am J Ophthalmol. 2000;129:715–722.
[CrossRef] [PubMed]ZhouQ, WeinrebRN. Individualized compensation of anterior segment birefringence during scanning laser polarimetry. Invest Ophthalmol Vis Sci. 2002;43:2221–2228.
[PubMed]LeitgebRA, HitzenbergerCK, FercherAF. Performance of fourier domain vs. time domain optical coherence tomography. Opt Express. 2003;11:889–894.
[CrossRef] [PubMed]De BoerJF, CenseB, ParkHB, PierceMC, TearneyGJ, BoumaBE. Improved signal-to-noise ratio in spectral-domain compared with time-domain optical coherence tomography. Opt Lett. 2003;28:2067–2069.
[CrossRef] [PubMed]NassifNA, CenseB, ParkBH, et al. In vivo high resolution video rate spectral domain optical coherence tomography. Opt Express. 2004;12:367–376.
[CrossRef] [PubMed]JiaoS, KnightonR, HuangX, GregoryG, PuliafitoCA. Simultaneous acquisition of sectional and fundus ophthalmic images with spectral-domain optical coherence tomography. Opt Express. 2005;13:444–452.
[CrossRef] [PubMed]Schmidt-ErfurthU, LeitgebRA, MichelsS, et al. Three-dimensional ultrahigh resolution optical coherence tomography of macular diseases. Invest Ophthalmol Vis Sci. 2005;46:3393–3402.
[CrossRef] [PubMed]GötzingerE, PircherM, HitzenbergerCK. High speed spectral domain polarization sensitive optical coherence tomography of the human retina. Opt Express. 2005;13:10217–10229.
[CrossRef] [PubMed]