Brainstem-Evoked Transcription of Defensive Genes After Spinal Cord Injury

The spinal cord after injury shows altered transcription in numerous genes. We tested in a pilot study whether the nucleus raphé magnus, a descending serotonergic brainstem region whose stimulation improves recovery after incomplete spinal cord injury (SCI), can influence these transcriptional changes. Rats received 2 h of low-frequency electrical stimulation in the raphé magnus 3 days after an impact contusion at segment T8. Comparison groups lacked injuries or activated stimulators or both. Immediately following stimulation, spinal cords were extracted, their RNA transcriptome sequenced, and differential gene expression quantified. Confirming many previous studies, injury primarily increased inflammatory and immune transcripts and decreased those related to lipid and cholesterol synthesis and neuronal signaling. Stimulation plus injury, contrasted with injury alone, caused significant changes in 43 transcripts (39 increases, 4 decreases), all protein-coding. Injury itself decreased only four of these 43 transcripts, all reversed by stimulation, and increased none of them. The non-specific 5-HT7 receptor antagonist pimozide reversed 25 of the 43 changes. Stimulation in intact rats principally caused decreases in transcripts related to oxidative phosphorylation, none of which were altered by stimulation in injury. Gene ontology (biological process) annotations comparing stimulation with either no stimulation or pimozide treatment in injured rats highlighted defense responses to lipopolysaccharides and microorganisms, and also erythrocyte development and oxygen transport (possibly yielding cellular oxidant detoxification). Connectivity maps of human orthologous genes generated in the CLUE database of perturbagen-response transcriptional signatures showed that drug classes whose effects in injured rats most closely resembled stimulation without pimozide include peroxisome proliferator-activated receptor agonists and angiotensin receptor blockers, which are reportedly beneficial in SCI. Thus the initial transcriptional response of the injured spinal cord to raphé magnus stimulation is upregulation of genes that in various ways are mostly protective, some probably located in recently arrived myeloid cells.

Introduction

Mechanical SCI causes numerous molecular, cellular and structural changes that vary in duration and onset, both at the site of injury and beyond (Dulin et al., 2015; O’Shea et al., 2017). These include changes in gene transcription that collectively and individually could influence recovery via protein synthesis and downstream molecules. It is therefore of interest to find ways to modify influential transcriptional changes by enhancing or opposing them according to their benefit or harm. Here we explore whether activity in a raphé nucleus of the hindbrain, the NRM, regulates transcription after SCI. The NRM is one of several ventral medullary nuclei that send overlapping serotonergic projections (as well as non-serotonergic projections) directly to the dorsal and ventral spinal cord, others being the nucleus raphe obscurus, nucleus raphe pallidus and the laterally adjacent paragigantocellularis and gigantocellularis nuclei (Kwiat and Basbaum, 1992; Gautier et al., 2017).

Previously, we showed that prolonged electrical stimulation of the NRM (for several days to weeks) significantly enhanced recovery from various autonomic, sensory, and motor deficits of contusional cervical or thoracic SCI in rats, while improving myelination and serotonergic innervation near the lesion (Hentall and Burns, 2009; Hentall and Gonzalez, 2012; Vitores et al., 2018). A shorter period of NRM stimulation (2 h) was found to restore injury-depleted levels in the spinal cord of cAMP and its downstream targets, phosphorylated cAMP response element-binding protein, and pPKA (Carballosa-Gonzalez et al., 2014), which cause neurotrophic effects via alterations in gene expression (Hannila and Filbin, 2008; Carballosa-Gonzalez et al., 2014; Batty et al., 2017). Stimulation of the midbrain’s median and dorsal raphé, which send ascending projections throughout the forebrain, analogously yielded partial restoration of anatomic and functional deficits following TBI in rats, as well as normalizing cAMP levels, thus establishing the generality of the idea that raphé nuclei provide central links for restorative feedback. In mice, prolonged NRM stimulation attenuated signs of experimental autoimmune encephalomyelitis, a model of MS, including cellular pathology and cytokine upregulation (Madsen et al., 2017), which further supports this idea. Our general aim in the studies cited above has been to achieve lasting, clinically significant recovery with a period of prolonged electrical DBS by evoking neurotrophic or protective effects from activated neuronal systems. In contrast, most studies of DBS for neurodegenerative conditions such as SCI, TBI, and MS have focused on the acute reversal of overt problematic signs or symptoms (e.g., pain, paralysis, altered mental status), without necessarily remedying the underlying pathology (Roy and Aziz, 2014; Shin et al., 2014; Chari et al., 2017).

Here, we examine transcriptional changes produced by NRM stimulation in rats with thoracic (T8) SCI. Stimulation was applied for 2 h to the NRM on the third day after a weight-drop SCI at segment T8. At this post-injury stage, many molecular pathways that promote plasticity and inflammation are near their peak (David and Kroner, 2011; Mao et al., 2016). The relatively brief period of stimulation, lasting just long enough for initial transcriptional effects to emerge, was followed by immediate extraction of various portions of the spinal cord and standard processing for RNA sequence analysis. Both injured and non-injured controls were studied with and without stimulation (n ≥ 3). Additional injured rats (n = 2) received stimulation after being pretreated with the non-specific 5-HT7 antagonist pimozide, to assess possible participation of serotonin release from the axon terminals of NRM neurons. The results proved surprising. Our working hypothesis was that stimulation in injured rats would induce neurotrophic transcriptional effects and predominantly reverse injury-produced changes. Instead it altered a small number of genes that are mostly concerned with inflammation and erythrocyte formation, very few of which were altered by injury alone. These findings shed new light on the brainstem’s descending modulatory influence on endogenous processes in healthy and injured spinal cord, and also identified molecules with potential for improving outcomes after SCI.

Materials and Methods

Animal Surgical Procedures and Treatments

All experiments were performed in accordance with the guidelines of the NIH Guide for the Care and Use of Laboratory Animals, and were approved by the University of Miami Miller School of Medicine Institutional Animal Care and Use Committee. Female young-adult Sprague-Dawley rats (220–240 g, 10–12 weeks old), obtained from Harlan Sprague-Dawley, Inc. (Indianapolis, IN, United States), were anesthetized with intraperitoneal ketamine (50 mg/kg) plus xylazine (10 mg/kg) and mounted in a stereotaxic head holder. Following a midline back incision and lateral dissection of the spinous musculature, a T8 laminectomy was performed and a moderate bilateral contusion injury was caused with a NYU-MASCIS Impactor, which applied a force by letting a 10 g rod of 2 mm diameter drop 12.5 mm centered on the midline of the spinal cord (Kearney et al., 1988). Some control animals received a laminectomy without the contusion injury. Dissected muscle layers were sutured and the skin was closed with wound clips. Animals recovered on a 37°C heating blanket. The opioid buprenorphine (0.01 mg/kg bid, subcutaneous) was given daily for analgesia and the antibiotic gentamycin (0.01 mg/kg bid, subcutaneous) to prevent infection. Bladder volume was checked daily and manually emptied.

At 72 h after the injury or sham injury, the animals were anesthetized with isoflurane (1.2% in oxygen) by face mask after induction in a glass chamber and mounted in a stereotaxic holder. A single monopolar stimulating microelectrode was placed in the brainstem by making a linear skin incision and drilling a 1.8 mm craniostomy on the midline, 2.2 mm caudal to the interaural line. To apply stimulation, a tungsten microelectrode (AC impedance 0.5 megohm, diameter 0.13 mm) was inserted into the midline NRM at the stereotaxic coordinates 2.2 mm caudal and 10 mm ventral to the interaural line (flat-skull orientation). Accuracy was confirmed by observing bilateral facial twitching 1.2–1.5 mm above the target, which disappeared once the electrode reached the NRM. This localization method was validated histologically in a previously published study (Hentall and Gonzalez, 2012). Monopolar stimulation, consisting of 8-Hz trains of cathodal pulses (30 μA, 1 ms), was given for 2 h in a 5 min, 50% duty cycle, ending in the stimulation phase. Electrodes were also placed in sham and injured control animals that did not receive the 2 h of stimulation. One injured group received pimozide by intraperitoneal injection of 1 mg/kg, 1 h prior to stimulation onset. Immediately after the stimulation period, the subdural spinal cord (including pia-arachnoid) was rapidly dissected and transversely sectioned into three 0.5 cm blocks encompassing the C4–C6, T7–T9, and L1–L3 spinal segments. The blocks were immediately placed in liquid nitrogen and stored at −80°C. Animals were euthanatized by isoflurane overdose. No adverse effects of stimulation or other treatment were seen.

