|
|
||||||||
1From Eurolens Research, Department of Optometry and Neuroscience, University of Manchester Institute of Science and Technology, Manchester, United Kingdom; the 2Department of Mathematics, University of Manchester, Manchester, United Kingdom; and the 3Department of Medicine, Manchester Royal Infirmary, Manchester, United Kingdom.
| Abstract |
|---|
|
|
|---|
METHODS. The cornea and collected images of the subbasal nerve plexus of 18 diabetic patients (stratified into mild, moderate, and severe neuropathic groups using conventional clinical measures of neuropathy) and 18 age-matched nondiabetic control subjects were scanned, and a novel mathematical paradigm was applied to quantify the extent of nerve tortuosity, which was termed the tortuosity coefficient (TC).
RESULTS. TC was significantly different between the four clinical groups (F3 = 12.2, P < 0.001). It was significantly greater in the severe neuropathic group than in control subjects (P < 0.003) and in the mild (P < 0.004) and moderate (P < 0.01) neuropathic groups. TC did not correlate significantly with the age (r = -0.003, P > 0.05), duration of diabetes (r = -0.219, P > 0.05), or hemoglobin A1c (HbA1c; r = 0.155, P > 0.05) of diabetic patients.
CONCLUSIONS. Corneal confocal microscopy allows rapid, noninvasive in vivo evaluation of corneal nerve tortuosity. This morphologic abnormality relates to the severity of somatic neuropathy and may reflect an alteration in the degree of degeneration and regeneration in diabetes.
Quantitative sensory tests of thermal and pain perception are proposed to assess small-fiber damage.2 3 4 However, we have recently shown no relationship between quantitative sensory tests and small myelinated or unmyelinated fiber damage and repair.5 Alternative, more accurate measures of nerve fiber damage and repair include nerve biopsy with electron microscopy6 and ex vivo confocal microscopy of skin biopsy specimens,7 but both are invasive procedures.
The cornea represents one of the most densely innervated tissues of the body.8 9 Corneal innervation provides protective and trophic functions10 11 12 13 for corneal repair in relation to disease, trauma, or surgery.14 Defining alterations in the corneal nerves has been limited. We have recently used corneal confocal microscopy to quantify corneal nerve morphology in normal subjects15 and have developed this application to show that alterations in fiber density and branching relate to the severity of somatic neuropathy in diabetic patients.16
Corneal nerves course through the stroma which is composed of collagen and substances such as fibronectin and proteoglycans.17 18 These substances are known to be upregulated in diabetes19 and influence axonal outgrowth and regeneration.20 Much of our knowledge on nerve regeneration is based on experiments after peripheral sciatic nerve crush which have demonstrated increased tortuosity of regenerating nerves particularly in older animals.21 Any direct comparison between a peripheral and cranial nerve must be interpreted with caution, as the regenerative response may differ in the two sites. In the present study, we used corneal in vivo confocal microscopy to quantify corneal nerve tortuosity and relate it to the severity of somatic diabetic neuropathy.
| Methods |
|---|
|
|
|---|
Neuropathic Severity Evaluation
All patients underwent a clinical history and neurologic examination to rule out any other cause of neuropathy or previous corneal abnormality. The neuropathy disability score (NDS) was based on a clinical scoring system obtained from a neurologic examination that defined abnormalities of vibration perception threshold (VPT) using a tuning fork, pin-prick perception, and temperature perception threshold (TPT), as well as the presence or absence of ankle reflexes, producing a score ranging from 0 to 10.22 23 Quantitative vibration and thermal assessment were performed with a sensory evaluator (Computer Aided Sensory Evaluator IV [CASE IV]; WR Medical Electronics Co., Stillwater, MN).24 An electromyogram (EMG) machine (model MS92a; Medelec Ltd., Old Woking Surrey, UK) was used to determine the peroneal motor nerve conduction velocity (PMNCV; in meters per second).
Patients were stratified into those with mild (NDS < 3; VPT < 25; PMNCV > 35), moderate (NDS, 36, VPT 2535; PMNCV, 3035), or severe (NDS > 6; VPT > 35; PMNCV < 30) neuropathy.
Confocal Microscopy
Patients were examined with a corneal in vivo confocal microscope (model P4, Confoscan; Tomey) in accordance with our established protocol.15 16 They were instructed to gaze straight ahead. A fixation target was attached to the chin and head rest to facilitate steady fixation. One eye of each subject was selected at random for examination, and several scans of the entire depth of the central cornea were made to acquire satisfactory images of all corneal layers providing three-dimensional images with a lateral resolution of approximately 1 to 2 µm and final image size of 768 pixels x 576 x 3 pixels. Three good-quality images of the subbasal nerve plexus were available for investigation in all diabetic patients and control subjects. The subbasal nerve plexus layer is of particular relevance for defining neuropathic changes because it is the location of the main nerve plexus that supplies the overlying corneal epithelium. The investigator who examined the cornea with the confocal microscope and who undertook morphometric measurements of the images of Bowmans layer (PK) was masked with respect to the severity of neuropathy in the diabetic patients.
Image Processing
Digital images of the subbasal nerve plexus layer were processed using image-processing software (Scion Corp., Frederick, MD). The red-green-blue (RGB) color images were converted to 8-bit indexed color. A grip pen (Intuos; Wacom Technology Corp., Vancouver, WA) was used to trace manually in black, one nerve fiber at a time, along its axis, by selecting the pencil tool in the Toolbox window. Each image was then thresholded to 255 and the current gray-scale image was converted to binary, by setting pixels that had been highlighted by thresholding to black (255) and all other pixels to white (0), resulting in images where nerve fibers appeared black against a white background. Once the processed images had been saved in TIFF format, the TC was calculated with a computer program function that was created for this purpose (MatLab; The Mathworks, Natick, MA). The average TC was calculated for all nerve fibers, but not nerve branches, in all three images of each subject. When nerve fibers exhibited a branching pattern, then only the thickest branch was considered to be a continuation of the nerve fiber. The width of each branch was calculated by averaging three measurements of the diameter of the nerve branch.
TC Computation
Using the MatLab built-in function "im2double," we converted the image to an array (matrix) of numbers. The elements of the matrix were either zeros (background) or ones (nerve fiber). The coordinates of the nerve were the indices of the "nonzero" entries in the matrix, which were returned by the MatLab built-in function "find." A straight line that connected the end points of the nerve fiber was plotted, and the image was translated and rotated to the origin, to align the straight line with the x-axis.
The computation of TC of corneal nerves was based on the approach presented in a previous study,25 where the researchers proposed a quantitative index for evaluating arterial tortuosity, based on the second differences of the coordinates of the vessel midline.
In the present study we calculated TC for corneal nerves based on a series of simple mathematical calculations. Each corneal nerve was represented as the graph of a function. The derivative of a function f at a point x is a measure of the rate at which that function is changing as (one of) its independent variables change. This corresponds to the slope of the tangent to the graph of the function at that point. If we increase x by a small amount, dx, we can calculate f(x + dx).
We first considered equally spaced points xj on the straight line that connected the ends of the nerve. The approximation of the first derivative is given by the difference of two consecutive points on the nerve, divided by the step size (dx). The second derivative is calculated as the difference of two consecutive values of the first derivative, divided by the step size. The step size dx is the distance between the projections on the x-axis of two consecutive pixels of the nerve fiber. The value of dx is constant and equal to 1 pixel, because the number of columns in the matrix that have a nonzero entry is always the same as the number of x coordinates of the nerve fiber. The following equations give an approximation of the first and second derivatives in the interval (xjxj+1), respectively:
![]() |
![]() |
The first and second derivatives are squared and added. The sum is multiplied by the length of the interval (xj xj+1), to estimate the change in the direction of the nerve, within that interval. The sum of all the values is obtained and the square root taken. Once all the quantities have been computed, TC is calculated by the following formula:
![]() |
To test the validity of this approach, TC was calculated for four simple functions: (1) f(x) = sin(x); (2) f(x) = sin(2x); (3) f(x) = sin(4x); and (4) f(x) = x.
The TCs obtained for each function were (1) 2.5066, (2) 7.9265, (3) 29.2292 and (4) 0, respectively. This analysis verified that higher TCs are obtained for curves of greater tortuosity (frequency), whereas the TC for a straight line [f(x) = x] equals zero. The same TCs are obtained when the graphs of these functions are rotated at various angles, indicating that TC is independent of the angle of the nerve axis.
Statistical Analysis
A univariate analysis of variance (U-ANOVA) was conducted to compare the tortuosity of corneal nerves for the four clinical groups. Where differences within the clinical groups were established at P = 0.05 level, post hoc analysis was conducted using the least significant difference (LSD) test. Spearmans correlation coefficient was computed to test for significant associations between the TC and age, duration of diabetes, and HbA1c of the diabetic patients. Correlation was set to be significant at P = 0.05.
| Results |
|---|
|
|
|---|
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
Corneal nerves have protective and trophic functions in the cornea.10 11 12 13 Anatomically, they extend from the ophthalmic division of the trigeminal nerve through the anterior ciliary nerves entering the middle third of the stroma to form the subbasal epithelial plexus anterior to Bowmans layer and finally innervate the basal and superficial epithelial cell layer. Anatomic and immunohistological studies confirm the presence of catecholaminergic, adrenergic, and primarily nociceptive C fibers.8 9 10 11 12 13 These nerves respond primarily to noxious mechanical, thermal, and chemical stimuli; for example, application of topical capsaicin results in a 70% reduction in corneal nerve fiber density.29 Furthermore, recent studies in mutant mice in which TrkA-the high-affinity receptor for nerve growth factor (NGF) has been inactivated, demonstrate a marked reduction in response to mechanical, thermal, and chemical noxious stimuli and the number of nerve terminals in the cornea.30
After LASIK, the number of subbasal and stromal nerve fiber bundles decreases by 90% and, during the first year, reinnervation occurs but the number remains less than half of that before LASIK.31 These findings are of particular relevance to diabetic patients; Rosenberg et al.28 demonstrated a reduction in corneal nerve bundles and related it to loss of corneal sensation and severity of neuropathy in patients with type 1 diabetes. We have recently refined and extended these observations by demonstrating a significant reduction in corneal nerve fiber density suggestive of enhanced degeneration, together with a reduction in branching, suggestive of limited regeneration, which relates to measures of somatic neuropathy in diabetic patients.16 Corneal epithelial metabolism, cell adhesion, and wound healing depend on adequate corneal innervation.32 This may explain the significantly higher risk of development of postoperative epithelial complications and poorer refractive results in diabetic patients who undergo LASIK.14
The mechanisms governing corneal nerve integrity and hence their structure are potentially complex. In the corneal stroma, physical structures such as collagen, fibronectin, and proteoglycans19 33 as well as a number of growth factors including TGF-ß,34 fibroblast growth factor,35 and NGF36 have been shown to regulate nerve fiber damage and repair. This may be relevant, as many of these factors are upregulated in diabetes.19 The morphologic features of corneal nerve fiber degeneration and regeneration remain to be clearly delineated. However, recent studies have demonstrated a reduction in total nerve fiber and branch density, which has been related to loss of both somatic16 and corneal28 sensation. With regard to regeneration, sciatic nerve crush experiments have demonstrated increased tortuosity of regenerating nerves, particularly in older animals.21 Thus, increased tortuosity may represent a morphologic marker of nerve regeneration. The present work demonstrates increased tortuosity of corneal nerve fibers, which is independent of age, duration of diabetes, or glycemic control in diabetic patients with increasing severity of somatic neuropathy. Caution is advised on the functional and clinical relevance of this finding in relation to corneal sensation, especially with the small number of patients studied. Nevertheless, these observations provide further support for a significant impact of diabetes on corneal nerve integrity.
| Footnotes |
|---|
Submitted for publication June 23, 2003; revised October 2, 2003; accepted October 23, 2003.
Disclosure: P. Kallinikos, None; M. Berhanu, None; C. ODonnell, None; A.J.M. Boulton, None; N. Efron, None; R.A. Malik, 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: Rayaz A. Malik, Department of Medicine, Manchester Royal Infirmary, Oxford Road, Manchester M13 9WL, UK; rayaz.a.malik{at}man.ac.uk.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S. Mehra, M. Tavakoli, P. A. Kallinikos, N. Efron, A. J.M. Boulton, T. Augustine, and R. A. Malik Corneal Confocal Microscopy Detects Early Nerve Regeneration After Pancreas Transplantation in Patients With Type 1 Diabetes Diabetes Care, October 1, 2007; 30(10): 2608 - 2612. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Quattrini, M. Tavakoli, M. Jeziorska, P. Kallinikos, S. Tesfaye, J. Finnigan, A. Marshall, A. J.M. Boulton, N. Efron, and R. A. Malik Surrogate Markers of Small Fiber Damage in Human Diabetic Neuropathy Diabetes, August 1, 2007; 56(8): 2148 - 2154. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Greenstein, M. Tavakoli, M. Mojaddidi, A. Al-Sunni, G. Matfin, and R. A Malik Review: Microvascular complications: evaluation and monitoring relevance to clinical practice, clinical trials, and drug development The British Journal of Diabetes & Vascular Disease, July 1, 2007; 7(4): 166 - 171. [Abstract] [PDF] |
||||
![]() |
M. Tavakoli, P. A. Kallinikos, N. Efron, A. J.M. Boulton, and R. A. Malik Corneal Sensitivity Is Reduced and Relates to the Severity of Neuropathy in Patients With Diabetes Diabetes Care, July 1, 2007; 30(7): 1895 - 1897. [Full Text] [PDF] |
||||
![]() |
R. L. Niederer, D. Perumal, T. Sherwin, and C. N. J. McGhee Corneal Innervation and Cellular Changes after Corneal Transplantation: An In Vivo Confocal Microscopy Study Invest. Ophthalmol. Vis. Sci., February 1, 2007; 48(2): 621 - 626. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Saghizadeh, A. A. Kramerov, J. Tajbakhsh, A. M. Aoki, C. Wang, N.-N. Chai, J. Y. Ljubimova, T. Sasaki, G. Sosne, M. R. J. Carlson, et al. Proteinase and Growth Factor Alterations Revealed by Gene Microarray Analysis of Human Diabetic Corneas Invest. Ophthalmol. Vis. Sci., October 1, 2005; 46(10): 3604 - 3615. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |