Critical signaling pathways during Wallerian degeneration of peripheral nerve

Qiong Cheng, Ya-xian Wang, Jun Yu, Sheng Yi

Key Laboratory of Neuroregeneration of Jiangsu and Ministry of Education, Co-innovation Center of Neuroregeneration, Nantong University, Nantong, Jiangsu Province, China

Critical signaling pathways during Wallerian degeneration of peripheral nerve

Qiong Cheng#, Ya-xian Wang#, Jun Yu, Sheng Yi＊

Key Laboratory of Neuroregeneration of Jiangsu and Ministry of Education, Co-innovation Center of Neuroregeneration, Nantong University, Nantong, Jiangsu Province, China

Graphical Abstract

Wallerian degeneration is a critical biological process that occurs in distal nerve stumps aer nerve injury. To systematically investigate molecular changes underlying Wallerian degeneration, we used a rat sciatic nerve transection model to examine microarray analysis outcomes and investigate significantly involved Kyoto Enrichment of Genes and Genomes (KEGG) pathways in injured distal nerve stumps at 0, 0.5, 1, 6, 12, and 24 hours, 4 days, 1, 2, 3, and 4 weeks aer peripheral nerve injury. Bioinformatic analysis showed that only one KEGG pathway (cytokine-cytokine receptor interaction) was significantly enriched at an early time point (1 hour post-sciatic nerve transection). At later time points, the number of enriched KEGG pathways initially increased and then decreased.ree KEGG pathways were studied in further detail: cytokine-cytokine receptor interaction, neuroactive ligand-receptor interaction, and axon guidance. Moreover, temporal expression patterns of representative differentially expressed genes in these KEGG pathways were validated by real time-polymerase chain reaction. Taken together, the above three signaling pathways are important aer sciatic nerve injury, and may increase our understanding of the molecular mechanisms underlying Wallerian degeneration.

nerve regeneration; Wallerian degeneration; sciatic nerve transection; peripheral nerve regeneration; microarray; bioinformatic analysis; Kyoto Enrichment of Genes and Genomes; signaling pathway; cytokine-cytokine receptor interaction; neuroactive ligand-receptor interaction; axon guidance; neural regeneration

Introduction

Wallerian degeneration is an important degenerative process that occurs in distal nerve stumps in response to nerve fiber injury (Coleman, 2005; Coleman and Freeman, 2010). It occurs in both the central nervous system and peripheral nervous system, although it normally occurs within 1 to 2 weeks aer injury in the peripheral nervous system but does not occur until a few months or even years after injury to the central nervous system (Grif fin et al., 1992; George and Grif fin, 1994). Timely occurrence of Wallerian degeneration in the peripheral nervous system may contribute to axonal regeneration due to clearance of myelin debris and growth inhibitors, and subsequent establishment of a regenerative microenvironment (Avellino et al., 1995; Vargas and Barres, 2007). However, delayed Wallerian degeneration in the central nervous system may hinder axonal regeneration. Accordingly, it is believed that Wallerian degeneration plays a key role in nerve regeneration (Lunn et al., 1989; Brown et al., 1991, 1992).

Considering the importance of Wallerian degeneration, numerous studies have been performed to identify underlying biological changes. These studies show that macrophages, monocytes, and Schwann cells work together to remove axon and myelin debris, and clear a path for subsequent axonal regrowth and nerve regeneration (Geuna et al., 2009; Sta et al., 2014).ese morphological and genetic studies have identified many central factors, including nicotinamide mononucleotide adenylyltransferase 2 (Coleman and Freeman, 2010; Gilley and Coleman, 2010; Gilley et al., 2013). Furthermore, emerging high-throughput studies have been performed to decipher global molecular changes. For example, in our previous study, we jointly re-annotated and re-analyzed microarray data (Yao et al., 2012, 2013) using bioinformatic tools including Euclidean distance calculation, hierarchical clustering, principle component analysis, gene ontology analysis, Kyoto Enrichment of Genes and Genomes (KEGG) analysis, and Ingenuity Pathway Analysis. Altogether, we obtained an integrated global view of genetic changes in injured distal nerve stumps (Yu et al., 2016; Yi et al., 2017). In particular, KEGG analysis outcomes identified pathways withP-values less than 0.05, indicating they may play critical roles in Wallerian degeneration (Yi et al., 2017). Taking the importance of signaling pathways into account, in the current study, we examined in detail these significantly enriched pathways in distal nerve stumps at various time points following sciatic nerve transection. Our aim was to achieve greater insight into dynamic molecular changes underlying Wallerian degeneration and identify critical biological processes for treatment of peripheral nerve repair and regeneration.