RNAseq Sample Preparation

Cervical, thoracic and lumbar regions were combined for RNA analysis. Preparation and sequencing of RNA libraries was carried out by personnel of the John P. Hussman Institute for Human Genomics at the Center for Genome Technology, University of Miami Miller School of Medicine. RNA isolated from rat spinal cords was extracted using Trizol and purified using the QIAGEN (Venlo, Netherlands) RNeasy kit, by following the kit’s standard protocol. Samples were analyzed with an Agilent Bioanalyzer (Santa Clara, CA, United States). All analyzed samples passed minimum quality and quantity threshold required for RNA sequencing, including RNA integrity score above 7.0. Mean concentration of total RNA were 0.43 μg/μl (SD: 0.15); the mean purity (260/280 ratio) was 2.06 (SD: 0.041) and was > 2.02 in all individuals. From each sample, 500 ng of total RNA was provided to an Illumina TruSeq Stranded Total RNA Library Prep Kit with Ribo-Zero (San Diego, CA, United States) to create ribosomal RNA-depleted sequencing libraries with a unique barcode. Sequencing was performed to >25 million raw reads in a single-end 75 bp sequencing run on the Illumina NextSeq500.

Production of RNAseq Data

Raw sequence data was processed by the on-instrument Real Time Analysis software (v.2.7.7) transferred to de-multiplexed FASTQ files with the Illumina-supplied scripts in BCL2FASTQ software (v2.17). The quality of the reads was determined with FASTQC software¹ for per base sequence quality, duplication rates, and overrepresented k-mers. Illumina sequencing adapters were trimmed from the ends of the reads using the Trim Galore! Package², then aligned for Rattus norvegicus (Rnor_5.0) with the STAR aligner (v2.5.0a) (Dobin et al., 2013). Gene count quantification for aligned reads was performed using the GeneCounts function within STAR against the Ensembl gene build pipeline (release 79).

Data Analysis

Experimental groups, to which rats were assigned in random order, were as follows: group A, sham injury and sham stimulation; group B, SCI and sham stimulation; group C, SCI and stimulation; group D, SCI with stimulation and pimozide; group E, sham injury with stimulation. Many subjects had also provided material for a previously published assay of cAMP, cAMP response element-binding protein and protein kinase A (Carballosa-Gonzalez et al., 2014); exceptions were one subject in group E and both subjects in group D, whose samples were contemporaneous with the rest but had not previously been assayed for publication. Within 5 experimental groups, there were 10 possible contrasts, of which 6 were analyzed (Table 1). These contrasts are named in the text by the positively changed factor followed in parentheses by the constant background. Thus, ±st(SCI) refers to the contrast of stimulation (st) with no stimulation in injured rats, ±SCI(st) contrasts injury with no injury in stimulated rats, ±pim(st/SCI) contrasts pimozide’s presence with its absence in stimulated, injured rats; ±st(0), ±SCI(0), and ±st/SCI(0) respectively contrast stimulation, injury or both with their absence (in rats with no other interventions).

TABLE 1

Table 1. Groups and contrasts analyzed.

Differential gene expression analysis was performed on the gene count data with edgeR software (Robinson et al., 2010). Gene counts were normalized against total aligned reads to generate counts-per-million for each gene in each sample. Given the relatively small sample sizes per group, the exact test implemented in EdgeR was used to determine differential expression, producing a false discovery rate (FDR) p-value. Significant differential expression was decided by EdgeR with joint criteria of FDR < 0.05 and a greater than 100% change: abs(log₂(FC)) > 1, where FC is FC. For comparison, output from two other DGE software packages, DESeq2 and baySeq, was also examined (Hardcastle and Kelly, 2010; Love et al., 2014). Numbers of significant increases and decreases obtained in different contrasts were analyzed in 2 × 2 contingency tables with the two-tailed Chi-square test.

Pathway and network analyses were performed with Ingenuity Pathway Analysis (QIAGEN), using as input the lists of genes with FDR < 0.05, as determined by EdgeR. Gene ontology (biological process) analysis was performed for each set of genes meeting EdgeR cutoff criteria of abs(log₂(FC)) > 1 and FDR < 0.05 with the DAVID³. Statistical overrepresentation was determined with the default settings in DAVID. Gene names used in tables were official gene symbols given by DAVID, supplemented by Gene from the United States’ NCBI.

Human orthologs of rat genes flagged in the various EdgeR contrasts were obtained from the Ensembl 2018 website (Zerbino et al., 2018). The resulting human gene list was applied to the L1000 gene expression database of the CLUE software platform (Subramanian et al., 2017). This platform produced a connectivity map that was analyzed with the Query tool to find classes of perturbagens (drugs) whose expression signatures most closely resembled those yielded by the experiments.

Data Availability

The high-throughput sequence data is available from the GEO repository of the NCBI, accession number GSE133093. Other data relevant to the study’s conclusions are available on request from the corresponding author.

Results

Numbers Flagged as Upregulated or Downregulated

A total of 32,494 transcripts with unique rat gene identification (Ensembl) numbers were analyzed. These included 1,563 micro RNAs, 1,471 small RNAs, and 22,019 experimentally confirmed protein-coding RNAs; further details are given in the Supplementary Figure 1. For all six contrasts, the DGE software applications EdgeR and DESEq2 were of roughly equal leniency, and showed mean overlap of 75% for selected transcripts. The baySeq application was considerably more stringent, flagging on average 45% of those flagged by the other two software applications (Supplementary Figure 1). The following presentation is limited to results obtained with EdgeR. In the key contrasts ±st(SCI) and ±st(0), neither DESeq2 nor baySeq found any significant genes that were not flagged by edgeR.

The contrasts ±SCI(0), ±st/SCI(0) and ±SCI(st), all of which pair intact rats with rats injured 3 days previously, revealed a large number of differentially expressed transcripts. These numbers differed significantly (pairwise Chi-squared contingency tables, p < 0.0001). Contrast ±SCI(0) showed fewer differentially transcribed genes (n = 982) than ±st/SCI(0) (n = 1150), and ±SCI(st) showed the most differentially expressed transcripts (n = 1374). The relatively brief interventions of 2 h of stimulation or 3 h of pimozide’s presence caused far fewer significant alterations in transcription. The numbers of genes affected in the contrast ±pim(st/SCI) (n = 81) was significantly different (p = 0.0009) from the numbers for ±st(SCI) (n = 43), but the latter was not significantly different (p = 0.10) from the contrast ±st(0) (n = 28). With the brief interventions, there was also less variation in the difference between counts for the paired conditions, as reflected in their standard deviations (Table 1). All three interventions (injury, stimulation, pimozide) consistently produced more increases than decreases in transcript counts (Table 1): pimozide [contrast ±pim(st/SCI)] revealed 64% increases, injury alone [contrast ±SCI(0)] revealed 94% increases and stimulation after injury [contrast ±st(SCI)] revealed 91% increases. Differing kinetics of decay and synthesis may in some cases have contributed to this imbalance (see section “Discussion”).

Genes Affected by Stimulation: Responses to Injury Alone

The contrast showing the effect of stimulation in injury, contrast ±st(SCI), was of greatest importance for the goals of this study. It found that 43 transcripts, all protein-coding, had changed significantly (Table 2). When one of these 43 transcripts was flagged in another contrast involving injury, the direction of the effect was consistent: for example, in the opposite direction of the effects from injury alone [contrast ±SCI(0)] and in the same direction as the effect of injury plus stimulation compared with untreated rats [±st/SCI(0)] and as the effect of injury versus sham on stimulated rats [±SCI(st)] (Table 2). However, it was notable that only 9% of injury-produced transcriptional changes [contrast ±SCI(0)], all downregulation, were reversed by NRM stimulation [contrast ±st(SCI)].

TABLE 2

Table 2. Genes altered by stimulation in injured animals, selected by contrast ±st(SCI).

