Novel laser model allows us to perform accurate and reproducible optic nerve transection
To develop a new surgical model, we used the Leica LMD7 laser microdissection platform and Xla.Tg (tubb2b:mapt-GFP)Amaya transgenic X. laevis line to visualize the optic nerve cut and the subsequent axonal regrowth. X. laevis tadpoles were used at developmental stages 43–46 because at these stages the tadpoles are transparent and suitable for live-fluorescent imaging (Fig. 1A). Tadpoles were anesthetized and loaded onto custom surgical beds (Fig. 1A’). The left optic nerve of each tadpole was transected with a 349 nm UV cutting laser at 5× magnification, 300 µm distal to the retina (using the “measure” tool in the Leica LMD7 software). The cut can be confirmed by a gap in the fluorescent optic nerve visualization (Fig. 1B, 1C). The tadpoles were immediately removed, recovered in 1X MMR media, and used in later experiments. Our survival rate with this model is 97.3% over many experiments, and nerves can be cut at precisely the same distance (Fig. 1C). Following our novel laser transection injury, we observed with longitudinal confocal imaging that axons distal to the injury site degenerate from day 0–3 post laser (PL), followed by regeneration of the optic nerve from the globe-proximal cut end, which is observed from the third day after optic nerve transection (Fig. 1D). After 14 days, the regrown optic nerve is anatomically and morphologically indistinguishable from the naïve (uninjured) optic nerve (Fig. 1D).
Pre-existing RGCs do not die and new RGCs are not generated from retinal progenitors after optic nerve transection
It is known that damage to axons can lead to RGC death, and that dying RGCs in amphibian retinas can trigger RGC differentiation from progenitors for replacement[10–15]. Thus, we explored whether in our model damage to axons leads to death of pre-existing RGCs, and if axonal regeneration is the result of sprouting of axons from newly differentiated RGCs. Since amphibian Müller glia and RPE undergo a reprogramming process and generate a pool of proliferating progenitor cells, we examined these cell types as a potential source of newborn RGCs after optic nerve injury in the model. To this end, the left optic nerves were cut (the right optic nerves were used as controls). Before collecting control and experimental retinas 1, 2, 3, 6, and 10 days after injury, the tadpoles swam for 4 hours (pulse) in EdU-containing tadpole buffer (1X MMR). The retinas were then collected for analysis of proliferating cells. After 1, 2, and 3 days after injury, there was staining for EdU in the ciliary marginal zone, the region of the tadpole retina which is responsible for normal retinal growth, which was equal in the control as well as the injured retina (n = 4, Fig. 2A, 2B). We found no EdU staining in the ciliary marginal zone in experimental and control tadpoles 10 days after damage (Fig. 2A). We also did not detect any staining in those areas of the experimental and control retinas where Müller glia and RPE are located (Fig. 2A). Since we found dividing cells only in the ciliary marginal zone, progenitors from this zone could be responsible for the appearance of newborn RGCs (Fig. 2B). To test this hypothesis, we used the following design. The optic nerve was cut and after 0, 3 or 7 days the tadpoles were placed in EdU- containing tadpole buffer for 3 days. We expected that such a period would be enough for cells differentiated from progenitors to take their place in the ganglion cell layer of the retina. We found only a small number of EdU-positive newborn cells in the ganglion cell layer, which did not differ between experimental and control retinas (n = 4, Fig. 2C). Next, we examine cell death in retinal slices using a TUNEL assay. To this end, the left optic nerves were cut, and the retinas were collected after 1, 2, 3, 6, and 10 days. The right optic nerves were used as controls. We also used TUNEL-positive controls (Fig. 2D). We found only a small number of TUNEL-positive cells (maximum 3 cells per retinae and sometimes no cells at all, n = 4) in both experimental and control retinas (Fig. 2E).
RGC axons sprout in all directions during an early stage of optic nerve regeneration, but the plasticity and mobility of RGC axons allows them to restore a full-fledged optic nerve at a late stage
We were unsure if the presence of a chemoattractant gradient could still be expected in our tadpoles at this stage of development. If a chemoattractant gradient were present, then the RGC axons should grow in a certain direction and not in all possible directions. To this end, we cut the optic nerve and observe how it regenerates in a time-dependent manner. To examine the patterns of axon growth and remodeling, we conducted Sholl analysis at 3, 6, and 14 days after injury (Fig. 3A). Axons were semi-automatically traced with the auto-snapping feature in Neuroanatomy/SNT plugin and then the Sholl analysis plugin was used to assess branching complexity by measuring the frequency at which the traced axons intersect concentric circles spaced at regular intervals from the origin (in our case, from the cut site). At 100 µm from the cut site, the number of intersections at day 3 was 11 ± 1 (n = 5 regrowing nerves), for day 6 it was 12 ± 1 (n = 5), and for day 14 it was at 2 ± 1 (n = 10). However, while the shape at day 3 and 6 after injury indicates many axons branching out from near the cut site, there are fewer branch intersection points the farther out from the cut site you go (Fig. 3A). We also noticed that there are many axons during these periods of time that form loops or randomly change their direction of movement (Fig. 3B). All these data indicate that a chemoattractant gradient is not present during an early stage of optic nerve regeneration. By day 14, however, there is only one line and the initial peak of intersections is gone (Fig. 3A). To better understand the process of formation of a full-fledged optic nerve after injury, we quantified the thickness of the regrowing nerve where it meets the brain compared to the contralateral nerve (Fig. 3C). At day 1 after injury, the percentage of the contralateral nerve was 55 ± 35% (n = 3, a large variation was observed in the sizes of the degenerating distal segment), and this was observed to be the degenerating distal segment. At day 2 the average was 1 ± 1% (n = 5) and the average of day 3 was 8 ± 3% (n = 6). At day 4 the average was 19 ± 6% (n = 6), and at day 6 the average was 28 ± 6% (n = 6). At day 10 the average was 70 ± 17% (n = 8). At day 14 the average was 87 ± 11% of the contralateral (n = 10). One would expect that in the absence of a gradient at a later stage, the emerging new optic nerve would be less organized. However, our data indicate that the new optic nerve does not differ from the control (contralateral) optic nerve at 14 days after injury (Fig. 3D). These results suggest that at a later stage, attraction occurs from axons that have established contact with the brain. This attraction may be caused by the secretion of chemoattractants from these axons and the formation of a gradient. This may also be caused by the fact that the axons which do not make it to the correct place in the brain are either rerouted or trimmed back and regrown.
Regrown nerve resembles the contralateral uninjured nerve in 3D spatial positioning, resulting in functional vision restoration
After observing consistent regrowth and remodeling of axons from pre-existing RGCs, we wanted to see if this robust optic nerve regeneration with our model resulted in useful vision. First, we checked if the regrown nerve eventually resembles the contralateral nerve in 3D spatial positioning. We performed depth-coded imaging and noticed that at 14 days, when the nerve most consistently resembles the contralateral, the depth of the regrown nerve was approximately the same and the point at which the nerve intersects with the brain was also in the same orthogonal plane (Fig. 4A, 4B). We designed an analysis to observe volume colocalization with the contralateral nerve overlayed on the mirror plane in 3D space using Imaris image analysis software (Fig. 4C). This was not a flawless system, because animals may not be biologically perfectly symmetrical with their left and right side, but since the tadpole is rapidly growing, this internal control was better than a pre-injury and post-injury overlay analysis. We used a naïve age-matched cohort to examine how much volume ratio overlap is typical in the left and right optic nerves with this imaging method, and if there were any differences when either one or both nerves are cut and regrown. There were no statistical differences at day 14 between single cut, double cut or naïve cohorts (Fig. 4D). These results indicate that even when both optic nerves are cut and there is no chemoattractant gradient early on, RGC axons have internal information that allows them to establish contact with the correct part of the brain. This statement is confirmed by our behavioral tests, according to which tadpoles practically restore vision on the 14th day after optic nerve injury (Fig. 4E-4G). We took advantage of the fact that sighted X. laevis tadpoles prefer a white background when given the choice between a white and black background (Fig. 4E) [16]. For all behavioral tests, we performed bilateral optic nerve transections to ensure complete blindness. We recorded swimming videos from naïve, as well as 3, 7, 14 and 21 days after injury. The naïve, uninjured tadpoles showed a clear preference for the white side of the tank (Fig. 4F). They spent an average of 88 ± 3% (n = 12) of the time on the light side of the tank. At day 3 and day 7 after injury, tadpoles exhibited no preference and swam back and forth across the black and white sides of the tank relatively equally (day 3: 48 ± 4%, n = 8; day 7: 45 ± 3%, n = 8; Fig. 4F, 4G). At day 14 and 21, however, the tadpoles again showed a preference for the white side of the tank (day 14: 75 ± 5%, n = 8; day 21 73 ± 9%, n = 5). The naïve percentage was significantly different (p < 0.0001) from both day 3 and day 7 post injury. However naive vs. day 14 and naïve vs day 21 were not significantly different. This is consistent with tadpoles losing vision after double injury, being blind at day 3 and day 7 and having regained vision by day 14 and day 21. Considering these results and the data presented above, we can propose the following model. Injury to the optic nerve results in the sprouting of axons from pre-existing RGCs in all directions. The RGC axons, upon reaching the brain, can recognize whether this part is correct or incorrect. If the axons reach the correct part of the brain, they strengthen contact with it and possibly create a gradient of attracting other axons. If the axons reach the wrong part of the brain, they are either rerouted or trimmed back and regrown. There were a few exceptions to this rule, and we have noticed on occasion axons which grew to the wrong location and strengthened contact to the incorrect location, such as towards another nerve, but these cases were very rare (Fig. 5).
The optic nerve transection has a slight effect on the retina at the molecular level
Our morphological observations of proliferation and death of cells in the retina following optic nerve transection indicated the absence of significant changes compared to control retinas. However, how do these observations fit with changes in the retina at the molecular level? To answer this question, we performed RNA-seq analysis of the control and experimental retinas 3 days after the optic nerve transection. We chose this time point because at 3 days post optic nerve transection we observed significant axonal sprouting, which would be reflected in changes in expression in the retina. Due to the small quantity of material, we pooled 8 retinas to prepare one RNA sample (Fig. 6A). In total, we prepared 4 experimental and 4 control (retinas from the contralateral eyes) RNA samples that were used in the RNA-seq analysis. Since the retina contains many different cell types, we used deep sequencing to detect changes in gene expression in small populations of retinal cell types, including RGCs. To this end, 44,930,137 ± 2,304,700 fragments were sequenced on average per library among which 40,233,947 ± 2,096,169 fragments were uniquely mapped to the mouse genome. Our results indicate that counts per million (CPM) distributions do not differ between experimental and control retinas (Fig. 6B). While the value of the correlation coefficient (0.987), principal component analysis (PCA), heatmap and volcano plots show differences in expression levels between experimental and control retinas, the number of genes that showed a two-fold change in expression was small (Fig. 6C-6F, Supplementary Data S1). The number of genes whose expression was twice as high (log2 fold change [log2FC] ≥ 1) in the experimental retinas compared to the controls was 44 (FDR < 0.1); the number of genes whose expression was reduced more than twofold (log2FC ≤ − 1) in the experimental retinas compared to the control retinas was 38 (Supplementary Data S1). The gene ontology (GO) enrichment analysis revealed that genes with increased expression are involved in basic metabolic processes or in processes associated with the activity of neurons and glial cells (Fig. 6G, Supplementary Data S1). Among the genes with increased expression, gap43 (log2FC = 1.43, FDR = 0.014) and uchl1 (log2FC = 1.32, FDR = 0.073) are worth noting, the role of which in the process of axon regeneration after damage to the optic nerve has already been shown (Fig. 6E, Supplementary Data S1)[9]. It is interesting to note that the expression of clu (clusterin), the activity of which is characteristic of reactive astrocytes and Muller glia, was increased (log2FC = 1.25, FDR = 0.095). Meanwhile, Muller glia marker glul (glutamine synthetase) expression was reduced (log2FC=-1.23, FDR = 0.029). We also noticed reduced expression of a RGC marker rbpms2 (log2FC=-1.11, FDR = 0.014). However, we were unable to identify significant processes in which genes with reduced expression are involved. It should also be noted that we did not detect activity of signaling cascades associated with inflammation and cell death. Overall, our results could be described as that regeneration of RGC axons occurs in the resting state of the retina.
There is no significant inflammatory response in the optic nerve after injury, but the activity of signaling cascades involved in cell proliferation is increased
The lack of significant changes in gene expression in the experimental retinas compared to control retinas prompted us to examine changes in gene expression in the optic nerves at the same time point (3 days post injury). To this end, tadpole nerves were lasered on the left side and 3 days later (the time of axon branching) left (experimental) and right (control) optic nerves were dissected. Due to the very low amount of RNA in nerves, we pooled 15 nerves per sample, 3 samples per condition - either lasered or control (Fig. 7A). These samples were used for RNA-seq analysis. While CPM distributions did not differ between experimental and control optic nerves, we found a significant difference in gene expression between experimental and control optic nerves as inferred from the value of the correlation coefficient (0.968), principal component analysis (PCA), and volcano plot (Fig. 7B-7E). We found that the expression of 457 genes was two-fold or higher in experimental optic nerves compared to control optic nerves, while the expression of 276 genes was two-fold or lower (FDR < 0.1, Fig. 7F, Supplementary Data S1). The lists of these genes were used in the GO enrichment analysis. We found that a significant number of genes upregulated in damaged optic nerves were involved in cell proliferation (Fig. 7G, Supplementary Data S2). This result may be explained by increased expression of genes involved in gliogenesis (e.g., GO:0014013 regulation of gliogenesis [FDR = 0.009]), which may be associated with the restoration of the optic nerve structure. Genes whose expression is reduced are mainly involved in synaptogenesis (Fig. 7G, Supplementary Data S2). The reduced expression of these genes may be explained in part by the fact that new axons have not yet established contact with the brain at this time point (3 days post injury). It was surprising to us that among the biological processes with a high enrichment FDR, there were none that were involved in the inflammatory response. Several processes involved in the inflammatory response can be found with low FDR (e.g., GO:0006954 inflammatory response [FDR = 0.011], GO:0002534 cytokine production involved in inflammatory response [FDR = 0.025], GO:0050728 negative regulation of inflammatory response [FDR = 0.035], Supplementary Data S2). However, the genes involved in these processes and having increased expression in the damaged optic nerves play a supporting rather than a key role (Supplementary Data S2). The absence of a toxic inflammatory response in the optic nerve after injury likely promotes axonal regeneration.