Materials and Methods

Rat sciatic nerve transection

A total of 66 adult, 2-month-old, male Sprague-Dawley rats weighing 180 to 220 g were obtained from the Experimental Animal Center, Nantong University, China (animal license No. SCXK [Su] 2014-0001 and SYXK [Su] 2012-0031). All animal surgery procedures were performed in accordance with Institutional Animal Care guidelines of Nantong University, and ethically approved by the Administration Committee of Experimental Animals, Jiangsu Province, China (No. 20160229-007).e experiment followed the National Institutes of Health Guide for the Care and Use of Laboratory Animals (NIH Publications No. 8023, revised 1978), and“Consensus Author Guidelines on Animal Ethics and Welfare” by the International Association for Veterinary Editors (IAVE).

Figure 1 Significant differentially expressed genes in axon guidance signaling.

Table 2 Significant differentially expressed genes in cytokine-cytokine receptor interaction

Table 1 List of significantly involved Kyoto Enrichment of Genes and Genomes at various time points during Wallerian degeneration identified by bioinformatic analysis

Rat sciatic nerve transection was performed, as previously described (Yu et al., 2016; Yi et al., 2017). Briefly, Sprague-Dawley rats were equally and randomly divided into eleven groups with 6 rats in each group. Following anesthetization, rat hair was shaved and the surgical area cleansed with 75% ethanol. An incision was made on the lateral aspect of the mid-thigh of the rat lehind limb, the sciatic nerve was lifted, and a 10-mm long segment was removed from the middle of the femur. At 0.5, 1, 6, 12, and 24 hours, 4 days, and 1, 2, 3, and 4 weeks post-sciatic nerve transection, rats were sacrificed by decapitation and distal nerve stumps were collected. Sham-operated rats (rats with sciatic nerves exposed but uninjured) were used as controls and designated as 0 hour post-sciatic nerve transection.

RNA extraction and microarray analysis

RNA samples were extracted from distal nerve stumps using Trizol reagent (Life Technology, Carlsbad, CA, USA). Remaining DNA was removed using RNeasy spin columns (Qiagen, Valencia, CA, USA). Purified RNA samples were then quantified using a NanoDrop ND-1000 spectrophotometer (Infinigen Biotechnology Inc., City of Industry, CA, USA).

Microarray analysis was performed in accordance with previous studies (Yao et al., 2012, 2013). Briefly, an Affymetrix GeneChip Oven 640 and Gene Array Scanner 3000 (Affymetrix, Santa Clara, CA, USA) were used to obtain microarray outcomes.ese outcomes were then analyzed by the R soware platform (v.2.13.0) and limma (linear regression model) package (Ritchie et al., 2015).

Bioinformatic analysis

Expression levels of mRNAs at 0.5, 1, 6, 12, and 24 hours, 4 days, and 1, 2, 3, and 4 weeks post-sciatic nerve transection were compared with those at 0 hour post-nerve transection. Genes with fold changes > 2 or < -2 (absolute value of log2 fold change > 1) and adjustedP-values < 0.05 were considered to be differentially expressed. Differentially expressed genes were then systematically analyzed using Database for Annotation, Visualization, and Integrated Discovery to enrich significant involved KEGG pathways.

Quantitative real time-polymerase chain reaction

RNA samples (0.5 μg) were reverse transcribed to cDNA using the Prime-Script Reagent Kit (TaKaRa, Dalian, China) for subsequent amplification. Quantitative real time-polymerase chain reaction (RT-PCR) was then performed using SYBR Green Premix Ex Taq (TaKaRa) with specific primer pairs on an Applied Biosystems Stepone real-time PCR System (Applied Biosystems, Foster City, CA, USA). Primer pair sequences were: Eph receptor A4 (Epha4), (forward) 5′-CGC CGT AGT ATC AGT GGG TG-3′ and (reverse) 5′-GTC TGT TCG GTA CTG GCT CA-3′; met proto-oncogene (Met), (forward) 5′-CGC TGC AGG CTG TGG ATT TA-3′and (reverse) 5′-GGT GAA ATG TGC TGT GCG AG-3′; interleukin 11 (Il11), (forward) 5′-CCG ACT GGA ACG GCT ACT TC-3′ and (reverse) 5′-GAC GAT GTC GAT GGT GGC TT-3′; prostaglandin E receptor 2 (Ptger2), (forward) 5′-TTC TAT GGC GGA GAC GG-3′ and (reverse) 5′-GGT CCC ACT TTT CCT TTC GGG-3′; and glyceraldehyde 3-phosphate dehydrogenase (Gapdh), (forward) 5′-CCT TCA TTG ACC TCA ACT ACA TG-3′ and (reverse) 5′-CTT CTC CAT GGT GAA GAC-3′. Ct values were obtained for mRNA quantification using the ΔΔCt method andGapdhas the reference gene.