The effect of stimulation on injury could also be shown less directly, by comparing genes flagged as altered by stimulation plus injury, in contrast ±st/SCI(0), with those flagged by injury alone, in contrast ±SCI(0). Some genes (n = 315) showed altered transcription after stimulation plus injury, in contrast ±st/SCI(0), but were not affected by injury alone, in contrast ±SCI(0). Most genes (n = 835) were flagged by both contrasts, and some (n = 147) were flagged by ±SCI(0) but not ±st/SCI(0). As seen in Table 2, of the 43 genes selected by contrast ±st(SCI), 15 were flagged by ±st/SCI(0) but not by ±SCI(0), 3 were flagged by ±SCI(0) but not by ±st/SCI(0) and one gene (Mmmp13) was flagged by both contrasts. This numerical imbalance confirms that stimulation caused changes in transcription mainly in genes that were unaffected by injury alone.

Effect of Pimozide on the Response to Stimulation

Pimozide significantly altered 58% of the genes affected by stimulation after injury. Despite the possibility of bias for increases over decreases described above and the low sample number (n = 2 rats), genes flagged in contrast ±pim(st/SCI) showed a consistent direction of effect on genes flagged in contrast ±st(SCI) (Table 2). That is, upregulation was reversed in 24 genes and downregulation was reversed in one gene.

Effect of Injury on the Response to Stimulation

Among rats that received stimulation, comparison of injured with intact groups, in contrast ±SCI(st), yielded further insight. A large number (n = 1374) of genes was flagged by ±SCI(st), most of which (36 of 43) were also flagged by contrasting stimulation with no stimulation, ±st(SCI) (Table 2). Injury thus profoundly alters the response to stimulation. Confirming this conclusion, the transcriptional response of intact rats to stimulation, as determined by contrast ±st(0), yielded a list of 28 flagged genes (Table 3), none of which was among those significantly affected by stimulation in injured rats [±st(SCI)]. Many genes flagged by ±st(0) were not protein-coding: 10 were genes for small RNAs. Unlike in injured rats, significant effects of stimulation in intact rats were in the same direction as the effects (when present) of injury and pimozide (Table 3).

TABLE 3

Table 3. Genes selected by contrast ±st(0) (effect of stimulation in non-injured rats).

Genes of Interest

Individual genes of greatest interest for understanding the effect of NRM stimulation on SCI are those flagged by all of the three contrasts that involved stimulation (without pimozide) in injury: ±st(SCI), ±SCI(st) and ±st/SCI(0). Twelve of these genes were increased by stimulation: Mmp13, Defa5, Ppbp, Mpo, Rhd, S100a8, S100a9, Camp, Klf1, Robo-1, Timd2, and Fcnb. Of these 12 genes, 7 also passed the considerably more restrictive filtering criterion of FDR < 0.00001, and all except one of these (Camp) was reversed by pimozide (Table 2). Three genes that were flagged in these three contrasts were decreased by stimulation: Slc17a7, Gbp6 and the uncharacterized AABR07043200.1, of which only the decrease in Gbp6 was reversed by pimozide. Figure 1 shows the normalized counts for all 5 experimental conditions of the key characterized genes that were upregulated or down-regulated.

FIGURE 1

Figure 1. Normalized RNA molecular counts from individual rats for 14 characterized genes identified as prominently up-regulated (n = 12) or down-regulated (n = 2) by stimulation. The included genes were those flagged by all three contrasts ±st(SCI), ±SCI(st) and ±st/SCI(0), as listed in Table 2. The vertical axis is scaled logarithmically and log zero values (equal to –∞) are presented separately off-scale.

Few sRNA and miRNA genes underwent significantly altered transcription under any experimental condition. Mir100 was down-regulated in the contrasts ±SCI(0) and ±SCI(st). Mir3593 was down-regulated in the contrast ±SCI(st). Finally Mir6326, Mir675, and Mir21 were up-regulated in the contrast ±SCI(st).

Inferred Pathways and Biological Processes

Injury primarily increased inflammatory and immune transcripts and decreased those related to lipid and cholesterol synthesis and neuronal signaling, confirming previous reports (Bareyre and Schwab, 2003; Chamankhah et al., 2013; Shi et al., 2017). This was seen in the top canonical pathways determined with Ingenuity Pathway Analysis (Table 4) and in DAVID gene ontologies (biological process) for the contrast ±SCI(0) (some of which are listed in Table 5). Stimulation after injury had a major effect on pathways concerned with defense responses to lipopolysaccharides and various microorganisms (fungi, bacteria, yeast, viruses) and related inflammatory and immune responses (Tables 4, 5). Many of these same pathways were among the principal targets of pimozide. Also prominent among the pathways affected by stimulation in injury were erythrocyte development, oxygen transport and heme biosynthesis (Tables 4, 5). On the other hand, stimulation without injury, contrast ±st(0), had its largest effects on the canonical pathways for oxidative phosphorylation and mitochondrial dysfunction, with little influence on inflammatory and immune processes (Table 4). The DAVID gene ontologies (biological process) for contrast ±st(0) found that only two terms reached (weak) significance: “response to hyperoxia” and “regulation of apoptotic process.”

TABLE 4

Table 4. Top canonical pathways indicated by IPA for the effects of injury alone [±SCI(0)], of stimulation after injury [±st(SCI)], stimulation on intact rats [±st(0)], and of pimozide on injured, stimulated rats [±pim(st/SCI)].

TABLE 5

Table 5. Gene Ontology terms (biological process, GOTERM_BP_DIRECT) found to be enriched in by stimulation after injury [±st(SCI)] or its modification by pimozide [±pim(st/SCI)].

Resemblances of Expression Patterns to Perturbation With Drug Classes

The L1000 gene expression database of the CLUE software platform (clue.io) (Subramanian et al., 2017) was used to ascertain drug classes whose actions on cell cultures resembled to some quantifiable degree the effect of NRM stimulation on SCI. Human orthologs of the rat genes flagged in the two EdgeR contrasts ±st(SCI) and ±pim(st/SCI) were obtained from the Ensembl 2018 website (Zerbino et al., 2018). From the contrast ±st(SCI), which had flagged 43 rat genes, 37 human orthologous genes could be identified, and from the contrast ±pim(st/SCI), which originally yielded 81 rat genes, 65 human orthologous genes were identified (Supplementary Table 1). These two human gene lists were input to the CLUE Query tool to probe the database connectivity map, constructed from genetic changes in various human cell culture lines in response to diverse drugs. The classes of perturbagens (drugs) whose expression signatures most closely resembled those obtained in the experiments are listed in Table 6. The similarity score for the ±pim(st/SCI) contrast was reversed so that the absence of the drug pimozide was designated as the positive effect. High similarity (>98%) was reached in both contrasts by four perturbagen classes: PPAR agonists, angiotensin receptor antagonists, glucocorticoid receptor agonists, and nucleoporin loss of function (LOF).

TABLE 6

Table 6. Perturbagens (pharmacologic and genetic) whose expression signatures most closely resemble the effect of stimulation after injury and its modification by pimozide.

Discussion

General Observations and Technical Considerations

The effect of the NRM, or indeed of any brainstem system, on transcription in the spinal cord has not previously been reported. Here we show that transcriptional changes occurred in the first 2 h of NRM stimulation in a small fraction of the total genome. A longer period of stimulation probably could have produced more changes, as seen with the 3-day injury period and even with the 3-h presence of pimozide. The effect of the NRM stimulation differed completely between injured and normal spinal cord. The changes due to stimulation were for the most part not reversals (or enhancements) of those caused by injury alone, and thus implicate a distinct set of genes. In SCI, stimulation affected primarily transcription in genes with defensive immune functions. Other notable changes were in genes involved in erythrocyte development, such as constituent molecules of hemoglobin. There was a paucity of neuron-specific effects. These general findings lead to a somewhat revised and more complex model of the NRM’s repair effect, with defensive and protective processes preceding (and perhaps permitting or evoking) the regenerative benefits of stimulation.

The validity of the findings and their interpretation must take into account the low numbers in each group. Low numbers are not uncommon in RNAseq studies, typically imposed by costs and by other external factors. A study of biological replicate numbers with various DGE software applications (Schurch et al., 2016), found that EdgeR, DESeq2, and baySeq can all produce satisfactory and mutually consistent conclusions with as few as three replicates, given a < 5% false discovery rate and a fold change of ±50% [or abs(log2(FC)) > 0.58]. Here, we applied a more restrictive criterion for fold change of ±100% [abs(log₂(FC)) > 1] along with the < 5% false positive rate, and all groups except pimozide treatment consisted of > = 3 rats. Further strengthening the conclusions, the core list of altered genes (Figure 1, see section “Genes of Interest”) was derived from three paired contrasts that utilized all but two of the studied animals (n = 15). Thus the early changes in gene transcription found with NRM stimulation in injured animals can be specified with reasonable certainty. However, the findings strictly apply to female rats only, under conditions of treatment with buprenorphine, gentamycin and halothane, and the possibility remains of different genes emerging from contrasts made under other conditions or in other types of subject.

