|Year : 2020 | Volume
| Issue : 3 | Page : 503-511
Protein microarray analysis of cytokine expression changes in distal stumps after sciatic nerve transection
Xiao-Qing Cheng1, Xue-Zhen Liang2, Shuai Wei1, Xiao Ding1, Gong-Hai Han1, Ping Liu1, Xun Sun1, Qi Quan1, He Tang1, Qing Zhao1, Ai-Jia Shang3, Jiang Peng4
1 Beijing Key Laboratory of Regenerative Medicine in Orthopedics, Institute of Orthopedics, Chinese PLA General Hospital, Beijing, China
2 Beijing Key Laboratory of Regenerative Medicine in Orthopedics, Institute of Orthopedics, Chinese PLA General Hospital, Beijing; The First Clinical Medical College, Shandong University of Traditional Chinese Medicine, Jinan, Shandong Province, China
3 Co-innovation Center of Neuroregeneration, Nantong University, Nantong, Jiangsu Province; Department of Neurosurgery, Chinese PLA General Hospital, Beijing, China
4 Beijing Key Laboratory of Regenerative Medicine in Orthopedics, Institute of Orthopedics, Chinese PLA General Hospital, Beijing; Co-innovation Center of Neuroregeneration, Nantong University, Nantong, Jiangsu Province, China
|Date of Submission||05-Dec-2018|
|Date of Acceptance||08-Mar-2019|
|Date of Web Publication||26-Sep-2019|
Co-innovation Center of Neuroregeneration, Nantong University, Nantong, Jiangsu Province; Department of Neurosurgery, Chinese PLA General Hospital, Beijing
Beijing Key Laboratory of Regenerative Medicine in Orthopedics, Institute of Orthopedics, Chinese PLA General Hospital, Beijing; Co-innovation Center of Neuroregeneration, Nantong University, Nantong, Jiangsu Province
Source of Support: This study was supported by the National Key Research & Development Program of China, No. 2017YFA0104702 (to AJS), and the
National Basic Research Program of China (973 Program), No. 2014CB542201 (to JP), Conflict of Interest: None
A large number of chemokines, cytokines, other trophic factors and the extracellular matrix molecules form a favorable microenvironment for peripheral nerve regeneration. This microenvironment is one of the major factors for regenerative success. Therefore, it is important to investigate the key molecules and regulators affecting nerve regeneration after peripheral nerve injury. However, the identities of specific cytokines at various time points after sciatic nerve injury have not been determined. The study was performed by transecting the sciatic nerve to establish a model of peripheral nerve injury and to analyze, by protein microarray, the expression of different cytokines in the distal nerve after injury. Results showed a large number of cytokines were up-regulated at different time points post injury and several cytokines, e.g., ciliary neurotrophic factor, were downregulated. The construction of a protein-protein interaction network was used to screen how the proteins interacted with differentially expressed cytokines. Kyoto Encyclopedia of Genes and Genomes pathway and Gene ontology analyses indicated that the differentially expressed cytokines were significantly associated with chemokine signaling pathways, Janus kinase/signal transducers and activators of transcription, phosphoinositide 3-kinase/protein kinase B, and notch signaling pathway. The cytokines involved in inflammation, immune response and cell chemotaxis were up-regulated initially and the cytokines involved in neuronal apoptotic processes, cell-cell adhesion, and cell proliferation were up-regulated at 28 days after injury. Western blot analysis showed that the expression and changes of hepatocyte growth factor, glial cell line-derived neurotrophic factor and ciliary neurotrophic factor were consistent with the results of protein microarray analysis. The results provide a comprehensive understanding of changes in cytokine expression and changes in these cytokines and classical signaling pathways and biological functions during Wallerian degeneration, as well as a basis for potential treatments of peripheral nerve injury. The study was approved by the Institutional Animal Care and Use Committee of the Chinese PLA General Hospital, China (approval number: 2016-x9-07) in September 2016.
Keywords: cytokines; distal stump; gene ontology; Kyoto Encyclopedia of Genes and Genomes pathway; peripheral nerve injury; protein microarray; protein-protein interaction network; Wallerian degeneration
|How to cite this article:|
Cheng XQ, Liang XZ, Wei S, Ding X, Han GH, Liu P, Sun X, Quan Q, Tang H, Zhao Q, Shang AJ, Peng J. Protein microarray analysis of cytokine expression changes in distal stumps after sciatic nerve transection. Neural Regen Res 2020;15:503-11
|How to cite this URL:|
Cheng XQ, Liang XZ, Wei S, Ding X, Han GH, Liu P, Sun X, Quan Q, Tang H, Zhao Q, Shang AJ, Peng J. Protein microarray analysis of cytokine expression changes in distal stumps after sciatic nerve transection. Neural Regen Res [serial online] 2020 [cited 2020 Feb 17];15:503-11. Available from: http://www.nrronline.org/text.asp?2020/15/3/503/266062
Xiao-Qing Cheng, Xue-Zhen Liang
These authors contributed equally to this work.
Chinese Library Classification No. R447; R363; R741
| Introduction|| |
Although the peripheral nerve has remarkable regenerative abilities after injury, it is difficult to recover completely from long-term nerve defects following peripheral nerve injury (PNI) (Raimondo et al., 2011). After PNI, Wallerian degeneration, a complicated process involving distal axonal degeneration and myelin breakdown, takes place immediately (Geuna et al., 2009). Subsequently, multiple macrophages and monocytes migrate to remove the axon and myelin debris. Schwann cells proliferate to form bands of Büngner as a bridge to the defective peripheral nerve system. They secrete a large number of chemokines, cytokines and other trophic factors and extracellular matrix molecules, which form a favorable microenvironment for peripheral nerve regeneration (Frostick et al., 1998; Chen et al., 2007). This microenvironment is one of the major factors of regenerative success (Webber and Zochodne, 2010). Therefore, it is important to investigate the key molecules and regulators affecting nerve regeneration after PNI.
In early studies, microarrays were widely used to characterize differentially expressed genes in the distal nerve stump following PNI (Pan et al., 2017). Temporal expression patterns of upregulated and downregulated genes during Wallerian degen-eration were validated (Yu et al., 2016; Yi et al., 2017). RNA-sequencing was also performed to evaluate comprehensive transcriptomic expression and identify numerous differentially expressed microRNAs at different time points during peripheral nerve regeneration (Yi et al., 2015). These studies have helped elucidate global gene expression changes involved in peripheral nerve repair. However, in many cases during peripheral nerve regeneration, multiple cellular events occur after a gene is translated into a protein. Several studies have suggested that some proteins have significant effects on many biological processes during peripheral nerve repair and Wallerian degeneration, including immune responses, macrophage recruitment and Schwann cell reprogramming (van der Laan et al., 1996; da Costa et al., 1997; Siebert et al., 2000; Takahashi et al., 2007; Clements et al., 2017). Nevertheless, the integrated relationships between proteins involved in PNI and recovery are not yet clear.
In this study, we aimed to investigate the dynamic differential expression of 67 cytokines using a protein micro-array to achieve greater insight into the relative pathways or networks during Wallerian degeneration of injured sciatic nerve. We used bioinformatic analyses (protein-protein interaction (PPI) network, Kyoto Encyclopedia of Genes and Genomes Pathway (KEGG) pathway and gene ontology (GO) analysis).
| Materials and Methods|| |
In this study, eighty male Sprague-Dawley rats (specific pathogen free level) weighing 200–250 g, aged 7–8 weeks, were obtained from the Experimental Animal Research Center at the Chinese PLA General Hospital, China. The rats were housed in a temperature-controlled environment and fed water and food ad libitum. The Institutional Animal Care and Use Committee of the Chinese PLA General Hospital, China approved the experimental procedures (approval No. 2016-x9-07) in September 2016.
Rats were randomly divided into five groups: control and 4 periods post injury. Briefly, rats were anesthetized by intraperitoneal injection of pentobarbital (30 mg/kg; Sigma-Aldrich, St. Louis, MO, USA). The sciatic nerve was exposed via incision on the lateral aspect of the right hind limb and was excised through the middle site of the exposed sciatic nerve (Yu et al., 2012). At 1, 7, 14, and 28 days after nerve transection, the rats were sacrificed by decapitation. The distal nerve stumps were removed and stored at −80°C until further use. As a control, rats in the 0 day group underwent sham surgery of the right sciatic nerves.
Protein microarray analysis
Protein samples were extracted from the distal nerve stumps of sciatic nerves of the rats and lysed in radioimmuno-precipitation assay buffer lysis buffer containing a protease inhibitor cocktail (Pulilai, Beijing, China). Subsequently, protein concentrations were determined using the Pierce Bicinchoninic acid assay Kit (Thermo Fisher, Waltham, MA, USA).
Microarray analysis was performed using the Rat Cytokine Array 67 (RayBiotech, Guangzhou, Guangdong Province, China), as described previously (Luck et al., 2017). Briefly, the microarray was incubated with 100 μL sample diluent at room temperature for 30 minutes to block the slides. After removing the diluent, the array was completely covered with 100 μL of the protein sample (2 mg/mL) and incubated at 4°C overnight. Biotinylated antibody cocktail (80 μL) was added to each well at room temperature for 2 hours. Cy3 equivalent dye-conjugated streptavidin (80 μL) was added to each well, and the wells were covered with aluminum foil and placed in a dark room to avoid light exposure at room temperature for 1 hour. The signals were visualized using the LuxScan 10K scanner (CapitaBio, Beijing, China) at 532 nm wavelength.
The expression levels of proteins at 1, 7, 14, and 28 days after sciatic nerve transection were compared with those in the control group. Proteins with an expression fold change > 2 or < −2 and adjusted P-value < 0.05 were considered significantly differentially expressed. The differentially expressed proteins, Venn diagram and Principal Component analysis were performed by R packages (http://www.R-project.org) named gplots, VennDiagram and scatterplot3d, respectively.
The GeneMANIA database (http://genemania.org) is a gene and protein analysis tool designed to predict PPIs. The various proteins were mapped using GeneMANIA to evaluate the interactive relationships among them. The PPI networks were then constructed using Cytoscape software (http://www.cytoscape.org).
Proteins selected from the PPI networks were systematically analyzed using the Database for Annotation, Visu-alization and Integrated Discovery (DAVID; http://www.david.niaid.nih.gov) program to identify significantly enriched KEGG pathways and GO categories (Dennis et al., 2003). The KEGG pathways and GO categories were performed by R packages named ggplot2.
Western blot analysis
Sciatic nerves were rinsed in cold PBS and lysed on ice in radioimmunoprecipitation assay buffer containing a pro-tease inhibitor cocktail (Pulilai), and the resulting tissue lysates were mixed with sample buffer and boiled at 95°C for 5 minutes. Equal amounts of protein from each sample were subjected to 10–15% sodium dodecyl sulfate-polyacrylamide gel electrophoresis and transferred to polyvinylidene fluoride membranes (Pulilai). The membranes were blocked in 5% nonfat dry milk at 4°C for 1 hour and incubated with rabbit anti-hepatocyte growth factor (HGF) antibody (1:1000; ab83760, Abcam, Cambridge, UK), rabbit anti-glial cell line-derived neurotrophic factor (GDNF) antibody (1:1000; ab18956, Abcam), rabbit anti-ciliary neurotrophic factor (CNTF) antibody (1:1000; ab46172, Abcam) or rabbit anti-β-actin antibody (1:1000; Proteintech, Chicago, IL, USA) at 4°C overnight. These were followed by the appropriate secondary antibody, donkey-anti-rabbit-HRP (1:5000, Pulilai) at room temperature for 1 hour. The membranes were developed using an enhanced chemiluminescence substrate (Thermo Fisher). Measurement of the protein band intensities was conducted using ImageJ software (National Institutes of Health, Bethesda, MD, USA).
All data are expressed as the mean ± standard error of the mean (SEM). Significant differences among data were determined by one-way analysis of variance, followed by Tukey’s post hoc test. All analyses were performed using GraphPad Prism 7.0 (GraphPad Software, Inc., San Diego, CA, USA) and a P-value < 0.05 was designated as indicating statistical significance.
| Results|| |
Differentially expressed proteins in the distal nerve stump following sciatic nerve transection
To gain a better understanding of the microenvironment of the distal nerve stump, we examined the expression patterns of 67 proteins in injured sciatic nerves at different time points using the Rat Cytokine Array 67 [Table 1]. Proteins with a fold change in expression > 2 or < −2 and an adjusted P < 0.05 were defined as differentially expressed. A full list of the differentially expressed proteins at each time point evaluated is displayed in a heatmap [Figure 1]A, [Figure 1]B, [Figure 1]C, [Figure 1]D. Compared with the control group (0 day after PNI), the expression of nearly 20% of the total proteins increased at 1 day after PNI, including some chemokines and interleukins (ILs). The number of upregulated proteins increased to 33% of the total proteins at 7 days after PNI. In addition to chemokines and ILs, Notch 1/2 and Neuropilin 1/2, which are related to some classical pathways, increased in expression starting at 7 days after PNI. At 14 days after PNI, no ILs, except for IL-22, showed any remarkable expression changes. At 28 days after PNI, the number of upregulated proteins had decreased slightly. From day 1 to day 28 after PNI, only a few proteins, including CNTF, were downregulated. The numbers of overlapping differentially expressed proteins at 1, 7, 14, 28 days post injury are displayed in a Venn diagram [Figure 1]E. Principal component analysis based on the expression values of the 67 proteins demonstrated five completely independent clusters: control and the other time points after injury [Figure 1]F.
|Table 1: The 67 cytokine proteins included on the RayBiotech Rat Cytokine Array 67|
Click here to view
|Figure 1: Overview of differentially expressed cytokines at the distal stump following sciatic nerve transection.|
(A–D) Heatmap and hierarchical clustering at 1 (A), 7 (B), 14 (C), 28 (D) days. Protein expression was increased (red) and decreased (green) compared with the control (0 day after sciatic nerve transection). (E) Venn diagram displaying the number of different proteins among the four groups of sciatic nerves at each time point after sciatic nerve transection. (F) Principal component analysis of the differentially expressed cytokines at different time points following sciatic nerve transection. PC: Principal component; PCA: principal component analysis.
Click here to view
Protein-protein interaction network of the differentially expressed proteins
Many proteins were temporally differentially expressed in the injured sciatic nerve, suggesting that those cytokines have a dramatic effect on the nerve microenvironment after PNI. To explore this further, PPI networks were constructed by uploading the up- and downregulated proteins into GeneMANIA. In addition to the differentially expressed proteins, other chemotactic and inflammatory factors were selected within these networks, such as CC chemokine ligand (CCL17), chemokine (C-X-C motif) ligand (CXCL13), the tissue inhibitor of metalloproteinases-3 (TIMP3) and others. At 1 day after PNI, the differentially expressed proteins were mainly related to chemotactic family proteins (CCL2, CCL3, CCL7, CCL17, CXCL2, CXCL3, CXCL6, CXCL9, CXCL11 and CXCL13) [Figure 2]A. At 7 days after PNI, the proteins in the PPI included GFR alpha2, GFR alpha3, GFR alpha4, IL9, IL19, IL20, tyrosine-protein kinase 2 (JAK 2), Jag1, Jag2, TIMP3, TIMP4, VEGFB, GDNF, Neurturin (NRTN) [Figure 2]B. At 14 days after PNI, CCL1, CCL2, CCL5, CCL7, CCL15, CXCL1, CXCL2, intercellular adhesion molecule 2 (ICAM2), ICAM4, ICAM5, GFR alpha1, GDNF family receptor (GFR) alpha2, GFR alpha4, Notch-3, Galectin-1, Galectin-2, Galectin-5, Galectin-12 were selected in the PPI network [Figure 2]C. At 28 days after PNI, GFR alpha2, GFR alpha3, GFR alpha4, Notch-3, VEGFB, GDNF, NRTN, placenta growth factor, TIMP2, TIMP3, TIMP4, ICAM1, ICAM2, ICAM4, ICAM5, Galectin-12, GTPase Hras, Jag2 interacted with the differentially expressed proteins [Figure 2]D. In addition, there were several connections between the differentially expressed proteins and the selected proteins. Co-expression characteristics and described physical interactions accounted for most of the aforementioned targets and their interacting proteins. Other results, including shared protein domains, co-localization and predictions are shown in [Figure 2].
|Figure 2: Networks of protein-protein interactions of the differentially expressed proteins in the distal nerve stump following sciatic nerve transection.|
(A–D) Differentially expressed proteins at 1 (A), 7 (B), 14 (C), 28 (D) days. The sizes of the spots in the network represent the weights of the interaction involving the protein. The red spots represent the upregulated proteins and green spots the downregulated proteins at each time point after sciatic nerve transection. The blue spots represent the selected proteins interact with the differentially expressed proteins. The colors of the the connecting lines represent the different interactions between protein and protein. The thickness of the connecting lines represents the score of each connection.
Click here to view
Kyoto Encyclopedia of Genes and Genomes enriched analysis of integrated metabolic pathways following sciatic nerve injury
The proteins from the PPI networks were correlated with integrated metabolic pathways using KEGG analysis in DAVID. The top 10 canonical pathways associated with differentially expressed cytokines at each time point were displayed and analyzed using a P-value threshold of 0.05 [Figure 3]. Cytokine-cytokine receptor interactions and rheumatiod arthritis were activated throughout the entire post-injury period. Notably, the chemokine signaling pathways, Janus kinase/signal transducers and activators of transcription (JAK/STAT) and phosphoinositide 3-kinase/protein kinase B (PI3K/Akt), and leishmainasis were drastically stimulated at 1–7 days after PNI, while the notch signaling pathway and dorso-ventral axis formation were enriched at 14–28 days after PNI. Besides, malaria and influenze A were activited at 1–14 days after PNI. A full list of the canonical pathways and their involved molecules identified at 1, 7, 14, and 28 days after PNI is provided in Additional Table 1.
|Figure 3: Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enriched among the proteins in the protein-protein interaction network in the distal nerve stump following sciatic nerve transection.|
(A–D) Differentially expressed proteins at 1 (A), 7 (B), 14 (C) and 28 (D) days. The top 10 KEGG pathways with P < 0.05 are listed. The x-axis represents the gene ratio, defined as the ratio of the numbers of differentially expressed proteins to all proteins annotated in the KEGG pathway using Database for Annotation, Visualization, and Integrated Discovery (DAVID; http://www.david.niaid.nih.gov).
Click here to view
Gene ontology analysis of biological processes following sciatic nerve injury
To further explore the effects of cytokines following sciatic nerve injury, GO analysis was performed in DAVID to analyze biological processes. The top 10 biological processes (P < 0.05) associated with these proteins are displayed in [Figure 4]. At 1–7 days after PNI, biological processes related to inflammation, the immune response, and cell chemotaxis were activated. At 14 days after PNI, even more cell chemotaxis processes, as well as cell response to IL-1, interferon-gamma were activated, including neutrophil, monocyte, lymphocyte, and eosinophil chemotaxis. At 28 days after PNI, neuronal apoptotic processes, cell-cell adhesion, and cell proliferation were enriched. However, some biological processes were negatively regulated at different time after PNI. For example, the negative regulation of apoptosis first increased at 1–7 days after PNI, and then progressively declined at 14 days after PNI, before increasing again at 28 days after PNI. Membrane protein ectodomain proteolysis was inhibited from 7–28 days after PNI. Catalytic activity reduced at 7 and 28 days after PNI. A full list of the biological processes and involved molecules at 1, 7, 14, and 28 days after PNI is provided in Additional Table 2. More information on the highly enriched cellular components associated with PNI associated with the external side of the plasma membrane, extracellular matrix, immunological synapse, and receptor complex are shown in Additional Table 3.
|Figure 4: Gene ontology (GO) biological processes of those proteins enriched in the protein-protein interaction network in the distal nerve stump following sciatic nerve transection relative to all proteins.|
(A–D) Differentially expressed proteins at 1 (A), 7 (B), 14 (C), 28 (D) days. The top 10 GO biological processes with P < 0.05 are listed. The x-axis represents the gene ratio, defined as the ratio of the numbers of differentially expressed proteins to all proteins annotated in the GO biological processes using Database for Annotation, Visualization, and Integrated Discovery (DAVID; http://www.david.niaid.nih.gov).
Click here to view
The relationships between proteins and functional terms in the distal nerve stump following sciatic nerve transection
The identified cytokines seem to affect many functional pathways and PPI networks via direct and indirect interactions. Therefore, to visualize smaller high-dimensional data subsets, the relationships between functional terms and differentially expressed proteins were assessed using the GOplot package in R (Walter et al., 2015) to integrate the expression data with the results of these analyses. Monocyte chemoattractant protein-1 (MCP-1/CCL2) regulated upon activation of normal T cell expressed and secreted (RANTES)/CCL5 and Galectin-3 that are related to monocyte and neutrophil chemotaxis, whereas CNTF and GDNF are specific to neuronal apoptotic processes. Notch1, Notch2, and Galectin-3 may participate in cell differentiation during nerve regeneration [Figure 5]. Additionally, inflammatory responses involve many cytokines, including CCL2, CCL17, CXCL2, CXCL9, CXCL11, CXCL13, NGFR, GDNF, IL10, IL20, and IL19 [Figure 5].
|Figure 5: The relationship between the differentially expressed proteins and functional terms (biological processes and canonical pathways) in the distal nerve stump following sciatic nerve transection.|
Multiple biological processes and canonical pathways are labeled by different colors in circles on the right. The proteins linked to each term are listed in circles on the left.
Click here to view
Western blot verification of the protein array results in the distal nerve stump following sciatic nerve transection
Western blot analysis was performed to validate the temporal expression patterns of the differentially expressed growth factors in the nerves identified in the microarray analyses. Western blot analysis confirmed that CNTF expression levels decreased markedly, whereas GDNF and HGF were increased after PNI [Figure 6]A & [Figure 6]B, similar to the results of microarray analyses [Figure 6]C.
|Figure 6: Western blot analyses of the differentially expressed growth factors in the distal nerve stump following sciatic nerve transection.|
(A) Protein bands of HGF, CNTF and GDNF were measured by Western blot analysis. β-Actin was used as the reference protein. (B) Quantitative result of HGF, CNTF and GDNF expression. (C) The fold changes of HGF, CNTF and GFR alpha-1 were measured by protein microarray. Data represent the mean ± SEM (n = 4) from three independent experiments, statistically analyzed by one-way analysis of variance followed by Tukey’s post hoc test. *P < 0.05, **P < 0.01, vs. control group (ctrl; 0 day after nerve injury). CNTF: Ciliary neurotrophic factor; GDNF: glial cell line-derived neurotrophic factor; HGF: hepatocyte growth factor.
Click here to view
| Discussion|| |
Wallerian demyelination is the most typical cause of PNI, resulting in a series of complicated cellular responses and molecular mechanisms. Numerous studies have demonstrated that PNI induces cytokine production in immune and non-immune cells at sites distal to the nerve lesion (Karanth et al., 2006; Kiguchi et al., 2017). Such cytokines are closely related to Wallerian demyelination and participate in peripheral nerve regeneration (Rotshenker, 2011; He et al., 2016; Lin et al., 2019). In the current study protein microarray and bioinformatic analyses were performed to examine the detailed kinetic changes in cytokine production in distal nerve stumps following sciatic nerve injury.
In this study, the 67 cytokines on the protein microarray comprise growth factors, chemotaxis factors, and other proteins. By screening differentially expressed cytokines after PNI, we discovered that some growth factors may be critical for sciatic nerve injury and regeneration. For instance, previous studies reported that GFR alpha-1 is expressed in myelinated peripheral nerves and the neuromuscular junction, exerting its effects on motor neurons by interacting with GDNF (Hase et al., 1999; Rosich et al., 2017). In line with previous findings, we detected upregulation of GFR alpha-1 protein at 7, 14, and 28 days after PNI. Additionally, HGF also increased 1.5-fold at 28 days post injury. HGF has been shown to promote the migration and proliferation of Schwann cells and increase the expression of neurotrophic factors and inflammatory cytokines such as GDNF and tumor necrosis factor-α (Ko et al., 2018). In contrast, we found that CNTF was downregulated over most of the post-injury period studied. The PPI network suggested that CNTF regulates neuronal apoptosis via the JAK/STAT signaling pathway. All of these results highlight the central roles of growth factors in nerve regeneration.
Chemokine factors have been identified as important modulators of peripheral nerve regeneration (Taskinen and Röyttä, 2000). MCP-1/CCL2, a CCL family member, is expressed at low levels under basal conditions and is upregulated rapidly and markedly in Schwann cells and neurons (Schreiber et al., 2001; Tanaka et al., 2004; Niemi et al., 2013). In this study, we found that MCP-1 expression increased 10-fold from 1 to 28 days after sciatic nerve injury. Consequently, the injury-induced increase in MCP-1 may lead to recruitment of inflammatory monocytes and macrophages to nerves via the Toll-like receptor-4 or STAT3-dependent signaling pathways (Niemi et al., 2013). Moreover, other CCL family members, including RANTES/CCL5, cutaneous T-cell-attracting chemokine (CTACK)/CCL27 were also upregulated at different time points after PNI. RANTES/CCL5 is another important chemokine that exhibits strong chemoattractant activity towards monocytes and leukocytes, inducing immune cell migration and protection of neurons, either directly or indirectly (Raport et al., 1996; Tillie-Leblond et al., 2000; Tokami et al., 2013; Solga et al., 2015). CTACK was reported to accelerate skin regeneration via specific chemokine-receptor interactions (Inokuma et al., 2006), suggesting an important role of CTACK in nerve regeneration.
Notch proteins (Notch-1–4) are transmembrane receptors that regulate cellular processes, including cell proliferation, apoptosis, and angiogenesis (Bolós et al., 2007; Fortini, 2009; Kopan and Ilagan, 2009). During nerve regeneration, notch signaling mediates the differentiation of adipose-derived stem cells into Schwann-like cells and of monocytes into macrophages (Ohishi et al., 2001; Kingham et al., 2009). In our study, both Notch-1 and -2 were upregulated and may play important roles in the differentiation of the nervous system. It is also worth noting that Neuropilin-1 was upregulated from 1 to 28 days and Neuropilin-2 within 7 days after PNI. These proteins are closely related to axonal guidance, angiogenesis, and motor neuron migration. Triggering receptor expressed on myeloid cells 1 (TREM1) is another protein involved in immune response signaling pathways (Collins et al., 2009; Kingham et al., 2009) and activation of TREM1 may increase the secretion of pro-inflammatory cytokines (Walter, 2016). TREM1 continually increased over the time period we examined and was related to neutrophil chemotaxis. These results suggest that a cytokine network is involved in the kinetics of macrophage recruitment and nerve removal of damaged nerves.
Bioinformatic analyses are efficient methods for interpreting proteomic or genomic information. PPI network bioinformatic data are used to predict protein functionality within sequence homology clusters (Athanasios et al., 2017). Chemotactic factors, immune factors, and other cytokines were identified in our PPI network. Thymus and activation-regulated chemokine (TARC)/CCL17 accelerates fibroblast migration and wound healing (Kato et al., 2011). Rac GTPase, together with mammalian T-cell lymphoma invasion, metastasis factor 1 (TIAM-1) and CDC42, has been shown to mediate axon guidance (Demarco et al., 2012). These proteins may participate in nerve regeneration directly or indirectly. Additionally, the KEGG pathways resulted in the cytokine-cytokine receptor interactions keeping active during Wallerian degeneration. Other enriched pathways included the JAK/STAT, PI3K/Akt, and chemokine signaling pathways. Furthermore, our study demonstrated a diverse array of biological processes, including chemotaxis, inflammatory and immune responses, cell migration, cell proliferation, apoptosis, and angiogenesis, that were significantly activated in the distal nerve stump following sciatic nerve transection. Based on the relationships between cytokines and our bioinformatic data, further in-depth studies are needed to determine the detailed mechanism of peripheral nerve regeneration.
In summary, we analyzed global changes in cytokine expression patterns at the distal nerve stump following PNI using protein microarray analysis. Although our results did not elucidate the mechanism of PNI, bioinformatic analysis enabled us to gain a comprehensive view of cytokine expression changes with time. They also show the relationships of these cytokines with canonical pathways, biological functions, and networks during Wallerian degeneration. Overall, our study may help identify the potential clinic treatments for PNI.
Author contributions: Study design: XQC, XZL, AJS, JP; experiments performing: XQC, GHH, SW, XD, PL, XS, QQ, HT; data analysis: XQC, XZL; manuscript writing: XQC; manuscript revising: QZ, AJS, JP. All authors approved the final version of the paper.
Conflicts of interest: The authors declare that there is no conflict of interests regarding the publication of this paper.
Financial support: This study was supported by the National Key Research & Development Program of China, No. 2017YFA0104702 (to AJS), and the National Basic Research Program of China (973 Program), No. 2014CB542201 (to JP). The funders had no involvement in the study design; data collection, analysis, and interpretation; paper writing; or decision to submit the paper for publication.
Institutional review board statement: This study was approved by the Institutional Animal Care and Use Committee of the Chinese PLA General Hospital (approval No. 2016-x9-07) in September 2016.
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 Table 1: All canonical pathways and involved molecules at 1, 7, 14, and 28 days after sciatic nerve injury[Additional file 1].
Additional Table 2: All biological function categories and involved molecules at 1, 7, 14, and 28 days after sciatic nerve injury[Additional file 2].
Additional Table 3: All cellular component categories and involved molecules at 1, 7, 14, and 28 days after sciatic nerve injury[Additional file 3].
Funding: This study was supported by the National Key Research & Development Program of China, No. 2017YFA0104702 (to AJS), and the National Basic Research Program of China (973 Program), No. 2014CB542201 (to JP).
| References|| |
Athanasios A, Charalampos V, Vasileios T, Ashraf GM (2017) Protein-protein interaction (PPI) network: recent advances in drug discovery. Curr Drug Metab 18:5-10.
Bolós V, Grego-Bessa J, de la Pompa JL (2007) Notch signaling in development and cancer. Endocr Rev 28:339-363.
Chen ZL, Yu WM, Strickland S (2007) Peripheral regeneration. Annu Rev Neurosci 30:209-233.
Clements MP, Byrne E, Camarillo Guerrero LF, Cattin AL, Zakka L, Ashraf A, Burden JJ, Khadayate S, Lloyd AC, Marguerat S, Parrinello S (2017) The wound microenvironment reprograms Schwann cells to invasive mesenchymal-like cells to drive peripheral nerve regeneration. Neuron 96:98-114.e7.
Collins CE, La DT, Yang HT, Massin F, Gibot S, Faure G, Stohl W (2009) Elevated synovial expression of triggering receptor expressed on myeloid cells 1 in patients with septic arthritis or rheumatoid arthritis. Ann Rheum Dis 68:1768-1774.
da Costa CC, van der Laan LJ, Dijkstra CD, Bruck W (1997) The role of the mouse macrophage scavenger receptor in myelin phagocytosis. Eur J Neurosci 9:2650-2657.
Demarco RS, Struckhoff EC, Lundquist EA (2012) The Rac GTP exchange factor TIAM-1 acts with CDC-42 and the guidance receptor UNC-40/DCC in neuronal protrusion and axon guidance. PLoS Genet 8:e1002665.
Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA (2003) DAVID: Database for annotation, visualization, and integrated discovery. Genome Biol 4:P3.
Fortini ME (2009) Notch signaling: the core pathway and its posttranslational regulation. Dev Cell 16:633-647.
Frostick SP, Yin Q, Kemp GJ (1998) Schwann cells, neurotrophic factors, and peripheral nerve regeneration. Microsurgery 18:397-405.
Geuna S, Raimondo S, Ronchi G, Di Scipio F, Tos P, Czaja K, Fornaro M (2009) Chapter 3: Histology of the peripheral nerve and changes occurring during nerve regeneration. Int Rev Neurobiol 87:27-46.
Hase A, Suzuki H, Arahata K, Akazawa C (1999) Expression of human GFR alpha-1 (GDNF receptor) at the neuromuscular junction and myelinated nerves. Neurosci Lett 269:55-57.
He XZ, Wang W, Hu TM, Ma JJ, Yu CY, Gao YF, Cheng XL, Wang P (2016) Peripheral nerve repair: theory and technology application. Zhongguo Zuzhi Gongcheng Yanjiu 20:1044-1050.
Inokuma D, Abe R, Fujita Y, Sasaki M, Shibaki A, Nakamura H, McMillan JR, Shimizu T, Shimizu H (2006) CTACK/CCL27 accelerates skin regeneration via accumulation of bone marrow-derived keratinocytes. Stem Cells 24:2810-2816.
Karanth S, Yang G, Yeh J, Richardson PM (2006) Nature of signals that initiate the immune response during Wallerian degeneration of peripheral nerves. Exp Neurol 202:161-166.
Kato T, Saeki H, Tsunemi Y, Shibata S, Tamaki K, Sato S (2011) Thymus and activation-regulated chemokine (TARC)/CC chemokine ligand (CCL) 17 accelerates wound healing by enhancing fibroblast migration. Exp Dermatol 20:669-674.
Kiguchi N, Kobayashi D, Saika F, Matsuzaki S, Kishioka S (2017) Pharmacological regulation of neuropathic pain driven by inflammatory macrophages. Int J Mol Sci 18:E2296.
Kingham PJ, Mantovani C, Terenghi G (2009) Notch independent signalling mediates Schwann cell-like differentiation of adipose derived stem cells. Neurosci Lett 467:164-168.
Ko KR, Lee J, Lee D, Nho B, Kim S (2018) Hepatocyte growth factor (HGF) promotes peripheral nerve regeneration by activating repair schwann cells. Sci Rep 8:8316.
Kopan R, Ilagan MX (2009) The canonical Notch signaling pathway: unfolding the activation mechanism. Cell 137:216-233.
Lin YF, Xie Z, Zhou J, Yin G, Lin HD (2019) Differential gene and protein expression between rat tibial nerve and common peroneal nerve during Wallerian degeneration. Neural Regen Res 14:2183-2191.
Luck C, DeMarco VG, Mahmood A, Gavini MP, Pulakat L (2017) Differential regulation of cardiac function and intracardiac cytokines by rapamycin in healthy and diabetic rats. Oxid Med Cell Longev 2017:5724046.
Niemi JP, DeFrancesco-Lisowitz A, Roldan-Hernandez L, Lindborg JA, Mandell D, Zigmond RE (2013) A critical role for macrophages near axotomized neuronal cell bodies in stimulating nerve regeneration. J Neurosci 33:16236-16248.
Ohishi K, Varnum-Finney B, Serda RE, Anasetti C, Bernstein ID (2001) The Notch ligand, Delta-1, inhibits the differentiation of monocytes into macrophages but permits their differentiation into dendritic cells. Blood 98:1402-1407.
Pan B, Shi ZJ, Yan JY, Li JH, Feng SQ (2017) Long non-coding RNA NONMMUG014387 promotes Schwann cell proliferation after peripheral nerve injury. Neural Regen Res 12:2084-2091.
Raimondo S, Fornaro M, Tos P, Battiston B, Giacobini-Robecchi MG, Geuna S (2011) Perspectives in regeneration and tissue engineering of peripheral nerves. Ann Anat 193:334-340.
Raport CJ, Gosling J, Schweickart VL, Gray PW, Charo IF (1996) Molecular cloning and functional characterization of a novel human CC chemokine receptor (CCR5) for RANTES, MIP-1beta, and MIP-1alpha. J Biol Chem 271:17161-17166.
Rosich K, Hanna BF, Ibrahim RK, Hellenbrand DJ, Hanna A (2017) The effects of glial cell line-derived neurotrophic factor after spinal cord injury. J Neurotrauma 34:3311-3325.
Rotshenker S (2011) Wallerian degeneration: the innate-immune response to traumatic nerve injury. J Neuroinflammation 8:109.
Schreiber RC, Krivacic K, Kirby B, Vaccariello SA, Wei T, Ransohoff RM, Zigmond RE (2001) Monocyte chemoattractant protein (MCP)-1 is rapidly expressed by sympathetic ganglion neurons following axonal injury. Neuroreport 12:601-606.
Siebert H, Sachse A, Kuziel WA, Maeda N, Bruck W (2000) The chemokine receptor CCR2 is involved in macrophage recruitment to the injured peripheral nervous system. J Neuroimmunol 110:177-185.
Solga AC, Pong WW, Kim KY, Cimino PJ, Toonen JA, Walker J, Wylie T, Magrini V, Griffith M, Griffith OL, Ly A, Ellisman MH, Mardis ER, Gutmann DH (2015) RNA sequencing of tumor-associated microglia reveals Ccl5 as a stromal chemokine critical for neurofibromatosis-1 glioma growth. Neoplasia 17:776-788.
Takahashi K, Prinz M, Stagi M, Chechneva O, Neumann H (2007) TREM2-transduced myeloid precursors mediate nervous tissue debris clearance and facilitate recovery in an animal model of multiple sclerosis. PLoS Med 4:e124.
Tanaka T, Minami M, Nakagawa T, Satoh M (2004) Enhanced production of monocyte chemoattractant protein-1 in the dorsal root ganglia in a rat model of neuropathic pain: possible involvement in the development of neuropathic pain. Neurosci Res 48:463-469.
Taskinen HS, Röyttä M (2000) Increased expression of chemokines (MCP-1, MIP-1alpha, RANTES) after peripheral nerve transection. J Peripher Nerv Syst 5:75-81.
Tillie-Leblond I, Hammad H, Desurmont S, Pugin J, Wallaert B, Tonnel AB, Gosset P (2000) CC chemokines and interleukin-5 in bronchial lavage fluid from patients with status asthmaticus. Potential implication in eosinophil recruitment. Am J Respir Crit Care Med 162:586-592.
Tokami H, Ago T, Sugimori H, Kuroda J, Awano H, Suzuki K, Kiyohara Y, Kamouchi M, Kitazono T; REBIOS Investigators (2013) RANTES has a potential to play a neuroprotective role in an autocrine/paracrine manner after ischemic stroke. Brain Res 1517:122-132.
van der Laan LJ, Ruuls SR, Weber KS, Lodder IJ, Dopp EA, Dijkstra CD (1996) Macrophage phagocytosis of myelin in vitro determined by flow cytometry: phagocytosis is mediated by CR3 and induces production of tumor necrosis factor-alpha and nitric oxide. J Neuroimmunol 70:145-152.
Walter J (2016) The triggering receptor expressed on myeloid cells 2: a molecular link of neuroinflammation and neurodegenerative diseases. J Biol Chem 291:4334-4341.
Walter W, Sánchez-Cabo F, Ricote M (2015) GOplot: an R package for visually combining expression data with functional analysis. Bioinformatics 31:2912-2914.
Webber C, Zochodne D (2010) The nerve regenerative microenvironment: early behavior and partnership of axons and Schwann cells. Exp Neurol 223:51-59.
Yi S, Tang X, Yu J, Liu J, Ding F, Gu X (2017) Microarray and qPCR analyses of wallerian degeneration in rat sciatic nerves. Front Cell Neurosci 11:22.
Yi S, Zhang H, Gong L, Wu J, Zha G, Zhou S, Gu X, Yu B (2015) Deep sequencing and bioinformatic analysis of lesioned sciatic nerves after crush injury. PLoS One 10:e0143491.
Yu B, Zhou S, Wang Y, Qian T, Ding G, Ding F, Gu X (2012) miR-221 and miR-222 promote Schwann cell proliferation and migration by targeting LASS2 after sciatic nerve injury. J Cell Sci 125:2675-2683.
Yu J, Gu X, Yi S (2016) Ingenuity pathway analysis of gene expression profiles in distal nerve stump following nerve injury: insights into wallerian degeneration. Front Cell Neurosci 10:274.
C-Editor: Zhao M; S-Editors: Yu J, Li CH; L-Editors: Yu J, Song LP; T-Editor: Jia Y
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6]