Skip to main content


Transcriptional mechanisms underlying sensitization of peripheral sensory neurons by Granulocyte-/Granulocyte-macrophage colony stimulating factors

Article metrics



Cancer-associated pain is a major cause of poor quality of life in cancer patients and is frequently resistant to conventional therapy. Recent studies indicate that some hematopoietic growth factors, namely granulocyte macrophage colony stimulating factor (GMCSF) and granulocyte colony stimulating factor (GCSF), are abundantly released in the tumor microenvironment and play a key role in regulating tumor-nerve interactions and tumor-associated pain by activating receptors on dorsal root ganglion (DRG) neurons. Moreover, these hematopoietic factors have been highly implicated in postsurgical pain, inflammatory pain and osteoarthritic pain. However, the molecular mechanisms via which G-/GMCSF bring about nociceptive sensitization and elicit pain are not known.


In order to elucidate G-/GMCSF mediated transcriptional changes in the sensory neurons, we performed a comprehensive, genome-wide analysis of changes in the transcriptome of DRG neurons brought about by exposure to GMCSF or GCSF. We present complete information on regulated genes and validated profiling analyses and report novel regulatory networks and interaction maps revealed by detailed bioinformatics analyses. Amongst these, we validate calpain 2, matrix metalloproteinase 9 (MMP9) and a RhoGTPase Rac1 as well as Tumor necrosis factor alpha (TNFα) as transcriptional targets of G-/GMCSF and demonstrate the importance of MMP9 and Rac1 in GMCSF-induced nociceptor sensitization.


With integrative approach of bioinformatics, in vivo pharmacology and behavioral analyses, our results not only indicate that transcriptional control by G-/GMCSF signaling regulates a variety of established pain modulators, but also uncover a large number of novel targets, paving the way for translational analyses in the context of pain disorders.


Pain is one of the most severe and common symptoms of a variety of cancers and is a primary determinant of the poor quality of life in cancer patients. In a large number of clinical cases, cancer-associated pain, particularly the neuropathic component thereof, is resistant to conventional therapeutics or their application is severely limited owing to the widespread side effects. Because many types of carcinomas and sarcomas metastasize to skeletal bones, they are associated with spontaneous pain, hyperalgesia and allodynia. As potential mechanisms, tumor-derived factors, such as NGF[1], endothelins[24], amongst others, have been studied, which either directly activate nociceptive nerves or sensitize them towards sensory stimuli[5, 6].

Several types of non-hematopoietic tumors secrete hematopoietic colony stimulating factors, which act on myeloid cells and tumor cells[7]. In a recent study, we demonstrated that receptors and signaling mediators of granulocyte- and granulocyte-macrophage colony stimulating factors (G-/GMCSF) are also broadly expressed on sensory nerves in mouse models of bone metastases as well as in human biopsies of pancreatic adenocarcinoma[8]. Using animal models of bone metastases which closely mimic the nature and progression of cancer pain in humans, we reported that GCSF and GMCSF directly act on receptors on diverse DRG neurons to subserve important functions in the generation of pain hypersensitivity in tumor-affected regions[8]. Importantly, behavioral, electrophysiological and biochemical experiments demonstrated sensitization of sensory nerves towards thermal and mechanical stimuli as well as an increase in neurotransmitter release upon exposure to G-/GMCSF. By adapting RNAi methodology in vivo, we demonstrated that a specific loss of GMCSFRα in DRG led to a reduction in bone tumor–evoked pain without interfering with the tumor growth, indicating that GMCSF signaling in peripheral nerves contributes substantially to cancer pain[8]. Recent studies on post-surgical pain and inflammatory pain also point to a key role for these cytokines[912].

G-/GMCSF activates the JAK family of receptor tyrosine kinases, which unfolds its activity by not only regulating enzymes and target proteins within its local milieu, but importantly also by activating the STAT family of transcription factors, which subsequently dimerize and translocate to the cell nucleus to regulate gene expression[13]. Albeit we have reported local, acute activation of the ERK Kinase as well as PI3 Kinase in sensory nerves upon a short-term exposure to G-/GMCSF, nothing is known so far about the nature of genes regulated transcriptionally in DRG neurons upon exposure to G/GMCSF. However, long-term transcriptional mechanisms of G/GMCSF action are arguably of even greater importance in pathophysiological states involving chronic, continual release of G/GMCSF, such as tumor-affected tissues, rheumatoid arthritis, amongst others[14, 15]. Addressing precise mechanisms via which the G-/GMCSF-JAK-STAT pathway elicits long-term nociceptive sensitization is thus important for understanding mechanisms of cancer pain and other chronic disorders associated with G-/GMCSF release.

In lieu of the attractive therapeutic opportunities offered by these findings, we aimed to elucidate cellular targets of G/GMCSFR in DRG neurons, particularly with respect to transcriptional regulation. Not only did we find a variety of known, established ‘pain-related’ mediators to be transcriptional targets of G-/GMCSF, but also several protein-protein interaction hubs were observed to be under G-/GMCSF regulation in sensory neurons via detailed bioinformatics analyses. Behavioral and pharmacological analyses on 4 of the emerging targets confirmed that Rac1 and Matrix metallopeptidase 9 (MMP9) contribute to GMCSF-induced nociceptive sensitization. These integrative approaches advance our understanding of chronic pain mechanisms and hold promise in the development of novel therapeutic approaches.

Materials and methods

Animal usage

All animal usage procedures were in accordance with ethical guidelines laid down by the International Association of the Study of Pain and the local governing body (Regierungspräsidium Karlsruhe). All behavioral measurements were done in awake, unrestrained, age-matched adult (more than 2 months-old) C57/Bl6 mice. Mice were housed in plastic cages, with ambient temperature and a 12 h diurnal light cycle. Food and water were provided ad libitum.

Sensory neuronal cultures and G-/GMCSF treatment

Adult DRG neuronal cultures were prepared following the protocol explained previously[8]. Briefly, neuronal cells isolated from adult wild type mice were seeded on Poly-L-Lysine coated cover slips and maintained in F12 Media (Sigma) supplemented with 15% Amino Acids (Gibco), 10% bovine serum (Invitrogen), 1% Penicillin/Streptomycin (Gibco), 0.5% L-Glutamine (Gibco) and Nerve Growth Factor (100 mg ml-1, Roche). 4 days old enriched adult neuronal cultures were starved of growth factors and serum for 4 h. At the end of 4 h, starving culture medium was replaced with medium containing 0.5% Fetal bovine serum (Life Technologies, 10270106) together with either 1× PBS (vehicle) or 200 ng/mL of murine GMCSF (Peprotech, 315–03) or 200 ng/mL of murine GCSF (Peprotech 250–05) dissolved in 1×PBS. Neurons were left in the incubator for 24 h. Each treatment was performed in triplicate culture wells (n = 3) to test biological variability. At the end of 24 h, total RNA was isolated and used for microarray expression or qRT-PCR analysis.

RNA isolation from cultured sensory neurons and DRGs

Total RNA from cultured sensory neurons treated with murine GMCSF or GCSF or PBS for 24 h was isolated using mirVana™ miRNA Isolation Kit (Ambion, AM 1561) following manufacturer’s instructions and dissolved in 20 μl of nuclease-free water. Purification steps were performed using RNAse-free DNAse kit (Qiagen, 79254) following manufacturer’s instructions. RNA concentration was determined using the NanoDrop spectrophotometer (NanoDrop Technologies, Wilmington, Germany) and the quality of total RNA was checked by gel analysis using the total RNA Nanochip assay on an Agilent 2100 Bioanalyzer (Agilent Technologies GmbH, Waldbronn, Germany). Only samples with RNA index values greater than 7 were selected for mRNA profiling. 200 ng of total RNA from each biological sample was used as starting material for mRNA expression analysis.

For in vivo testing, lumbar DRGs L3, L4 and L5 were collected at 25 h, 36 and 48 h after bilateral intraplantar application of 20 ng murine GMCSF and flash frozen in liquid nitrogen. Total RNA was isolated and processed following the same protocol explained above for cultured sensory neurons.

Microarray expression, networking and gene ontology analysis

The mRNA profiling was performed on polyadenylated RNA using Illumina mouse sentrix-6 chips. cDNA library preparation, hybridization and scanning steps were performed by employing in-house standardized protocols and including stringent positive and negative controls at each step at the genomics and proteomics core facility, German Cancer Research Centre (DKFZ, Heidelberg, Germany). The array intensity data were imported into Beadstudio ver. 3 from Illumina and the quantile array normalization method was employed to account for intra- and inter-array variations in expression intensities within each experimental group[16]. Magnitude of induction or repression of individual transcript was compared over vehicle-treated samples. To be able to understand the magnitude of regulation at transcript level, we first converted probe-level signals to transcript-level signals by using in house developed Perl scripts and applying following criteria – i) Take the average fold-change if all the probes for one transcript showed the same direction (positive or negative) of regulation, ii) Discard the transcripts for which different probes showed different directions of regulation iii) Take the regulation value from the majority of probes if only one probe out of several probes is showing different expression signals.

All the networking analyses of the expression data were performed using MetaCore™ software in which a network is built around an initial list of “seed nodes,” which can originate from the uploaded experiment, or be manually assembled, or else be automatically converted by MetaCore™ from a list of genes. For the gene ontology enrichment analyses, we used software called bioCompendium ( developed at the European Molecular Biology Laboratory, Heidelberg, Germany.

GMCSF and inhibitors application in vivo

Murine GMCSF was purchased from Peprotech, dissolved in 1× PBS of physiological PH and 20 ng was applied into the intraplantar surface of adult C57/Bl6 mice unilaterally for 4 times at 8 h intervals. Inhibitors for MMP-9 (CAS 1177749-58-4) and Rac1 (1090893-12-1) were purchased from Calbiochem and dissolved in 10% and 50% Dimethyl Sulfoxide (DMSO, Applichem), respectively. Calpain 1/2 inhibitor (Calpain inhibitor III, A3672, 0250) was purchased from Sigma-Aldrich and dissolved in 20% DMSO. TNFα inhibitor was purchased from Pfizer (Enbrel®)) and diluted in 1× PBS. One hour after the last GMCSF dosage application, different groups of mice received 0.15 pmoles, 1.5 nmoles, 10 nmoles or 100 pmoles of MMP-9 or Rac1 or Calpain 1/2 or TNFα inhibitors, respectively, in 10 μl volume of vehicle into the same paw into which multiple dosages of GMCSF were applied. BSA dissolved in 1× PBS was used as vehicle control for TNFα inhibition experiments. Mechanical hyperalgesia was recorded after 4 and 8 h after the last GMCSF dosage application while thermal hypersensitivity was recorded after 5 and 9 h of the last GMCSF dosage application.

Mechanical and thermal pain behavioral tests

Mice were habituated to the experimental setup in at least 2 separate sessions within the week preceding the time of behavioral testing. The observer was fully blinded to the identity of the groups in all behavioral tests. To measure mechanical sensitivity, animals were placed on an elevated wire grid and the plantar hind paw was stimulated using calibrated von Frey monofilaments of 0.07 g, 0.16 g, 0.4 g and 1.0 g strength (Bioseb, France). Paw withdrawal was recorded as a positive response. Data is expressed as percentage of frequency of response over 5 stimulations and data from representative filament is shown in this manuscript. For thermal nociceptive testing, radiant heat was applied using Hargreaves’ apparatus (Ugo Basile, Italy) to the plantar surface of the hind paw, until mice retract it sharply. The time taken to retract (paw withdrawal latency) the hind paw was recorded. A cut-off of 15 seconds heat exposure was followed in order to avoid any potential damage to the tissue.

Quantification of mRNA expression

We used NanoString-nCounter™ based gene quantification method to validate microarray expression data. Probes specifically targeting the desired gene of interest were obtained from Nanostring Technologies, USA and analyses were performed at the nCounter core facility of the Medical Faculty of Heidelberg, Heidelberg University, Germany. Two hundred ng of total RNA were used to analyze the expression of diverse target genes, using 5 housekeeping genes, namely Clathrin, heavy polypeptide (Cltc), Glyceraldehyde-3-phosphate dehydrogenase (Gapdh), glucuronidase beta (Gusb), Hypoxanthine guanine phosphoribosyl transferase (Hprt) and Tubulin, beta 5 class I (Tubb5), as internal controls. Expression of target genes was analyzed by comparing treated and control samples. Fold-change of test gene was expressed as arithmetic average value over all 5 housekeeping genes.

Taqman assays (Life Technologies, USA) were used for QRTPCR-based quantification of Rac1 (assay ID Mm01201653_mH), Calpain2 (assay ID Mm00486669_m1), MMP9 (assay ID Mm00600164_g1) and TNFα (assay ID Mm00443260_g1). 20 ng of total RNA was used to prepare the cDNA using random primers from the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, 4368814) following manufacturer’s instructions. Four μl of prepared cDNA were PCR amplified in each reaction using mRNA-specific primers and TaqMan® Universal Master Mix II, (Applied Biosystems, 4440040) following manufacturer’s instructions on Chromo 4 detection system (BioRad, USA). The expression level of the target mRNA was normalized to the expression of Glyceraldehyde 3-phosphate dehydrogenase (GAPDH, assay ID 4352932E, Applied biosystems). Each mRNA was amplified from triplicate samples and Ct values were recorded. Fold-change in the mRNA expression in vehicle- or GMCSF-treated sensory neuronal cultures was calculated using D DC T method[17] which measures the relative change in expression of a mRNA from treatment to control compared to the reference gene.

Data analysis

All data are presented as mean ± standard error of the mean (S.E.M.), Two-tailed Student’s t-test or the Analysis of Variance (ANOVA) for repeated measures followed by post-hoc Fisher’s LSD test was utilized to determine statistically significant differences (p < 0.05), unless mentioned otherwise for a particular experiment.


GMCSF-mediated changes in the gene expression repertoire in sensory neurons

To investigate transcriptional expression changes caused by exposure to GMCSF or GCSF application at a genome-wide level, we performed a genome-wide gene profiling screen from cultured DRG neurons derived from adult mice. Neuron-enriched cultures were starved of growth factors and serum for 4 h and treated with GMCSF or GCSF (200 ng/mL in PBS) or vehicle (PBS) in medium containing 0.5% serum for 24 h. Total RNA isolated from 3 such independent experiments was subjected to quality control as described under methods and processed in profiling experiments using Illumina prespotted arrays. cDNA library preparation, hybridization and scanning steps were performed by employing in-house standardized protocols, and including stringent positive and negative controls at each step. The array intensity data were imported into Beadstudio ver. 3 from Illumina and the quantile array normalization method was employed to correct for systematic differences between arrays which do not represent a biological variation of interest between experimental groups[18]. Magnitude of fold-induction or fold-repression of individual transcript was compared over vehicle-treated samples, as described in our previous studies[19].

Data acquired from the array provided expression signals from 46090 probes targeting 30723 genome-wide transcripts. Out of 30723 transcripts arrayed 15833 and 16882 showed significant fold change in their expression in GMCSF-treated and GCSF-treated samples, respectively, as compared to vehicle-treated samples (P ≤ 0.05, two-tailed t-test assuming equal variance, Benjamini and Hochberg false discovery rate correction, n = 3 independent profiling experiments) (Figure 1 and Additional file1: Table S1). Because both GMCSF and GCSF act in a pronociceptive manner, we then studied commonly regulated transcripts and found that 3898 transcripts showed significant upregulation with GMCSF as well as with GCSF as stimuli for DRG neurons (Figure 1B). These included several genes which have been implicated in nociceptive modulation, such as chemokine (C-C motif) ligand 2 and 3 (Ccl2 and Ccl3)[20, 21], transient receptor potential cation channel, subfamily V, member 1 (TRPV1), a molecular sensor and transducer for heat, protons and algogens[22], amongst others (a few examples are shown in Figure 1C). Moreover, 9254 transcripts were commonly downregulated upon exposure with GMCSF as well as with GCSF. These also included several pain-related known genes, such as the voltage-dependent calcium channel subunit aplha2/delta1 (Cacna2d1)[23], the AMPA receptor interacting protein, GRIP1[24, 25], amongst others (Figure 1C). Interestingly, however, 421 genes showed reciprocal regulation upon exposure to GMCSF or GCSF (Figure 1B), e.g. the nociceptive modulatory chemokine (C-C motif) ligand (Ccl5)[26] and the matrix metalloprotease 3 (MMP3)[27, 28] were significantly upregulated after GMCSF exposure, but downregulated upon exposure to GCSF (Figure 1C). In general, genes encoding nociceptive modulating pronociceptive chemokines and cytokines appeared more strongly regulated by GMCSF signaling than by GCSF signaling in DRG neurons (some examples in Figure 1C).

Figure 1

GMCSF- or GCSF- mediated gene pool in the peripheral sensory neurons. (A) Heat map representation of significantly-regulated transcripts showing more than two-fold change (two-tailed t-test assuming equal variance, P with Benjamini and Hochberg False Discovery Rate correction <0.05) following chronic exposure to GMCSF or GCSF in sensory neurons of the DRG. The color scale represents quantile normalized hybridization intensities for each gene. (B) Venn diagram representing the number of genes commonly or differentially regulated by GMCSF and GCSF exposure in DRG neurons. (C) Selected pain-related genes that were significantly regulated (two-tailed t-test assuming equal variance, P with Benjamini and Hochberg False Discovery Rate correction <0.05). Genes commonly (green text) or differently (red text) regulated by GMCSF and GCSF stimulation in sensory neurons are shown. Gene regulation confirmed by qRT-PCR analysis is highlighted in bold text.

To test the validity of the microarray data reported above, we performed quantitative measurements of the expression of several candidate regulated genes using Nanostring-nCounter technique. Amongst putatively regulated transcripts in the GMCSF-mediated gene pool, we tested 100 genes quantitatively and found that 78 genes were regulated as predicted by profiling data, which included up-regulated genes such as chemokine (C-C motif) ligand 5 (Ccl5), interleukin 1 alpha (Il1a), Ccl3 (Figure 2A) and down-regulated genes such as Cacna2d1, synapsin II (Syn2), amongst others (Figure 2B). Along the same lines, GCSF-mediated regulation of several genes, including pain-related genes such as calcitonin/calcitonin-related polypeptide, alpha, transcript variant 1(Calca), Ccl3, and fibroblast growth factor 7 (Fgf7) amongst several others could be confirmed using the same PCR-based methodology (Figure 2C,2D). These results thus validate the results obtained with the microarray expression arrays via an independent method.

Figure 2

