Skip to main content

Expression profiling of genes modulated by minocycline in a rat model of neuropathic pain

Abstract

Background

The molecular mechanisms underlying neuropathic pain are constantly being studied to create new opportunities to prevent or alleviate neuropathic pain. The aim of our study was to determine the gene expression changes induced by sciatic nerve chronic constriction injury (CCI) that are modulated by minocycline, which can effectively diminish neuropathic pain in animal studies. The genes associated with minocycline efficacy in neuropathic pain should provide insight into the etiology of neuropathic pain and identify novel therapeutic targets.

Results

We screened the ipsilateral dorsal part of the lumbar spinal cord of the rat CCI model for differentially expressed genes. Out of 22,500 studied transcripts, the abundance levels of 93 transcripts were altered following sciatic nerve ligation. Percentage analysis revealed that 54 transcripts were not affected by the repeated administration of minocycline (30 mg/kg, i.p.), but the levels of 39 transcripts were modulated following minocycline treatment. We then selected two gene expression patterns, B1 and B2. The first transcription pattern, B1, consisted of 10 mRNA transcripts that increased in abundance after injury, and minocycline treatment reversed or inhibited the effect of the injury; the B2 transcription pattern consisted of 7 mRNA transcripts whose abundance decreased following sciatic nerve ligation, and minocycline treatment reversed the effect of the injury. Based on the literature, we selected seven genes for further analysis: Cd40, Clec7a, Apobec3b, Slc7a7, and Fam22f from pattern B1 and Rwdd3 and Gimap5 from pattern B2. Additionally, these genes were analyzed using quantitative PCR to determine the transcriptional changes strongly related to the development of neuropathic pain; the ipsilateral DRGs (L4-L6) were also collected and analyzed in these rats using qPCR.

Conclusion

In this work, we confirmed gene expression alterations previously identified by microarray analysis in the spinal cord and analyzed the expression of selected genes in the DRG. Moreover, we reviewed the literature to illustrate the relevance of these findings for neuropathic pain development and therapy. Further studies are needed to elucidate the roles of the individual genes in neuropathic pain and to determine the therapeutic role of minocycline in the rat neuropathic pain model.

Background

Studies conducted in recent years have shown that multiple endogenous factors initiate and regulate neuropathic pain, and neuroimmune interactions play important roles in this process, which may be one reason for the loss of efficacy for many analgesics [15]. Recent reports suggest that inhibitors of microglial activation and cytokine synthesis, including minocycline, pentoxyfylline and propentofylline, may significantly inhibit the development of neuropathic pain in animal models [610]. The effects of these substances are a result of the decreased secretion of numerous cytokines due to the reduced activation of the microglia or other cells [6, 7, 9, 11]. Administration of these inhibitors also enhances the effect of analgesic drugs, such as morphine [7, 1214].

Our previous data highlighted the importance of immune response- and microglia activation-related genes in the development of neuropathic pain in the spinal dorsal horn and their involvement in the persistence of its symptoms [4, 7, 14, 15]. Therefore, the aim of the present study was to determine which of the genes that are altered after sciatic nerve injury are modulated by minocycline. Minocycline is a strong modulator of the neuroimmune response and readily permeates the blood–brain barrier [12, 1621]. Minocycline is effective at reducing allodynia in animal models, and therefore, it appears to be a promising substance that could be used in analgesia, as confirmed by experimental studies [7, 15, 22]. This drug is also effective in the treatment of neurodegenerative diseases in a clinical setting [23]; however, it has not been used in the treatment of neuropathic pain. The neuroprotective properties of minocycline have been observed in spinal cord injury models [24], hypoxia models [25, 26] and in Huntington’s and Parkinson’s diseases in the clinical setting [20, 21, 25, 26]. Minocycline was shown to alleviate neuropathic pain symptoms by inhibiting the activation of the microglia within the spinal cord, leading to the reduced release of pro-inflammatory agents, which often exhibit pronociceptive properties [9, 27]. Additionally, minocycline is believed to promote the maintenance of the blood–brain barrier integrity by reducing the expression of chemokine receptors and metalloproteinases, as well as by reducing the production of reactive oxygen species [28]. By affecting the degree of activation of the microglia, minocycline enhances the analgesic effects of morphine, which has been shown to be less effective in neuropathic pain models [7, 1214, 29].

Microarray analysis for global gene expression indicated that both neuropathic and inflammatory pain are associated with a dramatic shift in the regulation of many genes in the spinal cord and dorsal root ganglia DRG [30, 31].

Therefore, the aim of the present study was to determine which of the genes that were altered after sciatic nerve injury are modulated by minocycline. As DNA microarray experiments allowed us to identify numerous transcripts that are important for the development of neuropathy and are possible targets for drug therapy, we planned the experiments to study the molecular basis of the inhibitory effects of minocycline on the development of neuropathic pain. We performed behavioral studies and microarray screening for genes in the lumbar section of the rat spinal cord in the rat model of neuropathic pain (chronic constriction injury to the sciatic nerve, CCI), and we identified genes that are associated with minocycline efficacy in neuropathic pain. The changes in the abundance of some transcripts in the ipsilateral dorsal part of the lumbar spinal cord after sciatic nerve ligation and the modulatory role of minocycline were confirmed by quantitative PCR analysis of samples independent from those used for the microarray analysis, as well as DRG samples from the same group of rats subjected to CCI.

Results

The effect of minocycline i.p. administration on the development of mechanical allodynia and thermal hyperalgesia in CCI-exposed rats

Unilateral, loose ligation of the sciatic nerve led to the development of symptoms typical of neuropathic pain, such as allodynia (Figure 1A) and hyperalgesia (Figure 1B). In the von Frey test, strong tactile allodynia on the paw ipsilateral to the injury was observed on day seven after CCI; at this time, the ipsilateral paw responded to a stimulation of 14.8 ± 0.7 g (Figure 1A), compared to the reactions of the hind paws of naive rats to 26.0 ± 0.7 g. The strongest hyperalgesia was observed on the seventh day in the cold plate test (Figure 1B). At this time, the ipsilateral paw reacted after 6.2 ± 0.7 s (Figure 1B), compared to reaction after 27.3 ± 2.4 s in the naive rats.

Figure 1
figure 1

The effect of minocycline on the development of allodynia and hyperalgesia 7 days after CCI. Influence of preemptive and repeated administration of minocycline (30 mg/kg; i.p.; 16 h and 1 h before CCI and then for 7 days twice daily) on the development of allodynia (A; von Frey test) and hyperalgesia (B; cold plate test) seven days after CCI in rats. The data are presented as the mean ± S.E.M (10–12 rat per group). Allodynia and hyperalgesia were assessed 60 min after drug administration. The inter-group differences were analyzed using an ANOVA and Bonferroni’s multiple comparison test; **p < 0.01 and ***p < 0.001 indicate a significant difference when compared to the control (naïve rats), and ###p < 0.001 indicates a significant difference when compared to the V-CCI rats (ANOVA, Bonferroni’s test). Naive (N), vehicle (V; water for injection), minocycline (MC).

Administration of minocycline led to a significant reduction in the above described symptoms. Preemptive and repeated treatment with minocycline (twice daily; 30 mg/kg i.p.) significantly attenuated the allodynia to 18.4 ± 1.5 g (Figure 1A) and the hyperalgesia to 19.3 ± 2.9 s (Figure 1B) on day seven after CCI.

Profiling the gene expression alterations following minocycline i.p. administration to CCI-exposed rats

