The Expanded Bead Size of Corneal C-Nerve Fibers Visualized by Corneal Confocal Microscopy Is Associated with Slow Conduction Velocity of the Peripheral Nerves in Patients with Type 2 Diabetes Mellitus
This study aims to establish the corneal nerve fiber (CNF) morphological alterations in a large cohort of type 2 diabetic patients and to investigate the association between the bead size, a novel parameter representing composite of accumulated mitochondria, glycogen particles, and vesicles in CNF, and the neurophysiological dysfunctions of the peripheral nerves. 162 type 2 diabetic patients and 45 healthy control subjects were studied in detail with a battery of clinical and neurological examinations and corneal confocal microscopy. Compared with controls, patients had abnormal CNF parameters. In particular the patients had reduced density and length of CNF and beading frequency and increased bead size. Alterations in CNF parameters were significant even in patients without neuropathy. The HbA1c levels were tightly associated with the bead size, which was inversely related to the motor and sensory nerve conduction velocity (NCV) and to the distal latency period of the median nerve positively. The CNF density and length positively correlated with the NCV and amplitude. The hyperglycemia-induced expansion of beads in CNF might be a predictor of slow NCV in peripheral nerves in type 2 diabetic patients.
Diabetic neuropathy is a major complication of diabetes that may affect the sensory, autonomic, and motor nerves . Mitochondria play a pivotal role in controlling nerve function by their fusion, fission, and trafficking [2, 3] and their morphological and functional anomalies play a pivotal role in developing diabetic neuropathy [4, 5]. Corneal confocal microscopy (CCM) has been used to diagnose small C-fiber neuropathy [6, 7] and visualizes composites of the accumulated mitochondria, glycogen particles, and vesicles  as a bead . In diabetic mice the presence of collections of small mitochondria and the formation of markedly enlarged mitochondria has been reported by Schmidt et al. . Using IENFs of patients with diabetic neuropathy, Hamid et al.  observed larger mitochondrial signals composed of clustered smaller, potentially dysfunctional mitochondria or of larger predegenerative mitochondria due to a disruption in axonal transport. In small fiber neuropathy patients with relatively preserved IENFs, mitochondrial loss as judged by the significant reduction of mitochondrial oxidative phosphorylation was found . However, there has been no report to associate the bead size (BS) of corneal nerve fibers (CNFs) and dysfunctions of the peripheral nerves in diabetic patients, because the size of an individual bead is too small to be assessed. We already reported that enlarging and smoothing of beads by S-Spline Max algorithm make it feasible to determine the size of an individual bead .
The aim of the present study is to associate the BS in CNFs determined by subsequent image processing with the neurophysiological dysfunctions of the peripheral nerves in patients with type 2 diabetes, elucidating whether the BS of CNFs might predict dysfunctions in the peripheral nerves.
2. Research Design and Methods
162 Japanese patients with type 2 diabetes and 45 healthy control subjects (HbA1c < 5.7% (38.8 mmol/mol), fasting plasma glucose < 5.5 mmol/L or casual postprandial plasma glucose < 7.7 mmol/L) were enrolled in this study at the period between May 2013 and January 2015 at the Ishibashi Clinic, Hiroshima, Japan. Subjects were excluded from the study if they consumed alcohol > one unit/day, had neurological symptoms due to neurological disorders other than diabetic neuropathy or had been diagnosed with the neuropathy due to another cause, or wore contact lenses or had history of refractive surgery or anterior segment trauma. The healthy volunteers were recruited from the general population. Written informed consent was obtained according to the Declaration of Helsinki. The Ethics Committee of the Ishibashi Clinic approved the protocol of the present research.
2.2. Laboratory Data
The HbA1c levels were converted to National Glycohemoglobin Standardization Program (NGSP) units by adding 0.4% to the measured values  and subsequently converted to the International Federation of Clinical Chemistry values by using the equation . The serum creatinine, lipid levels, the urinary albumin/creatinine ratio (ACR), and the estimated glomerular filtration rate (eGFR) were also determined.
2.3. Neurophysiological Examinations
Electrophysiology and nerve conduction velocity (NCV) studies were performed using an electromyography instrument (NeuroPak S1, NIHON KOHDEN, Tokyo, Japan). The motor conduction velocity (MCV, median nerve) and sensory conduction velocity (SCV, ulnar and sural nerve), action potential amplitudes of these three nerves, and the distal latency period of the median nerve in the wrist segment were determined.
The vibration perception threshold (VPT) was measured at the left medial malleolus using a biothesiometer (Biomedical Instruments, Newbury, OH, USA). The average of eight readings has been included. The warm and cold perception thresholds (PTs) at the dorsum of the foot were determined using a thermal stimulator that was controlled by a Peltier element (Intercross-200, Intercross Co., Tokyo, Japan). To assess the cardiovagal function of the autonomic nervous system, the coefficient of variation of R-R intervals () was calculated from the R-R intervals of 200 samples on an electrocardiogram.
2.4. Assessment of Neuropathy
Severity of neuropathy and neurological deficits were assessed using the Neuropathy Disability Score (NDS)  which includes evaluation of vibration, pin prick, and temperature perception as well as the presence or absence of ankle reflexes to establish the severity of neuropathy: NDS 0–2, no neuropathy; NDS 3–5, mild neuropathy; NDS 6–8, moderate neuropathy; and NDS 9-10, severe neuropathy.
All study subjects were examined using Heidelberg Retina Tomograph (HRT III, Heidelberg Engineering, Heidelberg, Germany) . The technicians who were performing CCM examination were blinded to patient groups and the signs and symptoms of patients. The examined eye was anesthetized by instilling one drop of 0.4% benoxinate hydrochloride (Santen Pharmaceutical Co., Osaka, Japan). Comfort Gel (Dr. Mann Pharma, Berlin, Germany) was applied to the lens, and a disposable sterilized TomoCap (Heidelberg Engineering GmbH, Dossenheim, Germany) was mounted on the holder to cover the objective lens. After applying Comfort Gel to the TomoCap, the lens was slowly advanced until the gel touched the cornea. More than 50 images of the subbasal nerve plexus were captured using a section mode, and we analyzed at least 6 high-clarity images per subject for the quantification of the following parameters to define changes in the CNFs: (1) CNF density (CNFD): the total number of major nerve fibers/mm2 of corneal tissue; (2) CNF length (CNFL): the total length of all nerve fibers (mm/mm2); (3) corneal nerve branch density (CNBD): the number of branches emanating from all major nerve trunks/mm2 of corneal tissue; (4) tortuosity grade (TG); (5) frequency/0.1 mm of beading (BF); and (6) BS determined after enlarging 5 times and smoothing the original image of CCM using S-Spline Max algorithm (PhotoZoom Pro 4, Gungle Inc., Tokyo, Japan). The pixel numbers of 120 beads were counted using Photoshop Elements 8.0 (Adobe Systems Inc., San Jose, CA, USA) and averaged.
Except for the TG and BS, all measurements were performed using ImageJ (Texelcraft, Tokyo, Japan); the TG was measured by the grading system of Oliveira-Soto and Efron .
2.6. Statistical Analyses
The post hoc analysis of sample power revealed that, with the use of a one-sided ANOVA (significance of 0.05) and Kruskal-Wallis test, the present study population provided statistical power ranged from 0.82 to 0.99. All statistical analyses were performed using the SPSS medical package (SPSS version 19, Chicago, IL, USA). All values are presented as the mean ± SEM. All data sets were tested for the normality using the Shapiro-Wilk test. For normally distributed variables, the comparisons between controls and total diabetic patients or subgroups stratified by NDS were made using one-way ANOVA (for continuous variables) and the -test (for categorical variables) followed by Bonferroni corrections. For nonnormally distributed variables, the Kruskal-Wallis test was applied with subsequent Mann-Whitney test and Bonferroni corrections. The correlations between the CNF measures and the clinical factors or the results of neurophysiological tests were assessed by the Spearman rank correlation test. The sensitivity and specificity of CCM measures in differentiating between control subjects and patients without neuropathy or between patients without neuropathy and with neuropathy were assessed using receiver operating characteristic (ROC) analysis. was considered significant.
3.1. Characteristics of Healthy Controls and Patients with Type 2 Diabetes
The clinical characteristics and detailed assessment of diabetic neuropathy in diabetic patients and controls are summarized in Table 1. The gender ratio and age were similar between controls, total diabetic patients, and their subgroups. The BMI of all patients and the mild neuropathy subgroup was higher than that in controls. The blood pressure in all patients was higher than that in control subjects. The angiotensin receptor blocker or angiotensin-converting enzyme inhibitor was prescribed more frequently for the total patient group and the moderate neuropathy subgroup than for controls. The HbA1c levels in the total patient group and all diabetic subgroups were higher than those in controls, while no difference was found between diabetic subgroups. The LDL-cholesterol level in the total patient group and a mild neuropathy subgroup was elevated compared with that in controls. Statins were prescribed more frequently for the total diabetic patients than for the controls. The HDL-cholesterol level in the total patient group and patients without neuropathy and with moderate neuropathy was lower than that in controls. The triglycerides level in the total patient group and subgroups with mild or moderate neuropathy was increased compared with that in controls. The ACR in the whole patient group and all diabetic subgroups was significantly higher than controls. There were significant differences in the mean NDS between all subgroups of diabetic patients stratified by NDS.
3.2. Neurological Examinations
Except for the MCV of the median nerve, there existed no difference in the neurophysiological tests between controls and the subgroup without neuropathy. The MCV, amplitude, and distal latency period of the median nerve and the SCV and amplitude of ulnar and sural nerve deteriorated relative to the severity of neuropathy (Table 1). VPT increased relative to the severity of neuropathy. and temperature PTs in the moderate neuropathy subgroup were impaired compared with those in controls or patients without neuropathy (Table 1).
3.3. Corneal Nerves Morphological Parameters in Control Subjects and Type 2 Diabetic Patients
In the cohort of patients without neuropathy, the CNFD, CNFL, CNBD, and BF were markedly reduced, and TG and BS were increased compared with controls. The further mild deterioration of CNFD and CNFL was found in the moderate neuropathy subgroup and BS in patients with severe neuropathy (Table 2). BS in patients without neuropathy expanded compared with that of controls and further expansion was found only in patients with severe neuropathy (Figure 1).
Figure 2 illustrated the representative beads in original CCM images (a), when the images were simply magnified 5 times by Photoshop (b), and when the images were enlarged 5 times with smoothing by the S-Spline Max algorithm (c) in a control subject (1), a patient without diabetic neuropathy (2), a patient with mild diabetic neuropathy (3), a patient with moderate diabetic neuropathy (4), or a patient with severe diabetic neuropathy (5). The images that were simply magnified by Photoshop were not suitable to determine the pixel numbers accurately. In comparison, the images processed by the S-Spline Max algorithm were clearly demarcated, and the pixel numbers were easily determined. Compared with the BS of controls, those in diabetic patients appeared to expand. The intrarater variability (averaged CV) in the measurement of the pixel numbers of bead in 10 control subjects was %.
According to the ROC curves for the four CCM parameters between controls and patients without neuropathy (NDS < 3) (Figure 3 and Table 3), AUC was 0.847 for CNFD, 0.830 for CNFL, 0.883 for BF, and 0.997 for BS. Because the AUC, sensitivity, and specificity of BS are the best among parameters of CCM, it is considered that BS is the most reliable marker of CNF alteration in type 2 diabetes before developing neuropathy (Table 3). On the other hand, AUC between patients without neuropathy (NDS < 3) and with neuropathy (NDS > 3) did not seem to be useful in detecting the presence of neuropathy: AUC was 0.610 for CNFD, 0.585 for CNFL, 0.538 for BF, and 0.602 for BS (Table 3).
3.4. Correlations between CCM Parameters and Clinical Factors or Neurophysiological Tests
There was significant correlation between HbA1c and all CNF parameters except TG and BF. NDS was associated with CNFD and CNFL inversely and with BS positively (Table 4). The BS had the robust inverse association with MCV and the positive association with the distal latency period of the median nerve and inversely related with the SCV of the ulnar and sural nerve. On the other hand, CNFD and CNFL had the modest association with the NCV as well as the amplitudes of the peripheral nerves. The weak but significant relationships between the CNFD, CNFL, and CNBD and or cold PT were found, while BS had the positive relationship with a VPT (Table 4).
Diabetic neuropathy is a major complication of diabetes and leads to significant morbidity and mortality . The chronic hyperglycemia-induced changes in mitochondrial structure, function, and dynamics may lead to the development of diabetic neuropathy [2, 3, 17]. The increase in the number of small mitochondria due to excessive fission in the dorsal root fibers of db/db mouse with diabetic neuropathy has been reported . In a study with Akita diabetic mouse with diabetic neuropathy, sympathetic ganglia show the accumulation of minute mitochondria . In patients with diabetic neuropathy the larger mitochondrial signals composed of clustered smaller, potentially dysfunctional mitochondria or of larger predegenerative mitochondria in IENFs due to a disruption in axonal transport were observed . In the small fiber neuropathy patients in early stage, mitochondrial loss was reported . Since an invasive skin or nerve biopsy is mandatory for the morphological assessment of mitochondrial number and size, it could not be performed as a clinical routine procedure.
To our knowledge, this is the largest cohort of type 2 diabetic patients that have been studied in detail for standard clinical examinations and CCM, and the new parameters including BS have been studied for the first time.
The CNF pathology assessed by CCM has been proposed as a surrogate marker for a small fiber neuropathy in patients with diabetes [6, 7], because the ROC curves for the CCM parameters (CNFD, CNFL, and CNBD) for NDS > 3 revealed good AUC, sensitivity, and specificity  and because the CCM analysis using a conventional and novel algorithm to reconstruct CNF images had revealed the significant correlations between CNFD, CNFL, or CNBD and the NCV, VPT, temperature PT, heart rate variability, or the clinical severity of diabetic neuropathy . The ROC analysis of CCM parameters in the present study revealed that the AUC, sensitivity, and specificity of CCM measures in differentiating between control subjects and patients without neuropathy (NDS < 3) were excellent and among the parameters of CNF morphology, BS had the best AUC, sensitivity, and specificity.
Using the thin section of a freshly isolated human cornea the electron microscopy can identify beads in CNFs as a composite of the accumulated mitochondria, glycogen particles, and vesicles [8, 21] playing some role in controlling CNF functions  with similar ultrastructure to nociceptor terminals. The present study revealed that BF decreased and BS expanded in patients with diabetes compared with control subjects. In the present study BS inversely associated with CNFD (269, ), CNFL (.212, ), and CNBD (229, ), and BF had positive relationship with CNFD (, ), CNFL (, ), and CNBD (, ). Therefore, altered beading structures may have relationship with the loss of main nerve fibers as well as branches. Although the increased irregularity in the periodicity of CNF beading has been reported in streptozotocin- (STZ-) diabetic rat , alteration in the BS and in composition of bead remained to be clarified. The axonal trafficking of mitochondria is controlled by cytoskeletons, motor proteins, and ATP fuel supply . In STZ-diabetic rats, axonal neurofilaments in peripheral nerve were deleted , and sciatic level of kinesin 5B motor protein which is involved in the axonal transport of mitochondria is changed . These results suggested that changes in the distribution of mitochondria in CNFs might occur resulting in altered number and size of beads. However, we could not determine the exact mitochondrial area in bead using mitochondria marker, because the present study was an in vivo human study using CCM. Therefore, the contribution of mitochondrial number and size to the BS remained to be determined. In patients with diabetic neuropathy  and the spontaneously diabetic BB-Wistar-rat , mitochondrial accumulation of glycogen particles was observed. In experimental diabetic animals, the loss of vesicles in presynapse  and autonomic nerve endings  was reported. Therefore, the alterations in components of bead other than mitochondria might influence the density and size of bead.
The present research found the significant relationships between parameters of main CNF and branch and the NCV and amplitudes of the peripheral nerves, heart rate variability, and a cold PT. Although we and others previously reported the preliminary findings of lower BF of CNFs in diabetic patients [30–32], the relationship between the BS and neurophysiological tests had never been investigated. The beads can be clearly demonstrated along CNFs by CCM. However, the individual bead in the original CCM image is too small for assessing its area. The enlarging and smoothing of the original CCM image by S-Spline Max algorithm enable us to determine the size of an individual bead in CNFs in spite of the uncertain accuracy and exactness of this method.
Based on the current results, CCM can detect CNF alteration in patients without clinical evidence of neuropathy compared with those of the controls. According to ROC analysis BS had the largest AUC with the best sensitivity and specificity compared with CNFD, CNFL, and BF. These results indicated that the diabetes-induced expansion of BS and changes in other morphological parameters of CCM might have predictive value for the dysfunctions of the peripheral nerves. The expansion of beads occurred in patients without the neuropathy, getting larger in patients with severe neuropathy. The BS had good negative relationship with MCV and SCV of the peripheral nerves and directly associated with a distal latency period of the median nerve. In the present study the BS correlated strongly to HbA1c levels. In myelinated murine axon, the oxidative stress alters the external morphology and reduces the transport of mitochondria at the nodes of Ranvier. These mitochondrial changes expand from the node of Ranvier bidirectionally . In STZ-diabetic mice, there is a significant increase in synaptic delay compared with control mice, and the accumulation of degenerated mitochondria of presynaptic axon of the neuromuscular junction was observed . However, because there has been no report investigating the morphological changes in mitochondria of CNFs in patients with diabetes, we could not assume that the bead in the CNFs is a surrogate marker of the size of mitochondrial area in the peripheral nerves in diabetic patients.
We acknowledge limitations to the present study, which may affect the interpretation of the results. First, although we measured the BS of CNF in patients with type 2 diabetes, we did not determine whether the number of mitochondria in an expanded bead in CNFs increased or not. Furthermore, alteration in other components of bead could influence BS. Second, we assessed BS after enlarging and smoothing of original bead. However, the accuracy and preciseness of this method were not established. The future improvement in a resolution of CCM apparatus will make it possible to assess the BS directly using an original CCM image. Lastly, the potential bias determining BF and BS by human errors could not be ruled out. The automatic analyzing system for beading as already developed for other corneal morphological parameters of CNFs  would eliminate this type of errors.
In conclusion, the expansion of beads of the CNFs in patients with type 2 diabetes occurred before the development of neuropathy and was related to the slow NCV and the prolonged distal latency period of the peripheral nerves. The BS in the CNFs has a predictive value for developing the slow NCV of the peripheral nerves in patients with type 2 diabetes. However, elucidating the mechanisms of the expansion of bead in CNFs was beyond the scope of the present study.
Fukashi Ishibashi and Mitra Tavakoli are the guarantors of this work and as such had full access to all data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis and interpretation.
This study received no financial support. The authors have reported no competing interests.
Fukashi Ishibashi designed the study, researched data, and wrote the entire paper. Rie Kojima and Miki Taniguchi performed a CCM examination and neurophysiological tests. Aiko Kosaka and Harumi Uetake gathered the clinical and laboratory data and statistically analyzed all data. Mitra Tavakoli advised on the statistical analysis, interpreted the results, and reviewed and revised the whole paper.
L. J. Müller, L. Pels, and G. F. J. M. Vrensen, “Ultrastructural organization of human corneal nerves,” Investigative Ophthalmology and Visual Science, vol. 37, no. 4, pp. 476–488, 1996.View at: Google Scholar
R. E. Schmidt, C. A. Parvin, and K. G. Green, “Synaptic ultrastructural alterations anticipate the development of neuroaxonal dystrophy in sympathetic ganglia of aged and diabetic mice,” Journal of Neuropathology and Experimental Neurology, vol. 67, no. 12, pp. 1166–1186, 2008.View at: Publisher Site | Google Scholar
F. Ishibashi, “Glycemic state has real-time impact on beading size and frequency while requiring several years to influence nerve fibers by corneal confocal microscopy,” Diabetes, vol. 61, supplement 1, p. A148, 2012.View at: Google Scholar
A. Kashiwagi, M. Kasuga, E. Araki et al., “International clinical harmonization of glycated hemoglobin in Japan: from Japan Diabetes Society to National Glycohemoglobin Standardization Program values,” Journal of Diabetes Investigation, vol. 3, no. 1, pp. 39–40, 2012.View at: Publisher Site | Google Scholar
S. Yagihashi, M. Kamijo, and K. Watanabe, “Reduced myelinated fiber size correlates with loss of axonal neurofilaments in peripheral nerve of chronically streptozotocin diabetic rats,” The American Journal of Pathology, vol. 136, no. 6, pp. 1365–1373, 1990.View at: Google Scholar
M. Rahmati, R. Gharakhanlou, M. Movahedin et al., “Treadmill training modifies KIF5B motor protein in the STZ-induced diabetic rat spinal cord and sciatic nerve,” Archives of Iranian Medicine, vol. 18, no. 2, pp. 94–101, 2015.View at: Google Scholar
F. Ishibashi, A. Kawasaki, E. Yamanaka, A. Kosaka, and H. Uetake, “Morphometric features of corneal epithelial basal cells, and their relationship with corneal nerve pathology and clinical factors in patients with type 2 diabetes,” Journal of Diabetes Investigation, vol. 4, no. 5, pp. 492–501, 2013.View at: Publisher Site | Google Scholar
F. Ishibashi, R. Kojima, A. Kawasaki, E. Yamanaka, A. Kosaka, and H. Uetake, “Correlation between sudomotor function, sweat gland duct size and corneal nerve fiber pathology in patients with type 2 diabetes mellitus,” Journal of Diabetes Investigation, vol. 5, no. 5, pp. 588–596, 2014.View at: Publisher Site | Google Scholar
M. A. Dabbah, J. Graham, I. N. Petropoulos, M. Tavakoli, and R. A. Malik, “Automatic analysis of diabetic peripheral neuropathy using multi-scale quantitative morphology of nerve fibres in corneal confocal microscopy imaging,” Medical Image Analysis, vol. 15, no. 5, pp. 738–747, 2011.View at: Publisher Site | Google Scholar