Nanostring-nCounter-based validation of selected genes those are up-regulated by GMCSF (A), down-regulated by GMCSF (B), up-regulated by GCSF (C) and down-regulated by GCSF (D) in the sensory neurons. Fold-change expression in the genes is expressed as arithmetic average over 5 housekeeping genes namely Cltc, Gapdh, Gusb, Hprt and Tubb5 in all panels. * P < 0.05, one-way ANOVA followed by Fisher’s LSD Post-hoc analysis.

In a next step, to understand systems level interactions in the GMCSF- or GCSF-mediated gene pools, we performed a direct-interactions analysis using Metacore software[19, 29]. When we applied this to all significantly regulated transcripts following the criteria explained above for Figure 1, it yielded too dense a network to allow meaningful interpretations (data not shown). Therefore, we stringently filtered out the transcripts which showed at least 4-fold up- or down-regulation upon exposure to GMCSF (thus arriving at 661 transcripts) or GCSF (611 transcripts). Of these, only 467 GMCSF-target genes and 454 GCSF-target genes were well annotated with known higher level mapping in Metacore and were used for the direct-network analysis. The network map generated by the genego direct-interaction network analysis tool revealed a dense network of genes in the GMCSF-target pool with 3 major nodal points namely, two transcription factors, E26 avian leukemia oncogene 1, 5' domain, transcript variant 2 (Ets1), Hypoxia inducible factor 1 alpha subunit (Hif1a) and a metallo-protease, namely Mmp9 (Figure 3). These 3 nodal points are intensively related to many kinases such as mitogen-activated protein kinase 3 (MAPK3), generic binding proteins such as Synapsin, Ras super family members such as Rac1, receptors like Toll-like receptor 2 encoding gene (Tlr2), all of which are either directly or indirectly implicated in nociceptive mechanisms. Similarly, the direct-interaction network for the GCSF-mediated gene pool also revealed a densely connected network with genes encoding the key posttranslational sumoylation protein (Sumo1), the cyclin-dependent kinase inhibitor 1A (Cdkn1a), CREB binding protein (Crebbp or CBP), calpain 2 (Capn2), MAPK3 and the RhoGTPase Rac1 (Rac1) serving as major nodal points. These nodes are intensively linked to genes encoding Calmodulin 2 (Calm2), the Transient Receptor Family channel V1 (Trpv1 or capsaicin receptor), Actin-modulatory protein profilin 1 (Pfn1), among several others (Figure 4). These results indicate that GMCSF- and GCSF-signaling interlinks transcriptional and post-translational modification mechanisms to key nociceptive modulatory proteins. We further performed direct interaction network analysis on genes that were commonly regulated following GMCSF or GCSF exposure. These commonly-regulated networks revealed shared nodal points such as Rac1, mitogen-activated protein kinase 3 (Mapk3), among others (Additional file2: Figure S1 and Additional file3: Figure S2).

Figure 3

Direct-interactions-network analysis on GMCSF-mediated gene pool with fold-change between +4 and −4 as compared to control-treated sensory neurons and P-BH < 0.05 ( t -test, P with Benjamini and Hochberg False Discovery Rate < 0.05). Genes upregulated and downregulated in GMCSF-dependent manner are marked with red and blue circles, respectively. Please see Additional file 3: Figure S2 for information on legends.

Figure 4

Direct-interactions-network analysis on GCSF-mediated gene pool with fold change between +4 and −4 and as compared to control-treated sensory neurons and P-BH < 0.05 ( t -test, P with Benjamini and Hochberg False Discovery Rate < 0.05). Genes upregulated and downregulated in GCSF-dependent manner are marked with red and blue circles, respectively. Please see Additional file 3: Figure S2 for information on legends.

A gene ontology enrichment analysis on the same subsets of GMCSF- or GCSF-target genes performed using the bioCompendium software revealed that a major proportion of transcripts (> 50%) show protein-binding activity. Interestingly, cytokine activity and chemokine-receptor binding categories were found to be represented in the molecular function enrichment analysis on GMCSF-target pool in DRG neurons (Additional file4: Figure S3), consistent with our observation of high regulation levels of several nociception-related cytokines and chemokines.

In the next step, by using the same subsets of significantly regulated GMCSF- or GCSF-modulated genes as explained for Figure 1A, we performed a network analysis in which networks are built on the basis of relationships and interactions contained in the MetaCore™ Database. Interestingly, the network which emerged with a top rank from the gene pool of GMCSF-targets in sensory neurons revealed that the classical signaling cascade consisting of JAK kinases (JAK1 and JAK2) and STAT transcription factors, STAT1 and STAT3 are tightly linked to Tumor necrosis factor-alpha (TNF-alpha) and its receptor TNF-R1, both of which were observed to be directly regulated by GMCSF in sensory neurons in our profiling analyses (Additional file5: Figure S4). Moreover, a link to NF-kappa-I Kappa B signaling, which has also been implicated in sensory neurons[30], was also apparent (Additional file5: Figure S4). These findings further indicate a close link between GMCSF-induced transcriptional control and induction of key nociceptive modulators, such as TNF-alpha.

Functional significance of GM-/GCSF-regulated gene pool in GM-/GCSF-induced nociceptive hypersensitivity

Finally, to functionally validate our results on GMCSF- and GCSF-associated genes, we selected protein products of a set of four candidate genes from different functional classes and with functional relevance to pain modulation, namely the RhoGTPase Rac1, the matrix metallopeptidase 9 (MMP9), a chemokine TNF-alpha and a generic protease calpain 2. To confirm GMCSF-mediated modulation of these four genes, we compared their mRNA expression in the total RNA isolated from the DRG neuronal cultures following chronic treatment with GMCSF or vehicle, i.e. a similar paradigm as with the expression array screening. Analysis of results confirmed GMCSF-mediated robust upregulation of Rac1, MMP9, TNFα and Calpain 2 as compared to vehicle-treated samples (Figure 5A, *P < 0.05, one-way ANOVA followed by Fisher’s LSD Post-hoc analysis).

Figure 5

QRTPCR-based validation of GM-CSF-induced regulation of the expression of Rac1, MMP9, Calpain2 and TNFα in vitro as well as in vivo . Gene regulation at 24 h following the start of GMCSF exposure in cultured sensory neurons is shown in (A) and in DRGs at different time points following intraplantar GMCSF injection is shown in (B). * denotes P ≤ 0.05 as compared to vehicle group, One-Way ANOVA followed by Fisher’s LSD Post-hoc analysis, n = 3 mice per group. (C) Schematic representation of the protocol followed to investigate changes in GMCSF-mediated mechanical and thermal hypersensitivity upon inhibition of Rac1 or MMP9 or Calpain 2 or TNFα in vivo.

In previous studies, we have analyzed short-term effects of acute exposure to GCSF and GMCSF[8]. However, in order to mimic chronic clinical conditions which are associated with longer exposure to G-/GMCSF and to match the course of the following behavioral experiments with the time frame of our gene regulation studies, we administered multiple dosages of 20 ng murine GMCSF as described in the scheme shown in Figure 5C and inhibitors were applied one hour after the last GMCSF dosage application. Mechanical sensitivity was recorded upon ipsilateral plantar application of graded von Frey filaments at 3 h and 7 h post injection and thermal sensitivity was recorded using Hargreaves’ apparatus at 4 and 8 h following inhibitor application. For in vivo confirmation of GMCSF-mediated regulation of Rac1, MMP9, Calpain2 and TNFα, DRGs were collected at 24, 36 and 48 h after the first GMCSF dosage application. QRTPCR-based expression analysis confirmed GMCSF-mediated up-regulation of Rac1, MMP9, Calpain 2 and TNFα in the DRGs of GMCSF-treated mice as compared to the vehicle-treated group of mice (Figure 5B, *P < 0.05, one-way ANOVA followed by Fisher’s LSD Post-hoc analysis).

Peripheral Rac1 activation is required for GMCSF-mediated mechanical and thermal hyperalgesia

