Abstract
Neural stem cells (NSCs) can generate new neurons to repair brain injury and central nervous system disease by promoting neural regeneration. MicroRNAs (miRNAs) involve in neural development, brain damage, and neurological diseases repair. Recent reports show that several miRNAs express in NSCs and are important to neurogenesis. Neurites play a key role in NSC-related neurogenesis. However, the mechanism of NSC neurite generation is rarely studied. We surprisingly noticed that the neurites increased after bone morphogenetic protein (BMP) treatment in rat NSCs. This process was accompanied by the dynamic change of miRNA-29. Then we discovered that miR-29a regulated neural neurites in rat hippocampus NSCs. Overexpression of miR-29a reduced the cell soma area and promoted the neurite outgrowth of NSCs. Cell soma area became small, whereas the number of neurite increased. Moreover, neurite complexity increased dramatically, with more primary and secondary branches after miR-29a overexpression. In addition, miR-29a overexpression still maintained the stemness of NSCs. Besides, we identified that miR-29a can promote the neurite outgrowth by targeting extracellular matrix-related genes like Fibrillin 1 (Fbn1), Follistatin-like 1 (Fstl1), and laminin subunit gamma 2 (Lamc2). These findings may provide a novel role of miR-29a to regulate neurite outgrowth and development of NSCs. We also offered a possible theoretical basis to the migration mechanism of NSCs in brain development and damage repair.
Introduction
Neural stem cells (NSCs) have self-renewing ability and can differentiate into neurons, astrocytes, and oligodendrocytes in mammalian brain [1]. They are important to neuronal plasticity, neurogenesis, and neuronal integration [2,3]. Neurogenesis initiates from NSCs and produces new functional neurons. Therefore, neurogenesis plays an important role both in embryonic neural development and adult brain plasticity [4,5]. In adult brain, neurogenesis is regulated by several key factors and signaling mechanisms [5].
MicroRNAs (miRNAs) regulate gene expression at the post-translational levels by either inhibiting messenger RNA (mRNA) translation or inducing mRNA degradation in many biological processes [6,7]. miRNAs are enriched in central nervous system and play pivotal roles during neurogenesis and brain development [8,9]. Reports show that miRNAs involve in neurogenesis, cell fate determination, and cell migration [10]. For example, miR-9 regulates NSC proliferation and differentiation by targeting nuclear receptor TLX (nuclear receptor subfamily 2 group E member 1, NR2E1) and other transcription factors [11,12]. The interplay between epigenetic machinery and miRNAs associates with several cellular mechanism of NSCs fate determination [13].
A key step of neurogenesis is neuronal maturation, which includes dendritic and axonal growth, synaptogenesis, neuronal and synaptic pruning, and so on [14,15]. miRNAs have been implicated in neuronal maturation [16]. miR-9 contributes to neural cells and regulates NSC proliferation and differentiation. In addition, miR-9 overexpression promotes axonal branching and reduces the axonal outgrowth [17]. miR-21 exhibits neurites outgrowth [18]. miR-431 promotes axonal outgrowth and induces axonal regeneration [19]. Not only axons, miRNAs also play an important role in dendrites. miR-137 regulates the dendritic morphology by targeting a ubiquitin ligase named Mind bomb one (Mib1) [16]. Overexpression of miR-541 represses synapsin-I and reduces the neurite extensions [20]. miR-132 is reported to promote dendritic growth and remodeling [21]. These findings illustrate the important role of miRNAs in neurite outgrowth and branching.
Extracellular matrix (ECM) is a complex network of material such as proteins and polysaccharides, provides structural, adhesive, and biochemical signaling support to cells. The ECM is composed of three molecule types: (1) structural proteins, that is, glycoproteic collagens (Col) and nonglycosylated proteins elastins (Eln); (2) specialized glycoproteins, such as fibrillin (Fbn), fibronectin (FN), and laminin (Lam); and (3) heavily glycosylated proteoglycans (PoG), such as chondroitin, heparan, keratan sulfate proteoglycans, and hyaluronic acid [22]. Signaling from the ECM influences regulation of gene expression [22]. Cell surface receptors transduce signals into cells from ECM, which regulate diverse cellular functions, such as survival, growth, migration, and differentiation [23].
miR-29 family suppress many ECM genes, including collagens, elastin, and fibrillins [24,25]. Family members' miR-29a/29b/29c share the same seed sequence, and they are regulated by various transcriptional regulators and signaling pathways [26]. They have emerged as a key modulator of ECM homeostasis [27]. miR-29b overexpression is associated with aging vasculature [24] and aneurysm development [28,29]. Reduced miR-29b is noted in pulmonary fibrosis by upregulation of ECM genes [30]. miR-29 inhibition enhances elastin expression with elastin haploinsufficiency [31]. miR-29a has been demonstrated the anti-tumor functions as a tumor suppressor in chronic lymphocytic leukemia and lung cancer [32], hepatocellular carcinoma [33], pancreatic cancer cells [34], pediatric acute myeloid leukemia [35], papillary thyroid carcinoma [36], and myeloma [37]. However, miR-29a also promotes tumorigenesis in breast cancer [38], nasopharyngeal carcinoma [27], and oral squamous cell carcinoma [39]. miR-29a also induces apoptosis in a p53-dependent way [40]. Thus, miR-29a might have different functions in different cells.
Moreover, ECM regulation by miR-29a has been implicated in fibrosis of many organs including heart [41,42] and liver [43,44]. Reports point that miR-29a highly expresses in mice cerebrum, and can trigger upregulation axon branching by directly repressing doublecortin (DCX) [45 –47]. miR-29a also increases neurite outgrowth by decreasing PTEN (phosphatase and tensin homolog) expression in PC12 cells [48]. Yes-associated protein (YAP) promotes the neurite outgrowth by targeting the promoter of miR-29a in N2a cells [49].
During the maturation of miRNAs, the pre-miRNA hairpin is cleaved to generate two strands. In most miRNAs, the more stable strand is functional and less stable strand tends to degradation. Only one strand is overly dominant, although in some cases, both strands could be functional [6]. Based on miRBase analysis (
In this study, we observed that neurites of rat NSCs increased after bone morphogenetic protein (BMP) treatment. This process was accompanied by the dynamic change of miRNA-29a. Then we investigated the function of miR-29a in this process. We showed that overexpression of miR-29a promoted the neurite outgrowth of NSCs. Moreover, neurite complexity increased dramatically, with multiple primary and secondary branches. Then we found that miR-29a significantly influenced many ECM-related genes, including the collagen family gene Col3a1 (collagen, type III, alpha 1); the laminin-encoding gene laminin subunit gamma 2 (Lamc2); the fibrillin family gene Fbn1, Sparc (secreted acidic cysteine rich glycoprotein), and Follistatin-like 1 (Fstl1). Similar neurite outgrowth displayed after several typical ECM genes knock down. Together, our results suggested that miR-29a may play an important role in NSC neurite outgrowth and development by targeting ECM.
Materials and Methods
Cell culture
Rat NSC cells were cultured under the following conditions: Dulbecco's modified Eagle's medium (DMEM)/F12 (Gibco, Grand Island, NY) containing 1 × N2 supplement (Gibco), 1 × GlutaMAX (Gibco), 1 × penicillin/streptomycin (Gibco), and 20 ng/mL FGF-2 (HumanZyme, Chicago, IL). When cell density was up to 70%–80%, we prepared to passage. Trypsin (0.05%; Gibco) was used to digest NSCs, and trypsin inhibitor (Gibco) to stop the trypsin. Then cells were seeded at the appropriate density after counting.
Poly-l -ornithine and laminin-coated plates
The plates were coated with poly-
BMP-treated NSC cells
The final concentration of BMP-4 (HumanZyme) was 50 ng/mL. BMP-4 was added to the culture medium after passage. We did not change medium or passage during treatment. Discard BMP-4 according to the required time points.
NSC cells transfection
The miR-29a precursor sequences were cloned into FG12 vector (No. 14884; Addgene) to construct recombinant miR-29a overexpression plasmid (miR-29 OX), with scrambled miRNA control plasmid (Control). After 48 h of BMP-4 treatment, the transfection of plasmids was performed in NSCs using electroporation following the manufacturer's instructions (Pulse Generator CUY21EDIT II; BEX), and the concentration of plasmids was 1 μg/μL. After transfection, NSCs were seeded into six-well plates at a density of 1.5 × 105/well. For reverse transcription-quantitative real time polymerase chain reaction (RT-qPCR) assay and RNA sequencing (RNA-Seq), NSCs were treated by 50 ng/mL BMP-4 for 48 h, then miR-29a-3p mimics and mimics control were transfected into NSCs by Lipofectamine 3000 (Invitrogen, Carlsbad, CA), respectively, according to the instructions. miR-29a-3p mimics were double-stranded small RNA molecules, designed based on the mature sequence of rat miR-29a-3p, synthesized by GenePharma Co. (Shanghai, China). For gene knockdown experiments, specific small interfering RNA (siRNA) and control were transfected into NSCs by Lipofectamine 3000 using the same method. Each set of experiments was repeated at least three times.
RT-qPCR assay
Total RNA was isolated from NSCs using TRIzol reagent (Invitrogen). miR-29 and mRNAs were reversely transcribed into complementary DNA (cDNA) using TaqMan™ Advanced miRNA cDNA Synthesis Kit (Thermo Fisher) or Reverse Transcription System Kit (Takara, Dalian, China). RT-qPCR reactions were conducted using RT-qPCR (SYBR Green) Kit (Tiangen, Beijing, China) on an Applied Biosystems 7500 Real-time PCR System (Applied Biosystems, Foster City, CA). All samples were amplified in technical triplicates. The relative quantifications of miR-29 and interested mRNAs were performed by the 2−ΔΔCt method. U6 and GAPDH were used as endogenous controls to normalize the data of miR-29 and mRNAs, respectively. The primers of interested mRNAs were as follows: Sox2 [SRY (sex determining region Y)-box transcription factor 2] forward primer: 5′-CGGCGGCAACCAGAAGAACAG-3′, reverse primer: 5′-CGCTTGGCCTCGTCGATGAAC-3′; Nes (Nestin) forward primer: 5′-ATGAGGAAGGAGCAGAGTCAGGAG-3′, reverse primer: 5′-CAGCACCTCTCAAGCCATCATCC-3′; Fstl1 forward primer: 5′-GTGGCAGTAATGGCAAGACCTACC-3′, reverse primer: 5′-CGGTTAGCCTGATAGCAGACAACG; Col3a1 forward primer: 5′-GACACGCTGGTGCTCAAGGAC, reverse primer: 5′-GTTCGCCTGAAGGACCTCGTTG; Lamc2 forward primer: 5′-AGACACGCTCAACACATTGGAAGG-3′, reverse primer: 5′-AGGTGGAGGTGGTTGCTCTGTC-3′; Sparc forward primer: 5′-AAGCTCCACCTGGACTACATCGG-3′, reverse primer: 5′-CGCTTCTCGTTCTCGTGGATCTTC-3′; Fbn1 forward primer: 5′-TGCCGCATATCTCCTGACCTCTG-3′, reverse primer: 5′-TAGCCTTCGTCACACTTGCATTCG-3′; Cdk6 (cyclin-dependent kinase 6) forward primer: 5′-TGACAGACATCGACGAGCTAGGC-3′, reverse primer: 5′-GCTGGACGACAGGTGAGAATGC-3′; Mmp24 (matrix metallopeptidase 24) forward primer: 5′-GAGGTAGAGCGGCGGAAGGAG-3′, reverse primer: 5′-CCACAGCCACAGCATTCACAGAG-3′; Col5a3 (Collagen, type V, alpha 3) forward primer: 5′-TCAGGTGACCACAGGCACTCTATC-3′, reverse primer: 5′-TTGATGGTGGCTGCTGTTGTCTG-3′. For gene knockdown experiments, we also used RT-qPCR to quantify silencing efficiency. Methods and primers were as described previously. All primers were synthesized by Sangon Biotech Co. (Shanghai, China).
Immunostaining assay
Cells were fixed with 4% PFA for 15 min at room temperature, and then washed with PBS for 3 × 10 min. After blocking with 3% bovine serum albumin (BSA) solution in PBS containing 0.3% Triton X-100 and 0.05% Tween-20 at room temperature for 1 h, cells were incubated with primary antibodies at 4°C overnight. The dilution buffer of antibodies was the same as blocking solution and the dilution ratios were as follows: Sox2 (1:2,000; Santa Cruz, Dallas, TX), neuroepithelial stem cell protein Nestin (Nes; 1:2,000; Abcam, Cambridge, MA), marker of proliferation Ki67 (1:1,000; Thermo Fisher), DCX (1:500; Santa Cruz), glial fibrillary acidic protein (GFAP, 1:1,000; DakoCytomation, Carpinteria, CA), and microtubule-associated protein 2 (MAP2, 1:1,000; Synaptic Systems, Gottingen, Germany). On the second day, primary antibodies were discarded and cells were washed with PBS containing 0.05% Tween-20 for 3 × 10 min. Thereafter, cells were incubated with second antibodies at room temperature for 2 h. The dilution ratio of fluorescent-labeled second antibodies (Jackson ImmunoResearch Laboratories, West Grove, PA) was 1:1,000. Finally, nuclei were stained with DAPI (1:1,000; Sigma) for 10 min. Then the solution was discarded and cells were washed with PBS containing 0.05% Tween-20 for 3 × 10 min. If coverslips were used, the coverslips were picked up carefully and sealed with mounting media (Thermo Fisher).
RNA sequencing
NSCs were plated in 100-mm plates and treated with 50 ng/mL BMP-4 for 48 h. Then we transfected cells with miR-29a mimics and control, respectively, by Lipofectamine 3000 according to the instructions. After 48-h transfection (up in BMP conditions), RNA was isolated, and then total RNA was treated by mRNA enrichment or rRNA removal. Purified mRNA was fragmented into small pieces, and then first-strand cDNA was generated using random hexamer-primed reverse transcription, followed by a second-strand cDNA synthesis. Afterward, A-Tailing Mix and RNA Index Adapters were added by incubating to end repair. The cDNA fragments obtained from previous step were amplified by PCR, and products were purified. Then the double-stranded PCR products were heated, denatured, and circularized by the splint oligo sequence to get the final library. This single-strand circle DNA was sequenced on BGIseq500 platform (BGI Co.). Sample processing and sequencing were performed by BGI Co. in three separate experiments.
The raw data obtained by sequencing were called raw reads. Low-quality reads was filtered out and the so-called clean reads were then aligned to the reference genome by HISAT (Hierarchical Indexing for Spliced Alignment of Transcripts) Software [50], following by new transcript prediction, single nucleotide polymorphism (SNP) and insertion-deletion (InDel), and differential splicing gene detection. After obtaining a new transcript, the new transcript was added with potential coding proteins to the reference gene sequence to form a complete reference sequence for alignment using Bowtie2 [51], and then calculated the gene and transcript expression level by RSEM [52]. Finally, differentially expressed genes between different samples were detected according to requirements. In addition, cluster analysis and functional enrichment analysis can perform on the differentially expressed genes.
Gene ontology enrichment and pathway analysis
We performed pathway and network analysis of differentially expressed genes in Metascape (
Sholl analysis
Sholl analysis was performed by a plugin for ImageJ. It creates a series of concentric circles around the cell soma, and counts how many times the cell intersects with the circumference of these circles [53,54]. User guide provided by Tiago Ferreira and Wayne Rasband of ImageJ (
Statistical analysis
Images were taken under the inverted fluorescence microscope (ECLIPSE Ti; Nikon) or confocal microscope (TCS SP5 II; Leica). Cells were quantified blind using the ImageJ software. Then, the data were collected from ImageJ and calculated by Excel or SPSS. Data were given as mean ± standard deviation (SD). Error bars were defined as SD. Statistical comparison of data sets was performed by two-tailed Student's t-test or one-way analysis of variance (ANOVA). t-Tests were used in Figs. 5 and 6. ANOVA was used in Figs. 3, 4, and 7. P < 0.05 was considered to indicate significant difference. Charts were organized by Adobe Illustrator CS5.
Results
Identification of rat hippocampus NSCs after BMP treatment
NSCs were derived from adult rat hippocampus SGZ (subgranular zone). NSCs were given as a gift from Dr. Fred H. Gage, Salk Institute. In the presence or absence of BMP, NSCs presented two different morphologies. After 50 ng/mL BMP-4 treatment for 48 h, we were very surprised to find an increase in protrusions from NSC cell body, with more neurites and more interaction between neighboring cells compared with control (Fig. 1A).

Rat hippocampus NSCs after BMP treatment.
To identify these NSC cells with more neurites, we detected them by Sox2 and Nestin, two classic NSC cell markers. More than 90% cells were positive for Sox2 and Nestin (Fig. 1B). This suggested that they were NSCs.
To continuously observe these NSC cells with more neurites, we selected different time points for photographing; 12, 36, 72, and 96 h were arranged after BMP-4 treatment. During the first 12 h, the neurites begin to form and cell–cell connection started to establish. In the following few days, the neurites continued to grow and became more complicated (Fig. 1C, upper). Intercellular connections were also getting closer (Fig. 1C, upper). At 96 h, we discarded the old medium and changed fresh medium without BMP-4. Similarly, we selected a series of time points to observe the morphology of NSCs after BMP discarding. The neurite growth was almost stopped and some of the protrusions began to degrade over time (Fig. 1C, lower).
Dynamic change of miR-29 in NSCs
The growth and disappearance of neurites involved in large-scale ECM synthesis and degradation. Multiple genes and signaling pathways were related to this biological process. miR-29 was a known ECM-associated miRNA, which targeted a large number of ECM-related genes [55]. To investigate whether miR-29 involved in this neurite outgrowth process, we selected the corresponding time points after 50 ng/mL BMP-4 treatment and examined the expression pattern of miR-29 in rat hippocampus NSCs. Our results demonstrated that the expression level of miR-29-3p in NSCs was higher than 5p (Fig. 2). Therefore, we decided to put more attention on 3p and chose miR-29-3p as the representative of miR-29 family members. miR-29-3p and 5p expressions showed a similar trend after BMP treatment. During the first 12 h (B0-B12), the expression of miR-29 decreased sharply, ∼60% down (Fig. 2). Consistent with this process, the protrusions of NSC cells began to form. After the initial decline, miR-29 expression slightly rebound. Later, the tendency began to decline again after reaching another high point in ∼48 h (B48). In 96 h, we discarded the old medium and changed the new medium without BMP (day 0). Then miR-29 expression stopped the downward trend and began to gradually rebound (Fig. 2).

Expression pattern of miR-29 after BMP treatment in rat hippocampus NSCs. Expression pattern of miR-29-3p and 5p during the dynamic process of BMP treatment or BMP withdrawn. Data were shown as mean ± SD. SD, standard deviation.
Overexpression of miR-29a significantly altered the cell soma area of rat NSCs
To further study the effect of miR-29 family to rat NSCs, we constructed pre-miR-29a overexpression plasmid and a control plasmid, then transfected them into BMP-treated NSCs separately by electroporation. On the first day after electroporation (24 h), NSCs began to express green fluorescence in both the overexpression group and the control group; on the third day (72 h), the luminescence intensity of GFP reached a peak; but after day 5, the expression of fluorescence had a decreasing trend because it was transient transfection. Some cells died and cell density also decreased (Fig. 3A). After careful comparison of GFP-positive cells, it was found that cell soma areas in the overexpression group were significantly smaller than those in the control group (Fig. 3A). We calculated the fluorescence pixels of the cell body using ImageJ, and compared the relative size of the cell soma area by using this unit. After counting thousands of cells continuously from the first day to the sixth day after electroporation, results indicated that the cell soma areas in miR-29a overexpression group were significantly reduced (Fig. 3B). In addition, the nucleoli of miR-29a overexpressed NSCs were shrunk, rounder, and smoother (Fig. 3A).

Overexpression of miR-29a altered the cell soma area in rat NSCs.
Effects of miR-29a on neurite outgrowth and complexity of rat NSCs
Overexpression of miR-29a not only reduced the cell soma area, but also increased the number of cell processes. Moreover, neighboring cells were widely linked and connected (Fig. 4A). Neurite morphology was often analyzed by Sholl analysis and by counting the number of neurites and branches [56]. Therefore, we investigated the number and complexity of neurites. GFP+ cells were randomly selected from the second day to the fourth day after electroporation. Then we quantified the number of processes that extended directly from the cell soma, and branches extended from these processes, including primary and secondary branches. The primary and secondary branches represented the complexity of neurites. Processes that extended directly from the cell soma of overexpressed miR-29a NSCs were greater than the control group, and this difference became more obvious over time (day 2, 5.13 vs. 3.62; day 3, 7.11 vs. 4.36; day 4, 8.26 vs. 4.97) (Fig. 4B).

miR-29a overexpression promoted neurite outgrowth and branching of rat NSCs.
Statistical results verified that the complexity of neurites in the overexpression group was much higher. On the second day, the average primary branch of miR-29a overexpressed NSCs was 2.34, and the secondary branch was 1.21, whereas the control group had only 0.54 primary branches and almost had no secondary branches (average 0.044). On the third day, the primary and secondary branches in the overexpressed group were 3.15 and 1.78, respectively, which were significantly higher than the previous day. Although the control group was also increased, the average number was still very few, with only 0.85 primary branches and 0.21 secondary branches per cell on average. On the fourth day after overexpression of miR-29a, the average number of branches per cell was almost double compared with the second day: 4.52 primary branches and 2.81 secondary branches, whereas the control group had 1.41 primary branches and 0.19 secondary branches (Fig. 4C).
Sholl analysis results also showed the comparison of complexity in another perspective. Both the sum and mean intersections that overexpressed miR-29a cells were more than that of the control. On day 2, the sum intersections between neurites of miR-29a overexpressed NSCs and the concentric circles were 30.1, the mean intersections of each circumference were 5.39; whereas the sum intersections were 15 and mean intersections were 3.59 in the control group. On day 3, the sum intersections of overexpression group were 36.6 and in the control group it was 16.7, whereas the mean intersections were 10.09 and 4.45 in the two groups, respectively. The same situation also appeared on the fourth day: the overexpression group had 38.7 sum intersections and 9.97 mean intersections, which were more than that of the control group (17.3 sum intersections and 4.84 mean intersections) (Fig. 4D). Taken together, we demonstrated that miR-29a overexpression significantly increased neurite outgrowth and complexity.
miR-29a overexpression maintains the characteristics of NSCs
To confirm the characteristics of transfected cells, we detected the expression of Sox2 and Nestin. We found that both Sox2 and Nestin were expressed in NSC cells that transfected with either the miR-29a overexpression plasmid (Fig. 5A) or the control plasmid (Fig. 5B). It also demonstrated that overexpression of miR-29a did not alter the stem cell characteristics of NSCs.

Stem cell characteristics identification and proliferation analysis to NSCs. Cells transfected with control plasmid
Next, we selected several morphological markers to further confirm the effects of miR-29a overexpression, including DCX, GFAP, and MAP2. Both miR-29a overexpression and control NSC cells were DCX positive, and the expression pattern of DCX was consistent with GFP (Supplementary Fig. S1A). However, few GFP+ cells showed GFAP positive in the two groups (Supplementary Fig. S1B). This result was also confirmed by Confocal (Supplementary Fig. S1C). In addition, MAP2 was hardly expressed in NSCs most probably because it is a mature neuron dendritic marker (data were not shown).
Then, we tried to find out whether miR-29a overexpression influences the proliferation ability of NSCs. Cells were taken for immunostaining from the second day to the fourth day after transfection, respectively. Among them, cell proliferation marker Ki67 labeled proliferating cells and GFP indicated transfected cells. As a result, Ki67-positive cells were significantly decreased, from 45% on the second day to <20% on the fourth day (Fig. 5C), suggesting that NSCs were less active over culturing time. Moreover, the proportion of proliferating cells in the two groups showed no significant difference, indicating that overexpression of miR-29a did not affect the proliferation ability of NSCs (Fig. 5D).
The expression level of miR-29a target ECM-related genes
miR-29 targeted large number of ECM protein-encoding genes, including the collagen family, elastin, laminin, fibrillin, and so on, or genes related to ECM protein synthesis and degradation. Therefore, we selected several typical ECM genes to detect their expression level under miR-29 overexpression. We also tested Sox2 and Nestin and the cell cycle-related gene Cdk6. Because the efficiency of electrotransfection was ∼30%–40%, and low efficiency produced some interference with the quantification of RNA, we tried to use small molecule mimics to achieve the purpose of overexpression. As mentioned previously, miR-29a-3p was more abundant in rat NSCs, so we synthesized miR-29a-3p small molecule mimics and transfected them into BMP-treated NSCs by Lipofectamine 3000 for RT-qPCR assay.
It is very interesting that the expression level of ECM-related genes, which miR-29 targeted, displayed different changes. Overexpression of miR-29a-3p significantly suppressed those ECM-related genes, including the collagen family gene Col3a1 (Fig. 6A), the laminin-encoding gene Lamc2 (Fig. 6B), the fibrillin family gene Fbn1 (Fig. 6C), and Sparc (Fig. 6D), which encoded protein involved in ECM synthesis. Fstl1 (Fig. 6E), another interested gene also showed significant downregulation. But not all genes were downregulated. Col5a3 (Fig. 6F) showed no change in expression. Mmp24 was significantly increased (Fig. 6G). However, there was no significant difference in Sox2 and Nestin between the miR-29a-3p overexpression and control (Fig. 6H). Moreover, Cdk6, an important regulator of cell cycle also showed no significant change after miR-29a-3p overexpression (Fig. 6I). These results indicated that miR-29a overexpression maintained the characteristics and proliferative capacity of NSC cells, which is consistent with the results given in Fig. 5.

Relative expression level of ECM genes in miR-29a overexpression NSCs. After 48 h of miR-29a overexpression, Col3a1
Knockdown experiments of several typical miR-29a target ECM-related genes
miR-29a overexpression altered the neurite outgrowth of NSCs, so we wanted to investigate whether downregulating certain genes may affect the neurite outgrowth. We selected three typical ECM-related genes—Fbn1, Fstl1, Lamc2. They were downregulated by miR-29 overexpression. Then we designed three specific siRNA for each of them to observe the effect of specific gene to NSC neurites. For morphological analysis, we stained DCX and GFAP as before.
First of all, we detected the efficiency of siRNA silencing (Fig. 7A). siRNA had different silencing efficiency compared with the control. For each gene, at least two siRNAs achieved >80% silencing efficiency (Fig. 7A). We could see that DCX-positive cells were increased after siRNA transfection compared with the control (Fig. 7B). As shown previously, there were few GFAP-positive cells in the control group. However, after siRNA transfection, there were varying increases of GFAP-positive cells and cell processes (Fig. 7C). The statistical results indicated that the most significant change was Fbn1-siRNA group, with a higher proportion of GFAP-positive cells in NSCs. In addition, proportion of GFAP-positive cells had also enhanced both in the Fstl1-siRNA and Lamc2-siRNA groups compared with the control (Fig. 7D). In general, knockdown of certain ECM-related genes, which miR-29a targeted, might also achieve similar results with miR-29a overexpression in different degrees.

Knockdown experiments of several typical miR-29a target ECM-related genes.
RNA-Seq analysis of miR-29a regulated mRNAs in BMP-treated rat NSCs
To identify potential mRNA targets of miR-29a in rat NSC cells, cells were transfected with miR-29a-3p mimics and control, followed by RNA-Seq of total RNA. Gene expression Venn diagram is given in Fig. 8A. We could see that 17,336 genes expressed in both groups, and 625 genes specifically expressed in miR-29a overexpression group, whereas only 754 genes expressed in the control group (Fig. 8A). FPKM values of the differential genes of each comparison group were clustered, and the heat map is given in Fig. 8B. We calculated the FPKM values ratio of the control (numerator) to the miR-29a OX group (denominator), and filtered out genes with a value not <1.5, it can be considered as a specific downregulation. Conversely, if the ratio of the miR-29a OX (numerator) to the control (denominator) was not <1.5, it could be considered as a specific upregulation.

KEGG Pathway and GO Biological Process Enrichment Analysis of RNA-Seq data.
We analyzed the transcriptome data from RNA-seq to identify network processes, pathway maps, and GO (processes) between miR-29a versus control data sets. The top 20 GO cellular processes downregulated or upregulated are given in Fig. 8C and D. Between two images, the gene clustering with specific downregulation was more interesting (Fig. 8C). Cell adhesion-related processes were downregulated such as cell–substrate adhesion, regulation of response to external stimulus, and positive regulation of cell adhesion (Fig. 8C). Branch morphology-related processes were upregulated like morphogenesis of a branching epithelium and reproductive structure development (Fig. 8D). These processes involve a large number of ECM synthesis and regulation. They may widely take part in intercellular communication, cell junction, cell migration, and intercellular signal transduction.
Discussion
In this study, we first identified the dynamic change of miR-29a after BMP-4 treatment and withdrew BMP-4 in rat NSCs. Then we found that there was a negative correlation between miR-29a expression and NSCs neurite outgrowth. We found that miR-29a overexpression increased the neurite number of NSCs significantly. Moreover, the complexity of neurites also increased, with multiple primary and secondary branches. Cell soma area of miR-29a overexpressed NSCs was dramatically reduced. Meanwhile, miR-29a overexpression did not alter the stem cell characteristics of NSCs. Finally, we also identified several ECM-related genes that downregulated under miR-29a-3p overexpression. Then siRNA were designed targeting three typical genes, Fbn1, Fstl1, and Lamc2. Silencing any of Fbn1, Fstl1, and Lamc2 could promote neurite outgrowth of NSCs.
Neurogenesis is dynamically regulated by a number of extracellular signals including BMPs [57]. BMPs can act as short-range morphogens and profoundly affect adult NSC proliferation and differentiation [5,58]. BMP negatively regulates neurogenesis by promoting NSCs astroglial fate commitment and quiescence [3,59,60]. Certain levels of BMP-4 combination with FGF-2 reversibly induced the quiescent state of NSCs [3]. We introduced BMP-4 into NSCs, and were surprised to find the different morphology of neurites. We also found the role of miR-29 in neurite outgrowth.
TGF-β/BMP signaling is involved in a vast majority of cellular processes and is fundamentally important throughout life. In addition, the crosstalk have been recognized between BMP/TGF-β signaling and many other major signaling pathways such as MAPK, Wnt, Hedgehog, Notch, miR-29, and FGF [61,62]. TGF-β1 is the effective target genes of miR-29 and meanwhile, miR-29b decreased the expression of TGF-β1 [62]. These reports reveal the complex feedback between BMP and miR-29. However, the regulation network between them and the function in neurites still need further investigation. In this study, we only focused on the role of miR-29a in neurite outgrowth within the presence of BMP.
miR-29 family is associated with tumorigenesis, fibrosis, and ECM-related protein regulation [26,63]. Reports have shown that miR-29 family is upregulated in the adult cortex and striatum [64]. These findings indicate that miR-29 may play an unknown role in neural development and neurogenesis. In 2014, Shin et al. demonstrated that miR-29b plays a pivotal role in fetal mouse neurogenesis by regulating ICAT-mediated Wnt/β-catenin signaling [65]. In 2018, Shi et al. showed that miR-29a promotes the neuronal differentiation and decreases the astrocytes differentiation of rat embryonic NSCs by targeting PTEN [66]. Our data show that overexpression of miR-29a neither changed the stem cell characteristics of adult rat hippocampus NSCs, nor affected its proliferation ability. These data suggested that miR-29a may have the different function as miR-29b in NSCs. In addition, the effects of miR-29a in adult and embryonic NSCs need more attention.
As a famous tumorigenesis regulator, the predicted target gene of miR-29 contains a large number of ECM-associated genes. miR-29a upregulation triggers axon branching by directly repressing DCX [46]. Zou et al. found that miR-29a can promote neurite outgrowth of PC12 cells by decreasing PTEN [48]. Another article also reported similar result. YAP-regulated miR-29a promotes the neurite outgrowth by decreasing PTEN in N2a cells [49]. These findings indicated that miR-29a is an important regulator of neurite outgrowth, and also suggest that miR-29a may be functional in NSCs and neurogenesis.
By analyzing the strand that is complementary to the RNA of interest, we focused on miR-29a-3p, which potentially targeted a large number of ECM genes. Our results also confirmed that miR-29a-3p is more abundant than 5p in rat NSCs. We demonstrated that miR-29a overexpression significantly reduced the cell soma area and increased neurite outgrowth. Furthermore, neurites became more complex with multiple primary and secondary branches. RNA-Seq analysis of miR-29a-regulated mRNAs verified that the specific downregulation gene clusters are mostly related to cell adhesion, cell–substrate interaction, and extracellular structure. These genes may widely take part in intercellular communication, cell junction, cell migration, and intercellular signal transduction. This novel finding may be of benefit to study the interaction of NSCs and niche during development.
At the same time, we tried to knockdown miR-29a by using specific inhibitor (synthesized by GenePharma Co.). Previous reports point that brain-specific knockdown of miR-29 results in neuronal cell death in mice [67]. However, we had not yet observed the visible differences in rat NSCs (data not shown). One possible reason is that the expression level of miR-29a has already been downregulated in BMP-treated NSC cells (as shown in Fig. 2), so the inhibitor we used may not achieve any effect. Anyway, it needs further investigation.
During brain development and neurogenesis, the migration of newborn neural cells plays a pivotal role. Migration is controlled by a complex dynamic molecular regulatory network. Our findings provided a new sight that miRNAs might affect neural cells migration by regulating neurite outgrowth and branching. In addition, the function of miR-29 in radial glial cells was also worth noting. As famous “scaffold” that guide neuronal migration, radial glial cells also serve as progenitors of neurons and glia [68 –70]. Radial glial cells also help to direct axonal and dendritic process outgrowth, and then regulate synaptic development and function [68 –70]. Does miR-29 function in radial glial cells? Whether it regulates the process of glial cells? What is the role of miR-29 in cell migration regulatory network? In future study, we will pay more attention to these questions and try to explore the intrinsic mechanism of miR-29 function.
Conclusions
During brain development and neurogenesis, the migration of new born neuronal cells plays an essential role and is critical to neocortex formation. The process of cell migration is also important to brain injury repair. These biological processes are partly controlled by neurite outgrowth and branching in NSCs. However, the regulatory mechanism is not sufficient and the functional study of miR-29 in NSC is limited. In this study, we investigated the function of miR-29a in BMP-treated NSCs. First, we demonstrated that miR-29a reduced cell soma area of rat NSCs and promoted neurite outgrowth and branching. Then, we found that NSCs of overexpressed miR-29a still maintain the stem cell characters. Finally, we identified that miR-29a regulated the neurite outgrowth of NSCs by downregulating the ECM-related genes. This novel finding may provide a new insight for studying the interaction of NSCs and niche during migration, and also offered a possible theoretical basis to the migration mechanism of NSCs in brain development and brain damage repair.
Footnotes
Author Disclosure Statement
No competing financial interests exist.
Funding Information
The study was supported by National Natural Science Foundation of China grant (81901031), Natural Science Foundation of Shanghai (19ZR1445400), General Program Shanghai Municipal Health and Family Planning Commission grant (201740091), Natural Science Foundation of Shanghai (16ZR1430600), National Key R&D Program of China (2017YFA0104100), National Natural Science Foundation of China grants (31571058, 31471029) and the National Major Scientific and Technological Special Project for Significant New Drugs Development (2018ZX09201002-005).
Supplementary Material
Supplementary Figure S1
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