The microarray method was used to study the effect of minocycline on the expression changes caused by injury to the sciatic nerve (Figure 2). The analysis involved three groups of animals: N – naïve, control animals; V-CCI – animals that received water for injection seven days after sciatic nerve ligation; and MC-CCI – animals that received repeated intraperitoneal administration of minocycline (30 mg/kg i.p.; 16 h and 1 h before nerve damage and then twice daily for seven days after CCI (Figure 2). The DNA microarray method (Illumina RatRef-12 V1) was used to analyze the changes in mRNA levels in the lumbar segment of the spinal cord. The results of the microarray analysis are presented in the form of a heat map. Statistical analysis (ANOVA) revealed that, out of 22,500 tested mRNA transcripts, 93 transcripts displayed changes in their relative abundance following CCI, at an ANOVA threshold value of p < 0.001. Of these, 54 transcripts were not affected by minocycline administration, while the levels of the other 39 transcripts were modulated following minocycline administration. The results are presented as two main gene transcription patterns: A and B (N vs. V-CCI vs. MC-CCI) (Figure 2A, B).

Figure 2
figure 2

Profiling the gene expression alterations following minocycline i.p. administration to CCI-exposed rats. Gene expression changes 7 days after nerve injury and administration of minocycline (30 mg/kg; i.p.; 16 h and 1 h before CCI and then for 7 days twice daily). The statistical analysis was performed at an ANOVA threshold value of p < 0.001 (N vs. V-CCI vs. MC-CCI). The results of the cDNA microarray analysis are presented in the form of a heat map. The individual columns represent the respective microarrays, while the individual rows represent the regulated genes. The intensity of the color reflects the relative abundance of the transcript and is proportional to the standard deviation from the in-row average, as indicated in the legend below the heat map. (A) The gene pattern represents changes in the abundance levels of 54 mRNA transcripts induced in the lumbar spinal segment on day 7 of neuropathic pain; these abundance levels of these genes were not affected by the addition of minocycline (30 mg/kg; 16 h and 1 h before sciatic nerve damage and then twice daily for 7 days). (B) Gene pattern B represents changes in the abundance levels of 39 mRNA transcripts induced in the lumbar spinal segment on day 7 of neuropathic pain; the abundance levels of these genes were modulated by repeated intraperitoneal administration of minocycline (30 mg/kg; 16 h and 1 h before sciatic nerve damage and then twice daily for 7 days). (C) The percentage analysis of the alterations in the CCI-induced gene expression profile in the spinal cord.

Transcripts not modulated by minocycline administration

Statistical analysis (ANOVA) revealed that 58% (54 transcripts) of the 93 examined transcripts (RGD1560961, Ltc4s, LOC308350, Btk, Ms4a6a, LOC681932, Ptprc, RGD1566043, Irf1, Angpt4, Psmb4, Ifi47, RGD1359108, Ccr5, Tifa, RGD1562655, Ifngr1, Cd14, Jak3, Lgals3bp, Hla-dma, RGD1564553, Fam105a, Irf5, Plek, Rac2, C4-2, St14, Myo1f, Adap2, Mefv, Apobec1, LOC497841, Pik3ap1, Lyl1, Lgmn, C1qa, C1qb, Clec4a3, Glipr1, Coro1a, C1qc, Plac8, Tspo, Laptm5, Arhgap9, Tbxas1, Fcgr2b, LOC498276, Ncf1, Tlr7, Hhex, Tlr2, and Sart2) were characterized by an increase in abundance following CCI, and minocycline administrations did not affect these changes (Figure 2A i C -pattern A). These tests were conducted at an ANOVA statistical threshold of p < 0.001 (N vs. V-CCI vs. MC-CCI).

Transcripts modulated by minocycline administration

Statistical analysis (ANOVA) revealed that 93 transcripts showed changes in their relative abundance following sciatic nerve ligation. Of these, the levels of 39 (42%) transcripts were modulated following repeated injections of minocycline, and these transcripts were described as having transcription pattern B (Figure 2B, C).

Transcription pattern B

Transcription pattern B1 – 10.75% of 93 transcripts

Ten transcripts (Kmo, Pcdha3, Clec7a, Lep, Apobec3b, Fam22f, Card11, Slc7a7, Cd40, and Etaa1) were characterized by an increase in mRNA abundance after CCI, and minocycline administrations reduced the level of these changes.

Transcription pattern B2 – 7.53% of 93 transcripts

Seven transcripts (Rwdd3, LOC364991, Gimap5, RGD1305604, Gpr110, LOC287028, and HeyI) were characterized by a decrease in mRNA abundance following CCI, and minocycline administrations reduced the level of these changes.

Transcription pattern B3 – 8.60% of 93 transcripts

Eight transcripts (RGD1559604, Ms4a6a, Btg2, LOC498606, Wdr45, Fn1, Abt1, and Pla1a) were characterized by an increase in mRNA abundance following CCI, and minocycline administrations enhanced the level of these changes.

Transcription pattern B4 – 15.05% of 93 transcripts

Fourteen transcripts (Ldhb, LOC501111, LOC688240, Ebp, LOC499941, Nkx2-2, Dhcr7, Cnp, Rtn4, Tnnc2, Slc22a8, LOC290549, Car4, and Idi1) were characterized by high baseline mRNA levels, which decreased following CCI and continued to decrease as a result of minocycline administrations.

qPCR validation of the mRNA transcript changes in the spinal cord and DRG following minocycline i.p. administration seven days after CCI

Genes from the B1 transcription pattern

Genes similarly regulated by minocycline in the spinal cord and DRG

Cd40

In the dorsal lumbar spinal cord, we observed compared to naïve rats the upregulation of monocyte marker Cd40 mRNA (1.0 ± 0.02 vs. 1.3 ± 0.009) using microarray analysis (Figure 3A). Microarray analysis of gene expression for T-cell (Cd3g, Cd3e, Cd3d, CD4, and CD8), B-cells (CD19) and NK-cells (CD335) markers suggest that there is no activation or infiltration of those cells into the spinal cord (Table 1). Minocycline significantly diminished the spinal level of Cd40 mRNA (from 1.3 ± 0.009 to 1.1 ± 0.02).

Figure 3
figure 3

Genes from the B1 transcription pattern - similarly spinally and DRG regulated by minocycline. Comparison of the CCI-induced changes in transcript levels of genes with expression pattern B1 (Cd40, Clec7a and Apobec3b) in the lumbar spinal cord using microarray (A, D, G) and qPCR (B, E, H) analysis; transcript levels in the DRG (C, F, I) and the modulation by repeated intraperitoneal administration of minocycline (30 mg/kg; i.p. 16 h and 1 h before sciatic nerve damage and then twice daily for 7 days) were also examined. The statistical analysis was performed using an ANOVA and Bonferroni’s test. **p < 0.01 and ***p < 0.001 indicate a significant difference when compared to the control (naïve rats); #p < 0.05, ##p < 0.01 and ###p < 0.001 indicate a significant difference when compared to the V-CCI rats. Naive (N), vehicle (V; water for injection), minocycline (MC).

Table 1 The influence of minocycline on lymphocyte and monocyte cell markers 7 days after CCI

qPCR analysis (Figure 3B) confirmed the changes detected by microarray. In the spinal cord, the upregulation of Cd40 mRNA (1.0 ± 0.05 vs. 3.4 ± 0.2) was observed compared to the naïve rats. Minocycline significantly diminished the spinal level of Cd40 mRNA (from 3.4 ± 0.2 to 2.5 ± 0.3). In the DRG (Figure 3C), the upregulation of Cd40 mRNA (1.0 ± 0.04 vs. 1.4 ± 0.04) was observed compared to the naïve rats and minocycline significantly diminished the level of Cd40 (from 1.4 ± 0.04 to 1.17 ± 0.05).

Clec7a

Using microarray analysis, in the dorsal lumbar spinal cord, we observed compared to naïve rats the upregulation of Clec7a mRNA (1.0 ± 0.009 vs. 1.9 ± 0.06) (Figure 3D) and minocycline significantly diminished the level of Clec7a mRNA (from 1.9 ± 0.06 to 1.5 ± 0.04).

qPCR analysis (Figure 3E) confirmed the changes detected by microarray. In the spinal cord, the upregulation of Clec7a mRNA (1.0 ± 0.009 vs. 13.4 ± 2.00) was observed compared to the naïve rats, and treatment with minocycline significantly diminished the level of Clec7a mRNA (from 13.4 ± 2.00 to 7.2 ± 1.1). In the DRG (Figure 3F), the upregulation of Clec7a mRNA (1.0 ± 0.13 vs. 1.8 ± 0.12) was observed compared to the naïve rats, and repeated treatment with minocycline significantly diminished the level of Clec7a mRNA (from 1.8 ± 0.12 to 1.14 ± 0.17).

Apobec3b

In the dorsal lumbar spinal cord, we observed compared to naïve rats the upregulation of Apobec3b mRNA (1.0 ± 0.008 vs. 1.3 ± 0.03) using microarray analysis (Figure 3G) and repeated treatment with minocycline significantly diminished the level of Apobec3b mRNA (from 1.35 ± 0.03 to 1.17 ± 0.007).

qPCR analysis (Figure 3H) confirmed the changes detected by microarray. In the spinal cord, the upregulation of Apobec3b mRNA (1.0 ± 0.2 vs. 2.8 ± 0.4) was observed compared to the naïve rats, and repeated treatment with minocycline significantly diminished the level of Apobec3b mRNA (from 2.8 ± 0.4 to 1.8 ± 0.002). In the DRG (Figure 3I), Apobec3b mRNA (1.0 ± 0.2 vs. 4.5 ± 0.3) was upregulated compared to the naïve rats, and repeated treatment with minocycline significantly diminished the level of Apobec3b in the DRG (from 4.5 ± 0.3 to 2.3 ± 0.4).

Genes differentially regulated by minocycline in the spinal cord and DRG

Slc7a7

Using microarray analysis (Figure 4A), we observed compared to naïve rats the spinal upregulation of Slc7a7 mRNA (1.0 ± 0.01 vs. 1.7 ± 0.03) and repeated treatment with minocycline diminished the level of Slc7a7 mRNA (from 1.7 ± 0.03 to 1.5 ± 0.07).

Figure 4
figure 4

Genes from the B1 transcription pattern - differentially spinally and DRG regulated by minocycline. Comparison of the changes in the transcripts levels of genes with expression pattern B1 (Slc7a7 and Fam22f) in the lumbar spinal cord using microarray (A, D) and qPCR (B, E) analysis; transcript levels in the DRG (C, F) under neuropathic pain and the modulation by repeated intraperitoneal administration of minocycline (30 mg/kg; i.p. 16 h and 1 h before sciatic nerve damage and then twice daily for 7 days) were also examined using qPCR analysis. The statistical analysis was performed using an ANOVA and Bonferroni’s test. ***p < 0.001 indicates a significant difference when compared to the control (naïve rats); #p < 0.05, ##p < 0.01 and ###p < 0.001 indicate a significant difference when compared to the V-CCI rats. Naive (N), vehicle (V; water for injection), minocycline (MC).

qPCR analysis (Figure 4B) confirmed the changes detected by microarray, but only in the spinal cord samples. In the spinal cord, Slc7a7 mRNA was upregulated (1.0 ± 0.04 vs. 3.3 ± 0.2) compared to the naïve rats and repeated treatment with minocycline significantly diminished the level of Slc7a7 mRNA (from 3.3 ± 0.2 to 2.1 ± 0.2). In the DRG (Figure 4C), the downregulation of Slc7a7 mRNA was observed (1.0 ± 0.1 vs. 0.6 ± 0.02), compared to naïve rats, and repeated treatment with minocycline increased the level of Slc7a7 (from 0.6 ± 0.02 to 0.8 ± 0.04).

Fam22f

In the dorsal lumbar spinal cord, we observed compared to naïve rats the upregulation of Fam22f mRNA (1.0 ± 0.01 vs. 1.9 ± 0.06) using microarray analysis (Figure 4D) and repeated treatment with minocycline diminished the level of Fam22f mRNA (from 1.9 ± 0.06 to 1.5 ± 0.04).

qPCR analysis (Figure 4E) confirmed the changes detected by microarray, but only in the spinal cord. In the spinal cord, Fam22f mRNA (1.0 ± 0.02 vs. 2.5 ± 0.3) was upregulated compared to the naïve rats, and repeated treatment with minocycline significantly diminished the spinal level of Fam22f mRNA (from 2.4 ± 0.3 to 1.4 ± 0.2). In contrast, in the DRG (Figure 4F), the downregulation of Fam22f mRNA (1.0 ± 0.1 vs. 0.7 ± 0.03) was observed compared to the naïve rats, and repeated treatment with minocycline increased (from 0.7 ± 0.03 to 0.8 ± 0.03) the level of Fam22f.

Genes from the B2 transcription pattern

Genes similarly regulated by minocycline in the spinal cord and DRG

Rwdd3

In the dorsal lumbar spinal cord, we observed compared to naïve rats a slight downregulation of Rwdd3 mRNA (1.0 ± 0.004 vs. 0.8 ± 0.006) using microarray analysis (Figure 5A), and repeated treatment with minocycline resulted in an increase (from 0.8 ± 0.006 to 0.9 ± 0.008) in the level of Rwdd3 mRNA.

Figure 5
figure 5

Genes from the B2 transcription pattern - similarly spinally and DRG regulated by minocycline. Comparison of the changes in the transcript levels of genes with expression pattern B2 (Rwdd3 and Gimap5) in the lumbar spinal cord using microarray (A, D) and qPCR (B, E) analysis; transcript levels in the DRG (C, F) under neuropathic pain and the modulation by repeated intraperitoneal administration of minocycline (30 mg/kg; i.p. 16 h and 1 h before sciatic nerve damage and then twice daily for 7 days) were also examined using qPCR. The statistical analysis was performed using an ANOVA and Bonferroni’s test. **p <0.01 and ***p <0.001 indicate a significant difference when compared to the control (naïve rats); #p < 0.05, ##p < 0.01 and ###p < 0.001 indicate a significant difference when compared to the V-CCI rats. Naive (N), vehicle (V; water for injection), minocycline (MC).

qPCR analysis (Figure 5B) confirmed the changes detected by microarray. In the spinal cord, Rwdd3 mRNA (1.0 ± 0.006 vs. 0.8 ± 0.05) was downregulated compared to the naïve rats, and repeated treatment with minocycline significantly increased the level of Rwdd3 mRNA (from 0.8 ± 0.05 to 1.0 ± 0.08). In the DRG (Figure 5C), Rwdd3 mRNA (1.0 ± 0.006 vs. 0.75 ± 0.02) was downregulated compared to the naïve rats, and repeated treatment with minocycline did not influence the level of Rwdd3 (0.75 ± 0.02 vs. 0.72 ± 0.01).

Gimap5

In the dorsal lumbar spinal cord, we observed compared to naïve rats the downregulation of Gimap5 mRNA (1.0 ± 0.01 vs. 0.7 ± 0.003) using microarray analysis (Figure 5D). Chronic treatment with minocycline decreased (from 0.7 ± 0.003 to 0.8 ± 0.02) the level of Gimap5 mRNA.

qPCR analysis (Figure 5E) confirmed the changes detected by microarray. In the spinal cord, the downregulation of Gimap5 mRNA (1.0 ± 0.06 vs. 0.7 ± 0.04) was observed compared to the naïve rats, and repeated treatment with minocycline significantly increased the level of Gimap5 mRNA (from 0.7 ± 0.04 to 0.9 ± 0.08). In the DRG (Figure 5F), Gimap5 mRNA was downregulated (1.0 ± 0.06 vs. 0.5 ± 0.04) compared to the naïve rats, and treatment with minocycline increased (from 0.5 ± 0.04 to 0.7 ± 0.04) the level of Gimap5.

Discussion

The widely accepted definition of neuropathic pain is “pain arising as a direct consequence of a lesion or disease affecting the somatosensory system” [32]. Neuropathic pain is very challenging to manage because of the heterogeneity of its etiologies, symptoms and underlying mechanisms, as well as our limited understanding of its development and progression [33]. Hypersensitivity to thermal and mechanical nociceptive stimuli in animals was observed on day 7 following sciatic nerve injury in our experiments, and this result is in agreement with our previous studies, as well as others [7, 15, 3438]. Recently, accumulating evidence has suggested that glial cell activation and neuroinflammation are critical for the development and maintenance of persistent pain [1, 4, 5, 14, 39]. The results obtained after administration of minocycline in neuropathic pain models suggest that it has a therapeutic potential; however, it is well known that the its mechanism of action is not selective. In our previously published studies, we demonstrated an increase in the activation of monocytes on day 7 after sciatic nerve injury in the lumbar spinal cord and/or in the DRG [7, 15]. Other studies revealed microglia activation on day 2 following sciatic nerve injury, with its highest activation being observed between days 7 and 10 [1, 7, 15, 40, 41]; therefore, we choose day 7 for our microarray analysis. Although it is evident that the main cause of the attenuation of neuropathy is the inhibition of microglial activation, minocycline also acts through many other targets. Our studies allowed us to isolate a series of genes whose expression profiles followed several patterns after CCI and minocycline treatment. As shown by our studies, 93 out of 22,500 studied transcripts undergo abundance level changes following sciatic nerve injury; of these, 54 transcripts were not affected by the repeated administration of minocycline, while strong modulation was observed in 39 transcripts. These 39 genes are potential targets for minocycline and could be interesting from a therapeutic viewpoint.

The results obtained from microarray analysis revealed that the abundance levels of 54 transcripts increased following sciatic nerve ligation, but these genes were not affected by the repeated administration of minocycline; these genes are shown in Figure 2A in the results section. The roles that many of these genes play in neuropathy are not known; therefore, these genes will be the subjects of our future research. In the present study, we validated the expression changes of seven of the most interesting genes from the 39 genes whose abundance levels increased after sciatic nerve injury and were modulated by repeated administration of minocycline using qPCR analysis of spinal cord samples, and we further studied the mRNA levels of these genes in DRG samples from the same experimental scheme. Furthermore, we selected two gene expression patterns to study in detail: B1 and B2. The B1 transcription pattern contained 10 transcripts; the mRNA abundance levels of these genes increased following injury, and minocycline administration reversed or inhibited the effect of the injury. The B2 transcription pattern contained 7 transcripts; the mRNA abundance levels of these genes decreased after sciatic nerve ligation, and minocycline administration reversed the effect of the injury. According to the literature, we selected the following genes for qPCR analysis: Cd40, Clec7a, Apobec3b, Slc7a7, and Fam22f from pattern B1 and Rwdd3 and Gimap5 from pattern B2.

Cd40, Clec7a and Apobec3b are genes from the B1 transcription pattern that are upregulated by sciatic nerve injury in the spinal cord and DRG and are diminished by minocycline

CD40, a 48 kDa cell surface tumor necrosis factor (TNF) family receptor, has been shown to be upregulated in microglia upon activation in both in vitro and in vivo studies. CD40 is also expressed by a wide variety of cells, such as neurons, dendritic cells, microglia, B cells, macrophages, keratinocytes, endothelial cells, thymic epithelial cells, fibroblasts and various tumor cells [42]. It is known that CD40-mediated microglia activation contributes to disease progression in a variety of neuroinflammatory diseases, such as multiple sclerosis, Alzheimer’s disease, and cerebral ischemia, leading to the production of a wide array of cytokines, chemokines, matrix metalloproteinases and neurotoxins [43, 44]. After CD40 is bound, numerous signaling pathways are activated, leading to changes in gene expression and function. The interaction between CD40 and CD154 appears to be critical for a productive immune response, the upregulation of various costimulatory molecules (ICAM-1, VCAM-1, E-selectin, LFA-3, B7.1, B7.2, and class II MHC) and the production of numerous cytokines/chemokines (IL-1, IL-6, IL-8, IL-10, IL-12, TNF-α, MIP-1α, and MCP-1). Thus, signaling through CD40 in macrophages/microglia induces a number of soluble mediators that have important functional roles in the CNS (central nervous system). Macrophages/microglia have been shown to express CD40 in patients with multiple sclerosis [45], and CD40-deficient mice fail to develop these diseases [46, 47]. The strategies used to attenuate inflammatory responses within the CNS by inhibiting the activation of macrophages and microglia (by suppressing CD40 expression) may be beneficial for a growing number of neuroinflammatory diseases. Others have shown that TGF-β and IL-4 inhibit the IFN-γ-induced CD40 expression in microglia [48, 49]. In 2012, Cao et al. [50] showed that CD40 plays an important role in leukocyte infiltration into the lumbar spinal cord after L5 spinal nerve transection. Studies using a CD40 neutralizing Ab suggest that CD40 is required early on in order to promote the maintenance of injury-induced mechanical hypersensitivity. These data are consistent with our data, which showed that intraperitoneal administration of minocycline inhibited the development of mechanical allodynia and thermal hyperalgesia in parallel with the observed downregulation of CD40.

Clec7a, also called Dectin-1, was recently identified as the most important receptor for beta-glucan; Clec7a is a type II transmembrane protein that binds beta-1,3 and beta-1,6 glucans. Clec7a is primarily expressed by cells of myeloid origin, including monocytes, macrophages, microglia, neutrophils, most subsets of dendritic cells and a subset of T cells, B cells, mast cells, and eosinophils. Our studies using qPCR demonstrated that mRNA for Clec7a is expressed in rat primary microglial cells cultures (date not shown). Clec7a can recognize an unidentified endogenous ligand on T cells and may act as a costimulatory molecule; it can also induce a variety of cellular responses, including phagocytosis, respiratory burst and cytokine production [51]. Recently, Salazar-Aldrete et al. [52] showed that monocytes from patients with systemic lupus erythematosus and rheumatoid arthritis exhibited decreased expression of Clec7a in parallel with the enhanced synthesis of proinflammatory cytokines. Recent studies have highlighted the importance of Clec7a in anti-fungal immunity in both mice and humans and have suggested the possible involvement of this receptor in the control of mycobacterial infections [53]. Using microarray and qPCR analyses, we have shown that Clec7a mRNA is upregulated in the spinal cord and DRG at seven days after CCI-induced neuropathic pain and is downregulated after intraperitoneal administration of minocycline in the spinal cord and DRG.

Apobec3b is overexpressed in many tumors, including tumors found in the breast, cervix, bladder, lung, head and neck [5457]. The overexpression of Apobec3b in animal models confirmed its tumor-type specificity [58]. It is not known which cells expressed Apobec3b, but after observing that minocycline downregulates the injury-increased level of Apobec3b mRNA, we hypothesize that microglia cells may be one type of cell that can express Apobec3b. Further studies are needed to clarify the role of this gene in the development of neuropathy. We observed a very similar pattern of Apobec3b transcription, as well as its regulation by minocycline, in the spinal cord and dorsal root ganglia, which suggest that microglia as well as other cell types, such as macrophages and leukocytes, can be the source of Apobec3b [58].

Slc7a7 and Fam22f are genes from the B1 transcription pattern that are upregulated by sciatic nerve injury in the spinal cord, are downregulated in the DRG and are modulated by minocycline in both structures

Slc7a7 (solute carrier family 7, amino acid transporter light chain, y + L system, member 7) is an important amino acid transporter responsible for the sodium-independent influx/efflux of cationic and large neutral amino acids across the membrane [59]. Deregulation of the amino acid transporter Slc7a7 is involved in multiple types of cancer, including gliobastoma, non-small cell lung cancer and multiple myeloma [6062]. Studies conducted in recent years have demonstrated that diseases associated with Slc7a7 also include lysinuric protein intolerance and cystinuria [63]. There are no reports detailing the roles Slc7a7 plays during neuropathic pain, inflammation or neurodegenerative disorders. Our study suggested that Slc7a7 may play a role after nerve injury. In the present study, using microarray and qPCR analyses, we showed that the level of Slc7a7 mRNA is upregulated in the spinal cord seven days after CCI-induced neuropathic pain but is downregulated in the DRG. Repeated administration of minocycline reversed both of these changes. The regulation of this gene by minocycline in the spinal cord suggests it is present in microglia/macrophages. However, the involvement of the mononuclear phagocyte system appears to play a crucial role in the development of immunological complications, as was previously suggested by Barilli et al. [64]. The clinical significance of Slc7a7 expression in pain therapy needs to be clarified.

Fam22f’s role in neuropathic pain and/or neurodegenerative disorders needs to be described, as there is no information on this gene in the literature. Until now, it was only known that Fam22f belongs to the Fam22 family of transmembrane proteins, which span from one side of the membrane to the other side of the membrane. We tested this gene in the present study, and we demonstrated that Fam22f mRNA is upregulated in the spinal cord seven days after CCI-induced neuropathic pain but is downregulated in the DRG. Repeated administration of minocycline reversed both of these changes. The contrasting expression changes in the spinal cord and DRG suggest the participation of various cell types that express the gene in the tested structures after injury, but the role of minocycline in the expression of this gene remains to be explained.

Rwdd3 and Gimap5 are genes from the B2 transcription pattern that are downregulated by sciatic nerve injury but are increased after minocycline treatment in the spinal cord and DRG

Rwdd3 belongs to the RWD domain family and was first described in 1996 by Bonaldo [65]. Recently, in 2013, Chi-Cheng Huang et al. [66] identified Rwdd3 as a gene associated with breast cancer risk. Some authors have suggested a role for Rwdd3 in neurosensory hearing loss and cellular stress [67]. Recently, Bergmann et al. [68] examined the association between the single nucleotide polymorphism rs2296308 in Rwdd3 and the development of neuropathy in paclitaxel-treated cancer patients, as was suggested by Schneider et al. [69]; however, this study was unable to confirm such a correlation. Our studies suggest the importance of Rwdd3 in the development of neuropathic pain, as we observed the CCI-induced upregulation of this gene, as well as its downregulation by chronic minocycline treatment.

Gimap5 refers to GTPase of immunity-associated nucleotide binding protein 5. Gimap5 is one of seven members of the Gimap family, which has been shown to be integral to T cell survival and development. These small GTPases regulate proapoptotic and antiapoptotic T cell pathways [7073], as well as thymocyte maturation and differentiation [72, 74]. There is a lack of information concerning the role or even regulation of Gimap5 during neuropathic pain; however, its important role in the function of the immune system suggests its importance. Our microarray and qPCR data strongly suggest it plays a significant role in the development of neuropathic pain. Our study using qPCR revealed the absence of Gimap5 mRNA in rat primary microglial cell culture (data not shown). Until now, it has been known that a loss of Gimap5 function causes T cell lymphopenia in rats due to the near complete loss of post-thymic peripheral CD8 T cells, which triggers a lethal autoimmune disease [75]. It has also been shown that Gimap5 knockout mice lack peripheral NK cells and CD8+ T cells and exhibit dynamic changes in immune homeostasis, as demonstrated by the progressive loss of CD4+ T cells and B cells and the development of neutrophilia [76, 77]. A frameshift mutation in Gimap5 is a prerequisite for the development of spontaneous type 1 diabetes in rats [72, 73]. In 2011, Moralejo et al. [78] suggested that further elucidating the role of the Gimap5 in T cell survival, development and/or activation would aid in our understanding of the pathways involved in the onset of spontaneous diabetes mellitus type 1 and may ultimately uncover the pathways leading to development of this disease in humans. The mechanisms underlying the pro-survival function of Gimap5 in T cells have not been elucidated. Recently, some authors have shown that Gimap5-deficiency in T cells impairs Ca2+ entry via the plasma membrane channels [77]. Gimap5sph/sph mice, an ENU germline mutant with a missense mutation in Gimap5, showed a progressive loss of the peripheral lymphocyte populations and developed spontaneous colitis, resulting in early mortality [79]. Genetic aberrations in Gimap5 have been linked to lymphopenia and the loss of immunological tolerance. Gimap5 is essential for maintaining lymphocyte quiescence and immunological tolerance. Its role during neuropathy needs further explanation and will be studied in our laboratory in the near future. In the present study, we have shown that Gimap5 mRNA is downregulated in the spinal cord and DRG seven days after CCI-induced neuropathic pain. We also showed that intraperitoneal administration of minocycline inhibited the development of neuropathic pain symptoms and in parallel we observed that the level of Gimap5 mRNA back to the level measured in naïve rats on the spinal cord level and DRG.

Conclusions

To summarize, the DNA microarray method allowed us to choose 93 transcripts from 22,500 studied transcripts due to their abundance level changes after sciatic nerve injury, as well as 39 (42%) transcripts that were additionally modulated by repeated administration of minocycline. Further studies are needed to elucidate the roles of the individual genes in the development of neuropathic pain and to determine which genes altered by minocycline may be relevant to its analgesic action in neuropathic pain. The gene analysis in the present study offers the first step towards future research into genes that can be modulated and may be good targets for effective and safe therapies for the treatment of neuropathic pain.

Materials and methods

Animals

Male Wistar rats (300–350 g) from Charles River (Hamburg, Germany) were housed in cages lined with sawdust under a standard 12/12 h light/dark cycle (lights on at 08:00 h), with ad libitum access to food and water. All efforts were made to minimize animal suffering and to reduce the number of animals used in this study. All experiments were performed according to the recommendations of the International Association for the Study of Pain (IASP) [80] and the National Institutes of Health Guide for the Care and Use of Laboratory Animals, and these experiments were approved by the local Bioethics Committee (Krakow, Poland).

Surgical preparations

Chronic constriction injury (CCI) was performed according to Bennett and Xie [81]. The right sciatic nerve was exposed under sodium pentobarbital anesthesia (60 mg/kg; i.p.). Four ligatures (4/0 silk) were made around the nerve, distal to the sciatic notch with 1 mm spacing, until a brief twitch in the respective hind limb was observed. After surgery, all rats developed symptoms of long-lasting neuropathic pain, such as allodynia and hyperalgesia.

Drug administration and experimental scheme

Minocycline hydrochloride (30 mg/kg; Sigma, Schnelldorf, Germany) was dissolved in water for injection and administered preemptively by i.p. injection 16 and 1 h before CCI, then twice daily for seven days, as previously described [7, 9]. This administration schedule was used because systemic microglia inhibitors attenuate the activation of microglia more efficiently when the inhibitor is injected before injury [9, 27, 82]. The control groups received the vehicle (water for injection) on the same schedule (Scheme 1).

Scheme 1
scheme 1

Drug administration and experimental schedule during 7 days after CCI.

Drug administration and experimental schedule during 7 days after CCI.

According to the scheme, naive animals were subjected to all procedures (behavioral tests and tissue collection) parallel to the V-CCI and MC-CCI groups, with exception of CCI surgery procedure including anaesthesia. The influence of minocycline administration on gene expression in the ipsilateral dorsal part of the lumbar spinal cord (L4–L6) was studied using DNA microarray analysis. Injury-induced changes and the influence of minocycline were verified using the qPCR method, and the gene expression patterns in the ipsilateral DRG (L4–L6) were also examined.

Behavioral tests

Tactile allodynia (von Frey test)

Allodynia was measured using an automatic von Frey apparatus (Dynamic Plantar Aesthesiometer Cat. No. 37400, Ugo Basile, Italy). The animals were placed in plastic cages with wire net floors 5 min before the experiment. The strengths of the von Frey stimuli used in our experiments ranged from 0.5 to 26 g. The von Frey test began with the lowest filament and increased in order; the filament was applied to the midplantar surface of the hind paw, and measurements were taken automatically, as described previously [7, 15]. The ipsilateral paw was tested, and the mean value was calculated. There was almost no response to the highest strength (26 g) in the naive animals. Therefore, a line was drawn at this value. For the CCI-exposed rats, the significantly different reactions of the paws between the CCI-exposed and naive rats were recorded.

Cold hyperalgesia (cold plate test)

Hyperalgesia was assessed using the cold plate test (Cold/Hot Plate Analgesia Meter No. 05044, Columbus Instruments, USA) as previously described [7, 15]. The temperature of the cold plate was kept at 5°C, and the cut-off latency was 30 s. The animals were placed on the cold plate, and the time until the hind paw was lifted was recorded. In the naive rat group, the reaction of the first hind paw to be lifted was measured. In the rats subjected to nerve injury, the ipsilateral paw reacted first.

Biochemical tests

Tissue collection and RNA isolation

Ipsilateral and contralateral fragments of the dorsal part of the lumbar (L5–L6) spinal cord and the ipsi- and contralateral DRG (L5–L6) were removed immediately after decapitation on day 7 after CCI. The tissue samples were placed in individual tubes containing the tissue storage reagent RNAlater (Qiagen Inc.) and were stored at −70°C for RNA isolation. Total RNA was extracted using the TRIzol reagent (Invitrogen), as previously described [83]. The RNA concentration was measured using a NanoDrop ND-1000 Spectrometer (NanoDrop Technologies), and RNA quality was determined by chip-based capillary electrophoresis using an RNA 6000 Nano LabChip Kit and an Agilent Bioanalyzer 2100 (Agilent) according to the manufacturer’s instructions.

Quantitative reverse transcriptase polymerase chain reaction (qPCR)

Reverse transcription was performed on 2 μg of total RNA using Omniscript reverse transcriptase (Qiagen Inc.) at 37°C for 60 min. RT reactions were carried out in the presence of an RNAse inhibitor (rRNAsin, Promega) and an oligo (dT16) primer (Qiagen Inc.). cDNA was diluted 1:10 with H2O, and for each reaction, ~ 50 ng of cDNA synthesized from the total RNA of an individual animal was used for the quantitative real-time PCR (qPCR) reaction. qPCR was performed using Assay-On-Demand TaqMan probes according to the manufacturer’s protocol (Applied Biosystems), and the reactions were run on an iCycler device (BioRad, Hercules). The following TaqMan primers and probes were used: Rn01527838_g1 (Hprt, hypoxanthine guanine rat hypoxanthine guanine phosphoribosyl transferase); Rn01522736_m1 (Rwdd3); Rn01772952_m1 (Fam22f); Rn00580189_m1 (Slc7a7); Rn01423590_m1 (Cd40); Rn00595553_m1 (Gimap5); Rn01505455_m1 (Apobec3b); and Rn01459401_m1 (Clec7a). The amplification efficiency for each assay (between 1.7 and 2) was determined by running a standard dilution curve. The cycle threshold values were calculated automatically by the iCycler IQ 3.0 software using the default parameters. RNA abundance was calculated as 2-(threshold cycle). HPRT transcript levels do not significantly change in rats exposed to CCI [15] and, therefore, served as an adequate housekeeping gene.

Microarray analysis

A starting amount of 200 ng of high quality total RNA was used to generate cDNA and cRNA using the Illumina TotalPrep RNA Amplification Kit (Illumina Inc., San Diego, CA, USA) according to the our previous studies [84, 85]. The obtained cDNA served as a template for in vitro transcription with T7 RNA polymerase and biotin UTP to generate multiple copies of biotinylated cRNA. Each cRNA sample (1.5 μg) was hybridized overnight to a RatRef-12 V1 BeadChip array (Illumina); subsequently, the chips were washed, dried and scanned using the BeadArray Reader (Illumina). Raw microarray data were generated using BeadStudio v3.0 (Illumina). Samples from 2 rats were pooled for each microarray, and 3 biological replicates were used for each experimental point. Microarray quality control was performed using the BeadArray R package v1.10.0. The following parameters were checked: number of outliers, number of beads and percentage of detected probes. After background subtraction, the data were normalized using quantile normalization and were then log2-transformed. The obtained signal was taken as the measure of mRNA abundance derived from the level of its gene expression. All statistical analyses were performed using the R software version 2.11.1.

Data analysis

The mean ± SEM of the behavioral data are presented in grams and seconds, and each group contained 6–15 rats. The results of the experiments were evaluated using one-way analysis of variance (ANOVA). The data are presented as the mean ± SEM. *p < 0.05, **p < 0.01 and ***p < 0.001 indicate a significant difference when compared to the control group (naïve rats); #p < 0.05 and ###p < 0.001 indicate a significant difference when compared to the vehicle-treated CCI-exposed rats.

The results of the qPCR analyses are presented as the fold change compared to the control group (naive rats) and were calculated for the ipsilateral sides of the spinal cords, as well as for the DRGs of the CCI-exposed rats. The quantitative qPCR analysis data are presented as the mean ± SEM and represent the normalized averages derived from the threshold cycle in the qPCR analysis of 4 to 10 samples for each group. The inter-group differences were analyzed using an ANOVA, followed by Bonferroni’s multiple comparison test. *p < 0.05, **p < 0.01 and ***p < 0.001 indicate a significant difference when compared to the control group (naïve rats); #p < 0.05 and ###p < 0.001 indicate a significant difference when compared to the vehicle-treated CCI-exposed rats.

Microarray analyses were performed in three groups: naïve, CCI-exposed and minocycline-treated CCI-exposed rats. The data are presented as fold changes compared to the naïve rats in the ipsilateral dorsal lumbar spinal cord and DRG. The inter-group differences were analyzed using ANOVAs, followed by the calculation of the false discovery rate (FDR). **p < 0.01 and ***p < 0.001 indicate significant differences when compared to the naïve rats. #p < 0.05, ##p < 0.01 and ###p < 0.001 indicate significant differences when compared to the CCI-treated group.

Abbreviations

CCI:

Chronic constriction injury to the sciatic nerve

CNS:

Central nervous system

DRG:

Dorsal root ganglia

N:

Naive

MC:

Minocycline

IASP:

International Association for the Study of Pain

V:

vehicle.

References

  1. Austin PJ, Moalem-Taylor G: The neuroimmune balance in neuropathic pain: involvement of inflammatory immune cells, immune-like glial cells and cytokines. J Neuroimmunol 2010, 229: 26–50.

    Article  CAS  PubMed  Google Scholar 

  2. DeLeo JA, Yezierski RP: The role of neuroinflammation and neuroimmune activation in persistent pain. Pain 2001, 90: 1–6.

    Article  CAS  PubMed  Google Scholar 

  3. Lee ML, Sciorra LJ: Partial monosomy of the long arm of chromosome 11 in a severely affected child. Ann Genet 1981, 24: 51–53.

    CAS  PubMed  Google Scholar 

  4. Mika J, Zychowska M, Popiolek-Barczyk K, Rojewska E, Przewlocka B: Importance of glial activation in neuropathic pain. Eur J Pharmacol 2013, 716: 106–119.

    Article  CAS  PubMed  Google Scholar 

  5. Watkins LR, Maier SF: Glia: a novel drug discovery target for clinical pain. Nat Rev Drug Discov 2003, 2: 973–985.

    Article  CAS  PubMed  Google Scholar 

  6. Lundblad R, Ekstrom P, Giercksky KE: Pentoxifylline improves survival and reduces tumor necrosis factor, interleukin-6, and endothelin-1 in fulminant intra-abdominal sepsis in rats. Shock 1995, 3: 210–215.

    Article  CAS  PubMed  Google Scholar 

  7. Mika J, Osikowicz M, Makuch W, Przewlocka B: Minocycline and pentoxifylline attenuate allodynia and hyperalgesia and potentiate the effects of morphine in rat and mouse models of neuropathic pain. Eur J Pharmacol 2007, 560: 142–149.

    Article  CAS  PubMed  Google Scholar 

  8. Neuner P, Klosner G, Schauer E, Pourmojib M, Macheiner W, Grünwald C, Knobler R, Schwarz A, Luger TA, Schwarz T: Pentoxifylline in vivo down-regulates the release of IL-1 beta, IL-6, IL-8 and tumour necrosis factor-alpha by human peripheral blood mononuclear cells. Immunology 1994, 83: 262–267.

    CAS  PubMed Central  PubMed  Google Scholar 

  9. Raghavendra V, Tanga F, Rutkowski MD, DeLeo JA: Anti-hyperalgesic and morphine-sparing actions of propentofylline following peripheral nerve injury in rats: mechanistic implications of spinal glia and proinflammatory cytokines. Pain 2003, 104: 655–664.

    Article  CAS  PubMed  Google Scholar 

  10. Sweitzer SM, Schubert P, DeLeo JA: Propentofylline, a glial modulating agent, exhibits antiallodynic properties in a rat model of neuropathic pain. J Pharmacol Exp Ther 2001, 297: 1210–1217.

    CAS  PubMed  Google Scholar 

  11. Sweitzer S, Martin D, DeLeo JA: Intrathecal interleukin-1 receptor antagonist in combination with soluble tumor necrosis factor receptor exhibits an anti-allodynic action in a rat model of neuropathic pain. Neuroscience 2001, 103: 529–539.

    Article  CAS  PubMed  Google Scholar 

  12. Cui Y, Liao XX, Liu W, Guo RX, Wu ZZ, Zhao CM, Chen PX, Feng JQ: A novel role of minocycline: attenuating morphine antinociceptive tolerance by inhibition of p38 MAPK in the activated spinal microglia. Brain Behav Immun 2008, 22: 114–123.

    Article  CAS  PubMed  Google Scholar 

  13. Hains BC, Waxman SG: Activated microglia contribute to the maintenance of chronic pain after spinal cord injury. J Neurosci 2006, 26: 4308–4317.

    Article  CAS  PubMed  Google Scholar 

  14. Mika J, Osikowicz M, Rojewska E, Korostynski M, Wawrzczak-Bargiela A, Przewlocki R, Przewlocka B: Differential activation of spinal microglial and astroglial cells in a mouse model of peripheral neuropathic pain. Eur J Pharmacol 2009, 623: 65–72.

    Article  CAS  PubMed  Google Scholar 

  15. Mika J, Rojewska E, Makuch W, Przewlocka B: Minocycline reduces the injury-induced expression of prodynorphin and pronociceptin in the dorsal root ganglion in a rat model of neuropathic pain. Neuroscience 2010, 165: 1420–1428.

    Article  CAS  PubMed  Google Scholar 

  16. Amin AR, Attur MG, Thakker GD, Patel PD, Vyas PR, Patel RN, Patel IR, Abramson SB: A novel mechanism of action of tetracyclines: effects on nitric oxide synthases. Proc Natl Acad Sci U S A 1996, 93: 14014–14029.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  17. Aronson AL: Pharmacotherapeutics of the newer tetracyclines. J Am Vet Med Assoc 1980, 176: 1061–1068.

    CAS  PubMed  Google Scholar 

  18. Colovic M, Caccia S: Liquid chromatographic determination of minocycline in brain-to-plasma distribution studies in the rat. J Chromatogr B Analyt Technol Biomed Life Sci 2003, 791: 337–343.

    Article  CAS  PubMed  Google Scholar 

  19. Tikka T, Fiebich BL, Goldsteins G, Keinanen R, Koistinaho J: Minocycline, a tetracycline derivative, is neuroprotective against excitotoxicity by inhibiting activation and proliferation of microglia. J Neurosci 2001, 21: 2580–2588.

    CAS  PubMed  Google Scholar 

  20. Tikka TM, Koistinaho JE: Minocycline provides neuroprotection against N-methyl-D-aspartate neurotoxicity by inhibiting microglia. J Immunol 2001, 166: 7527–7533.

    Article  CAS  PubMed  Google Scholar 

  21. Wu DC, Jackson-Lewis V, Vila M, Tieu K, Teismann P, Vadseth C, Choi DK, Ischiropoulos H, Przedborski S: Blockade of microglial activation is neuroprotective in the 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine mouse model of Parkinson disease. J Neurosci 2002, 22: 1763–1771.

    CAS  PubMed  Google Scholar 

  22. Yong VW, Wells J, Giuliani F, Casha S, Power C, Metz LM: The promise of minocycline in neurology. Lancet Neurol 2004, 3: 744–751.

    Article  PubMed  Google Scholar 

  23. Blum D, Chtarto A, Tenenbaum L, Brotchi J, Levivier M: Clinical potential of minocycline for neurodegenerative disorders. Neurobiol Dis 2004, 17: 359–366.

    Article  CAS  PubMed  Google Scholar 

  24. Stirling DP, Khodarahmi K, Liu J, McPhail LT, McBride CB, Steeves JD, Ramer MS, Tetzlaff W: Minocycline treatment reduces delayed oligodendrocyte death, attenuates axonal dieback, and improves functional outcome after spinal cord injury. J Neurosci 2004, 24: 2182–2190.

    Article  CAS  PubMed  Google Scholar 

  25. Yrjanheikki J, Keinänen R, Pellikka M, Hökfelt T, Koistinaho J: Tetracyclines inhibit microglial activation and are neuroprotective in global brain ischemia. Proc Natl Acad Sci U S A 1998, 95: 15769–15774.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  26. Yrjanheikki J, Tikka T, Keinänen R, Goldsteins G, Chan PH, Koistinaho J: A tetracycline derivative, minocycline, reduces inflammation and protects against focal cerebral ischemia with a wide therapeutic window. Proc Natl Acad Sci U S A 1999, 96: 13496–13500.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  27. Ledeboer A, Gamanos M, Lai W, Martin D, Maier SF, Watkins LR, Quan N: Involvement of spinal cord nuclear factor kappaB activation in rat models of proinflammatory cytokine-mediated pain facilitation. Eur J Neurosci 2005, 22: 1977–1986.

    Article  PubMed  Google Scholar 

  28. Mishra MK, Dutta K, Saheb SK, Basu A: Understanding the molecular mechanism of blood–brain barrier damage in an experimental model of Japanese encephalitis: correlation with minocycline administration as a therapeutic agent. Neurochem Int 2009, 55: 717–723.

    Article  CAS  PubMed  Google Scholar 

  29. Speth C, Dierich MP, Gasque P: Neuroinvasion by pathogens: a key role of the complement system. Mol Immunol 2002, 38: 669–679.

    Article  CAS  PubMed  Google Scholar 

  30. Rodriguez Parkitna J, Korostynski M, Kaminska-Chowaniec D, Obara I, Mika J, Przewlocka B, Przewlocki R: Comparison of gene expression profiles in neuropathic and inflammatory pain. J Physiol Pharmacol 2006, 57: 401–414.

    CAS  PubMed  Google Scholar 

  31. Vega-Avelaira D, Géranton SM, Fitzgerald M: Differential regulation of immune responses and macrophage/neuron interactions in the dorsal root ganglion in youngand adult rats following nerve injury. Mol Pain 2009, 5: 70.

    Article  PubMed Central  PubMed  Google Scholar 

  32. Treede RD, Jensen TS, Campbell JN, Cruccu G, Dostrovsky JO, Griffin JW, Hansson P, Hughes R, Nurmikko T, Serra J: Neuropathic pain: redefinition and a grading system for clinical and research purposes. Neurology 2008, 70: 1630–1635.

    Article  CAS  PubMed  Google Scholar 

  33. Beniczky S, Tajti J, Tímea Varga E, Vécsei L: Evidence-based pharmacological treatment of neuropathic pain syndromes. J Neural Transm 2005, 112: 735–749.

    Article  CAS  PubMed  Google Scholar 

  34. Backonja MM, Miletic G, Miletic V: The effect of continuous morphine analgesia on chronic thermal hyperalgesia due to sciatic constriction injury in rats. Neurosci Lett 1995, 196: 61–64.

    Article  CAS  PubMed  Google Scholar 

  35. Hunter JC, Gogas KR, Hedley LR, Jacobson LO, Kassotakis L, Thompson J, Fontana DJ: The effect of novel anti-epileptic drugs in rat experimental models of acute and chronic pain. Eur J Pharmacol 1997, 324: 153–160.

    Article  CAS  PubMed  Google Scholar 

  36. Makuch W, Mika J, Rojewska E, Zychowska M, Przewlocka B: Effects of selective and non-selective inhibitors of nitric oxide synthase on morphine- and endomorphin-1-induced analgesia in acute and neuropathic pain in rats. Neuropharmacology 2013, 75: 445–457.

    Article  CAS  PubMed  Google Scholar 

  37. Mika J, Rojewska E, Makuch W, Korostynski M, Luvisetto S, Marinelli S, Pavone F, Przewlocka B: The effect of botulinum neurotoxin A on sciatic nerve injury-induced neuroimmunological changes in rat dorsal root ganglia and spinal cord. Neuroscience 2011, 175: 358–366.

    Article  CAS  PubMed  Google Scholar 

  38. Obara I, Parkitna JR, Korostynski M, Makuch W, Kaminska D, Przewlocka B, Przewlocki R: Local peripheral opioid effects and expression of opioid genes in the spinal cord and dorsal root ganglia in neuropathic and inflammatory pain. Pain 2009, 141: 283–291.

    Article  CAS  PubMed  Google Scholar 

  39. Watkins LR, Milligan ED, Maier SF: Spinal cord glia: new players in pain. Pain 2001, 93: 201–205.

    Article  CAS  PubMed  Google Scholar 

  40. Kreutzberg GW: Microglia: a sensor for pathological events in the CNS. Trends Neurosci 1996, 19: 312–318.

    Article  CAS  PubMed  Google Scholar 

  41. Marchand F, Perretti M, McMahon SB: Role of the immune system in chronic pain. Nat Rev Neurosci 2005, 6: 521–532.

    Article  CAS  PubMed  Google Scholar 

  42. Schonbeck U, Libby P: The CD40/CD154 receptor/ligand dyad. Cell Mol Life Sci 2001, 58: 4–43.

    Article  CAS  PubMed  Google Scholar 

  43. Garlichs CD, Kozina S, Fateh-Moghadam S, Tomandl B, Stumpf C, Eskafi S, Raaz D, Schmeier A, Yilmaz A, Ludwig J, Neundorfer B, Daniel WG: Upregulation of CD40–CD40 ligand (CD154) in patients with acute cerebral ischemia. Stroke 2003, 34: 1412–1418.

    Article  CAS  PubMed  Google Scholar 

  44. Tan J, Town T, Paris D, Mori T, Suo Z, Crawford F, Mattson MP, Flavell RA, Mullan M: Microglial activation resulting from CD40–CD40L interaction after b-amyloid stimulation. Science 1999, 286: 2352–2355.

    Article  CAS  PubMed  Google Scholar 

  45. Gerritse K, Laman JD, Noelle RJ, Aruffo A, Ledbetter JA, Boersma WJA, Claassen E: CD40–CD40 ligand interactions in experimental allergic encephalomyelitis and multiple sclerosis. Proc Natl Acad Sci U S A 1996, 93: 2499–2504.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  46. Becher B, Durell BG, Miga AV, Hickey WF, Noelle RJ: The clinical course of experimental autoimmune encephalomyelitis and inflammation is controlled by the expression of CD40 within the central nervous system. J Exp Med 2001, 193: 967–974.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  47. Grewal I, Foellmer H, Grewal K, Xu J, Hardardottir F, Baron J, Janeway C, Flavell R: Requirement for CD40 ligand in costimulation induction, T cell activation, and experimental allergic encephalomyelitis. Science 1996, 273: 1864–1867.

    Article  CAS  PubMed  Google Scholar 

  48. Nguyen V, Walker WS, Benveniste EN: Post-transcriptional inhibition of CD40 gene expression in microglia by transforming growth factor-b. Eur J Immunol 1998, 28: 2537–2548.

    Article  CAS  PubMed  Google Scholar 

  49. Nguyen VT, Benveniste EN: Critical role of TNF-a and NFkB in IFN-c-induced CD40 expression in microglia/macrophages. J Biol Chem 2002, 277: 13796–13803.

    Article  CAS  PubMed  Google Scholar 

  50. Cao L, Beaulac H, Eurich A: Differential lumbar spinal cord responses among wild type, CD4 knockout, and CD40 knockout mice in spinal nerve L5 transection-induced neuropathic pain. Mol Pain 2012, 8: 88.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  51. Brown G: Dectin-1: a signalling non-TLR pattern-recognition receptor. Nat Rev Immunol 2006, 6: 33–43.

    Article  CAS  PubMed  Google Scholar 

  52. Salazar-Aldrete C, Galán-Díez M, Fernández-Ruiz E, Niño-Moreno P, Estrada-Capetillo L, Abud-Mendoza C, Layseca-Espinosa E, Baranda L, González-Amaro R: Expression and function of dectin-1 is defective in monocytes from patients with systemic lupus erythematosus and rheumatoid arthritis. J Clin Immunol 2013, 33: 368–377.

    Article  CAS  PubMed  Google Scholar 

  53. Viriyakosol S, Jimenez Mdel P, Gurney MA, Ashbaugh ME, Fierer J: Dectin-1 is required for resistance to coccidioidomycosis in mice. MBio 2013, 4: e00597-e00612.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  54. Burns MB, Lackey L, Carpenter MA, Rathore A, Land AM, Leonard B, Refsland EW, Kotandeniya D, Tretyakova N, Nikas JB, Yee D, Temiz NA, Donohue DE, McDougle RM, Brown WL, Law EK, Harris RS: APOBEC3B is an enzymatic source of mutation in breast cancer. Nature 2013, 494: 366–370.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  55. Lada AG, Dhar A, Boissy RJ, Hirano M, Rubel AA, Rogozin IB, Pavlov YI: AID/APOBEC cytosine deaminase induces genome-wide kataegis. Biol Direct 2012, 18: 7–47.

    Google Scholar 

  56. Roberts SA, Lawrence MS, Klimczak LJ, Grimm SA, Fargo D, Stojanov P, Kiezun A, Kryukov GV, Carter SL, Saksena G, Harris S, Shah RR, Resnick MA, Getz G, Gordenin DA: An APOBEC cytidine deaminase mutagenesis pattern is widespread in human cancers. Nat Genet 2013, 45: 970–976.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  57. Taylor BJ, Nik-Zainal S, Wu YL, Stebbings LA, Raine K, Campbell PJ, Rada C, Stratton MR, Neuberger MS: DNA deaminases induce break-associated mutation showers with implication of APOBEC3B and 3A in breast cancer kataegis. Elife 2013, 16: e00534.

    Google Scholar 

  58. Kuong KJ, Loeb LA: APOBEC3B mutagenesis in cancer. Nat Genet 2013, 45: 964–965.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  59. He L, Vasiliou K, Nebert DW: Analysis and update of the human solute carrier (SLC) gene superfamily. Hum Genomics 2009, 3: 195–206.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  60. Agnelli L, Forcato M, Ferrari F, Tuana G, Todoerti K, Walker BA, Morgan GJ, Lombardi L, Bicciato S, Neri A: The reconstruction of transcriptional networks reveals critical genes with implications for clinical outcome of multiple myeloma. Clin Cancer Res 2011, 17: 7402–7412.

    Article  CAS  PubMed  Google Scholar 

  61. Fan S, Zhao Y, Li X, Du Y, Wang J, Song X, Zhou F, Chen H, Chen G, Zhao Y, Mao Y, Lan Q: Genetic variants in SLC7A7 are associated with risk of glioma in a Chinese population. Exp Biol Med 2013, 238: 1075–1081.

    Article  Google Scholar 

  62. Xie L, Song X, Yu J, Guo W, Wei L, Liu Y, Wang X: Solute carrier protein family may involve in radiation-induced radioresistance of non-small cell lung cancer. J Cancer Res Clin Oncol 2011, 137: 1739–1747.

    Article  CAS  PubMed  Google Scholar 

  63. Borsani G, Bassi MT, Sperandeo MP, De Grandi A, Buoninconti A, Riboni M, Manzoni M, Incerti B, Pepe A, Andria G, Ballabio A, Sebastio G: SLC7A7, encoding a putative permease-related protein, is mutated in patients with lysinuric protein intolerance. Nat Genet 1999, 21: 297–301.

    Article  CAS  PubMed  Google Scholar 

  64. Barilli A, Rotoli BM, Visigalli R, Bussolati O, Gazzola GC, Gatti R, Dionisi-Vici C, Martinelli D, Goffredo BM, Font-Llitjós M, Mariani F, Luisetti M, Dall’Asta V: Impaired phagocytosis in macrophages from patients affected by lysinuric protein intolerance. Mol Genet Metab 2012, 105: 585–589.

    Article  CAS  PubMed  Google Scholar 

  65. Bonaldo MF, Lennon G, Soares MB: Normalization and subtraction: two approaches to facilitate gene discovery. Genome Res 1996, 6: 791–806.

    Article  CAS  PubMed  Google Scholar 

  66. Huang CC, Tu SH, Lien HH, Jeng JY, Huang CS, Huang CJ, Lai LC, Chuang EY: Concurrent gene signatures for han chinese breast cancers. PLoS One 2013, 8: e76421.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  67. Lee KM, Kang BS, Lee HL, Son SJ, Hwang SH, Kim DS, Park JS, Cho HJ: Spinal NF-kappaB activation induces COX-2 upregulation and contributes to inflammatory pain hypersensitivity. Eur J Neurosci 2004, 19: 3375–3381.

    Article  PubMed  Google Scholar 

  68. Bergmann TK, Vach W, Feddersen S, Eckhoff L, Gréen H, Herrstedt J, Brosen K: GWAS-based association between RWDD3 and TECTA variants and paclitaxel induced neuropathy could not be confirmed in Scandinavian ovarian cancer patients. Acta Oncol 2013, 52: 871–874.

    Article  CAS  PubMed  Google Scholar 

  69. Schneider B, Li K, Miller K, Flockhart DA, Radovich M, Han Cock BA: Genetic associations with taxane-induced neuropathy by a genome-wide association study (GWAS) in E5103 [abstracts 1000]. ASCO Meeting 2011. http://meetinglibrary.asco.org/content/78499–102

    Google Scholar 

  70. Carter C, Dion C, Schnell S, Coadwell WJ, Graham M, Hepburn L, Morgan G, Hutchings A, Pascall JC, Jacobs H, Miller JR, Butcher GW: A natural hypomorphic variant of the apoptosis regulator Gimap4/IAN1. J Immunol 2007, 179: 1784–1795.

    Article  CAS  PubMed  Google Scholar 

  71. Dalberg U, Markholst H, Hornum L: Both Gimap5 and the diabetogenic BBDP allele of Gimap5 induce apoptosis in T cells. Int Immunol 2007, 19: 447–453.

    Article  CAS  PubMed  Google Scholar 

  72. Nitta T, Nasreen M, Seike T, Goji A, Ohigashi I, Miyazaki T, Ohta T, Kanno M, Takahama Y: IAN family critically regulates survival and development of T lymphocytes. PLoS Biol 2006, 4: e103.

    Article  PubMed Central  PubMed  Google Scholar 

  73. Nitta T, Takahama Y: The lymphocyte guard-IANs: regulation of lymphocyte survival by IAN/GIMAP family proteins. Trends Immunol 2007, 28: 58–65.

    Article  CAS  PubMed  Google Scholar 

  74. Filen JJ, Filén S, Moulder R, Tuomela S, Ahlfors H, West A, Kouvonen P, Kantola S, Björkman M, Katajamaa M, Rasool O, Nyman TA, Lahesmaa R: Quantitative proteomics reveals GIMAP family proteins 1 and 4 to be differentially regulated during human T helper cell differentiation. Mol Cell Proteomics 2009, 8: 32–44.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  75. Pandarpurkar M, Wilson-Fritch L, Corvera S, Markholst H, Hornum L, Greiner DL, Mordes JP, Rossini AA, Bortell R: Ian4 is required for mitochondrial integrity and T cell survival. Proc Natl Acad Sci U S A 2003, 100: 10382–10387.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  76. Aksoylar HI, Lampe K, Barnes MJ, Plas DR, Hoebe K: Loss of immunological tolerance in Gimap5-deficient mice is associated with loss of Foxo in CD4+ T cells. J Immunol 2012, 188: 146–154.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  77. Chen XL, Serrano D, Mayhue M, Wieden HJ, Stankova J, Boulay G, Ilangumaran S, Ramanathan S: GTPase of the immune-associated nucleotide-binding protein 5 (GIMAP5) regulates calcium influx in T-lymphocytes by promoting mitochondrial calcium accumulation. Biochem J 2013, 449: 353–364.

    Article  CAS  PubMed  Google Scholar 

  78. Moralejo DH, Fuller JM, Rutledge EA, Van Yserloo B, Ettinger RA, Jensen R, Osborne W, Kwitek A, Lernmark A: BB rat Gimap gene expression in sorted lymphoid T and B cells. Life Sci 2011, 89: 748–754.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  79. Chen Y, Yu M, Dai X, Zogg M, Wen R, Weiler H, Wang D: Critical role for Gimap5 in the survival of mouse hematopoietic stem and progenitor cells. J Exp Med 2011, 208: 923–935.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  80. Zimmermann M: Ethical guidelines for investigations of experimental pain in conscious animals. Pain 1983, 16: 109–110.

    Article  CAS  PubMed  Google Scholar 

  81. Bennett GJ, Xie YK: A peripheral mononeuropathy in rat that produces disorders of pain sensation like those seen in man. Pain 1988, 33: 87–107.

    Article  CAS  PubMed  Google Scholar 

  82. Sorkin L, Svensson CI, Jones-Cordero TL, Hefferan MP, Campana WM: Spinal p38 mitogen-activated protein kinase mediates allodynia induced by first-degree burn in the rat. J Neurosci Res 2009, 87: 948–955.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  83. Chomczynski P, Sacchi N: Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem 1987, 162: 156–159.

    Article  CAS  PubMed  Google Scholar 

  84. Korostynski M, Piechota M, Dzbek J, Mlynarski W, Szklarczyk K, Ziolkowska B, Przewlocki R: Novel drug-regulated transcriptional networks in brain reveal pharmacological properties of psychotropic drugs. BMC Genomics 2013, 14: 606.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  85. Piechota M, Korostynski M, Solecki W, Gieryk A, Slezak M, Bilecki W, Ziolkowska B, Kostrzewa E, Cymerman I, Swiech L, Jaworski J, Przewlocki R: The dissection of transcriptional modules regulated by various drugs of abuse in the mouse striatum. Genome Biol 2010, 11: R48.

    Article  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by MNiSzW/NCN grant 2011/03/B/NZ4/00042 and statutory funds.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Joanna Mika.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

ER - participated in the design of the study, behavioral and qPCR studies, analysis and interpretation of data, drafting of the manuscript. MK - microarray studies and analysis, data interpretation and discussion of results. RP - discussion of microarrays results, decision on the form of the manuscript BP - conception and design, develop and decision on the content of the manuscript JM - conception and design, supervised the experiments, drafting of the manuscript, elaboration and discussion of the whole manuscript. All authors have read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( https://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rojewska, E., Korostynski, M., Przewlocki, R. et al. Expression profiling of genes modulated by minocycline in a rat model of neuropathic pain. Mol Pain 10, 47 (2014). https://doi.org/10.1186/1744-8069-10-47

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1744-8069-10-47

Keywords