Sphingosine-1-phosphate signaling regulates the ability of Müller glia to become neurogenic, proliferating progenitor-like cells

Abstract

The purpose of these studies is to investigate how Sphingosine-1-phosphate (S1P) signaling regulates glial phenotype, dedifferentiation of Müller glia (MG), reprogramming into proliferating MG-derived progenitor cells (MGPCs), and neuronal differentiation of the progeny of MGPCs in the chick retina. We found that S1P-related genes are highly expressed by retinal neurons and glia, and levels of expression were dynamically regulated following retinal damage. Drug treatments that activate S1P receptor 1 (S1PR1) or increase levels of S1P suppressed the formation of MGPCs. Conversely, treatments that inhibit S1PR1 or decrease levels of S1P stimulated the formation of MGPCs. Inhibition of S1P receptors or S1P synthesis significantly enhanced the neuronal differentiation of the progeny of MGPCs. We report that S1P-related gene expression in MG is modulated by microglia and inhibition of S1P receptors or S1P synthesis partially rescues the loss of MGPC formation in damaged retinas missing microglia. Finally, we show that TGFβ/Smad3 signaling in the resting retina maintains S1PR1 expression in MG. We conclude that the S1P signaling is dynamically regulated in MG and MGPCs in the chick retina, and activation of S1P signaling depends, in part, on signals produced by reactive microglia.

Figure 1:. Patterns of expression of S1P-related genes.

scRNA-seq was used to identify patterns of expression of S1P-related factors among retinal cells with the data presented in UMAP (a,b,d-h, k,l) or dot plots (c,j,m,n). Aggregate scRNA-seq libraries were generated for cells from (i) control retinas and retinas 24, 48, and 72 h after NMDA-treatment, retinas treated with two or three doses of insulin and FGF2, and retinas treated with insulin, FGF2 and NMDA (a–f). MG were bioinformatically isolated and analyzed from the large aggregate library (g-i), and control retinas and retinas 3, 12 and 48 h after NMDA (k–l). UMAP-ordered cells formed distinct clusters of neuronal cells, resting MG, early activated MG, activated MG and MGPCs based on distinct patterns of gene expression (see methods and Figure S1). UMAP heatmap plots illustrate patterns and levels of expression of S1P receptors S1PR1 and S1PR3 (d), S1P metabolism and transport genes SPHK1, SGPL1, and SPNS2 (e), and ceramide metabolism genes ASAH1, CERS5, and CERS6 (f) illustrate levels and patterns of expression across the whole retina. Dot plots illustrate relative levels of expression (heatmap) and percent expression (dot size) in MG, activated MG, and MGPCs (j,m), whole retina, and other cell types (n) in different UMAP clusters. Significance of difference was determined by using a Wilcox rank sum with Bonferoni correction. Abbreviations: MG, Müller glia; NMDA, N-methyl-D-aspartate; UMAP, uniform manifold approximation and projection.

Figure 2:. Fluorescence in situ hybridization (FISH) for S1PR1, S1PR3 and SPHK1.

Retinas were obtained from undamaged, saline-treated eyes and eyes injected with NMDA at 3, 48 or 72 hrs after treatment. Retinal sections were labeled with antibodies to Sox2 (red) or FISH probes to S1PR1 (green puncta; b), S1PR3 (cyan puncta; d), or SPHK1 (green puncta; f). Dot plots illustrate the average expression (heatmap) and percent expressed (dot size) for S1PR1, S1PR3 and SPHK1 in MG from different treatments and replicates (rep1 or rep2) (a,c,e). Hollow arrows indicate MG nuclei labeled for Sox2 alone and solid arrows indicate MG nuclei associated with numerous FISH puncta. Calibration bars panels b, d and f represent 50 μm. Areas indicated by cyan or yellow are enlarged 2-fold in adjacent panels. Abbreviations: ONL – outer nuclear layer, INL – inner nuclear layer, IPL – inner plexiform layer, GCL – ganglion cell layer, ns – not significant. NMDA, N-methyl-D-aspartate.

Figure 3:. Activation of cell signaling pathways and microglia by S1P.