In our profiling analysis and quantitative PCR analysis, Rac1 expression increased significantly in DRG neurons following a 24 h-long exposure to GMCSF or GCSF. In spinal neurons, regulation of Rac1 activity is known to impact dendritic spine morphology and density as well as pain hypersensitivity following spinal cord injury[31]. However, Rac1 has not been addressed in peripheral sensory neurons in the context of nociceptive modulation in sensory neurons. To address whether Rac1 is upregulated at 24 h after exposure to GMCSF contributes to GMCSF-evoked nociceptive hypersensitivity, we selected dosage of the Rac1-specific inhibitor, NSC23766[32], based on the concentration used by Tan et. al[31] in rats. Extrapolating this concentration to mice and to account for the dilution factor in the CSF, we selected approx. 10 times lesser concentration (1.5 nmoles) in the current study. We divided mice treated with GMCSF over 24 h into 2 groups – one received a single intraplantar injection of NSC23766, a specific Rac1 inhibitor (1.5 nmoles in 10 μl of 50% DMSO) and the other group received a single intraplantar injection of vehicle (10 μl of 50% DMSO) 1 h after the last plantar treatment with GMCSF; GMCSF-mediated mechanical and thermal hypersensitivity was analyzed approx. 3 h and 7 h after inhibitor or vehicle application in both groups. Whereas mice injected with vehicle showed significant mechanical hypersensitivity to 0.16 g of von Frey force as compared to vehicle at 3 h as well as 7 h after the inhibitor application (* P = 0.003 at 3 h and 0.007 at 7 h, One-Way ANOVA with repeated measures followed by Fisher’s LSD Post-hoc analysis, n = 6 mice per group), mice injected with the Rac1 inhibitor did not show any significant deviation from basal response frequencies to 0.16 g force (P > 0.05 as compared to basal; n = 6 mice; † P < 0.01 between groups, One-Way ANOVA with repeated measures followed by Fisher’s LSD Post-hoc analysis, Figure 6A and Additional file6: Table S2). Along the same lines, vehicle-treated mice showed a significant decrease in withdrawal response latencies to plantar heat (i.e. thermal hyperalgesia) as compared to basal values (* P = 0.003 and 0.002 at 4 and 8 h post inhibitor application, respectively, One-way ANOVA followed by Fisher’s LSD Post-hoc analysis, Figure 6B and Additional file7: Table S3). In contrast, Rac1-inhibitor-treated mice did not show thermal hyperalgesia at 4 h after inhibitor application; furthermore, thermal hyperalgesia at both time points tested after inhibitor treatment was significantly reduced as compared to vehicle-treated mice († P = 0.033 at 4 h and P = 0.019 at 8 h between vehicle- and Rac1-inhibitor-treated mice, One-way ANOVA with repeated measures followed by Fisher’s LSD Post-hoc analysis, Figure 6B Additional file7: Table S3). As a control to rule out systemic effects of Rac1 inhibitor, we injected inhibitor or vehicle into the paw contralateral to the paw injected with GMCSF – this treatment failed to block GMCSF-induced mechanical and thermal hyperalgesia in the paw ipsilateral to the GMCSF-injected paw (Additional file8: Figure S5 panels A and C, *P < 0.05 following GMCSF administration, One-Way ANOVA with repeated measures followed by Fisher’s LSD Post-hoc analysis, n = 6 mice). In the same paradigm, mechanical or thermal response frequencies were unaltered as compared to the basal readings in the paw contralateral to the GMCSF-injected paw (Additional file8: Figure S5 panels B and G, P > 0.05). Similarly, injection of the same dosage of Rac1 inhibitor unilaterally into the intraplantar surface in the absence of GMCSF treatment did not affect mechanical and thermal response frequencies when tested up to 7 h post-injection (Additional file9: Figure S6, panels A and B). Thus, these results indicate that locally activated Rac1 specifically contributes to both mechanical and thermal hypersensitivity induced by a prolonged peripheral exposure to GMCSF.

Figure 6

In vivo validation of GM-CSF transcriptional targets in DRG neurons in the context of nociceptive sensitization evoked by long-term exposure to GMCSF. Changes in GMCSF-mediated mechanical hypersensitivity upon inhibition of Rac1 (A), MMP9 (C), Calpain-2 (E) or TNFα (G) as compared to corresponding vehicle-treated mice are shown. Response frequency to the von Frey filament at 0.16 g force is represented on the Y-axis. Changes in GM-CSF-mediated thermal hypersensitivity upon inhibition Rac1 (B), MMP9 (D), Calpain-2 (F) or TNFα (H) as compared to corresponding vehicle-treated mice are shown. Withdrawal latency in seconds to calibrated radiant heat is represented on the Y-axis. * denotes P ≤ 0.05 as compared to basal values, denotes P ≤ 0.05 relative to corresponding vehicle treated group, One-Way ANOVA with repeated measures followed Fisher’s LSD Post-hoc analysis , n = 6 mice per group.

Peripheral MMP9 activation is required for ongoing nociceptive sensitization

In the current study, expression of Mmp9, the gene encoding the extracellular matrix protease MMP9, increased by about 5-fold in DRG following GMCSF stimulation, but not following GCSF stimulation. Given that MMP9 has been shown to participate in inflammatory as well as neuropathic pain in a peripheral as well as a spinal context[3335], we were interested in addressing whether GMCSF-induced upregulation of MMP9 was functionally linked to GMCSF-evoked exaggeration of mechanical and thermal sensitivity. Using the scheme shown in Figure 5C and described in detail above for the Rac1-associated experiments, we administered a single dose of 0.15 pmoles of a potent MMP9 inhibitor in 10 μl of 10% DMSO (vehicle) to the plantar surface 1 h after the last plantar administration of GMCSF. We selected the MMP9 inhibitor dosage based on its high potency (IC50 = 5nM) and its reported intrathecal dosage to attenuate CFA-mediated mechanical allodynia in rats[36]. Upon peripheral MMP9 inhibition, we observed a complete abrogation of GMCSF-evoked mechanical hypersensitivity to 0.16 g of von Frey force as well as thermal hyperalgesia at 3–4 h after MMP9 inhibitor application (Figure 6C and Additional file6: Table S2 *P < 0.05 as compared to basal values; † P < 0.05 as compared between inhibitor and vehicle groups; One-way ANOVA with repeated measures followed by Fisher’s LSD Post-hoc analysis). This effect on mechanical and thermal hyperalgesia was partially or fully lost, respectively, at 7–8 h after inhibitor application, reflecting the duration of action of a single dose of the MMP9-inhibitor at the low dose utilized in this study (Figure 6D and Additional file7: Table S3). Similar to the experiments described above with Rac1 inhibition, we observed that injection of the MMP9 inhibitor in the paw contralateral to the paw injected with GMCSF did not significantly influence GMCSF-mediated mechanical and thermal hyperalgesia in the paw ipsilateral to the GMCSF-injected paw (Additional file8: Figure S5 panels B and D, *P < 0.05 following GM-CSF administration as compared to basal values, One-way ANOVA with repeated measures followed by Fisher’s LSD Post-hoc analysis) nor induced hyperalgesia in the paw contralateral to the GMCSF-injected paw (Additional file8: Figure S5, panels F and H, P > 0.05). Furthermore, injection of the MMP9 inhibitor in the absence of GM-CSF did not affect nociceptive sensitivity (Additional file9: Figure S6 panels C and D, P > 0.05 as compared to basal values, One-way ANOVA with repeated measures followed by Fisher’s LSD Post-hoc analysis). These results indicate that similar to Rac1, peripheral MMP9 activation is important for ongoing nociceptive sensitization that develops upon a prolonged exposure to GMCSF.

Peripheral calpain2 and TNFα are dispensable in the mediation of GMCSF-mediated hyperalgesia

In the current study, Capn2, which encodes the calcium-dependent cysteine proteases calpain 2, emerged as a gene showing a high level of induction, common to both GMCSF and GCSF, in sensory neurons. Drugs which inhibit both Calpain 1 and Calpain 2 have been shown to inhibit mechanical hyperalgesia following inflammation or in response to some mediators[37, 38]. Therefore, we also tested the potential of a peripherally administered Calpain 1/2 inhibitor (Calpain inhibitor III, 10 nmol in 10 μl of 10% DMSO) to modulate GMCSF-mediated mechanical and thermal hyperalgesia using the paradigm described above and inhibitor concentrations that have been shown to be effective in previous studies on inflammatory pain[37, 38]. Our analysis showed, however, that GMCSF-induced hypersensitivity was not significantly different between mice receiving Calpain1/2 inhibitor or vehicle (10 μl of 10% DMSO) in all tests and at all time points tested (Figure 6E and F; Additional file6: Table S2 and Additional file7: Table S3; * P < 0.05 as compared to basal values; One-way ANOVA with repeated measures followed by Fisher’s LSD Post-hoc analysis).

In the current study, tumor necrosis factor alpha (TNFα) was upregulated by 10-fold following GMCSF stimulation but only moderately following GCSF stimulation in microarray or QRTPCR experiments. TNFα inhibition in the spinal cord has been reported to be protective against neuropathic pain[3941] and systemically-applied TNFα inhibitor protects mice from tumor-induced thermal hyperalgesia[42]. Moreover, a specific TNFα inhibitor, namely Etanercept (trade name Enbrel®), has been suggested to be efficacious against autoimmune diseases such as Crohn’s disease[43] and is in clinical practice for the treatment of several peripheral inflammatory diseases such as rheumatoid arthritis[44]. Therefore, we sought out to investigate the potential of peripherally applied TNFα decoy receptor, Etanercept, choosing a dose of 100 pmol, which is slightly higher than the dose reported to reduce neuropathic hypersensitivity[40]. We observed that GMCSF-induced mechanical and thermal hyperalgesia were not significantly different between groups treated with vehicle (BSA) or Etanercept (Figure 6G and H, Additional file6: Table S2 and Additional file7: Table S3; * P < 0.05 following GM-CSF administration as compared to basal values, One-way ANOVA with repeated measures followed by Fisher’s LSD Post-hoc analysis). Thus, albeit TNFα is upregulated in peripheral sensory neurons following GM-CSF administration, it does not appear to directly contribute to GM-CSF-induced nociceptive hypersensitivity.


We and others have previously demonstrated that the hematopoietic growth factors, GMCSF and GCSF, sensitize sensory nerves directly via activation of receptors located on sensory nerves and contribute significantly to cancer-associated pain[8]. Several recent studies have extended these results to acute and chronic pain in the context of post-operative pain, inflammatory pain and osteoarthritic pain[912]. Thus, understanding the genomic program induced in sensory neurons by this set of key cytokines is crucial for gaining insights into the pathophysiological role of G-/GMCSF, as well as in developing therapeutic options. To understand the molecular mediators via which the GMCSF receptor is exerting its protective pronociceptive affects, we performed a genome-wide expression screen in cultured sensory neurons following G-/GMCSF ligand treatments. Using in-depth systems level in silico analysis, and in vivo pharmacology to functionally test selected candidate genes, we present G-/GMCSF-mediated transcriptome change and its importance in modulating GMCSF-mediated hyperalgesia in the sensory neurons.