Pimozide is an antipsychotic drug that acts as an antagonist with a high affinity for both 5-HT7 receptors and for dopamine D2 subtype receptors (Roth et al., 1994). In the two subjects tested with pimozide, the drug consistently blocked and never enhanced the effects of stimulation in injury on individual genes. The descending pathways from the NRM are strongly serotonergic whereas dopaminergic projections are not prominent in the spinal cord and their cell bodies are distant to the site of stimulation (van Dijken et al., 1996). Thus the results with pimozide suggest that the effect of NRM stimulation after SCI was mediated in many genes by serotonin release, while a role for dopamine release is unlikely. The relatively rapid effect of pimozide in reversing gene expression induced by 2 h of NRM stimulation in SCI is likely to have occurred at the level of post-synaptic 5-HT7 receptors in the spinal cord. These receptors have been implicated in anti-inflammation, repair, and neuroplasticity at various CNS sites including the spinal cord (Kvachnina et al., 2005; de las Casas-Engel et al., 2013; Carballosa-Gonzalez et al., 2014; Di Pilato et al., 2014; Volpicelli et al., 2014; Fields and Mitchell, 2015). However, an action of pimozide at other CNS sites or other organs, which then indirectly causes a block of the influence of NRM stimulation on RNA synthesis in the spinal cord appears unlikely.

Measurements made after the 2 h time window of stimulation or after 3 days of injury represent snapshots of dynamic processes, which may include transient or long-lasting or even multiphasic changes in messenger RNA level. Changes in transcription rate are thought to determine the direction of most effects, because degradation rates and post-transcriptional processing are constant for most individual genes (Rabani et al., 2011, 2014), although inflammatory and immune signaling genes and targets of NFκB signaling are more likely to show modulated degradation (Rabani et al., 2011). The half-life of decay is typically slower than the transcription rate, which suggests the possibility of bias in favor of detecting increases over decreases. The half-life of messenger RNA in mammalian cell cultures averages 7–9 h, with a lower median, e.g., 4–6 h, and regulators of transcription and signal transduction tend to have shorter half-lives (Friedel et al., 2009; Schwanhausser et al., 2011). Therefore, given typical half-lives, the finding that 94% of changes were increases after the 3 day post-injury period probably reflects true differences in the direction of transcription rates. However, the finding that 91% of changes after the 2 h stimulation window were increases is alternatively explained by the slowness of decay.

Background Studies

This work is part of a series exploring the central control of endogenous repair mechanisms. It is proposed that the various brainstem raphé nuclei, including the NRM, form the central links of feedback loops controlling repair. In our model, raphé neurons receive sensory, chemical and central-state signals indicative of recent injury (e.g., pain, vestibular input, hypothermia, circulating cytokines). Their evoked electrical activity then causes axonal release of serotonin and co-release of neuropeptide transmitters from terminals that, collectively, innervate almost the entire CNS (Molliver, 1987; Hornung, 2003). Trophic and protective cellular effects can be produced by these various neuropeptides: thyrotropin releasing hormone (Daimon et al., 2013), substance P (Jiang et al., 2012; Kim et al., 2015), galanin (Hobson et al., 2008), met-enkephalin (Narita et al., 2006; Campbell et al., 2012) as well as by serotonin (Azmitia, 2007; Trakhtenberg and Goldberg, 2012), all acting through non-synaptic volume transmission (Ridet et al., 1993; Agnati et al., 1995; Borroto-Escuela et al., 2015). As in the response to infection which is frequently concurrent with injury, an endogenous repair system from the perspective of natural selection should promiscuously exploit (facilitate or depress) multiple molecular and cellular mechanisms.

The main experimental evidence for this model is that some days or weeks of stimulation in a raphé nucleus starting in the first days following incomplete contusional SCI or fluid-pressure TBI or induction of MS-like signs was found to improve behavioral signs over several weeks (Hentall and Gonzalez, 2012; Carballosa Gonzalez et al., 2013; Madsen et al., 2017). Concomitant histological changes seen at chronic endpoints following stimulation include increased myelination and serotonergic innervation around the injury after T8 SCI (Hentall and Gonzalez, 2012), increased calcitonin gene related peptide after C5 SCI (Vitores et al., 2018), restored volume loss in neocortex after TBI (Carballosa Gonzalez et al., 2013), and lessened cellular pathology and cytokine production in experimental autoimmune encephalomyelitis (Madsen et al., 2017).

The mechanisms leading to these long-term changes are likely to involve complexly interwoven cellular and molecular processes. To better understand the proximate causes of improved recovery with NRM stimulation, this and some other studies have focused on more acute, post-injury endpoints and shorter stimulation periods. An immunostaining study of 3-day C5 injury after 2 days of NRM stimulation showed changes in numbers of inflammatory cell types and a transition from neural precursors to radial glia that facilitate differentiation (Jermakowicz et al., 2019). In rats with 3-day T8 contusions, 2 h of NRM stimulation, as in the present protocol, restored normal levels of injury-depleted cAMP, phosphorylated cAMP response binding element (pCREB), and pPKA (Carballosa-Gonzalez et al., 2014). The same study showed that NRM stimulation had no significant effects on levels of cAMP, pPKA or pCREB in intact animals, but pimozide in intact animals with or without stimulation lowered cAMP to injured levels. Reversal of cAMP depletion has also been shown in the neocortex and hippocampus of rats with 3-day old TBI following 3 days of electrical stimulation in the median raphé nucleus (Carballosa Gonzalez et al., 2013).

Interpretation of Main Findings

Cyclic AMP can promote axonal regeneration following SCI via elevated pPKA and pCREB, which subsequently increases expression of various trophic genes, although a non-genetic route via pPKA can also have beneficial effects on the cytoskeleton via the small G protein, Rho (Lehmann et al., 1999; Hannila and Filbin, 2008). Given the previously demonstrated NRM-evoked recovery of cAMP after T8 SCI (Carballosa-Gonzalez et al., 2014), it was unsurprising that transcriptional effects of stimulation were seen. However, it was not expected that the main increases would be in defensive and pro-inflammatory genes, for example, Ctsg, Np4, RatNP-3b, Defa5, Elane, Mpo, S100a8, S100a9, and Ppbp. The proteins S100A8 and S100A9 are frequently combined as the antimicrobial dimer calprotectin (Striz and Trebichavsky, 2004). They probably arrive in the injured spinal cord as constituents of myeloid cells, particularly neutrophils (Fleming et al., 2006; Mitchell et al., 2008), which are thought to be deleterious to recovery (McCreedy et al., 2018), although S100A8 and S100A9 can also be induced in macrophages and microglia by neuropathic states (Abe et al., 1999). While there is some indirect evidence that S100A9 promotes recovery from SCI (Roet et al., 2013), but pro-inflammatory effects are more likely the dominant initial action.

Other genes upregulated by stimulation in SCI may also have deleterious effects. Mpo has been shown by gene knockout to exacerbate secondary injury after SCI (Kubota et al., 2012). Ctsg is a pronociceptive mediator in the spinal cord (Liu et al., 2015), as is carbonic anhydrases (gene Car1) (Asiedu et al., 2010). Neutrophil elastase, which is released by activated neutrophils, is reported to be a key mediator of secondary pathogenesis in SCI (Semple et al., 2015). On the other hand, matrix metalloproteinase 13 is a matrix-degrading enzyme released from monocytes that is important in functional recovery after SCI (Shechter et al., 2011). Of the several microRNAs tagged in these studies, Mir21, which was upregulated in the contrast ±SCI(st), is of greatest interest. Mir21 is increased by exercise after SCI and improves recovery in rats, probably by regulating the pro-regenerative PTEN/mTOR pathway, since it lowers PTEN messenger RNA and raises mTOR messenger RNA (Liu et al., 2012). Few of the genes affected by stimulation after SCI (Table 2) are particularly associated with the CNS. One exception was the downregulated solute carrier family 17 member 7 (Slc17a7), also known as Vglut1, whose protein transports glutamate into synaptic vesicles (Shigeri et al., 2004); its levels in the spinal cord are reported to be decreased after peripheral nerve damage (Rotterman et al., 2014; Wang et al., 2016). Acute decreases in Vglut1 (Slc17a7) with stimulation, as observed here, could lead to lower net glutamate release, which is likely to reduce damage due to over-stimulation of neurons (Krzyzanowska et al., 2014), although recovery ultimately requires restoration of Vglut1 (Burnside et al., 2018). A second exception was Robo-1, which was increased by stimulation, whose protein product an important player in pioneer longitudinal axon guidance (Kim et al., 2011).

