|Year : 2021 | Volume
| Issue : 7 | Page : 1258-1265
Cortical transcriptome analysis after spinal cord injury reveals the regenerative mechanism of central nervous system in CRMP2 knock-in mice
Ayaka Sugeno1, Wenhui Piao2, Miki Yamazaki3, Kiyofumi Takahashi4, Koji Arikawa4, Hiroko Matsunaga4, Masahito Hosokawa5, Daisuke Tominaga6, Yoshio Goshima7, Haruko Takeyama8, Toshio Ohshima MD, PhD 9
1 Laboratory for Molecular Brain Science, Department of Life Science and Medical Bioscience, Graduate School of Advanced Science and Engineering, Waseda University; Computational Bio Big-Data Open Innovation Laboratory (CBBD-OIL), National Institute of Advanced Industrial Science and Technology (AIST), Tokyo, Japan
2 Laboratory for Molecular Brain Science, Department of Life Science and Medical Bioscience, Graduate School of Advanced Science and Engineering, Waseda University, Tokyo, Japan
3 Biomolecular Engineering Laboratory, Department of Life Science and Medical Bioscience, Graduate School of Advanced Science and Engineering, Waseda University; Computational Bio Big-Data Open Innovation Laboratory (CBBD-OIL), National Institute of Advanced Industrial Science and Technology (AIST), Tokyo, Japan
4 Research Organization for Nano and Life Innovation, Waseda University, Tokyo, Japan
5 Biomolecular Engineering Laboratory, Department of Life Science and Medical Bioscience, Graduate School of Advanced Science and Engineering; Research Organization for Nano and Life Innovation; Institute for Advanced Research of Biosystem Dynamics, Waseda Research Institute for Science and Engineering, Graduate School of Advanced Science and Engineering, Waseda University, Tokyo, Japan
6 Computational Bio Big-Data Open Innovation Laboratory (CBBD-OIL), National Institute of Advanced Industrial Science and Technology (AIST), Tokyo; Cellular and Molecular Biotechnology Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Japan
7 Department of Molecular Pharmacology and Neurobiology, Graduate School of Medicine, Yokohama City University, Yokohama, Japan
8 Biomolecular Engineering Laboratory, Department of Life Science and Medical Bioscience, Graduate School of Advanced Science and Engineering, Waseda University; Computational Bio Big-Data Open Innovation Laboratory (CBBD-OIL), National Institute of Advanced Industrial Science and Technology (AIST); Research Organization for Nano and Life Innovation; Institute for Advanced Research of Biosystem Dynamics, Waseda Research Institute for Science and Engineering, Graduate School of Advanced Science and Engineering, Waseda University, Tokyo, Japan
9 Laboratory for Molecular Brain Science, Department of Life Science and Medical Bioscience; Institute for Advanced Research of Biosystem Dynamics, Waseda Research Institute for Science and Engineering, Graduate School of Advanced Science and Engineering, Waseda University, Tokyo, Japan
|Date of Submission||29-Jun-2020|
|Date of Decision||08-Sep-2020|
|Date of Acceptance||10-Oct-2020|
|Date of Web Publication||12-Dec-2020|
Laboratory for Molecular Brain Science, Department of Life Science and Medical Bioscience; Institute for Advanced Research of Biosystem Dynamics, Waseda Research Institute for Science and Engineering, Graduate School of Advanced Science and Engineering, Waseda University, Tokyo
Source of Support: This work was supported by Grants-in-Aid for Scientific Research on Priority Areas from The Ministry of Education, Culture, Sports, Science and
Technology (No. 26430043; to TO), Conflict of Interest: None
Recent studies have shown that mutation at Ser522 causes inhibition of collapsin response mediator protein 2 (CRMP2) phosphorylation and induces axon elongation and partial recovery of the lost sensorimotor function after spinal cord injury (SCI). We aimed to reveal the intracellular mechanism in axotomized neurons in the CRMP2 knock-in (CRMP2KI) mouse model by performing transcriptome analysis in mouse sensorimotor cortex using micro-dissection punching system. Prior to that, we analyzed the structural pathophysiology in axotomized or neighboring neurons after SCI and found that somatic atrophy and dendritic spine reduction in sensorimotor cortex were suppressed in CRMP2KI mice. Further analysis of the transcriptome has aided in the identification of four hemoglobin genes Hba-a1, Hba-a2, Hbb-bs, and Hbb-bt that are significantly upregulated in wild-type mice with concomitant upregulation of genes involved in the oxidative phosphorylation and ribosomal pathways after SCI. However, we observed substantial upregulation in channel activity genes and downregulation of genes regulating vesicles, synaptic function, glial cell differentiation in CRMP2KI mice. Moreover, the transcriptome profile of CRMP2KI mice has been discussed wherein energy metabolism and neuronal pathways were found to be differentially regulated. Our results showed that CRMP2KI mice displayed improved SCI pathophysiology not only via microtubule stabilization in neurons, but also possibly via the whole metabolic system in the central nervous system, response changes in glial cells, and synapses. Taken together, we reveal new insights on SCI pathophysiology and the regenerative mechanism of central nervous system by the inhibition of CRMP2 phosphorylation at Ser522. All these experiments were performed in accordance with the guidelines of the Institutional Animal Care and Use Committee at Waseda University, Japan (2017-A027 approved on March 21, 2017; 2018-A003 approved on March 25, 2018; 2019-A026 approved on March 25, 2019).
Keywords: CNS regeneration; cortex; CRMP2; hemoglobin; metabolic pathway; spinal cord injury; spine; transcriptome
|How to cite this article:|
Sugeno A, Piao W, Yamazaki M, Takahashi K, Arikawa K, Matsunaga H, Hosokawa M, Tominaga D, Goshima Y, Takeyama H, Ohshima T. Cortical transcriptome analysis after spinal cord injury reveals the regenerative mechanism of central nervous system in CRMP2 knock-in mice. Neural Regen Res 2021;16:1258-65
|How to cite this URL:|
Sugeno A, Piao W, Yamazaki M, Takahashi K, Arikawa K, Matsunaga H, Hosokawa M, Tominaga D, Goshima Y, Takeyama H, Ohshima T. Cortical transcriptome analysis after spinal cord injury reveals the regenerative mechanism of central nervous system in CRMP2 knock-in mice. Neural Regen Res [serial online] 2021 [cited 2021 Jan 28];16:1258-65. Available from: http://www.nrronline.org/text.asp?2021/16/7/1258/301035
Chinese Library Classification No. R459.9; R364; R741
| Introduction|| |
Mammalian central nervous system (CNS) has limited ability to regenerate and recover after injury. Spinal cord injury (SCI) is one such typical injury which is characterized by its semi-permanent impairments in sensorimotor functions below the spinal cord lesion. Though new treatment regimen for SCI by stem cell transplantation is gradually developing (Nakamura and Okano, 2013), promising approach for regenerating CNS with intrinsic factors is yet to be revealed.
Collapsin response mediator protein 2 (CRMP2) is one such candidate protein which is abundantly expressed in the adult mammalian brain (Wang and Strittmatter, 1996). CRMP2 is a member of the CRMP family proteins (CRMP1-5) which were identified as downstream molecules from Semaphorin 3A - an axon guidance cue (Schmidt and Strittmatter, 2007). CRMP2 has been reported to bind to microtubules, kinesin or dynein dependent vesicles, transport proteins, and interact with calcium channels (Kawano et al., 2005; Kimura et al., 2005; Wang et al., 2009; Rahajeng et al., 2010; Lin et al., 2011). Through these complex interactions, CRMP2 regulates synaptic transmission (Brittain et al., 2009; Wang et al., 2009), controls microtubule dynamics (Zheng et al., 2018), and modulates vesicle trafficking (Arimura et al., 2009; Rahajeng et al., 2010). The activity of CRMP2 is controlled by phosphorylation and Cdk5 is one of the major upstream kinases that phosphorylates CRMP2 at Ser522. Recently, in a mouse model, inhibition of phosphorylation at Ser522 (CRMP2KI) induces axon elongation after SCI along with reduced inflammation and microtubule stabilization at the spinal cord lesion site (Nagai et al., 2016). CRMP2KI mice showed partial recovery in the lost sensorimotor function subsequently and the phosphorylation inhibition of CRMP2 at Ser522 was suggested to have a therapeutic effect on SCI pathophysiology (Nagai et al., 2016).
The pathophysiology of SCI was also observed in the brain sensorimotor cortex. The axotomized neurons and neighborhood neurons significantly shrink and the number of synaptic spines reduces after SCI (Barron et al., 1988; Kim et al., 2006; Carter et al., 2008; Ghosh et al., 2012). These structural change is known to contribute to the rewiring of neuronal circuits within the cortex for compensating the lost sensorimotor functions (Ghosh et al., 2010). The effect of such changes in the cortex on the entire pathophysiology of SCI is not clearly understood and it has also not been investigated in the context of CRMP2 knock-in (CRMP2KI) mouse models.
In order to clear out such cortex pathology after SCI, we analyzed the structural pathophysiology and mRNA transcription in the cortex of CRMP2KI mice.
| Materials and Methods|| |
Animals and surgical procedures
Female wild-type (WT) and CRMP2KI (Crmp2KI/KI) mice (Yamashita et al., 2012) were crossed with YFP-H (thy1-Yellow Fluorescent Protein Mouse line H) or GFP-M (thy1-Green Fluorescent Protein Mouse line M) mice (Feng et al., 2000). Offspring at the age of 6–9 weeks, weighing 22–28 g, were used in the study.The fluorescent signal was used as a landmark to demarcate the region of interest–cortex layer V. Animals were randomly selected for sham or SCI models. For the SCI model, dorsal transection at T7–8 spinal cord was performed as previously described (Hill et al., 2009; Nagai et al., 2015). In brief, mice were anesthetized by inhalation with 2–5% isoflurane (DS Pharma Animal Health, Japan) and a T7–8 laminectomy was performed with 1.5-mm-deep dorsal transection. After surgery, manual bladder evacuation was done every day. All animals could easily access food and water. SCI mice who showed 0–1 Basso Mouse Scale (BMS) score (Basso et al., 2006) at 1 day after SCI were used for the experiments. For RNA sequencing, sham-operated model was made in order for precise analysis. In the condition of sham-operated model, mice were kept under anesthesia for 15–20 min. All animals were group-housed (N = 2 to 4 per cage) in plastic cages and maintained on a 12/12-hour light/dark cycle with the room temperature (RT) maintained at 25 ± 2°C. All these experiments were performed in accordance with the guidelines of the Institutional Animal Care and Use Committee at Waseda University, Japan (2017-A027 approved on March 21, 2017; 2018-A003 approved on March 25, 2018; 2019-A026 approved on March 25, 2019).
Soma size analysis by YFP signal
Mice were anesthetized using diethyl ether (Wako, Osaka, Japan) and transcardially perfused with 4% paraformaldehyde (PFA) (nacalai tesque, Kyoto, Japan) in phosphate buffer saline (pH 7.4), one or four weeks after surgery. Brains were dissected out from the animals and post-fixed overnight. They were continuously soaked in PBS, 10% sucrose and 20% sucrose for 12–24 hours and embedded in 2:1 mixture of Tissue Tek® optimal cutting temperature (OCT) compound (Sakura Finetek Japan, Osaka, Japan) and 20% sucrose. Frozen brains were sectioned into 30 µm sections by cryostat (Leica biosystems, Wetzlar, Germany). Slices containing bregma-posterior 0.5 mm, which has been reported to be sensorimotor cortex of the hindlimbs in mice by retrograde tracer (Ghosh et al., 2012), were observed by confocal microscope (FV-1000, Olympus, Tokyo, Japan) in 1.5 µm Z-stacks. Soma size of YFP-positive neurons in the cortex layer V between 500 and 1500 µm from the midline was measured with ImageJ software ver 1.52a (NIH, Bethesda, MD, USA) (Schneider et al., 2012) by encircling the soma.
Golgi-stain and spine analysis
Female mice at the age of 6 weeks were used for spine analysis for enhanced precision. SliceGolgi Kit (Bioenno Tech, LLC, Santa Ana, CA, USA) was used for Golgi staining technique. All the procedures were performed following SliceGolgi Kit instructions except the impregnation period. In brief, dissected brains fixed in fixative solution overnight at 4°C. Hindlimb sensorimotor cortex was sliced (200 µm) by vibratome (Micro Slicer DTK-1000, D.S.K, Kyoto, Japan) and Golgi-stained with impregnation period of 2.5–3 days. Apical dendrites of pyramidal neurons in cortex layer V were selectively observed with bright-field microscope (BX51, Olympus) within 7 days after Golgi-staining. The observed dendrites were selected randomly after confirming proper staining of the dendrites as well as ensuring the dendrites did not overlap with neighboring dendrites for the ease of analysis. All the spines and spine-like protrusions on the apical dendrites were counted with 60× objective lens (UPLFLN 60X, Olympus) and the length of the dendrites was measured with ImageJ software (Schneider et al., 2012).
Mice brain sections were prepared as described above. From bregma to bregma posterior 1.0 mm, the brains were sectioned into 15 μm. The slices were washed with PBS for 30 minutes, blocked with 3% horse serum in 0.1% PBSTr (PBS with 0.01% Triton X-100) for 10 minutes, 3% horse serum in 0.01% PBSTr for 30 minutes and incubated with anti-Hbα rabbit monoclonal IgG (1:50, Cat# SN70-09, Invitrogen, Carlsbad, CA, USA) overnight at 4°C. Slides were washed three times for 10 minutes each with 0.01% PBSTr the following day and incubated with Alexa Fluor 594 (1:1000, Abcam, Tokyo, Japan) rabbit IgG for 1 hour at room temperature. Following 3 times of 0.01% PBSTr washing, slices were incubated for 30 min with Hoechst 33258 (Fujifilm, Tokyo, Japan). Subsequently, the sections were washed with 0.01% PBSTr three times and mounted with FluoromountTM aqueous mounting medium (MERCK). Images were acquired by confocal laser scanning microscope in 1.5 µm Z-stacks (FV1000, Olympus).
Tissue dissection for RNA sequencing
Female mice at the age of 6 weeks were used for enhanced precision. The mice brains were dissected out 1 week after sham or SCI operation. Mice bodies were laid on ice to prevent degradation of RNA during the dissection process. Fresh brains were temporarily stored in PBS on ice, embedded in SCEM (Super Cryoembedding Medium, SECTION-LAB Co. Ltd., Hiroshima, Japan) using liquid nitrogen and stored at –80°C until sectioning. Fresh brains were sectioned at 30 µm thickness coronally with Cryostat (CM1860, Leica). The slices were mounted on FrameSlides PPS-Membrane 4.0 µm (Leica) and fixed immediately by air-drying. Dried-mounted slices were stored at room temperature or –20°C until use. Air-drying was conducted every time when the slices were taken out from –20°C freezer in order to prevent RNA degradation by condensation. Micro-dissection punching (Yoda et al., 2017) was performed for sampling cortex layer V, whose pyramidal neurons project axons to the spinal cord. The standard area of punching was determined by the YFP signal of YFP-H lines (Feng et al., 2000). Hindlimb sensorimotor cortex was also selected from micro-dissected samples as described previously (Ghosh et al., 2012). We confirmed the area similarly by Fast Blue tracer as well (Additional Figure 1).
Library construction and RNA sequencing
cDNA libraries were made after Poly (A)+ RNA extraction as per an available protocol with Smart-seq2 (Picelli et al., 2014). RNA quality was confirmed with RIN > 7 using bulk whole brain slice sample. Sequencing libraries were constructed by Nextera XT DNA library prep kit (Illumine®, San Diego, CA, USA). Constructed libraries were sequenced by paired-end way with Illumina MiSeq (Illumine®). Adapter sequences were trimmed with FLEXBAR (Dodt et al., 2012) and trimmed sequence reads were aligned to the Ensembl mouse reference genome (GRCm38 ver.92) using HISAT2 (Kim et al., 2015) with the default parameters. Annotated gene expression data was normalized as transcripts per million (TPM) value.
Differential gene expression and pathway analysis
The samples with high portion of mitochondrion genes (TPM ≥ 35%) and low number of detected genes (≤ 10,000 genes) were removed from analysis. Genes expressing more than 5 as averaged TPM value and detected in more than 4 samples from each condition (WT + sham, WT + SCI, CRMP2KI + sham, and CRMP2KI + SCI) were statistically analyzed as shown in statistics section. For pathway enrichment analysis, Gene Ontology (GO, http://www.geneontology.org), Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.jp/kegg/), Reactome (https://reactome.org/), and Hallmark gene sets from MSigDB (Liberzon et al., 2015) databases from Molecular Signatures Database (MSigDB, https://www.gsea-msigdb.org/gsea/msigdb/collections.jsp) (Liberzon et al., 2011) were used with clusterProfiler (Yu et al., 2012), ReactomePA (Yu and He, 2016), and msigdbr R packages. Fold change (FC) was calculated using averaged TPM value. FC evaluation for finding differentially expressed genes (DEGs) (Top DEGs) was performed by averaging intensity-based z-score. Concretely, genes were separated into several groups by averaged TPM from sham samples and z-score was calculated from the FC values within the group. This procedure was performed repeatedly with different grouping threshold and z-scores were averaged. Genes showing the highest or the lowest averaged z-score were defined to be the top DEGs.
Quantitative real-time PCR validation
cDNA libraries for RT-qPCR validation were constructed with the same procedure as RNA sequencing library construction. SYBR Green PCR Master Mix (Applied Biosystems, Foster City, CA, USA) was used for qPCR. Reactions were carried out in triplicates with 10 ng of cDNA and 200 nM of each primer in a final volume of 20 μL. The reaction protocol was 95°C for 1 minute; 40 amplification cycles at 95°C for 15 seconds, 50°C for 30 seconds, 72°C for 1 minute. Relative quantification of gene expression was analyzed with ΔΔCt method and the internal control ACTB, whose mRNA expression did not change statistically in the RNA sequencing experiments. Genes and primers used for validation are as follows: ACTB-F (TTG TGT AAG GTA AGG TGT GC), ACTB-R (GGT TGA GGT GTT GAG GCA G), Hba-F (ATG CCC ACA AGC TGC GT), Hba-R (GTG CTC ACA GAG GCA AGG AA), Hbb-F (GCC CAG CAC AAT CAC GA), and Hbb-R (TGC CTT TAA CGA TGG CCT GA). In order to calculate relative expression ratio of Hb genes, ACTB was used as a housekeeping control for qPCR analysis.
For the comparison of soma size distribution analysis, Kolmogorov-Smirnov test was performed with GraphPad Prism software version 6.0b (GraphPad Software Inc., San Diego, CA, USA). Statistical differences of spine analysis were calculated with an unpaired two-tailed Student’s t-test. DEGs were defined using Wilcoxon test. Adjusted P-values were calculated with Benjamini-Hochberg method. Genes and pathways with P-value < 0.05 were defined as significantly changed or enriched.
| Results|| |
CRMP2KI suppresses somatic atrophy in sensorimotor cortex after spinal cord injury
A previous study demonstrated that neurons of the corticospinal tract showed somatic atrophy after SCI in rat and mice (Barron et al., 1988; Carter et al., 2008). In mice, cortical layer V neurons undergo severe atrophy at 4 weeks after SCI (Carter et al., 2008). In order to check whether these changes occur upon the inhibition of phosphorylation of CRMP2, we examined the size of the sensorimotor neurons in cortex layer V observing YFP signal. In WT mice, sensorimotor neurons showed severe somatic atrophy at 1 week after SCI and this atrophy continued until 4 weeks after SCI [Figure 1]A and [Figure 1]B. However, in CRMP2KI mice, sensorimotor neurons shrank more mildly at 1 week after SCI and recovered to their original size at 4 weeks after SCI [Figure 1]A and [Figure 1]B. Between WT and CRMP2KI mice, significant difference was observed for both 1 week and 4 weeks after SCI (P < 0.05; [Figure 1]C.
|Figure 1: Constitutively activated CRMP2 suppresses soma atrophy in corticospinal neurons after SCI.|
(A) YFP-labeled Layer V neurons (green) show severe atrophy 1 week after SCI and partially recovered 3 weeks later. N = 3, 4, 3, 3, 4, 3 for wild-type (WT), WT + SCI 1 wk, WT + SCI 4 wk, CRMP2KI, CRMP2KI + SCI 1 wk, CRMP2KI + SCI 4 wk, respectively. Scale bars: 20 μm. (B) Distribution of soma size of YFP-labeled layer V after SCI in WT (above) and CRMP2KI (below) mice (Kolmogorov-Smirnov test, *P < 0.05, **P < 0.01, ***P < 0.001). (C) Distributions of soma size were compared between WT and CRMP2KI mice at 1 and 4 weeks after injury. CRMP2: Collapsin response mediator protein 2; KI: knock-in; SCI: spinal cord injury.
Click here to view
CRMP2KI suppresses postsynaptic spine reduction in sensorimotor cortex after spinal cord injury
Neurons of the corticospinal tract have fewer post-synaptic spines after SCI (Ghosh et al., 2012; Kim et al., 2006). In order to reveal the effect of phosphorylation inhibition of CRMP2 on synaptic pathology, we analyzed the spine density. The location of analysis was apical dendrites of the pyramidal neurons. The number of spines decreased 1 week after SCI in WT mice [Figure 2]. However, CRMP2KI mice did not show any changes in spine density [Figure 2]B. CRMP2KI mice had originally fewer spines on the pyramidal neurons in cortex layer V at 7 weeks.
|Figure 2: Constitutively activated CRMP2 suppresses spine decrease in CST neurons after SCI.|
(A) Golgi-stained layer V neurons after SCI in WT and CRMP2 KI mice. Scale bars: 50 μm. (B) Spine density of layer V neurons was calculated from apical dendrites (n = 17, 25, 25, 30 for WT, WT + SCI 1 wk, CRMP2KI, CRMP2KI + SCI 1 wk, respectively). Data are presented as the mean ± SEM. ***P < 0.001 (Student’s t-test). CRMP2: Collapsin response mediator protein 2; KI: knock-in; SCI: spinal cord injury; WT: wild-type.
Click here to view
Cortical transcriptome profiles in all comparisons
We found the structural pathologies in sensorimotor cortex were suppressed in synaptic and somatic levels by the phosphorylation inhibition of CRMP2. However, the intracellular response in axotomized neurons is not fully understood. To reveal the detailed pathophysiology in axotomized neurons of the cortex, we performed RNA sequencing for micro-dissected cortex layer V [Figure 3]A, 1 week after SCI in WT and CRMP2KI mice. For the precise detection of DEGs, we performed statistical analysis for each pair of groups instead of multiple comparisons. We identified 967 genes in WT mice and 811 genes in CRMP2KI mice as DEGs after SCI. The mutation effect of CRMP2 phosphorylation inhibition was analyzed using the sham model samples of WT and CRMP2KI mice. A total of 1674 DEGs were found in CRMP2KI mice compared with WT mice. PCA analysis was performed with all DEGs’ TPM data and micro-dissected samples were separated well reflecting experimental conditions [Figure 3]B. We found that the number of overlapped DEGs did not account for the high percentage of whole DEGs [Figure 3]C.
|Figure 3: RNA sequencing analysis of micro-dissected samples.|
RNA sequencing analysis was performed for micro-dissected samples obtained from the hindlimb sensorimotor cortex layer V. Each condition has 5–7 micro-dissected samples from 3, 4 mice. (A) Mouse cortex slice of bregma posterior 0.5 mm. YFP fluorescence was used to mark the dissection regions. Arrows indicate the micro-dissected regions of 100 μm diameter. (B) Principal component analysis (PCA) of micro-dissected samples. Transcripts per million (TPM) data of all differentially expressed genes (DEGs) found in three comparisons, i.e., SCI vs. sham in wild-type (WT), SCI vs. sham in CRMP2KI, CRMP2KI sham vs. WT sham, were used. (C) The number of overlapped DEGs found in three comparisons was shown in Venn diagram. CRMP2: Collapsin response mediator protein 2; KI: knock-in; SCI: spinal cord injury.
Click here to view
Differential gene expression by the phosphorylation inhibition of CRMP2
We analyzed the mutation effect of CRMP2 phosphorylation inhibition. We detected 1674 DEGs. Among them, 1176 genes were upregulated, while 498 genes were downregulated (Additional Figure 2). Pathway enrichment analysis revealed that many pathways related to intracellular metabolism and neuronal functions were differentially expressed in CRMP2KI sensorimotor cortex [Table 1] and Additional Table 1 and 2). Further, genes involved in the cytosolic or mitochondrial translation and oxidative phosphorylation, proteasome, spliceosome pathways were significantly upregulated (P < 0.05, q < 0.05), and fatty acid metabolism, lysosome, neuronal system, synaptic pathway genes were downregulated (P < 0.05, q < 0.2). In addition to the oxidative phosphorylation and ribosomal pathways in mitochondria, genes regulating mitochondrion organization were also upregulated (Additional Figure 2C). GO analysis also displayed that neuronal specific structural or functional terms such as axon, and vesicles were enriched among downregulated DEGs (Additional Figure 2D). This suggests that the inhibition of phosphorylation of CRMP2 might dynamically change the intracellular metabolic system and neuronal functions. We also analyzed the top 20 significant DEGs (Additional Tables 3 and 4). Various genes whose function is restricted to cytoskeletal protein binding like Nrcam, Map1b, and Vapa were greatly downregulated.
|Table 1: Enriched biological pathways of differentially expressed genes after SCI in WT and CRMP2KI mice|
Click here to view
Differential gene expression after spinal cord injury in WT and CRMP2KI mice
We detected 967 DEGs after SCI in WT mice. Of them, 886 DEGs showed upregulation. Pathway analysis revealed that cytosolic or mitochondrial translation and oxidative phosphorylation pathway genes were significantly enriched among the upregulated genes ([Table 1] and Additional table 5). From the limited number of 81 downregulated genes, we could not detect any enriched pathways. Interestingly, hemoglobin genes Hba-a1, Hba-a2, Hbb-bs, and Hbb-bt were remarkably upregulated and became the top 4 upregulated genes (FC > 6.9) after SCI [Table 2].
|Table 2: Top 10 significantly upregulated DEGs after spinal cord injury in WT mice|
Click here to view
On the other hand, in CRMP2KI mice, we detected a total of 811 DEGs in which 224 were upregulated; while 587 were downregulated after SCI (Additional Figure 3). Upregulated cluster consisted mainly of channel and transporter activity genes [Table 1] and Additional figure 3C). From the downregulated genes, neuronal system and cytosolic translation pathway genes were specifically enriched. Genes regulating synaptic activity and organization of axon, myelin sheath, vesicles, and synaptic membranes were all downregulated ([Figure 4]A and Additional Table 6). Notably, gliogenesis, glial cell differentiation, and developmental genes were also enriched in the downregulated DEGs after SCI in CRMP2KI mice (Additional Table 6). Hemoglobin genes were not differentially expressed in CRMP2KI mice because of the broad distribution of TPM values.
|Figure 4: Differential gene expression analysis after SCI.|
(A) Cellular component GO enrichment analysis of downregulated genes after SCI in CRMP2KI mice. Many neuronal specific genes were enriched in addition to ribosomal pathway genes. P. adust: Adjusted P-value calculated with Benjamini-Hochberg method. (B) Heatmap using the differentially expressed genes (DEGs) after SCI in wild-type (WT) mice. The clusters containing characteristic enriched pathways and hemoglobin genes have been shown alongside. (n = 5–7, N = 3,4). CRMP2: Collapsin response mediator protein 2; KI: knock-in; SCI: spinal cord injury.
Click here to view
As summarized in [Table 1], oxidative phosphorylation pathway genes were upregulated in both CRMP2KI mutation and SCI. In order to compare these DEGs expression level in other condition, heatmap was generated for DEGs after SCI in WT mice and pathway analysis was performed for several clusters [Figure 4]B. We found that oxidative phosphorylation pathway and mitochondrial component genes were expressed at a high level in both CRMP2KI + sham and CRMP2KI + SCI condition compared with WT + SCI condition.
Hemoglobin expression in sensorimotor cortex after SCI
In transcriptome analysis, four hemoglobin genes were remarkably upregulated after spinal cord injury in WT mice [Table 2]. These hemoglobin expressions were validated with RT-qPCR at first. cDNA libraries were prepared from micro-dissected sensorimotor cortex. Both Hba and Hbb genes were significantly upregulated after SCI in WT mice [Figure 5]A. However, there was no significant difference between sham and SCI in CRMP2KI mice. These results were consistent with the transcriptome analysis obtained from RNA sequencing.
|Figure 5: Hemoglobin (Hb) was upregulated after SCI and increased at protein level.|
(A) mRNA levels of Hba and Hbb were validated by RT-qPCR from the micro-dissected samples 1 week after SCI (n = 6, 6, 4, 4 for wild-type (WT) sham, WT SCI, CRMP2KI sham, CRMP2KI SCI, respectively). *P < 0.05 (Student’s t-test). (B) Immunostaining of Hbα in hindlimb sensorimotor cortex in WT mice. The ratio of Hbα positive cells was analyzed with respect to Hoechst staining (3–7 micro-dissected samples were obtained from five mice). Scale bars: 100 μm. Data are presented as the mean ± SEM. CRMP2: Collapsin response mediator protein 2; KI: knock-in; SCI: spinal cord injury.
Click here to view
Next, we performed immunostaining of Hbα in WT mice to analyze if this hemoglobin upregulation could occur at the protein level in the sensorimotor cortex after SCI. The ratio of Hbα positive cells increased after SCI but not significantly (P = 0.063) unlike our expectation [Figure 5]B.
| Discussion|| |
The inhibition of CRMP2 phosphorylation has been proven to be effective for the diagnosis of SCI by this study. Previous study showed that the phosphorylation inhibition of CRMP2 suppresses glial scar formation and inflammation at the spinal cord lesion, which leads to axonal regeneration (Nagai et al., 2016). In this study, we found that spine reduction and somatic atrophy were suppressed in the sensorimotor cortex in CRMP2KI mice. Spine reduction is considered to regulate the rewiring of neuronal circuits as compensation for lost sensorimotor function and brain activity map changed in rat SCI model (Ghosh et al., 2010, 2012; Zhang et al., 2015). We show that CRMP2KI mice have the same spine density as the originally fewer density, 1 week after SCI. According to the previous study (Jin et al., 2016), fewer spines were also observed in CA1 pyramidal neurons of CRMP2KI mice at the age of 5 weeks. Focusing on the effect of SCI, the phosphorylation inhibition of CRMP2 suppressed further spine reduction by SCI. Although it is necessary to take into account that there might be relevance between suppression of spine reduction after SCI and the original fewer spines in CRMP2KI mice, the phosphorylation inhibition of CRMP2 suppressed further post-synaptic spine reduction by SCI from the primary number. This supports the idea that it could be possible in CRMP2KI mice, there is only a mild alteration of electrical activity in the somatosensory cortex, with concomitant suppression of somatic atrophy and spine reduction. Simultaneously, somatic atrophy gave an indication of neuronal stress after SCI. Significant somatic atrophy is observed 4 weeks after SCI in WT mice, but somatic size in CRMP2KI regained the original size on the same time scale. This result suggests that phosphorylation inhibition of CRMP2 contribute to suppressing the somatic atrophy of axotomized neurons and facilitating the recovery toward original size after SCI. This size recovery is also observed in the mice treated with ChABC which have retrograde neuroprotective effect after SCI (Carter et al., 2008). This gave us the possibility that the phosphorylation inhibition of CRMP2 also have similar neuroprotective effect in axotomized neurons. Considering the sensorimotor structural pathology after SCI, the phosphorylation inhibition of CRMP2 contributes to suppressing SCI pathology not only at the spinal cord lesion site; but also in the intracellular response of CNS neurons.
We also performed RNA sequencing analysis in micro-dissected sensorimotor cortex to reveal the detailed intracellular response in axotomized neurons after SCI in WT and CRMP2KI mice. Each cortical transcriptome profile has different characteristics when subjected to SCI, CRMP2 mutation and both. We observed that genes involved in the cytosolic, mitochondrial translation and oxidative phosphorylation pathways were significantly upregulated after SCI. Previous studies showed that genes involved in similar pathways like translation, RNA splicing, and oxidative phosphorylation were differentially expressed within a week after SCI in the cortex (Jaerve et al., 2012; Baek et al., 2017). This showed that our transcriptome analysis with micro-dissected cortex is consistent with the previous studies. One of these studies also suggested that impaired oxidative phosphorylation and mitochondrial dysfunction occur after SCI (Baek et al., 2017). We also found upregulation of mitochondrion component genes. Oxidative phosphorylation pathway takes place in mitochondria and impaired oxidative phosphorylation induces oxidative stress by generating reactive oxygen species (ROS) in cells generally. ROS and oxidative stress have been considered to be a hallmark of SCI and progressive SCI pathology at spinal cord lesion sites (Xu et al., 2005; Jia et al., 2012). In addition, one transcriptome analysis in cortex layer V also detected differential expression in response to oxidative stress (Kruse et al., 2011). Our findings regarding upregulated genes and enriched pathways after SCI in WT also suggests the oxidative stress as well. Notably, we further detected 4 hemoglobin (Hb) genes remarkably increased after SCI in WT mice, which was also validated by qPCR. Recently, Hb was found to be expressed in CNS tissues, including cortical neurons. Also, some studies on neurodegenerative diseases such as Alzheimer’s disease, Multiple Sclerosis, and Parkinson’s disease suggested the relevance between Hb and their pathophysiology (Biagioli et al., 2009; Richter et al., 2009; Schelshorn et al., 2009; Brown et al., 2016; Freed and Chakrabarti, 2016; Altinoz et al., 2019). In the field of SCI studies, this is the first report that indicates the involvement of Hb in SCI pathophysiology. The role of Hb in nervous system has not been fully understood and continuous studies were needed for clear explanation. However Hb is suggested to be involved in O2 homeostasis and mitochondrial homeostasis by buffering protons by other studies (Biagioli et al., 2009; Freed and Chakrabarti, 2016; Altinoz et al., 2019). Though we could not detect statistically significant upregulation of Hb at the protein level, Hb transcription upregulation and simultaneously upregulated Hb functions, oxidative phosphorylation gave us some idea that the involvement of oxidative stress or impairment of mitochondrial homeostasis in axotomized neurons after SCI.
In this study, we reveal for the first time the mRNA transcription profile in sensorimotor cortex in CRMP2KI in sham-operated and SCI condition. Completely different pathways were up- or down-regulated after SCI. All the enriched pathways except the pathways involved in translation were related to typical functions of nervous system and most of the enriched pathways were downregulated after SCI. The point to be emphasized is that genes regulating synaptic compartments and functions were strongly downregulated simultaneously with the suppression of spine reduction after SCI. In addition, similar synaptic pathways were downregulated by the genetic modification of CRMP2 originally. This means axotomized neurons in CRMP2KI have quite lower mRNAs regulating synaptic functions. Although further studies are needed for clarifying the concrete mechanism, it is likely that downregulation of synaptic genes contributes to the suppression of spine reduction after SCI in CRMP2KI and might contribute to the improvement of SCI pathophysiology as well. We also found an interesting observation that glial cell differentiation and myelin sheath pathways were downregulated after SCI in CRMP2KI mice. Besides, some studies showed that microglial activation appeared in cortex (Wu et al., 2014a, b; Jure and Labombarda, 2017). We could not detect the pathway of microglial activation after SCI in WT mice, but downregulation of glial cells differentiation pathway in CRMP2KI indicates that there might be some changes in the response of glial cells including microglia, astrocytes and oligodendrocytes upon the phosphorylation inhibition of CRMP2. Indeed, the glial scar formation and inflammation at spinal cord lesion was suppressed in CRMP2KI mice (Nagai et al., 2016). Our transcriptome analysis in sensorimotor cortex indirectly supports the different response after SCI in glial cells in CRMP2KI. In addition, axon elongation was also confirmed after SCI in CRMP2KI (Nagai et al., 2016) mice. The downregulation of genes regulating myelin sheath could be involved in suppressed retrograde axon degeneration in CRMP2KI mice. To summarize, the downregulation of these pathways in CRMP2KI may be directly or indirectly involved in pathophysiology of SCI, spine reduction, glial scar formation, and axon degeneration.
It is also important that we found several metabolic intracellular pathways and synaptic compartments are differentially expressed in CRMP2KI mice. Genes involved in translation, oxidative phosphorylation and proteasome pathways were upregulated and fatty acid metabolism, lysosome pathways, and synaptic compartments were downregulated. These changes in the metabolic pathways might come from protein-protein interactions of CRMP2. Previously, CRMP2 interactome analysis indicated that CRMP2 interacts with ten proteins that regulate energy metabolism and nine proteins that regulate protein metabolism including ribosomal subunit (Martins-de-Souza et al., 2015). The interactions between these proteins and mutated CRMP2 is not clear, but the phosphorylation inhibition of CRMP2 could change these interaction levels as well and affect the metabolic system overall. Our transcriptome results are partly consistent with and partly contradict the previous proteome study interestingly. The protein levels of molecules regulating oxidative phosphorylation, proteasome and fatty acid metabolism pathways were actually decreased in CRMP2KI mice (Nakamura et al., 2018). It is possible that the phosphorylation inhibition of CRMP2 broadly induce the improvements in SCI pathophysiology via intracellular metabolic pathways. In addition, a previous phosphoproteome analysis in CRMP2KI mice indicated proteins regulating synaptic functions had different phosphorylation profile as well (Nakamura et al., 2018). Our results that genes involved in regulating synaptic compartments were differentially expressed in CRMP2KI mice might be caused by this alteration of phosphorylation in synaptic proteins. We also found that genes regulating oxidative phosphorylation and mitochondrial components were even more upregulated in both CRMP2KI condition than SCI in WT mice. It is possible that the phosphorylation inhibition of CRMP2 might already have similar condition as observed after SCI in WT mice.
In conclusion, our study showed many possibilities of SCI pathophysiology in axotomized neurons in WT and CRMP2KI mice. Further, the transcriptome profile of CRMP2KI was revealed in detail. We confirmed that CRMP2KI suppressed the structural pathophysiology of spine reduction and somatic atrophy after SCI in cortex. In addition, translation and oxidative phosphorylation pathways were upregulated in cortex after SCI. Interestingly, fatty acid metabolism, lysosome, and neuronal system relating pathways were downregulated originally by phosphorylation inhibition of CRMP2 though translation, oxidative phosphorylation, and proteasome pathways were increased. After SCI in this CRMP2 mutation, neuronal system, synaptic, translation, glial cell differentiation, and myelin sheath pathways were decreased though only channel activity function was upregulated. This transcriptome analysis indicated (1) CRMP2KI seems to influence the whole metabolic system in CNS such as oxidative phosphorylation, fatty acid metabolism, and proteasome pathways, (2) these metabolic systems might contribute to the SCI pathophysiology in axotomized neurons as well and (3) response to SCI in CRMP2KI mice seems to involve synaptic activity and glial cell differentiation in addition to the putatively effective pathways such as microtubules stabilization. Hb involvement in SCI pathology is also suggested for the first time in this study. However in this study we could not cut in the functional level. In spite of the suggestions from mRNA transcriptome, locations that each protein works and protein interaction after SCI are not directly revealed. Cell-to-cell interaction in cortex is also unclear. At this point, our methods have some limitations.
The intracellular mechanisms underlying for axon elongation and functional recovery in CRMP2KI or Hb functions in axotomized neurons are still unclear in detail. In order to clarify the concrete mechanism on CNS regeneration, the complicated SCI pathophysiology and pathophysiological phenomena in CRMP2KI need to be further investigated.
Acknowledgments: We thank Sachiyo Aburatani, Hisashi Anbutsu and Yoriko Uotani for support in OIL, and support by the Platform Project for Supporting Drug Discovery and Life Science Research (Basis for Supporting Innovative Drug Discovery and Life Science Research (BINDS)) from AMED under Grant Number JP20am0101104.
Author contributions: AS, WP, MY, and KT performed experiments. AS, DT, HM, KA, and MH analyzed RNA sequencing data. TO and YG developed CRMP2KI mice. HT and TO supervised the overall project. AS and TO wrote the manuscript.
Conflicts of interest: The authors declare no conflict of interest.
Financial support: This work was supported by Grants-in-Aid for Scientific Research on Priority Areas from The Ministry of Education, Culture, Sports, Science and Technology (No. 26430043; to T.O)
Institutional review board statement: All these experiments were performed in accordance with the guidelines of the Institutional Animal Care and Use Committee at Waseda University, Japan (2017-A027 approved on March 21, 2017; 2018-A003 approved on March 25, 2018; 2019-A026 approved on March 25, 2019).
Copyright license agreement: The Copyright License Agreement has been signed by all authors before publication.
Data sharing statement: Datasets analyzed during the current study are available from the corresponding author on reasonable request.
Plagiarism check: Checked twice by iThenticate.
Peer review: Externally peer reviewed.
Additional Figure 1: Fast blue tracer positive region of sensorimotor cortex[Additional file 1].
Additional Figure 2: Transcriptome analysis between WT and CRMP2KI mice[Additional file 2].
Additional Figure 3: Transcriptome analysis after spinal cord injury in CRMP2KI[Additional file 3].
Additional Table 1: Enriched pathways of upregulated DEGs by the phosphorylation deficient of CRMP2 (P < 0.05, q < 0.05)[Additional file 4].
Additional Table 2: Enriched pathways of downregulated DEGs by the phosphorylation deficient of CRMP2 (P < 0.05, q < 0.2)[Additional file 5].
Additional Table 3: Top20 significantly upregulated DEGs by the phosphorylation deficient of CRMP2[Additional file 6].
Additional Table 4: Top20 significantly downregulated DEGs by the phosphorylation deficient of CRMP2[Additional file 7].
Additional Table 5: Enriched pathways of upregulated DEGs after spinal cord injury in WT mice (P < 0.05, q < 0.05)[Additional file 8].
Additional Table 6: Enriched pathways of downregulated DEGs after spinal cord injury in CRMP2KI mice (P < 0.05, q < 0.05)[Additional file 9].
Funding: This work was supported by Grants-in-Aid for Scientific Research on Priority Areas from The Ministry of Education, Culture, Sports, Science and Technology (No. 26430043; to TO).
| References|| |
Altinoz MA, Guloksuz S, Schmidt-Kastner R, Kenis G, Ince B, Rutten BPF (2019) Involvement of hemoglobins in the pathophysiology of Alzheimer’s disease. Exp Gerontol 126:110680.
Arimura N, Kimura T, Nakamuta S, Taya S, Funahashi Y, Hattori A, Shimada A, Ménager C, Kawabata S, Fujii K, Iwamatsu A, Segal RA, Fukuda M, Kaibuchi K (2009) Anterograde transport of TrkB in axons is mediated by direct interaction with Slp1 and Rab27. Dev Cell 16:675-686.
Baek A, Cho SR, Kim SH (2017) Elucidation of gene expression patterns in the brain after spinal cord injury. Cell Transplant 26:1286-1300.
Barron KD, Dentiger MP, Popp AJ, Makes R (1988) Neurons of layer Vb of rat sensorimotor cortex atrophy but do not die after thoracic cord transection. J Neuropathol Exp Neurol 47:62-74.
Basso DM, Fisher LC, Anderson AJ, Jakeman LB, McTigue DM, Popovich PG (2006) Basso Mouse Scale for locomotion detects differences in recovery after spinal cord injury in five common mouse strains. J Neurotrauma 23:635-659.
Biagioli M, Pinto M, Cesselli D, Zaninello M, Lazarevic D, Roncaglia P, Simone R, Vlachouli C, Plessy C, Bertin N, Beltrami A, Kobayashi K, Gallo V, Santoro C, Ferrer I, Rivella S, Beltrami CA, Carninci P, Raviola E, Gustincich S (2009) Unexpected expression of alpha- and beta-globin in mesencephalic dopaminergic neurons and glial cells. Proc Natl Acad Sci U S A 106:15454-15459.
Brittain JM, Piekarz AD, Wang Y, Kondo T, Cummins TR, Khanna R (2009) An atypical role for collapsin response mediator protein 2 (CRMP-2) in neurotransmitter release via interaction with presynaptic voltage-gated calcium channels. J Biol Chem 284:31375-31390.
Brown N, Alkhayer K, Clements R, Singhal N, Gregory R, Azzam S, Li S, Freeman E, McDonough J (2016) Neuronal hemoglobin expression and its relevance to multiple sclerosis neuropathology. J Mol Neurosci 59:1-17.
Carter LM, Starkey ML, Akrimi SF, Davies M, McMahon SB, Bradbury EJ (2008) The Yellow Fluorescent Protein (YFP-H) mouse reveals neuroprotection as a novel mechanism underlying chondroitinase ABC-mediated repair after spinal cord injury. J Neurosci 28:14107-14120.
Dodt M, Roehr JT, Ahmed R, Dieterich C (2012) FLEXBAR-flexible barcode and adapter processing for next-generation sequencing platforms. Biology (Basel) 1:895-905.
Feng G, Mellor RH, Bernstein M, Keller-Peck C, Nguyen QT, Wallace M, Nerbonne JM, Lichtman JW, Sanes JR (2000) Imaging neuronal subsets in transgenic mice expressing multiple spectral variants of GFP. Neuron 28:41-51.
Freed J, Chakrabarti L (2016) Defining a role for hemoglobin in Parkinson’s disease. NPJ Parkinsons Dis 2:16021.
Ghosh A, Haiss F, Sydekum E, Schneider R, Gullo M, Wyss MT, Mueggler T, Baltes C, Rudin M, Weber B, Schwab ME (2010) Rewiring of hindlimb corticospinal neurons after spinal cord injury. Nat Neurosci 13:97-104.
Ghosh A, Peduzzi S, Snyder M, Schneider R, Starkey M, Schwab ME (2012) Heterogeneous spine loss in layer 5 cortical neurons after spinal cord injury. Cereb Cortex 22:1309-1317.
Hill RL, Zhang YP, Burke DA, DeVries WH, Zhang Y, Magnuson DSK, Whittemore SR, Shields CB (2009) Anatomical and functional outcomes following a precise, graded, dorsal laceration spinal cord injury in C57BL/6 Mice. J Neurotrauma 26:1.
Jaerve A, Kruse F, Malik K, Hartung HP, Müller HW (2012) Age-dependent modulation of cortical transcriptomes in spinal cord injury and repair. PLoS One 7:49812.
Jia Z, Zhu H, Li J, Wang X, Misra H, Li Y (2012) Oxidative stress in spinal cord injury and antioxidant-based intervention. Spinal Cord 50:264-274.
Jin X, Sasamoto K, Nagai J, Yamazaki Y, Saito K, Goshima Y, Inoue T, Ohshima T (2016) Phosphorylation of CRMP2 by Cdk5 regulates dendritic spine development of cortical neuron in the mouse hippocampus. Neural Plast 2016:6790743.
Jure I, Labombarda F (2017) Spinal cord injury drives chronic brain changes. Neural Regen Res 12:1044-1047.
Kawano Y, Yoshimura T, Tsuboi D, Kawabata S, Kaneko-Kawano T, Shirataki H, Takenawa T, Kaibuchi K (2005) CRMP-2 Is Involved in Kinesin-1-Dependent Transport of the Sra-1/WAVE1 Complex and Axon Formation. Mol Cell Biol 25:9920-9935.
Kim BG, Dai HN, McAtee M, Vicini S, Bregman BS (2006) Remodeling of synaptic structures in the motor cortex following spinal cord injury. Exp Neurol 198:401-415.
Kim D, Langmead B, Salzberg SL (2015) HISAT: A fast spliced aligner with low memory requirements. Nat Methods 12:357-360.
Kimura T, Arimura N, Fukata Y, Watanabe H, Iwamatsu A, Kaibuchi K (2005) Tubulin and CRMP-2 complex is transported via Kinesin-1. J Neurochem 93:1371-1382.
Kruse F, Bosse F, Vogelaar CF, Brazda N, Küry P, Gasis M, Müller H (2011) Cortical gene expression in spinal cord injury and repair: insight into the functional complexity of the neural regeneration program. Front Mol Neurosci 4:26.
Liberzon A, Birger C, Thorvaldsdóttir H, Ghandi M, Mesirov JP, Tamayo P (2015) The molecular signatures database hallmark gene set collection. Cell Syst 1:417-425.
Liberzon A, Subramanian A, Pinchback R, Thorvaldsdóttir H, Tamayo P, Mesirov JP (2011) Molecular signatures database (MSigDB) 3.0. Bioinformatics 27:1739-1740.
Lin PC, Chan PM, Hall C, Manser E (2011) Collapsin response mediator proteins (CRMPs) are a new class of microtubule-associated protein (MAP) that selectively interacts with assembled microtubules via a taxol-sensitive binding interaction. J Biol Chem 286:41466-41478.
Martins-de-Souza D, Cassoli JS, Nascimento JM, Hensley K, Guest PC, Pinzon-Velasco AM, Turck CW (2015) The protein interactome of collapsin response mediator protein-2 (CRMP2/DPYSL2) reveals novel partner proteins in brain tissue. Proteomics Clin Appl 9(9-10):817-831.
Nagai J, Kitamura Y, Owada K, Yamashita N, Takei K, Goshima Y, Ohshima T (2015) Crmp4 deletion promotes recovery from spinal cord injury by neuroprotection and limited scar formation. Sci Rep 5:8269.
Nagai J, Owada K, Kitamura Y, Goshima Y, Ohshima T (2016) Inhibition of CRMP2 phosphorylation repairs CNS by regulating neurotrophic and inhibitory responses. Exp Neurol 277:283-295.
Nakamura H, Takahashi-Jitsuki A, Makihara H, Asano T, Kimura Y, Nakabayashi J, Yamashita N, Kawamoto Y, Nakamura F, Ohshima T, Hirano H, Tanaka F, Goshima Y (2018) Proteome and behavioral alterations in phosphorylation-deficient mutant collapsin response mediator protein 2 knock-in mice. Neurochem Int 119:207-217.
Nakamura M, Okano H (2013) Cell transplantation therapies for spinal cord injury focusing on induced pluripotent stem cells. Cell Res 23:70-80.
Picelli S, Faridani OR, Björklund ÅK, Winberg G, Sagasser S, Sandberg R (2014) Full-length RNA-seq from single cells using Smart-seq2. Nat Protoc 9:171-181.
Rahajeng J, Giridharan SSP, Naslavsky N, Caplan S (2010) Collapsin response mediator protein-2 (Crmp2) regulates trafficking by linking endocytic regulatory proteins to dynein motors. J Biol Chem 285:31918-31922.
Richter F, Meurers BH, Zhu C, Medvedeva VP, Chesselet MF (2009) Neurons express hemoglobin α- and β-chains in rat and human brains. J Comp Neurol 515:538-547.
Schelshorn DW, Schneider A, Kuschinsky W, Weber D, Krüger C, Dittgen T, Bürgers HF, Sabouri F, Gassler N, Bach A, Maurer MH (2009) Expression of hemoglobin in rodent neurons. J Cereb Blood Flow Metab 29:585-595.
Schmidt EF, Strittmatter SM (2007) The CRMP family of proteins and their role in Sema3A signaling. Adv Exp Med Biol 600:1-11.
Schneider CA, Rasband WS, Eliceiri KW (2012) NIH Image to ImageJ: 25 years of image analysis. Nat Methods 9:671-675.
Wang LH, Strittmatter SM (1996) A family of rat CRMP genes is differentially expressed in the nervous system. J Neurosci 16:6197-6207.
Wang Y, Brittain JM, Wilson SM, Khanna R (2009) Emerging roles of collapsin response mediator proteins (CRMPs) as regulators of voltage-gated calcium channels and synaptic transmission. Commun Integr Biol 3:172-175.
Wu J, Stoica BA, Luo T, Sabirzhanov B, Zhao Z, Guanciale K, Nayar SK, Foss CA, Pomper MG, Faden AI (2014a) Isolated spinal cord contusion in rats induces chronic brain neuroinflammation, neurodegeneration, and cognitive impairment. Involvement of cell cycle activation. Cell Cycle 13:2446-2458.
Wu J, Zhao Z, Sabirzhanov B, Stoica BA, Kumar A, Luo T, Skovira J, Faden AI (2014b) Spinal cord injury causes brain inflammation associated with cognitive and affective changes: Role of cell cycle pathways. J Neurosci 34:10989-11006.
Xu W, Chi L, Xu R, Ke Y, Luo C, Cai J, Qiu M, Gozal D, Liu R (2005) Increased production of reactive oxygen species contributes to motor neuron death in a compression mouse model of spinal cord injury. Spinal Cord 43:204-213.
Yamashita N, Ohshima T, Nakamura F, Kolattukudy P, Honnorat J, Mikoshiba K, Goshima Y (2012) Phosphorylation of CRMP2 (collapsin response mediator protein 2) is involved in proper dendritic field organization. J Neurosci 32:1360-1365.
Yoda T, Hosokawa M, Takahashi K, Sakanashi C, Takeyama H, Kambara H (2017) Site-specific gene expression analysis using an automated tissue micro-dissection punching system. Sci Rep 7:4325.
Yu G, He QY (2016) ReactomePA: An R/Bioconductor package for reactome pathway analysis and visualization. Mol Biosyst 12:477-479.
Yu G, Wang LG, Han Y, He QY (2012) clusterProfiler: an R package for comparing biological themes among gene clusters. Omi A J Integr Biol 16:284-287.
Zhang K, Zhang J, Zhou Y, Chen C, Li W, Ma L, Zhang L, Zhao J, Gan W, Zhang L, Tang P (2015) Remodeling the dendritic spines in the hindlimb representation of the sensory cortex after spinal cord hemisection in mice. PLoS One 10:1-16.
Zheng Y, Sethi R, Mangala LS, Taylor C, Goldsmith J, Wang M, Masuda K, Karaminejadranjbar M, Mannion D, Miranda F, Herrero-Gonzalez S, Hellner K, Chen F, Alsaadi A, Albukhari A, Fotso DC, Yau C, Jiang D, Pradeep S, Rodriguez-Aguayo C, et al. (2018) Tuning microtubule dynamics to enhance cancer therapy by modulating FER-mediated CRMP2 phosphorylation. Nat Commun 9:476.
C-Editors: Zhao M, Li CH; T-Editor: Jia Y
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5]
[Table 1], [Table 2]