One of the most interesting observations of the present study was that many chemokines such as Ccl3, Ccl5, insulin-like growth factor 1 (Igf1), Vascular endothelial growth factor A (Vegfa), transcript variant 2 and TNFα, which are known to play a significant role in pain modulation[20, 26, 45, 46] are regulated in a G-/GMCSF dependent manner in sensory neurons. It was previously reported that the lack of CCL5 results in decreased hypersensitivity in partial sciatic nerve ligation model of neuropathic pain. Further, infiltration of pro-inflammatory cytokines such as tumor necrosis factor-α, [is this supposed to be alpha?] interleukin [IL]-1b, IL-6, and interferon-c was significantly reduced in the damaged nerves while that of anti-inflammatory cytokines such as IL-4 and IL-10 was significantly increased in the injured nerves following the global CCL5 loss in mice[26]. TNFα was shown to have a significant contribution in the development of pain sensitivity following peripheral nerve injury.

Rac1 is a small GTPase ubiquitously expressed in neurons and other cell types. Spinal cord neurons, astrocytes and oligodendrocytes express high levels of Rac1 mRNA and its expression is elevated following spinal cord injury (SCI) in rats[47]. Intrathecally applied specific Rac1 inhibitor ameliorates SCI-induced changes in spine morphology, attenuates injury-induced hyper excitability of wide-dynamic range neurons, and increases SCI-mediated pain thresholds[31]. Rac1 activity regulates dendritic spine morphology in hippocampal neurons through actin cytoskeleton reorganization and promotes clustering of glutamate AMPA receptors in spines[4850]. In the current study, Rac1 expression increased by about 4–5 fold in sensory neurons following exposure to GMCSF or GCSF stimulation, suggesting that this may work as a shared mechanism between GMCSF and GCSF stimuli to modulate functional as well as structural changes in sensory neurons. Indeed, we have previously shown that knocking down GMCSF receptor-alpha in sensory neurons of the DRG attenuates tumor-induced structural remodeling and sprouting of peripheral sensory terminals[8]. Further, inhibiting Rac1 activity in the peripheral terminals resulted in complete protection from GMCSF-induced mechanical hypersensitivity and partial protection from thermal hypersensitivity. These effects of peripherally applied Rac1 inhibitor were devoid of a systemic component and did not come about via an alteration of basal mechanical and thermal sensitivity. These results underline the importance of Rac1 in the periphery in the modulation of nociceptive stimuli, in the context of both mechanical and thermal modalities.

MMP9 is under the transcriptional control of Ets1 and belongs to a family of extracellular proteases called matrix metalloproteinases (MMPs)[51]. MMPs play a critical role in neuroinflammation through the cleavage of the extracellular matrix (ECM) proteins, cytokines, and chemokines[5254]. Consistent with previously reported changes in MMP9 expression in nervous tissue following damage to peripheral nerves[33, 35, 52, 55, 56], we observed a 5-fold increase in MMP9 expression following GMCSF stimulus in sensory neurons. Intrathecally administered MMP9 inhibitors have been shown to block mechanical allodynia associated with peripheral inflammation[36] as well as nerve injury[35]. Functional experiments performed in the current study show that MMP9 inhibition in the periphery leads to significant protection from GMCSF-mediated mechanical hypersensitivity in the absence of a systemic influence or without affecting basal nociception. Interestingly, these findings thus indicate that MMP9 exerts pronociceptive effects not only in the central terminals but also at the peripheral nerve endings. Our findings further suggest that MMP9 may be linked to the role of GMCSF in tumor-associated pain, since both inflammatory and neuropathic mechanisms contribute to cancer pain. However, peripheral blockade of MMP9 did not affect GMCSF-mediated thermal hyperalgesia, which is consistent with a lack of its role in inflammatory thermal hyperalgesia.

Calpains are a family of ubiquitously expressed calcium-dependent cysteine proteases and have been implicated in several cellular processes such as proliferation, apoptosis and differentiation[5759], pathological conditions including neuronal plasticity[60], neuronal cell death[61]. Proteolytic targets of calpain include cytoskeletal proteins such as neurofilament, spectrin, and membrane proteins such as growth factor receptors[62, 63], cytokines[64, 65] and transcription factors[66]. In the current study, long-term GMCSF or GCSF exposure led to robust transcriptional increase in the expression of calpain 2 and 7, but not other calpain family members, in sensory neurons. However, inhibiting the calpain 2 protease by using a specific calpain I/II inhibitor in the periphery did not have any impact on modulating GMCSF-induced mechanical and thermal hyperalgesia. This lack of phenotypic changes is not due to a lack of efficacy of the inhibitor since intrathecal application of comparable concentrations has been shown to significantly reduce zymosan-mediated thermal hypersensitivity in rats[37] or lysophosphatidic acid (LPA)-mediated thermal hyperalgesia in mice[38]. Thus, our results indicate that calpain is dispensable for GMCSF-induced nociceptive sensitization in the periphery; however, they do not rule out that GMCSF-mediated induction of calpain expression may be modulating other functions and processes in the DRG which were not studied here.

TNFα is a proinflammatory chemokine which was previously studied intensively in the context of nociceptive modulation[67, 68]. Importantly, intraperitoneal application of TNFα decoy receptor etanercept relieved mice from tumor-mediated hyperalgesia[42]. However, it was also reported that intrathecally and intraperitoneally applied etanercept protects the mice from diabetic neuropathy-induced mechanical hyperalgesia but intraplantar application of the same inhibitor showed inefficient in protecting the mice from diabetic neuropathy-induced hypersensitivity[40]. Consistent with this observation, we observed that peripheral inhibition of TNFα is not sufficient to abrogate the nociceptive stimulus-mediated sensitivity. Taken together, these observations suggest that TNFα is recruited downstream of other tumor-associated mediators.

Other than the genes directly regulated by G-/GMCSF transcriptionally, our systems level network analysis revealed several other pathways which could be directly connected to the G-/GMCSF-mediated genes, such as IKK-NF-κB pathway. This classical, canonical pathway involves TNFα or IL-1β stimuli via respective receptors leading to the activation of an inhibitor of the NF-κB (I-κB) kinase complex, consisting of the regulatory subunit I-κB kinases (IKKs). This pathway is crucial for the activation of innate immunity and inflammation[69]. Functionally, increased NF-κB levels in DRG neurons following sciatic nerve crush[70] or upregulation of NF- κB transactivation following sciatic nerve transection[71] have been reported. Pharmacological intervention at several nodes of the NF-κB pathway has been shown to modulate nociceptive responses[69]. Inactivation of NF-κB specifically in primary sensory neurons of DRG via Cre-LoxP-mediated deletion of IKKB inhibitor kappa B kinase beta (IKKβ) has been shown to have a protective effect on nerve injury-mediated hyperalgesia[72]. Interestingly, the expression of TNFα and IL-1β robustly increased following G-/GMCSF stimuli in the current study, implying a clearer activation of NF-κB pathway in sensory neurons. Another molecular pathway which emerged to be very closely connected to G-/GMCSF-mediated transcriptome in sensory neurons is caspase signaling. Caspases are a family of proteases and play an important role in mediating programmed cell death following different noxious stimuli[73, 74]. Peripheral nerve injury promotes neuronal cell death in the spinal dorsal horn[7577] and arresting this nerve injury-induced neuronal loss in the spinal dorsal horn by blocking caspase activity reduces neuropathy-induced pain sensitivity[78]. Caspase signaling pathways differentially contribute to neuropathy-induced and TNF-mediated pain behaviors[79]. Our results on network analysis indicate that GMCSF signaling may be interlinked with TNF-alpha and caspase signaling in DRG neurons.

Thus, the striking changes we report in the transcription of many pain-related ion channels, chemokines, growth factors and proteases among several other classes of genes in DRG neurons following prolonged exposure to G-/GMCSF imply that G-/GMCSF signaling is a trigger point for activation of multiple pain modulatory pathways and that blocking the G-/GMCSF signaling may be very effective in alleviating a broad set of pain disorders.


In summary, the current study demonstrates genome-wide transcriptome changes following chronic G-/GMCSF stimulus in the sensory neurons. Utilizing state-of-the art in silico systems level analysis, this study not only reveals that several key pain-related genes to be transcriptional targets of G-/GMCSF signaling, but also provides novel insights into network interactions with several other novel candidate genes. Using in vivo pharmacology, we provide the importance of peripheral MMP9 and Rac1 signaling in inhibiting GMCSF-mediated mechanical and thermal hypersensitivity. Thus, with integrative approach of genomics, bioinformatics, in vivo pharmacology and behavioral analyses, this study advances the understanding of nociceptive mechanisms in sensory neurons and provides a basis for further pursuing G-/GMCSF signaling in therapeutic treatment of pain disorders.



Dorsal root ganglion


Granulocyte macrophage colony stimulating factor


Granulocyte colony stimulating factor