Another interesting group of genes unexpectedly upregulated by stimulation in injury is related to erythrocyte development and oxygen transport: Klf1, Rhd, Rhag, Hbb, Hba, Hba2, Hbb-b1. The last four produce components of hemoglobin, which is constitutively expressed in the brain and spinal cord (Biagioli et al., 2009; Richter et al., 2009; Russo et al., 2013), and may act protectively by sequestering oxygen and free radicals (Xie and Yang, 2016). Indeed, “cellular oxidant detoxification” also featured among biological processes enriched by stimulation after SCI (Table 5). The key regulator of erythropoiesis, erythropoietin (EPO), is known to promote recovery and regeneration in SCI (Gassmann et al., 2003; Matis and Birbilis, 2009). There were no significant changes in expression of the Epo gene in this study, but it is possible that the increase in erythropoietic genes was due to increased erythropoiten release (Hasselblatt et al., 2006). Both the protein and Epo messenger RNA are reported to be increased by serotonin in the mouse hippocampus (Choi and Son, 2013). Further study of how erythropoiten synthesis and release responds at different time points post-injury and during NRM stimulation appears to be warranted.

The prevalence of altered transcription in exogenous cells of the injured spinal cord requires further experimental investigation. Among various possibilities to test, arriving myeloid cells such as neutrophils could be more attracted to the spinal cord’s injury site after they have been altered in some way by NRM stimulation. Alternatively, transcription in exogenous cells, once arrived, could be influenced (perhaps toward a pro-regeneration phenotype) by the NRM stimulation or by local sequelae of that stimulation. Classic synapses are not necessary for the effects of the NRM stimulation, since both serotonin and co-released neuropeptides participate in non-synaptic volume transmission in the spinal cord (Ridet et al., 1993). A growing body of evidence implicates serotonin in both neural and humoral aspects of immune control. Serotonin and its receptors are present on numerous cells of the adaptive and innate immune systems and have been implicated in macrophage activation (Jiang et al., 2009; Krabbe et al., 2012; Baganz and Blakely, 2013).

Conclusion

The present study adds important new details to the general concept that the raphé system senses and alleviates stress or injury. The effects of NRM stimulation on gene transcription were very different in injured and intact animals, mirroring our prior study of the effects of NRM stimulation on levels of cAMP, PKA, and CREB (Carballosa-Gonzalez et al., 2014). The earliest changes in gene transcription brought about by NRM stimulation were seen to be confined mainly to defensive and erythropoietic functions. Whereas a week or more of NRM stimulation facilitates repair (Hentall and Burns, 2009; Hentall and Gonzalez, 2012; Madsen et al., 2017; Vitores et al., 2018), the shorter 2-h period used here intensified transcription of predominantly inflammatory genes. Although the products of these genes are often assumed to be harmful, immune responses cannot necessarily be simply parsed into harmful and beneficial processes (Doty et al., 2015; Li et al., 2018). Conceivably, the inflammatory effects of substances such as S100A8 and S100A9, however harmful they may be initially, are necessary preliminaries with net benefits to the overall repair process, and lead later to neural-specific effects. Consistent with these ideas, three of the four top perturbagen classes that resembled the effect of 2 h of stimulation (Table 6) are known to produce improvements in SCI: PPARs agonists (Park et al., 2007), angiotensin receptor antagonists (Guler et al., 2010) and glucocorticoid receptor agonists, although this last group’s benefits in SCI are more controversial (Akhtar et al., 2009). Furthermore, all three classes inactivate inflammatory pathways, such as those controlled by NF-κB (Bekhbat et al., 2017; Labandeira-Garcia et al., 2017; Villapol, 2018). The present study thus tends to confirm the therapeutic promise of these drug classes. The study also provides further support for interim DBS of the NRM or, more practically from the surgical standpoint, of its midbrain periaqueductal gray input, as a therapy for SCI that may ultimately be superior to single or combined drugs in its simplicity and safety of application and the range of signs and symptoms improved (Chari et al., 2017).

Data Availability Statement

The datasets generated for this study can be found in the National Center for Biotechnology Information (NCBI), GEO Accession Number GSE133093.

Ethics Statement

The animal study was reviewed and approved by the Institutional Animal Care and Use Committee of the University of Miami, Miller School of Medicine.

Author Contributions

IH designed the experiments and wrote the manuscript. WJ managed the RNA sequencing and analysis. MC-G and AV performed animal surgery, stimulation, and tissue extraction.

Funding

This work was funded in part by Robert J. Kleberg, Jr. and Helen C. Kleberg Foundation.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Dr. Anthony J. Griswold for valuable help with the high-throughput sequencing and for advice on its analysis.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncel.2019.00510/full#supplementary-material

Abbreviations

cAMP, cyclic adenosine monophosphate; DAVID, Database for Annotation, Visualization and Integrated Discovery; DBS, deep brain stimulation; DGE, differential gene expression; EPO, erythropoietin; FC, fold change; MS, multiple sclerosis; NCBI, National Center for Biotechnology Information; NRM, nucleus raphé magnus; pCREB, phosphorylated cAMP response binding element; pim, pimozide; PPAR, peroxisome proliferator-activated receptor; pPKA, phosphorylated protein kinase A; SCI, spinal cord injury; st, stimulation; TBI, traumatic brain injury

Footnotes

1. ^ http://www.bioinformatics.babraham.ac.uk/projects/fastqc/
2. ^ www.bioinformatics.babraham.ac.uk/projects/trim_galore
3. ^ https://david.ncifcrf.gov/tools.jsp

References

Abe, M., Umehara, F., Kubota, R., Moritoyo, T., Izumo, S., and Osame, M. (1999). Activation of macrophages/microglia with the calcium-binding proteins MRP14 and MRP8 is related to the lesional activities in the spinal cord of HTLV-I associated myelopathy. J. Neurol. 246, 358–364. doi: 10.1007/s004150050363

PubMed Abstract | CrossRef Full Text | Google Scholar

Agnati, L. F., Bjelke, B., and Fuxe, K. (1995). Volume versus wiring transmission in the brain: a new theoretical frame for neuropsychopharmacology. Med. Res. Rev. 15, 33–45. doi: 10.1002/med.2610150104

PubMed Abstract | CrossRef Full Text | Google Scholar

Akhtar, A. Z., Pippin, J. J., and Sandusky, C. B. (2009). Animal studies in spinal cord injury: a systematic review of methylprednisolone. Altern. Lab. Anim. 37, 43–62. doi: 10.1177/026119290903700108

PubMed Abstract | CrossRef Full Text | Google Scholar

Asiedu, M., Ossipov, M. H., Kaila, K., and Price, T. J. (2010). Acetazolamide and midazolam act synergistically to inhibit neuropathic pain. Pain 148, 302–308. doi: 10.1016/j.pain.2009.11.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Azmitia, E. C. (2007). Serotonin and brain: evolution, neuroplasticity, and homeostasis. Int. Rev. Neurobiol. 77, 31–56. doi: 10.1016/S0074-7742(06)77002-7

CrossRef Full Text | Google Scholar

Baganz, N. L., and Blakely, R. D. (2013). A dialogue between the immune system and brain, spoken in the language of serotonin. ACS Chem. Neurosci. 4, 48–63. doi: 10.1021/cn300186b

PubMed Abstract | CrossRef Full Text | Google Scholar

Bareyre, F. M., and Schwab, M. E. (2003). Inflammation, degeneration and regeneration in the injured spinal cord: insights from DNA microarrays. Trends Neurosci. 26, 555–563. doi: 10.1016/j.tins.2003.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Batty, N. J., Fenrich, K. K., and Fouad, K. (2017). The role of cAMP and its downstream targets in neurite growth in the adult nervous system. Neurosci. Lett. 652, 56–63. doi: 10.1016/j.neulet.2016.12.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Bekhbat, M., Rowson, S. A., and Neigh, G. N. (2017). Checks and balances: the glucocorticoid receptor and NFkB in good times and bad. Front. Neuroendocrinol. 46:15–31. doi: 10.1016/j.yfrne.2017.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Biagioli, M., Pinto, M., Cesselli, D., Zaninello, M., Lazarevic, D., Roncaglia, P., et al. (2009). Unexpected expression of alpha- and beta-globin in mesencephalic dopaminergic neurons and glial cells. Proc. Natl. Acad. Sci. U.S.A. 106, 15454–15459. doi: 10.1073/pnas.0813216106