Statistical analysis

Data were presented as the mean ± SEM. Differences between groups were tested by one-way analysis of variance using GraphPad Prism 6.0 (GraphPad Soware, Inc., San Diego, CA, USA).P-values < 0.05 were considered statistically significant.

Results

Overview of significantly involved KEGG pathways during Wallerian degeneration

Our previous study showed that many KEGG pathways are significantly involved in Wallerian degeneration (Yi et al., 2017). To further study these activated KEGG pathways, we analyzed enriched KEGG pathways at various time points post-sciatic nerve transection. KEGG pathways withP-values less than 0.05 are listed inTable 1.

At 0.5 hour post-sciatic nerve transection, no KEGG pathway was significantly involved. At 1 hour post-nerve transection, only one KEGG pathway, cytokine-cytokine receptor interaction was significantly activated. At 6 hours postnerve transection, hematopoietic cell lineage, neuroactive ligand-receptor interaction, Jak-STAT signaling pathway, and mitogen-activated protein kinases signaling pathway were also activated. At 12 hours post-nerve transection, more KEGG pathways were involved including those involved in cancer, calcium signaling pathway, axon guidance, insulin signaling pathway, and Prion diseases. At later time points, the number of significantly involved KEGG pathways decreased. Some of the above-mentioned pathways (e.g., neuroactive ligand-receptor interaction) remained highly activated. Some novel KEGG pathways emerged at longer time points post-nerve transection, such as cell adhesion molecules (CAMs) and tight junctions.

Cytokine-cytokine receptor interaction