Phosphate buffered saline.


  1. 1.

    Jimenez-Andrade JM, Mantyh WG, Bloom AP, Ferng AS, Geffre CP, Mantyh PW: Bone cancer pain. Ann N Y Acad Sci 2010, 1198: 173–181. 10.1111/j.1749-6632.2009.05429.x

  2. 2.

    Davar G: Endothelin-1 and metastatic cancer pain. Pain Med 2001, 2: 24–27. 10.1046/j.1526-4637.2001.002001024.x

  3. 3.

    Grant K, Loizidou M, Taylor I: Endothelin-1: a multifunctional molecule in cancer. Br J Cancer 2003, 88: 163–166. 10.1038/sj.bjc.6700750

  4. 4.

    Stosser S, Agarwal N, Tappe-Theodor A, Yanagisawa M, Kuner R: Dissecting the functional significance of endothelin A receptors in peripheral nociceptors in vivo via conditional gene deletion. Pain 2009, 148: 206–214.

  5. 5.

    Mantyh PW: Cancer pain and its impact on diagnosis, survival and quality of life. Nat Rev Neurosci 2006, 7: 797–809.

  6. 6.

    Stosser S, Schweizerhof M, Kuner R: Hematopoietic colony-stimulating factors: new players in tumor-nerve interactions. J Mol Med (Berl) 2011, 89: 321–329. 10.1007/s00109-010-0697-z

  7. 7.

    Hamilton JA: Colony-stimulating factors in inflammation and autoimmunity. Nat Rev Immunol 2008, 8: 533–544. 10.1038/nri2356

  8. 8.

    Schweizerhof M, Stosser S, Kurejova M, Njoo C, Gangadharan V, Agarwal N, Schmelz M, Bali KK, Michalski CW, Brugger S, et al.: Hematopoietic colony-stimulating factors mediate tumor-nerve interactions and bone cancer pain. Nat Med 2009, 15: 802–807. 10.1038/nm.1976

  9. 9.

    Cook AD, Pobjoy J, Sarros S, Steidl S, Durr M, Lacey DC, Hamilton JA: Granulocyte-macrophage colony-stimulating factor is a key mediator in inflammatory and arthritic pain. Ann Rheum Dis 2013, 72: 265–270. 10.1136/annrheumdis-2012-201703

  10. 10.

    Cook AD, Pobjoy J, Steidl S, Durr M, Braine EL, Turner AL, Lacey DC, Hamilton JA: Granulocyte-macrophage colony-stimulating factor is a key mediator in experimental osteoarthritis pain and disease development. Arthritis Res Ther 2012, 14: R199. 10.1186/ar4037

  11. 11.

    Leah E: Experimental arthritis: GM-CSF mediates pain and disease in a mouse model of osteoarthritis. Nat Rev Rheumatol 2012, 8: 634.

  12. 12.

    Onuora S: Pain: granulocyte-macrophage colony-stimulating factor required for inflammatory and arthritic pain. Nat Rev Rheumatol 2012, 8: 499. 10.1038/nrrheum.2012.139

  13. 13.

    Rawlings JS, Rosler KM, Harrison DA: The JAK/STAT signaling pathway. J Cell Sci 2004, 117: 1281–1283. 10.1242/jcs.00963

  14. 14.

    Boyd TD, Bennett SP, Mori T, Governatori N, Runfeldt M, Norden M, Padmanabhan J, Neame P, Wefes I, Sanchez-Ramos J, et al.: GM-CSF upregulated in rheumatoid arthritis reverses cognitive impairment and amyloidosis in Alzheimer mice. J Alzheimers Dis 2010, 21: 507–518.

  15. 15.

    Cornish AL, Campbell IK, McKenzie BS, Chatfield S, Wicks IP: G-CSF and GM-CSF as therapeutic targets in rheumatoid arthritis. Nat Rev Rheumatol 2009, 5: 554–559. 10.1038/nrrheum.2009.178

  16. 16.

    Chudin E, Kruglyak S, Baker SC, Oeser S, Barker D, McDaniel TK: A model of technical variation of microarray signals. J Comput Biol 2006, 13: 996–1003. 10.1089/cmb.2006.13.996

  17. 17.

    Fu WJ, Hu J, Spencer T, Carroll R, Wu G: Statistical models in assessing fold change of gene expression in real-time RT-PCR experiments. Comput Biol Chem 2006, 30: 21–26. 10.1016/j.compbiolchem.2005.10.005

  18. 18.

    Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 2003, 19: 185–193. 10.1093/bioinformatics/19.2.185

  19. 19.

    Simonetti M, Hagenston AM, Vardeh D, Freitag HE, Mauceri D, Lu J, Satagopam VP, Schneider R, Costigan M, Bading H, Kuner R: Nuclear calcium signaling in spinal neurons drives a genomic program required for persistent inflammatory pain. Neuron 2013, 77: 43–57. 10.1016/j.neuron.2012.10.037

  20. 20.

    Quick ML, Mukherjee S, Rudick CN, Done JD, Schaeffer AJ, Thumbikat P: CCL2 and CCL3 are essential mediators of pelvic pain in experimental autoimmune prostatitis. Am J Physiol Regul Integr Comp Physiol 2012, 303: R580-R589. 10.1152/ajpregu.00240.2012

  21. 21.

    Zhang N, Rogers TJ, Caterina M, Oppenheim JJ: Proinflammatory chemokines, such as C-C chemokine ligand 3, desensitize mu-opioid receptors on dorsal root ganglia neurons. J Immunol 2004, 173: 594–599.

  22. 22.

    Julius D, Basbaum AI: Molecular mechanisms of nociception. Nature 2001, 413: 203–210. 10.1038/35093019

  23. 23.

    Tulleuda A, Cokic B, Callejo G, Saiani B, Serra J, Gasull X: TRESK channel contribution to nociceptive sensory neurons excitability: modulation by nerve injury. Mol Pain 2011, 7: 30. 10.1186/1744-8069-7-30

  24. 24.

    Dong H, O’Brien RJ, Fung ET, Lanahan AA, Worley PF, Huganir RL: GRIP: a synaptic PDZ domain-containing protein that interacts with AMPA receptors. Nature 1997, 386: 279–284. 10.1038/386279a0

  25. 25.

    Tao YX: Dorsal horn alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptor trafficking in inflammatory pain. Anesthesiology 2010, 112: 1259–1265. 10.1097/ALN.0b013e3181d3e1ed

  26. 26.

    Liou JT, Yuan HB, Mao CC, Lai YS, Day YJ: Absence of C-C motif chemokine ligand 5 in mice leads to decreased local macrophage recruitment and behavioral hypersensitivity in a murine neuropathic pain model. Pain 2012, 153: 1283–1291. 10.1016/j.pain.2012.03.008

  27. 27.

    Catrina AI, Lampa J, Ernestam S, af Klint E, Bratt J, Klareskog L, Ulfgren AK: Anti-tumour necrosis factor (TNF)-alpha therapy (etanercept) down-regulates serum matrix metalloproteinase (MMP)-3 and MMP-1 in rheumatoid arthritis. Rheumatology (Oxford) 2002, 41: 484–489. 10.1093/rheumatology/41.5.484

  28. 28.

    Green MJ, Gough AK, Devlin J, Smith J, Astin P, Taylor D, Emery P: Serum MMP-3 and MMP-1 and progression of joint damage in early rheumatoid arthritis. Rheumatology (Oxford) 2003, 42: 83–88. 10.1093/rheumatology/keg037

  29. 29.

    Ekins S, Nikolsky Y, Bugrim A, Kirillov E, Nikolskaya T: Pathway mapping tools for analysis of high content data. Methods Mol Biol 2007, 356: 319–350.

  30. 30.

    Lee KM, Kang BS, Lee HL, Son SJ, Hwang SH, Kim DS, Park JS, Cho HJ: Spinal NF-kB activation induces COX-2 upregulation and contributes to inflammatory pain hypersensitivity. Eur J Neurosci 2004, 19: 3375–3381. 10.1111/j.0953-816X.2004.03441.x

  31. 31.

    Tan AM, Stamboulian S, Chang YW, Zhao P, Hains AB, Waxman SG, Hains BC: Neuropathic pain memory is maintained by Rac1-regulated dendritic spine remodeling after spinal cord injury. J Neurosci 2008, 28: 13173–13183. 10.1523/JNEUROSCI.3142-08.2008

  32. 32.

    Gao Y, Dickerson JB, Guo F, Zheng J, Zheng Y: Rational design and characterization of a Rac GTPase-specific small molecule inhibitor. Proc Natl Acad Sci USA 2004, 101: 7618–7623. 10.1073/pnas.0307512101

  33. 33.

    Chattopadhyay S, Myers RR, Janes J, Shubayev V: Cytokine regulation of MMP-9 in peripheral glia: implications for pathological processes and pain in injured nerve. Brain Behav Immun 2007, 21: 561–568. 10.1016/j.bbi.2006.10.015

  34. 34.

    Ji RR, Xu ZZ, Wang X, Lo EH: Matrix metalloprotease regulation of neuropathic pain. Trends Pharmacol Sci 2009, 30: 336–340. 10.1016/

  35. 35.

    Kawasaki Y, Xu ZZ, Wang X, Park JY, Zhuang ZY, Tan PH, Gao YJ, Roy K, Corfas G, Lo EH, Ji RR: Distinct roles of matrix metalloproteases in the early- and late-phase development of neuropathic pain. Nat Med 2008, 14: 331–336. 10.1038/nm1723

  36. 36.

    Kular L, Rivat C, Lelongt B, Calmel C, Laurent M, Pohl M, Kitabgi P, Melik-Parsadaniantz S, Martinerie C: NOV/CCN3 attenuates inflammatory pain through regulation of matrix metalloproteinases-2 and −9. J Neuroinflammation 2012, 9: 36. 10.1186/1742-2094-9-36

  37. 37.

    Kunz S, Niederberger E, Ehnert C, Coste O, Pfenninger A, Kruip J, Wendrich TM, Schmidtko A, Tegeder I, Geisslinger G: The calpain inhibitor MDL 28170 prevents inflammation-induced neurofilament light chain breakdown in the spinal cord and reduces thermal hyperalgesia. Pain 2004, 110: 409–418. 10.1016/j.pain.2004.04.031

  38. 38.

    Xie W, Uchida H, Nagai J, Ueda M, Chun J, Ueda H: Calpain-mediated down-regulation of myelin-associated glycoprotein in lysophosphatidic acid-induced neuropathic pain. J Neurochem 2010, 113: 1002–1011. 10.1111/j.1471-4159.2010.06664.x

  39. 39.

    Abram SE, Yi J, Fuchs A, Hogan QH: Permeability of injured and intact peripheral nerves and dorsal root ganglia. Anesthesiology 2006, 105: 146–153. 10.1097/00000542-200607000-00024

  40. 40.

    Dogrul A, Gul H, Yesilyurt O, Ulas UH, Yildiz O: Systemic and spinal administration of etanercept, a tumor necrosis factor alpha inhibitor, blocks tactile allodynia in diabetic mice. Acta Diabetol 2011, 48: 135–142. 10.1007/s00592-010-0237-x

  41. 41.

    Shubayev VI, Myers RR: Anterograde TNF alpha transport from rat dorsal root ganglion to spinal cord and injured sciatic nerve. Neurosci Lett 2002, 320: 99–101. 10.1016/S0304-3940(02)00010-1

  42. 42.

    Constantin CE, Mair N, Sailer CA, Andratsch M, Xu ZZ, Blumer MJ, Scherbakov N, Davis JB, Bluethmann H, Ji RR, Kress M: Endogenous tumor necrosis factor alpha (TNFalpha) requires TNF receptor type 2 to generate heat hyperalgesia in a mouse cancer model. J Neurosci 2008, 28: 5072–5081. 10.1523/JNEUROSCI.4476-07.2008

  43. 43.

    Sandborn WJ, Hanauer SB, Katz S, Safdi M, Wolf DG, Baerg RD, Tremaine WJ, Johnson T, Diehl NN, Zinsmeister AR: Etanercept for active Crohn’s disease: a randomized, double-blind, placebo-controlled trial. Gastroenterology 2001, 121: 1088–1094. 10.1053/gast.2001.28674

  44. 44.

    Lovell DJ, Giannini EH, Reiff A, Cawkwell GD, Silverman ED, Nocton JJ, Stein LD, Gedalia A, Ilowite NT, Wallace CA, et al.: Etanercept in children with polyarticular juvenile rheumatoid arthritis. Pediatric Rheumatology Collaborative Study Group. N Engl J Med 2000, 342: 763–769. 10.1056/NEJM200003163421103

  45. 45.

    Cuellar JM, Montesano PX, Carstens E: Role of TNF-alpha in sensitization of nociceptive dorsal horn neurons induced by application of nucleus pulposus to L5 dorsal root ganglion in rats. Pain 2004, 110: 578–587. 10.1016/j.pain.2004.03.029

  46. 46.

    Czeschik JC, Hagenacker T, Schafers M, Busselberg D: TNF-alpha differentially modulates ion channels of nociceptive neurons. Neurosci Lett 2008, 434: 293–298. 10.1016/j.neulet.2008.01.070

  47. 47.

    Erschbamer MK, Hofstetter CP, Olson L: RhoA, RhoB, RhoC, Rac1, Cdc42, and Tc10 mRNA levels in spinal cord, sensory ganglia, and corticospinal tract neurons and long-lasting specific changes following spinal cord injury. J Comp Neurol 2005, 484: 224–233. 10.1002/cne.20471

  48. 48.

    Nakayama AY, Harms MB, Luo L: Small GTPases Rac and Rho in the maintenance of dendritic spines and branches in hippocampal pyramidal neurons. J Neurosci 2000, 20: 5329–5338.

  49. 49.

    Tashiro A, Minden A, Yuste R: Regulation of dendritic spine morphology by the rho family of small GTPases: antagonistic roles of Rac and Rho. Cereb Cortex 2000, 10: 927–938. 10.1093/cercor/10.10.927

  50. 50.

    Tashiro A, Yuste R: Regulation of dendritic spine motility and stability by Rac1 and Rho kinase: evidence for two forms of spine motility. Mol Cell Neurosci 2004, 26: 429–440. 10.1016/j.mcn.2004.04.001

  51. 51.

    Woessner JF Jr: The family of matrix metalloproteinases. Ann N Y Acad Sci 1994, 732: 11–21. 10.1111/j.1749-6632.1994.tb24720.x

  52. 52.

    Rosenberg GA: Matrix metalloproteinases in neuroinflammation. Glia 2002, 39: 279–291. 10.1002/glia.10108

  53. 53.

    Parks WC, Wilson CL, Lopez-Boado YS: Matrix metalloproteinases as modulators of inflammation and innate immunity. Nat Rev Immunol 2004, 4: 617–629. 10.1038/nri1418

  54. 54.

    Yong VW: Metalloproteinases: mediators of pathology and regeneration in the CNS. Nat Rev Neurosci 2005, 6: 931–944. 10.1038/nrn1807

  55. 55.

    Demestre M, Wells GM, Miller KM, Smith KJ, Hughes RA, Gearing AJ, Gregson NA: Characterisation of matrix metalloproteinases and the effects of a broad-spectrum inhibitor (BB-1101) in peripheral nerve regeneration. Neuroscience 2004, 124: 767–779. 10.1016/j.neuroscience.2003.12.037

  56. 56.

    Leppert D, Ford J, Stabler G, Grygar C, Lienert C, Huber S, Miller KM, Hauser SL, Kappos L: Matrix metalloproteinase-9 (gelatinase B) is selectively elevated in CSF during relapses and stable phases of multiple sclerosis. Brain 1998, 121(Pt 12):2327–2334.

  57. 57.

    Li G, Iyengar R: Calpain as an effector of the Gq signaling pathway for inhibition of Wnt/beta -catenin-regulated cell proliferation. Proc Natl Acad Sci USA 2002, 99: 13254–13259. 10.1073/pnas.202355799

  58. 58.

    Neumar RW, Xu YA, Gada H, Guttmann RP, Siman R: Cross-talk between calpain and caspase proteolytic systems during neuronal apoptosis. J Biol Chem 2003, 278: 14162–14167. 10.1074/jbc.M212255200

  59. 59.

    Saito Y, Saido TC, Sano K, Kawashima S: The calpain-calpastatin system is regulated differently during human neuroblastoma cell differentiation to Schwannian and neuronal cells. FEBS Lett 1994, 353: 327–331. 10.1016/0014-5793(94)01075-7

  60. 60.

    Vanderklish PW, Krushel LA, Holst BH, Gally JA, Crossin KL, Edelman GM: Marking synaptic activity in dendritic spines with a calpain substrate exhibiting fluorescence resonance energy transfer. Proc Natl Acad Sci USA 2000, 97: 2253–2258. 10.1073/pnas.040565597

  61. 61.

    Chera B, Schaecher KE, Rocchini A, Imam SZ, Ray SK, Ali SF, Banik NL: Calpain upregulation and neuron death in spinal cord of MPTP-induced parkinsonism in mice. Ann N Y Acad Sci 2002, 965: 274–280.

  62. 62.

    Baliova M, Betz H, Jursky F: Calpain-mediated proteolytic cleavage of the neuronal glycine transporter, GlyT2. J Neurochem 2004, 88: 227–232.

  63. 63.

    Chan SL, Mattson MP: Caspase and calpain substrates: roles in synaptic plasticity and cell death. J Neurosci Res 1999, 58: 167–190. 10.1002/(SICI)1097-4547(19991001)58:1<167::AID-JNR16>3.0.CO;2-K

  64. 64.

    Banik NL, Matzelle D, Gantt-Wilford G, Hogan EL: Role of calpain and its inhibitors in tissue degeneration and neuroprotection in spinal cord injury. Ann N Y Acad Sci 1997, 825: 120–127. 10.1111/j.1749-6632.1997.tb48421.x

  65. 65.

    Kampfl A, Posmantur RM, Zhao X, Schmutzhard E, Clifton GL, Hayes RL: Mechanisms of calpain proteolysis following traumatic brain injury: implications for pathology and therapy: implications for pathology and therapy: a review and update. J Neurotrauma 1997, 14: 121–134. 10.1089/neu.1997.14.121

  66. 66.

    Shields DC, Banik NL: Putative role of calpain in the pathophysiology of experimental optic neuritis. Exp Eye Res 1998, 67: 403–410. 10.1006/exer.1998.0537

  67. 67.

    Leung L, Cahill CM: TNF-alpha and neuropathic pain–a review. J Neuroinflammation 2010, 7: 27. 10.1186/1742-2094-7-27

  68. 68.

    Zheng W, Ouyang H, Zheng X, Liu S, Mata M, Fink DJ, Hao S: Glial TNFalpha in the spinal cord regulates neuropathic pain induced by HIV gp120 application in rats. Mol Pain 2011, 7: 40. 10.1186/1744-8069-7-40

  69. 69.

    Niederberger E, Geisslinger G: The IKK-NF-kappaB pathway: a source for novel molecular drug targets in pain therapy? FASEB J 2008, 22: 3432–3442. 10.1096/fj.08-109355

  70. 70.

    Fernyhough P, Smith DR, Schapansky J, Van Der Ploeg R, Gardiner NJ, Tweed CW, Kontos A, Freeman L, Purves-Tyson TD, Glazner GW: Activation of nuclear factor-kappaB via endogenous tumor necrosis factor alpha regulates survival of axotomized adult sensory neurons. J Neurosci 2005, 25: 1682–1690. 10.1523/JNEUROSCI.3127-04.2005

  71. 71.

    Pollock G, Pennypacker KR, Memet S, Israel A, Saporta S: Activation of NF-kappaB in the mouse spinal cord following sciatic nerve transection. Exp Brain Res 2005, 165: 470–477. 10.1007/s00221-005-2318-6

  72. 72.

    Kanngiesser M, Haussler A, Myrczek T, Kusener N, Lim HY, Geisslinger G, Niederberger E, Tegeder I: Inhibitor kappa B kinase beta dependent cytokine upregulation in nociceptive neurons contributes to nociceptive hypersensitivity after sciatic nerve injury. J Pain 2012, 13: 485–497. 10.1016/j.jpain.2012.02.010

  73. 73.

    Degterev A, Boyce M, Yuan J: A decade of caspases. Oncogene 2003, 22: 8543–8567. 10.1038/sj.onc.1207107

  74. 74.

    Ferri KF, Kroemer G: Organelle-specific initiation of cell death pathways. Nat Cell Biol 2001, 3: E255-E263. 10.1038/ncb1101-e255

  75. 75.

    Moore JD, Rothwell NJ, Gibson RM: Involvement of caspases and calpains in cerebrocortical neuronal cell death is stimulus-dependent. Br J Pharmacol 2002, 135: 1069–1077. 10.1038/sj.bjp.0704538

  76. 76.

    Sugimoto T, Bennett GJ, Kajander KC: Transsynaptic degeneration in the superficial dorsal horn after sciatic nerve injury: effects of a chronic constriction injury, transection, and strychnine. Pain 1990, 42: 205–213. 10.1016/0304-3959(90)91164-E

  77. 77.

    Whiteside GT, Munglani R: Cell death in the superficial dorsal horn in a model of neuropathic pain. J Neurosci Res 2001, 64: 168–173. 10.1002/jnr.1062

  78. 78.

    Scholz J, Broom DC, Youn DH, Mills CD, Kohno T, Suter MR, Moore KA, Decosterd I, Coggeshall RE, Woolf CJ: Blocking caspase activity prevents transsynaptic neuronal apoptosis and the loss of inhibition in lamina II of the dorsal horn after peripheral nerve injury. J Neurosci 2005, 25: 7317–7323. 10.1523/JNEUROSCI.1526-05.2005

  79. 79.

    Joseph EK, Levine JD: Caspase signalling in neuropathic and inflammatory pain in the rat. Eur J Neurosci 2004, 20: 2896–2902. 10.1111/j.1460-9568.2004.03750.x