PubMed Abstract | CrossRef Full Text | Google Scholar

Borroto-Escuela, D. O., Agnati, L. F., Bechter, K., Jansson, A., Tarakanov, A. O., and Fuxe, K. (2015). The role of transmitter diffusion and flow versus extracellular vesicles in volume transmission in the brain neural-glial networks. Philos. Trans. R. Soc. Lond. B Biol. Sci. 370:20140183. doi: 10.1098/rstb.2014.0183

PubMed Abstract | CrossRef Full Text | Google Scholar

Burnside, E. R., De Winter, F., Didangelos, A., James, N. D., Andreica, E. C., Layard-Horsfall, H., et al. (2018). Immune-evasive gene switch enables regulated delivery of chondroitinase after spinal cord injury. Brain 141, 2362–2381. doi: 10.1093/brain/awy158

PubMed Abstract | CrossRef Full Text | Google Scholar

Campbell, A. M., Zagon, I. S., and McLaughlin, P. J. (2012). Opioid growth factor arrests the progression of clinical disease and spinal cord pathology in established experimental autoimmune encephalomyelitis. Brain Res. 1472, 138–148. doi: 10.1016/j.brainres.2012.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Carballosa Gonzalez, M. M., Blaya, M. O., Alonso, O. F., Bramlett, H. M., and Hentall, I. D. (2013). Midbrain raphe stimulation improves behavioral and anatomical recovery from fluid-percussion brain injury. J. Neurotrauma. 30, 119–130. doi: 10.1089/neu.2012.2499

PubMed Abstract | CrossRef Full Text | Google Scholar

Carballosa-Gonzalez, M. M., Vitores, A., and Hentall, I. D. (2014). Hindbrain raphe stimulation boosts cyclic adenosine monophosphate and signaling proteins in the injured spinal cord. Brain Res. 1543, 165–172. doi: 10.1016/j.brainres.2013.11.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Chamankhah, M., Eftekharpour, E., Karimi-Abdolrezaee, S., Boutros, P. C., San-Marina, S., and Fehlings, M. G. (2013). Genome-wide gene expression profiling of stress response in a spinal cord clip compression injury model. BMC Genomics 14:583. doi: 10.1186/1471-2164-14-583

PubMed Abstract | CrossRef Full Text | Google Scholar

Chari, A., Hentall, I. D., Papadopoulos, M. C., and Pereira, E. A. (2017). Surgical neurostimulation for spinal cord injury. Brain Sci. 7:E18. doi: 10.3390/brainsci7020018

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, M., and Son, H. (2013). Effects of serotonin on erythropoietin expression in mouse hippocampus. Exp. Neurobiol. 22, 45–50. doi: 10.5607/en.2013.22.1.45

PubMed Abstract | CrossRef Full Text | Google Scholar

Daimon, C. M., Chirdon, P., Maudsley, S., and Martin, B. (2013). The role of Thyrotropin Releasing Hormone in aging and neurodegenerative diseases. Am. J. Alzheimers Dis. 1. doi: 10.7726/ajad.2013.1003

PubMed Abstract | CrossRef Full Text | Google Scholar

David, S., and Kroner, A. (2011). Repertoire of microglial and macrophage responses after spinal cord injury. Nat. Rev. Neurosci. 12, 388–399. doi: 10.1038/nrn3053

PubMed Abstract | CrossRef Full Text | Google Scholar

de las Casas-Engel, M., Dominguez-Soto, A., Sierra-Filardi, E., Bragado, R., Nieto, C., Puig-Kroger, A., et al. (2013). Serotonin skews human macrophage polarization through HTR2B and HTR7. J. Immunol. 190, 2301–2310. doi: 10.4049/jimmunol.1201133

PubMed Abstract | CrossRef Full Text | Google Scholar

Di Pilato, P., Niso, M., Adriani, W., Romano, E., Travaglini, D., Berardi, F., et al. (2014). Selective agonists for serotonin 7 (5-HT7) receptor and their applications in preclinical models: an overview. Rev. Neurosci. 25, 401–415. doi: 10.1515/revneuro-2014-0009

PubMed Abstract | CrossRef Full Text | Google Scholar

Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., et al. (2013). STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21. doi: 10.1093/bioinformatics/bts635

PubMed Abstract | CrossRef Full Text | Google Scholar

Doty, K. R., Guillot-Sestier, M. V., and Town, T. (2015). The role of the immune system in neurodegenerative disorders: adaptive or maladaptive? Brain Res. 1617, 155–173. doi: 10.1016/j.brainres.2014.09.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Dulin, J. N., Antunes-Martins, A., Chandran, V., Costigan, M., Lerch, J. K., Willis, D. E., et al. (2015). Transcriptomic approaches to neural repair. J. Neurosci. 35, 13860–13867. doi: 10.1523/JNEUROSCI.2599-15.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Fields, D. P., and Mitchell, G. S. (2015). Spinal metaplasticity in respiratory motor control. Front. Neural. Circuits 9:2. doi: 10.3389/fncir.2015.00002

PubMed Abstract | CrossRef Full Text | Google Scholar

Fleming, J. C., Norenberg, M. D., Ramsay, D. A., Dekaban, G. A., Marcillo, A. E., Saenz, A. D., et al. (2006). The cellular inflammatory response in human spinal cords after injury. Brain 129(Pt 12), 3249–3269. doi: 10.1093/brain/awl296

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedel, C. C., Dolken, L., Ruzsics, Z., Koszinowski, U. H., and Zimmer, R. (2009). Conserved principles of mammalian transcriptional regulation revealed by RNA half-life. Nucleic Acids Res. 37:e115. doi: 10.1093/nar/gkp542

PubMed Abstract | CrossRef Full Text | Google Scholar

Gassmann, M., Heinicke, K., Soliz, J., and Ogunshola, O. O. (2003). Non-erythroid functions of erythropoietin. Adv. Exp. Med. Biol. 543, 323–330. doi: 10.1007/978-1-4419-8997-0_22

PubMed Abstract | CrossRef Full Text | Google Scholar

Gautier, A., Geny, D., Bourgoin, S., Bernard, J. F., and Hamon, M. (2017). Differential innervation of superficial versus deep laminae of the dorsal horn by bulbo-spinal serotonergic pathways in the rat. IBRO Rep. 2, 72–80. doi: 10.1016/j.ibror.2017.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Guler, A., Sahin, M. A., Ucak, A., Onan, B., Inan, K., Oztas, E., et al. (2010). Protective effects of angiotensin II type-1 receptor blockade with olmesartan on spinal cord ischemia-reperfusion injury: an experimental study on rats. Ann. Vasc. Surg. 24, 801–808. doi: 10.1016/j.avsg.2010.03.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Hannila, S. S., and Filbin, M. T. (2008). The role of cyclic AMP signaling in promoting axonal regeneration after spinal cord injury. Exp. Neurol. 209, 321–332. doi: 10.1016/j.expneurol.2007.06.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Hardcastle, T. J., and Kelly, K. A. (2010). baySeq: empirical Bayesian methods for identifying differential expression in sequence count data. BMC Bioinformatics 11:422. doi: 10.1186/1471-2105-11-422

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasselblatt, M., Ehrenreich, H., and Siren, A. L. (2006). The brain erythropoietin system and its potential for therapeutic exploitation in brain disease. J. Neurosurg. Anesthesiol. 18, 132–138. doi: 10.1097/00008506-200604000-00007

PubMed Abstract | CrossRef Full Text | Google Scholar

Hentall, I. D., and Burns, S. B. (2009). Restorative effects of stimulating medullary raphe after spinal cord injury. J. Rehabil. Res. Dev. 46, 109–122.

PubMed Abstract | Google Scholar

Hentall, I. D., and Gonzalez, M. M. (2012). Promotion of recovery from thoracic spinal cord contusion in rats by stimulation of medullary raphe or its midbrain input. Neurorehabil. Neural. Repair. 26, 374–384. doi: 10.1177/1545968311425178

PubMed Abstract | CrossRef Full Text | Google Scholar

Hobson, S. A., Bacon, A., Elliot-Hunt, C. R., Holmes, F. E., Kerr, N. C., Pope, R., et al. (2008). Galanin acts as a trophic factor to the central and peripheral nervous systems. Cell. Mol. Life Sci. 65, 1806–1812. doi: 10.1007/s00018-008-8154-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Hornung, J. P. (2003). The human raphe nuclei and the serotonergic system. J. Chem. Neuroanat. 26, 331–343. doi: 10.1016/j.jchemneu.2003.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Jermakowicz, W. J., Sloley, S. S., Dan, L., Vitores, A., Carballosa-Gautam, M. M., and Hentall, I. D. (2019). Cellular changes in injured rat spinal cord following electrical brainstem stimulation. Brain Sci. 9:E124. doi: 10.3390/brainsci9060124

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, M. H., Chung, E., Chi, G. F., Ahn, W., Lim, J. E., Hong, H. S., et al. (2012). Substance P induces M2-type macrophages after spinal cord injury. Neuroreport 23, 786–792. doi: 10.1097/WNR.0b013e3283572206

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, X., Wang, J., Luo, T., and Li, Q. (2009). Impaired hypothalamic-pituitary-adrenal axis and its feedback regulation in serotonin transporter knockout mice. Psychoneuroendocrinology 34, 317–331. doi: 10.1016/j.psyneuen.2008.09.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Kearney, P. A., Ridella, S. A., Viano, D. C., and Anderson, T. E. (1988). Interaction of contact velocity and cord compression in determining the severity of spinal cord injury. J. Neurotrauma. 5, 187–208. doi: 10.1089/neu.1988.5.187

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, K. T., Kim, H. J., Cho, D. C., Bae, J. S., and Park, S. W. (2015). Substance P stimulates proliferation of spinal neural stem cells in spinal cord injury via the mitogen-activated protein kinase signaling pathway. Spine J. 15, 2055–2065. doi: 10.1016/j.spinee.2015.04.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, M., Roesener, A. P., Mendonca, P. R., and Mastick, G. S. (2011). Robo1 and Robo2 have distinct roles in pioneer longitudinal axon guidance. Dev. Biol. 358, 181–188. doi: 10.1016/j.ydbio.2011.07.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Krabbe, G., Matyash, V., Pannasch, U., Mamer, L., Boddeke, H. W., and Kettenmann, H. (2012). Activation of serotonin receptors promotes microglial injury-induced motility but attenuates phagocytic activity. Brain Behav. Immun. 26, 419–428. doi: 10.1016/j.bbi.2011.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Krzyzanowska, W., Pomierny, B., Filip, M., and Pera, J. (2014). Glutamate transporters in brain ischemia: to modulate or not? Acta Pharmacol. Sin. 35, 444–462. doi: 10.1038/aps.2014.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kubota, K., Saiwai, H., Kumamaru, H., Maeda, T., Ohkawa, Y., Aratani, Y., et al. (2012). Myeloperoxidase exacerbates secondary injury by generating highly reactive oxygen species and mediating neutrophil recruitment in experimental spinal cord injury. Spine 37, 1363–1369. doi: 10.1097/BRS.0b013e31824b9e77

PubMed Abstract | CrossRef Full Text | Google Scholar

Kvachnina, E., Liu, G., Dityatev, A., Renner, U., Dumuis, A., Richter, D. W., et al. (2005). 5-HT7 receptor is coupled to G alpha subunits of heterotrimeric G12-protein to regulate gene transcription and neuronal morphology. J. Neurosci. 25, 7821–7830. doi: 10.1523/JNEUROSCI.1790-05.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Kwiat, G. C., and Basbaum, A. I. (1992). The origin of brainstem noradrenergic and serotonergic projections to the spinal cord dorsal horn in the rat. Somatosens Mot. Res. 9, 157–173. doi: 10.3109/08990229209144768

PubMed Abstract | CrossRef Full Text | Google Scholar

Labandeira-Garcia, J. L., Rodriguez-Perez, A. I., Garrido-Gil, P., Rodriguez-Pallares, J., Lanciego, J. L., and Guerra, M. J. (2017). Brain renin-angiotensin system and microglial polarization: implications for aging and neurodegeneration. Front. Aging Neurosci. 9:129. doi: 10.3389/fnagi.2017.00129

PubMed Abstract | CrossRef Full Text | Google Scholar

Lehmann, M., Fournier, A., Selles-Navarro, I., Dergham, P., Sebok, A., Leclerc, N., et al. (1999). Inactivation of Rho signaling pathway promotes CNS axon regeneration. J. Neurosci. 19, 7537–7547. doi: 10.1523/jneurosci.19-17-07537.1999

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, T., Peng, M., Yang, Z., Zhou, X., Deng, Y., Jiang, C., et al. (2018). 3D-printed IFN-gamma-loading calcium silicate-beta-tricalcium phosphate scaffold sequentially activates M1 and M2 polarization of macrophages to promote vascularization of tissue engineering bone. Acta Biomater. 71, 96–107. doi: 10.1016/j.actbio.2018.03.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, G., Detloff, M. R., Miller, K. N., Santi, L., and Houle, J. D. (2012). Exercise modulates microRNAs that affect the PTEN/mTOR pathway in rats after spinal cord injury. Exp. Neurol. 233, 447–456. doi: 10.1016/j.expneurol.2011.11.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X., Tian, Y., Meng, Z., Chen, Y., Ho, I. H., Choy, K. W., et al. (2015). Up-regulation of cathepsin G in the development of chronic postsurgical pain: an experimental and clinical genetic study. Anesthesiology 123, 838–850. doi: 10.1097/ALN.0000000000000828

PubMed Abstract | CrossRef Full Text | Google Scholar

Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Madsen, P. M., Sloley, S. S., Vitores, A. A., Carballosa-Gautam, M. M., Brambilla, R., and Hentall, I. D. (2017). Prolonged stimulation of a brainstem raphe region attenuates experimental autoimmune encephalomyelitis. Neuroscience 346, 395–402. doi: 10.1016/j.neuroscience.2017.01.037

PubMed Abstract | CrossRef Full Text | Google Scholar

Mao, Y., Mathews, K., and Gorrie, C. A. (2016). Temporal response of endogenous neural progenitor cells following injury to the adult rat spinal cord. Front. Cell Neurosci. 10:58. doi: 10.3389/fncel.2016.00058

PubMed Abstract | CrossRef Full Text | Google Scholar

Matis, G. K., and Birbilis, T. A. (2009). Erythropoietin in spinal cord injury. Eur. Spine J. 18, 314–323. doi: 10.1007/s00586-008-0829-0

PubMed Abstract | CrossRef Full Text | Google Scholar

McCreedy, D. A., Lee, S., Sontag, C. J., Weinstein, P., Olivas, A. D., Martinez, A. F., et al. (2018). Early targeting of L-selectin on leukocytes promotes recovery after spinal cord injury, implicating novel mechanisms of pathogenesis. eNeuro 5:ENEURO.0101-18.2018. doi: 10.1523/ENEURO.0101-18.2018

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitchell, K., Yang, H. Y., Tessier, P. A., Muhly, W. T., Swaim, W. D., Szalayova, I., et al. (2008). Localization of S100A8 and S100A9 expressing neutrophils to spinal cord during peripheral tissue inflammation. Pain 134, 216–231. doi: 10.1016/j.pain.2007.10.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Molliver, M. E. (1987). Serotonergic neuronal systems: what their anatomic organization tells us about function. J. Clin. Psychopharmacol. 7(6 Suppl.), 3S–23S.

PubMed Abstract | Google Scholar

Narita, M., Kuzumaki, N., Miyatake, M., Sato, F., Wachi, H., Seyama, Y., et al. (2006). Role of delta-opioid receptor function in neurogenesis and neuroprotection. J. Neurochem. 97, 1494–1505. doi: 10.1111/j.1471-4159.2006.03849.x

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Shea, T. M., Burda, J. E., and Sofroniew, M. V. (2017). Cell biology of spinal cord injury and repair. J. Clin. Invest. 127, 3259–3270. doi: 10.1172/JCI90608

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, S. W., Yi, J. H., Miranpuri, G., Satriotomo, I., Bowen, K., Resnick, D. K., et al. (2007). Thiazolidinedione class of peroxisome proliferator-activated receptor gamma agonists prevents neuronal damage, motor dysfunction, myelin loss, neuropathic pain, and inflammation after spinal cord injury in adult rats. J. Pharmacol. Exp. Ther. 320, 1002–1012. doi: 10.1124/jpet.106.113472

