Abstract
Optic neuropathies comprise a diverse group of disorders that ultimately lead to retinal ganglion cell (RGC) degeneration. Despite varying etiologies, these conditions share a conserved pathological progression: axonal damage in the optic nerve triggers progressive RGC degeneration. Understanding species-specific differences in neuronal resilience is critical for identifying key survival mechanisms and potential neuroprotective targets. In this study, we compare RGC densities and survival rates following optic nerve crush (ONC) in three vertebrate models—mice, zebrafish, and killifish—under standardized experimental conditions. Transcriptomic analysis confirmed that, similar to RBPMS in mice, Rbpms2 serves as a pan-RGC marker in zebrafish and killifish. Using these markers, we reveal significant species-specific differences in RGC density, with fish species exhibiting over a 5-fold higher density than mice at equivalent life stage. Killifish also show an age-dependent decline in RGC density. Furthermore, we identify distinct injury responses across species: mice undergo rapid degeneration, losing ∼80% of their RGCs by day 14 after ONC; zebrafish maintain full RGC retention for 2 weeks before experiencing a loss of ∼12%; and killifish display a biphasic response to ONC, with young adults retaining two-thirds of their RGCs by day 21, while older fish exhibit a more pronounced second wave of RGC loss, ultimately preserving just over half of their RGCs by 21 days after injury. These findings highlight fundamental differences in neuroprotective capacity among species, providing a comparative framework to uncover molecular mechanisms governing RGC survival and to identify therapeutic strategies for treating optic neuropathies and neurodegeneration across diverse pathologies.
1 Introduction
Neuronal loss is a hallmark of both acute and chronic neurodegenerative conditions, including those affecting the visual system. Optic neuropathies encompass diverse disorders, from glaucomatous neurodegeneration to ischemic and traumatic optic nerve injuries. Despite differing etiologies, these conditions converge on a common pathological sequence: initial damage of retinal ganglion cell (RGC) axons in the optic nerve in turn triggering progressive RGC death within the retina (; ). Being the retina’s sole output neurons, loss of RGCs disrupts visual signaling to the brain, resulting in irreversible blindness. A key challenge in treating optic neuropathies—and neurodegenerative diseases broadly—is preventing neuronal degeneration and enhancing resilience to injury. Thus, elucidating the molecular and cellular mechanisms underlying RGC vulnerability is essential for advancing neuroprotective strategies.
To this end, model organisms, and the differences among them, are essential for unraveling RGC biology and pathophysiology in retinal disorders. RGC density varies across species and age, affecting visual acuity and disease susceptibility (). Mice are widely used due to their genetic similarity to humans and well-established models of glaucomatous injury, exhibiting a comparable central nervous system (CNS) injury response (). Zebrafish provide insights into neurorepair due to their remarkable regenerative abilities, which are maintained throughout adulthood (; ). Killifish are by now a well-accepted gerontology model, with retinal aging mirroring key hallmarks of human aging, including increased oxidative stress, stem cell exhaustion, gliosis and inflammaging (; ). Furthermore, albeit teleost fish, killifish lose their regenerative capabilities and fail to recover from CNS injury at old age (; ), thereby recapitulating mammalian-like phenotypes after injury (). Therefore, they offer a valuable model for studying age-related changes in RGC vulnerability, as shown by increased stress response systems and reduced expression of bdnf (brain-derived neurotropic factor) in the aged retina (). Thus, species-specific differences in injury responses may uncover key factors influencing CNS neuronal survival and, in turn, identify novel neuroprotective targets that could form the foundation for new therapeutic strategies.
Although RGC survival after optic nerve crush (ONC) has been studied singularly in mice, zebrafish, and killifish, a comprehensive, cross-species analysis using whole-mount (WM) retinas for zebrafish and killifish has yet to be performed. Unlike the use of narrow field-of-view micrographs, entire retinal WMs rule out regional differences in injury response that might skew the general outcome. This is particularly relevant for killifish, which, due to their fast-growing nature, undergo significant retinal stretching during adulthood (), potentially leading to overestimation of RGC loss when using spatially pre-selected samples like cryosections (). Furthermore, a major limitation in fish models has been suboptimal RGC labeling methods. In zebrafish, RGC survival has largely been assessed using Tg(isl2b:reporter) lines (; ; ). Recent single cell RNAseq data however, disputes the validity of isl2b (ISL LIM homeobox 2b) as a pan-RGC marker in zebrafish, as its expression is neither uniform nor ubiquitous across RGC subtypes (). In killifish, RGCs have been labeled via retrograde biocytin tracing (; ), but this can lead to incomplete labeling of RGCs when the procedure is not carried out correctly. Furthermore, RGC axons are also marked by retrograde biocytin tracing, occluding the view of the ganglion cell layer in regions of high axonal density, such as in the central retina. Conversely, Rbpms2 (RNA-binding protein with multiple splicing 2) has emerged as a pan-RGC marker in both species based on RNA sequencing data (; ), but has yet to be used to assess RGC densities and/or survival rates following ONC in these teleost models.
In this study, we perform optic nerve crush (ONC) in three key model organisms—mice, zebrafish, and killifish—to directly compare RGC survival across species. By standardizing experimental conditions within a single laboratory, we minimize inter-laboratory variability, strengthening cross-species comparisons. Using pan-RGC markers, we reveal species-specific RGC densities and injury responses. The studied fish species exhibit significantly higher RGC densities compared to mice of an equivalent life stage, and RGC density declines with increasing age in killifish. Furthermore, murine RGCs exhibit low resilience to ONC, zebrafish preserve RGC numbers for 2 weeks before a delayed degenerative phase, and killifish display a biphasic cell loss pattern affected by age. As such, our findings provide a foundation for future studies exploring molecular mechanisms underlying differential RGC resilience across vertebrates.
2 Materials and methods
2.1 Transcriptomic data
To display the expression pattern of RGC markers across the retina, previously published scRNAseq datasets were used: (1) for mouse, (2) for zebrafish, and (3) for killifish. All three datasets were used “as is” and the original dimensionality reduction projections and clustering were used for plotting. For zebrafish, only the non-injured cells of the NMDA dataset of Hoang et al. were used.
2.2 Animal housing
Mice (Mus musculus) of the C57Bl/6N strain (Charles River Laboratories, France) were housed under standard laboratory conditions (12 h light/12 h dark cycle, 21°C, 50% humidity) (). Water was accessible at all times and they were fed ad libitum. All experiments were carried out using 10 weeks-old mice of both sexes.
Zebrafish (Danio rerio) of the AB wildtype strain were raised and maintained under standard laboratory conditions as described (), specifically at 28°C with a conductivity of 650 μS and pH of 7.5 on a 14 h light/10 h dark cycle. They were fed twice a day with a mixture of dry food and brine shrimp (Artemia Salina nauplii, Ocean Nutrition). All experiments were conducted on 21 weeks-old adult zebrafish of similar size, including both males and females.
African turquoise killifish (Nothobranchius furzeri) were raised and housed at standard laboratory conditions (), i.e., at 28°C with a conductivity of 600 μS and pH of 7.0 on a 12 h light/dark cycle. Fish were fed twice a day with a combination of brine shrimp and Chironomidae larvae. All experiments made use of 6- and 18 weeks-old killifish of the GRZ-AD inbred strain. Only female fish were used to minimize sex difference as male and female killifish differ substantially in size and age trajectory, potentially altering both cell densities and survival properties.
All animal experiments received approval from the KU Leuven Animal Ethics Committee and were conducted in strict compliance with the European Communities Council Directive of 20 October 2010 (2010/63/EU).
2.3 Optic nerve crush
In mice, a unilateral ONC was performed as previously described (; ). Briefly, animals were anesthetized through intraperitoneal injection of ketamine (75 mg/kg body weight, Nimatek, Eurovet) and medetomidine (1 mg/kg, Domitor, Pfizer). Next, the optic nerve was exposed after an incision of the conjunctiva and crushed 1 mm from the optic nerve head using a Dumont #7 cross-action forceps (Fine Science Tools) for 10 s. After ONC, anesthesia was reversed using atipamezole (1 mg/kg, Antisedan, Pfizer). Additionally, local analgesia (oxybuprocaïne 0.4%, Unicaïne) was applied to the eye before surgery, and antibiotic ointment (tobramycin 0.3%, Tobrex) was applied afterward.
Detailed protocols to unilaterally crush both zebrafish and killifish optic nerves have been previously reported (; ). Briefly, fish were anesthetized using 0.03% Tris-buffered tricaine (MS-222, Sigma-Aldrich). The connective tissue surrounding the eye was removed and the eyeball was lifted out of its orbit, exposing both the optic nerve and ophthalmic artery. Using a Dumont #7 cross-action forceps (Fine Science Tools), the optic nerve was crushed at a distance of 0.5 mm from the optic nerve head for 10 s, without damaging the ophthalmic artery. Both the distance from the optic nerve head and the duration of the crush were kept constant for young adult zebrafish, young adult and old killifish. After the procedure, the fish were returned to system water to recover.
2.4 Retrograde labeling of fish retinal ganglion cells using biocytin
Retrograde labeling of both teleost RGCs was performed as described in detail for both zebrafish (), and killifish (). Fish were anesthetized and the eye was lifted from its socket as described above. The optic nerve was completely cut at a distance of 500 μm from the optic nerve head. A gel foam, drenched with a saturated biocytin solution, was placed at the proximal end of the optic nerve. The fish were awakened in a recovery tank for 3 h to allow passive retrograde transport of the tracer. Next, fish were euthanized, and tissues were collected (see section “2.5 Tissue sampling and processing”).
2.5 Tissue sampling and processing
All mice were killed with an overdose of sodium pentobarbital (60 mg/kg, Dolethal, Vetoquinol), while fish were euthanized using an overdose of tricaine (0.1% Tris-buffered tricaine, MS-222). Animals were subsequently transcardially perfused using 0.9% NaCl (mice) or phosphate buffered saline (PBS, fish) and 4% paraformaldehyde (PFA) as previously described (Wu et al., 2021; ; ). Eyes were fixated in 4% PFA for 1 h after enucleation and subsequently washed three times with PBS and stored in storage buffer (0.4M NaN3 in PBS) until further use. Retinas were dissected and whole-mounted, fixated for an additional hour in 4% PFA, washed three times in PBS and stored in storage buffer until further use.
Next, murine and fish retinal WMs were immunohistochemically stained for RBPMS (rabbit anti-RBPMS, PhosphoSolutions, AB_249225, 1:250) or Rbpms2 (rabbit anti-Rbpms2, Abcam, Ab181098, 1:200), respectively. WMs were permeabilized by washing steps in 0.5% Triton X-100 in PBS (PBST 0.5%) and a 15 min freeze-thaw step at −80°C. Samples were blocked for 2 h using pre-immune donkey serum (PID, 1:5 in PBST 2%). Hereafter, retinal WMs were incubated with the primary antibody overnight in PBS with 2% Triton X-100 and 10% PID. Tissues were washed several times with PBST 0.5% before an incubation of 2 h with the secondary antibody (Alexa-647-conjugated donkey anti-rabbit IgG, 1:200, Thermofisher). Retinas from fish that were retrogradely traced using biocytin were incubated with streptavidin-Alexa488 for 2 h. All steps were carried out at room temperature. Retinas were mounted on glass slides with anti-fading mounting reagent Mowiol (10%, Sigma-Aldrich).
Finally, mosaic images of entire retinal WMs were imaged using a wide-field epifluorescent microscope (Leica DM6) containing a HC PL FLUOTAR L 20x/0.40 CORR objective (resolution = 2.17 pixels/μm).
2.6 Retinal whole-mount analysis
2.6.1 Establishment of automated counting models
The deep learning-based cell counting models for zebrafish and killifish RGCs were derived from RGCode () via transfer-learning. Briefly, frames of fixed size (177 by 177 μm) were obtained from Rbpms2-stained retinal WMs, sampling the central, mid and peripheral retina equally. The dataset was divided into training and testing datasets, and the cells were annotated by two independent expert counters. Zebrafish training dataset contained 48 frames, while the testing dataset was composed out of 36. For killifish, 36 and 24 frames were employed for the training and testing, respectively. Both counters annotated all testing frames, while the training ones were equally split between the two counters. Transfer-learning was performed starting from the weights of the original RGCode neural network, as described previously (). The performance of the newly-trained models was evaluated by linear-regression analysis against the human-annotated frames of the testing datasets. If the slope and coefficient of determination of the model were comparable to those between human counters, the performance was deemed satisfactory. Performance metrics of the final models are reported in Supplementary Figure 1.
2.6.2 Analysis
The number and density of RGCs per retina, together with the retinal area, were obtained by running the RGCode2 pipeline. Murine retinas were analyzed using the original RGCode model, while for zebrafish and killifish the newly-trained models were used. For retinal segmentation, the original RGCode model was used for all organisms, as its performance was satisfactory across all three animal species. RGCs density across the whole retina was obtained by dividing the total number of RGCs by the retinal area. Isodensity maps of the retina were generated as probability density functions via gaussian kernel density estimation (KDE). The KDE was generated from the centroids of the detected cells, with a bandwidth of 100 μm, and scaled to represent cells per square millimeter.
2.7 Statistical analysis
All data analyses were performed on raw, unsaturated micrographs. For visualization, some images were contrast-enhanced by adjusting the white point, applying identical, linear enhancements across comparable images. Statistical analyses included ANOVA, t-tests and U-tests, as described in the figure legends. The median with Kruskal–Wallis ANOVA and Mann-Whitney U-tests were used when data failed the Shapiro–Wilk normality test; otherwise, the mean with Welch ANOVA was chosen. All data processing, plotting, statistical analyses, were performed in Python using pandas, seaborn, matplotlib, and dabest (). A p-value < 0.05 was considered significant.
3 Results
3.1 Validation of Rbpms2 as a pan-RGC marker in the teleost retina
From the wide array of available murine pan-RGCs markers (), only a few are amenable for automated cell counting. Among these, RBPMS has gradually replaced BRN3A (Pou4f1, POU class 4 homeobox 1) over the last decade as the golden standard marker for evaluating RGC numbers (; ; ; ; ; Figures 1A–C). In contrast, rbpms2b has only recently been reported as a potential pan-RGC marker in zebrafish and killifish based on single cell RNA sequencing (; ; Figures 1D, E, G, H). While isl2b has been traditionally used for labeling zebrafish RGCs, it is not ubiquitously expressed in RGCs of zebrafish or killifish (; ; ; Figures 1F, I). To validate Rbpms2 as a pan-RGC marker for teleosts, we performed retrograde labeling assays known to label all RGCs in zebrafish and killifish. We observed near-complete co-labeling of biocytin and Rbpms2 in the zebrafish (Figure 1J) and killifish (Figure 1K) retina. However, biocytin labeling was heterogenous, with intensities ranging from faint to bright, whereas Rbpms2 provided uniform labeling of all RGCs. Additionally, Rbpms2 specifically marked RGC somas, unlike biocytin, which also labeled axons, masking the retinal ganglion cell layer and complicating quantification.
FIGURE 1
Thus, both RBPMS and Rbpms2 serve as reliable markers for assessing RGC densities and survival in murine and teleost retinas, respectively, and can be leveraged for automated cell counting platforms such as RGCode (
3.2 Comparative analysis of RGC density in the retina of mice, zebrafish and killifish
RGC densities were compared across species, including young adult mice (C57Bl/6N, 10 weeks-old), young adult zebrafish (AB, 21 weeks-old), and African turquoise killifish (GRZ-AD), a teleost gerontology model. In killifish, both young adult (6 weeks-old) and old (18 weeks-old) age groups were analyzed to assess age-related changes in RGC density. To this end, we developed RGCode2, an expansion of RGCode (
Scaled images of retinal WMs from mice, zebrafish and killifish reveal significant differences in retinal size (Figures 2A–D), confirmed by retinal area analysis based on automated segmentation by RGCode2 (Figure 2I). Notably, killifish exhibit substantial retinal expansion between 6 and 18 weeks (Figures 2C, D, I), consistent with pervious findings (
FIGURE 2

Different retinal ganglion cell (RGC) density in retinas of adult mice, zebrafish and killifish. (A–D) Scaled representation of retinas from young adult mice [10 weeks-old, (A)], young adult zebrafish [21 weeks-old, (B)], young adult killifish [6 weeks-old, (C)] and old killifish [18 week-old, (D)]. Mice have larger retinas than zebrafish and young killifish, while the retina of aged killifish is considerably larger than that of their younger counterparts. Scale bar 1 mm. (E–H) Representative micrographs of RGCs labeled with RBPMS [mouse, (E)] or Rbpms2 (fish), sampled from the temporal retina. Young adult zebrafish (F) and killifish (G) show a comparable density, higher than the one of old killifish (H) and young adult mice. Moreover, RGCs from the fish species are considerably smaller than the ones of mice. Scale bar 25 μm. (I) Automated quantification of the area of retinal WMs, revealing that unlike young adult fish, old killifish approach the size of murine retinas. (J) Automated quantification of RGC numbers in retinal WMs. Mice exhibit the lowest RGC count, with approximately 45,000 cells. In contrast, young fish possess around 70,000 RGCs. Aged killifish have the highest count, reaching approximately 125,000, nearly twice as many as young killifish. (K) Automated quantification of RGC density in retinal WMs. RGC density is considerably lower in mice compared to fish species. Notably, old killifish exhibit a significantly reduced RGC density compared to young adult fish. D, dorsal; N, nasal; RGCs, retinal ganglion cells; T, temporal; V, ventral; WMs, whole-mount.
On average, 10 weeks-old mouse retinas contain 44,499 ± 470 RGCs (Figure 2J), as previously reported (
RGC densities were calculated to reduce the impact of dissection artifacts and function as a more reliable interspecies comparison metric. RGC densities also serve as a proxy for visual acuity since it determines the degree of spatial detail that can be transmitted to the brain. Young adult zebrafish (21 weeks-old) and killifish (6 weeks-old) exhibit similar RGC densities of 17,390 ± 301 and 17,769 ± 272 RGCs/mm2, respectively, which are over five times higher than in young adult mice (3,179 ± 26 RGCs/mm2, Figures 2E–G, K). In old killifish (18 weeks-old), density declines to 14,041 ± 255 RGCs/mm2 (Figures 2H, K), but remains more than four times higher than in 10 weeks-old mice. Scaled isodensity projections (Supplementary Figure 2) confirm previously published density maps of mice (
Although the retina is highly conserved across vertebrates in terms of anatomy, including its layered structure and cell types, RGC densities vary greatly between species. These differences highlight the importance of cross-species comparisons of RGC function and responses to stimuli or injury, providing valuable insights into neural network dynamics and injury mechanisms.
3.3 Differential RGC susceptibility to optic nerve injury in canonical vertebrate models
Extensive literature describes a striking difference in intrinsic survival capacity of RGCs between mammals and teleost fish (Zou et al., 2013;
Based on a previously published RGC survival curve in mice (
FIGURE 3

Differential retinal ganglion cell (RGC) susceptibility to optic nerve injury in canonical vertebrate models. (A) Experimental timeline for the RGC survival experiment in young adult mice and zebrafish. RGC survival was evaluated at 7 and 14 days post-ONC in mice and at 7, 14, and 21 dpi in zebrafish. With assets from BioRender.com. (B) Representative micrographs of RPBMS-stained murine retinal WMs following ONC injury show substantial RGC loss at 7 dpi, which becomes even more pronounced by 14 dpi. Scale bar 25 μm. (C) Representative micrographs of Rpbms2-stained zebrafish retinal WMs following ONC injury. ONC injury leads to no appreciable loss of RGCs at 7 and 14 dpi, but a minor loss can be observed at 21 dpi. Scale bar 25 μm. (D) Quantification of RGC survival in adult murine WMs. ONC leads to the loss of over 50% of RGCs at 7 dpi, and a further one until 14 dpi, when only about 20% of the RGCs remain. (E) Quantification of RGC survival in adult zebrafish WMs. There is no significant loss of RGCs within the first 2 weeks after ONC. A small but significant decrease is measured at 21 dpi, with approximately 10% of the RGCs lost. Data from two independent experiments, presented as percentages relative to the median of uninjured retinas and presented as median ± 25–75th CI. Kruskal-Wallis ANOVA. P-values reported within the figure. CI, confidence interval; dpi, days post injury; ONC, optic nerve crush; RGCs, retinal ganglion cells; WMs, whole-mounts.
To conclude, we confirm the distinct resilience profiles of these two canonical vertebrate models, with mice displaying minimal neuroprotection and zebrafish exhibiting a significantly higher neuroprotective capacity.
3.4 Biphasic RGC loss in killifish after optic nerve crush injury
To assess the impact of aging on RGC survival, we used the African turquoise killifish, a well-established model in gerontology. Accordingly, we examined RGC survival in young adult (6 weeks-old) and old (18 weeks-old) female killifish (Figure 4A).
FIGURE 4

Biphasic retinal ganglion cell (RGC) loss in killifish after optic nerve crush injury. (A) Representative image of a young adult (6 weeks-old) and old killifish (18 weeks-old) and experimental timeline for the RGC survival experiment, where RGC survival is evaluated at 4, 7, 14, and 21 days following ONC injury. Scale bar: 1 mm. (B) Representative micrographs of Rbpms2-stained WMs of young adult (6 weeks-old) and old killifish (18 weeks-old). For both ages, an appreciable loss of RGCs is evident at 7 dpi with further loss at 21 dpi, when compared to uninjured age-matched control fish. (C) Quantification of RGC survival in adult killifish WMs shows a first wave of RGC loss at 4 dpi, with 20% of the RGCs lost in both age groups, and this loss remains steady through 7 dpi. A second wave of loss is observed at 14 dpi, with older fish losing more RGCs (50%) compared to young fish (40%). No further loss is detected at 21 dpi in either age group. (D) Quantification of RGC survival per retinal quadrant in adult killifish WMs reveals no significant differences in inter-quadrant RGC loss after ONC by 21 dpi for both young adult (6 weeks-old) and old (18 weeks-old) fish. Data from two independent experiments, presented as percentages relative to the median of their uninjured age-matched control and presented as median ± 25–75th CI. Two-way Kruskal-Wallis ANOVA with pairwise Mann-Whitney U tests (C), One-way Kruskal-Wallis ANOVA and post hoc Mann-Whitney U test with Bonferroni correction [(D), 06 weeks], One-way Welch ANOVA and post hoc Games-Howell test [(D), 18 weeks]. p-values reported within the figure for significant differences. CI, confidence interval; DN, dorsonasal; dpi, days post injury; DT, dorsotemporal; ONC, optic nerve crush; RGCs, retinal ganglion cells; VN, ventronasal; VT, ventrotemporal; WMs, whole-mounts.
The killifish retina, however, undergoes rapid expansion during early life stages, which is driven by a dynamic balance between cell addition and tissue stretching (
In contrast to zebrafish, both young adult and old killifish exhibit a rapid loss of RGCs following ONC, with approximately 79% (IQR: 75%, 82%) and 83% (IQR: 75%, 88%) of RGCs surviving at 4 dpi, respectively, a level that remains stable until 7 dpi (Figures 4B, C). This initial phase of RGC loss mirrors the pattern observed in mice (
These findings indicate that while both age groups exhibit a biphasic pattern of RGC loss, aging is a key factor determining the overall extent of degeneration, primarily due to the increased RGC loss observed during the second wave.
4 Discussion and conclusion
In this study, we validated the use of Rbpms2 as a pan-RGC marker to study ganglion cell numbers in both zebrafish and killifish. Additionally, we present an updated version of the automated RGC counting platform, RGCode (
Despite the conserved cytoarchitecture of the vertebrate retina, from lampreys to humans (
RGC survival in the mouse retina has been extensively studied, demonstrating an acute sigmoidal monophasic loss pattern during the first 3 months following injury (
Unlike mice, zebrafish exhibit strong neuroprotective properties, showing no RGC loss until 14 dpi and only a minor reduction in RGC density by 21 dpi. This loss is unlikely to result directly from injury (primary degeneration), as RGC transcriptional profiles at 14dpi closely resemble those of uninjured cells (Zhang et al., 2025). As such, it is more conceivable that this moderate loss of RGCs reflects secondary degeneration. Immune system overactivation is one plausible cause, as immune cell abundance peaks at 7 dpi, but only returns to baseline levels by 28 dpi (
Killifish exhibit strong CNS regenerative capacity in early adulthood, but this ability declines with age, approaching a mammalian-like phenotype at old age (
In summary, cross-species comparisons of RGC survival following ONC reveal distinct resilience patterns. Mice exhibit minimal neuroprotection, whereas zebrafish display robust neuroprotective capacities. Killifish, with their biphasic RGC loss profile, potentially offer a unique model to study both intrinsic injury susceptibility (primary degeneration) and cell loss driven by reduced electrical activity and/or detrimental factors in the retinal microenvironment (secondary degeneration). As such, these findings provide a framework for elucidating molecular mechanisms that regulate RGC survival and for identifying potential therapeutic targets to enhance neuroprotection, which could aid in the treatment of optic neuropathies, and by extension of neurodegenerative disorders and traumatic CNS injuries.
Statements
Data availability statement
RGCode2 is available on the public GitLab repository https://gitlab.com/NCDRlab/rgcode2, other data will be made available upon request.
Ethics statement
The animal study was approved by KU Leuven Animal Ethics Committee. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
JDDS: Writing – original draft, Writing – review and editing, Methodology. AZ: Methodology, Writing – review and editing, Writing – original draft. P-JS: Writing – original draft, Funding cquisition, Methodology, Writing – review and editing. LMo: Funding acquisition, Supervision, Conceptualization, Writing – review and editing. LMa: Funding acquisition, Writing – review and editing, Visualization, Writing – original draft, Supervision, Conceptualization, Methodology. SB: Methodology, Supervision, Conceptualization, Writing – original draft, Visualization, Funding acquisition, Writing – review and editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. P-JS, SB, and LMa hold a personal Research Foundation Flanders (FWO, Belgium) fellowship (1114525N, 1165020N, and 1S42720N, respectively). Experiments are financially supported by KU Leuven (C14/22/074, KAC14/22/074 and KA16-00745) and FWO (G082221N, G092222N).
Acknowledgments
We would like to thank Simon Buys, Arnold Van Den Eynde, Evelien Herinckx, Véronique Brouwers for the daily animal maintenance and environmental control; Iene Kemps, Dr. Lien Andries, Dr. Marie Claes and Marijke Christiaens for their technical support.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2025.1596464/full#supplementary-material
Supplementary Table 1Mean retinal ganglion cell (RGC) count, retinal area and RGC density of non-injured retinas per species and age, along with their relative standard error of the mean (SEM).
Supplementary Table 2Median retinal ganglion cell (RGC) count, retinal area and RGC density per species, age and timepoint after optic nerve crush (ONC), along with their relative 25th and 75th percentiles. For killifish, the values of age-matched control were included as well.
Supplementary Table 3Results of Kruskal-Wallis two-way ANOVA between injured young adult and old killifish. Multiple comparisons statistics and p-values are reported.
Abbreviations
AMC, age-matched controls; CNS, central nervous system; DN, dorsonasal; dpi, days post-injury; DT; dorsotemporal; FWO, fonds wetenschappelijk onderzoek; KDE, kernel density estimation; IQR, interquartile range; Isl2b, ISL LIM Homeobox 2; ONC, optic nerve crush; PBS, phosphate buffered saline; PFA, paraformaldehyde; PID, pre-immune donkey serum; Pou4f1, POU Class 4 Homeobox 1; RBPMS, RNA-binding protein with multiple splicing; Rbpms2, RNA-binding protein with multiple splicing 2; SEM, standard error of mean; VN, ventronasal; VT, ventrotemporal; WM, whole-mount.
References
1
BadenT. (2024a). Ancestral photoreceptor diversity as the basis of visual behaviour.Nat. Ecol. Evol.83374–386. 10.1038/s41559-023-02291-7
2
BadenT. (2024b). The vertebrate retina: A window into the evolution of computation in the brain.Curr. Opin. Behav. Sci.57:101391. 10.1016/J.COBEHA.2024.101391
3
BadenT.EulerT.BerensP. (2020). Understanding the retinal basis of vision across species.Nat. Rev. Neurosci.215–20. 10.1038/s41583-019-0242-1
4
BeckersA.BergmansS.Van DyckA.MoonsL. (2023). Analysis of axonal regrowth and dendritic remodeling after optic nerve crush in adult zebrafish.Methods Mol. Biol.2636163–190. 10.1007/978-1-0716-3012-9_9/FIGURES/2
5
BeckersA.Van DyckA.BollaertsI.Van houckeJ.LefevereE.AndriesL.et al (2019). An antagonistic axon-dendrite interplay enables efficient neuronal repair in the adult zebrafish central nervous system.Mol. Neurobiol.563175–3192. 10.1007/s12035-018-1292-5
6
BeckersA.VanhunselS.Van DyckA.BergmansS.MasinL.MoonsL. (2021). Injury-induced autophagy delays axonal regeneration after optic nerve damage in adult zebrafish.Neuroscience47052–69. 10.1016/J.NEUROSCIENCE.2021.07.009
7
BergmansS.NoelN. C. L.MasinL.HardingE. G.KrzywańskaA. M.De SchutterJ. D.et al (2024). Age-related dysregulation of the retinal transcriptome in African turquoise killifish.Aging Cell23:e14192. 10.1111/acel.14192
8
BergmansS.SerneelsP. J.MasinL.MoonsL. (2023a). Tissue stretching is a confounding factor for the evaluation of neurodegeneration in the fast-ageing killifish.Biogerontology24403–419. 10.1007/s10522-023-10026-1
9
BergmansS.VanhunselS.Van HouckeJ.MariënV.ArckensL.MoonsL. (2023b). Morphological analysis of aging hallmarks in the central nervous system of the fast-aging african turquoise killifish.Cold Spring Harb. Protoc.2023739–754. 10.1101/pdb.prot107827
10
BiehlmaierO.NeuhaussS. C. F.KohlerK. (2001). Onset and time course of apoptosis in the developing zebrafish retina.Cell Tissue Res.306199–207. 10.1007/S004410100447/METRICS
11
BoiaR.RuzafaN.AiresI. D.PereiroX.AmbrósioA. F.VecinoE.et al (2020). Neuroprotective strategies for retinal ganglion cell degeneration: Current status and challenges ahead.Int. J. Mol. Sci.21:2262. 10.3390/IJMS21072262
12
ChenS.LathropK. L.KuwajimaT.GrossJ. M. (2022). Retinal ganglion cell survival after severe optic nerve injury is modulated by crosstalk between Jak/Stat signaling and innate immune responses in the zebrafish retina.Development149:dev199694. 10.1242/dev.199694
13
ClaesM.De GroefL.MoonsL. (2019). Target-Derived neurotrophic factor deprivation puts retinal ganglion cells on death row: Cold hard evidence and caveats.Int. J. Mol. Sci.20:4314. 10.3390/IJMS20174314
14
ClaesM.MoonsL. (2022). Retinal ganglion cells: Global number, density and vulnerability to glaucomatous injury in common laboratory mice.Cells11:2689. 10.3390/CELLS11172689
15
De GroefL.DekeysterE.GeeraertsE.LefevereE.StalmansI.Salinas-NavarroM.et al (2016). Differential visual system organization and susceptibility to experimental models of optic neuropathies in three commonly used mouse strains.Exp. Eye Res.145235–247. 10.1016/j.exer.2016.01.006
16
FleischV. C.FraserB.AllisonW. T. (2011). Investigating regeneration and functional integration of CNS neurons: Lessons from zebrafish genetics and other fish species.Biochim. Biophys. Acta - Mol. Basis Dis.1812364–380. 10.1016/j.bbadis.2010.10.012
17
GuoL.DavisB. M.RavindranN.GalvaoJ.KapoorN.HaamediN.et al (2020). Topical recombinant human Nerve growth factor (rh-NGF) is neuroprotective to retinal ganglion cells by targeting secondary degeneration.Sci. Rep.10:3375. 10.1038/s41598-020-60427-2
18
HahnJ.MonavarfeshaniA.QiaoM.KaoA. H.KölschY.KumarA.et al (2023). Evolution of neuronal cell classes and types in the vertebrate retina.Nature624415–424. 10.1038/s41586-023-06638-9
19
HoJ.TumkayaT.AryalS.ChoiH.Claridge-ChangA. (2019). Moving beyond P values: Data analysis with estimation graphics.Nat. Methods16565–566. 10.1038/s41592-019-0470-3
20
HoangT.WangJ.BoydP.WangF.SantiagoC.JiangL.et al (2020). Gene regulatory networks controlling vertebrate retinal regeneration.Science370:eabb8598. 10.1126/science.abb8598
21
JohnsonE. C.GuoY.CepurnaW. O.MorrisonJ. C. (2009). Neurotrophin roles in retinal ganglion cell survival: Lessons from rat glaucoma models.Exp. Eye Res.88808–815. 10.1016/J.EXER.2009.02.004
22
JonesK. A.MaltbyS.PlankM. W.KlugeM.NilssonM.FosterP. S.et al (2018). Peripheral immune cells infiltrate into sites of secondary neurodegeneration after ischemic stroke.Brain. Behav. Immun.67299–307. 10.1016/J.BBI.2017.09.006
23
KanedaM.NagashimaM.NunomeT.MuramatsuT.YamadaY.KuboM.et al (2008). Changes of phospho-growth-associated protein 43 (phospho-GAP43) in the zebrafish retina after optic nerve injury: A long-term observation.Neurosci. Res.61281–288. 10.1016/J.NEURES.2008.03.008
24
KölschY.HahnJ.SappingtonA.StemmerM.FernandesA. M.HelmbrechtT. O.et al (2021). Molecular classification of zebrafish retinal ganglion cells links genes to cell types to behavior.Neuron109645–662.e9. 10.1016/j.neuron.2020.12.003
25
Levkovitch-VerbinH.QuigleyH. A.MartinK. R. G.ZackD. J.PeaseM. E.ValentaD. F. (2003). A model to study differences between primary and secondary degeneration of retinal ganglion cells in rats by partial optic nerve transection.Invest. Ophthalmol. Vis. Sci.443388–3393. 10.1167/IOVS.02-0646
26
LiJ.ChoiJ.ChengX.MaJ.PemaS.SanesJ. R.et al (2024). Comprehensive single-cell atlas of the mouse retina.iScience27:109916. 10.1016/j.isci.2024.109916
27
LukowskiS. W.LoC. Y.SharovA. A.NguyenQ.FangL.HungS. S.et al (2019). A single-cell transcriptome atlas of the adult human retina.EMBO J.38:e100811. 10.15252/embj.2018100811
28
MangrumW. I.DowlingJ. E.CohenE. D. (2002). A morphological classification of ganglion cells in the zebrafish retina.Vis. Neurosci.19767–779. 10.1017/S0952523802196076
29
MariënV.Van HouckeJ.ArckensL. (2022). Intracardial perfusion of the African turquoise killifish.2–7. 10.17504/protocols.io.b2ryqd7w
30
MarquesI. J.LupiE.MercaderN. (2019). Model systems for regeneration: Zebrafish.Dev.1461–13. 10.1242/dev.167692
31
MasinL.ClaesM.BergmansS.CoolsL.AndriesL.DavisB. M.et al (2021). A novel retinal ganglion cell quantification tool based on deep learning.Sci. Rep.11:702. 10.1038/s41598-020-80308-y
32
Nadal-NicolásF. M.Galindo-RomeroC.Lucas-RuizF.Marsh-AmstrongN.LiW.Vidal-SanzM.et al (2023). Pan-retinal ganglion cell markers in mice, rats, and rhesus macaques.Zool. Res.44226–248. 10.24272/J.ISSN.2095-8137.2022.308
33
O’Hare DoigR. L.BartlettC. A.MaghzalG. J.LamM.ArcherM.StockerR.et al (2014). Reactive species and oxidative stress in optic nerve vulnerable to secondary degeneration.Exp. Neurol.261136–146. 10.1016/J.EXPNEUROL.2014.06.007
34
OyinboC. A. (2011). Secondary injury mechanisms in traumatic spinal cord injury: A nugget of this multiply cascade.Acta Neurobiol. Exp.71281–299. 10.55782/ANE-2011-1848
35
ParkE.VelumianA. A.FehlingsM. G. (2004). The role of excitotoxicity in secondary mechanisms of spinal cord injury: A review with an emphasis on the implications for white matter degeneration.21754–774. 10.1089/0897715041269641
36
Ramon y CajalS. (1893). La Rétine des Vertébrés.La Cellule9119–257.
37
Ramos-CejudoJ.WisniewskiT.MarmarC.ZetterbergH.BlennowK.de LeonM. J.et al (2018). Traumatic brain injury and Alzheimer’s disease: The cerebrovascular link.EBioMedicine2821–30. 10.1016/j.ebiom.2018.01.021
38
RheaumeB. A.JereenA.BolisettyM.SajidM. S.YangY.RennaK.et al (2018). Single cell transcriptome profiling of retinal ganglion cells identifies cellular subtypes.Nat. Commun.9:2759. 10.1038/s41467-018-05134-3
39
RodriguezA. R.de Sevilla MüllerL. P.BrechaN. C. (2014). The RNA binding protein RBPMS is a selective marker of ganglion cells in the mammalian retina.J. Comp. Neurol.5221411–1443. 10.1002/CNE.23521
40
Sánchez-MigallónM. C.Valiente-SorianoF. J.Salinas-NavarroM.Nadal-NicolásF. M.Jiménez-LópezM.Vidal-SanzM.et al (2018). Nerve fibre layer degeneration and retinal ganglion cell loss long term after optic nerve crush or transection in adult mice.Exp. Eye Res.17040–50. 10.1016/J.EXER.2018.02.010
41
ShvetcovA.RuitenbergM. J.DelerueF.GoldW. A.BrownD. A.FinneyC. A. (2023). The neuroprotective effects of estrogen and estrogenic compounds in spinal cord injury.Neurosci. Biobehav. Rev.146:105074. 10.1016/J.NEUBIOREV.2023.105074
42
SkarieJ.ColleryR. F.LinkB. A. (2015). A transgenic zebrafish model to study retinal ganglion cell axon regeneration.Invest. Ophthalmol. Vis. Sci.56:4966.
43
SkarieJ. M.NickellsR. W. (2016). Mechanisms of retinal ganglion cell survival in zebrafish after axon injury.Invest. Ophthalmol. Vis. Sci.57:2555.
44
SmolenP.DashP. K.RedellJ. B. (2023). Traumatic brain injury-associated epigenetic changes and the risk for neurodegenerative diseases.Front. Neurosci.17:1259405. 10.3389/FNINS.2023.1259405
45
StirlingD. P.CumminsK.Wayne ChenS. R.StysP. (2014). Axoplasmic reticulum Ca2+ release causes secondary degeneration of spinal axons.Ann. Neurol.75220–229. 10.1002/ANA.24099
46
TapiaM. L.Nascimento-dos-SantosG.ParkK. K. (2022). Subtype-specific survival and regeneration of retinal ganglion cells in response to injury.Front. Cell Dev. Biol.10:956279. 10.3389/FCELL.2022.956279
47
TranN. M.ShekharK.WhitneyI. E.JacobiA.BenharI.HongG.et al (2019). Single-Cell profiles of retinal ganglion cells differing in resilience to injury reveal neuroprotective genes.Neuron1041039–1055.e12. 10.1016/j.neuron.2019.11.006
48
Van houckeJ.BollaertsI.GeeraertsE.DavisB.BeckersA.Van HoveI.et al (2017). Successful optic nerve regeneration in the senescent zebrafish despite age-related decline of cell intrinsic and extrinsic response processes.Neurobiol. Aging601–10. 10.1016/j.neurobiolaging.2017.08.013
49
Van houckeJ.GeeraertsE.VanhunselS.BeckersA.NoterdaemeL.ChristiaensM.et al (2019). Extensive growth is followed by neurodegenerative pathology in the continuously expanding adult zebrafish retina.Biogerontology20109–125. 10.1007/s10522-018-9780-6
50
Van houckeJ.MariënV.ZandeckiC.VanhunselS.MoonsL.AyanaR.et al (2021). Aging impairs the essential contributions of non-glial progenitors to neurorepair in the dorsal telencephalon of the Killifish Nothobranchius furzeri.Aging Cell20:e13464. 10.1111/acel.13464
51
VanhunselS.BergmansS.BeckersA.EtienneI.Van BergenT.De GroefL.et al (2022a). The age factor in optic nerve regeneration: Intrinsic and extrinsic barriers hinder successful recovery in the short-living killifish.Aging Cell21:e13537. 10.1111/acel.13537
52
VanhunselS.BergmansS.MoonsL. (2022b). Killifish switch towards mammalian-like regeneration upon aging.Aging142924–2925. 10.18632/aging.203995
53
VanhunselS.BergmansS.BeckersA.EtienneI.Van houckeJ.SeuntjensE.et al (2021). The killifish visual system as an in vivo model to study brain aging and rejuvenation.NPJ Aging Mech. Dis.7:22. 10.1038/s41514-021-00077-4
54
VanhunselS.BergmansS.MoonsL. (2023). The optic nerve crush injury paradigm in african turquoise killifish to study axonal regeneration in an aged environment.Cold Spring Harb. Protoc.2023872–884. 10.1101/pdb.prot107828
55
VaradarajanS. G.HunyaraJ. L.HamiltonN. R.KolodkinA. L.HubermanA. D. (2022). Central nervous system regeneration.Cell18577–94. 10.1016/J.CELL.2021.10.029
56
WiseP. M.DubalD. B.WilsonM. E.RauS. W.BöttnerM. (2001). Minireview: Neuroprotective effects of estrogen—new insightsinto mechanisms of action.Endocrinology142969–973. 10.1210/ENDO.142.3.8033
57
WuJ.CaiY.WuX.YingY.TaiY.HeM. (2021). Transcardiac perfusion of the mouse for brain tissue dissection and fixation.Bio-Protocol11:e3876. 10.21769/bioprotoc.3876
58
YangN.YoungB. K.WangP.TianN. (2020). The susceptibility of retinal ganglion cells to optic nerve injury is type specific.Cells9:677. 10.3390/CELLS9030677
59
ZhangA.BergmansS.DyckA.Van, BeckersA.MoonsL.et al (2025). Successful axonal regeneration is driven by evolutionarily conserved metabolic reprogramming.bioRxiv [Preprint]10.1101/2025.02.12.637870
60
ZouS.TianC.GeS.HuB. (2013). Neurogenesis of retinal ganglion cells is not essential to visual functional recovery after optic nerve injury in adult zebrafish.PLoS One8:e57280. 10.1371/JOURNAL.PONE.0057280
Summary
Keywords
retina, retinal ganglion cell, survival, optic nerve crush, mouse, killifish, zebrafish, neuroprotection
Citation
De Schutter JD, Zhang A, Serneels P-J, Moons L, Masin L and Bergmans S (2025) Differential retinal ganglion cell resilience to optic nerve injury across vertebrate species. Front. Neurosci. 19:1596464. doi: 10.3389/fnins.2025.1596464
Received
19 March 2025
Accepted
28 April 2025
Published
20 May 2025
Volume
19 - 2025
Edited by
Pete A. Williams, Karolinska Institutet (KI), Sweden
Reviewed by
Caridad Galindo-Romero, University of Murcia, Spain
Sabrina Reinehr, Ruhr-University Bochum, Germany
Updates

Check for updates
Copyright
© 2025 De Schutter, Zhang, Serneels, Moons, Masin and Bergmans.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lieve Moons, lieve.moons@kuleuven.be
†These authors have contributed equally to this work and share first authorship
‡These authors have contributed equally to this work and share last authorship
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.