Download references


The authors are grateful to Rose LeFaucheur for secretarial work and to Dunja Baumgartl-Ahlert for technical assistance. This work was supported by grants from the Association for International Cancer Research and the Sander Foundation to RK. RK is a principal investigator in the Excellence Cluster ‘CellNetworks’ of Heidelberg University.

Author information

Correspondence to Kiran Kumar Bali or Rohini Kuner.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

KB performed the major portion of the experiments, analyzed data and wrote the manuscript. VV and PG contributed to experiments and performed data analysis. VPS and RS contributed to the bioinformatics part of the study. RK designed and supervised the study and wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Table S1: Significantly regulated genes following GMCSF or GCSF chronic stimulus in the sensory neurons with the 2 fold cut-off and pBH < 0.05 (two-tailed t-test assuming equal variance, P with Benjamini and Hochberg False Discovery Rate correction <0.05). (XLSX 442 KB)

Additional file 2: Figure S1: Symbols used to represent different functional classes of protein to represent direct-interactions networks in Figures 3 and 4 and network interactions in Additional file 3: Figure S2 and Additional file 4: Figure S3. (PDF 293 KB)

Additional file 3: Figure S2: Direct-interactions-network analysis on the gene pool commonly up- and down-regulated following GMCSF or GCSF exposure in the sensory neurons. Gene pool with fold-change between +4 and −4 as compared to control-treated sensory neurons and P-BH < 0.05 (t-test, P with Benjamini and Hochberg False Discovery Rate < 0.05). Genes upregulated and downregulated in GMCSF-dependent manner are marked with red and blue circles respectively. Please see Additional file 2: Figure S1 for information on legends. (PDF 176 KB)


