Fig 2: Two distinct clusters in ventral lumbar spinal motoneurons post-SNC(A) Overview of the investigated timepoints and workflow of the experiment.(B) Pearson’s correlations among samples captured based on expressed genes (gene with average FPKM >1 was kept).(C) Gene number identified per condition. The number above the boxplot represented the sample number.(D) Histogram of AUC score of neuronal markers reflecting the high activity of neuronal characteristics.(E) Principal component analysis (PCA) visualization of samples distribution.(F) UMAP clustering of samples. Two distinct clusters were shown.(G) The sample proportion of each group from different injury timepoints in two clusters.(H) Expression profiles of top 20 differential expression genes in either cluster 1 or cluster 2. (I) Immunostaining of ANXA3 (red) in rat ventral lumbar spinal cord at different injured timepoints post-SNC. Neuron and nuclear were labeled with NeuN (green) and DAPI (blue). Scale bar: 100 µm.

Fig 3: Immunohistochemistry and positive rate statistics. (A1-D2) show the specific expression of ANXA3 with SORL1 in AS and controls. (E, F) shows the positive rate statistics of all immunohistochemical images of SORL1 in AS and control groups. “*” indicates p < 0.05, “***” indicates p < 0.001.

Fig 4: Immune correlation lollipop plot. (A) shows the immune cells with significant difference of ANXA3 in AS and control. (B) shows the immune cells with significant difference of SORL1 in AS and control.

Fig 5: Drug sensitivity analysis. The graph shows the drug sensitivity analysis of ANXA3 and SORL1 with various drugs. If Cor > 0, it means the gene is positively correlated with drug sensitivity, the higher the gene expression, the stronger the sensitivity of this drug; conversely, Cor < 0 means the gene is negatively correlated with drug sensitivity, the higher the gene expression, the worse the sensitivity of this drug.