Identification and characterisation of RAs from the postnatal rat cortex
To avoid neural progenitor cell (NPC) contamination, only cerebral cortex regions where NPCs were not abundant were used. Cultured cells were expanded in medium containing TGF-β1, which stimulated astrocytes to become RAs [29]. For about one week, until the cells reached 90% confluence, cytarabine was added to remove proliferating oligodendrocytes and microglia. After purification, most cells were typically star-shaped or polygonal, producing branching processes from the cell bodies (Figure 1A). More than 95% of the cells expressed the RA markers GFAP and S100 without DCX or MAP2 positive neuronal contamination (Figures 1B-1D). This suggested that the cultured cells were mainly RAs.
Characterisation of RAs overexpressing miR-124
To create RAs overexpressing miR-124, we transfected RAs with a miR-124 mimic (miR-124) or a negative control (NC-124). After a 48-h transfection, cell bodies in the miR-124 group became smaller and their branching processes slenderer compared with those in the NC-124 group (Figure 2A). qRT-PCR analysis demonstrated that miR-124 expression increased. The highest transfection efficiency was with 50 nM miR-124 for 48 h (about 23.4 times higher than in the control) (Figures 2B and 2C). As expected, miR-124 overexpression reduced mRNA levels of the RA markers GFAP and S100 (Figure 2D). Additionally, GFAP and S100 protein expression was also downregulated after miR-124 overexpression (Figures 2E and 2F), which was consistent with the mRNA levels. These results suggested that we had successfully constructed miR-124 overexpressing RAs and that miR-124 overexpression suppressed the original astrocytic characteristics of RAs.
RAs synergistically treated with miR-124, ruxolitinib, SB203580, and forskolin acquired neuronal properties
Although miR-124 inhibited RA characteristics, it alone could not convert RAs into neurons. We then treated RAs overexpressing miR-124 with ruxolitinib, SB203580, and forskolin (Figure 3A). This synergistic combination dramatically converted RA morphology into neuron-like shapes after a 7-day treatment. The induced neuronal cell (iN) bodies became smaller and more compact and formed slender axon-like branching structures; no significant morphological changes were observed in the NC-124 control group (Figure 3B). Moreover, miR-124 group cells had smaller territory sizes and longer branching processes (Figures 3C and 3D). Immunostaining demonstrated that early neuronal markers DCX and TUJ1 and mature neuronal markers MAP2 and NEUN were detectable in miR-124 group on the 7th day of induction (Figures 3E-3G); the neurosynaptic junction marker SYN1 was detectable on the 21st day (Figure 3H). Based on neuronal marker expression and cell morphologies, the conversion efficiencies for DCX, TUJ1, MAP2, NEUN, and SYN1 were approximately 41.5%, 64.4%, 38.7%, 87.3%, and 0.7%, respectively (Figure 3M). With regard to neuronal subtypes, our iNs were predominantly CHAT+ (25.9%) cholinergic and VGLUT1+ (22.3%) glutamatergic, occasionally TH+ (6.1%) dopaminergic, but rarely GAD67+ (1.4%) inhibitory neurons on the 7th day of induction (Figures 3I-3L and 3N). These data suggested that rat cortical RAs acquired a neuronal fate after synergistic miR-124 and small molecules treatment, mainly forming cholinergic and glutamatergic neurons.
MiR‑124, ruxolitinib, SB203580, and forskolin were indispensable for RA-to-neuron conversion
To explore whether the same induction effects could be achieved using fewer components, different sets of miR-124 and small molecules were used for induction. The small molecules alone and miR-124 combined with at least one small molecule promoted the generation of DCX-positive neurons at day 7, but conversion efficiency was the highest when miR-124 and the three small molecules were added synergistically (Figures 4A and 4B). MAP2 immunostaining results were also consistent (Figures 4C and 4D). As JAK/STAT3, P38 MAPK, and ERK1/2 signalling is involved in RA migration, astrogliosis, and proliferation [30-32], increased phosphorylation levels would enhance astrocyte activation. Therefore, we assessed the phosphorylation levels in the iNs. Western blot results at day 7 of induction revealed that phosphorylation levels were most downregulated when miR-124 and the three small molecules were added synergistically (Figures 4E and 4F); this also reduced RA reactivity and benefited neuronal transdifferentiation the most. These results suggested that the synergistic effects of miR-124 and the three small molecules best promoted RA-to-neuron conversion.
Gene expression profile of iNs was closer to that of cortical neurons than of RAs
To examine the differences between RAs, iNs, and rat neurons at the genetic level, we first isolated rat cortical neurons and maintained them in a defined, serum-free NM. At day 10 after purification, the cells exhibited a typical neuronal morphology (Figure 5A). Most cells were immune-positive for TUJ1 and MAP2 (Figure 5B), suggesting that the isolated cells were mainly neurons. As miR-124 expression gradually increases during neuronal differentiation and maturation [33], we assessed miR-124 expression in iNs. At day 7 of induction, total miR-124 expression was 23.9 times higher than that in control group (Figure 5C). Considering that exogenous miR-124 decreases gradually (Figure 2C), the total miR-124 overexpression was mainly caused by increased endogenous miR-124, further indicating RA-to-neuron conversion. To identify the gene expression profiles of RAs, iNs (day 7), and rat neurons (day 10), RNA-seq analyses were performed. The number of sDEGs (TPM ≥ 5 in at least one sample, |FoldChange|> 2, q-value < 0.05) between iNs and RAs (3065) was more than that between iNs and rat neurons (2695) (Figure 5D), indicating fewer differences between the latter two. Consistently, hierarchical cluster analysis indicated that the gene expression profile of iNs was more similar to that of rat neurons than of RAs (Figure 5E), similar to the qRT-PCR results (Figure 5F). These results suggested that RAs silenced their astrocytic transcriptional programs and adopted a neuronal identity during induction.
As the induction efficiency was incomplete, we were interested in whether prolonged induction time or environmental factors would affect the fate of the remaining RAs. Therefore, we cultured the 21-day iNs in AM and NM. On the third day of continuous culturing, no morphological changes were observed. On the 11th day, cell morphology in NM was still unchanged, similar to that of rat neurons. However, cell numbers in AM increased, and the morphology was similar to that of RAs (Additional file 1: Figure S1A), which might be because serum exposure allowed the remaining RAs to proliferate. These results indicated that prolonged exposure and a neuronal environment could not change the fate of non-neuronal transformed RAs. Subsequently, we assessed the gene expression of these cells, RAs, and rat neurons. The genetic pattern of cells in AM resembled that of RAs and that of cells in NM resembled that of rat neurons (Additional file 1: Figure S1B).
Transcriptome analysis of iNs and RAs
To explore the molecular and cellular changes of iNs compared with their initial cells, we further analysed the transcriptome of iNs and RAs. At day 7 of induction, the transcriptome scatterplot analysis revealed that, of the 3065 sDEGs, 1820 were upregulated and 1245 were downregulated (Figure 6A). The expression of RA-enriched genes (e.g., Gfap, Aldh1l1, Vim, Nes, Slit2, Axin2, and Mmp2) was decreased, while that of neuron-enriched genes (e.g., Dcx, Tuj1, Syt1, Kcnc1, Map2, Syn1, and Ascl1) was significantly increased. Additionally, expression of genes involved in brain development or cholinergic neuron differentiation (e.g., Foxg1, Bcl11b, and Myt1l) was also significantly upregulated. Further, qRT-PCR results confirmed expression patterns of the selected genes, which were consistent with the RNA-seq results (Figure 6B). The GO analysis (Biological process, Cellular component, and Molecular function) showed that the upregulated genes were significantly enriched in neuron projection morphogenesis, neuron differentiation, axonogenesis, and synaptic vesicle (Figure 6C), and the downregulated genes were significantly enriched in cell proliferation, cell migration, cell-substrate adhesion, and proteoglycan binding (Figure 6D), suggesting that the astrogliosis and scar formation properties of RAs had weakened, while neuronal characteristics had increased. Additionally, KEGG analysis revealed that upregulated cyclic guanosine monophosphate-dependent protein kinase G and cyclic adenosine monophosphate (cAMP) signalling, soluble NSF attachment protein receptor interactions in vesicular transport, and other upregulated signals were mainly involved in neuron development, axon formation, and neurotransmitter release (Figure 6E), while downregulated MAPK signalling, cell adhesion molecules, PI3K-AKT signalling, and other downregulated signals were related to astrogliosis and glial scarring of RAs (Figure 6F), further pointing to RA-to-neuron conversion. Together, the transcriptional regulatory network of RAs was reprogrammed toward that of neuronal lineage under our induction.
MiR-124 regulated neuronal differentiation of RAs via the SOX9-NFIA-HES1 axis
To clarify the molecular mechanism involved in the regulatory role of miR-124, we focused on the NOTCH downstream effectors SOX9 and NFIA, characteristic RA TFs [34], as well as HES1, whose oscillations are related to neurogenesis [35]. MiR-124 overexpression notably repressed SOX9, NFIA, and HES1 expression (Figures 7A and 7B), indicating the regulatory effect of miR-124 on them. Considering that Sox9 is a target of miR-124 [36], we explored the relationship between SOX9, NFIA, and HES1 in RAs. After transfection with si-Sox9 (Additional file 1: Figure S2A), Sox9 mRNA and protein levels were significantly reduced (Figures 7C-7E). Moreover, NFIA and HES1 expression was also downregulated with reductions in SOX9 levels (Figures 7D and 7E), suggesting that SOX9 regulated NFIA and HES1. Additionally, a physical association between SOX9 and NFIA confirmed by co-immunoprecipitation further indicated that SOX9 directly downregulated NFIA (Figure 7F). Subsequently, we transfected RAs with si-Nfia (Additional file 1: Figure S2B), and Nfia mRNA level and protein expression were remarkably decreased (Figures 7G-7I). Moreover, HES1 also decreased upon reduced levels of NFIA (Figures 7H and 7I), which reduces HES1 expression by binding to the Hes1 promoter region [37], suggesting the direct repression of NFIA on HES1. Altogether, miR-124 could downregulate HES1 expression by targeting the SOX9-NFIA-HES1 axis in RAs. To explore whether this axis was related to neurodifferentiation, we transfected RAs with si-Hes1 (Additional file 1: Figure S2C), which significantly downregulated Hes1 mRNA expression (Figures 7J-7L). Furthermore, the mRNA expression of Ascl1, a gene related to neuronal differentiation, was significantly elevated by Hes1 silencing (Figure 7J), revealing the negative regulation of Hes1 in neuronal differentiation. Simultaneously, we detected the effect of Sox9 or Nfia inhibition on Ascl1. The results showed that Sox9 or Nfia inhibition could also increase Ascl1 mRNA expression (Additional file 1: Figure S2D), but not as much as Hes1 inhibition, further indicating the indispensable role of Hes1 in this regulatory axis. These results suggested that in RAs, miR-124 inhibited NOTCH signalling by targeting the SOX9-NFIA-HES1 axis, thereby promoting RA-to-neuron conversion.
HES1 was involved in the RA-to-neuron conversion co-regulated by miR-124, ruxolitinib, SB203580, and forskolin
To decipher whether the SOX9-NFIA-HES1 axis operated during the induction after the small molecules were added, we first induced the aforementioned si-Hes1 RAs with the three small molecules for 3 days, and compared their HES1 expression with that in RAs, RAs treated with the three small molecules only, and RAs treated with both miR-124 and the three small molecules. Western blotting showed that HES1 expression was significantly suppressed in the miR-124 and small molecules group (Figures 8A and 8B). Immunostaining for DCX demonstrated that conversion efficiency in the miR-124 and small molecules group was better than those in other groups (Figures 8E [a-d] and 8F). Given that the miR-124 and small molecules group had the lowest HES1 expression but the best conversion efficiency; we speculated whether HES1 expression determined differences in conversion efficiency. We then applied DAPT, a NOTCH pathway inhibitor, to further suppress HES1 expression. For better effects, we also prolonged induction time to 7 days. As expected, DAPT application further suppressed HES1 expression in the miR-124 and small molecules group (Figures 8C and 8D), but conversion efficiency was still lower than that in the miR-124 and small molecules group without added DAPT (Figures 8E [e-h] and 8G), demonstrating that apart from suppressing HES1 expression, the combination might also regulate neuronal conversion by other means. These results suggested that HES1 expression was involved in RA-to-neuron conversion co-regulated by the miR-124 and small molecule cocktail, but this was not the only mechanism.