Retinas were obtained from eyes injected with vehicle or S1P at P7 and P8, and harvested 4hrs after the last injection. Retinal sections were labeled with antibodies to Sox2 (red) and pS6 (green), pStat3 (green), or pERK1/2 (a), or antibodies to CD45 (green) or DNA stain (Draq5; red) (c). Arrows indicate the nuclei of MG. Calibration bars panels a and c represent 50 μm. Areas indicated by cyan or yellow are enlarged 2-fold in adjacent panels. The histograms in b illustrate the mean (bar ± SD) fluorescence intensity sum for pS6, pStat3, pERK1/2 and cFos in MG. The histograms in d illustrate the mean (bar ± SD) numbers of CD45⁺ microglia in the GCL-ONL, at the INL, or whole retina. Each dot represents one biological replicate and blue lines connect replicates from control and treated retinas from one individual. Significance of difference (p-values) was determined by using a paired t-test. Abbreviations: ONL – outer nuclear layer, INL – inner nuclear layer, IPL – inner plexiform layer, GCL – ganglion cell layer, ns – not significant.

Figure 4:. Effects of different inhibitors to S1P receptor, synthesizing enzymes and degrading enzymes on the formation of proliferating MGPCs.

Schematic diagram of the receptors and enzymes involved in S1P-signaling (a). The diagram includes the different agonists and antagonists that were applied to target the different components of the S1P pathway (a). Eyes were injected with vehicle or drug at P6 and P7, NMDA ± drug at P8, EdU ± drug at P9, EdU at P10 and retinas harvested at P11. Retinal sections were labeled for EdU-incorporation and antibodies to Sox2 (b) or Sox9, neurofilament, and phospho-histone H3 (k). Arrows indicate MG, small double arrows indicate presumptive EdU-labeled microglia, and hollow arrow-heads indicate presumptive proliferating NIRG cells. The calibration bar represents 50 μm (b,k). Histograms (c-j, m-p) illustrate the mean (bar ± SD), each dot represents one biological replicate, and blue lines connect counts from control and treated retinas from one individual. Significance of difference (p-values) was determined by using a paired t-test. Abbreviations: ONL – outer nuclear layer, INL – inner nuclear layer, IPL – inner plexiform layer, GCL – ganglion cell layer, ns – not significant, NMDA, N-methyl-D-aspartate, EdU- 5-Ethynyl-2-deoxyuridine.

Figure 5:. Effects of Sphk1 inhibitor on on S1P levels and activation of cell signaling pathways

Eyes were injected with vehicle or PF543 at P7, NMDA ± PF543 at P8, EdU ± PF543 P9, and retinas harvested 4hrs after the last injection. Sections of the retina were labeled for DNA (Draq5; red), ATF3 (green), pS6 (green), Sox2 (red), pSmad1/5/9 (green) (a). Draq5 and Sox2 are included as a partial-field overlay (red, a). Arrows indicate the nuclei of MG. Calibration bar represents 50 μm (a). Eyes were injected with vehicle, PF543, or NMDA at P7, treated with vehicle or PF543 at P8, and harvested 4 hours later (e). LC/MS was used to quantify levels of S1P in retinas treated with saline, NMDA, PF543, NMDA + PF543, or NMDA + SGPL-in-31. Histograms (b-f) illustrate the mean (bar ± SD), each dot represents one biological replicate and blue lines connect counts from control and treated retinas from one individual. Significance of difference (p-values) in b, c, and d was determined by using a paired t-test. Significance of difference (p-values) in e and f was determined using a one-way ANOVA with Šidák correction for multiple comparisons. Abbreviations: ONL – outer nuclear layer, INL – inner nuclear layer, IPL – inner plexiform layer, GCL – ganglion cell layer, ns – not significant, NMDA, N-methyl-D-aspartate, EdU- 5-Ethynyl-2-deoxyuridine.

Figure 6:. Inhibition of S1P synthesis and S1P receptors stimulates neuronal differentiation from MGPCs.

Eyes were injected with NMDA at P7, EdU at P9, vehicle or drug at P10 and P11, and retinas harvested at P12. Retinas sections were labeled for EdU (red; a,c,e,g) and antibodies to HuC/D (green; c), calretinin (green; e) or visinin (green; g). Arrows indicate the nuclei of regenerated neurons. Calibration bars in panels a, c, e and g represent 50 μm. Histograms represent the mean (bar ± SD) and each dot represents one biological replicate (blue lines connect data points from control and treated retinas from the same individual) (b,d,f). Significance of difference (p-values) was determined by using a paired t-test. Abbreviations: ONL – outer nuclear layer, INL – inner nuclear layer, IPL – inner plexiform layer, GCL – ganglion cell layer, ns – not significant, NMDA, N-methyl-D-aspartate, EdU- 5-Ethynyl-2-deoxyuridine.

