|
|
||||||||
1From the Department of Ophthalmology, Brown University Medical School, Providence, Rhode Island; and the 2Departments of Ophthalmology and 3Neuroscience, Mount Sinai School of Medicine, New York, New York.
| Abstract |
|---|
|
|
|---|
METHODS. Fifteen-month-old DBA2/NNia mice were labeled retrogradely with fluorogold. Retinas were flat mounted and imaged in their entirety using an epifluorescence microscope with a motorized stage. Digital maps of the retinal wholemounts were constructed to automatically count and establish spatial coordinates for RGCs over the entire retina. RGC size and NND were determined from these maps.
RESULTS. RGC counts in the group of 15-month-old DBA/2NNia animals ranged from 22,330 to 92,157 cells per retina. Mean RGC cell size per retina ranged from 22.35 to 35.64 µm2 and correlated linearly with total RGC counts. NND distribution histograms were compared for retinas with variable degrees of RGC loss. The distribution of NNDs in each retina was skewed toward larger distance values in more affected retinas. In partially damaged retinas, areas with severe pathology coincided with areas of maximal loss of large RGCs, and areas of preserved RGCs correlated with larger cell sizes.
CONCLUSIONS. Damaged retinas have a smaller mean cell size, indicating preferential loss of larger RGCs or size reduction of surviving cells. NND analysis of the RGC population in a retina is a useful measure of glaucomatous RGC loss. The skewed NND distribution of surviving RGCs and the finding that RGC loss correlates with a shift/amplitude change in the mode of the histogram and its tail suggests two different patterns of RGC loss possibly attributable to different pathologic processes in glaucomatous DBA/2 mice.
More recently, murine models of glaucomatous optic neuropathy have also become of interest to some investigators.11 12 13 14 15 16 Highly inbred mouse strains have a more restricted biologic variability that, together with spontaneous or gene-induced development of high IOP,11 12 17 18 19 makes them desirable glaucoma models for investigation. Other potentially important examples for glaucoma research are mice with gene knockouts of components of the apoptotic cell death program of neurons,20 21 which is thought to be the primary mechanism for retinal ganglion cell (RGC) loss in glaucoma. A critical requirement in most studies of mouse glaucoma models is the quantification of RGC somata or their axons.
Automated counting methods to quantify RGC loss in the glaucomatous retina have certain advantages over manual counting based on sampling.22 The assessment of RGCs in virtually 100% of the retinal area provides increased counting precision and allows the investigation of spatial and cell size patterns of surviving RGCs in a retina with glaucoma. In this study, we examined the detailed spatial arrangement and the cell size of RGCs in murine retinas with different degrees of glaucomatous RGC loss. We used an automated counting method to quantify the glaucomatous focal RGC loss encountered in aged DBA/2NNia mice.16 The spontaneous development of secondary glaucoma in this mouse strain is partially age and sex dependent and results from gene mutations causing iris atrophy and pigment dispersion,12 together with an immunologic component,23 that is not yet fully understood. Therefore, individual eyes in a nonselected, mixed-sex population of DBA/2 mice of the same age can exhibit a range of retinal pathology. We used such a group of mice at 15 months of age to study the patterns of RGC survival in retinas with varying degrees of glaucomatous neuropathy. We created data libraries of the studied retinas and digitally constructed spatial RGC maps to extend the analysis from cell counts to morphologic and topographic associations of surviving RGCs. To further characterize the previously reported spatial non-uniformity of RGC loss,16 we investigated the use of nearest neighbor distance (NND) for every RGC in a retina as a measure of focal RGC density. NND is defined as the distance between the centroid of each RGC and the centroid of its closest neighbor. The centroid is defined as the moment center for each object and can be calculated for any object, symmetric or asymmetric. It corresponds to its geometric center for a symmetric object. NND and the spatial coordinates of the centroid for each object in a binary image are directly provided by a plug-in of Photoshop 5.5 (Adobe, San Jose, CA).24 Consequently, for each cell, the NND approximates the radius of a circular retinal area occupied by that cell alone. NND is thus inversely related to RGC density and is independent of the constraints inherent in density analysis of RGCs calculated from sampling frames. Therefore, the NND approach has potential advantages for quantifying and characterizing the overall RGC loss in a retina compared with gross RGC counts estimated from sampling or frame cell density analysis.
| Methods |
|---|
|
|
|---|
Retrograde Labeling
Mice were anesthetized by intraperitoneal administration of xylazine (10.8 mg/kg), acepromazine (1.2 mg/kg), and ketamine (54.0 mg/kg). Their skulls were exposed, and holes 2 mm in diameter were drilled bilaterally at points 2.92 mm posterior and 0.5 mm lateral to the bregma to expose the occipital cortex.26 27 Under constant direct observation, the occipital cortex overlying the superior colliculus (SC) was gently aspirated and the SC was exposed. A piece of gel foam (Pharmacia & Upjohn, Kalamazoo, MI) soaked in 5% aqueous fluorogold (Fluorochrome, Denver, CO) was applied to each SC. The gel foam was covered with antibiotic ointment, and the overlying skin was sutured.
Flat Mounting/Digital Imaging
After 7 days, to allow for the retrograde labeling of retinal ganglion cells, the animals were humanely killed by transcardial perfusion with 4% paraformaldehyde in 0.1 M phosphate-buffered saline under general anesthesia. Eyes were immediately enucleated; retinas were dissected from the globes and prepared as flattened wholemounts after relaxing incisions were performed, preserving orientation. Retinal wholemounts were imaged on an epifluorescence microscope (Axioskop 2; Zeiss, Thornwood, NY) equipped with a digital camera (Monospot; Diagnostic Instrumentation Inc., Sterling Heights, MI) and a motorized stage (Biopoint; Ludl Electronic Products, Ltd., Hawthorne, NY), allowing the sequential acquisition of adjacent nonoverlapping frames that covered the total area of each retina.22 Imaging was performed with a 10x objective, NA 0.50 mm (Zeiss), in 30 frames per retina arranged in a reproducible 5 x 6 grid. Images from 8 of 10 flat-mounted retinas of adequate quality were further analyzed for construction of NND maps.
Digital Processing
Commercially available software (Photoshop 5.5 [Adobe]; Image Tool 2.0 [University of Texas Health Science Center at San Antonio, San Antonio, TX]) was used to convert the fluorescent digital retinal images to corresponding black-and-white images used for the subsequent analysis, as previously described.16 22
To validate the automated RGC counting method for this imaging setup, automated counts were compared with those obtained from images counted manually by two investigators in a masked fashion. More than 146,000 manually counted RGCs from another set of flat-mounted retinal images were used for this purpose. The size threshold criteria for black objects representing RGCs versus background were determined by calculating the object size in pixels of the smallest among a large number of RGCs positively identified on the original images as RGCs. This size threshold coincided with the inflection point of the object size distribution in the set of manually counted frames, as previously reported.16 22 This inflection point separated the bimodal object size distribution, with the first mode representing the noise and the latter representing RGCs. In addition, the determined size threshold provided the highest correlation between manual and automated counts (R2 = 0.9631) (Fig. 1) . By using this algorithm, total RGC counts were determined for the eight retinas studied in detail.
|
Finally, from the same data libraries, we calculated a mean object size, reflected by its area, for all objects (RGCs) and correlated this value with the total RGC counts. Pixel size for an object was converted (Photoshop 5.5; Adobe) to a circular area of equivalent diameter providing, in this case, an estimate of RGC size (function: measure all in Photoshop 5.5; Adobe).
Statistical Analysis
Commercial programs (NCSS 2000 and PASS 6.0 [NCSS]; Excel 2000 [Microsoft Corp., Redmond, WA]) were used to perform all statistical analyses.
| Results |
|---|
|
|
|---|
Nearest Neighbor Distance
Mean NND values correlated well with total cell counts for all retinas (R2 = 0.9621; 95% confidence interval, 0.9710.987; P < 0.001) (Fig. 2) .
|
1%-2%) with high NND values (Fig. 3A) . Therefore, all objects whose NNDs exceeded 20 arbitrary units or 21.09 µm were included in the last bin for illustration purposes (Fig. 3A , histogram; Fig. 3B , cumulative plot). This percentage of cells with NND greater than 20 arbitrary units differed significantly between preserved retinas that maintained approximately 75% of their RGCs16 (total RGCs > 64,000; n = 3) and severely affected retinas (total RGCs < 64,000; n = 5) (P < 0.001; Student t test) and was inversely related to the total count of surviving RGCs in senescent DBA/2NNia murine retinas, as shown in Figure 4 (R2 = 0.8958).
|
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
As illustrated here, automated counting methods may allow for easier extraction of detailed information about RGCs than sampling and manual counting routines. Information that can be obtained from digital images includes not only total RGC counts per retina but also estimates of local density, cell size, and shape information that allows mapping of RGC loss and survival anywhere in the retina.
NND analysis provides an alternative measure to quantify RGC loss and has the capability to differentiate more uniformly distributed cell loss from the patchy, non-uniform type of loss typical of glaucoma. NND has conceptual homologies to the inverse of the focal cell density calculated for each individual cell and thus does not depend on averaging effects when cells are counted in a representative area (frame) of the retina. As with density, NND defined as centroid-to-centroid distance is affected to some extent by cell size. Despite this limitation, we reasoned that NND might be better able to differentiate between focal and more uniformly distributed RGC loss within each retina, thus providing a better way to model the pattern of RGC loss and allowing comparisons between retinas. Diffuse loss would manifest as a shift of the NND distribution to the right, whereas focal loss would appear as an increase in the proportion of cells with very large NNDs. The NND distribution analysis indicates that both types of loss appear to take place simultaneously in the DBA/2 mouse model of glaucomatous optic neuropathy, suggesting that there may be two pathologic processes leading to RGC depletion in this model of glaucoma. Our findings validate the use of NND as a useful tool for detecting focal and diffuse RGC loss. Examining the effect of the neuroprotective agents/inhibitors of the apoptotic pathway on the distribution characteristics of the NND, as shown in Figure 3 , could provide a means to differentiate effects on diffuse or focal RGC loss. Finally, NND analysis can be used to compare different glaucoma models and to provide indications of the type of RGC pathology these models cause in relation to the duration and severity of the insult. An additional important advantage is that NND analysis for a group of retinas provides an overall measure of disease that is less dependent on biologic variability and potential inaccuracies than the single number of a total RGC count derived from manual or automated counting and may prove itself to be more sensitive in detecting early glaucomatous loss. Therefore, the increased statistical power for comparisons of NND distributions may significantly increase our ability to screen for effective treatments that preserve RGCs.
Cell size analysis is an additional capability offered by the construction of digital libraries of all the retrogradely labeled objects (RGCs) within each retina. The large number of RGCs counted by automated methods can overcome shortcomings in accuracy of the absolute cell size compared with more precise estimates derived from a few RGCs in limited areas of the retina captured under high magnification. In the present study, occasional underestimation or overestimation of RGC size because of depth-of-field variations in the captured images or merged cells appears to be random and does not make cell size comparison among retinas any less valid. This conclusion is based on the finding that correlations between cell size and cell loss in the data presented were not affected when the databases were filtered from grossly asymmetric objects representing artifacts or merging of more than one cell (data not shown). The filtering process resulted in the removal of less than 1.5% of the objects and only minimally affected the absolute values of cell sizes and their distribution.
Mean cell size, estimated from the area of the imaged RGCs, in affected retinas of senescent DBA/2 animals appears to be 15% to 20% smaller than in undamaged retinas and results from relative depletion of larger RGCs compared with smaller ones or from cell shrinkage. Larger cells also appear to be located on islands of preserved RGCs within affected retinas. This finding is in agreement with reports that ocular hypertension is associated with an overall size reduction in RGCs in monkey, cat, and rat glaucoma models (Cordeiro MF, et al. IOVS 2004;45:ARVO E-Abstract 1114; Guo L, et al. IOVS 2004;45:ARVO E-Abstract 2153),28 29 30 31 32 and has recently been reported to also occur in a rat glaucoma model.33 The degree of overall cell size reduction in the DBA/2 retina is within the range reported for Macaca fascicularis monkeys after 6 to 14 weeks of ocular hypertension.34 In this study, we present new evidence of spatial correlation between larger cell size and areas of RGC preservation. This result implies the converse, namely that in those areas of a glaucomatous retina in which extensive RGC loss occurs, the largest cells die, or, if they survive, they shrink. Our data, however, could not indicate whether smaller RGCs are less susceptible to injury or whether RGCs in the glaucomatous retina tend to shrink before cell death. Nevertheless our findings show that cell size distribution analysis can provide another measure of the extent of disease in the glaucomatous retina.
In summary, our data confirm the usefulness of automated RGC counting as a method for evaluating glaucomatous retinal damage. Measures of RGC distribution, such as NND, and RGC size can improve our understanding of the pathophysiology of the disease and increase the power of experimental studies. RGC distribution characteristics, such as NND, and RGC cell size complement RGC total cell number as a method of quantifying RGC cell loss in a glaucoma model and offer a more complete description of the disease in the RGC layer. Application of these additional measures of RGC loss may help overcome some of the problems caused by the variability encountered in many glaucoma models.
| Acknowledgements |
|---|
| Footnotes |
|---|
Submitted for publication July 22, 2005; revised November 19, 2005; accepted March 13, 2006.
Disclosure: T. Filippopoulos, None; J. Danias, None; B. Chen, None; S.M. Podos, None; T.W. Mittag, 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: Theodoros Filippopoulos, Department of Ophthalmology, Brown University Medical School, APC 7th, 593 Eddy Street, Providence, RI 02903; theodoros_filippopoulos{at}brown.edu.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
Y. Zhou, O. Grinchuk, and S. I. Tomarev Transgenic Mice Expressing the Tyr437His Mutant of Human Myocilin Protein Develop Glaucoma Invest. Ophthalmol. Vis. Sci., May 1, 2008; 49(5): 1932 - 1939. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. P. Buckingham, D. M. Inman, W. Lambert, E. Oglesby, D. J. Calkins, M. R. Steele, M. L. Vetter, N. Marsh-Armstrong, and P. J. Horner Progressive Ganglion Cell Degeneration Precedes Neuronal Loss in a Mouse Model of Glaucoma J. Neurosci., March 12, 2008; 28(11): 2735 - 2744. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Scholz, T. Buder, S. Seeber, E. Adamek, C.-M. Becker, and E. Lutjen-Drecoll Dependency of Intraocular Pressure Elevation and Glaucomatous Changes in DBA/2J and DBA/2J-Rj Mice Invest. Ophthalmol. Vis. Sci., February 1, 2008; 49(2): 613 - 621. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Soto, E. Oglesby, B. P. Buckingham, J. L. Son, E. D. O. Roberson, M. R. Steele, D. M. Inman, M. L. Vetter, P. J. Horner, and N. Marsh-Armstrong Retinal Ganglion Cells Downregulate Gene Expression and Lose Their Axons within the Optic Nerve Head in a Mouse Glaucoma Model J. Neurosci., January 9, 2008; 28(2): 548 - 561. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. R. Howell, R. T. Libby, T. C. Jakobs, R. S. Smith, F. C. Phalan, J. W. Barter, J. M. Barbay, J. K. Marchant, N. Mahesh, V. Porciatti, et al. Axons of retinal ganglion cells are insulted in the optic nerve early in DBA/2J glaucoma J. Cell Biol., December 31, 2007; 179(7): 1523 - 1537. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |