Peripheral nerve injury is accompanied by chronic transcriptome-wide changes in the mouse prefrontal cortex
Molecular Pain volume 9, Article number: 21 (2013)
Peripheral nerve injury can have long-term consequences including pain-related manifestations, such as hypersensitivity to cutaneous stimuli, as well as affective and cognitive disturbances, suggesting the involvement of supraspinal mechanisms. Changes in brain structure and cortical function associated with many chronic pain conditions have been reported in the prefrontal cortex (PFC). The PFC is implicated in pain-related co-morbidities such as depression, anxiety and impaired emotional decision-making ability. We recently reported that this region is subject to significant epigenetic reprogramming following peripheral nerve injury, and normalization of pain-related structural, functional and epigenetic abnormalities in the PFC are all associated with effective pain reduction.
In this study, we used the Spared Nerve Injury (SNI) model of neuropathic pain to test the hypothesis that peripheral nerve injury triggers persistent long-lasting changes in gene expression in the PFC, which alter functional gene networks, thus providing a possible explanation for chronic pain associated behaviors.
SNI or sham surgery where performed in male CD1 mice at three months of age. Six months after injury, we performed transcriptome-wide sequencing (RNAseq), which revealed 1147 differentially regulated transcripts in the PFC in nerve-injured vs. control mice. Changes in gene expression occurred across a number of functional gene clusters encoding cardinal biological processes as revealed by Ingenuity Pathway Analysis. Significantly altered biological processes included neurological disease, skeletal muscular disorders, behavior, and psychological disorders. Several of the changes detected by RNAseq were validated by RT-QPCR and included transcripts with known roles in chronic pain and/or neuronal plasticity including the NMDA receptor (glutamate receptor, ionotropic, NMDA; grin1), neurite outgrowth ( roundabout 3; robo3), gliosis ( glial fibrillary acidic protein; gfap), vesicular release (synaptotagmin 2; syt2), and neuronal excitability ( voltage-gated sodium channel, type I; scn1a).
This study used an unbiased approach to document long-term alterations in gene expression in the brain following peripheral nerve injury. We propose that these changes are maintained as a memory of an insult that is temporally and spatially distant from the initial injury.
Peripheral nerve injury can result in a multitude of changes within an organism, including motor dysfunction, pain and associated cognitive and emotional comorbidities. While acute pain related to an injury is protective and normally resolves, chronic pain can be detrimental to the overall wellbeing and functioning of the individual. Indeed, occasionally the “memory” of injury, in the form of chronic pain, persists long after the initial recovery phase and becomes difficult to reverse. This is due, in part, to changes in anatomy and function that take place in the peripheral as well as the central nervous system. These changes can occur at many levels: individual molecules, synapses, cellular function, and network activity . It is therefore not surprising that injury is often accompanied by local or systemic alterations in gene expression. To date, several reports have identified strong links between injury and transcriptional changes in the peripheral nervous system , blood  and in the brain . However, the full profile of transcriptional changes that accompany chronic pain in response to peripheral injury is unknown.
The fact that peripheral injury results in chronic behavioral changes suggests that transient exposure to injury is chronically embedded in the transcription programming within the central nervous system, resulting in altered phenotypic behaviors. In the current study, we tested this hypothesis by focusing on the transcriptional changes that occur in the prefrontal cortex (PFC) of chronically neuropathic mice six months after the induction of neuropathy. This area was chosen based on evidence indicating its involvement in pain modulation [5, 6] and epigenetic differences that accompany peripheral nerve injury . In humans with back pain, reversible pathological changes in both cortical thickness and functional activation have been shown in the PFC . In animal models of neuropathic pain, the PFC undergoes synaptic re-organization as early as 8 days following nerve injury , and shows reductions in grey matter 5 months post-injury . Since the PFC has also been implicated in depression and anxiety , common co-morbidities of chronic pain, transcriptional changes in brain regions such as the PFC can also provide an explanation for these co-morbidities.
While brain-specific transcription associated with neuropathic pain and injury has been previously investigated, few studies have looked at the effects at time points longer than ~ 1 month post-injury, whereas individual patients often suffer for many years following the initial injury. This distinction is critical as the development of co-morbid conditions, such as anxiety-like behavior in the rat, takes multiple months to develop . Thus, early time points might not fully incorporate the impact of long term chronic pain on CNS plasticity. Furthermore, previous studies have been limited by microarray-based technologies that are inherently biased by probe density/design and limited to coding mRNA transcripts . We therefore investigated the long-term (6 months post-injury) transcriptional changes induced by peripheral injury in the prefrontal cortex using whole transcriptome sequencing (RNAseq) in a mouse model of chronic pain (Spared Nerve Injury, SNI). Furthermore, we employed bioinformatic tools to identify functional gene clusters that were altered in the prefrontal cortex.
The aim of this work is to provide a comprehensive and unbiased look at the molecular correlates of peripheral nerve injury and catalog long-term transcriptional changes that may play a role in the pathologies associated with injury and pain. We believe that our findings will help shed light on the “signature” of painful neuropathy in the brain at both the molecular and network levels.
Peripheral injury is accompanied by behavioral signs of neuropathic pain six months post-injury
The persistence of nerve injury-induced hypersensitivity to mechanical and cold stimuli and injury-related motor impairment were confirmed six months following SNI (data not shown; mechanical thresholds (grams) = 0.20±0.05 in SNI vs 0.82±0.07 in controls, p<0.0001; acetone-evoked behaviors (seconds) = 2.9±0.4 vs. 0.4±0.05 in controls, p>0.0001; motor impairment (latency to fall from accelerating rotarod) = 76±11 in SNI vs. 225±21 in controls, p>0.0001; n=10/group).
Peripheral injury is accompanied transcriptomic changes in the prefrontal cortex six months post-injury
Six months following SNI, significant changes in the expression levels of 1147 different transcripts were identified in the prefrontal cortex (Figure 1, Additional file 1: Table S1). Considering the genomic base pair representation of exonic (3.2%), intronic (35.5%) and intergenic (61.3%) elements in mouse (Figure 1A), the PFC SNI-associated transcriptome was equally partitioned between coding RNA accounting for 40% (exonic) and 60% non-coding (intergenic+intronic) base pairs of the transcriptome (Figure 1B). Within differentially expressed transcripts in SNI vs. sham control animals, the largest transcriptional changes were observed in protein-coding exons (63.4%) with the rest made up of non-coding RNA. Within non-coding RNAs, the largest observed changes were in non-translated transcripts (retained introns and processed transcripts) followed by classes of microRNAs, transcribed processed pseudogenes, lincRNAs and nonsense mediated decay (Figure 1C).
The transcripts with the largest log2fold changes in RNA expression in response to SNI (summarized in Table 1) included several transcripts with previously uncharacterized roles in the brain. Within this subset, specific transcripts with a role in subcellular structure (actbl, col11a2, krt12, krt20) and development pathways (lrat, aldh3a1, crb1) show some of the greatest transcriptional changes. These transcripts showed large transcriptional differences from control ranging between 4.20 and 5.20 log2fold change for repressed genes and between 2.61 and 5.60 log2fold for induced genes.
We further classified differentially expressed transcripts in SNI animals that have brain-specific functions (Summarized in Table 2). Examples include down-regulation of neurotransmitter channel and receptor subunits gabrg1/3, clca1, slc14a1 and drd2 and the astrocyte marker gfap while sodium channel subunit scn1a, NMDA subunit grin1, and promoters of neuronal growth xlr4b, robo3, prc, and cux1 were all up-regulated in PFC from injured mice. Specific transcripts identified with asterisks in Tables 1 and 2 were further validated with RT-qPCR (labeled with asterisks) and summarized in Figures 2A-D, 3B-C, 4B-C and 5B. Within these validated genes, robo3, scn1a, grin1, xlr4b, krt20, syt2, and lbp showed marked induction following SNI while gfap and clca1 showed marked repression (unpaired 2-tailed t-test, p<0.05, n=8). Whereas the genes in Figure 2 were selected for validation due to an interesting role in the CNS, Figures 4, 5, and 6 highlight validated genes within the context of an identified functional gene cluster.
Peripheral nerve injury results in changed brain-specific transcriptional programs
Ingenuity Pathway Analysis (IPA) identified specific networks that were dysregulated six months post-SNI in the prefrontal cortex. Biological functions with a cutoff p-value < 0.05 were considered statistically significant. We identified the following biological functions specific to peripheral nerve injury: neurological disease (p-value= 2.02E-21-1.46E-02), behavior (p-value= 2.68E-12-1.31E-02), psychological disorder (p-value=7.80E-11- 1.46E-02), protein synthesis (p-value=1.70E-10 to 4.23E-03), and nervous system development and function (p-value=5.93E-10 to 1.42E-02) (Figure 6, Additional file 2: Table S2). In our screen, genetic disorders and skeletal muscular disorders emerged but due to an unclear role in CNS we did not perform further validation within these pathways.
With each biological function comprised of hundreds of functional gene clusters we then focused on transcriptional clusters previously associated with neurological function. Within this scope, peripheral injury resulted in up-regulation of pathways involved in cellular growth and proliferation (Figure 3A), molecular transport and neurological disease (Figure 4A) and neuronal development (Figure 5A). We validated the induction of representative genes included in these clusters by QRT-PCR: robo3 and krt20 in the cellular growth and proliferation pathway (Figure 3B-C), scn1a and syt2 in the molecular transport and neurological disease pathway (Figure 4B-C) and grin1 in the neuronal development pathway (Figure 5B).
Our results delineate for the first time a transcriptomic signature in the prefrontal cortex resulting from peripheral injury six months prior. Interestingly, both coding and non-coding transcripts are altered.
The coding transcripts include both genes that were previously implicated in the pathology associated with neuronal plasticity as well as genes with yet an unknown role in brain function, neuronal plasticity or in chronic pain. We further mapped the functional gene pathways whose transcription was altered, identified specific clusters involved in neuronal plasticity and validated candidate genes within these pathways. Genes known to play a role in brain structure and function that were differentially expressed and validated in response to peripheral injury were: clca1, syt2, grin1, scn1a, krt20, xlr4b, gfap, lbp and robo3. Considering the importance of the PFC to chronic pain and its associated co-morbidities, these broad, functionally-relevant changes in the transcriptome of the PFC provides a possible substrate for the long-term systemic effects of peripheral injury and may elucidate transcriptional mechanisms of supraspinal pathologies associated with chronic pain.
Coding/noncoding transcriptomic changes
While annotated coding transcripts account for about a third (n=209000) of the long-term differential transcriptome in response to SNI, the remaining differentially-expressed transcripts were non-coding unannotated RNA (intergenic+intronic) (n=435000) (Figure 1B). Within differential exonic transcripts, we further classified all annotated genes into their respective RNA classes (Figure 1C). Approximately 60% of exonic elements represent protein-coding genes and the remaining 40% demonstrate a wide array of noncoding RNA with a previously uncharacterized role in peripheral nerve injury (Additional file 1: Table S1). The role of noncoding RNA has been discussed extensively [16, 17] and it is known to be involved in a variety of functions ranging from translation , splicing  to transcriptional regulation . The fact that SNI induces non-coding (both annotated and un-annotated) transcripts in the brain that could be detected long after the initial injury suggests that they might be playing a functional role in the brain response to peripheral injury. However, the mechanisms involved remain unknown at this stage. Nevertheless, our data provide further support to the emerging idea that genome-function in the brain involves more than the commonly studied protein coding gene sequences.
Annotated SNI-associated transcriptional differences were categorized into biological functions using ingenuity pathway analysis and we identified statistically significant affected gene pathways pertinent to neurological disease, behavior and psychological disorders (Figure 6 Additional file 2: Table S2). The SNI-associated transcriptome of the PFC appears particularly relevant to several of the co-morbidities associated with peripheral nerve injury such as poor sleep, anxiety and depression  as the PFC is highly implicated in these conditions [22–24]. Considering the identified molecular pathways (Figure 6), chronic neuropathic pain may alter higher order biological functions that may mediate these behavioral disorders.
Genes involved in chronic injury and altered neuronal growth and proliferation in the PFC
The up-regulation of robo3, krt20, xlr4b in response to peripheral injury and overall pathway analysis suggests altered regulation of gene networks involved in neuronal cellular growth (Figures 2 & 3, Additional file 2: Table S2). Previous work using the same model of peripheral nerve injury reported altered neuronal growth in injured animals characterized by longer basal dendrites with more branches in cortical layer 2/3 pyramidal neurons in injured vs. control animals . The cell adhesion membrane protein ROBO3, which is up-regulated in the injured mice, has been shown to be essential for axonal guidance in drosophila  and is a regulator of cortical interneuron growth in mice  and serotonergic neuronal differentiation . While not previously considered in the context of peripheral injury, its role in regulating morphology and function of neurons is consistent with a possible contribution to cortical remodeling in response to injury (Figure 3B).
xlr4b, a non-coding RNA, has been shown to regulate chromatin remodeling as well as modulate cognitive defects in a mouse model of Turner’s Syndrome . Upstream regulators of xlr4b, cux1 and cux2 control the number and maturity of dendritic spines in cortical neurons in cortical layer 2/3 pyramidal neurons . Interestingly, in addition to xlr4b, cux1 is up-regulated (4 fold) in the SNI animals, supporting the hypothesis that this gene regulatory pathway is activated (Additional file 1: Table S1 and Figure 3A). clca1, which is down-regulated in the PFC in response to peripheral nerve injury (Figure 2A), encodes a calcium chloride channel. Calcium chloride channels are highly conserved structurally but have diverse expression patterns and physiological functions [30, 31]. CLCA1 has been implicated in neuronal death , suggesting that CLCA1 might play role in neuronal reorganization in the PFC.
Lastly, while krt20 is not directly implicated in neuronal growth, as a soft keratin and intermediary filament, its expression is ubiquitous. krt20 is up-regulated in the PFC (Figure 3C) in response to peripheral nerve injury and it is considered a member of functional pathways involved in cell cycle and proliferation (IPA analysis). Taken together, these data are consistent with alterations in regulation of cell growth and proliferation as a chronic response to peripheral nerve injury. However, direct experiments need to test whether these transcriptional changes in the PFC are causally involved in altering neuronal or associated cell growth in the PFC.
Molecular transport and neurological disease
Chronic pain changes the brain at various levels: changes in cortical grey matter [33, 34], white matter , and overall cortical excitability  have all been reported. Recent studies have shown that chronic pain is linked to changes in synaptic structure and function in the primary somatosensory cortex , hippocampus , and anterior cingulate cortex . It is therefore not surprising that several neurological disorders are co-morbid with chronic pain. Previous reports have shown chronic pain is associated with impaired decision-making , anxiety, depression and sleep disorders . In the SNI model, the PFC exhibits signs of damage following peripheral injury [9, 10]. Considering that the PFC is shown to regulate higher order emotional processing, transcriptional changes in the PFC can identify putative candidates mediating known co-morbid pathologies seen with chronic injury and/or pain.
Mutations in sodium channels have been previously implicated in pain and in the modulation of neuropathic and inflammatory pain [42, 43]. In the present study, IPA-identified pathways reveal that the gene encoding sodium channel scn1a, which is a member of the neurological disease functional gene cluster, is up-regulated in the SNI group (Figure 4C). Familial missense mutations in scn1a have been linked to epilepsy , migraines  and brain structure in aging individuals . Interestingly, these familial missense mutations can result in increased sodium channel activity; for example, the D188V mutation (linked to epilepsy) results in impaired inactivation, and the Q1489K mutation (linked to familial hemiplegic migraine) results in changes in channel gating which include accelerated recovery from inactivation and increased persistent current . Thus, the functional result of these missense mutations may be analogous to the impact of nerve injury-induced up-regulation of the same channel.
This is the first demonstration of long-term overexpression of scn1a in the prefrontal cortex following peripheral injury. Nevertheless our data is consistent with previous reports showing up-regulation of SCN3A, for example, in dorsal root ganglia following axotomy . The data is consistent with the hypothesis that scn1a over-expression (and to a lesser degree scn1b, scn2b, and scn4b (Additional file 1: Table S1)) is involved in increased cortical excitability  associated with chronic pain.
Synaptotagmin II (syt2) is up-regulated in our model in response to peripheral injury (Figure 3B) which is consistent with its role in neurological function and disease. Synatoptoagmins are Ca2+ sensors involved in vesicular trafficking and exocytosis and mediate vesicular release important for neurotransmission . Syt2 was previously reported to be involved in vesicular GABAergic transmission , and may therefore also be involved in pain-mediated alterations in GABA and glutamate receptor activation known to accompany models of peripheral injury .
LBP, a mediator of innate immunity , is up-regulated in the PFC of animals with peripheral nerve injury. Interestingly, LBP-KO mice show increased spine density and abnormal spine morphology, indicating a role for LBP role in synaptic pruning . We speculate that its role in TLR4 signaling may link it to the neuroinflammation and depression induced by treatment with lipopolysaccharide . Altered LBP expression in the PFC in response to SNI might be playing similar roles in the gene networks, possibly mediating cortical neuroinflammation in animals with peripheral nerve injury.
Nervous system development and function
One of the functional gene pathways that is highly enriched by genes that are differentially expressed six months after nerve injury is the pathway involved in nervous system development and functions (Figure 5, Additional file 2: Table S2). One member of this cluster which is up-regulated by SNI is grin1 (NMDA receptor subunit 1). NMDA receptors are linked to neuroplastic changes including long-term potentiation  and implicated in psychiatric disorders [57–59] and anxiety [60, 61].
Within the nervous system development and function cluster we found a marked decrease in the expression of gfap, a marker for astrocytes. Astrocytes are implicated in synapse maintenance, secretion of neurotrophins, and overall function in maintaining neuroplasticity. Previous work has shown that peripheral nerve injury is associated with increased levels of gfap in the spinal and medullary dorsal horns [62, 63], anterior cingulate cortex [64–66], and the periaqueductal gray, during early and intermediate time points after injury . However, there is evidence showing gfap is down-regulated in the PFC in psychological disorders [67, 68]. Furthermore, pain-related pathological changes in neuroinflammation in the spinal cord, such as microglial activation, may occur in the absence of parallel dysregulation in supraspinal structures , and glial activation following spinal cord injury is both spatially and temporally regulated. Thus, the well documented increases in astroctye activation in the spinal cord following nerve injury  may not be mirrored in the PFC and/or gfap could be up-regulated during the early and middle stages of injury but down-regulated at very chronic stages once the co-morbidities have appeared.
We demonstrate broad changes in gene expression in the mouse prefrontal cortex six months after peripheral nerve injury, illustrating a long-term impact of a peripheral injury on brain genome function. The use of RNAseq allowed for an unbiased picture of the transcriptomic changes involved in chronic neuropathy in the PFC. Furthermore, the ingenuity pathway analysis revealed functional gene networks that were significantly implicated including networks involved in brain development and function. These reported changes and their functional analysis will generate hypotheses on the molecular mechanisms that mediate chronic pain and its co-morbidities that will need to be tested in future experiments.
These wide-spread changes in gene expression in the PFC are consistent with our previous report demonstrated significant genome-wide changes in global DNA methyation that we predicted would result in the dysregulation of hundreds of genes. Thus, epigenetic mechanisms that embed the transient peripheral injury into a long-term programmatic change in gene function in the brain may be contribution to these mechanisms .
Although we demonstrate a causal relationship between peripheral injury and a transcriptome change six months later, it is unknown whether these are the same changes that occurred in the brain at the time of injury or whether (and more probably) a cascade of gene expression alterations led to the chronic profile observed in our study. Furthermore, it is unknown if similar changes occur in other spinal or supraspinal regions. Further experiments are required at multiple time point and in additional brain regions, including those not implicated in pain signaling, before the significance of these findings can be fully understood. Finally, there is currently no evidence that these patterns are the cause of chronic pain or its associated co-morbidities. Nevertheless, our study strongly supports the plausibility that long-term changes in gene expression in the CNS are involved in chronic pain and its associated behaviors, and generates hypotheses on the genes and functional gene networks that might be involved.
8–10 week-old male CD1 mice (Charles River, St-Constant, QC, Canada) were used. Animals were housed in ventilated polycarbonate cages and received water and standard laboratory rodent diet (Harlan Teklad, soy-free diet 2020X) ad libitum. All experiments were approved by the Animal Care Committee at McGill University, and conformed to the ethical guidelines of the Canadian Council on Animal Care and the guidelines of the Committee for Research and Ethical Issues of the International Association for the Study of Pain published in PAIN, 16 (1983) 109–110. All surgery was performed under isoflurane anesthesia, and all efforts were made to minimize suffering.
Induction of nerve injury
Neuropathy was induced using the spared nerve injury model. Under deep anesthesia, an incision was made on the lateral surface of the thigh through the muscle, exposing the three terminal branches of the sciatic nerve: the sural, common peroneal and tibial nerves. The common peroneal and the tibial nerves were tightly ligated with 6.0 silk (Ethicon) and sectioned distal to the ligation. The sural nerve was left intact. Sham surgery involved exposing the nerve without damaging it [71, 72].
Calibrated monofilaments (Stoelting Co., Wood Dale, IL) were applied to the plantar surface of the hindpaw and the 50% threshold to withdraw (grams) was calculated as previously described . The stimulus intensity ranged from 0.008 g to 4 g.
A modified version of the acetone drop test  was used: total duration of acetone evoked behaviors (flinching, licking or biting) was measured for 1 minute after acetone (~25 μl) was applied to the plantar surface of the hindpaw.
The accelerating rotarod assay was used (IITC Life Science Inc., Woodland Hills, CA) with the mouse adapter . The task includes a speed ramp from 0 to 30 rpm over 60 s, followed by an additional 240 s at the maximal speed.
Animals were sacrificed 6 months after nerve injury or sham surgery by decapitation following isoflurane anesthesia. Anatomical regions were defined according to the stereotaxic coordinates (rostral–caudal, medial–lateral and dorsal–ventral from bregma) by Paxinos and Franklin . The prefrontal cortex (right and left; +1 to +3, -1 to +1, 0 to −2.5) was extracted, frozen on dry ice and stored at −80 C until use.
RNA extraction and RNA sequencing
RNA was extracted with Trizol (Invitrogen) according to the manufacturer’s protocol and treated with DNAse. cDNA libraries and RNA sequencing was performed by Genome Quebec on the Illumina Genome Analyzer IIX following Illumina guidelines.
ug of RNA was used for cDNA conversion using random hexamers (Roche Molecular Biochemicals) according to manufacturer’s guidelines. Expression primers were designed over overlapping exons with Autoprime.de  or specific to single exons that were differentially expressed in RNAseq and are shown in Additional file 3: Table S3. All reactions for all subjects (n=6-8 per group) were performed on the LightCycler 480 (Roche) in triplicate and statistical significance was determined as P<0.05 using one-tailed t-tests. Quantitative PCR was amplified with a pre-incubation at 95°C for 10 minutes followed by 45 cycles of [95°C for 10 seconds, 60°C for 10 seconds, 72°C for 10 sec] followed by 10 minutes of 72°C. Expression was measured relative to GAPDH expression, which demonstrated stable unchanged expression from RNAseq results.
Two animals were sequenced per condition with 36 bp reads with single end sequencing. RNAseq reads were aligned to the genome (mm9) using Top Hat [59, 78]. Between 36 and 40 million reads aligned to unique genomic locations. Differential expression was determined from read counts per transcript, exon and 1000 bp genomic partition using Bioconductor  package DESeq  with default parameter settings. Transcript annotations were obtained from Ensembl version 63.37 (http://www.ensembl.org/). Concordance of read counts per 1000bp genomic partition between pairs of samples ranged between 0.88 to 0.99.
Ingenuity Software was used to perform whole pathway analysis in the identification of affected networks and their relationship to each other based on the differential expression between SNI and Sham treated mice. Briefly, our data set containing gene identifiers and corresponding expression values was uploaded into the application. Each identifier was mapped to its corresponding object in the Ingenuity® Knowledge Base. Differentially expressed genes, called Network Eligible molecules, were overlaid onto a global molecular network developed from information contained in the Ingenuity Knowledge Base. Networks of Network Eligible Molecules were then algorithmically generated based on their connectivity. Pathways presented were chosen from following candidate gene validation (as indicated with asterisks in Table 2). Right-tailed Fisher’s exact test was used to calculate a p-value determining the probability that each pathway, biological function and/or disease assigned to that data set is due to chance alone.
Calcium-activated chloride channel regulator 1
Central nervous system
Glial fibrillary acidic protein
- GRIN 1:
Glutamate receptor, ionotropic, N-methyl D-aspartate 1
Ingenuity pathway analysis
Lipopolysaccharide binding protein
Polymerase chain reaction
Roundabout, axon guidance receptor, homolog 3
Sodium channel, voltage -gated, type I, alpha subunit
Spared nerve injury
X-linked lymphocyte-regulated 4B.
Kuner R: Central mechanisms of pathological pain. Nat Med 2010, 16: 1258–1266. 10.1038/nm.2231
Xiao HS, Huang QH, Zhang FX, Bao L, Lu YJ, Guo C, Yang L, Huang WJ, Fu G, Xu SH: Identification of gene expression profile of dorsal root ganglion in the rat peripheral axotomy model of neuropathic pain. Proc Natl Acad Sci USA 2002, 99: 8360–8365. 10.1073/pnas.122231899
Grace PM, Hurley D, Barratt DT, Tsykin A, Watkins LR, Rolan PE, Hutchinson MR: Harnessing pain heterogeneity and RNA transcriptome to identify blood-based pain biomarkers: a novel correlational study design and bioinformatics approach in a graded chronic constriction injury model. J Neurochem 2012, 122: 976–994. 10.1111/j.1471-4159.2012.07833.x
Poh KW, Yeo JF, Stohler CS, Ong WY: Comprehensive gene expression profiling in the prefrontal cortex links immune activation and neutrophil infiltration to antinociception. J Neurosci 2012, 32: 35–45. 10.1523/JNEUROSCI.2389-11.2012
Tracey I, Bushnell MC: How neuroimaging studies have challenged us to rethink: is chronic pain a disease? J Pain 2009, 10: 1113–1120. 10.1016/j.jpain.2009.09.001
Schweinhardt P, Bushnell MC: Pain imaging in health and disease–how far have we come? J Clin Invest 2010, 120: 3788–3797. 10.1172/JCI43498
Tajerian M, Alvarado S, Millecamps M, Vachon P, Crosby C, Bushnell MC, Szyf M, Stone LS: Peripheral nerve injury is associated with chronic, reversible changes in global DNA methylation in the mouse prefrontal cortex. PLoS One 2013, 8: e55259. 10.1371/journal.pone.0055259
Seminowicz DA, Wideman TH, Naso L, Hatami-Khoroushahi Z, Fallatah S, Ware MA, Jarzem P, Bushnell MC, Shir Y, Ouellet JA, Stone LS: Effective treatment of chronic low back pain in humans reverses abnormal brain anatomy and function. J Neurosci 2011, 31: 7540–7550. 10.1523/JNEUROSCI.5280-10.2011
Metz AE, Yau HJ, Centeno MV, Apkarian AV, Martina M: Morphological and functional reorganization of rat medial prefrontal cortex in neuropathic pain. Proc Natl Acad Sci USA 2009, 106: 2423–2428. 10.1073/pnas.0809897106
Seminowicz DA, Laferriere AL, Millecamps M, Yu JS, Coderre TJ, Bushnell MC: MRI structural brain changes associated with sensory and emotional function in a rat model of long-term neuropathic pain. NeuroImage 2009, 47: 1007–1014. 10.1016/j.neuroimage.2009.05.068
Augustinsson LE, Sullivan L, Sullivan M: Physical, psychologic, and social function in chronic pain patients after epidural spinal electrical stimulation. Spine (Phila Pa 1976) 1986, 11: 111–119. 10.1097/00007632-198603000-00002
Metzker ML: Sequencing technologies - the next generation. Nat Rev Genet 2010, 11: 31–46. 10.1038/nrg2626
Batten ML, Imanishi Y, Maeda T, Tu DC, Moise AR, Bronson D, Possin D, Van Gelder RN, Baehr W, Palczewski K: Lecithin-retinol acyltransferase is essential for accumulation of all-trans-retinyl esters in the eye and in the liver. J Biol Chem 2004, 279: 10422–10432.
Shiao T, Tran P, Siegel D, Lee J, Vasiliou V: Four amino acid changes are associated with the Aldh3a1 locus polymorphism in mice which may be responsible for corneal sensitivity to ultraviolet light. Pharmacogenetics 1999, 9: 145–153.
Bulgakova NA, Knust E: The Crumbs complex: from epithelial-cell polarity to retinal degeneration. J Cell Sci 2009, 122: 2587–2596. 10.1242/jcs.023648
Mercer TR, Gerhardt DJ, Dinger ME, Crawford J, Trapnell C, Jeddeloh JA, Mattick JS, Rinn JL: Targeted RNA sequencing reveals the deep complexity of the human transcriptome. Nat Biotechnol 2012, 30: 99–104.
Wood SH, Craig T, Li Y, Merry B, de Magalhaes JP: Whole transcriptome sequencing of the aging rat brain reveals dynamic RNA changes in the dark matter of the genome. Age (Dordr) 2012.
Wilson BA, Masel J: Putatively noncoding transcripts show extensive association with ribosomes. Genome Biol Evol 2011, 3: 1245–1252. 10.1093/gbe/evr099
Kishore S, Stamm S: The snoRNA HBII-52 regulates alternative splicing of the serotonin receptor 2C. Science 2006, 311: 230–232. 10.1126/science.1118265
Yao MJ, Chen G, Zhao PP, Lu MH, Jian J, Liu MF, Yuan XB: Transcriptome analysis of microRNAs in developing cerebral cortex of rat. BMC Genomics 2012, 13: 232. 10.1186/1471-2164-13-232
Nicholson B, Verma S: Comorbidities in chronic neuropathic pain. Pain Med 2004, 5(Suppl 1):S9-S27.
Choi JM, Padmala S, Pessoa L: Impact of state anxiety on the interaction between threat monitoring and cognition. NeuroImage 2012, 59: 1912–1923. 10.1016/j.neuroimage.2011.08.102
Drevets WC: Functional anatomical abnormalities in limbic and prefrontal cortical structures in major depression. Prog Brain Res 2000, 126: 413–431.
Drummond SP, Brown GG, Gillin JC, Stricker JL, Wong EC, Buxton RB: Altered brain response to verbal learning following sleep deprivation. Nature 2000, 403: 655–657. 10.1038/35001068
Zlatic M, Landgraf M, Bate M: Genetic specification of axonal arbors: atonal regulates robo3 to position terminal branches in the Drosophila nervous system. Neuron 2003, 37: 41–51. 10.1016/S0896-6273(02)01131-5
Barber M, Di Meglio T, Andrews WD, Hernandez-Miranda LR, Murakami F, Chedotal A, Parnavelas JG: The role of Robo3 in the development of cortical interneurons. Cereb Cortex 2009, 19(Suppl 1):i22–31.
Couch JA, Chen J, Rieff HI, Uri EM, Condron BG: robo2 and robo3 interact with eagle to regulate serotonergic neuron differentiation. Development 2004, 131: 997–1006. 10.1242/dev.00962
Raefski AS, O'Neill MJ: Identification of a cluster of X-linked imprinted genes in mice. Nat Genet 2005, 37: 620–624. 10.1038/ng1567
Cubelos B, Sebastian-Serrano A, Beccari L, Calcagnotto ME, Cisneros E, Kim S, Dopazo A, Alvarez-Dolado M, Redondo JM, Bovolenta P: Cux1 and Cux2 regulate dendritic branching, spine morphology, and synapses of the upper layer neurons of the cortex. Neuron 2010, 66: 523–535. 10.1016/j.neuron.2010.04.038
Gruber AD, Elble RC, Ji HL, Schreur KD, Fuller CM, Pauli BU: Genomic cloning, molecular characterization, and functional analysis of human CLCA1, the first human member of the family of Ca2+−activated Cl- channel proteins. Genomics 1998, 54: 200–214. 10.1006/geno.1998.5562
Hartzell C, Putzier I, Arreola J: Calcium-activated chloride channels. Annu Rev Physiol 2005, 67: 719–758. 10.1146/annurev.physiol.67.032003.154341
Wahl AS, Buchthal B, Rode F, Bomholt SF, Freitag HE, Hardingham GE, Ronn LC, Bading H: Hypoxic/ischemic conditions induce expression of the putative pro-death gene Clca1 via activation of extrasynaptic N-methyl-D-aspartate receptors. Neuroscience 2009, 158: 344–352. 10.1016/j.neuroscience.2008.06.018
Moayedi M, Weissman-Fogel I, Salomons TV, Crawley AP, Goldberg MB, Freeman BV, Tenenbaum HC, Davis KD: Abnormal gray matter aging in chronic pain patients. Brain Res 2012, 1456: 82–93.
As-Sanie S, Harris RE, Napadow V, Kim J, Neshewat G, Kairys A, Williams D, Clauw DJ, Schmidt-Wilcke T: Changes in regional gray matter volume in women with chronic pelvic pain: a voxel-based morphometry study. Pain 2012, 153: 1006–1014. 10.1016/j.pain.2012.01.032
Szabo N, Kincses ZT, Pardutz A, Tajti J, Szok D, Tuka B, Kiraly A, Babos M, Voros E, Bomboi G: White matter microstructural alterations in migraine: a diffusion-weighted MRI study. Pain 2012, 153: 651–656. 10.1016/j.pain.2011.11.029
Li W, Dai S, An J, Li P, Chen X, Xiong R, Liu P, Wang H, Zhao Y, Zhu M: Chronic but not acute treatment with caffeine attenuates traumatic brain injury in the mouse cortical impact model. Neuroscience 2008, 151: 1198–1207. 10.1016/j.neuroscience.2007.11.020
Kim SK, Nabekura J: Rapid synaptic remodeling in the adult somatosensory cortex following peripheral nerve injury and its association with neuropathic pain. J Neurosci 2011, 31: 5477–5482. 10.1523/JNEUROSCI.0328-11.2011
Kim SK, Eto K, Nabekura J: Synaptic structure and function in the mouse somatosensory cortex during chronic pain: in vivo two-photon imaging. Neural Plast 2012, 2012: 640259.
Eto K, Wake H, Watanabe M, Ishibashi H, Noda M, Yanagawa Y, Nabekura J: Inter-regional contribution of enhanced activity of the primary somatosensory cortex to the anterior cingulate cortex accelerates chronic pain behavior. J Neurosci 2011, 31: 7631–7636. 10.1523/JNEUROSCI.0946-11.2011
Apkarian AV, Sosa Y, Krauss BR, Thomas PS, Fredrickson BE, Levy RE, Harden RN, Chialvo DR: Chronic pain patients are impaired on an emotional decision-making task. Pain 2004, 108: 129–136. 10.1016/j.pain.2003.12.015
Gore M, Sadosky A, Stacey BR, Tai KS, Leslie D: The burden of chronic low back pain: clinical comorbidities, treatment patterns, and health care costs in usual care settings. Spine (Phila Pa 1976) 2012, 37: E668–677. 10.1097/BRS.0b013e318241e5de
Liu M, Wood JN: The roles of sodium channels in nociception: implications for mechanisms of neuropathic pain. Pain Med 2011, 12(Suppl 3):S93–99.
Dib-Hajj SD, Cummins TR, Black JA, Waxman SG: Sodium channels in normal and pathological pain. Annu Rev Neurosci 2010, 33: 325–347. 10.1146/annurev-neuro-060909-153234
Escayg A, MacDonald BT, Meisler MH, Baulac S, Huberfeld G, An-Gourfinkel I, Brice A, LeGuern E, Moulard B, Chaigne D: Mutations of SCN1A, encoding a neuronal sodium channel, in two families with GEFS+2. Nat Genet 2000, 24: 343–345. 10.1038/74159
Dichgans M, Freilinger T, Eckstein G, Babini E, Lorenz-Depiereux B, Biskup S, Ferrari MD, Herzog J, van den Maagdenberg AM, Pusch M, Strom TM: Mutation in the neuronal voltage-gated sodium channel SCN1A in familial hemiplegic migraine. Lancet 2005, 366: 371–377. 10.1016/S0140-6736(05)66786-4
Meier S, Demirakca T, Brusniak W, Wolf I, Liebsch K, Tunc-Skarka N, Nieratschker V, Witt SH, Matthaus F, Ende G: SCN1A Affects Brain Structure and the Neural Activity of the Aging Brain. Biol Psychiatry 2012, 72: 677–683. 10.1016/j.biopsych.2012.03.017
Eijkelkamp N, Linley JE, Baker MD, Minett MS, Cregg R, Werdehausen R, Rugiero F, Wood JN: Neurological perspectives on voltage-gated sodium channels. Brain 2012, 135: 2585–2612. 10.1093/brain/aws225
Hains BC, Saab CY, Klein JP, Craner MJ, Waxman SG: Altered sodium channel expression in second-order spinal sensory neurons contributes to pain after peripheral nerve injury. J Neurosci 2004, 24: 4832–4839. 10.1523/JNEUROSCI.0300-04.2004
Cheah CS, Yu FH, Westenbroek RE, Kalume FK, Oakley JC, Potter GB, Rubenstein JL, Catterall WA: Specific deletion of NaV1.1 sodium channels in inhibitory interneurons causes seizures and premature death in a mouse model of Dravet syndrome. Proc Natl Acad Sci USA 2012, 109: 14646–14651. 10.1073/pnas.1211591109
Petrenko AG, Perin MS, Davletov BA, Ushkaryov YA, Geppert M, Sudhof TC: Binding of synaptotagmin to the alpha-latrotoxin receptor implicates both in synaptic vesicle exocytosis. Nature 1991, 353: 65–68. 10.1038/353065a0
Bragina L, Fattorini G, Giovedi S, Melone M, Bosco F, Benfenati F, Conti F: Analysis of Synaptotagmin, SV2, and Rab3 Expression in Cortical Glutamatergic and GABAergic Axon Terminals. Front Cell Neurosci 2011, 5: 32.
Ji G, Neugebauer V: Pain-related deactivation of medial prefrontal cortical neurons involves mGluR1 and GABA(A) receptors. J Neurophysiol 2011, 106: 2642–2652. 10.1152/jn.00461.2011
Muta T, Takeshige K: Essential roles of CD14 and lipopolysaccharide-binding protein for activation of toll-like receptor (TLR)2 as well as TLR4 Reconstitution of TLR2- and TLR4-activation by distinguishable ligands in LPS preparations. Eur J Biochem 2001, 268: 4580–4589. 10.1046/j.1432-1327.2001.02385.x
Wei L, Simen A, Mane S, Kaffman A: Early life stress inhibits expression of a novel innate immune pathway in the developing hippocampus. Neuropsychopharmacology 2012, 37: 567–580. 10.1038/npp.2011.239
Garate I, Garcia-Bueno B, Madrigal JL, Bravo L, Berrocoso E, Caso JR, Mico JA, Leza JC: Origin and consequences of brain Toll-like receptor 4 pathway stimulation in an experimental model of depression. J Neuroinflammation 2011, 8: 151. 10.1186/1742-2094-8-151
Chen PE, Errington ML, Kneussel M, Chen G, Annala AJ, Rudhard YH, Rast GF, Specht CG, Tigaret CM, Nassar MA: Behavioral deficits and subregion-specific suppression of LTP in mice expressing a population of mutant NMDA receptors throughout the hippocampus. Learn Mem 2009, 16: 635–644. 10.1101/lm.1316909
Sakurai K, Toru M, Yamakawa-Kobayashi K, Arinami T: Mutation analysis of the N-methyl-D-aspartate receptor NR1 subunit gene (GRIN1) in schizophrenia. Neurosci Lett 2000, 296: 168–170. 10.1016/S0304-3940(00)01599-8
Tarabeux J, Kebir O, Gauthier J, Hamdan FF, Xiong L, Piton A, Spiegelman D, Henrion E, Millet B, Fathalli F: Rare mutations in N-methyl-D-aspartate glutamate receptors in autism spectrum disorders and schizophrenia. Transl Psychiatry 2011, 1: e55. 10.1038/tp.2011.52
Hamdan FF, Gauthier J, Araki Y, Lin DT, Yoshizawa Y, Higashi K, Park AR, Spiegelman D, Dobrzeniecka S, Piton A: Excess of de novo deleterious mutations in genes associated with glutamatergic systems in nonsyndromic intellectual disability. Am J Hum Genet 2011, 88: 306–316. 10.1016/j.ajhg.2011.02.001
Labrie V, Clapcote SJ, Roder JC: Mutant mice with reduced NMDA-NR1 glycine affinity or lack of D-amino acid oxidase function exhibit altered anxiety-like behaviors. Pharmacol Biochem Behav 2009, 91: 610–620. 10.1016/j.pbb.2008.09.016
Furuse T, Wada Y, Hattori K, Yamada I, Kushida T, Shibukawa Y, Masuya H, Kaneda H, Miura I, Seno N: Phenotypic characterization of a new Grin1 mutant mouse generated by ENU mutagenesis. Eur J Neurosci 2010, 31: 1281–1291. 10.1111/j.1460-9568.2010.07164.x
Zhang X, Wang J, Zhou Q, Xu Y, Pu S, Wu J, Xue Y, Tian Y, Lu J, Jiang W, Du D: Brain-derived neurotrophic factor-activated astrocytes produce mechanical allodynia in neuropathic pain. Neuroscience 2011, 199: 452–460.
Piao ZG, Cho IH, Park CK, Hong JP, Choi SY, Lee SJ, Lee S, Park K, Kim JS, Oh SB: Activation of glia and microglial p38 MAPK in medullary dorsal horn contributes to tactile hypersensitivity following trigeminal sensory nerve injury. Pain 2006, 121: 219–231. 10.1016/j.pain.2005.12.023
Chen FL, Dong YL, Zhang ZJ, Cao DL, Xu J, Hui J, Zhu L, Gao YJ: Activation of astrocytes in the anterior cingulate cortex contributes to the affective component of pain in an inflammatory pain model. Brain Res Bull 2012, 87: 60–66. 10.1016/j.brainresbull.2011.09.022
Lu Y, Zhu L, Gao YJ: Pain-related aversion induces astrocytic reaction and proinflammatory cytokine expression in the anterior cingulate cortex in rats. Brain Res Bull 2011, 84: 178–182. 10.1016/j.brainresbull.2010.12.007
Kuzumaki N, Narita M, Hareyama N, Niikura K, Nagumo Y, Nozaki H, Amano T, Suzuki T: Chronic pain-induced astrocyte activation in the cingulate cortex with no change in neural or glial differentiation from neural stem cells in mice. Neurosci Lett 2007, 415: 22–27. 10.1016/j.neulet.2006.12.057
Miguel-Hidalgo JJ, Waltzer R, Whittom AA, Austin MC, Rajkowska G, Stockmeier CA: Glial and glutamatergic markers in depression, alcoholism, and their comorbidity. J Affect Disord 2010, 127: 230–240. 10.1016/j.jad.2010.06.003
Toro CT, Hallak JE, Dunham JS, Deakin JF: Glial fibrillary acidic protein and glutamine synthetase in subregions of prefrontal cortex in schizophrenia and mood disorder. Neurosci Lett 2006, 404: 276–281. 10.1016/j.neulet.2006.05.067
Zhang F, Vadakkan KI, Kim SS, Wu LJ, Shang Y, Zhuo M: Selective activation of microglia in spinal cord but not higher cortical regions following nerve injury in adult mouse. Mol Pain 2008, 4: 15. 10.1186/1744-8069-4-15
Svensson CI, Brodin E: Spinal astrocytes in pain processing: non-neuronal cells as therapeutic targets. Mol Interv 2010, 10: 25–38. 10.1124/mi.10.1.6
Decosterd I, Woolf CJ: Spared nerve injury: an animal model of persistent peripheral neuropathic pain. Pain 2000, 87: 149–158. 10.1016/S0304-3959(00)00276-1
Bourquin AF, Suveges M, Pertin M, Gilliard N, Sardy S, Davison AC, Spahn DR: Decosterd I: Assessment and analysis of mechanical allodynia-like behavior induced by spared nerve injury (SNI) in the mouse. Pain 2006, 122(14):e11–14.
Chaplan SR, Bach FW, Pogrel JW, Chung JM, Yaksh TL: Quantitative assessment of tactile allodynia in the rat paw. J Neurosci Methods 1994, 53: 55–63. 10.1016/0165-0270(94)90144-9
Choi Y, Yoon YW, Na HS, Kim SH, Chung JM: Behavioral signs of ongoing pain and cold allodynia in a rat model of neuropathic pain. Pain 1994, 59: 369–376. 10.1016/0304-3959(94)90023-X
Jones BJ, Roberts DJ: The quantiative measurement of motor inco-ordination in naive mice using an acelerating rotarod. J Pharm Pharmacol 1968, 20: 302–304. 10.1111/j.2042-7158.1968.tb09743.x
Paxinos G, Franklin K: The mouse brain in stereotaxic coordinates. Waltham Massachusetts: Academic Press; 2004.
Wrobel GKF: Lichter: Autoprime: Selecting primers for expressed sequences. Genome Biol 2004, 5: 11.
Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 2009, 25: 1105–1111. 10.1093/bioinformatics/btp120
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol 2004, 5: R80. 10.1186/gb-2004-5-10-r80
Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol 2010, 11: R106. 10.1186/gb-2010-11-10-r106
This work was funded by CIHR MOP-102586 and, a FRSQ bourse de chercheur-boursier to LSS, a grant to MS from the CIHR MOP-42411 a gift from the Sackler McGill Program in Psychobiology and Epigenetics (MS) and a grant from the Canadian Institute for Advanced Research (MS). MT was supported by a studentship from The Louise and Alan Edwards Foundation.
The authors declare that they have no competing interests.
SA and MT Overall experimental design, collecting the mouse behavioral data, performing RNA data collection and analysis and drafting the manuscript. MM Induction of nerve injury, collection of brain tissue Suderman: Bioinformatic analysis of RNAseq data. LSS Overall design of the experiments, supervision of the experiments and writing of the manuscript. Szyf: Overall design of the experiments, supervision of the experiments and writing of the manuscript. All authors have read and approved the final manuscript.
Sebastian Alvarado, Maral Tajerian contributed equally to this work.
Electronic supplementary material
Additional file 2: Table S2: Ingenuity Pathway Analysis of altered biological functions in Sham vs. SNI PFC. (XLSX 18 KB)
Authors’ original submitted files for images
About this article
Cite this article
Alvarado, S., Tajerian, M., Millecamps, M. et al. Peripheral nerve injury is accompanied by chronic transcriptome-wide changes in the mouse prefrontal cortex. Mol Pain 9, 21 (2013). https://doi.org/10.1186/1744-8069-9-21
- Neuropathic pain
- Chronic pain
- RNA sequencing
- Noncoding RNA
- Neuronal plasticity