Mathematically modeling the involvement of axons in Leber's hereditary optic neuropathy

Leber's hereditary optic neuropathy (LHON), a mitochondrial disease, has clinical manifestations that reflect the initial preferential involvement of the papillomacular bundle (PMB). The present study seeks to predict the order of axonal loss in LHON optic nerves using the Nerve Fiber Layer Stress Index (NFL-S(I)), which is a novel mathematical model. Optic nerves were obtained postmortem from four molecularly characterized LHON patients with varying degrees of neurodegenerative changes and three age-matched controls. Tissues were cut in cross-section and stained with p-phenylenediamine to visualize myelin. Light microscopic images were captured in 32 regions of each optic nerve. Control and LHON tissues were evaluated by measuring axonal dimensions to generate an axonal diameter distribution map. LHON tissues were further evaluated by determining regions of total axonal depletion. A size gradient was evident in the control optic nerves, with average axonal diameter increasing progressively from the temporal to nasal borders. LHON optic nerves showed an orderly loss of axons, starting inferotemporally, progressing centrally, and sparing the superonasal region until the end. Values generated from the NFL-S(I) equation fit a linear regression curve (R(2) = 0.97; P < 0.001). The quantitative histopathologic data from this study revealed that the PMB is most susceptible in LHON, supporting clinical findings seen early in the course of disease onset. The present study also showed that the subsequent progression of axonal loss within the optic nerve can be predicted precisely with the NFL-S(I) equation. The results presented provided further insight into the pathophysiology of LHON.