Kinetics of neurodegeneration based on a risk-related biomarker in animal model of glaucoma
© Hayashi et al.; licensee BioMed Central Ltd. 2013
Received: 18 May 2012
Accepted: 4 January 2013
Published: 18 January 2013
Neurodegenerative diseases including Parkinson’s and Alzheimer’s diseases progress slowly and steadily over years or decades. They show significant between-subject variation in progress and clinical symptoms, which makes it difficult to predict the course of long-term disease progression with or without treatments. Recent technical advances in biomarkers have facilitated earlier, preclinical diagnoses of neurodegeneration by measuring or imaging molecules linked to pathogenesis. However, there is no established “biomarker model” by which one can quantitatively predict the progress of neurodegeneration. Here, we show predictability of a model with risk-based kinetics of neurodegeneration, whereby neurodegeneration proceeds as probabilistic events depending on the risk.
We used five experimental glaucomatous animals, known for causality between the increased intraocular pressure (IOP) and neurodegeneration of visual pathways, and repeatedly measured IOP as well as white matter integrity by diffusion tensor imaging (DTI) as a biomarker of axonal degeneration. The IOP in the glaucomatous eye was significantly increased than in normal and was varied across time and animals; thus we tested whether this measurement is useful to predict kinetics of the integrity. Among four kinds of models of neurodegeneration, constant-rate, constant-risk, variable-risk and heterogeneity models, goodness of fit of the model and F-test for model selection showed that the time course of optic nerve integrity was best explained by the variable-risk model, wherein neurodegeneration kinetics is expressed in an exponential function across cumulative risk based on measured IOP. The heterogeneity model with stretched exponential decay function also fit well to the data, but without statistical superiority to the variable-risk model. The variable-risk model also predicted the number of viable axons in the optic nerve, as assessed by immunohistochemistry, which was also confirmed to be correlated with the pre-mortem integrity of the optic nerve. In addition, the variable-risk model identified the disintegrity in the higher-order visual pathways, known to underlie the transsynaptic degeneration in this disease.
These findings indicate that the variable-risk model, using a risk-related biomarker, could predict the spatiotemporal progression of neurodegeneration. This model, virtually equivalent to survival analysis, may allow us to estimate possible effect of neuroprotection in delaying progress of neurodegeneration.
KeywordsDiffusion tensor Neurodegenerative mechanisms Biomarker Kinetic model Glaucoma
Neurodegenerative diseases are slowly progressive and intractable disorders of the nervous system. The temporal profile of neurodegenerative diseases shows a pattern of high variability across patients in terms of symptoms, neuropathology and neuroimaging findings [1–3]. This variability, as well as their slowly progressive nature, poses difficulty in assessing treatment efficacy . For any current and future drug that potentially delays disease progression, it is important to know whether and how much it delays deterioration caused by neurodegenerative diseases. Early detection of disease is now becoming possible by measurement of neurochemical [4, 5] and neuroimaging  biomarkers specifically related to the pathogenic events (for review, see [7, 8]. Although several biomarker models have been proposed to illustrate the relationship between the biomarker cascade and symptomatic/clinical stages [1, 2], there is scarce evidence that biomarkers can be used to quantitatively interpret the neurodegenerative process. To this end, it is necessary to establish a model for the kinetics of neurodegeneration using a biomarker, and to predict the time courses for testing the long-term efficacy of treatment.
Glaucoma, a leading cause of the adult-onset blindness, is known to involve a pathology in which axons of retinal ganglion neurons are mechanically injured by a chronic increase in intra-ocular pressure resulting from dysregulated aqueous fluid circulation . In this disease, neurodegeneration occurs not only in the retinal ganglion neuronal cell bodies and their axons in the optic nerves/tracts, but also transsynaptically in the lateral geniculate nucleus, optic radiation, and visual cortex [14–16]. The causality between increased IOP and glaucomatous neurodegeneration is well established based on several observations: 1) the higher the IOP, the more accelerated is the progression of glaucoma ; 2) therapies that alleviate IOP prevent the progress of the disease ; and 3) experimental animals given a treatment that increases IOP by blocking aqueous-fluid absorption exhibit typical visual-field deficits and retinal pathologies similar to human patients [19, 20]. Although its primary cause is different, the glaucoma shares close similarities with Parkinson’s and Alzheimer’s disease in several pathological findings  and is recently considered as one of neurodegenerative diseases that need neuroprotective therapies .
Here, using diffusion-weighted magnetic resonance imaging, we test whether the causal model explains the time course of neurodegeneration in the living brains of experimental animals with an established neurodegeneration model. The technique of diffusion-weighted magnetic resonance imaging has recently enabled in vivo measurements of white matter microstructure, thus allowing the quantitative, longitudinal assessment of neurodegeneration in various diseases including Alzheimer’s [23, 24] and Parkinson’s diseases [25, 26], amyotrophic lateral sclerosis [27, 28] and glaucoma . By applying the diffusion tensor model , diffusion-weighted data could be used to calculate fractional anisotropy (FA), a scale that expresses anisotropic diffusion motion and is proven to be correlated with the density of viable neuronal axons if conditions permit (for review, see ). We used macaque monkeys (Macaca fascicularis) with glaucoma induced by laser photocoagulation of the trabecular meshwork , which is known to result in a chronic increase in IOP and a pathology that mimics glaucoma, degeneration of axons of retinal ganglion neurons.
The animals with glaucoma showed increased IOP in the affected eye (Figure 1B and Additional file 1: Table S1) during the follow-up period from 33 to 168 days after laser coagulation that blocks aqueous fluid absorption. The mean IOP values in the glaucomatous eye was higher than those in the baseline (54.9 ± 4.2 vs. 23.2 ± 0.8 mmHg, paired-T test, P < 0.005) or than the contralateral (22.8 ± 0.73, T test, P < 0.05, Additional file 1: Table S1). In addition, as has been seen in our previous study , values of IOP in the glaucomatous eye were significantly variable across time (analysis of covariance, F 17,19 = 4.34, P < 0.05) and subject (F 1,19 = 5.35, P < 0.05) (see Figure 1B), which led us to assume that the variability in IOP fluctuates with the rate of progression of degeneration.
Parameter estimates in the constant-rate and risk-based models of primary neurodegeneration using FA values at the optic nerve
F 1,9 = 38.6*, P < 0.005
F 1,10 = 37.9*, P < 0.0005
F 1,8 = 3.5**, N.S.
F 1,9 = 0.6**, N.S.
The relevant statistics were also performed with voxel-based analysis of FA images, because this might help in finding a better model if it identified significant voxels located in known pathways of degeneration, such as the optic nerve or tracts. When both the cumulative risk and the post-operative period were entered as regressors in the statistical model, the FA images revealed a cluster located in the visual pathways, including the optic nerve and tracts, and the sagittal stratum, which were significant determinants in the coefficient of cumulative risk (Figure 1E, Additional file 2: Table S2), but not in the post-operative period. In addition, the decrease in FA was specific to the laser photocoagulation treatment; such a time- or risk-dependent decrease was not found in the FA values of the contralateral optic nerve (Additional file 3: Figure S1). Because the sagittal stratum includes the optic radiation, containing axons of transsynaptic neurons, we considered that the degree of secondary neurodegeneration could be also approximated as a linear function of the degree of primary neurodegeneration (see Appendix 1).
This was tested in detail by analysis with voxel-based regression and model selection. When voxel-based multi-regression analysis was performed using a regressor of the FA values at the primarily degenerated optic nerve, we found significant clusters for the contrast of primary-degeneration FA values at the bilateral optic radiation and the posterior callosum (Additional file 3: Figure S1C, Additional file 4: Table S3). These contain axons of the lateral geniculate and visual cortical neurons, respectively, which are known to be affected by transsynaptic neurodegeneration in glaucoma . These results also depict the pathways to the superior colliculus, which is known to receive direct inputs from retinal neurons  (Additional file 3: Figure S1C, Additional file 4: Table S3). Transsynaptic neurodegeneration expressed by FA values at the optic radiation had a linear relationship with primary neuronal death as measured by FA at the optic nerve (R 2 = 0.82, P < 0.0001, Additional file 5: Table S4), suggesting that this type of pathology has its basis in tight inter-neuronal interactions for cellular survival and supporting the hypothesized protective role of neurotropic factors in transsynaptic neurodegeneration . The kinetics of FA values in the transsynaptic degeneration site were also an exponential function of the cumulative risk (R 2 = 0.85 in ipsilateral and R 2 = 0.80 in contralateral, Additional file 3: Figure S1E, Additional file 5: Table S4) and fitted by the variable-risk model significantly better than by the constant-risk model (F 1,8 = 19.0, P < 0.005 in ipsilateral and F 1,8 = 14.1, P < 0.01 in contralateral, Additional file 5: Table S4). Thus, these findings indicated that transsynaptic degeneration proceeds both in a linear relationship with the primary one, and in an exponential function of the cumulative risk.
Parameter estimates in the risk-based kinetic models of primary neurodegeneration using the L/R ratio of density of NF at the optic nerve
F 1,3 = 22.4, P < 0.05
Using established experimental animals of neurodegeneration and in vivo neuroimaging data, we demonstrated the risk-dependent, stochastic nature of the kinetics of neurodegeneration, in which neurons probabilistically degenerate depending upon the product of time and amount of risk. Plausible exponential decay indicates that the earlier the disease risk is accumulated, the more rapidly neurodegeneration occurs, stressing the need for earlier detection using risk-related biomarkers. The current data showed that the risk-dependent decay spatially extended not only into the optic nerve but also into the remote white matter, consistent with the known widespread pathology of glaucoma in the visual pathways including the optic nerve, lamina cribrosa, lateral geniculate nucleus, and visual cortex [15, 16]. The IOP increase has been known not only to elicit stress on the ganglion neurons and their unmyelinated axons in the retina or/and optic head, but also to induce glial activation and expression of tumor necrosis factor α, which may induce transsynaptic degeneration in remote areas [16, 17]. The kinetics of the secondary transsynaptic neurodegeneration was approximated as a linear function of the kinetics of primary neurodegeneration, making the model effective when assessing degeneration in the brain with its enormous amount of inter-regional connectivity.
The current results provide three important findings as it relates to the kinetics of neurodegeneration. First, the results of our dynamic risk-based model for neurodegeneration confirm the stochastic nature of neurodegeneration in one-hit model, which was originally arisen from knowledge of inherited neurodegenerative diseases [11, 37]. In the inherited trinucleotide repeat diseases, the age of symptomatic onset decreases exponentially with increasing length of trinucleotide repeats . These diseases involve pathologies including initiation of nucleated polymerization, a rate-limiting thermodynamically unfavorable state, followed by rapid irreversible elongation to form fibrils [13, 38]. Because the probability of the initial nucleation is increased by longer trinucleotide repeats, the risk function could simply be approximated as a step function across time; thus, the model explains the length-dependent variability in the age onset [13, 37] and the exponential progress of neurodegeneration .
Second, the results extend the applicability of the risk-based stochastic model to the non-inherited, time or subject-variable neurodegenerative diseases. The variable-risk model better explained the variability of the data by 27% and 9% for quantitative and relative FA respectively (R 2 = 0.94 and 0.98) than did the constant-risk model (R 2 = 0.67 and 0.89, respectively). F-test for these nested models showed significant superiority of the variable-risk model (Table 1). This is important findings which indicate that the neurotoxic risk is not constant across time and subjects, even in typical animal model of neurodegeneration, which is all made by the same procedure. A majority of neurodegenerative diseases, such as Parkinson’s and Alzheimer’s diseases, do not involve trinucleotide repeats in the protein aggregation mechanisms and these diseases show significant between-subject variability in symptoms. This may be due to complicated pathomechanisms, such as the nonlinear nature of the risk [38, 39], or multiple co-risks such as aging [9, 40], microcirculation  and microglial activation . The current result showing the close relationship between pathogenic (IOP) and neurodegenerative (DTI) biomarkers is also consistent with a recent finding in Alzheimer’s disease patients , whereby the deposition of beta-amyloid, a presumed pathogenic factor, was strongly related to hippocampal atrophy in the very early stage. Moreover, the current neurodegeneration model is virtually equivalent to those applied in the survival analysis of the Cox proportional-hazards regression model with time-dependent covariates (see Appendix 2). The model can further be generalized by including multiple risk factors (Appendix 2), such as genetic vulnerability (e.g. Apo E4 allele in Alzheimer’s disease), age, sex, and microcirculation, using multiplicative combination to form the net risk function (Eq. 8). The model can also be used for estimating the effect of neuroprotection by testing the interaction between the intervention and the risk to the time course of neurodegeneration.
Third, our findings of transsynaptic neurodegeneration suggest that the regionally exposed risk can induce system-wide neurodegenerative changes. The IOP increase may only injure the ganglion neurons in the retina, but the risk-based kinetic model allowed us to detect the FA decrease in the optic radiations, the 1st order transsynaptic sites. The observations are well consistent with findings by recent studies in glaucomatous patients [29, 44]. The effect size may not be large in the higher order transsynaptic connections (as shown in Additional file 3: Figure S1); however, a specific spatial pattern covering the primary and secondary connections may increase diagnostic specificity for neurodegeneration that often proceeds in the specific neural circuit, e.g. visual pathways in case of glaucoma.
Although the heterogeneity model with stretched exponential decay (in Eq. 6) did fit well to our data, it did not significantly better explain our data than the risk-based model. The variability of the data was largely explained by the risk-based model (94-95%), which was improved only slightly by 1–2% when applied the heterogeneity model and F-test for model selection did not show significant difference between these two nested models (Table 1). Similarly, previous study showed that the stretched exponential model better fit to the survival curve of neurodegeneration using 16 kinds of data  and improved explanation of variability only by 0–11% (median = 4%), while the (constant-) risk-based model explained 72–98 % of data variability (median = 89%). Therefore, it is possible that the heterogeneity model may better explain the neuronal survival than the risk-based model, but the effect size seems to be relatively small as compared with models with less smaller number of parameters. Moreover, we consider that the model needs to be carefully assessed for its eligibility, particularly, for the parameter, γ. The value of γ may be influenced not only by biological factors (e.g. heterogeneity in the decay rate or late-stage biological compensatory response to the initial neurotoxic events) but also by measurement accuracy of the biomarker. Therefore, the actual factor that governs the variability of γ should be assessed systematically in each experimental condition. As for estimation of heterogeneity of decay rates, bi-exponential decay model could be better suited than the current model to evaluate the level of heterogeneity . Moreover, since adding the parameter (γ) in the model may not only increase goodness of fit but also bias the base parameter (β or r 0 ), it may need to be evaluated particularly when for example the study intervenes early-stage neuroprotection that may expect any change in β.
A few issues suggest the need for some caution when applying the current approach. First, to increase generalizability, the current model needs to be tested in other kinds of neurodegenerative diseases and animal models. While the origin of glaucomatous pathology is chronic physical stress to retinal ganglion cells due to continuously increased intraocular pressure, most of neurodegenerative pathologies in Parkinson’ disease or Alzheimer’s disease is primarily based on the biochemical events, such as abnormal protein aggregations . Therefore future studies should also address potential of the current model in various pathologies of neurodegeneration. Second, the predictability of the current model may depend on the measurements of pathogenic biomarker, such as initial time point or frequency. This effect is difficult to estimate without a priori knowledge of temporal variation in the pathogenic biomarker; thus pilot studies that estimate the potential variation may be required. In terms of diseases with abnormal protein aggregation, this effect may not be large because it takes years and even decades for the accumulation of proteins to occur. Third, the predictability of the model rests on the quality of biomarker data, and this may be different across diseases. For example, in Parkinson’s disease, fluorodopa PET can be used as a neurodegeneration biomarker because the uptake constant of fluorodopa has a direct linear relationship to the number of surviving dopaminergic neurons in the substantia nigra . A potent pathogenic role of microglial activation can be estimated by testing the predictability of the model using PET and 11C-PK11195, a marker of translocator protein [49, 50].
In summary, we have successfully explained the kinetics of glaucomatous neurodegeneration by a dynamic risk-based kinetic model and in vivo measurement of biomarkers. The results support that neurodegeneration occurs and proceeds as a probabilistic events based on the risk to neurons. Although many studies may be preferable in other kinds of neurodegenerative diseases to generalize the current kinetic model, the model may provide opportunities for predicting the time course of neurodegeneration in the living brain in the early stage of the disease, and may bridge the dissociation between heterogenous pathology in patients and standardized pathology in experimental animal models of neurodegeneration .
Kinetic models of neurodegeneration
where N 0 is the number of neurons at time zero. The function for risk, r(t), should be non-negative and could be substituted by an exponential function as follows:
where r 0 constitutes a component of baseline risk, while the other term, exp(β · p(t)), constitutes the risk. In the variable-risk model (Eq. 3), the risk for glaucoma, intraocular pressure (IOP), was incorporated into the p(t), which was calculated as measured IOP values minus 29 mmHg, the 95% confidence upper limit of IOP in the normal optic nerve. The optimal values for r 0 , β , N 0 were estimated to obtain the best fit of the model by a non-linear least square method using the values of IOP and FA in the glaucomatous optic nerves at each time point.
where r 0 is a rate constant for the decrease in the number of neurons, expressed as counts per unit of time. The initial condition was again substituted by N(0) = N 0 .
In this model, γ is a shape parameter such that if 0 < γ < 1, N(t) shows a typical form of stretched exponential decay across ∫ r(t)dt, while if γ = 1 it becomes identical to the original risk-based exponential decay as in Eq. 2, and if γ > 1 sigmoidal (or compressed exponential) decay. Mathematically, the function is also known as a Weibull function .
Eight macaque monkeys (Macaca fascicularis, body weight 4–5 kg, all male) were used; five animals were used to model glaucomatous pathologies, while three were control. The animal model of glaucoma was made using an established method of laser photocoagulation, as described previously . Glaucoma model animals were evaluated in detail for histological changes in microglial cells, astrocytes, and neurons in the lateral geniculate nucleus, and the results were published elsewhere . Visual field defects and pathological features of this animal model are also described in detail in our previous studies [19, 20, 54]. Before being enrolled into these experiments, all animals were confirmed by ophthalmoscopy to have no abnormalities in their ocular fundus. After glaucoma was induced, the animals were repeatedly followed up with measurements of intraocular pressure (IOP) and diffusion tensor imaging. All surgical and experimental procedures conformed to NIH guidelines for the care of experimental animals (National Institutes of Health Committee on Care and Use of Laboratory Animals, 1985). The study was approved by the Institutional Animal Care and Use Committee of RIKEN and the National Cardiovascular Center Research Institute.
Induction of experimental glaucoma
Elevated IOP was induced by applying argon blue/green laser photocoagulation burns to the trabecular meshwork of the left eye, with the right eye being used as an untreated control. For the laser treatment, the animals were anesthetized with an intramuscular injection of ketamine (8.75 mg/kg) plus xylazine (0.5 mg/kg). A single-mirror Goldmann lens filled with a physiological solution was placed on the eye to be treated. The argon laser was focused on the mid-portion of the trabecular meshwork, and a total of 150 laser-beam spots were applied around 360° (spot size 100 μm; power 1.0 W; exposure time 0.2 sec) using an argon laser photo-coagulator (Ultima 2000 SE; Coherent Inc., CA, USA) attached to a standard slit-lamp microscope (BQ 900; Haag-Streit, Köniz, Switzerland). The same laser treatment was repeatedly applied two weeks after the first treatment to maintain continuous IOP elevation, as has been confirmed previously. Potential complications of the glaucoma surgery, including retinal ischemia were not suspected by detailed histological evaluations .
The IOP was intermittently measured (over an interval of 3–28 days) in both eyes during the course of glaucoma in each animal using a calibrated applanation pneumotonometer (Model 30 Classic Pneumotonometer; Medtronic Solan, FL, USA). The measurement was performed under generalized anesthesia using intramuscular ketamine (8.75–10 mg/kg) and local anesthesia with 0.4% oxybuprocaine hydrochloride.
Diffusion tensor magnetic resonance images (DTI) were acquired using a 3-Tesla magnetic resonance image (MRI) scanner (Signa Horizon Lx VH3, General Electric Healthcare, Little Chalfont, Buckinghamshire, UK), which provides a maximum gradient strength of 40 mT m-1, rising in 268 μs using a customized eight-channel phased array receiver coil. During the DTI scanning, animals were intubated and deeply anesthetized using a gaseous anesthetic, 1.5% isoflurane, with their respiration assisted by a ventilator (Cato, Dräger, Germany). To maintain stable physiological conditions, we monitored the partial pressures of oxygen, carbon dioxide and isoflurane in the inspiratory and expiratory gases. We also intermittently performed gas analysis of the arterial blood that was collected from the tail artery, in which a 24-gauge needle was indwelled. The animal’s head was fixated to a customized acrylic retainer that was firmly attached to the receiver coil. The animal and the retainer were placed on the MRI gantry with the animal’s head centered in the bore of MRI scanner. After scanning a localizer for planning a field of view for subsequent scans, a higher-order shimming was performed using a spiral sequence that allows us to calculate the fitted field map including second-order components (for a total of 10 components), and to minimize inhomogeneity in the static magnetic field. Then, the DTI data were collected using a multi-shot spin-echo type EPI sequence (number of shot = 2, TR = 17000 ms, TE = 81.9 ms, FA = 90°, number of slice = 51) with an isotropic spatial resolution of 0.9 mm and with diffusion-weighted gradients of 27 directions (b value = 1000 s mm-2). The diffusion-weighted gradients were applied in a symmetrical, twice-refocused pulse to reduce eddy-current-induced distortion in the DTI images . The DTI data, a total of 24 volumes for three non-diffusion-weighted volumes and 21 diffusion-weighted volumes, were obtained by three separate scans, each consisting of one non-diffusion and seven diffusion-weighted volumes. Before each scan, we monitored the temperature of the gradient coil, and started the scan when the temperature was between 22° and 23°C, at which the shift and distortion across scans were minimized. If the temperature was less than that specified, we performed dummy scans to keep the temperature within range. The DTI data were scanned three times for each time point during the follow-up of glaucoma. The total time for obtaining DTI data was approximately 5 hours, depending on the temperature and computation time for the image reconstruction. We also scanned two gradient-echo sequences with different TE values (TR = 200 ms, TE1 = 4.4 ms, TE2 = 6.637 ms) and obtained the field map to be used for a post-process of reducing distortion of DTI images, which originates from the magnetic field inhomogeneities caused by magnetic susceptibility differences between neighboring tissues, such as air-bone or air-tissue.
Raw DTI data were corrected for image shift, rotation and distortion (with a linear transformation of 12 degrees of freedom) and for B0-inhomogeneity distortion (with non-linear warping calculated on the fieldmap) using the programs FLIRT and PRELUDE, respectively, both part of the Functional magnetic resonance imaging of the brain (FMRIB) Software Library (FSL) , developed by the Analysis Group of FMRIB Centre, University of Oxford. Estimated shifts for these corrections were concurrently applied to the original images to minimize reslicing noise in the corrected data. Then, the image for fractional anisotropy (FA) was computed by fitting a diffusion tensor model  to the corrected diffusion data using the program, FDT, part of the FSL.
For region of interest (ROI) analysis of FA based on the primary degeneration kinetic model, we delineated the ROI with a columnar shape (3 mm long, diameter of 2 mm, number of voxels of 12), centered at the center of the optic nerve, at a distance of 6 mm from the junction of the eyeball and the nerve. For each ROI, we obtained FA values averaged across three scans obtained at the same time point. For fitting the kinetic models of primary neurodegeneration, we used a quantitative value of FA or ratio of FA (rFA, expressed as a ratio to the FA of the contralateral optic nerve) in glaucomatous animals and determined optimum values for N 0 , r 0 and β. When we used rFA, N 0 was fixed to one. The integration of risk function, r(t), in Eq. 1 was calculated based on the trapezoid rule using discrete time point data for IOP. The goodness of fit was evaluated based on coefficients of determination (R 2 ) and the P-value computed in the fitting process. The F-test for model selection was performed for the two nested models between with constant- and variable-risk or between with variable-risk and heterogeneity. The fitting and F-test were computed using Pybld (http://www.mi.med.osaka-u.ac.jp/pybld/pybld.html) built in the language Python. We also confirmed that the goodness of fits in the stochastic variable-risk model was independent of the arbitrarily determined threshold (29 mmHg, 95% confidence upper limit of normal IOP). For presenting graphs, the optimized values determined at the fitting process for FA data, was used for N 0 , r 0 and β to calculate the cumulative risk, ∫exp(r(t))dt to be used in the x-axis of graphs (Figure 1D, 1F and Additional file 3: Figure S1E) and for making the regressor in the subsequent analysis of voxel-based statistics and of secondary neurodegeneration (see below). FA and rFA values in normal animals were used for the scatter plot, but not for fitting the kinetic models of neurodegeneration. For graph plotting, the cumulative risks of control (contralateral to glaucomatous side) and normal optic nerves were considered to be zero, which means that baseline risk, r 0 , was also zero.
For voxel-based analysis of FA images, we used tract-based spatial statistics (TBSS) , part of the FSL. First, FA images were brain-extracted using BET  and were aligned into a common space (a matrix of 137 × 167 × 100; a voxel size of 0.4 mm-cubic) using the nonlinear registration tool FNIRT, part of the FSL, which uses a b-spline representation of the registration warp field , followed by linear registration to the standard space of the anterior-posterior commissure line of the macaque brain . Next, the mean FA image was created and thinned to create a mean FA skeleton representing the centers of all tracts common to all data. The mean FA skeleton was created by thresholding the mean FA image at an FA value larger than 0.3. Each scan’s aligned FA data were then projected onto this skeleton and the resulting data fed into voxel-wise cross-subjects statistics for the linear regression analysis. The statistical analysis for estimation of the neurodegenerative model included analyzing the skeletonized FA images by testing significance of coefficients of post-operative time and the cumulative risk in the multi-regression statistical model. The analysis for the secondary neurodegenerative changes involved analyzing the skeletonized FA images by testing regressors of FA values in the affected optic nerve and of post-operative time. The predictor in the variable-risk model was made by using the optimized values for N 0 , r 0 and β determined from prior analysis of the quantitative FA in the primary degenerated area, the optic nerve (Table 1). For the statistical threshold, we applied threshold-free cluster-corrected P-values less than 0.05. Voxel-wise statistics were performed by permutation-based nonparametric inference using the program Randomise, part of the FSL. For post hoc analysis of the model of secondary neurodegeneration, a cubic ROI consisting of 27 voxels was placed with the center of the ROI located at the local maximum in the optic radiation in the contrast for the regressor of optic nerve FA. The model equation, Eq. 6, was fitted to the FA values in the areas of primary degeneration (glaucomatous optic nerve) and secondary degeneration (optic radiation), and the optimum values for k t and N B0 were determined. The model was also built based on the risk, by substituting N A (t) in Eq. 6 by Eq. 2, then was fitted to the FA data of the optic radiation to determine optimum values for N 0 , r 0 , and N B0 . An F-test was performed for model selection between constant- and variable-risk models, taking into account the additional parameter (β) in the latter.
The animals were analyzed by conventional histological evaluations, including toluidine blue and cresyl violet staining of retina, optic nerve and lateral geniculate nucleus. In addition, to correlate the optic nerve DTI changes with the degeneration of axons, optic nerve immunofluorescence staining was performed using SMI-31 antibody, a maker of neurofilament (NF), a constituent of neuronal axons. Animals were euthanized 33–168 days after glaucoma induction under deep general anesthesia (isoflurane in N2O and O2 by inhalation) and were perfused via the common carotid artery with 1 L of 0.9% saline containing 10 U/ml heparin at room temperature, followed by 1 L of 4% paraformaldehyde in 0.01 M phosphate-buffered saline (PBS; pH 7.4). Eyeballs with optic nerves were enucleated at the time of brain removal and 4% paraformaldehyde in PBS solution was injected into the vitreous body and postfixed by immersion in 4% paraformaldehyde in PBS for at least 1 week at 4°C. After immersion, a 3-mm long segment of optic nerve was cut out at 6 to 9 mm from the eyeball–optic nerve junction. The eyeballs and optic nerve segments were soaked in 10, 15, and 30% (w/v) sucrose in 0.1 M phosphate buffer, pH 7.4, at 4°C, for at least 24 h each, and then frozen in embedding compound (Tissue-Tek; Sakura Finetechnical Co.Ltd., Tokyo, Japan).
Coronal sections of the optic nerve segment were cut at 10-μm thickness and every 20th section was mounted onto the same slide glass until we obtained 20 slide glasses each containing four sections. Because cutting of the sections started from the end closer to the eyeball, the sections evaluated for counting should be 6–7 mm distant from the eyeball–nerve junction. Optic nerve sections were washed with 0.01 M PBS, preincubated with 10% normal goat serum in 0.01 M PBS for 30 min, and then incubated overnight at 4°C with mouse anti-SMI31 monoclonal antibody (NE1022; 1: 1000, Calbiochem, San Diego, CA, USA). They were washed with 0.01 M PBS and then incubated for 3 h at room temperature with Alexa Fluor 546 F(ab’)2 fragment of goat anti-mouse IgG (H + L) (1:1000 dilution) (Molecular Probes, Eugene, OR). For counting the number of SMI-31-positive axons, we randomly chose a single slide glass from the twenty, and identified five regions (nasal, temporal, superior, inferior, and central) in each of the four optic nerve sections. In the center of each region, we placed a sample volume of interest with a size of 218 × 164 × 10 μm (in x, y, and section thickness) under a microscope fitted with a 40× objective × 3 digital zoom; thus, a total of 0.18 mm2 per section (areal ratio of 3.6 ± 0.4 % as mean ± S.E.) was assessed for counting the number of SMI-31-positive (i.e. NF-positive) axons. Counting and photography were performed by a single observer who was blinded to the animals’ data and DTI results (Y. I). Thereby, for each optic nerve, we obtained 20 measurements of axon densities (five regions × four sections) and calculated the mean value that represents the NF-positive axon density of the corresponding optic nerve.
The obtained NF-positive axon density was used as an independent variable in the linear regression analysis of optic nerve FA. The ratio of NF density (between the glaucomatous and contralateral optic nerve) was used for fitting the risk-based kinetic models of neurodegeneration and determining the optimum values for r 0, and β. Model selection was also analyzed using an F-test between constant- and variable-risk models. Although the absolute NF density in monkey #1 (mild stage glaucomatous animal) was almost comparable with those in normal nerves (Additional file 6: Table S5), the ratio (to that in the contralateral non-affected optic nerve) was slightly decreased (0.86), similar to the axonal density (0.85) as evaluated by toluidine blue staining. In fact, toluidine blue-stained sections of this animal’s affected optic nerve revealed neurodegeneration and formation of myelin ovoids in the peripheral part of the optic nerve.
Model for secondary (transsynaptic) neurodegeneration
where k t is a constant coefficient. The value of k t should depend both on the strength of the trans-synaptic neurodegeneration and the connectivity between regions A and B, which may vary from region to region, or connectivity strength to the primary lesion in the brain. The differential equation (Eq. 7) indicates that the number of neurons in the area B (D B (t)) should decrease as a linear function of the primary degenerated area A (D A (t)) (the intercept of this linear function was substituted by N B0 ) as illustrated in Additional file 3: Figure S1A. Thus, if the neurodegeneration kinetics of the latter region is an exponential function of the cumulative risk (Eq. 2), the kinetics in region B should also be a decremental exponential function across the cumulative risk, but approach the plateau, N B0 , at infinite time (Additional file 3: Figure S1B).
General risk-based kinetic model of neurodegeneration
The general form of neurodegeneration kinetics, as formulated in Eq. 2 and 7, is virtually equivalent to those used in the survival analysis of the Cox hazard proportional regression model with time-dependent covariates [61, 62]. Thus, the treatment effect could be estimated by testing the significance of a covariate for the treatment or an interaction term between the pathogenic biomarker and the treatment.
We thank Akihide Yamamoto and Hiroshi Sato for their technical help, Kyoko Shioya for her generous support for handling animals, and Keizo Sugaya for his careful review of the manuscript. This study was supported by a Grant for Molecular Imaging Science Program from Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan; a KAKENHI Grant from MEXT [21613011 to T.H.], a Health Labour Sciences Research Grant from The Ministry of Health Labour and Welfare, Japan [H17-kokoro-025 to T.H.], and a Grant from the National Institute of Biomedical Innovation [10–05 to T.H.].
- DeKosky ST, Marek K: Looking backward to move forward: early detection of neurodegenerative disorders. Science. 2003, 302: 830-834. 10.1126/science.1090349.View ArticlePubMed
- Jack CR, Knopman DS, Jagust WJ, Shaw LM, Aisen PS, Weiner MW, Petersen RC, Trojanowski JQ: Hypothetical model of dynamic biomarkers of the Alzheimer’s pathological cascade. Lancet Neurol. 2010, 9: 119-128. 10.1016/S1474-4422(09)70299-6.PubMed CentralView ArticlePubMed
- Komarova NL, Thalhauser CJ: High Degree of Heterogeneity in Alzheimer’s Disease Progression Patterns. PLoS Comput Biol. 2011, 7: e1002251-10.1371/journal.pcbi.1002251.PubMed CentralView ArticlePubMed
- Fjell AM, Walhovd KB, Fennema-Notestine C, McEvoy LK, Hagler DJ, Holland D, Brewer JB, Dale AM: CSF biomarkers in prediction of cerebral and clinical change in mild cognitive impairment and Alzheimer’s disease. J Neurosci. 2010, 30: 2088-2101. 10.1523/JNEUROSCI.3785-09.2010.PubMed CentralView ArticlePubMed
- Perneczky R, Tsolakidou A, Arnold A, Diehl-Schmid J, Grimmer T, Förstl H, Kurz A, Alexopoulos P: CSF soluble amyloid precursor proteins in the diagnosis of incipient Alzheimer disease. Neurology. 2011, 77: 35-38. 10.1212/WNL.0b013e318221ad47.View ArticlePubMed
- Dickerson BC, Wolk DA: MRI cortical thickness biomarker predicts AD-like CSF and cognitive decline in normal adults. Neurology. 2012, 78: 84-90. 10.1212/WNL.0b013e31823efc6c.PubMed CentralView ArticlePubMed
- Shaw LM, Korecka M, Clark CM, Lee VM-Y, Trojanowski JQ: Biomarkers of neurodegeneration for diagnosis and monitoring therapeutics. Nat Rev Drug Discov. 2007, 6: 295-303. 10.1038/nrd2176.View ArticlePubMed
- Sherer TB: Biomarkers for parkinson’s disease. Sci Transl Med. 2011, 3: 79ps14-10.1126/scitranslmed.3002488.PubMed
- Drachman DA: Aging of the brain, entropy, and Alzheimer disease. Neurology. 2006, 67: 1340-1352. 10.1212/01.wnl.0000240127.89601.83.View ArticlePubMed
- Dunnett SB, Bjorklund A: Prospects for new restorative and neuroprotective treatments in Parkinson’s disease. Nature. 1999, 399: A32-A39.View ArticlePubMed
- Clarke G, Collins RA, Leavitt BR, Andrews DF, Hayden MR, Lumsden CJ, McInnes RR: A one-hit model of cell death in inherited neuronal degenerations. Nature. 2000, 406: 195-199. 10.1038/35018098.View ArticlePubMed
- Dubinsky JM, Kristal BS, Elizondo-Fournier M: On the probabilistic nature of excitotoxic neuronal death in hippocampal neurons. Neuropharmacology. 1995, 34: 701-711. 10.1016/0028-3908(95)00041-4.View ArticlePubMed
- Sugaya K, Matsubara S: Nucleation of protein aggregation kinetics as a basis for genotype-phenotype correlations in polyglutamine diseases. Mol Neurodegener. 2009, 4: 29-10.1186/1750-1326-4-29.PubMed CentralView ArticlePubMed
- Clark LG: The laminar organization and cell content of the lateral geniculate body in the monkey. J Anat Lond. 1941, 75: 419-433.
- Yoles E, Schwartz M: Degeneration of Spared Axons Following Partial White Matter Lesion: Implications for Optic Nerve Neuropathies. Exp Neurol. 1998, 153: 1-7. 10.1006/exnr.1998.6811.View ArticlePubMed
- Yucel Y, Gupta N: Glaucoma of the brain: a disease model for the study of transsynaptic neural degeneration. Prog Brain Res. 2008, 173: 465-478.View ArticlePubMed
- Kwon YH, Fingert JH, Kuehn MH, Alward WL: Primary open-angle glaucoma. N Engl J Med. 2009, 360: 1113-1124. 10.1056/NEJMra0804630.PubMed CentralView ArticlePubMed
- Leskea MC, Heijl A, Hyman L, Bengtsson B, Komaroff E: Factors for progression and glaucoma treatment: the Early Manifest Glaucoma Trial. Curr Opin Ophthalmol. 2004, 15: 102-106. 10.1097/00055735-200404000-00008.View Article
- Sasaoka M, Nakamura K, Shimazawa M, Ito Y, Araie M, Hara H: Changes in visual fields and lateral geniculate nucleus in monkey laser-induced high intraocular pressure model. Exp Eye Res. 2008, 86: 770-782. 10.1016/j.exer.2008.02.004.View ArticlePubMed
- Shimazawa M, Tomita G, Taniguchi T, Sasaoka M, Hara H, Kitazawa Y, Araie M: Morphometric evaluation of changes with time in optic disc structure and thickness of retinal nerve fibre layer in chronic ocular hypertensive monkeys. Exp Eye Res. 2006, 82: 427-440. 10.1016/j.exer.2005.08.001.View ArticlePubMed
- Gupta N, Yücel YH: Glaucoma as a neurodegenerative disease. Curr Opin Ophthalmol. 2007, 18: 110-114. 10.1097/ICU.0b013e3280895aea.View ArticlePubMed
- Danesh-Meyer HV: Neuroprotection in glaucoma: recent and future directions. Curr Opin Ophthalmol. 2011, 22: 78-86. 10.1097/ICU.0b013e32834372ec.View ArticlePubMed
- Douaud G, Jbabdi S, Behrens TEJ, Menke RA, Gass A, Monsch AU, Rao A, Whitcher B, Kindlmann G, Matthews PM, Smith S: DTI measures in crossing-fibre areas: increased diffusion anisotropy reveals early white matter alteration in MCI and mild Alzheimer’s disease. Neuroimage. 2011, 55: 880-890. 10.1016/j.neuroimage.2010.12.008.View ArticlePubMed
- Rose SE, Chen F, Chalk JB, Zelaya FO, Strugnell WE, Benson M, Semple J, Doddrell DM: Loss of connectivity in Alzheimer’s disease: an evaluation of white matter tract integrity with colour coded MR diffusion tensor imaging. J Neurol Neurosurg Psychiatr. 2000, 69: 528-530. 10.1136/jnnp.69.4.528.PubMed CentralView ArticlePubMed
- Yoshikawa K, Nakata Y, Yamada K, Nakagawa M: Early pathological changes in the parkinsonian brain demonstrated by diffusion tensor MRI. J Neurol Neurosurg Psychiatr. 2004, 75: 481-484. 10.1136/jnnp.2003.021873.PubMed CentralView ArticlePubMed
- Péran P, Cherubini A, Assogna F, Piras F, Quattrocchi C, Peppe A, Celsis P, Rascol O, Démonet J-F, Stefani A, Pierantozzi M, Pontieri FE, Caltagirone C, Spalletta G, Sabatini U: Magnetic resonance imaging markers of Parkinson’s disease nigrostriatal signature. Brain. 2010, 133: 3423-3433. 10.1093/brain/awq212.View ArticlePubMed
- Ellis CM, Simmons A, Jones DK, Bland J, Dawson JM, Horsfield MA, Williams SC, Leigh PN: Diffusion tensor MRI assesses corticospinal tract damage in ALS. Neurology. 1999, 53: 1051-1058. 10.1212/WNL.53.5.1051.View ArticlePubMed
- Sach M, Winkler G, Glauche V, Liepert J, Heimbach B, Koch MA, Büchel C, Weiller C: Diffusion tensor MRI of early upper motor neuron involvement in amyotrophic lateral sclerosis. Brain. 2004, 127: 340-350. 10.1093/brain/awh041.View ArticlePubMed
- Garaci FG, Bolacchi F, Cerulli A, Melis M, Spanò A, Cedrone C, Floris R, Simonetti G, Nucci C: Optic nerve and optic radiation neurodegeneration in patients with glaucoma: in vivo analysis with 3-T diffusion-tensor MR imaging. Radiology. 2009, 252: 496-501. 10.1148/radiol.2522081240.View ArticlePubMed
- Basser PJ, Mattiello J, LeBihan D: MR diffusion tensor spectroscopy and imaging. Biophys J. 1994, 66: 259-267. 10.1016/S0006-3495(94)80775-1.PubMed CentralView ArticlePubMed
- Beaulieu C: Diffusion MRI. From Quantitative Measurement to In Vivo Neuroanatomy. Edited by: Johansen-Berg H, Behrens T. 2009, Elsevier, London
- Quigley HA, Hohman RM: Laser energy levels for trabecular meshwork damage in the primate eye. Invest Ophthalmol Vis Sci. 1983, 24: 1305-1307.PubMed
- Clarke G, Lumsden CJ: Scale-free neurodegeneration: cellular heterogeneity and the stretched exponential kinetics of cell death. J Theor Biol. 2005, 233: 515-525. 10.1016/j.jtbi.2004.10.028.View ArticlePubMed
- Huber DL: Statistical model for stretched exponential relaxation in macroscopic systems. Phys Rev B. 1985, 31: 6070-6071. 10.1103/PhysRevB.31.6070.View Article
- Kondo Y, Takada M, Tokuno H, Mizuno N: Single retinal ganglion cells projecting bilaterally to the lateral geniculate nuclei or superior colliculi by way of axon collaterals in the cat. J Comp Neurol. 1994, 346: 119-126. 10.1002/cne.903460108.View ArticlePubMed
- Poltorak M, Freed WJ: Immunoreactive phosphorylated epitopes on neurofilaments in neuronal perikarya may be obscured by tissue preprocessing. Brain Res. 1989, 480: 349-354. 10.1016/0006-8993(89)90206-0.View ArticlePubMed
- Gusella JF, MacDonald ME: Molecular genetics: unmasking polyglutamine triggers in neurodegenerative disease. Nat Rev Neurosci. 2000, 1: 109-115.View ArticlePubMed
- Murphy RM: Kinetics of amyloid formation and membrane interaction with amyloidogenic proteins. Biochim Biophys Acta. 2007, 1768: 1923-1934. 10.1016/j.bbamem.2006.12.014.View ArticlePubMed
- Ashe KH, Zahs KR: Probing the biology of Alzheimer’s disease in mice. Neuron. 2010, 66: 631-645. 10.1016/j.neuron.2010.04.031.PubMed CentralView ArticlePubMed
- Douglas PM, Dillin A: Protein homeostasis and aging in neurodegeneration. J Cell Biol. 2010, 190: 719-729. 10.1083/jcb.201005144.PubMed CentralView ArticlePubMed
- Craft S: The role of metabolic disorders in Alzheimer disease and vascular dementia: two roads converged. Arch Neurol. 2009, 66: 300-305. 10.1001/archneurol.2009.27.PubMed CentralView ArticlePubMed
- Prinz M, Priller J, Sisodia SS, Ransohoff RM: Heterogeneity of CNS myeloid cells and their roles in neurodegeneration. Nat Neurosci. 2011, 14: 1227-1235.View ArticlePubMed
- Chételat G, Villemagne VL, Bourgeat P, Pike KE, Jones G, Ames D, Ellis KA, Szoeke C, Martins RN, O’Keefe GJ, Salvado O, Masters CL, Rowe CC: Relationship between atrophy and beta-amyloid deposition in Alzheimer disease. Ann Neurol. 2010, 67: 317-324.PubMed
- El-Rafei A, Engelhorn T, Wärntges S, Dörfler A, Hornegger J, Michelson G: A framework for voxel-based morphometric analysis of the optic radiation using diffusion tensor imaging in glaucoma. Magn Reson Imaging. 2011, 29: 1076-1087. 10.1016/j.mri.2011.02.034.View ArticlePubMed
- Assaf Y, Cohen Y: Non-mono-exponential attenuation of water and N-acetyl aspartate signals due to diffusion in brain tissue. J Magn Reson. 1998, 131: 69-85. 10.1006/jmre.1997.1313.View ArticlePubMed
- Jack CR, Wiste HJ, Vemuri P, Weigand SD, Senjem ML, Zeng G, Bernstein MA, Gunter JL, Pankratz VS, Aisen PS, Weiner MW, Petersen RC, Shaw LM, Trojanowski JQ, Knopman DS: Brain beta-amyloid measures and magnetic resonance imaging atrophy both predict time-to-progression from mild cognitive impairment to Alzheimer’s disease. Brain. 2010, 133: 3336-3348. 10.1093/brain/awq277.PubMed CentralView ArticlePubMed
- Brundin P, Melki R, Kopito R: Prion-like transmission of protein aggregates in neurodegenerative diseases. Nat Rev Mol Cell Biol. 2010, 11: 301-307. 10.1038/nrm2873.PubMed CentralView ArticlePubMed
- Pate BD, Kawamata T, Yamada T, McGeer EG, Hewitt KA, Snow BJ, Ruth TJ, Calne DB: Correlation of striatal fluorodopa uptake in the MPTP monkey with dopaminergic indices. Ann Neurol. 1993, 34: 331-338. 10.1002/ana.410340306.View ArticlePubMed
- Gerhard A, Pavese N, Hotton G, Turkheimer F, Es M, Hammers A, Eggert K, Oertel W, Banati RB, Brooks DJ: In vivo imaging of microglial activation with [11C](R)-PK11195 PET in idiopathic Parkinson’s disease. Neurobiol Dis. 2006, 21: 404-412. 10.1016/j.nbd.2005.08.002.View ArticlePubMed
- Shimazawa M, Ito Y, Inokuchi Y, Yamanaka H, Nakanishi T, Hayashi T, Ji B, Higuchi M, Suhara T, Imamura K, Araie M, Watanabe Y, Onoe H, Hara H: An alteration in the lateral geniculate nucleus of experimental glaucoma monkeys: In vivo positron emission tomography imaging of glial activation. PLoS One. 2012, 7: e30526-10.1371/journal.pone.0030526.PubMed CentralView ArticlePubMed
- Obeso JA: Modeling clinical features of neurodegeneration. Nat Med. 2010, 16: 1372-View ArticlePubMed
- Weon BM, Je JH: Theoretical estimation of maximum human lifespan. Biogerontology. 2009, 10: 65-71. 10.1007/s10522-008-9156-4.View ArticlePubMed
- Weibull W: A statistical distribution function of wide applicability. J Appl Mech. 1951, 18: 293-297.
- Ito Y, Shimazawa M, Chen Y-N, Tsuruma K, Yamashima T, Araie M, Hara H: Morphological changes in the visual pathway induced by experimental glaucoma in Japanese monkeys. Exp Eye Res. 2009, 89: 246-255. 10.1016/j.exer.2009.03.013.View ArticlePubMed
- Reese TG, Heid O, Weisskoff RM, Wedeen VJ: Reduction of eddy-current-induced distortion in diffusion MRI using a twice-refocused spin echo. Magn Reson Med. 2003, 49: 177-182. 10.1002/mrm.10308.View ArticlePubMed
- Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TE, Johansen-Berg H, Bannister PR, De Luca M, Drobnjak I, Flitney DE, Niazy RK, Saunders J, Vickers J, Zhang Y, De Stefano N, Brady JM, Matthews PM: Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage. 2004, 23 (Suppl 1): S208-S219.View ArticlePubMed
- Smith SM, Jenkinson M, Johansen-Berg H, Rueckert D, Nichols TE, Mackay CE, Watkins KE, Ciccarelli O, Cader MZ, Matthews PM, Behrens TE: Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data. Neuroimage. 2006, 31: 1487-1505. 10.1016/j.neuroimage.2006.02.024.View ArticlePubMed
- Smith SM: Fast robust automated brain extraction. Hum Brain Mapp. 2002, 17: 143-155. 10.1002/hbm.10062.View ArticlePubMed
- Rueckert D, Sonoda LI, Hayes C, Hill DL, Leach MO, Hawkes DJ: Nonrigid registration using free-form deformations: application to breast MR images. IEEE Trans Med Imaging. 1999, 18: 712-721. 10.1109/42.796284.View ArticlePubMed
- Bowden DM, Martin MF: Primate Brain Maps: Structure of the Macaque Brain: A Laboratory Guide with Original Brain Sections, Printed Atlas and Electronic Templates for Data and Schematics. 2000, Elsevier Science, Amsterdam, 1
- Cox DR: Regression Models and Life-Tables. J Statist Assoc Series B (Methodological). 1972, 34: 187-220.
- Fisher LD, Lin DY: Time-dependent covariates in the Cox proportional-hazards regression model. Annu Rev Public Health. 1999, 20: 145-157. 10.1146/annurev.publhealth.20.1.145.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.