- Open Access
An RNA-sequencing transcriptome of the rodent Schwann cell response to peripheral nerve injury
Journal of Neuroinflammation volume 19, Article number: 105 (2022)
The important contribution of glia to mechanisms of injury and repair of the nervous system is increasingly recognized. In stark contrast to the central nervous system (CNS), the peripheral nervous system (PNS) has a remarkable capacity for regeneration after injury. Schwann cells are recognized as key contributors to PNS regeneration, but the molecular underpinnings of the Schwann cell response to injury and how they interact with the inflammatory response remain incompletely understood.
We completed bulk RNA-sequencing of Schwann cells purified acutely using immunopanning from the naïve and injured rodent sciatic nerve at 3, 5, and 7 days post-injury. We used qRT-PCR and in situ hybridization to assess cell purity and probe dataset integrity. Finally, we used bioinformatic analysis to probe Schwann cell-specific injury-induced modulation of cellular pathways.
Our data confirm Schwann cell purity and validate RNAseq dataset integrity. Bioinformatic analysis identifies discrete modules of genes that follow distinct patterns of regulation in the 1st days after injury and their corresponding molecular pathways. These findings enable improved differentiation of myeloid and glial components of neuroinflammation after peripheral nerve injury and highlight novel molecular aspects of the Schwann cell injury response such as acute downregulation of the AGE/RAGE pathway and of secreted molecules Sparcl1 and Sema5a.
We provide a helpful resource for further deciphering the Schwann cell injury response and a depth of transcriptional data that can complement the findings of recent single cell sequencing approaches. As more data become available on the response of CNS glia to injury, we anticipate that this dataset will provide a valuable platform for understanding key differences in the PNS and CNS glial responses to injury and for designing approaches to ameliorate CNS regeneration.
Like the CNS, the peripheral nerve is protected within several layers of connective tissue and behind a diffusion barrier, the blood–nerve barrier. Within the epineurium, the peripheral nerve is composed of endoneurial fibroblasts, endothelial cells, pericytes, resident macrophages, Schwann cells, and efferent and afferent axons . Following breakdown of the blood–nerve barrier the canonical process of Wallerian degeneration ensues, recruiting additional cell types into the nerve. Neutrophils occupy the nerve 1–2 days after injury followed by invasion by hematogenous macrophages, which persist until approximately 2 weeks after injury  (Fig. 1A).
Understanding the mechanisms underlying the remarkable ability of the PNS to regenerate after injury has been a topic of intense interest for decades. The Schwann cell, a neural crest-derived glial cell type comprising 80% of the cells in the uninjured nerve has long been recognized as a key player in PNS regeneration, undergoing a remarkable cellular transformation after nerve injury and assuming multiple pro-regenerative functions during recovery of the injured nerve . While the Schwann cell transformation after nerve injury was long thought to represent de-differentiation to a more immature state, we now know that post-injury Schwann cells are distinct from their precursors. Genetically and morphologically, post-injury Schwann cells, or “Repair cells”, are unique, transdifferentiated cells that facilitate tissue rebuilding. Repair cells turn off genes for myelination and activate those for diverse features including upregulation of growth-promoting proteins (GDNF, artemin, BDNF and others), innate immune activation (secretion of cytokines TNFα, LIF, IL1α, IL1β, and MCP-1), structural reorganization, and myelin debris clearance [18, 20, 21]. Numerous recent studies have made considerable progress towards precisely defining the molecular changes behind the Schwann cell injury response, implicating key transcriptional regulators and signaling pathways in this transformation [1, 5, 18, 33].
Until recently, knowledge of global changes in Schwann cell gene expression during the 1st week after injury was largely based on data obtained from whole sciatic nerve before and after crush or transection injury [27, 32, 44, 49]. This data, however, does not allow for distinction between Schwann cell-specific transcriptomic changes and gene expression in other cell types within the nerve, particularly in hematogenous macrophages, which flood the nerve in large numbers beginning 3 days after injury . Purification of Schwann cells away from the other cell types of the peripheral nerve is essential to gain a clear view of the Schwann cell response to injury. Another valuable resource in the field has been microarray data from neonatal Schwann cells in which Erk signaling is activated to mimic injury . While also instructive, this dataset depicts Schwann cells purified at an early developmental timepoint, cultured in vitro, and activated in an in vitro environment in the absence of axon contact, and is thus unlikely to recapitulate all aspects of in vivo transdifferentiation. Circumventing the need for cell purification, several groups have recently performed single cell transcriptomic analysis of the rodent sciatic nerve following autoimmune and mechanical injury [22, 40, 46]. These studies provide exciting insight into heterogeneity of post-injury gene expression changes in Schwann cells and other cell peripheral nerve cell populations, but lack the sequencing depth to identify more subtle responses to injury. Thus, we view this work as a fundamental new resource that will complement the excellent studies above.
Here, we provide a high-resolution RNAseq transcriptome of Schwann cells acutely purified from the uninjured and injured rat sciatic nerve at 3, 5, and 7 days after nerve crush injury at postnatal day (P)18 . In addition, we provide RNAseq analysis of the whole sciatic nerve at the same timepoints. To our knowledge, this is the first deep sequencing RNAseq dataset of purified Schwann cells at across this timespan, a period known to involve dynamic changes in Schwann cell structure and function. Bioinformatics analysis of our dataset highlights dataset integrity and the importance of purifying Schwann cells away from other cell types in the injured nerve to identify Schwann cell-specific gene changes. We identify the most highly up- and down-regulated transcripts at each post-injury timepoint and further classify genes into discrete modules that follow distinct patterns of regulation with respect to time after PNS injury. Interrogation of the genes in these modules identifies aspects of the Schwann cell response to injury with potential for future study. We also compare immune roles of Schwann cells with those of invading macrophages to better differentiate injury responses in these two key cell types. To ensure widespread distribution of our dataset, we have incorporated these data in a user-friendly website previously created by Ben Barres’ lab (www.BrainRNAseq.org) that provides a simple and searchable platform for browsing gene expression data. Our dataset provides a helpful resource for further deciphering the Schwann cell injury response. Furthermore, as additional data become available on the response of CNS glia to injury, we anticipate that this dataset will provide a valuable platform for understanding key differences in the PNS and CNS glial responses to injury and for designing approaches to ameliorate CNS regeneration.
All procedures involving animals were conducted in conformity with Stanford University guidelines and are in compliance with national and state laws and policies according to APLAC protocol 10726.
Sciatic nerve crush injury
All surgical experiments were performed under 2.5% isoflurane on rats at postnatal day 18. Sciatic nerve crush injury was performed as previously described . Briefly, rats were anesthetized using isoflurane. Upper thigh was shaved and sterilized using isopropanol. A 1-cm incision was made using a scalpel and nerve was visualized via blunt dissection using forceps. The left sciatic nerve was crushed at mid-thigh for 5 s using forceps marked with sterile graphite to mark the crush site. Carprofen (5 mg/kg subcutaneous) was administered for analgesia. No sham operation was performed prior to the collection of naïve nerves.
Purification of Schwann cells from the naïve and injured sciatic nerve
10–20 sciatic nerves were used for the generation of each sample of purified cells. For injured nerves, sciatic nerve tissue was collected beginning proximally 1 mm distal to the graphite-marked crush site and extending distally along the tibial nerve branch as far as possible approaching the medial calcaneus. The proximal portion of the peroneal nerve branch was included when intact dissection was possible. Schwann cells were isolated from the entirety of the harvested nerve segment. For uninjured nerves, the anatomically equivalent nerve segment was isolated. Schwann cells were acutely purified using immunopanning from the injured or uninjured rat sciatic nerve according to Lutz . Briefly, following enzymatic digestion using collagenase/dispase, nerves were mechanically dissociated to create a single cell suspension. Myelin debris was removed using a nitex filter and BSA cushion. Sequential immunopanning first depleted unwanted cell populations (CD45 and Thy1) and then selected for Schwann cells (O4). Performing depletion of other cell populations prior to positive selection of Schwann cells ensures optimal Schwann cell purity. Yields per nerve ranged from 10 to 20,000 cells for uncrushed nerves to approximately 100,000 cells per nerve after nerve crush injury.
Immunocytochemistry of immunopanned cells
Cultured cells were fixed with 4% paraformaldehyde for 10 min at room temperature 24 h after immunopanning and then permeabilized and blocked with 10% goat serum, 0.2% TritonX-100. Cells were labeled using antibodies to goat Iba1 (Abcam) and mouse S100 (Sigma s2532) as well as Cellmask (ThermoFisher H34558) and DAPI.
RNA isolation from whole nerves
RNA from whole nerves was purified according to the RNeasy Fibrous Tissue kit. For injured nerves, only the nerve segment distal to the site of crush injury was included. For uninjured nerves the anatomically equivalent nerve segment was isolated. A minimum of 10 nerves was used for the generation of each RNA sample. RNA quality was verified through Bionalyzer analysis. All samples submitted for RNASeq had a RIN of 8 or above.
RNA isolation from purified Schwann cells
RNA from purified Schwann cells was isolated at the end of the immunopanning procedure by scraping the O4-positive panning dish according to the instructions of the RNeasy MICRO kit. RNA quality was verified through Bionalyzer analysis. All samples submitted for RNASeq had a RIN of 8 or above.
RNA-seq library construction and sequencing
Total RNA that passed QC standards (RIN > 8.0 as assessed by bioanalyzer) proceeded to library synthesis using the Ovation® RNA-seq system V2 (Nugen 7102). First and second-strand cDNA synthesis and SPIA amplification were performed following the manufacturer’s instructions, and cDNA was fragmented with a sonicator (Covaris S2) using the following parameters: duty cycle 10%, intensity 5, cycles/burst 100, time 5 min. We then used the Next Ultra RNA-seq library prep kit for Illumina (NEB E7530) and NEBNext® multiplex oligos for Illumina® (NEB E7335 E7500) to perform end repair, adaptor ligation, and 5–6 cycles of PCR enrichment according to manufacturer’s instructions. The quality of the libraries were then assessed by bioanalyzer and qPCR and high-quality libraries were sequenced by the Illumina NextSeq sequencer to obtain 150-bp pair-end reads to a minimum depth of 20 million reads per sample. All libraries passed standard QC metrics (average Q-scores > 35, mapping % > 80).
Bioinformatic processing of RNAseq data
FASTQ files were first groomed using FASTQ groomer and then mapped using TopHat2, which invokes Bowtie as an internal read mapper. The paired end option was selected and rat genome Rnor_6.0 was used as the reference genome. We then ran Cufflinks to assemble transcripts and estimate expression level as fragments per kilobase of transcript sequence per million mapped fragments (FPKM). All subsequent data analysis was performed using R software (https://www.r-project.org/). We deposited all raw .fastq files from the RNA-seq data in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO, Accession number GSE177037). Unsupervised hierarchical clustering of FPKM data was performed in R using the top 2000 most variable genes with the following settings: Linkage; average, ColumnPdistance; Spearman, standardized across rows.
Differential expression and pathway analysis
To assess the functional difference of genes post-nerve crush, we carried out pathway analysis on the top 2000 differentially expressed genes at each timepoint. FASTQ files were aligned using STAR  (v2.7.1) to the Rnor 6.0 reference genome and counts were generated using RSEM  (v1.3.1). Differentially expressed genes were determined using the DEseq2  (v1.26) package within the R environment. Differences were visualized as volcano plots using the EnhancedVolcano package  (v1.80) from Bioconductor. Pathway analysis was carried out using Bioconductor package pathfindR  with p < 0.01 for both the Gene Ontology-Molecular Function and KEGG databases.
Weighted gene co-expression network analysis (WGCNA)
WGCNA analysis was performed using the R package ‘‘WGCNA’’ [24, 25]. The power value within the blockwise consensus module was determined by examining the soft-threshold-mean-connectivity curve (p = 21). The minimum WGCNA module size was 30 genes. Modules with < 0.25 similarity were merged.
We designed primers spanning exon–exon junctions to avoid amplification of genomic DNA. The primers used are listed in Additional file 3. We performed 40 cycles of amplification for all samples. The specificity and efficiency of all primers were first validated using gel electrophoresis and qPCR. The determination of each gene’s CT in qPCR was performed in triplicate. When determining fold changes in gene expression across post-injury timepoints, the CT of each gene was normalized according to the CT of the housekeeping gene in the same sample.
In situ hybridization validation
Full-length rat cDNA expression plasmids for the genes analyzed were linearized and transcribed with T7, Sp6, or T3 RNA polymerases to generate fluorescein (FITC)-labeled single-stranded anti-sense riboprobes according to manufacturer’s instructions (11175025910; Roche) and purified over G-50 columns (GE Healthcare) (Additional file 4). Fresh frozen 12-µm-thick nerve cryosections from P18 rats were processed as described  for in situ hybridization and immunohistochemistry using a mouse S100 antibody (Sigma s2532).
Schwann cell RNA purification and RNAseq transcriptional profiling
To gain insight into the Schwann cell response to peripheral nerve injury, we generated an RNAseq resource for Schwann cells acutely purified from the naïve and injured rat sciatic nerve at 3, 5, and 7 days post-crush. For comparison, we also generated RNAseq data for the naïve rat sciatic nerve (Fig. 1B, D) (Additional file 1). We chose to examine gene changes over the first 7 days after mouse nerve crush injury because this time period encompasses the key events of Wallerian degeneration including axon degeneration, myelin debris formation, Schwann cell transdifferentiation, Schwann cell proliferation, the major phase of Schwann cell-mediated debris clearance, macrophage recruitment, and the onset of macrophage-mediated debris clearance . We aimed to capture the transcriptional changes associated with these events. The full dataset of purified Schwann cells and whole nerve segments can be found in Additional file 4.
To assess the purity of our O4-immunopanned Schwann cells we used immunocytochemistry and expression of cell type-specific genes. Immunocytochemistry of cells purified from uncrushed nerves using the Schwann cell-specific marker S100 demonstrated a > 95% pure population of Schwann cells with typical spindle-like morphology 24 h after immunopanning (Fig. 1C). Both O4 and S100 are broadly expressed in the Schwann cell lineage .
Examination of FPKM values for known cell type-specific genes in our whole nerve and purified Schwann cell RNAseq data revealed high enrichment for Schwann cell-specific transcripts and depletion of transcripts from potential contaminating cell types. Specifically, we demonstrated the depletion of S100 expressing perineurial glia and satellite glial transcripts in our purified Schwann cell samples (Additional file 9) [13, 39] In addition, we were able to deplete macrophage, endothelial cell, and fibroblast cell type-specific markers at all timepoints examined (Fig. 2A). Given the large influx of myeloid cells into the nerve after injury, we wanted to more closely examine the possibility that our Schwann cell samples post-injury might be contaminated by infiltrating hematogenous myeloid cells. Based on absolute FPKM, we confirmed that while myeloid markers indeed increase significantly in the whole nerve transcriptome post-injury, these were near completely depleted in our purified Schwann cell samples (Fig. 2B).
Supporting the notion that purification of Schwann cells away from the sciatic nerve is crucial for the analysis of the Schwann cell injury response in order to avoid confounding transcriptome changes in other cell types after nerve injury, unsupervised hierarchical clustering demonstrated that expression profiles of naïve whole nerve and naïve Schwann cells cluster together while profiles of injured whole nerve and injured Schwann cells form distinct clusters (Fig. 2C). As a further quality control measure, we examined genes known to be up- or down-regulated in Schwann cells in the injured nerve. As expected, our RNAseq profiles revealed persistent expression of lineage marker Sox10, a marked downregulation of myelin genes MBP and MPZ and upregulation of immature and non-myelinating Schwann cell marker, NGFR after nerve crush injury (Fig. 2D).
Validation of RNAseq results
Next, we examined the validity of our Schwann cell RNAseq transcriptome datasets by two orthologous methods. We began by identifying a list of transcripts that appeared from our RNAseq data to be highly regulated after nerve injury. To be selected, transcripts had to attain an FPKM of at least 30 FKPM at any timepoint, and were either highly upregulated (sevenfold or more increase in FPKM) or down-regulated (sevenfold or more decrease in FPKM) during the course of injury. First, we used qRT-PCR to validate a subset of these top differentially expressed Schwann cell genes in the week following nerve crush injury. We found consistent agreement between the qRT-PCR data and our RNAseq findings (Fig. 3). Following qRT-PCR validation, we used in situ hybridization to confirm Schwann cell origin for the transcriptional changes observed. We designed fluorescence in situ hybridization probes directed against a subset of the top up- and down-regulated transcripts (Additional file 2), selected based on (a) biological interest and (b) amenability to specific probe design. We then performed simultaneous immunohistochemistry with an antibody directed against the Schwann cell marker, S100 (Fig. 4A–D) to validate Schwann cell-specificity for each in situ probe.
Our in situ hybridization confirmed the anticipated changes in gene expression in S100-positive Schwann cells, with putative upregulated transcripts exhibiting no signal co-localization with S100 immunostaining in uninjured nerves and co-localization with S100-positive Schwann cells post-injury. Similarly, putative down-regulated transcripts (SPARCL1) exhibited signal co-localized with S100 immunostaining prior to injury and were no longer detectable in immunolabeled Schwann cells after injury (Fig. 4).
For a subset of up- and down-regulated transcripts, we observed parallel trends in gene expression in S100-negative cells, in particular macrophages. Therefore, we wanted to further confirm that the transcriptomic changes we denoted originated in Schwann cells, and also differentiate immune responses due to Schwann cell expression of injury-related genes from those mediated by CD45-expressing myeloid cells. We sequenced CD45-positive myeloid cells that were depleted from our sciatic nerve cell suspension during CD45-based immunopanning post-injury. This data revealed that although CD45-positive macrophages exhibit occasional parallel trends in gene expression compared with purified Schwann cells, the level of the expression changes for these genes in this infiltrating population (absolute FPKM) was typically at least an order of magnitude less than in purified Schwann cells (Fig. 4A–D). This finding indicates that Schwann cells relative to macrophages are the primary drivers of gene expression changes for these differentially expressed genes post-injury.
Analysis of Schwann cell-specific gene expression before and after peripheral nerve injury
We next used volcano plots to visualize global changes in Schwann cell gene expression at each post-crush timepoint relative to the naïve Schwann cell transcriptome (Fig. 5A–C) (calculated using the DESeq2 package). These comparisons revealed widespread transcriptomic changes (average = 5398, padj < 0.05) in Schwann cells at each time post-crush stage. The number of significant differentially expressed genes (DEGs) was relatively consistent at each timepoint (5495, 5744, and 4955 genes at days 3, 5, 7 after injury), although the composition of these genes differed throughout each stage. The full list of DEGs at each timepoint can be found in Additional file 5.
To identify the most relevant biological pathways and processes induced by nerve crush in Schwann cells, we performed Gene Ontogeny (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway analysis of the top 2000 most upregulated and downregulated genes at each post-crush timepoint (Fig. 5D–F, all genes with padj < 0.05) (Additional files 6 and 7). Among upregulated genes, to determine which biological categories exhibit stage-dependent changes, we quantified how many GO hits were present at each post-injury timepoints (Fig. 5G). While some categories were quite consistent throughout all timepoints (motility, phagocytosis), others displayed significant enrichment at either day 3 (cell division) or day 5 (cytoskeletal rearrangement). The most highly downregulated GO terms at days 3, 5, and 7 after injury consistently relate to extracellular matrix production and cell adhesion (Additional file 7). Our GO analysis of downregulated transcripts also highlighted a strong silencing of ion channel expression, likely reflecting the dismantling of highly specialized Schwann cell membrane regions in close proximity to axons [3, 7, 45] (Additional file 7).
Using WGCNA to identify co-regulated networks in the Schwann cell injury response
Finally, we wondered whether we could identify families of genes who displayed similar expression patterns across all post-injury timepoints within Schwann cells. While this analysis could include consistent up- and down-regulated trends throughout the duration of the injury, it could alternatively involve, for example, groups of genes that peak at 3 days post-injury and then decay, or those that remain low until before increasing at 7 days post-injury. Categorizing genes by these patterns allows for a more nuanced understanding of which pathways share temporal characteristics during injury. Therefore, we performed unsupervised grouping of Schwann cell-specific genes into ten distinct modules using weighted gene co-expression network analysis (WCGNA), each defined by a distinct pattern of gene expression over the examined timecourse (Fig. 6A, B). These clusters can be mined individually to discover potentially novel pathways and groups of Schwann cell genes regulated in distinct patterns after nerve injury. A full list of genes belonging to each module can be found in Additional file 8. To demonstrate the feasibility of clustering injury response genes in these groups, we sampled three modules (1, 5, and 8) that exhibit distinct temporal patterns of expression and asked which pathways were most abundant in each (Fig. 6C–E). Module 1 genes, which decline rapidly after injury, were enriched in genes associated with membrane organization and formation of membrane sub-domains, consistent with dismantling of myelin nodes and paranodes. This is in contrast to module 5 and 8 genes, which spike rapidly after injury and are most associated with metabolic and cell cycle processes, likely reflecting Schwann cell proliferation after injury.
Acute regulation of PI3K/Akt/mTOR and AGE/RAGE activation in the Schwann cell injury response
Our WCGNA analysis also highlighted the temporal regulation of specific molecular pathways that are of interest for further study. The PI3K/Akt pathway, implicated in multiple studies involving the Schwann cell response to injury, Schwann cell proliferation, and in the initiation of Schwann cell myelination, was grouped into module 8, a cluster of genes for which expression increases sharply by day 3 after injury and then decreases in expression intensity by day 5 and further by day 7 after injury [10, 35, 47]. WGCNA analysis also identified the AGE–RAGE pathway within module 7, a cluster of genes sharply downregulated at 3 days post-injury that gradually increase back towards baseline expression at days 5 and 7 after injury.
The transformation of the Schwann cell in response to peripheral nerve injury has long been recognized as a key feature of the peripheral nerve’s remarkable ability to regenerate and yet the molecular underpinnings of this transformation have remained incompletely understood. Here, we provide an RNAseq dataset on Schwann cells acutely purified from rodent peripheral nerves at multiple timepoints during the 1st week after nerve crush to provide insight into the transcriptional changes underlying the Schwann cell response to injury.
We used published immunopanning methods to purify Schwann cells away from other cell types in the nerve, including resident macrophages, hematogenous immune cells, and endoneurial fibroblasts at postnatal day 18. Several methodological points deserve brief discussion. First, O4 is highly expressed by all Schwann cells (both the myelinating and non-myelinating (Remak)) at rat postnatal day 18 . While it is plausible that some Schwann cells may not bind to our O4 coated panning for technical reasons such as incomplete tissue digestion, we have no reason to believe that we have excluded a biologically meaningful Schwann cell sub-population. Secondly, purification of viable Schwann cells from uninjured whole sciatic nerves is increasingly difficult as the nerve ages. Connective tissue becomes increasingly difficult to digest and myelinating cell processes are increasingly difficult to dissociate from axons without causing Schwann cell death. We chose postnatal day 18 as our initial timepoint because it represents a temporal stage at which both myelinating and unmyelinating Schwann cells are reaching maturity and yet can still be dissociated from the nerve in sufficient numbers [19, 29]. Nonetheless, potential subsequent maturational changes should be taken into account when interpreting our data as representative of the adult nerve. Finally, given that immunopanning is performed at room temperature, we do not anticipate significant changes in cellular gene profiles during the cell purification process, with the exception of immediate early genes. We have confirmed this by comparing transcriptomes of whole nerve with dissociated nerve (without immunopanning).
While recent single cell expression analyses contribute unprecedented granularity to our knowledge of Schwann cell gene expression, they lack the sequencing depth provided by bulk RNAseq to more accurately define and quantify the molecular changes occurring in Schwann cells after injury (~ 1000 genes per Schwann cell detected in Wolbert et al. vs ~ 8000 in this study). Thus, both approaches are valuable tools for understanding the Schwann cell response to injury and should be viewed as complementary to one another. Our initial analysis of our bulk RNAseq results underlines several familiar themes of Schwann cell transdifferentiation, as well as potential avenues for future study. These are outlined below.
There has long been consensus that Schwann cells distal to a site of peripheral nerve injury reorganize to form cellular columns known as Bungner’s bands. The importance of these cellular columns for successful nerve regeneration is also well-established, yet the mechanisms underlying the formation of these bands remains incompletely understood. Elegant studies using lineage tracing analysis have recently provided insight into the Schwann cell sub-type contributions and changes in Schwann cell structure that enable the formation of these columns . The molecular underpinnings of these changes, however, remain unclear. Our GO and KEGG analysis of highly upregulated transcripts revealed a large proportion of terms related to cytoskeletal rearrangement, tissue morphogenesis, and ECM receptor expression. Further mining of these transcriptional changes within our database could yield important insight into the regulation of molecular pathways that direct this aspect of the Schwann cell injury response.
Recent work has highlighted that in parallel with macrophages, Schwann cells contribute to debris clearance within the injured peripheral nerve [5, 15]. Accordingly, our KEGG analysis of genes highly upregulated after injury highlights endocytosis and lysosomal degradation as key upregulated cellular processes in the 1st days after injury. While macrophages express a broad repertoire of phagocytic receptors in the injured nerve, Schwann cell-mediated phagocytosis involves specific upregulation of TAM (Tyro, Axl, Mer) receptors Axl and Mertk . The involvement of lectin-mediated phagocytosis, specifically galectin-3, in Schwann cells debris clearance has also been suggested . The role of this pathway in peripheral nerve repair has however been debated, with galectin-3 knock-out mice exhibiting an unexpected enhanced capacity for peripheral nerve regeneration . Galactose metabolism was one of the most upregulated pathways revealed through our KEGG analysis, suggesting that this the involvement of this pathway in the Schwann cell response to injury may merit further study.
Extracellular matrix changes
Our transcriptional profiling of the post-injury sciatic nerve also highlights extracellular matrix changes. On the side of infiltrating myeloid cells, we observed a strong upregulation of glycan breakdown enzymes after injury. Two Schwann cell transcripts, Sparcl1 and Sema5a, are among the most significantly downregulated genes within the first week after injury. Sparcl1, or Hevin, a secreted non-neuronal protein is implicated in synapse initiation and remodeling during development and after injury . Most work to date on Sparcl1 has been in the CNS [6, 11, 23]. The magnitude of transcriptional regulation revealed in our RNAseq dataset suggests that this protein is likely playing a role in the Schwann cell injury response beyond synapse regulation at the neuromuscular junction. Sema5a, a thrombospondin repeat-containing axon guidance molecule is known to contribute to CNS axon outgrowth inhibition after injury . The effect of Schwann cell downregulation of Sema5a on PNS axon regeneration has not yet been studied.
Regulation of PI3K/Akt/mTOR and AGE/RAGE activation
A complex interplay between the Mek/ERK1/2-MAPK and PI3K/Akt/mTOR pathways is suggested in Schwann cell myelination and myelin pathologies . A dual role for the PI3K/Akt/mTOR pathway in Schwann cell development and myelination has been reported, with activation of the pathway required for the inhibition of premature Schwann cell differentiation during development, and then subsequently for positive promotion of myelin growth . Whether this pathway is differentially implicated in developmental myelination versus post-injury re-myelination has not been extensively studied. Our data suggest that activation of PI3K/Akt/mTOR following nerve injury is dramatic at the acute stage, but rapidly declines towards baseline levels by 1 week post-injury. Further study of the temporal regulation of this pathway could suggest an anticipated second wave of expression. Additionally, differential spatial regulation of the PI3K/Akt/mTor signaling between proximal and more distal portions of the regenerating nerve may elucidate differences between developmental and post-injury Schwann cell differentiation and myelination. In our CD45-positive myeloid cell sequencing results, we observed a significant upregulation of MAPK signaling, as reported following CNS injury .
The AGE/RAGE pathway has to date been most heavily studied in age-related vascular dysfunction where it has been found to mediate age-related vascular stiffening. The accumulation of AGEs is known to be accelerated in diabetes [9, 38]. Although the AGE receptor RAGE is expressed in Schwann cells, its role in the peripheral nerve has not yet been studied . Recent work shows that age-related changes in Schwann cells contribute to failure of the aged PNS to regenerate to the same extent as the young PNS . Our observation that a decrease in RAGE occurs in parallel with the Schwann cell reprogramming indicates that a silencing of this pathway is at least involved in, if not required for a successful Schwann cell response to injury. Such silencing may not be feasible in the aged and diabetic nerve, where AGEs have heavily accumulated. Failed downregulation of AGE–RAGE signaling could help to explain a defective Schwann cell response to injury in aged nerves and possibly a defective homeostatic nerve repair mechanism in diabetic neuropathy.
Relevance to human peripheral nerve injury
The ultimate value of findings in rodent models of nerve injury depends on their relevance to human disease. Data on the human Schwann cell response to injury are extremely scarce. Recent studies using degenerating nerve fascicles in vitro, however, demonstrate important parallels between rodent and human Schwann cells [31, 43]. Importantly, human Schwann cells appear to recapitulate key features of murine repair Schwann cells. Namely, human Schwann cells demonstrate adaptive cellular reprogramming after injury including cJUN upregulation and expression of immature Schwann cell markers L1CAM, NEUM and NGFR, downregulation of myelin-associated genes such as MBP and MYP0, and participation in debris clearance . Schwann cells in both species demonstrate a reduction in the Schwann cell injury response with increased age [5, 31]. Interestingly, work on the relative timing of early events in the human Schwann cell injury response indicates that the Schwann cell injury response in humans is delayed in comparison with mice . Specifically, human Schwann cells exhibit delayed myelin degradation, delayed transdifferentiation (slower/less robust silencing of differentiation genes and slower/less robust induction of repair genes) and reduced adaptation of lipid metabolism. The regenerative capacity of the rodent peripheral nerve (4.5 mm/day) is also known to exceed that of the human nerve (1–1.5 mm/day) and could be tightly linked to differences in the Schwann cell injury response . The underpinnings of these differences are not yet well understood. This recent work supports the relevance of work performed in murine models of peripheral nerve injury to human peripheral nerve repair. Additional studies in human Schwann cells are needed to deepen our understanding of the similarities and differences between Schwann cell injury responses in these two species.
A fascinating aspect of Schwann cell biology is the bifurcation of the Schwann cell lineage during development to form two distinct populations of mature Schwann cells, non-myelinating (Remak) cells and myelinating Schwann cells. Our understanding of the gene expression differences between Remak and myelinating Schwann cells has until very recently remained rudimentary, based on small handful of cell-subtype specific markers such as GFAP and myelin proteins . Using single cell sequencing and genetic labeling, Gerber and colleagues have very recently provided the first transcriptome data for the developing Schwann cell lineage that includes separation of the Remak and myelinating Schwann cell populations . For years, despite advances in other aspects of the Schwann cell response to injury, how the roles of these two Schwann cell sub-types compare during peripheral nerve regeneration has remained elusive. Recent studies have provided the first insights into the morphological contributions of myelinating and non-myelinating Schwann cells after nerve transection, demonstrating that both myelinating and non-myelinating Schwann cells have the capacity to change their shape and transform into elongated repair cells and that repair cells in turn can shorten dramatically following axon regeneration and form myelin internodes . Applying single cell sequencing technologies akin to the approach of Suter et al. to the injured nerve’s myelinating and Remak populations will build upon available transcriptomes, and notably on the data we present here to further elucidate the transcriptional basis of Schwann cell transdifferentiation, heterogeneity of repair cells, maintenance of the repair cell state, and Schwann cell re-differentiation.
The transcriptomic data presented here provide a valuable resource to peripheral nerve biologists to fuel novel hypotheses and uncover new pathways underlying the Schwann cell injury response. In addition to the molecules highlighted above, we anticipate the discovery of to date un-appreciated pathways and proteins using our dataset that will enable mechanistic advances in our understanding of how these remarkable cells transform after nerve injury. Confirmation of our transcription-level changes at the protein level followed by knock-out or overexpression experiments has the potential to yield answers to outstanding questions in the field of Schwann cell biology. For example, to further our understanding of why the PNS regenerates so well, and to generate new molecular handles for both understanding CNS regeneration failure and designing strategies to overcome it.
Availability of data and materials
To ensure widespread distribution of our RNAseq dataset, we have incorporated these data in a user-friendly website previously created by Ben Barres’ lab (www.BrainRNAseq.org) that provides a simple and searchable platform for browsing gene expression data.
Central nervous system
Fragments per kilobase of transcript sequence per million mapped fragments
Gene Expression Omnibus
Kyoto Encyclopedia of Genes and Genomes
National Center for Biotechnology Information
Peripheral nervous system
Arthur-Farraj PJ, Latouche M, Wilton DK, Quintes S, Chabrol E, Banerjee A, et al. c-Jun reprograms Schwann cells of injured nerves to generate a repair cell essential for regeneration. Neuron. 2012;75(4):633–47. https://doi.org/10.1016/j.neuron.2012.06.021.
Arthur-Farraj PJ, Morgan CC, Adamowicz M, Gomez-Sanchez JA, Fazal SV, Beucher A, et al. Changes in the coding and non-coding transcriptome and DNA methylome that define the schwann cell repair phenotype after nerve injury. Cell Rep. 2017;20(11):2719–34. https://doi.org/10.1016/j.celrep.2017.08.064.
Baker MD. Electrophysiology of mammalian Schwann cells. Prog Biophys Mol Biol. 2002;78(2–3):83–103. https://doi.org/10.1016/s0079-6107(02)00007-x.
Blighe K, Rana S, Lewis M. EnhancedVolcano: publication-ready volcano plots with enhanced colouring and labeling. R package version 1.10.0. 2021.
Brosius Lutz A, Chung WS, Sloan SA, Carson GA, Zhou L, Lovelett E, et al. Schwann cells use TAM receptor-mediated phagocytosis in addition to autophagy to clear myelin in a mouse model of nerve injury. Proc Natl Acad Sci USA. 2017;114(38):E8072-e8080. https://doi.org/10.1073/pnas.1710566114.
Chen S, Zou Q, Chen Y, Kuang X, Wu W, Guo M, et al. Regulation of SPARC family proteins in disorders of the central nervous system. Brain Res Bull. 2020;163:178–89. https://doi.org/10.1016/j.brainresbull.2020.05.005.
De Logu F, Nassini R, Materazzi S, Carvalho Gonçalves M, Nosi D, Rossi Degl’Innocenti D, et al. Schwann cell TRPA1 mediates neuroinflammation that sustains macrophage-dependent neuropathic pain in mice. Nat Commun. 2017;8(1):1887. https://doi.org/10.1038/s41467-017-01739-2.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21. https://doi.org/10.1093/bioinformatics/bts635.
Egaña-Gorroño L, López-Díez R, Yepuri G, Ramirez LS, Reverdatto S, Gugger PF, et al. Receptor for advanced glycation end products (RAGE) and mechanisms and therapeutic opportunities in diabetes and cardiovascular disease: insights from human subjects and animal models. Front Cardiovasc Med. 2020;7:37. https://doi.org/10.3389/fcvm.2020.00037.
Figlia G, Norrmén C, Pereira JA, Gerber D, Suter U. Dual function of the PI3K-Akt-mTORC1 axis in myelination of the peripheral nervous system. Elife. 2017. https://doi.org/10.7554/eLife.29241.
Gan KJ, Südhof TC. SPARCL1 promotes excitatory but not inhibitory synapse formation and function independent of neurexins and neuroligins. J Neurosci. 2020;40(42):8088–102. https://doi.org/10.1523/jneurosci.0454-20.2020.
Gaudet AD, Popovich PG, Ramer MS. Wallerian degeneration: gaining perspective on inflammatory events after peripheral nerve injury. J Neuroinflamm. 2011;8:110. https://doi.org/10.1186/1742-2094-8-110.
Gerber D, Pereira JA, Gerber J, Tan G, Dimitrieva S, Yángüez E, Suter U. Transcriptional profiling of mouse peripheral nerves to the single-cell level to build a sciatic nerve ATlas (SNAT). Elife. 2021. https://doi.org/10.7554/eLife.58591.
Goldberg JL, Vargas ME, Wang JT, Mandemakers W, Oster SF, Sretavan DW, Barres BA. An oligodendrocyte lineage-specific semaphorin, Sema5A, inhibits axon growth by retinal ganglion cells. J Neurosci. 2004;24(21):4989–99. https://doi.org/10.1523/jneurosci.4390-03.2004.
Gomez-Sanchez JA, Carty L, Iruarrizaga-Lejarreta M, Palomo-Irigoyen M, Varela-Rey M, Griffith M, et al. Schwann cell autophagy, myelinophagy, initiates myelin clearance from injured nerves. J Cell Biol. 2015;210(1):153–68. https://doi.org/10.1083/jcb.201503019.
Gomez-Sanchez JA, Pilch KS, van der Lans M, Fazal SV, Benito C, Wagstaff LJ, et al. After nerve injury, lineage tracing shows that myelin and remak Schwann cells elongate extensively and branch to form repair Schwann cells, which shorten radically on remyelination. J Neurosci. 2017;37(37):9086–99. https://doi.org/10.1523/jneurosci.1453-17.2017.
Ishii A, Furusho M, Bansal R. Mek/ERK1/2-MAPK and PI3K/Akt/mTOR signaling plays both independent and cooperative roles in Schwann cell differentiation, myelination and dysmyelination. Glia. 2021;69(10):2429–46.
Jessen KR, Arthur-Farraj P. Repair Schwann cell update: adaptive reprogramming, EMT, and stemness in regenerating nerves. Glia. 2019;67(3):421–37. https://doi.org/10.1002/glia.23532.
Jessen KR, Mirsky R. Schwann cells: early lineage, regulation of proliferation and control of myelin formation. Curr Opin Neurobiol. 1992;2(5):575–81. https://doi.org/10.1016/0959-4388(92)90021-c.
Jessen KR, Mirsky R. Negative regulation of myelination: relevance for development, injury, and demyelinating disease. Glia. 2008;56(14):1552–65. https://doi.org/10.1002/glia.20761.
Jessen KR, Mirsky R. The repair Schwann cell and its function in regenerating nerves. J Physiol. 2016;594(13):3521–31. https://doi.org/10.1113/jp270874.
Kalinski AL, Yoon C, Huffman LD, Duncker PC, Kohen R, Passino R, et al. Analysis of the immune response to sciatic nerve injury identifies efferocytosis as a key mechanism of nerve debridement. Elife. 2020. https://doi.org/10.7554/eLife.60223.
Kucukdereli H, Allen NJ, Lee AT, Feng A, Ozlu MI, Conatser LM, et al. Control of excitatory CNS synaptogenesis by astrocyte-secreted proteins Hevin and SPARC. Proc Natl Acad Sci USA. 2011;108(32):E440-449. https://doi.org/10.1073/pnas.1104977108.
Langfelder P, Horvath S. Fast R Functions for Robust Correlations and Hierarchical Clustering. J Stat Softw. 2012;46(11):i11.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 2008;9:559. https://doi.org/10.1186/1471-2105-9-559.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011;12:323. https://doi.org/10.1186/1471-2105-12-323.
Liu JH, Tang Q, Liu XX, Qi J, Zeng RX, Zhu ZW, et al. Analysis of transcriptome sequencing of sciatic nerves in Sprague-Dawley rats of different ages. Neural Regen Res. 2018;13(12):2182–90. https://doi.org/10.4103/1673-5374.241469.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. https://doi.org/10.1186/s13059-014-0550-8.
Lutz AB. Purification of Schwann cells. Cold Spring Harb Protoc. 2014;2014(12):1234–6. https://doi.org/10.1101/pdb.top073981.
Ma CH, Omura T, Cobos EJ, Latrémolière A, Ghasemlou N, Brenner WCJ, et al. Accelerating axonal growth promotes motor recovery after peripheral nerve injury in mice. J Clin Invest. 2011;121(11):4332–47. https://doi.org/10.1172/JCI58675.
Meyer Zu Reckendorf S, Brand C, Pedro MT, Hegler J, Schilling CS, Lerner R, et al. Lipid metabolism adaptations are reduced in human compared to murine Schwann cells following injury. Nat Commun. 2020;11(1):2123. https://doi.org/10.1038/s41467-020-15915-4.
Nagarajan R, Le N, Mahoney H, Araki T, Milbrandt J. Deciphering peripheral nerve myelination by using Schwann cell expression profiling. Proc Natl Acad Sci USA. 2002;99(13):8998–9003. https://doi.org/10.1073/pnas.132080999.
Napoli I, Noon LA, Ribeiro S, Kerai AP, Parrinello S, Rosenberg LH, et al. A central role for the ERK-signaling pathway in controlling Schwann cell plasticity and peripheral nerve regeneration in vivo. Neuron. 2012;73(4):729–42. https://doi.org/10.1016/j.neuron.2011.11.031.
Narciso MS, Mietto Bde S, Marques SA, Soares CP, Mermelstein Cdos S, El-Cheikh MC, Martinez AM. Sciatic nerve regeneration is accelerated in galectin-3 knockout mice. Exp Neurol. 2009;217(1):7–15. https://doi.org/10.1016/j.expneurol.2009.01.008.
Newbern JM, Li X, Shoemaker SE, Zhou J, Zhong J, Wu Y, et al. Specific functions for ERK/MAPK signaling during PNS development. Neuron. 2011;69(1):91–105. https://doi.org/10.1016/j.neuron.2010.12.003.
Painter MW, Brosius Lutz A, Cheng YC, Latremoliere A, Duong K, Miller CM, et al. Diminished Schwann cell repair responses underlie age-associated impaired axonal regeneration. Neuron. 2014;83(2):331–43. https://doi.org/10.1016/j.neuron.2014.06.016.
Reichert F, Saada A, Rotshenker S. Peripheral nerve injury induces Schwann cells to express two macrophage phenotypes: phagocytosis and the galactose-specific lectin MAC-2. J Neurosci. 1994;14(5 Pt 2):3231–45. https://doi.org/10.1523/jneurosci.14-05-03231.1994.
Senatus LM, Schmidt AM. The AGE–RAGE axis: implications for age-associated arterial diseases. Front Genet. 2017;8:187. https://doi.org/10.3389/fgene.2017.00187.
Tasdemir-Yilmaz OE, Druckenbrod NR, Olukoya OO, Dong W, Yung AR, Bastille I, Segal RA. Diversity of developing peripheral glia revealed by single-cell RNA sequencing. Dev Cell. 2021;56(17):2516-2535.e8. https://doi.org/10.1016/j.devcel.2021.08.005.
Toma JS, Karamboulas K, Carr MJ, Kolaj A, Yuzwa SA, Mahmud N, et al. Peripheral nerve single-cell analysis identifies mesenchymal ligands that promote axonal growth. eNeuro. 2020. https://doi.org/10.1523/eneuro.0066-20.2020.
Ulgen E, Ozisik O, Sezerman OU. pathfindR: an R package for comprehensive identification of enriched pathways in omics data through active subnetworks. Front Genet. 2019;10:858. https://doi.org/10.3389/fgene.2019.00858.
Wada R, Yagihashi S. Role of advanced glycation end products and their receptors in development of diabetic neuropathy. Ann N Y Acad Sci. 2005;1043:598–604. https://doi.org/10.1196/annals.1338.067.
Weiss T, Taschner-Mandl S, Bileck A, Slany A, Kromp F, Rifatbegovic F, Frech C, Windhager R, et al. Proteomics and transcriptomics of peripheral nerve tissue and cells unravel new aspects of the human Schwann cell repair phenotype. Glia. 2016;64(12):2133–53. https://doi.org/10.1002/glia.23045.
Welleford AS, Quintero JE, Seblani NE, Blalock E, Gunewardena S, Shapiro SM, et al. RNA sequencing of human peripheral nerve in response to injury: distinctive analysis of the nerve repair pathways. Cell Transplant. 2020;29:963689720926157. https://doi.org/10.1177/0963689720926157.
Wilson GF, Chiu SY. Ion channels in axon and Schwann cell membranes at paranodes of mammalian myelinated fibers studied with patch clamp. J Neurosci. 1990;10(10):3263–74. https://doi.org/10.1523/jneurosci.10-10-03263.1990.
Wolbert J, Li X, Heming M, Mausberg AK, Akkermann D, Frydrychowicz C, et al. Redefining the heterogeneity of peripheral nerve cells in health and autoimmunity. Proc Natl Acad Sci USA. 2020;117(17):9466–76. https://doi.org/10.1073/pnas.1912139117.
Wu W, Liu Y, Wang Y. Sam68 promotes Schwann cell proliferation by enhancing the PI3K/Akt pathway and acts on regeneration after sciatic nerve crush. Biochem Biophys Res Commun. 2016;473(4):1045–51. https://doi.org/10.1016/j.bbrc.2016.04.013.
Xu X, Gao W, Li L, Hao J, Yang B, Wang T, et al. Annexin A1 protects against cerebral ischemia-reperfusion injury by modulating microglia/macrophage polarization via FPR2/ALX-dependent AMPK-mTOR pathway. J Neuroinflamm. 2021;18(1):119. https://doi.org/10.1186/s12974-021-02174-3.
Yi S, Zhang H, Gong L, Wu J, Zha G, Zhou S, et al. Deep sequencing and bioinformatic analysis of lesioned sciatic nerves after crush injury. PLoS ONE. 2015;10(12):e0143491. https://doi.org/10.1371/journal.pone.0143491.
Zhang Y, Sloan SA, Clarke LE, Caneda C, Plaza CA, Blumenthal PD, et al. Purification and characterization of progenitor and mature human astrocytes reveals transcriptional and functional differences with mouse. Neuron. 2016;89(1):37–53. https://doi.org/10.1016/j.neuron.2015.11.013.
We would like to thank past Barres Lab members for helpful comments and ideas concerning this work.
This work was supported by the Dr. Miriam and Sheldon G. Adelson Medical Research Foundation (B.A.B.), Christopher and Dana Reeve Foundation (B.A.B.), NIH Grants EY11310, NS069375, 1S10RR02678001, and NIH S10 Shared Instrument Grant S10RR027431-01 (to B.A.B.), NIH Grant T32 HD007249 Developmental and Neonatal Biology Training Program (to A.B.L.), National Institute of Mental Health Grants T32GM007365, F30MH106261, R01MH125956 (to S.A.S) Bio-X Fellowship (to S.A.S.), the NARSAD Brain and Behavior Young Investigator Award (to S.A.S.) and the Stanford Medical Scientist Training Program and School of Medicine (A.B.L.). M.S.B received funding from an AHA/Allen Brain Health Award (19PABHI34580007) and Foundation Leducq (19CVD01).
Ethics approval and consent to participate
All procedures involving animals were conducted in conformity with Stanford University guidelines and are in compliance with national and state laws and policies.
Consent for publication
The authors declare no conflicts of interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
PCA of all RNA-seq samples. (A) PCA of all samples colored by day post-crush. Circles indicate whole nerve samples and triangles represent Schwann cell purified samples. (B) PCA of Schwann cell samples alone, colored by time point.
Generation of in situ hybridization probes. RNA probes used for in situ hybridization validation. Expected probe lengths are listed beneath each column.
qRT-PCR primer sequences (rat) used to validate RNAseq findings.
All FPKM data for uncrushed (day 0) and crushed (days 3, 5, 7 post-crush) samples. This includes whole nerve and purified Schwann cells. Also included is dissociated sciatic nerve to help identify genes that change from the dissociation process, itself.
All differentially expressed genes and their p-values (and p-adjusted) in Schwann cells at days 3, 5, and 7 post-crush, as compared to day 0 (uncrushed).
Table of enriched gene pathways in purified Schwann cells at days 3, 5, and 7. First two sheets show pathways across all timepoints and the subsequent sheets highlight statistics of each pathway separated by timepoints.
Table of downregulated gene pathways in purified Schwann cells at days 3, 5, and 7.
Module assignments (from the WGCNA) for each gene.
Analysis of markers of different S100-expressing cell populations.
About this article
Cite this article
Brosius Lutz, A., Lucas, T.A., Carson, G.A. et al. An RNA-sequencing transcriptome of the rodent Schwann cell response to peripheral nerve injury. J Neuroinflammation 19, 105 (2022). https://doi.org/10.1186/s12974-022-02462-6
- Schwann cell
- Repair cell
- Injury response
- Peripheral nerve