|Year : 2019 | Volume
| Issue : 9 | Page : 1583-1593
RNA sequencing screening of differentially expressed genes after spinal cord injury
Yi Li1, Ying Chen2, Xiang Li2, Jian Wu2, Jing-Ying Pan2, Ri-Xin Cai2, Ri-Yun Yang2, Xiao-Dong Wang MD 3
1 School of Biology & Basic Medical Sciences, Soochow University, Suzhou; Department of Histology and Embryology, Medical College, Nantong University, Nantong, Jiangsu Province, China
2 Department of Histology and Embryology, Medical College, Nantong University, Nantong, Jiangsu Province, China
3 Department of Histology and Embryology, Medical College, Nantong University; Jiangsu Key Laboratory of Neuroregeneration, Co-innovation Center of Neuroregeneration, Nantong University, Nantong, Jiangsu Province, China
|Date of Submission||05-Dec-2018|
|Date of Acceptance||01-Feb-2019|
|Date of Web Publication||9-May-2019|
Department of Histology and Embryology, Medical College, Nantong University; Jiangsu Key Laboratory of Neuroregeneration, Co-innovation Center of Neuroregeneration, Nantong University, Nantong, Jiangsu Province
Source of Support: This study was supported by the National Natural Science Foundation of China, No. 31570983 (to XDW); the Priority Academic
Program Development of Jiangsu Higher Education Institutes of China, Conflict of Interest: None
In the search for a therapeutic schedule for spinal cord injury, it is necessary to understand key genes and their corresponding regulatory networks involved in the spinal cord injury process. However, ad hoc selection and analysis of one or two genes cannot fully reveal the complex molecular biological mechanisms of spinal cord injury. The emergence of second-generation sequencing technology (RNA sequencing) has provided a better method. In this study, RNA sequencing technology was used to analyze differentially expressed genes at different time points after spinal cord injury in rat models established by contusion of the eighth thoracic segment. The numbers of genes that changed significantly were 944, 1362 and 1421 at 1, 4 and 7 days after spinal cord injury respectively. After gene ontology analysis and temporal expression analysis of the differentially expressed genes, C5ar1, Socs3 and CCL6 genes were then selected and identified by real-time polymerase chain reaction and western blot assay. The mRNA expression trends of C5ar1, Socs3 and CCL6 genes were consistent with the RNA sequencing results. Further verification and analysis of C5ar1 indicate that the level of protein expression of C5ar1 was consistent with its nucleic acid level after spinal cord injury. C5ar1 was mainly expressed in neurons and astrocytes. Finally, the gene Itgb2, which may be related to C5ar1, was found by Chilibot database and literature search. Immunofluorescence histochemical results showed that the expression of Itgb2 was highly consistent with that of C5ar1. Itgb2 was expressed in astrocytes. RNA sequencing technology can screen differentially expressed genes at different time points after spinal cord injury. Through analysis and verification, genes strongly associated with spinal cord injury can be screened. This can provide experimental data for further determining the molecular mechanism of spinal cord injury, and also provide possible targets for the treatment of spinal cord injury. This study was approved ethically by the Laboratory Animal Ethics Committee of Jiangsu Province, China (approval No. 2018-0306-001) on March 6, 2018.
Keywords: nerve regeneration; spinal cord injury; RNA sequencing; C5ar1; Itgb2; Socs3; CCL6; astrocytes; rats; neural regeneration
|How to cite this article:|
Li Y, Chen Y, Li X, Wu J, Pan JY, Cai RX, Yang RY, Wang XD. RNA sequencing screening of differentially expressed genes after spinal cord injury. Neural Regen Res 2019;14:1583-93
|How to cite this URL:|
Li Y, Chen Y, Li X, Wu J, Pan JY, Cai RX, Yang RY, Wang XD. RNA sequencing screening of differentially expressed genes after spinal cord injury. Neural Regen Res [serial online] 2019 [cited 2019 May 27];14:1583-93. Available from: http://www.nrronline.org/text.asp?2019/14/9/1583/255994
Yi Li, Ying Chen. These authors contributed equally to this work.
Chinese Library Classification No. R446; R394; R741
| Introduction|| |
Spinal cord injury (SCI) results in extremely high disability rate and serious consequences. Its pathophysiological changes are mainly divided into two stages: primary injury and secondary injury. The former refers to the damage caused by external force directly or indirectly to the spinal cord, including crush injury, contusion or transverse injury (Widrin, 2018). Although the primary injury is transient, it triggers more complex secondary injuries, including neuronal apoptosis, glial cell proliferation and inflammatory responses, which further aggravate and extend the severity of injuries (Noller et al., 2017). Eventually, the primary injury can lead to the loss of sensory and/or motor function below the plane of damage, as well as cause dysfunction of various organs and cause chronic neuropathic pain (Park et al., 2018). SCI not only reduces the quality of life of patients, but also brings great emotional and economic burdens to families and society. However, current therapies are not totally satisfactory, and there is an urgent need for more effective treatments.
A growing number of researchers believe that an understanding of the mechanisms underlying SCI and its secondary damage is an important route to developing a treatment strategy for SCI (Dalbayrak et al., 2015; Garcia et al., 2016; Wang et al., 2018). Therefore, analysis starting at the molecular level for the gene expression and its corresponding regulatory pathways in the SCI process is required. However, research in the past often directly selected individual genes for detection and analysis of regulatory pathways (Park et al., 2016; Alizadeh et al., 2017; Tapia et al., 2017; Deng et al., 2018). The emergence and development of gene sequencing technologies, such as cDNA chips and gene chip technology, have largely addressed this problem (Uchida et al., 2010; Chamankhah et al., 2013). However, initial sequencing technology is limited in its use owing to its low resolution and accuracy (Marioni et al., 2008). The emergence of second-generation sequencing technologies such as RNA sequencing (RNA-Seq) technology, with its unprecedented sensitivity and specificity, in spite of its high cost (Narrandes and Xu, 2018), allow for better genome-wide sequencing.
RNA-Seq technology, also known as whole transcriptome shotgun sequencing, is a genome-wide sequencing technology that increases the coverage of DNA sequences (Ozsolak and Milos, 2011; Van Verk et al., 2013; Yi et al., 2018). Second-generation high-throughput sequencing is performed by extracting all transcribed RNA in a biological sample and then reversely transcribing them into cDNAs. On this basis, the overlapping assembly of the fragments is carried out, so that transcripts can be obtained to observe post-transcriptional changes, gene fusion, and gene expression changes (Maher et al., 2009). A global understanding of the up-and-down gene expression in the biological sample is thus formed (Ingolia et al., 2012). In recent years, if the procedures are tightly controlled, the results of RNA-Seq are generally more accurate and reliable so it is widely used in the molecular mechanism studies of various complex diseases, such as tumors, brain microangiopathy, depression, amyotrophic lateral sclerosis and SCI (Duan et al., 2015; Shi et al., 2017; Barham et al., 2018; Herman et al., 2018; Le et al., 2018; Yuan et al., 2018; Wang et al., 2019).
To this end, we applied RNA-Seq technology to efficiently analyze the differences in gene expression at different stages after SCI, and verified and functionally analyzed the selected genes to find possible determinants in the damage process (Chen et al., 2013). Our findings provide data for further study of the molecular mechanisms of SCI, and provide a reference for developing treatment strategies for SCI.
| Materials and Methods|| |
A total of 49 specific-pathogen-free 8-week-old adult female Sprague-Dawley rats weighing 200–220 g were purchased from the Laboratory Animal Center of Nantong University in China (license No. SYXK (Su) 2017-0046). All the surgical interventions and postoperative animal care were performed according to the National Institutes of Health (USA) guidelines for care and use of laboratory animals (NIH Publication No. 85-23, revised 1996). This study was approved ethically by the Laboratory Animal Ethics Committee of Jiangsu Province, China (approval No. 2018-0306-001) on March 6, 2018.
SCI model establishment
Forty-nine adult female Sprague-Dawley rats were randomly divided into a sham group (n = 10) and an SCI group (n = 39). The rats in the SCI group were randomly and equally subdivided into three given time groups, 1, 4 and 7 days after SCI. Animal models were prepared in accordance with a previous method (Donnelly et al., 2012). Under anesthesia, a 1-cm incision was made at the eighth thoracic segment (T8) to open the lamina and expose the spinal cord layer by layer. In the sham group, no damage was inflicted on the spinal cord and the incision was sutured layer by layer. In the SCI group, after fixation of the spine, a spinal cord hitter (IH-0400, PSI, Fairfax, VA, USA) inflicted a 1.5 N impact to the exposed spinal cord, followed by layer by layer suturing. The model rats were postoperatively injected intraperitoneally with 5-mL lactated Ringer’s solution (Anhui BBCA Pharmaceutical Co., Hefei, China) to supplement the blood volume, and then penicillin (100,000 units) (Shandong Lukang Pharmaceutical Co., Jining, China) was given daily for 1 week, and the bladder was massaged twice daily for artificial urination until it was able to urinate spontaneously.
Postoperative rat behavior testing
Postoperative rat behavioral testing was performed using the Basso, Beattie & Bresnahan (BBB) locomotor rating scale (Basso et al., 1995). Preoperative training was performed as follows: the rats were placed in a round, open basin to walk freely, and kept in the central area as much as possible. The test time was 4 minutes. Its movement range, trunk movement and coordination of the knee, hip and ankle during walking were observed. At least two testers performed each test and the mean of scores taken.
RNA-Seq technology and bioinformatics analysis
After anesthesia, a 5-mm spinal cord segment, including the injured portion, was taken from each rat, and total RNA was extracted according to the conventional method. Agilent 2200 (Agilent Technologies, Santa Clara, CA, USA) was used to test the integrity of samples. A sample was used for the establishment of a transcriptome library if the RNA-integrity number was more than 8. Qualified samples were sent to RiboBio Corporation (Guangzhou, China) for RNA-Seq analysis. Gene ontology (GO) analysis and temporal expression analysis were performed. Genes related to SCI and post-injury repair were subsequently screened for biological properties. Genes related to SCI were found by Chilibot database (http://www.chilibot.net/) and literature searches.
Real-time polymerase chain reaction
The qualified RNA samples were used as templates to reverse-transcribe cDNA in real-time polymerase chain reaction (PCR) assay. Then cDNA (250 mg/L) 2 μL, Fast Start Universal SYBR Green Master 10 μL, Primer (sense) 0.6 μL, Primer (antisense) 0.6 μL, and quantified to 20 μL with ddH2O. Primer sequences are shown in [Table 1]. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was the internal control, and the untemplated wells were used as negative controls. The relative quantitative analysis used was the 2–ΔΔCt method.
Western blot assay
Spinal cord tissues were added to the lysate at a ratio of 10 mL/g for protein sample preparation. The protein concentration was determined by the bicinchoninic acid method. The samples were separated on a 10% sodium dodecyl sulfate polyacrylamide gel electrophoresis gel at 60–110 V for 2 hours, and electrotransferred onto polyvinylidene difluoride membranes at 300 mA for 90 minutes. The membranes were blocked with Tris-buffered saline with Tween 20 containing 5% skim milk powder for 2 hours, then incubated with the primary antibodies (rabbit-anti-C5ar1, 1:800, Proteintech Group Inc., Wuhan, China; rabbit-anti-GAPDH, 1:1000, Abcam, Cambridge, MA, USA) at 4°C overnight. The next day, they were rinsed, then incubated with the secondary antibody (horseradish peroxidase-labeled goat anti-rabbit IgG, 1:1000, Abcam) for 2 hours at room temperature. The protein bands werevisualized with enhanced chemiluminescence (Bio-Rad laboratories Inc., Hercules, CA, USA), and analyzed using a Chemi DocTMXRS+ gel imaging system (Bio-Rad laboratories Inc.), and Image Lab 5.1 software (Bio-Rad laboratories Inc.).
Under anesthesia, the animals were perfused with normal saline and 4% paraformaldehyde. The 5-mm spinal cord segments, including the injured portion, were excised. Frozen sections were then made, blocked with 0.5% Triton-X 100 and 3% bovine serum albumin, incubated with primary antibodies (rabbit anti-C5ar1, 1:400, Proteintech Group Inc.; rabbit anti-Itgb2, 1:400, Abcam; mouse anti-glial fibrillary acidic protein (GFAP), 1:500, Abcam; mouse anti-NF200, 1:500, Abcam; goat anti-GFAP, 1:500, Abcam), overnight at 4°C. After rinsing with phosphate buffer saline, the samples were incubated with secondary antibodies (Alexa Fluor® 488-labeled donkey anti-rabbit IgG antibody, 1:1000, Thermo Fisher, Waltham, MA, USA; Alexa Fluor® 568-labeled donkey anti-mouse IgG, 1:1000, Thermo Fisher; Alexa Fluor® 647-labeled donkey anti-goat IgG antibody, 1:1000, Thermo Fisher) at room temperature for 2 hours in the dark. The samples were rinsed with phosphate buffer saline and mounted with DAPI-containing sealant. They were observed using laser confocal microscope (SP7, Leica, Wetzlar, Germany) and recorded. No primary antibodies were used in controls.
All data are expressed as the mean ± SD. Both data analysis and image generation were performed by Graph Pad Prism 6 (GraphPad Software Co., San Diego, CA, USA). The BBB score data were analyzed by repeated-measure analysis of variance for homogeneity of analysis, and then analysis of variance was used for comparison between groups. PCR and western blot data were analyzed by one-way analysis of variance for homogeneity of variance, and then analysis of variance was used for comparison between groups. A value of P < 0.05 was considered statistically significant.
| Results|| |
Changes in locomotor function of SCI rats
The general state of SCI rats was confirmed based on their BBB score to ensure that the SCI animal model had been successfully prepared. Subsequent experimental procedures were then performed. Compared with the sham group, the BBB score was significantly lower at each time point after SCI. The hindlimb movement was completely lost at 1 day after injury, and the BBB score was 0. The score on the hindlimb function was still low at 4 days after injury, and slightly restored at 7 days after injury [Figure 1]. This indicates that the hindlimb function of SCI rats had some but limited self-repair within 1 week after SCI. These results indicated that the animal model of SCI had been prepared successfully.
|Figure 1: Basso, Beattie & Bresnahan locomotor rating scale for the hindlimb function in rats.|
Data are expressed as the mean ± SD (repeated-measure analysis of variance). *P < 0.05, vs, sham group. All experiments were conducted five times. SCI: Spinal cord injury.
Click here to view
Screening of differentially expressed genes after SCI
A large number of differentially expressed genes were screened by high-throughput sequencing at the genome-wide level using RNA-Seq technology. To further understand the specific distribution of differential transcripts, comparisons between the sham group and experiment groups were performed using Cuffdiff. The results of the Venn plot showed that compared with the sham group, the number of genes that changed significantly at 1, 4 and 7 days after SCI were 944, 1362 and 1421, respectively, of which 652 genes were differentially expressed at both 1 and 4 days, 648 genes at both 1 and 7 days, and 1170 genes at both 4 and 7 days [Figure 2].
|Figure 2: Venn plots of differentially expressed genes in the tissues from the sham group at 1, 4, and 7 days after spinal cord injury.|
Differentially expressed genes were screened by high-throughput sequencing at the genome-wide level using RNA-Seq technology.
Click here to view
GO functional analysis of differentially expressed genes was performed to further explore the function of the differentially expressed genes. The results showed that among the top 20 groups of GO enrichments, multiple biological processes related to damage and repair were involved, such as GO: 0006954 (inflammatory response), GO: 0050729 (positive regulation of inflammatory response), GO: 0009612 (response to mechanical stimulus); and GO: 0045766 (positive regulation of angiogenesis) [Figure 3].
|Figure 3: GO functional analysis of differentially expressed genes.|
The abscissa represents various terms, and the ordinate represents numbers of genes.
Click here to view
Temporal expression analysis of differentially expressed genes after SCI
Temporal expression analysis of differentially expressed genes at different time points was performed to screen for genes with a consistently significant change (fold change > 4) from 0 to 1 day, 1 to 4 days, and 4 to 7 days [Table 2]. Searching the Chilibot database, the genes C5ar1, Socs3, and CCL6 found most strongly associated with SCI or spinal cord development were further screened.
|Table 2: Genes closely linked to spinal cord injury or spinal cord development, which showed significant changes to expression after spinal cord injury|
Click here to view
Expression of some genes after SCI verified by real-time PCR
Real-time PCR results verified the changes of the above three genes at different stages after SCI. C5ar1 expression was upregulated 1 day after injury, peaked at 4 days and downregulated significantly at 7 days after SCI compared with the sham group (P < 0.05 or P < 0.01; [Figure 4]). The variations in Socs3 and CCL6 expression were the similar to that of C5ar1 between SCI and sham groups (P < 0.05 or P < 0.01; [Figure 4]).
|Figure 4: Relative gene expression of C5ar1, Socs3 and CCL6 in the spinal cord after SCI.|
Data (2–ΔΔCt ratio to GAPDH mRNA) are expressed as the mean ± SD (n = 10 in sham group, n = 13 in SCI groups at 1, 4, and 7 days) relative to sham group (mean = 1). *P < 0.05, **P < 0.01, vs. sham group (one-way analysis of variance). All experiments were conducted five times. Sham: Sham group; 1 d, 4 d and 7 d: 1, 4, 7 days after SCI; SCI: spinal cord injury.
Click here to view
Expression of C5ar1 protein in the spinal cord after SCI
We selected the gene C5ar1 from the verified differentially expressed genes because it is involved in the response to injury, inflammatory response and neuronal apoptosis. The expression trend of C5ar1 after SCI by western blot assay was observed. C5ar1 expression significantly increased after SCI compared with that in the sham group (P < 0.05; [Figure 5]). C5ar1 expression peaked at 4 days after SCI. This finding was consistent with the RNA-Seq sequencing and real-time PCR results.
|Figure 5: Relative protein expression of C5ar1 in the spinal cord after SCI.|
Data are expressed as the mean ± SD (n = 10 in sham group, n = 13 in SCI groups at 1, 4, and 7 days), showing C5ar1 expression by the optical density ratio of C5ar1 to its GAPDH in each group, *P < 0.05, vs. sham group (one-way analysis of variance). All experiments were conducted five times. Sham: Sham group; 1 d, 4 d and 7 d: 1, 4, 7 days after SCI; SCI: spinal cord injury.
Click here to view
Bioinformatics analysis of C5ar1 in the spinal cord after SCI
Further bioinformatics analysis was performed on C5ar1. Itgb2, CD55 and TNF were most strongly associated with C5ar1 by searching the Chilibot database and query literature. These genes were also found to be differentially expressed temporally in the RNA-Seq results. Itgb2 expression was upregulated after 1 day of SCI, peaked at 4 days and downregulated at 7 days. TNF expression was upregulated at 1 day after SCI and downregulated at 4 days, but upregulated again at 7 days. CD55 expression was upregulated at 1 day after SCI and downregulated at 4 days until unchanged at 7 days. Therefore, only the expression trends of Itgb2 and C5ar1 were highly consistent (references are shown in [Table 3]).
Localization and expression of C5ar1 and Itgb2 in spinal cord tissue after SCI
The cellular localization and expression of C5ar1 in spinal cord tissue were observed by immunofluorescence histochemistry. The astrocyte and neuronal markers GFAP and NF200 were used for staining and localization (Trojanowski et al., 1986; Middeldorp and Hol, 2011; van Bodegraven et al., 2019). C5ar1 was mainly found in astrocytes and neurons, and its expression level was differentially expressed with time after SCI. C5ar1 was expressed in the cells at 1 day after SCI; both the immunofluorescence expression and the number of positive cells increased at 4 days. At 7 days, the immunofluorescence expression of cells decreased slightly [Figure 6] & [Figure 7]). This result was highly consistent with the changing trends of results from RNA-Seq, real-time PCR and western blot assay.
|Figure 6: Immunofluorescence staining of C5ar1 and GFAP in spinal cord tissue of the SCI rats.|
Green fluorescence (DyLight 488) in A, E, I and M refers to C5ar1; red fluorescence (DyLight 568) in B, F, J and N refers to GFAP; blue fluorescence in C, G, K and O refers to DAPI labeled nuclei. D, H, L and P are merged figures. Arrows indicate double-positive cells. Sham: Sham group; 1 day, 4 days and 7 days: 1, 4, 7 days after SCI. Scale bars: 30 μm. Scale bars in the insets of H, L, and P: 10 μm. GFAP: Glial fibrillary acidic protein; SCI: spinal cord injury.
Click here to view
|Figure 7: Immunofluorescence staining of C5ar1 and NF200 in spinal cord tissue of the SCI rats.|
Green fluorescence (DyLight 488) in A, E, I and M refers to C5ar1; red fluorescence (DyLight 568) in B, F, J and N refers to NF200; blue fluorescence in C, G, K and O refers to DAPI labeled nuclei; and D, H, L and P are merged figures. Sham: Sham group; 1 day, 4 days and 7 days: 1, 4, 7 days after SCI. Arrows indicate double-positive cells. Scale bar: 30 μm. Scale bars in the insets of H, L, and P: 10 μm. NF200: Neurofilament-200; SCI: spinal cord injury.
Click here to view
The expression of Itgb2, which was highly consistent with the changes in C5ar1 expression, was mainly expressed in astrocytes of the spinal cord tissue after SCI [Figure 8].
|Figure 8: Immunofluorescence staining of Itgb2, GFAP and NF200 in spinal cord tissue of the SCI rats.|
Green fluorescence (DyLight 488) in A, E, I and M refers to Itgb2; cyan fluorescence (DyLight 647) in B, F, J and N refers to GFAP; red fluorescence (DyLight 568) in C, G, K and O refers to NF200; and D, H, L and P are merged figures. Sham: Sham group; 1 day, 4 days and 7 days: 1, 4, 7 days after SCI. Arrows indicate double-positive cells. Scale bars: 75 μm. Scale bars in the insets of H, L, and P: 10 μm. GFAP: Glial fibrillary acidic protein; NF200: neurofilament-200; SCI: spinal cord injury.
Click here to view
| Discussion|| |
In recent years, researchers have focused on analyzing the complex pathological processes of SCI and its molecular mechanisms to determine the effects of changes to the microenvironment on SCI repair. This approach attempted to find techniques that might protect nerve tissue cells by improving the microenvironment around SCI and enhancing axonal regeneration, such as drug therapy, genetic modification, cell transplantation or tissue engineering. The RNA-Seq technology used in this study investigated the differences in transcriptome levels at different stages after SCI. The results of this study showed that compared with the sham group, the number of differentially expressed genes at 1, 4 and 7 days after SCI were 944, 1362 and 1421, respectively. Of the differentially expressed genes those genes of interest were confirmed after screening by GO functional analysis and temporal expression analysis. Multiple biological processes associated with injury repair were enriched, suggesting that these genes may play an important role in SCI and post-injury repair processes. The differentially expressed transcripts between samples were selected by the difference multiples (log2|fold_change|>2) and significant levels (q value < 0.05). Subsequent n searching of the Chilibot database, selected those genes with persistent and significant changes (fold change > 4) that were strongly associated with SCI, some of which were verified by real-time PCR. The results showed that the differential expression of C5ar1, Socs3 and CCL6 at different stages after SCI was consistent with the results of RNA-Seq sequencing.
Socs3 is a feedback inhibitor of the JAK/STAT pathway and plays a key role in regulating the immune response of JAK/STAT-mediated inflammatory cytokines (Croker et al., 2008). The structural basis for this effect is that Socs3 contains a kinase inhibitory region that in combination with the JAK/STAT/gp130 complex can inhibit downstream signaling and inflammatory cytokines (La Manna et al., 2018). In the central nervous system, Socs3 is expressed on almost all types of cells (Baker et al., 2009). A previous study has shown that the Socs3 protein is involved in the regulation and progression of microglial immune responses and, therefore, regulation of its expression may be a target for drug development in the control of inflammation (Cianciulli et al., 2017). After SCI, Socs3 was found to not only regulate the inflammatory response, but also to be associated with the formation of glial scars. For example, Socs3 may prevent the formation of glial scars (Li et al., 2010), and can regulate the proliferation and differentiation of NG2 cells in glial scars (Hackett et al., 2016). It had also been reported that increased Socs3 level was strongly associated with neuronal death (Park et al., 2014). The Socs3-dependent JAK/STAT pathway played a key role in neuronal loss and even acted in regulating axon regeneration during secondary injury of SCI (Liu et al., 2015). Multiple genes in zebrafish after SCI, equivalent to mammalian stat3, socs3, atf3, mmp9 and sox11 were upregulated (Hui et al., 2014). CCL6 is a member of the CC chemokine family and is normally expressed on microglia cells without stimulation. Under normal physiological environment or abnormal pathological conditions, CCL6 mediates the main immune effector cells of the central nervous system, or interaction between microglia and astrocytes, by activating the PI3-kinase/Akt pathway (Kanno et al., 2005). RNA-Seq screening indicated that both Socs3 and CCL6 were related to pathophysiological reactions after SCI, especially inflammatory reactions. They could be used as potential targets to reduce inflammatory responses as well as to improve the repair environment after SCI.
C5ar1 (CD88) is the most important receptor for the complement component C5a, which is structurally a typical G-protein coupled receptor that signals primarily through Gαi and Gα16 (Monk et al., 2007). C5ar1 is widely expressed in inflammatory cells, such as granulocytes, monocytes and macrophages, as well as in the central nervous system, especially in glial cells (Woodruff et al., 2010). Immunofluorescence histochemical staining of this study showed that C5ar1 was mainly present in astrocytes and neurons, which is consistent with previous findings. The widespread expression of C5ar1 means that it may be involved in various activities. First, damage to the central nervous system activates the complement system, triggering a cascade of biological effects after C5a binds to C5ar1, including recruitment of granulocytes, activation of glial cell and neuronal signal transduction pathways (Guo et al., 2013b). C5ar1 and C5a are also effectors of the neuropathic complement pathway and have multi-effect inflammatory effects, which can stimulate nuclear factor-kappa β, activator protein-1, and adenosine monophosphate to participate in protein signal transduction. Activation of C5ar1 is associated with the onset and severity of epileptic seizures and subsequent neuronal degeneration (Benson et al., 2015). Pathological changes in the hSOD1G93A mouse model of amyotrophic lateral sclerosis can be improved when antagonists inhibit C5a-C5ar1 receptor signaling (Lee et al., 2017). C5a-C5ar1 is also associated with permeability of the blood-brain barrier (Lee et al., 2017) and cognitive loss of Alzheimer’s disease. After SCI, C5ar1 is highly expressed and specifically binds to C5a, inducing a chemotactic response that participates in neuronal apoptosis and subsequent inflammatory reactions. This process is inextricably linked to a series of secondary injuries (Kwon et al., 2004; Qiao et al., 2010; Thundyil et al., 2012; Guo et al., 2013a). It is generally believed that the expression level of C5ar1 in the central nervous system is simultaneously associated with the inflammatory state. Its expression produces a variety of inflammatory mediators during complement activation, such as interleukin-2, interleukin-6, and tumor necrosis factor. The GO analysis of this study showed that C5ar1 was involved in multiple biological processes, such as injury response, inflammatory reaction and neuronal apoptosis. C5ar1 was found to have temporal differential expression after SCI, as shown by the changes in immunofluorescence histochemical staining, real-time PCR and western blot assay. C5ar1 expression was upregulated at 1 day after SCI, peaked at 4 days, and downregulated at 7 days, but was still higher than that in the sham group. This trend was consistent with the RNA-Seq results. This suggests that C5ar1 begins to be expressed in the early stage after SCI and then increases its own protein level, providing a material basis for mediating secondary pathological changes following SCI.
Among the differential genes identified above, we were most interested in C5ar1 (CD88). Therefore, we performed further analysis on C5ar1. By searching the Chilibot database and searching relevant literature, we found that many genes, such as Itgb2, CD55 and TNF, may be associated with C5ar1, and the expression trends of Itgb2 and C5ar1 are highly consistent. Itgb2 expression was upregulated at 1 day after SCI, reached the highest peak at 4 days then decreased at 7 days. Using immunofluorescence histochemistry, Itgb2, localized to astrocytes in the spinal cord sections, also exhibited temporal differential expression after SCI. Itgb2 (CD18) belongs to a family of β2 integrins. It is a transmembrane glycoprotein distributed on the cell surface that plays an important role in cell adhesion, cell surface signaling transduction and immune response. In the models of cerebral hemorrhage injuries, the symptoms of acute brain injury and cerebral edema can be reduced and even lower mortality rates in CD18-deficient mice compared with wild mice (Titova et al., 2008a, 2008b). As reported by Gris et al. (2004), the use of monoclonal antibodies to block Itgb2 significantly reduced secondary injury after SCI and promoted the recovery of motor and sensory functions in SCI rats. These findings indicate that, after SCI, an increase in Itgb2 level will aggravate secondary damage. However, further investigations on the relationship between C5ar1 and Itgb2, as well as the expression of CD55 and TNF and their roles and relationship with C5ar1, are required after SCI.
In conclusion, the differentially expressed genes, which play important roles in a series of pathological reactions after SCI and post-injury repair process, could be screened by RNA-Seq technology, and subjected to bioinformatic analysis. C5ar1 and Itgb2, two prominent genes of particular interest, were selected from all the differential genes, and subjected to various experimental techniques to verify changes in their expression and analyze their functions. On the one hand, it provides a basis for further elucidation of the molecular mechanism of SCI; on the other hand, it suggests targets for the treatment of SCI. Although several differentially expressed genes were found, only a few were studied and confirmed as participants in SCI. A more comprehensive study may be conducted in the future.
Author contributions: Study design and experimental implementation: YL, YC, XDW; experimental implementation: XL, JW, JYP, RXC, RYY; fundraising: XDW. All authors approved the final version of the paper.
Conflicts of interest: There are no conflicts of interest.
Financial support: This study was supported by the National Natural Science Foundation of China, No. 31570983 (to XDW); the Priority Academic Program Development of Jiangsu Higher Education Institutes of China. The funders did not participate in data collection and analysis, article writing or submission.
Institutional review board statement: The studies were approved by the Lab Animal Ethical Committee of Jiangsu Province, China (approval No. 2018-0306-001) on March 6, 2018.
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.
Open peer reviewer: Li Yao, Wichita State University, USA.
Additional file: Open peer review report 1[Additional file 1].
Funding: This study was supported by the National Natural Science Foundation of China, No. 31570983 (to XDW); the Priority Academic Program Development of Jiangsu Higher Education Institutes of China.
| References|| |
Abe T, Hosur KB, Hajishengallis E, Reis ES, Ricklin D, Lambris JD, Hajishengallis G (2012) Local complement-targeted intervention in periodontitis: proof-of-concept using a C5a receptor (CD88) antagonist. J Immunol 189:5442-5448.
Alizadeh A, Dyck SM, Kataria H, Shahriary GM, Nguyen DH, Santhosh KT, Karimi-Abdolrezaee S (2017) Neuregulin-1 positively modulates glial response and improves neurological recovery following traumatic spinal cord injury. Glia 65:1152-1175.
Amara U, Kalbitz M, Perl M, Flierl MA, Rittirsch D, Weiss M, Schneider M, Gebhard F, Huber-Lang M (2010) Early expression changes of complement regulatory proteins and C5A receptor (CD88) on leukocytes after multiple injury in humans. Shock 33:568-575.
Baker BJ, Akhtar LN, Benveniste EN (2009) SOCS1 and SOCS3 in the control of CNS immunity. Trends Immunol 30:392-400.
Barham C, Fil D, Byrum SD, Rahmatallah Y, Glazko G, Kiaei M (2018) RNA-Seq analysis of spinal cord tissues from hPFN1(G118V) transgenic mouse model of ALS at pre-symptomatic and end-stages of disease. Sci Rep 8:13737.
Basso DM, Beattie MS, Bresnahan JC (1995) A sensitive and reliable locomotor rating scale for open field testing in rats. J Neurotrauma 12:1-21.
Benson MJ, Thomas NK, Talwar S, Hodson MP, Lynch JW, Woodruff TM, Borges K (2015) A novel anticonvulsant mechanism via inhibition of complement receptor C5ar1 in murine epilepsy models. Neurobiol Dis 76:87-97.
Chamankhah M, Eftekharpour E, Karimi-Abdolrezaee S, Boutros PC, San-Marina S, Fehlings MG (2013) Genome-wide gene expression profiling of stress response in a spinal cord clip compression injury model. BMC Genomics 14:583.
Chen K, Deng S, Lu H, Zheng Y, Yang G, Kim D, Cao Q, Wu JQ (2013) RNA-seq characterization of spinal cord injury transcriptome in acute/subacute phases: a resource for understanding the pathology at the systems level. PLoS One 8:e72567.
Cianciulli A, Calvello R, Porro C, Trotta T, Panaro MA (2017) Understanding the role of SOCS signaling in neurodegenerative diseases: Current and emerging concepts. Cytokine Growth Factor Rev 37:67-79.
Croker BA, Kiu H, Nicholson SE (2008) SOCS regulation of the JAK/STAT signalling pathway. Semin Cell Dev Biol 19:414-422.
Dalbayrak S, Yaman O, Yilmaz T (2015) Current and future surgery strategies for spinal cord injuries. World J Orthop 6:34-41.
Deng G, Gao Y, Cen Z, He J, Cao B, Zeng G, Zong S (2018) miR-136-5p regulates the inflammatory response by targeting the IKKbeta/NF-kappaB/A20 pathway after spinal cord injury. Cell Physiol Biochem 50:512-524.
Donnelly EM, Madigan NN, Rooney GE, Knight A, Chen B, Ball B, Kinnavane L, Garcia Y, Dockery P, Fraher J, Strappe PM, Windebank AJ, O’Brien T, McMahon SS (2012) Lentiviral vector delivery of short hairpin RNA to NG2 and neurotrophin-3 promotes locomotor recovery in injured rat spinal cord. Cytotherapy 14:1235-1244.
Duan H, Ge W, Zhang A, Xi Y, Chen Z, Luo D, Cheng Y, Fan KS, Horvath S, Sofroniew MV, Cheng L, Yang Z, Sun YE, Li X (2015) Transcriptome analyses reveal molecular mechanisms underlying functional recovery after spinal cord injury. Proc Natl Acad Sci U S A 112:13360-13365.
Garcia E, Aguilar-Cevallos J, Silva-Garcia R, Ibarra A (2016) Cytokine and growth factor activation in vivo and in vitro after spinal cord injury. Mediators Inflamm 2016:9476020.
Girke G, Kohl B, Busch C, John T, Godkin O, Ertel W, Schulze-Tanzil G (2014) Tenocyte activation and regulation of complement factors in response to in vitro cell injury. Mol Immunol 60:14-22.
Gris D, Marsh DR, Oatway MA, Chen Y, Hamilton EF, Dekaban GA, Weaver LC (2004) Transient blockade of the CD11d/CD18 integrin reduces secondary damage after spinal cord injury, improving sensory, autonomic, and motor function. J Neurosci 24:4043-4051.
Guo Q, Cheng J, Zhang H, Zhang J, Su B, Bian C, Lin S (2013a) Expressions of C5a and its receptor CD88 after spinal cord injury in C3-deficient mice. Scand J Immunol 77:224-229.
Guo Q, Cheng J, Zhang J, Su B, Bian C, Lin S, Zhong C (2013b) Delayed post-injury administration of C5a improves regeneration and functional recovery after spinal cord injury in mice. Clin Exp Immunol 174:318-325.
Hackett AR, Lee DH, Dawood A, Rodriguez M, Funk L, Tsoulfas P, Lee JK (2016) STAT3 and SOCS3 regulate NG2 cell proliferation and differentiation after contusive spinal cord injury. Neurobiol Dis 89:10-22.
Herman PE, Papatheodorou A, Bryant SA, Waterbury CKM, Herdy JR, Arcese AA, Buxbaum JD, Smith JJ, Morgan JR, Bloom O (2018) Highly conserved molecular pathways, including Wnt signaling, promote functional recovery from spinal cord injury in lampreys. Sci Rep 8:742.
Hui SP, Sengupta D, Lee SG, Sen T, Kundu S, Mathavan S, Ghosh S (2014) Genome wide expression profiling during spinal cord regeneration identifies comprehensive cellular responses in zebrafish. PLoS One 9:e84212.
Ingolia NT, Brar GA, Rouskin S, McGeachy AM, Weissman JS (2012) The ribosome profiling strategy for monitoring translation in vivo by deep sequencing of ribosome-protected mRNA fragments. Nat Protoc 7:1534-1550.
Kanno M, Suzuki S, Fujiwara T, Yokoyama A, Sakamoto A, Takahashi H, Imai Y, Tanaka J (2005) Functional expression of CCL6 by rat microglia: a possible role of CCL6 in cell-cell communication. J Neuroimmunol 167:72-80.
Kwon BK, Tetzlaff W, Grauer JN, Beiner J, Vaccaro AR (2004) Pathophysiology and pharmacologic treatment of acute spinal cord injury. Spine J 4:451-464.
La Manna S, Lee E, Ouzounova M, Di Natale C, Novellino E, Merlino A, Korkaya H, Marasco D (2018) Mimetics of suppressor of cytokine signaling 3: Novel potential therapeutics in triple breast cancer. Int J Cancer 143:2177-2186.
Le TT, Savitz J, Suzuki H, Misaki M, Teague TK, White BC, Marino JH, Wiley G, Gaffney PM, Drevets WC, McKinney BA, Bodurka J (2018) Identification and replication of RNA-Seq gene network modules associated with depression severity. Translational psychiatry 8:180.
Lee JD, Kumar V, Fung JN, Ruitenberg MJ, Noakes PG, Woodruff TM (2017) Pharmacological inhibition of complement C5a-C5a1 receptor signalling ameliorates disease pathology in the hSOD1G93A mouse model of amyotrophic lateral sclerosis. Br J Pharmacol 174:689-699.
Li LM, Li JB, Zhu Y, Fan GY (2010) Soluble complement receptor type 1 inhibits complement system activation and improves motor function in acute spinal cord injury. Spinal Cord 48:105-111.
Liu X, Williams PR, He Z (2015) SOCS3: a common target for neuronal protection and axon regeneration after spinal cord injury. Exp Neurol 263:364-367.
Maher CA, Kumar-Sinha C, Cao X, Kalyana-Sundaram S, Han B, Jing X, Sam L, Barrette T, Palanisamy N, Chinnaiyan AM (2009) Transcriptome sequencing to detect gene fusions in cancer. Nature 458:97-101.
Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y (2008) RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res 18:1509-1517.
Middeldorp J, Hol EM (2011) GFAP in health and disease. Prog Neurobiol 93:421-443.
Monk PN, Scola AM, Madala P, Fairlie DP (2007) Function, structure and therapeutic potential of complement C5a receptors. Br J Pharmacol 152:429-448.
Mueller M, Herzog C, Larmann J, Schmitz M, Hilfiker-Kleiner D, Gessner JE, Theilmeier G (2013) The receptor for activated complement factor 5 (C5aR) conveys myocardial ischemic damage by mediating neutrophil transmigration. Immunobiology 218:1131-1138.
Narrandes S, Xu W (2018) Gene expression detection assay for cancer clinical use. J Cancer 9:2249-2265.
Noller CM, Groah SL, Nash MS (2017) Inflammatory stress effects on health and function after spinal cord injury. Top Spinal Cord Inj Rehabil 23:207-217.
Ozsolak F, Milos PM (2011) RNA sequencing: advances, challenges and opportunities. Nat Rev Genet 12:87-98.
Park J, Decker JT, Smith DR, Cummings BJ, Anderson AJ, Shea LD (2018) Reducing inflammation through delivery of lentivirus encoding for anti-inflammatory cytokines attenuates neuropathic pain after spinal cord injury. J Control Release 290:88-101.
Park KW, Lin CY, Lee YS (2014) Expression of suppressor of cytokine signaling-3 (SOCS3) and its role in neuronal death after complete spinal cord injury. Exp Neurol 261:65-75.
Park KW, Lin CY, Benveniste EN, Lee YS (2016) Mitochondrial STAT3 is negatively regulated by SOCS3 and upregulated after spinal cord injury. Exp Neurol 284:98-105.
Qiao F, Atkinson C, Kindy MS, Shunmugavel A, Morgan BP, Song H, Tomlinson S (2010) The alternative and terminal pathways of complement mediate post-traumatic spinal cord inflammation and injury. Am J Pathol 177:3061-3070.
Shi LL, Zhang N, Xie XM, Chen YJ, Wang R, Shen L, Zhou JS, Hu JG, Lu HZ (2017) Transcriptome profile of rat genes in injured spinal cord at different stages by RNA-sequencing. BMC Genomics 18:173.
Sopp P, Werling D, Baldwin C (2007) Cross-reactivity of mAbs to human CD antigens with cells from cattle. Vet Immunol Immunopathol 119:106-114.
Tapia VS, Herrera-Rojas M, Larrain J (2017) JAK-STAT pathway activation in response to spinal cord injury in regenerative and non-regenerative stages of Xenopus laevis. Regeneration (Oxf) 4:21-35.
Thundyil J, Pavlovski D, Hsieh YH, Gelderblom M, Magnus T, Fairlie DP, Arumugam TV (2012) C5a receptor (CD88) inhibition improves hypothermia-induced neuroprotection in an in vitro ischemic model. Neuromolecular Med 14:30-39.
Titova E, Kevil CG, Ostrowski RP, Rojas H, Liu S, Zhang JH, Tang J (2008a) Deficiency of CD18 gene reduces brain edema in experimental intracerebral hemorrhage in mice. Acta Neurochir Suppl 105:85-87.
Titova E, Ostrowski RP, Kevil CG, Tong W, Rojas H, Sowers LC, Zhang JH, Tang J (2008b) Reduced brain injury in CD18-deficient mice after experimental intracerebral hemorrhage. J Neurosci Res 86:3240-3245.
Trojanowski JQ, Walkenstein N, Lee VM (1986) Expression of neurofilament subunits in neurons of the central and peripheral nervous system: an immunohistochemical study with monoclonal antibodies. J Neurosci 6:650-660.
Uchida K, Nakajima H, Hirai T, Yayama T, Chen KB, Kobayashi S, Roberts S, Johnson WE, Baba H (2010) Microarray analysis of expression of cell death-associated genes in rat spinal cord cells exposed to cyclic tensile stresses in vitro. BMC Neurosci 11:84.
van Bodegraven EJ, van Asperen JV, Robe PAJ, Hol EM (2019) Importance of GFAP isoform-specific analyses in astrocytoma. Glia doi: 10.1002/glia.23594.
Van Verk MC, Hickman R, Pieterse CM, Van Wees SC (2013) RNA-Seq: revelation of the messengers. Trends Plant Sci 18:175-179.
Wang J, Dean DC, Hornicek FJ, Shi H, Duan Z (2019) RNA sequencing (RNA-Seq) and its application in ovarian cancer. Gynecol Oncol 152:194-201.
Widrin C (2018) Scalp acupuncture for the treatment of motor function in acute spinal cord injury: a case report. J Acupunct Meridian Stud 11:74-76.
Woodruff TM, Ager RR, Tenner AJ, Noakes PG, Taylor SM (2010) The role of the complement system and the activation fragment C5a in the central nervous system. Neuromolecular Med 12:179-192.
Yi S, Wang XH, Xing LY (2018) Transcriptome analysis of adherens junction pathway-related genes after peripheral nerve injury. Neural Regen Res 13:1804-1810.
Yuan X, Wu Q, Liu X, Zhang H, Xiu R (2018) Transcriptomic profile analysis of brain microvascular pericytes in spontaneously hypertensive rats by RNA-Seq. Am J Transl Res 10:2372-2386.
P-Reviewer: Yao L; C-Editor: Zhao M; S-Editors: Yu J, Li CH; L-Editors: Dawes EA, Haase R, Qiu Y, Song LP; T-Editor: Jia Y
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7], [Figure 8]
[Table 1], [Table 2], [Table 3]