Cytokine-cytokine receptor interaction (rno04060) was the only KEGG pathway activated immediately post-sciatic nerve transection, and strongly activated thereaer at early stages (1, 6, and 12 hours) following nerve transection. A schematic network of the cytokine-cytokine receptor interaction was built based on the KEGG Orthology Database (http://www.genome.jp/kegg/ko.html), with differentially expressed genes labeled (Table 2). Many chemokines, cytokines, and their corresponding receptors, including C-X-C motif chemokine ligand 2 (Cxcl2), interleukin 6 (Il6), interleukin-11 (Il11), colony stimulating factor 3 receptor (Csf3r), colony stimulating factor 2 (Csf2), MET proto-oncogene, receptor tyrosine kinase (Met), interferon A (Ifna), interleukin 1A (Il1a), and interleukin 1R2 (Il1r2) were significantly up-regulated in injured distal nerve stumps. Only one gene, kinase insert domain receptor (Kdr), a gene that encodes vascular endothelial growth factor receptor, was down-regulated post-sciatic nerve transection.

Neuroactive ligand-receptor interaction

Neuroactive ligand-receptor interaction (rno04080) was another KEGG pathway activated at a relatively early time point (6 hours post-sciatic nerve transection). Further, this KEGG pathway was activated at all time points from 6 hours to 4 weeks post-nerve transection, suggesting that it is critical for Wallerian degeneration and subsequent nerve regeneration. Receptors of numerous neurotransmitters and mediators, such as histamine receptor (Htr), prostaglandin E2 receptor 1 (Ptger), adenosine receptor (Adora), and leukotriene B4 receptor (Ltb4r) were up-regulated. In contrast, dopamine receptor (Drd), neurotensin receptor (Ntsr), MAS proto-oncogene or mas-related G-protein coupled receptor A (Mas1), cysteinyl-leukotriene receptor (Cysltr), corticotropin-releasing hormone receptor (Crhr), metabotropic glutamate receptor (Grm), purinergic receptor P2X (P2rx), and leptin receptor (Lepr) were down-regulated (Table 3).

Axon guidance

Besides cytokine-cytokine receptor interaction and neuroactive ligand-receptor interaction, two KEGG pathways significantly enriched at various time points during Wallerian degeneration, axon guidance signaling pathway (rno04360) was a KEGG pathway activated at 12 hours post-nerve transection, and also studied in detail based on its importance forsuccessful nerve regeneration and functional reconstruction. Accordingly,Met, a previously identified activated gene in the cytokine-cytokine receptor interaction pathway, was also involved in the axon guidance pathway. In addition, the ephrin receptor,Epha, was also up-regulated while the netrin receptor,Unc-5, was down-regulated (Figure 1).

Table 3 Significant differentially expressed genes in neuroactive ligand-receptor interaction

Figure 2 Quantitative real time-polymerase chain reaction for mRNA expression of Epha4 (A), Met (B), Il11 (C), and Ptger2 (D).

RT-PCR verification of microarray results

RT-PCR experiments were performed to validate temporal expression patterns of representative differentially expressed genes involved in the KEGG pathways: cytokine-cytokine receptor interaction (Il11), neuroactive ligand-receptor interaction (Ptger), and axon guidance (EphaandMet). Consistent with our microarray analysis outcomes, RT-PCR confirmed thatEpha4,Met,Il11, andPtger2expression levels were robustly increased, especially at longer time points post-sciatic nerve transection (Figure 2).

Discussion

Wallerian degeneration helps form a penetrable pathway for axon regrowth, and consequently, is very important for subsequent nerve repair and regeneration (Bittner et al., 2016). Until now, many critical factors for Wallerian degeneration have been identified. However, Wallerian degeneration is a complex biological process mediated by a group of molecules instead of a single gene or protein.us, obtaining a global view of molecular changes during Wallerian degeneration is of utmost importance. Accordingly here, by combined use of microarray and bioinformatic analysis, we synthetically analyzed highly activated KEGG signaling pathways (cytokine-cytokine receptor interaction, neuroactive ligand-receptor interaction, and axon guidance signaling pathway) following sciatic nerve transection.

Consistent with previous observations (Li et al., 2013, 2014; Yao et al., 2013; Yi et al., 2017), the KEGG pathway, cytokine-cytokine receptor interaction, was significantly involved from an early stage post-sciatic nerve injury. Detailed study of the cytokine-cytokine receptor interaction KEGG pathway suggests that most differentially expressed genes are related to immunoregulatory, inflammatory, and host defense processes such asCxcl2,Il6,Il11,Csf3r,Csf2,Ifna,Il1a, andIl1r2. Indeed, Wallerian degeneration has been thought of as an innate-immune response to external injuries (Rotshenker, 2011). Our outcomes demonstrate that immune reactions and inflammatory responses are initiated at an acute phase post-nerve injury, and suggest that innate immune and inflammatory responses in Wallerian degeneration might be critical for successful nerve regeneration and functional reconstitution.

However, it is worth noting that although we obtained some knowledge of dynamic molecular changes and essential signaling pathways during Wallerian degeneration, it remains unclear which cell types mediate these dynamic changes. Schwann cells and macrophages are important cells in distal nerve stumps, and consequently these dynamic changes may occur in these cell types. In our future studies, single cell sequencing will be performed to further decipher temporal expression patterns of genes in each cell type, and the specific role of each cell type during Wallerian degeneration.

Taken together, we systematically investigated significantly enriched KEGG pathways in distal nerve stumps following sciatic nerve transection, and identified critical KEGG pathways for peripheral nerve repair and regeneration. Our results may help elucidate critical genes, signaling pathways, and biological processes during Wallerian degeneration, and might contribute to identification of potential treatments for peripheral nerve repair and regeneration.

Author contributions:SY conceived and designed the experiments, and wrote the paper. QC, YXW, JY, and SY performed the experiments and analyzed the data. QC and SY provided reagents/materials/analysis tools. All authors approved the final version of the paper.

Conflicts of interest:None declared.

Research ethics:The study protocol was approved by the Administration Committee of Experimental Animals, Jiangsu Province, China (No. 20160229-007).e experiment followed the National Institutes of Health Guide for the Care and Use of Laboratory Animals (NIH Publications No. 8023, revised 1978), and “Consensus Author Guidelines on Animal Ethics and Welfare” produced by the International Association for Veterinary Editors (IAVE).e article was prepared in accordance with the “Animal Research: Reporting of In Vivo Experiments Guidelines” (ARRIVE Guidelines).

Open access statement:

Contributor agreement:A statement of “Publishing Agreement” has been signed by an authorized author on behalf of all authors prior to publication.

Plagiarism check:This paper has been checked twice with duplication-checking soware ienticate.

Peer review:A double-blind and stringent peer review process has been performed to ensure the integrity, quality and significance of this paper.Open peer reviewer:Bogdan K Beirowski, State University of New York, USA.

Avellino AM, Hart D, Dailey AT, MacKinnon M, Ellegala D, Kliot M (1995) Differential macrophage responses in the peripheral and central nervous system during wallerian degeneration of axons. Exp Neurol 136:183-198.

Brown MC, Lunn ER, Perry VH (1991) Poor growth of Mammalian motor and sensory axons into intact proximal nerve stumps. Eur J Neurosci 3:1366-1369.

Brown MC, Lunn ER, Perry VH (1992) Consequences of slow Wallerian degeneration for regenerating motor and sensory axons. J Neurobiol 23:521-536.

Coleman M (2005) Axon degeneration mechanisms: commonality amid diversity. Nat Rev Neurosci 6:889-898.

Coleman MP, Freeman MR (2010) Wallerian degeneration, wld(s), and nmnat. Annu Rev Neurosci 33:245-267.

George R, Grif fin JW (1994) Delayed macrophage responses and myelin clearance during Wallerian degeneration in the central nervous system: the dorsal radiculotomy model. Exp Neurol 129:225-236.

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.

Gilley J, Coleman MP (2010) Endogenous Nmnat2 is an essential survival factor for maintenance of healthy axons. PLoS Biol 8:e1000300.

Gilley J, Adalbert R, Yu G, Coleman MP (2013) Rescue of peripheral and CNS axon defects in mice lacking NMNAT2. J Neurosci 33:13410-13424.

Griffin JW, George R, Lobato C, Tyor WR, Yan LC, Glass JD (1992) Macrophage responses and myelin clearance during Wallerian degeneration: relevance to immune-mediated demyelination. J Neuroimmunol 40:153-165.

Li M, Guo W, Zhang P, Li H, Gu X, Yao D (2013) Signal flow and pathways in response to early Wallerian degeneration after rat sciatic nerve injury. Neurosci Lett 536:56-63.

Li M, Zhang P, Guo W, Li H, Gu X, Yao D (2014) Protein expression profiling during wallerian degeneration aer rat sciatic nerve injury. Muscle Nerve 50:73-78.

Lunn ER, Perry VH, Brown MC, Rosen H, Gordon S (1989) Absence of wallerian degeneration does not hinder regeneration in peripheral nerve. Eur J Neurosci 1:27-33.

Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK (2015) limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res 43:e47.

Rotshenker S (2011) Wallerian degeneration: the innate-immune response to traumatic nerve injury. J Neuroinflammation 8:109.

Shin YH, Chung HJ, Park C, Jung J, Jeong NY (2014) Adenosine 5’-triphosphate (ATP) inhibits schwann cell demyelination during Wallerian degeneration. Cell Mol Neurobiol 34:361-368.

Sta M, Cappaert NL, Ramekers D, Baas F, Wadman WJ (2014) The functional and morphological characteristics of sciatic nerve degeneration and regeneration aer crush injury in rats. J Neurosci Methods 222:189-198.

Vargas ME, Barres BA (2007) Why is Wallerian degeneration in the CNS so slow? Annu Rev Neurosci 30:153-179.

Yao D, Li M, Shen D, Ding F, Lu S, Zhao Q, Gu X (2012) Gene expression profiling of the rat sciatic nerve in early Wallerian degeneration aer injury. Neural Regen Res 7:1285-1292.

Yao D, Li M, Shen D, Ding F, Lu S, Zhao Q, Gu X (2013) Expression changes and bioinformatic analysis of Wallerian degeneration aer sciatic nerve injury in rat. Neurosci Bull 29:321-332.

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.

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.

Copyedited by James R, de Souza M, Yu J, Li CH, Qiu Y, Song LP, Zhao M

How to cite this article: Cheng Q, Wang YX, Yu J, Yi S (2017) Critical signaling pathways during Wallerian degeneration of peripheral nerve. Neural Regen Res 12(6):995-1002.

＊Correspondence to: Sheng Yi, Ph.D., syi@ntu.edu.cn.

orcid:

0000-0003-1316-3370

(Sheng Yi)

10.4103/1673-5374.208596

Accepted: 2017-05-17