Additional file 4: Figure S3: Gene-ontology enrichment analysis on gene pools induced by GMCSF (A) and GCSF (B). Analysis was performed using bioCompendium online repository ( using the gene pool with fold-change between +2 and −2, as compared to control-treated sensory neurons and P-BH < 0.05 (t-test, P with Benjamini and Hochberg False Discovery Rate <0.05). Percent of G-/GM-CSF induced genes over total input number of genes significantly enriched are represented in both panels (P-BH < 0.05, as compared to mouse genome, hypergeometric statistical method). (PDF 244 KB)

Additional file 5: Figure S4: Top scored network obtained from the gene pool regulated by GM-CSF stimulus in sensory neurons, using the same set of genes explained in the Figure 1-A. Genes upregulated and downregulated in GMCSF-dependent manner are shown with red and blue circles respectively. Thick cyan lines indicate the fragments of canonical pathways. Please refer Additional file 2: Figure S1 for the symbols representing different generic classes of proteins. (PDF 231 KB)

Additional file 6: Table S2: Impact of inhibition of Rac1, MMP9, Calpain2 or TNFα on GMCSF-mediated mechanical hypersensitivity. Response frequencies to calibrated von Frey filaments of different strength at paws ipsilateral and contralateral to intraplantar GMCSF administration are shown as compared to corresponding vehicle-treated mice. * denotes P ≤ 0.05 as compared to basal values, denotes P ≤ 0.05 relative to corresponding vehicle-treated group, One-Way ANOVA with repeated measures followed Fisher’s LSD Post-hoc analysis , n = 6 mice per group. (DOCX 14 KB)

Additional file 7: Table S3: Impact of Inhibition of Rac1, MMP9, Calpain2 or TNFα on GMCSF-mediated thermal hypersensitivity. Withdrawal latency in seconds to calibrated radiant heat from the paws ipsilateral and contralateral to intraplantar GMCSF administration are shown as compared to corresponding vehicle-treated mice. * denotes P ≤ 0.05 as compared to basal values, denotes P ≤ 0.05 relative to corresponding vehicle treated group, One-Way ANOVA with repeated measures followed Fisher’s LSD Post-hoc analysis, n = 6 mice per group. (DOCX 12 KB)

Additional file 8: Figure S5: Systemic effects of intraplantar application of inhibitors of specific pathways on GMCSF-induced nociceptive sensitization. Changes in the GMCSF-mediated mechanical hypersensitivity in the paw ipsilateral to the GMCSF-injected paw following Rac1 (A) or MMP9 (B) inhibitors application in the paw contralaeral to GMCSF-injected paw or changes in the GM-CSF-mediated mechanical hypersensitivity in the paw contralateral to the GMCSF-injected paw following Rac1 (E) or MMP9 (F) inhibitors application in the paw contralaeral to GMCSF-injected paw as compared to corresponding vehicle-treated mice are shown. Response frequency to the von Frey filament at 0.16 g force is represented on the Y-axis. Changes in the GMCSF-mediated thermal hypersensitivity in the paw ipsilateral to the GMCSF-injected paw following Rac1 (C) or MMP9 (D) inhibitors or changes in the GMCSF-mediated thermal hypersensitivity in the paw contralateral to the GMCSF-injected paw following Rac1 (G) or MMP9 (H) inhibitors as compared to corresponding vehicle-treated group of mice. Withdrawal latency in seconds to calibrated radiant heat is represented. * denotes P ≤ 0.05 as compared to basal values, One-Way ANOVA with repeated measures followed Fisher’s LSD Post-hoc analysis, n = 6 mice per group. (PDF 105 KB)

Additional file 9: Figure S6: Effect of blockade of Rac1 or MMP9 on basal mechanical and thermal sensitivity. Changes in the perception of mechanical sensitivity in response to von Frey filaments of increasing strength in the paw ipsilateral to the Rac1 (A) or MMP9 (C) inhibitor application is compared to the response frequency before inhibitor application (basal). Percentage Response frequency to the von Frey filament at 0.16 g force is represented on the Y-axis. Changes in the response to thermal stimuli following Rac1 (B) or MMP9 (D) inhibitors as compared to basal reading. Withdrawal latency in seconds to calibrated radiant heat is represented. * denotes P ≤ 0.05 as compared to basal values, One-Way ANOVA with repeated measures followed Fisher’s LSD Post-hoc analysis , n = 6 mice per group. (PDF 100 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Bali, K.K., Venkataramani, V., Satagopam, V.P. et al. Transcriptional mechanisms underlying sensitization of peripheral sensory neurons by Granulocyte-/Granulocyte-macrophage colony stimulating factors. Mol Pain 9, 48 (2013) doi:10.1186/1744-8069-9-48

Download citation


  • Dorsal Root Ganglion
  • Sensory Neuron
  • Granulocyte Colony Stimulate Factor
  • Dorsal Root Ganglion Neuron
  • Granulocyte Macrophage Colony Stimulate Factor