PubMed Abstract | CrossRef Full Text | Google Scholar

Rabani, M., Levin, J. Z., Fan, L., Adiconis, X., Raychowdhury, R., Garber, M., et al. (2011). Metabolic labeling of RNA uncovers principles of RNA production and degradation dynamics in mammalian cells. Nat. Biotechnol. 29, 436–442. doi: 10.1038/nbt.1861

PubMed Abstract | CrossRef Full Text | Google Scholar

Rabani, M., Raychowdhury, R., Jovanovic, M., Rooney, M., Stumpo, D. J., Pauli, A., et al. (2014). High-resolution sequencing and modeling identifies distinct dynamic RNA regulatory strategies. Cell 159, 1698–1710. doi: 10.1016/j.cell.2014.11.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Richter, F., Meurers, B. H., Zhu, C., Medvedeva, V. P., and Chesselet, M. F. (2009). Neurons express hemoglobin alpha- and beta-chains in rat and human brains. J. Comp. Neurol. 515, 538–547. doi: 10.1002/cne.22062

PubMed Abstract | CrossRef Full Text | Google Scholar

Ridet, J. L., Rajaofetra, N., Teilhac, J. R., Geffard, M., and Privat, A. (1993). Evidence for nonsynaptic serotonergic and noradrenergic innervation of the rat dorsal horn and possible involvement of neuron-glia interactions. Neuroscience 52, 143–157. doi: 10.1016/0306-4522(93)90189-m

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616

PubMed Abstract | CrossRef Full Text | Google Scholar

Roet, K. C., Franssen, E. H., de Bree, F. M., Essing, A. H., Zijlstra, S. J., Fagoe, N. D., et al. (2013). A multilevel screening strategy defines a molecular fingerprint of proregenerative olfactory ensheathing cells and identifies SCARB2, a protein that improves regenerative sprouting of injured sensory spinal axons. J. Neurosci. 33, 11116–11135. doi: 10.1523/JNEUROSCI.1002-13.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Roth, B. L., Craigo, S. C., Choudhary, M. S., Uluer, A., Monsma, F. J. Jr., Shen, Y., et al. (1994). Binding of typical and atypical antipsychotic agents to 5-hydroxytryptamine-6 and 5-hydroxytryptamine-7 receptors. J. Pharmacol. Exp. Ther. 268, 1403–1410.

PubMed Abstract | Google Scholar

Rotterman, T. M., Nardelli, P., Cope, T. C., and Alvarez, F. J. (2014). Normal distribution of VGLUT1 synapses on spinal motoneuron dendrites and their reorganization after nerve injury. J. Neurosci. 34, 3475–3492. doi: 10.1523/JNEUROSCI.4768-13.2014

PubMed Abstract | CrossRef Full Text | Google Scholar

Roy, H. A., and Aziz, T. Z. (2014). Deep brain stimulation and multiple sclerosis: therapeutic applications. Mult. Scler. Relat. Disord. 3, 431–439. doi: 10.1016/j.msard.2014.02.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Russo, R., Zucchelli, S., Codrich, M., Marcuzzi, F., Verde, C., and Gustincich, S. (2013). Hemoglobin is present as a canonical alpha2beta2 tetramer in dopaminergic neurons. Biochim. Biophys. Acta 1834, 1939–1943. doi: 10.1016/j.bbapap.2013.05.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Schurch, N. J., Schofield, P., Gierlinski, M., Cole, C., Sherstnev, A., Singh, V., et al. (2016). How many biological replicates are needed in an RNA-seq experiment and which differential expression tool should you use? RNA 22, 839–851. doi: 10.1261/rna.053959.115

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwanhausser, B., Busse, D., Li, N., Dittmar, G., Schuchhardt, J., Wolf, J., et al. (2011). Global quantification of mammalian gene expression control. Nature 473, 337–342. doi: 10.1038/nature10098

PubMed Abstract | CrossRef Full Text | Google Scholar

Semple, B. D., Trivedi, A., Gimlin, K., and Noble-Haeusslein, L. J. (2015). Neutrophil elastase mediates acute pathogenesis and is a determinant of long-term behavioral recovery after traumatic injury to the immature brain. Neurobiol. Dis. 74, 263–280. doi: 10.1016/j.nbd.2014.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Shechter, R., Raposo, C., London, A., Sagi, I., and Schwartz, M. (2011). The glial scar-monocyte interplay: a pivotal resolution phase in spinal cord repair. PLoS One 6:e27969. doi: 10.1371/journal.pone.0027969

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, L. L., Zhang, N., Xie, X. M., Chen, Y. J., Wang, R., Shen, L., et al. (2017). Transcriptome profile of rat genes in injured spinal cord at different stages by RNA-sequencing. BMC Genomics 18:173. doi: 10.1186/s12864-017-3532-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Shigeri, Y., Seal, R. P., and Shimamoto, K. (2004). Molecular pharmacology of glutamate transporters, EAATs and VGLUTs. Brain Res. Brain Res. Rev. 45, 250–265. doi: 10.1016/j.brainresrev.2004.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Shin, S. S., Dixon, C. E., Okonkwo, D. O., and Richardson, R. M. (2014). Neurostimulation for traumatic brain injury. J. Neurosurg. 121, 1219–1231. doi: 10.3171/2014.7.JNS131826

PubMed Abstract | CrossRef Full Text | Google Scholar

Striz, I., and Trebichavsky, I. (2004). Calprotectin - a pleiotropic molecule in acute and chronic inflammation. Physiol. Res. 53, 245–253.

PubMed Abstract | Google Scholar

Subramanian, A., Narayan, R., Corsello, S. M., Peck, D. D., Natoli, T. E., Lu, X., et al. (2017). A next generation connectivity map: L1000 platform and the first 1,000,000 profiles. Cell 171, 1437–1452.e17. doi: 10.1016/j.cell.2017.10.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Trakhtenberg, E. F., and Goldberg, J. L. (2012). The role of serotonin in axon and dendrite growth. Int. Rev. Neurobiol. 106, 105–126. doi: 10.1016/B978-0-12-407178-0.00005-3

PubMed Abstract | CrossRef Full Text | Google Scholar

van Dijken, H., Dijk, J., Voom, P., and Holstege, J. C. (1996). Localization of dopamine D2 receptor in rat spinal cord identified with immunocytochemistry and in situ hybridization. Eur. J. Neurosci. 8, 621–628. doi: 10.1111/j.1460-9568.1996.tb01247.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Villapol, S. (2018). Roles of peroxisome proliferator-activated receptor gamma on brain and peripheral inflammation. Cell. Mol. Neurobiol. 38, 121–132. doi: 10.1007/s10571-017-0554-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Vitores, A. A., Sloley, S. S., Martinez, C., Carballosa-Gautam, M. M., and Hentall, I. D. (2018). Some autonomic deficits of acute or chronic cervical spinal contusion reversed by interim brainstem stimulation. J. Neurotrauma. 35, 560–572. doi: 10.1089/neu.2017.5123

PubMed Abstract | CrossRef Full Text | Google Scholar

Volpicelli, F., Speranza, L., di Porzio, U., Crispino, M., and Perrone-Capano, C. (2014). The serotonin receptor 7 and the structural plasticity of brain circuits. Front. Behav. Neurosci. 8:318. doi: 10.3389/fnbeh.2014.00318

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, H. S., Yu, G., Wang, Z. T., Yi, S. P., Su, R. B., and Gong, Z. H. (2016). Changes in VGLUT1 and VGLUT2 expression in rat dorsal root ganglia and spinal cord following spared nerve injury. Neurochem. Int. 99, 9–15. doi: 10.1016/j.neuint.2016.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, L. K., and Yang, S. H. (2016). Brain globins in physiology and pathology. Med. Gas Res. 6, 154–163. doi: 10.4103/2045-9912.191361

PubMed Abstract | CrossRef Full Text | Google Scholar

Zerbino, D. R., Achuthan, P., Akanni, W., Amode, M. R., Barrell, D., Bhai, J., et al. (2018). Ensembl 2018. Nucleic Acids Res. 46, D754–D761. doi: 10.1093/nar/gkx1098

PubMed Abstract | CrossRef Full Text | Google Scholar