- Research article
- Open Access
Comparison of gene expression profile of the spinal cord of sprouting-capable neonatal and sprouting-incapable adult mice
BMC Genomics volume 20, Article number: 619 (2019)
The regenerative ability of severed axons in the central nervous system is limited in mammals. However, after central nervous system injury, neural function is partially recovered by the formation of a compensatory neural circuit. In a mouse pyramidotomy model, axonal sprouting of the intact side of the corticospinal tract is observed in the spinal cord, and the axons make new synapses with the denervated side of propriospinal neurons. Moreover, this sprouting ability is enhanced in neonatal mice compared to that in adult mice. Myelin-associated molecules in the spinal cord or intrinsic factors in corticospinal neurons have been investigated in previous studies, but the factors that determine elevated sprouting ability in neonatal mice are not fully understood. Further, in the early phase after pyramidotomy, glial responses are observed in the spinal cord. To elucidate the basal difference in the spinal cord, we compared gene expression profiles of entire C4–7 cervical cord tissues between neonatal (injured at postnatal day 7) and adult (injured at 8 weeks of age) mice by RNA-sequencing. We also tried to identify discordant gene expression changes that might inhibit axonal sprouting in adult mice at the early phase (3 days) after pyramidotomy.
A comparison of neonatal and adult sham groups revealed remarkable basal differences in the spinal cord, such as active neural circuit formation, cell proliferation, the development of myelination, and an immature immune system in neonatal mice compared to that observed in adult mice. Some inflammation-related genes were selectively expressed in adult mice after pyramidotomy, implying the possibility that these genes might be related to the low sprouting ability in adult mice.
This study provides useful information regarding the basal difference between neonatal and adult spinal cords and the possible differential response after pyramidotomy, both of which are necessary to understand why sprouting ability is increased in neonatal mice compared to that in adult mice.
Central nervous system injury is devastating as it frequently induces persistent neurological deficits and there is still no effective therapy to restore function besides rehabilitation. Especially, injury to the corticospinal tract (CST) causes motor deficits, because it is the main pathway involved in voluntary movement in humans . The cell body of corticospinal neurons (CSNs) localize to the cerebral cortex and their axons pass through the ventral surface of the medulla oblongata (pyramid). The main CST axon crosses to the contralateral side and then passes through the lateral column, forming synapses with neurons in the spinal cord . In rodents, although the main CST passes through the dorsal column of the spinal cord instead of the lateral column, the CST is essentially conserved, and it is important for precise movement . For these reasons, mouse or rat CST injury models such as traumatic brain injury  or pyramidotomy  are often used to study central nervous system injury.
After central nervous system injury, truncated axons rarely regenerate in mammals. However, it is known that axonal branches from surviving neurons are newly formed and that a compensatory neural circuit is established, leading to partial functional recovery. For example, the axonal branch of the intact side of the CST grows into the denervated side of the cervical cord and forms synapses with the propriospinal neuron approximately 4 weeks after brain injury in mice .
Intriguingly, the sprouting ability of the CST depends on the age of the mice. In a mouse brain injury model, axonal sprouting gradually decreases from postnatal day 3 (P3) to P21 . In a mouse pyramidotomy model, the sprouting ability of P7 mice is higher than that of 2-month-old mouse . Previously, the absence of myelin was thought to be the cause of the enhanced sprouting ability observed in neonatal mice. Components in myelin called myelin-associated molecules are known to inhibit axonal growth . In neonatal mice, myelination is not completed, and a gradual decrease in sprouting correlates with myelination . However, the deletion of paired immunoglobulin-like receptor B, which is a high affinity receptor for myelin-associated molecules, does not enhance sprouting after brain injury in adult mice , suggesting that there are other factors that determine sprouting ability.
Intrinsic factors represent another explanation for the high sprouting ability in neonatal mice. For example, activity of the mechanistic target of rapamycin kinase pathway in CSN is increased in neonatal mice compared to that in adult mice, and deletion of phosphatase and tensin homolog, which is a negative regulator of the mechanistic target of rapamycin kinase pathway, enhances sprouting in adult mice . However, sprouting ability in mutant mice is not as high as that in neonatal mice. In addition to it, pharmacological inhibition of the mechanistic target of rapamycin kinase pathway does not affect sprouting in neonatal mice after pyramidotomy , suggesting that there are other factors that determine this ability.
The glial response has both beneficial (such as neurite outgrowth-enhancing) and detrimental (such as axon growth-inhibiting) effects after central nervous system injury . For example, expression of astroglial and inflammatory markers at the spinal cord lesion site are higher in old mice compared to young mice, and deletion of phosphatase and tensin homolog enhances axon regeneration distal to injury only in young mice , suggesting that altered glial response or inflammatory response affect age dependent decline of axon regenerative ability. In pyramidotomy model, microglial activation in the denervated side of the dorsal column of the spinal cord was reported 4 days after pyramidotomy in mice . Therefore, it is possible that the spinal cord response in the early phase after CST injury is important for axonal sprouting, and the difference in this response might be another factor that determines the enhanced sprouting ability in neonatal mice.
Here, to elucidate the basal difference of the spinal cord between neonatal and adult mice, we compared gene expression profiles by RNA-sequencing (RNA-seq) using entire cervical cord tissues at C4-C7 levels. We also compared gene expression profile at the initial phase after pyramidotomy in an effort to identify discordant gene expression changes that inhibit axonal sprouting in adult mice. In a basal state, genes related to neurite circuit formation, cell proliferation, or the development of myelination were found to be upregulated in neonatal mice compared to expression in adult mice. A number of responses after pyramidotomy, such as the inflammatory response and axonal outgrowth, were common between neonatal and adult mice, but some genes related to the inflammatory response were selectively upregulated in adult mice. Although many transcriptomic analyses have been performed after CST injury in previous studies [12,13,14,15,16,17], none of them compared gene expression profiles between neonatal and adult spinal cords. Therefore, this is the first study to provide gene expression profiles of neonatal and adult spinal cords after pyramidotomy, which is necessary to understand why sprouting ability is enhanced in neonates.
First, to compare the sprouting ability between neonatal and adult mice after CST injury, we injured the left pyramid of the medulla oblongata (pyramidotomy or py) of postnatal day 7 (P7) neonatal mice or 8-week-old (8 W) adult mice. When we performed immunohistochemistry for protein kinase C gamma (PKCγ), which is a CST marker, the signal in the denervated side of the CST area of the cervical cord almost completely disappeared 4 weeks after injury in both neonatal and adult mice (Fig. 1a-d), indicating that the CST was successfully injured. Then, 2 weeks after injury, we injected anterograde tracer biotinylated or tetramethylrhodamine conjugated dextran amine into the intact side of the forelimb area of the motor cortex to label the axon of CSNs. Four weeks after injury, the number of tracer-positive axons sprouting into the denervated side of gray matter in the cervical cord was significantly higher in neonatal mice compared to that in adult mice (Fig. 1e-k), which is consistent with the results of a previous study .
Previously, it is reported that the microglial response in the spinal cord is observed 4 days after pyramidotomy in mice . To detect early changes in the gene expression profile, we assessed the cervical cord 1 day earlier than that report (i.e. 3 days after injury (dpi)), for RNA-seq (Fig. 2a). Although single cell analysis provides useful information, cell dissociation procedure might cause artifacts. Therefore, whole tissue RNA-seq using RNA directly extracted from fresh tissue is a good choice to grasp the picture of a phenomenon as a first step. The medulla oblongata was also harvested when we collected the cervical cord and the lesion was evaluated; only animals in which a lesion was successfully generated were selected for RNA-seq analysis (Fig. 2b, c). In neonatal mice, the lesion was confirmed by a complete loss of PKCγ signal in the left CST (Fig. 2b). In adult mice, PKCγ did not disappear at 3 dpi at the lesion site (Fig. 2c). Therefore, we confirmed the successful generation of a lesion in adult mice based on the abnormally high intensity of the PKCγ signal and disrupted CST structure (Fig. 2c).
Next, 13–18 million 75-bp single-end reads were produced for each sample with three biological replicates for each group (P7 sham, P7 pyramidotomy, 8 W sham, and 8 W pyramidotomy). First, we performed cluster analysis based on a distance matrix to examine the similarity between samples. For this, each replicate in the same group clustered together (Fig. 2d) supporting the validity of RNA-seq. Moreover, samples from animals of the same age, rather than the same surgical condition, clustered together, suggesting that age has a greater effect on gene expression than surgical procedure.
Then, we analyzed the base change between neonatal sham and adult sham groups. For this, 4637 genes were significantly upregulated in the P7 sham group compared to expression in the 8 W sham group, and 4616 genes were significantly downregulated (Fig. 3a). We chose some of the significantly up/downregulated genes and validated the results of RNA-seq by quantitative reverse transcription-polymerase chain reaction (qRT-PCR; Fig. 3b–e; Additional file 2B–W). Results of qRT-PCR were similar to those of RNA-seq, strongly supporting the validity of RNA-seq. Gene enrichment analyses revealed that several Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways such as cell cycle or axon guidance were significantly upregulated in neonatal mice (Fig. 3f, Additional file 4), whereas pathways including Toll-like receptor signaling were significantly downregulated (Fig. 3g). These results indicated a developmental change between the neonatal and adult spinal cord.
We also tried to compare gene expression profiles between neonatal pyramidotomy and neonatal sham groups (Fig. 4a). Although the number of genes that were downregulated in the neonatal pyramidotomy group compared to expression in the sham group (731) was higher than the number of genes that were upregulated (266), no gene was downregulated by more than 2-fold, whereas some genes showed more than 2-fold upregulation. When we compared gene expression profiles between adult pyramidotomy and adult sham groups, a similar tendency was observed (Fig. 4b). We considered that a small fold-change might be unreliable and thus focused on pathways that were activated in the pyramidotomy group. Most of the pathways enriched in the neonatal pyramidotomy group compared to expression in the neonatal sham group (Fig. 4c) were also enriched in the adult pyramidotomy group compared to that in the adult sham group (Fig. 4d). In fact, all top 10 pathways enriched in the P7 pyramidotomy group were also significantly enriched in the 8 W pyramidotomy group (Fig. 4c, d, Additional files 5 and 6), and top 10 pathways enriched in the 8 W pyramidotomy group were also significantly enriched in the P7 pyramidotomy group (Fig. 4c, d, Additional files 5 and 6), except for the cell cycle pathway (Fig. 4d). These results suggest that the response induced by pyramidotomy is very similar between neonatal and adult mice.
In an attempt to reveal any minor differences between neonatal and adult mice after pyramidotomy, we considered both surgery and age. To identify genes that were selectively up/downregulated in neonates after pyramidotomy, we explored genes that were significantly up/downregulated in the P7 pyramidotomy group compared to expression in the P7 sham and 8 W pyramidotomy groups. Both the number of differentially expressed genes (Fig. 5a, b) and the maximum fold-change (Fig. 5e, f) were smaller for the comparison of surgery groups (P7 pyramidotomy VS P7 sham) than for the comparison of age groups (P7 pyramidotomy VS 8 W pyramidotomy), again suggesting that differences in age have a greater effect on gene expression profiles than differences in surgical procedure. When we selected some of the genes that were significantly up/downregulated based on both surgery and age comparisons (hereafter referred to as P7 Py UP/DOWN genes; red dots in Fig. 5e, f), and examined the results by qRT-PCR (Additional file 2X-AA or AB-AG), we could not detect a significant difference, possibly due to too small of differences (less than 2-fold change based on RNA-seq data). Similarly, when we explored genes significantly up/downregulated in the 8 W pyramidotomy group compared to expression in the 8 W sham and P7 pyramidotomy groups, both the number of differentially expressed genes (Fig. 5c, d) and the maximum fold-change (Fig. 5g, h) were smaller for the comparison of surgery groups (8 W pyramidotomy VS 8 W sham) than for the comparison of age groups (8 W pyramidotomy VS P7 pyramidotomy), and we could not detect significant differences by qRT-PCR for genes that were significantly downregulated based on both surgery and age comparisons (hereafter referred to as 8 W Py DOWN genes; Additional file 2AH, AI). However, the genes significantly upregulated for both surgery and age comparisons (hereafter referred to as 8 W Py UP genes) showed relatively high fold-changes. Among them, two genes (Chemokine (C-C motif) ligand 6 (Ccl6) and Cd52 antigen (Cd52)) showed more than a 2-fold upregulation in the 8 W pyramidotomy groups compared to expression in the 8 W sham and P7 pyramidotomy groups. We then validated 8 W Py UP group selective expression by qRT-PCR (Fig. 5i–p) and performed gene enrichment analysis using the 8 W Py UP list (Fig. 5q). We found that KEGG pathways related to an inflammatory response (Staphylococcus aureus infection, pertussis, lysosome, and complement and coagulation cascade) were significantly enriched in this group. These results suggest that a specific inflammatory response is much more severe in adult mice than in neonatal mice after pyramidotomy.
To support this observation, we examined microglial/macrophage activation by immunohistochemistry for Ionized calcium binding adaptor molecule 1 (Iba1) at 3 dpi. In adult mice, Iba1+ cells accumulated at the denervated side of CST after pyramidotomy, but in neonatal mice, Iba1+ cells were not different between intact and injured mice (Fig. 6), suggesting that microglia/macrophage are activated at the denervated side of CST after pyramidotomy only in adult mice. To connect this observation to gene expression profile, we performed double staining of in situ hybridization for 8 W Py UP genes and immunohistochemistry for Iba1. We detected signals for Lymphocyte antigen 86 (Ly86) and macrophage expressed gene 1 (Mpeg1), although signal intensity between neonatal and adult groups are not comparable because sensitivity is different between neonatal and adult samples. For both of the genes, in situ hybridization signal was detected at the denervated side of CST in adult mice after pyramidotomy, whereas it was not different between intact and pyramidotomy group in neonatal mice (Figs. 7, 8 and Additional file 3), which is consistent to the result of RNA-seq and qRT-PCR (Fig. 5m-p). Iba1 signal was very weak in neonatal group, possibly because proteinase treatment decreased antigenicity of Iba1 in neonatal mice. Most of the signals for Ly86 and Mpeg1 in the denervated side of CST in 8 W Py group colocalized with Iba1 signal (Figs. 7n-p and 8n-p), suggesting that upregulation of immune related genes in 8 W Py group is at least partially explained by microglial/macrophage activation.
In this study, we compared gene expression profiles between neonatal and adult mice to identify the factors that determine enhanced sprouting ability in neonatal mice. The main finding of this study is that the difference in gene expression profiles between neonatal and adult mice was very large (Figs. 2d and 3a). The top 10 pathways that were significantly enriched in neonates could be categorized into four groups as follows: (1) protein digestion and absorption, extracellular molecule (ECM)-receptor interaction, focal adhesion, axon guidance, and Hippo signaling, (2) cell cycle and DNA replication, (3) phosphoinositide-3-kinase (PI3K)–Akt and Wnt signaling, and (4) the spliceosome (Fig. 3f). Obviously, proteins are absorbed in the intestine, and not the spinal cord; probably the reason why the protein absorption pathway was detected is that many collagen-encoding genes are included in this pathway, and elevated expression of these genes led to enrichment of this pathway. Therefore, probably five pathways in the first group essentially reflect the same phenomenon, specifically, active axon elongation and neural circuit construction in neonatal mice. Adhesion acts as a molecular clutch during axonal elongation, and the growth cone moves with the aid of guidance molecules . Axonal outgrowth is controlled by interactions between integrins and the ECM such as collagen, laminin, or tenascin . Hippo signaling is suggested to coordinate dendrite dynamics . Further, it was reported that synapse formation on Renshaw cells of the spinal cord increases until 2 weeks after birth [21, 22]. We collected the cervical cord of neonates at P10; therefore, integration of the neural circuit in the spinal cord might still be in progress at that time.
Cell cycle and DNA replication pathways might reflect active cell proliferation in the neonatal spinal cord. It has been reported that in P8 mice, approximately 10% of cells proliferate in the central canal lining of the spinal cord, whereas in adult mice, this value is less than 1% . These cells are thought to differentiate into glial cells.
PI3K–Akt and Wnt signaling pathways are known to be involved in the myelination or differentiation of oligodendrocyte precursor cells . A previous study from our laboratory reported that myelin can be detected in white matter as early as P8–P14 , suggesting that myelination is active in P10 mice used in the present study. Therefore, these pathways might reflect the differentiation of oligodendrocyte precursor cells and myelination in the neonatal spinal cord.
It was also reported that polypyrimidine tract-binding protein-associated splicing factor is highly expressed in differentiating neonatal neurons . It is possible that the spliceosome is important for the full maturation of neurons in the neonatal spinal cord. In addition to this, in adult mice, the oxidative phosphorylation pathway was activated (Fig. 3g). Oxidative phosphorylation is activated when neural stem cells differentiate into mature neurons . Therefore, these pathways might represent maturation of the neuron.
The other pathways enriched in adult mice were the lysosome, Toll-like receptor signaling, and complement and coagulation cascade pathways (Fig. 3g). The lysosome is an important component of phagocytic pathways in microglia , Toll-like receptors are important pattern recognition receptors in microglia , and the complement cascade is important for the inflammatory response and synapse elimination mediated by microglia . In the early developmental stages, microglia have an amoeboid form and they take on a complete ramified form at P28 . This suggests that the function of these cells is not completely developed in neonatal mice, and the aforementioned pathways might reflect the difference in microglial function.
Interestingly, some differentially expressed genes with prominent fold change have not been analyzed well. For example, among the genes that were significantly upregulated in the adult sham group, the fold-change of Ethanolamine phosphate phospholyase (Etnppl) was the largest. Etnppl encodes a protein that catalyzes the degradation of ethanolamine phosphate [31, 32], and its expression level in the brain was found to be higher in schizophrenia and bipolar patients compared to that in the normal population . This suggests that this gene is important for neural function; however, the function of Etnppl in vivo is not known. Interestingly, the expression level of this gene was quite low in neonates and mildly suppressed after pyramidotomy in adult mice, suggesting that its expression negatively correlates with axonal growth (Additional file 2 B, C). Analyzing this gene might reveal important aspects of neural development and neural circuit repair in future studies.
Second, we found that most differentially expressed genes after pyramidotomy were common between neonatal and adult mice (Fig. 4c, d), and among these, the top 10 pathways, except for the cell cycle pathway, were categorized into two groups as follows: (1) complement and coagulation cascades, nucleotide oligomerization domain (NOD)-like receptor signaling, cytokine-cytokine receptor interaction, osteoclast differentiation, Toll-like receptor signaling, chemokine signaling, phagosome, and nuclear factor (NF)-kappa B signaling, and (2) PI3K-Akt signaling, Janus kinase (Jak)-signal transducer and activator of transcription (STAT) signaling, ECM-receptor interaction, and focal adhesion (Fig. 4c, d). The first group comprises pathways related to the inflammatory response. Complement, pattern recognition receptors such as NOD-like or Toll-like receptors, cytokines, chemokines, and NF kappa B signaling function coordinately during the innate immune response , whereas osteoclast differentiation might reflect microglia/macrophage differentiation, because all of these cells are considered tissue macrophages and share the differentiation pathways . The phagosome represents phagocytosis by microglia/macrophages . Because these cells are known to be activated following axonal injury  and microglia accumulation is observed in the denervated side of the CST , these pathways probably reflect the activation of microglia on the degenerating CST.
PI3K, Jak–STAT signaling, ECM receptor interaction, and focal adhesion are involved in axonal sprouting or axonal growth. The deletion of phosphatase and tensin homolog, which is a negative regulator of PI3K–Akt signaling, enhances sprouting after pyramidotomy . The deletion of suppressor of cytokine signaling 3, which negatively regulates the Jak–STAT pathway, also enhances sprouting after pyramidotomy . These reports indicate that both pathways are important for sprouting. Axonal sprouting of the propriospinal neuron has also been suggested after spinal cord injury . Considering these reports, our results implicate axonal sprouting of the propriospinal neuron in both neonates and adults after pyramidotomy.
In contrast, the cell cycle was upregulated only in adults after pyramidotomy (Fig. 4d). There are two possibilities for this: (1) in neonates, cell proliferation is already active in the basal state (Fig. 3f), which masks the activation of cell proliferation after pyramidotomy, or (2) in neonates, cell proliferation is not activated after pyramidotomy. One possible type of cell that proliferates after pyramidotomy in adult mice is ependymal cells, because after spinal cord injury, ependymal cells were previously found to proliferate . Examining these possibilities and identifying proliferating cells after pyramidotomy in adult mice in future studies will deepen our understanding of the cell response after pyramidotomy and differences between neonates and adults.
Third, we found that pathways related to the inflammatory response showed significantly high expression in the adult pyramidotomy group compared to that in adult sham and neonatal pyramidotomy groups (Fig. 5q). Activation of Staphylococcus aureus infection and pertussis pathways probably does not mean that the mice were infected but rather reflect the inflammatory response. In the future, detailed mechanisms should be revealed not by pathway analyses, but by the analyses of individual candidate genes. Among the genes selectively expressed in the adult pyramidotomy group, only Ccl6 and Cd52 showed more than 2-fold upregulation (Fig. 5g, i-l). Ccl6 is considered the rodent homolog of human CCL23 [40, 41] and its expression is observed in some pathological conditions of the central nervous system. For example, it is expressed in neurons or infiltrating blood cells after ischemia in rats , or microglia/macrophages in a mouse model of experimental autoimmune encephalomyelitis ; moreover, CCL23 is expressed in infiltrating blood cells of the infarcted human brain . CCL6 has chemoattractive activity on microglia and astrocyte  in rats. It is possible that in adult mice, strong expression of Ccl6 results in microglia attraction, and that the activated microglia might suppress axonal sprouting, which is not observed in neonatal mice.
Cd52 is expressed in lymphocytes, natural killer cells, eosinophils, macrophages, and dendritic cells . It was shown that an antibody against CD52 can effectively deplete CD52-expressing cells for the treatment of multiple sclerosis (a humanized monoclonal antibody called Alemtuzumab is used)  or in a mouse model (experimental autoimmune encephalomyelitis) . Although the function of CD52 itself is not fully understood, this protein regulates the activation of T cells . Therefore, it is possible that CD52 activates microglia after pyramidotomy in adult mice.
Ly86 and Mpeg1 also showed adult pyramidotomy group selective upregulation in qRT-PCR, although they showed slightly less than 2-fold change in RNA-seq (Fig. 5g, m-p), and we confirmed their expression in microglia/macrophage at the denervated side of CST after pyramidotomy in adult mice by in situ hybridization (Figs. 7 and 8). Ly86, also known as MD-1, associates with Radioprotective 105, and is involved in lipopolysaccharide-induced B cell growth . Ly86 is selectively expressed in microglia in central nervous system , and its expression in macrophage is elevated by pathogen stimulation . Mpeg1, which is a member of perforin family and also known as Perforin-2 , is originally identified as a macrophage specific gene , and its promoter is used for gene expression in macrophage  or microglia  in transgenic zebrafish. Its expression in macrophage is also upregulated after infection . Therefore, upregulation of these genes might represent adult specific activation of microglia/macrophage, which was also observed in immunohistochemistry for Iba1 (Fig. 6). It is possible that the density of degenerated axons of CST is high in the dorsal column, therefore microglia/macrophage engulfing the degenerated axons is strongly activated there, and the above genes were preferentially upregulated in them.
For most cases, our RNA-seq and qRT-PCR results were very similar, but for some genes, careful interpretation is needed. We collected RNA-seq samples at Zeitgeber time 1–7, and the order of sampling was P7 pyramidotomy, P7 sham, 8 W pyramidotomy, and 8 W sham (Additional file 2A). Sampling of P7 groups were conducted on the same day, and that of 8 W groups were conducted on another same day. We found intergroup differences in the expression levels of several circadian clock genes based on RNA-seq data (Additional file 2AB-AG), and accordingly, the timing of sampling might affect the expression levels of some genes. Therefore, for qRT-PCR validation, we collected all samples at Zeitgeber time 6 (Additional file 2A) on different days, although difference of sampling date might also cause other noise. As expected, the differences in the expression levels of circadian clock genes were smaller compared to those based on RNA-seq data (Additional file 2AB-AG). Because the circadian rhythm affects many aspects of biological processes such as stem cell function , it is possible that other minor changes in gene expression that could not be validated by qRT-PCR might reflect differences in the timing of sampling. Nevertheless, the reproducibility of most RNA-seq data was validated by qRT-PCR, indicating that our RNA-seq data essentially reflect the true biological phenomenon.
An important limitation in this study is that we used both intact and injured side of the spinal cord for RNA-seq. Because our model is left side pyramidotomy, responses are different between the left and right sides of the spinal cord. For example, intact axons sprout into the denervated side , and microglia is activated only at the denervated side of CST . Therefore, gene expression profile is probably different between injured and denervated sides. In a preliminary experiment, we tried to separate intact and denervated sides under a microscope, but we found that rapid and reliable separation is difficult especially for fragile and small neonatal spinal cord. Particularly, microglia activation at the denervated side of CST is observed very close to midline , and even a small amount of contamination of different side of CST will severely affect the result of RNA-seq. From this reason, we did not separate the intact and the denervated sides, and used both sides of the spinal cord. Including both sides of the spinal cord might be one of the reasons of the modest change of gene expression in our RNA-seq result.
Modest change of gene expression might also be due to inefficient sample selection for RNA-seq. We found that PKCγ signal did not disappear in adult mice 3 days after injury (Fig. 2c), and abnormal high intensity signal was observed at the lesion site. This is not a non-specific binding of the antibody, because we did not observe it when we used normal serum IgG as a negative control (Additional file 1). PKCγ signal at the denervated side of CST disappeared in most of the samples when we collected them 4 weeks after injury (Fig. 1d), whereas the abnormal high intensity signal was observed in most of the samples when we collected them at 3 dpi (Fig. 2c), suggesting that the signal is not derived from PKCγ in intact CST. This signal might be caused by strong binding between the antibody and partially degraded PKCγ in the injured CST. It is possible that the degenerated CST (and PKCγ) was not completely cleared by microglia yet at 3 dpi in adult mice. Because distinguishing abnormal high intensity signal and normal signal is more difficult than distinguishing no signal and normal signal, it is possible that we included imperfectly injured adult mice as RNA-seq sample, and this obscured the difference of gene expression.
Recently, sprouting neuron-specific genes have been identified by collecting cortical neurons positive for retrograde tracer injected into the denervated side of spinal cord after pyramidotomy in adult mice . It would be interesting to compare their gene expression profile with that of neonatal corticospinal neuron in future studies. Combining it with our present study, our understanding about neonatal sprouting would be much deepened.
In this study, remarkable difference of the gene expression profile of the spinal cord between intact neonatal and adult mice was revealed. In the neonatal spinal cord, genes related to axonal growth, cell proliferation, and myelination were upregulated, whereas those related to the immune response were downregulated compared to expression in adult mice without injury. In an attempt to reveal gene expression difference after pyramidotomy, we found that some genes involved in the inflammatory response were selectively upregulated in adult mice after pyramidotomy. Our research provides useful information about the differences between neonates and adults, and further research will reveal why sprouting ability is increased in neonates compared to that in adults.
C57BL/6 J wild-type female mice were obtained from the domestic company (Japan SLC). For the neonatal group, the mice were injured at P7, and for the adult group, the mice were injured at 8 weeks of age. The mice were maintained with a 12-h day 12-h night cycle and fed ad libitum in specific pathogen free condition. All surgical manipulations were performed under complete aesthesis via the intraperitoneal injection of medetomidine, midazolam, and butorphanol for adult mice or the inhalation of isoflurane for neonatal mice. The experiments were not blinded or randomized, and statistical tests were not used to pre-determine the sample size or to assess whether the data met the assumption of the approach. Health status of the animals were not examined before experiments. When the tissues were harvested, the animals were sacrificed by intraperitoneal injection of overdose of medetomidine (1.2 mg/kg, Orion Pharma), midazolam (16 mg/kg, Astellas), and butorphanol (20 mg/kg, Meiji Seika Pharma).
Pyramidotomy was performed essentially as described previously . In brief, the ventral side of the neck was incised, the trachea and the esophagus were retracted, and an injury of ~ 0.25 mm in depth ~ 0.5 mm in width was made on the left of the basal artery. For the sham group, skin was sutured after revealing the basal artery.
Axon tracer injection
Axon tracer injection was performed essentially as described previously . In brief, 0.6 μl of 10% biotinylated (Invitrogen, D1956) or rhodamine-conjugated dextran amine (Invitrogen, D1817) in phosphate buffered saline was injected into three points of the cerebral cortex (1 mm right/0 mm anterior, 1 mm right/0.5 mm anterior, and 1.5 mm right/0.5 mm anterior to the bregma, depth = 0.6 mm, total = 1.8 μl). This area corresponds to the forelimb area of the motor cortex . Injection was performed 2 weeks after pyramidotomy and mice were sacrificed 2 weeks after injection (4 weeks after pyramidotomy).
Immunohistochemistry for PKCγ was performed essentially as described previously . In brief, the cervical cord (vertebra C4–C7 level) was harvested after perfusion with 4% paraformaldehyde in phosphate buffer and post-fixed with the same solution for 2–6 h; alternatively, the medulla oblongata was harvested without perfusion with paraformaldehyde and post-fixed with fixative overnight. The tissues were frozen and 10–20-μm thick slices were cut and mounted on the glass slide. The specimen was incubated with 400 ng/ml of an anti-PKCγ antibody (Santa Cruz Biotechnologies, sc-211), and then incubated with secondary antibody conjugated with Alexa Fluor 488 (Invitrogen, A11008) or Alexa Fluor 568 (Invitrogen, A11011). The same concentration of IgG from normal rabbit serum (Sigma, I5006) was used as a negative control. For visualization of biotinylated dextran amine, the specimen was incubated with 4 μg/ml of fluorescence-conjugated streptavidin (Invitrogen). For enhancement of rhodamine signal, 0.1 μg/ml of an anti-rhodamine antibody (Invitrogen, A-6397) and Alexa Fluor 568 conjugated secondary antibody (A11011) were used as primary and secondary antibody, respectively, for free floating slices. Then, 1 μg/ml of 4′,6-Diamidino-2-phenylindole, dihydrochloride (DAPI) was used for nuclear staining.
For Iba1, the cervical cord was post fixed for three over nights after perfusion with paraformaldehyde. 1 μg/ml of an anti-Iba1 antibody (Wako, 019–19741) was used.
Quantification of sprouting
Quantification of sprouting was performed essentially as described previously . In brief, using an image of a coronal section of the cervical cord of a biotinylated or rhodamine-conjugated dextran amine-injected mouse, a horizontal line was drawn from the central canal to the limb of the gray matter at the denervated side, and the line was divided into three parts of equal length by two vertical lines, V1 and V2 (from medial to lateral). Then, the numbers of axons crossing V2, V1, or a vertical line at the central canal (Mid), in the gray matter were counted. The numbers were normalized by the number of axons at the main CST of the same slice. 10–20 slices from 3 to 4 mice (3–6 slices per mouse) were used in each group. The statistical analysis was performed using Tukey’s honestly significant difference (HSD) test.
The cervical cord (vertebra level C4–C7) tissue was homogenized in TRIzol Reagent (Invitrogen) 3 dpi. The cervical cords of neonatal mice were harvested on the same day, and those of adult mice were harvested on the same day, but neonatal and adult samples were harvested on the different days. The homogenate was stored at − 20 °C during evaluation of the lesion by immunohistochemistry. Chloroform was added to the homogenate, and the aqueous phase was transferred to the column of an RNeasy Mini Kit (Qiagen). The aqueous phase from 1 to 3 mice were pooled, and applied to the same column. Then, RNA was extracted following the manufacturer’s protocol with on column digestion of genomic DNA using the RNase-Free DNase Set (Qiagen). The library for sequencing was constructed from total RNA using the TruSeq stranded mRNA sample prep kit (Illumina). For this, 13–18 million 75-bp single-end reads/sample were produced using HiSeq 2500 (Illumina).
The following analysis was performed essentially as described previously  with modifications. In brief, the reads were mapped to the mouse genome GRCm38  using HISAT2 [60, 61]. Read counts were obtained using HTSeq  with Mus_musculus.GRCm38.91.gtf  as the reference gene model. The following analyses were performed using R  with appropriate packages. The statistical test for differentially expressed genes and clustering analysis based on the Euclidean distance matrix with variance stabilizing transformations were performed using DESeq2  following the developer’s protocol. Normalized read counts were obtained using the “counts” function in DESeq2 with the option “normalized = TRUE”. This means that row read counts were normalized by size factors. Size factors indicate difference of read depth between each samples, and is obtained by calculating change of gene expression level from ideal reference (geometric mean of the samples), and take median of the value, assuming that median of the change is taken from a non-differentially expressed gene. The Wald test was used and an adjusted P < 0.05 was considered significant. The annotation was performed using the org. Mm.eg.db package . The data were plotted on graphs using ggplot2  and ggbeeswarm  packages. The pathway analyses were performed using gage  package for Figs. 3 and 4 using read count data of all genes as input or DAVID [69, 70] for Fig. 5 using the 8 W Py UP gene list as input following the developer’s protocols. For gage, a two-sample t-test was used and an adjusted P (q value) < 0.05 was considered statistically significant. For DAVID, an adjusted P < 0.05 based on a Fisher’s exact test adjusted using the Benjamini correction was considered statistically significant.
RNA was extracted as described in the RNA-seq analysis section except that homogenates were not pooled, and the timing of sampling was different (Additional file 2A). Total RNA was reverse transcribed using the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems). Real-time PCR was performed using Fast SYBR Green Master Mix (Applied Biosystems) and QuantStudio 7 Flex Real-Time PCR System (Applied Biosystems) following the manufacture’s protocol. cDNA samples corresponding to final RNA amount of 6–9 ng were applied to the MicroAmp Fast Optical 384-Well Reaction Plate (Applied Biosystems). The PCR condition was as follows: 95 °C for 20 s followed by 40 cycles of 95 °C for 3 s and 60 °C for 30 s. The primers used for qRT-PCR are listed in Additional file 7. The expression levels were calculated by the relative standard curve method using reverse-transcribed product from the P7 sham group or the 8 W pyramidotomy group as standards and normalized by the geometric mean of the expression levels of glyceraldehyde-3-phosphate dehydrogenase (Gapdh) and ribosomal protein S18 (Rps18, , setting the value in the P7 sham group as 1. The Tukey’s honestly significant difference (HSD) test was used for statistical analysis with a cutoff value of P < 0.05.
Double staining of in situ hybridization and immunohistochemistry
0.5–1 kbp length of in situ hybridization probes were designed for Ly86 and Mpeg1. For Mpeg1, 5 probes targeting different region of Mpeg1 mRNA were mixed to enhance the signal as described previously .in situ hybridization was performed essentially as described previously  with modifications. In brief, the cervical cord was fixed as described in immunohistochemistry for Iba1 section. 20 μm thick slices were cut, and the specimens were processed by free-floating method. The specimen was incubated with 0.2 μg/ml of proteinase K solution for 30 min at 37 °C following acetylation. Probe was diluted with hybridization buffer (750 mM NaCl, 75 mM sodium citrate, 50% formaldehyde, 2% blocking reagent (roche, 1096176), 0.1% 3-[(3-Cholamidopropyl) dimethylammonio] propanesulfonate, and 0.1% sodium dodecyl sulfate) at the concentration of 1 ng/μl each, and the specimen was incubated in the probe solution over night at 60 °C after pre-hybridization for 1 h. The specimen was incubated with 1:1000 dilution of alkaline phosphatase-conjugated anti-Digoxigenin antibody (roche, 11093274910) for 2-5 h after RNase treatment, following incubation with nitro blue tetrazolium chloride/5-bromo-4-chloro-3-indolyl-phosphate, toluidine-salt solution for two over nights. Then the specimen was mounted on a glass slide, and incubated with anti-Iba1 antibody solution as described in immunohistochemistry for Iba1 section after washing with methanol and ethanol.
Availability of data and materials
The RNA-seq datasets produced in this study are deposited in DNA DataBank of Japan (accession number DRA007590).
Ankyrin repeat and ubiquitin domain containing 1
ADP-ribosylation factor-like 4D
Complement component 3
Calcitonin/calcitonin-related polypeptide, alpha
Chemokine (C-C motif) ligand
Collagen, type V, alpha 3
D-amino acid oxidase
Ethanolamine phosphate phospholyase
Ionized calcium binding adaptor molecule 1
Janus kinase–signal transducer and activator of transcription
Potassium voltage-gated channel, subfamily G, member 4
Kyoto Encyclopedia of Genes and Genomes
Lymphocyte antigen 86
Macrophage expressed gene 1
- NF-kappa B:
Nuclear factor kappa B
- NOD-like receptor:
Nucleotide oligomerization domain-like receptor
Nuclear receptor subfamily 1, group D, member 1
Period circadian clock 1
Protein kinase C gamma
POU domain, class 3, transcription factor 1
Quantitative reverse transcription-polymerase chain reaction
Ribosomal protein S18
Sema domain, transmembrane domain (TM), and cytoplasmic domain, (semaphorin) 6A
Serine (or cysteine) peptidase inhibitor, clade B, member 1a
Tumor necrosis factor alpha induced protein 6
Topoisomerase (DNA) II alpha
- Tukey’s HSD test:
Tukey’s honestly significant difference test
WAP four-disulfide core domain 18
Welniarz Q, Dusart I, Roze E. The corticospinal tract: evolution, development, and human disorders. Dev Neurobiol. 2017;77:810–29.
Onyszchuk G, Al-Hafez B, He YY, Bilgen M, Berman NEJ, Brooks WM. A mouse model of sensorimotor controlled cortical impact: characterization using longitudinal magnetic resonance imaging, behavioral assessments and histology. J Neurosci Methods. 2007;160:187–96. https://doi.org/10.1016/j.jneumeth.2006.09.007.
Starkey ML, Barritt AW, Yip PK, Davies M, Hamers FPT, McMahon SB, et al. Assessing behavioural function following a pyramidotomy lesion of the corticospinal tract in adult mice. Exp Neurol. 2005;195:524–39.
Ueno M, Hayano Y, Nakagawa H, Yamashita T. Intraspinal rewiring of the corticospinal tract requires target-derived brain-derived neurotrophic factor and compensates lost function after brain injury. Brain. 2012;135:1253–67.
Omoto S, Ueno M, Mochio S, Takai T, Yamashita T. Genetic deletion of paired immunoglobulin-like receptor B does not promote axonal plasticity or functional recovery after traumatic brain injury. J Neurosci. 2010;30:13045–52. https://doi.org/10.1523/JNEUROSCI.3228-10.2010.
Liu K, Lu Y, Lee JK, Samara R, Willenberg R, Sears-Kraxberger I, et al. PTEN deletion enhances the regenerative ability of adult corticospinal neurons. Nat Neurosci. 2010;13:1075–81. https://doi.org/10.1038/nn.2603.
Buchli AD, Schwab ME. Inhibition of Nogo: a key strategy to increase regeneration, plasticity and functional recovery of the lesioned central nervous system. Ann Med. 2005;37:556–67.
Lee DH, Luo X, Yungher BJ, Bray E, Lee JK, Park KK. Mammalian target of Rapamycin’s distinct roles and effectiveness in promoting compensatory axonal sprouting in the injured CNS. J Neurosci. 2014;34:15347–55. https://doi.org/10.1523/JNEUROSCI.1935-14.2014.
Gaudet AD, Fonken LK. Glial cells shape pathology and repair after spinal cord injury. Neurotherapeutics. 2018;15:554–77.
Geoffroy CG, Hilton BJ, Tetzlaff W, Zheng B. Evidence for an age-dependent decline in axon regeneration in the adult mammalian central nervous system. Cell Rep. 2016;15:238–46. https://doi.org/10.1016/j.celrep.2016.03.028.
Leong SK, Ling EA, Fan DP. Glial reaction after pyramidotomy in mice and rats. Neurodegeneration. 1995;4:403–13.
Zai L, Ferrari C, Subbaiah S, Havton LA, Coppola G, Strittmatter S, et al. Inosine alters gene expression and axonal projections in neurons contralateral to a cortical infarct and improves skilled use of the impaired limb. J Neurosci. 2009;29:8187–97. https://doi.org/10.1523/JNEUROSCI.0414-09.2009.
Maier IC, Baumann K, Thallmair M, Weinmann O, Scholl J, Schwab ME. Constraint-induced movement therapy in the adult rat after unilateral corticospinal tract injury. J Neurosci. 2008;28:9386–403. https://doi.org/10.1523/JNEUROSCI.1697-08.2008.
Bareyre FM, Haudenschild B, Schwab ME. Long-lasting sprouting and gene expression changes induced by the monoclonal antibody IN-1 IN the adult spinal cord. J Neurosci. 2002;22:7097–110.
Siebert JR, Middelton FA, Stelzner DJ. Intrinsic response of thoracic propriospinal neurons to axotomy. BMC Neurosci. 2010;11:69. https://doi.org/10.1186/1471-2202-11-69.
Fink KL, López-Giráldez F, Kim IJ, Strittmatter SM, Cafferty WBJ. Identification of intrinsic axon growth modulators for intact CNS neurons after injury. Cell Rep. 2017;18:2687–701. https://doi.org/10.1016/j.celrep.2017.02.058.
Shi LL, Zhang N, Xie XM, Chen YJ, Wang R, Shen L, et al. Transcriptome profile of rat genes in injured spinal cord at different stages by RNA-sequencing. BMC Genomics. 2017;18:173. https://doi.org/10.1186/s12864-017-3532-x.
Vitriol EA, Zheng JQ. Growth cone travel in space and time: the cellular Ensemble of Cytoskeleton, adhesion, and membrane. Neuron. 2012;73:1068–80. https://doi.org/10.1016/j.neuron.2012.03.005.
Myers JP, Santiago-Medina M, Gomez TM. Regulation of axonal outgrowth and pathfinding by integrin-ecm interactions. Dev Neurobiol. 2011;71:901–23.
Emoto K. The growing role of the hippo-NDR kinase signalling in neuronal development and disease. J Biochem. 2011;150:133–41.
Alvarez FJ, Benito-Gonzalez A, Siembab VC. Principles of interneuron development learned from Renshaw cells and the motoneuron recurrent inhibitory circuit. Ann N Y Acad Sci. 2013;1279:22–31.
Siembab VC, Smith CA, Zagoraiou L, Berrocal MC, Mentis GZ, Alvarez FJ. Target selection of proprioceptive and motor axon synapses on neonatal V1-derived Ia inhibitory interneurons and Renshaw cells. J Comp Neurol. 2010;518:4675–701.
Alexovič Matiašová A, Ševc J, Tomori Z, Gombalová Z, Gedrová Š, Daxnerová Z. Quantitative analyses of cellularity and proliferative activity reveals the dynamics of the central canal lining during postnatal development of the rat. J Comp Neurol. 2017;525:693–707.
Gaesser JM, Fyffe-Maricich SL. Intracellular signaling pathway regulation of myelination and remyelination in the CNS. Exp Neurol. 2016;283:501–11. https://doi.org/10.1016/j.expneurol.2016.03.008.
Chanas-Sacré G, Mazy-Servais C, Wattiez R, Pirard S, Rogister B, Patton JG, et al. Identification of PSF, the polypyrimidine tract-binding protein-associated splicing factor, as a developmentally regulated neuronal protein. J Neurosci Res. 1999;57:62–73. https://doi.org/10.1002/(SICI)1097-4547(19990701)57:1<62::AID-JNR7>3.0.CO;2-Y.
Khacho M, Slack RS. Mitochondrial dynamics in the regulation of neurogenesis: from development to the adult brain. Dev Dyn. 2018;247:47–53.
Solé-Domènech S, Cruz DL, Capetillo-Zarate E, Maxfield FR. The endocytic pathway in microglia during health, aging and Alzheimer’s disease. Ageing Res Rev. 2016;32:89–103. https://doi.org/10.1016/j.arr.2016.07.002.
Kigerl KA, de Rivero Vaccari JP, Dietrich WD, Popovich PG, Keane RW. Pattern recognition receptors and central nervous system repair. Exp Neurol. 2014;258:5–16. https://doi.org/10.1016/j.expneurol.2014.01.001.
Presumey J, Bialas AR, Carroll MC. Complement system in neural synapse elimination in development and disease. Adv Immunol. 2017;135:53–79. https://doi.org/10.1016/bs.ai.2017.06.004.
Nayak D, Roth TL, McGavern DB. Microglia development and function. Annu Rev Immunol. 2014;32:367–402. https://doi.org/10.1146/annurev-immunol-032713-120240.
Veiga-da-Cunha M, Hadi F, Balligand T, Stroobant V, Van Schaftingen E. Molecular identification of hydroxylysine kinase and of ammoniophospholyases acting on 5-phosphohydroxy-L-lysine and phosphoethanolamine. J Biol Chem. 2012;287:7246–55. https://doi.org/10.1074/jbc.M111.323485.
Schiroli D, Cirrincione S, Donini S, Peracchi A. Strict reaction and substrate specificity of AGXT2L1, the human O-phosphoethanolamine phospho-lyase. IUBMB Life. 2013;65:645–50. https://doi.org/10.1002/iub.1178.
Shao L, Vawter MP. Shared gene expression alterations in schizophrenia and bipolar disorder. Biol Psychiatry. 2008;64:89–97. https://doi.org/10.1016/j.biopsych.2007.11.010.
Murphy K. Janeway’s Immunobiology. 8th ed. New York: Garland Science; 2011.
Stanley ER, Chitu V. CSF-1 receptor signaling in myeloid cells. Cold Spring Harb Perspect Biol. 2014;6.
Stoll G, Jander S. The role of microglia and macrophages in the pathophysiology of the CNS. Prog Neurobiol. 1999;58:233–47.
Jin D, Liu Y, Sun F, Wang X, Liu X, He Z. Restoration of skilled locomotion by sprouting corticospinal axons induced by co-deletion of PTEN and SOCS3. Nat Commun. 2015;6:8074. https://doi.org/10.1038/ncomms9074.
Filli L, Engmann AK, Zörner B, Weinmann O, Moraitis T, Gullo M, et al. Bridging the gap: a Reticulo-propriospinal detour bypassing an incomplete spinal cord injury. J Neurosci. 2014;34:13399–410. https://doi.org/10.1523/JNEUROSCI.0701-14.2014.
Meletis K, Barnabé-Heider F, Carlén M, Evergren E, Tomilin N, Shupliakov O, et al. Spinal cord injury reveals multilineage differentiation of ependymal cells. PLoS Biol. 2008;6:e182. https://doi.org/10.1371/journal.pbio.0060182.
Shibata K, Nomiyama H, Yoshie O, Tanase S. Genome diversification mechanism of rodent and Lagomorpha chemokine genes. Biomed Res Int. 2013;2013:856265. https://doi.org/10.1155/2013/856265.
Nomiyama H, Egami K, Tanase S, Miura R, Hirakawa H, Kuhara S, et al. Comparative DNA sequence analysis of mouse and human CC chemokine gene clusters. J Interf Cytokine Res. 2003;23:37–45. https://doi.org/10.1089/10799900360520432.
Simats A, García-Berrocoso T, Penalba A, Giralt D, Llovera G, Jiang Y, et al. CCL23: a new CC chemokine involved in human brain damage. J Intern Med. 2018;283:461–75.
Asensio VC, Lassmann S, Pagenstecher A, Steffensen SC, Henriksen SJ, Campbell IL. C10 is a novel chemokine expressed in experimental inflammatory demyelinating disorders that promotes recruitment of macrophages to the central nervous system. Am J Pathol. 1999;154:1181–91.
Kanno M, Suzuki S, Fujiwara T, Yokoyama A, Sakamoto A, Takahashi H, et al. Functional expression of CCL6 by rat microglia: a possible role of CCL6 in cell-cell communication. J Neuroimmunol. 2005;167:72–80.
Zhao Y, Su H, Shen X, Du J, Zhang X, Zhao Y. The immunological function of CD52 and its targeting in organ transplantation. Inflamm Res. 2017;66:571–8.
Guarnera C, Bramanti P, Mazzon E. Alemtuzumab: a review of efficacy and risks in the treatment of relapsing remitting multiple sclerosis. Ther Clin Risk Manag. 2017;13:871–9.
Turner MJ, Pang PT, Chretien N, Havari E, LaMorte MJ, Oliver J, et al. Reduction of inflammation and preservation of neurological function by anti-CD52 therapy in murine experimental autoimmune encephalomyelitis. J Neuroimmunol. 2015;285:4–12. https://doi.org/10.1016/j.jneuroim.2015.05.018.
Inohara N, Nuñez G. ML – a conserved domain involved in innate immunity and lipid metabolism. Trends Biochem Sci. 2002;27:219–21.
Rosenberg AB, Roco CM, Muscat RA, Kuchina A, Sample P, Yao Z, et al. Single-cell profiling of the developing mouse brain and spinal cord with split-pool barcoding. Science. 2018;360:176–82. https://doi.org/10.1126/science.aam8999.
Jennings RT, Odkhuu E, Nakashima A, Morita N, Kobayashi T, Yamai I, et al. Inflammatory responses increase secretion of MD-1 protein. Int Immunol. 2016;28:503–12.
Podack ER, Munson GP. Killing of microbes and Cancer by the immune system with three mammalian pore-forming killer proteins. Front Immunol. 2016;7:464. https://doi.org/10.3389/fimmu.2016.00464.
Spilsbury K, O’Mara MA, Wu WM, Rowe PB, Symonds G, Takayama Y. Isolation of a novel macrophage-specific gene by differential cDNA analysis. Blood. 1995;85:1620–9.
Ellett F, Pase L, Hayman JW, Andrianopoulos A, Lieschke GJ. mpeg1 promoter transgenes direct macrophage-lineage expression in zebrafish. Blood. 2011;117:e49–56. https://doi.org/10.1182/blood-2010-10-314120.
Svahn AJ, Graeber MB, Ellett F, Lieschke GJ, Rinkwitz S, Bennett MR, et al. Development of ramified microglia from early macrophages in the zebrafish optic tectum. Dev Neurobiol. 2013;73:60–71. https://doi.org/10.1002/dneu.22039.
Fields KA, McCormack R, de Armas LR, Podack ER. Perforin-2 restricts growth of chlamydia trachomatis in macrophages. Infect Immun. 2013;81:3045–54. https://doi.org/10.1128/IAI.00497-13.
Weger M, Diotel N, Dorsemans AC, Dickmeis T, Weger BD. Stem cells and the circadian clock. Dev Biol. 2017;431:111–23. https://doi.org/10.1016/j.ydbio.2017.09.012.
Tennant KA, Adkins DL, Donlan NA, Asay AL, Thomas N, Kleim JA, et al. The organization of the forelimb representation of the C57BL/6 mouse motor cortex as defined by intracortical microstimulation and cytoarchitecture. Cereb Cortex. 2011;21:865–76.
Tsujioka H, Kunieda T, Katou Y, Shirahige K, Fukazawa T, Kubo T. Interleukin-11 induces and maintains progenitors of different cell lineages during Xenopus tadpole tail regeneration. Nat Commun. 2017;8:495. https://doi.org/10.1038/s41467-017-00594-5.
Zerbino DR, Achuthan P, Akanni W, Amode MR, Barrell D, Bhai J, et al. Ensembl 2018. Nucleic Acids Res. 2018;46:D754–61.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357–60.
Pertea M, Kim D, Pertea GM, Leek JT, Salzberg SL. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat Protoc. 2016;11:1650–67. https://doi.org/10.1038/nprot.2016.095.
Anders S, Pyl PT, Huber W. HTSeq-A Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9.
R Core Team. R: A language and environment for statistical computing. 2018. https://www.r-project.org/.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550. https://doi.org/10.1186/s13059-014-0550-8.
Carlson M. org.Mm.eg.db: Genome wide annotation for Mouse. 2018. https://doi.org/10.18129/B9.bioc.org.Mm.eg.db.
Wickham H. ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York; 2016. https://ggplot2.tidyverse.org/.
Clarke E, Sherrill-Mix S. ggbeeswarm: Categorical Scatter (Violin Point) Plots. 2017. https://cran.r-project.org/package=ggbeeswarm.
Luo W, Friedman MS, Shedden K, Hankenson KD, Woolf PJ. GAGE: generally applicable gene set enrichment for pathway analysis. BMC Bioinformatics. 2009;10:161. https://doi.org/10.1186/1471-2105-10-161.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57. https://doi.org/10.1038/nprot.2008.211.
Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13. https://doi.org/10.1093/nar/gkn923.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3:RESEARCH0034. doi:https://doi.org/10.1186/gb-2002-3-7-research0034.
Watakabe A, Komatsu Y, Ohsawa S, Yamamori T. Fluorescent in situ hybridization technique for cell type identification and characterization in the central nervous system. Methods. 2010;52:367–74. https://doi.org/10.1016/j.ymeth.2010.07.003.
We acknowledge Dr. Daisuke Motooka and Dr. Daisuke Okuzaki (the NGS core facility of the Genome Information Research Center, Research Institute for Microbial Diseases, Osaka University) for library construction and sequencing for RNA-seq analysis.
This study was partly supported by JSPS KAKENHI, Grant Number JP18K15367 to H.T. and JP17H06178 to T.Y., and AMED-CREST (18gm1210005h0001) to T.Y. The funding bodies had no role in design of the study, collection, analysis, or interpretation of the data, or writing this manuscript.
All animal experiments were approved by the Animal Experimentation Committee of Osaka University (permission number 29–058-037).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Negative control of immunohistochemistry for PKCγ. A representative image of negative control of Fig. 2c. Normal serum IgG was applied instead of anti-PKCγ antibody. Red indicate non-specific binding of IgG, and blue indicate DAPI signal. Dorsal is at the top, ventral is bottom, left is to the left, and right is to the right. A scale bar: 500 μm. (TIF 2989 kb)
Comparison of RNA-seq and qRT-PCR. A, The tissues for RNA-seq or qRT-PCR were harvested at the indicated time points. B-AI, Comparison of the expression levels measured by RNA-seq (B, D, F, H, J, L, N, P, R, T, V, X, Z, AB, AD, AF, AH; n = 3) and qRT-PCR (C, E, G, I, K, M, O, Q, S, U, W, Y, AA, AC, AE, AG, AI; n = 5–6). Based on the RNA-seq data, genes upregulated (B-S) or downregulated (T-W) in the adult sham group compared to those in the neonatal sham group, genes in the P7 Py UP group (X-AA), P7 Py DOWN group (AB-AG), or 8 W Py DOWN group (AH, AI) were chosen. Vertical axes represent normalized read counts (for RNA-seq) or relative expression levels normalized to those of the geometric mean of Gapdh and Rps18, setting the value of the P7 sham group to 1 (for qRT-PCR). Mean ± S.E.M. * adjusted P < 0.05, Wald test (for RNA-seq) or P < 0.05, Tukey HSD test (for qRT-PCR). Wfdc18, WAP four-disulfide core domain 18; Serpinb1a, serine (or cysteine) peptidase inhibitor, clade B, member 1a; Dao, D-amino acid oxidase; Ankub1, ankyrin repeat and ubiquitin domain containing 1; Tnfaip6, tumor necrosis factor alpha induced protein 6; Arl4d, ADP-ribosylation factor-like 4D; Kcng4, potassium voltage-gated channel, subfamily G, member 4; Ctss, cathepsin S; Eln, elastin; Sema6a, sema domain, transmembrane domain (TM), and cytoplasmic domain, (semaphorin) 6A; Col5a3, collagen, type V, alpha 3; Nes, nestin; Calca, calcitonin/calcitonin-related polypeptide, alpha; Nr1d1, nuclear receptor subfamily 1, group D, member 1; Per1, period circadian clock 1; Pou3f1, POU domain, class 3, transcription factor 1. (TIF 1830 kb)
Negative control of in situ hybridization for Ly86 and Mpeg1. Representative images of negative controls for Figs. 7 and 8 using sense probes for Ly86 (A-D) or Mpeg1 (E-H) are shown. The spinal cord of sham (A, C, E, G) or pyramidotomy (B, D, F, H) group which were injured at P7 (A, B, E, F) or 8 W (C, D, G, H) are shown. Right sides of the spinal cords are marked with black pigment. Blue/purple color is non-specific binding of sense probe. Dorsal is at the top, ventral is bottom, left is to the left, right is to the right. Scale bars: 200 μm. Because floating slices shrank during hybridization, the edge of the slices were folded in most cases when they were mounted on glass slides. (TIF 9508 kb)
KEGG pathways significantly upregulated in the P7 sham group compared to the 8 W sham group. KEGG pathways significantly (adjusted P < 0.05, two-sample t-test) upregulated in the P7 sham group compared to the 8 W sham group are shown. (XLSX 10 kb)
KEGG pathways significantly upregulated in the P7 pyramidotomy group compared to the P7 sham group. KEGG pathways significantly (adjusted P < 0.05, two-sample t-test) upregulated in the P7 pyramidotomy group compared to the P7 sham group are shown. (XLSX 17 kb)
KEGG pathways significantly upregulated in the 8 W pyramidotomy group compared to the 8 W sham group. KEGG pathways significantly (adjusted P < 0.05, two-sample t-test) upregulated in the 8 W pyramidotomy group compared to the 8 W sham group are shown. (XLSX 16 kb)
List of primers used for qRT-PCR. Forward and reverse primers used for qRT-PCR are listed. (XLSX 10 kb)