Figure 7:. scRNA-seq of normal and damaged retinas with and without microglia.

Retinas were obtained from eyes injected with saline or clodronate liposomes at P6, saline or NMDA at P10, and tissues harvested at 24hrs after the last injection. UMAP ordering of cells is displayed for libraries of origin and distinct clusters of cells (a,b). Number of cells in each UMAP cluster and library of origin are listed in panel c. The identities of clusters of cells in UMAP plots were well-established markers (see methods). UMAP heatmap plots illustrate the patterns and levels of expression of S1PR1, S1PR2, S1PR3, SPHK1, SGPL1 and ASAH1 (d-i). Different cell types, including MG (j), amacrine cells (k), bipolar cells (ON and OFF cells; l), and retinal ganglion cells (RGCs; m) were bioinformatically isolated and dot plots were generated to assess levels of expression by treatment (library of origin). Dot plots illustrate the percentage of expressing cells (dot size) and significant (p<0.01) changes in expression levels (heatmap) for genes for cells from retinas treated with saline vs saline-clodronate and NMDA vs NMDA-clodronate. Significance of difference was determined by using a Wilcox rank sum with Bonferoni correction. Abbreviations: NMDA, N-methyl-D-aspartate; UMAP, uniform manifold approximation and projection.

Figure 8:. Rescue of MGPC proliferation in damaged retinas missing microglia

Arrows indicate MG, small double arrows indicate EdU-labeled microglia, and hollow arrow-heads indicate proliferating NIRG cells (a). The calibration bar represents 50 μm. Histograms (b-g) illustrate the mean (bar ± SD), each dot represents one biological replicate and blue lines connect replicates from control and treated retinas from one individual. Significance of difference (p-values) was determined by using a paired t-test. Abbreviations: ONL – outer nuclear layer, INL – inner nuclear layer, IPL – inner plexiform layer, GCL – ganglion cell layer, NIRG- non-astrocytic inner retinal glial cell, ns – not significant, EdU- 5-Ethynyl-2-deoxyuridine.

Figure 9:. Smad3 inhibitor application suppresses S1PR1 transcription

Retinas were obtained from eyes injected with vehicle or SIS3 at P7 and P8, and harvested at P9. Retinal sections were labeled with antibodies to Sox2 (red) and FISH probes to S1PR1 (green puncta; a). Solid arrows indicate MG nuclei associated with numerous FISH puncta. Calibration bar represents 50 μm. Histograms (c-j) illustrate the mean (bar ± SD), each dot represents one biological replicate and blue lines connect replicates from control and treated retinas from one individual. Significance of difference (p-values) was determined by using a paired t-test. Abbreviations: ONL – outer nuclear layer, INL – inner nuclear layer, IPL – inner plexiform layer, GCL – ganglion cell layer, MG- Müller glia, ns – not significant.

Figure 10:. Patterns of expression of S1P-related genes in human and zebrafish retinas.

scRNA-seq was used to identify patterns of expression of S1P-related factors among retinal cells with the data presented in UMAP plots. Aggregate snRNA-seq libraries were generated for human retinas cells at different ages ranging from 42 to 89 years of age (a-e), as described previously (Li et al., 2023). Heatmap plots illustrate the patterns and levels of expression for GLUL and RLBP1 (b), TGFB2 and HBEGF (c), S1PR1, S1PR2 and S1PR3 (d), and SGPL1 and SPHK1 (e) across retinal cells. Aggregate snRNA-seq libraries were generated for zebrafish retinas cells from eyes treated with vehicle, light damage or NMDA and harvest at 1.5, 2.25, 3, 4, 7 and 14 days after damage (f-j), as described previously (Lyu et al., 2023). Heatmap plots illustrate the patterns and levels of expression for glula and rlbp1a (g), ascl1a, top2a, neurod1, and neurod4 (h), s1pr1 and s1pr3 (i), and sgpl1, sphk1 and sphk2 (j) across retinal cells. Abbreviations: NMDA, N-methyl-D-aspartate; UMAP, uniform manifold approximation and projection.