Abstract
Neuropathic pain is a common neurobiological disease involving multifaceted maladaptations ranging from gene modulation to synaptic dysfunction, but the interactions between synaptic dysfunction and the genes that are involved in persistent pain remain elusive. In the present study, we found that neuropathic pain induced by the chemotherapeutic drug paclitaxel or L5 ventral root transection significantly impaired the function of GABAergic synapses of spinal dorsal horn neurons via the reduction of the GAD67 expression. We also found that mir-500 expression was significantly increased and involved in the modulation of GAD67 expression via targeting the specific site of Gad1 gene in the dorsal horn. In addition, knock-out of mir-500 or using mir-500 antagomir rescued the GABAergic synapses in the spinal dorsal horn neurons and attenuated the sensitized pain behavior in the rats with neuropathic pain. To our knowledge, this is the first study to investigate the function significance and the underlying molecular mechanisms of mir-500 in the process of neuropathic pain, which sheds light on the development of novel therapeutic options for neuropathic pain.
SIGNIFICANCE STATEMENT Neuropathic pain is a common neurobiological disease involving multifaceted maladaptations ranging from gene modulation to synaptic dysfunction, but the underlying molecular mechanisms remain elusive. The present study illustrates for the first time a mir-500-mediated mechanism underlying spinal GABAergic dysfunction and sensitized pain behavior in neuropathic pain induced by the chemotherapeutic drug paclitaxel or L5 ventral root transection, which sheds light on the development of novel therapeutic options for neuropathic pain.
Introduction
MicroRNAs (miRNAs) are an evolutionarily conserved class of small, noncoding RNAs that is involved in diverse biological processes including cellular proliferation and differentiation as negative regulators of gene expression (Filipowicz et al., 2008). Generally, miRNAs inhibit the translation of specific genes in a sequence-dependent manner via recognition of the 3′-UTR sequence of the target messenger RNAs (Bartel, 2009). Recently emerging studies suggested the potential role of miRNAs in the development and maintenance of hyperalgesia and/or allodynia induced by peripheral inflammation or nerve injury (Bai et al., 2007; Genda et al., 2013; Sengupta et al., 2013; Elramah et al., 2014), whereas the intracellular processes largely remained elusive. For instance, reduction of selective mature miRNAs (e.g., mir-124) in sensory neurons was reported in the rats with inflammatory muscle pain (Bai et al., 2007) or sciatic nerve injury (Wu et al., 2011), whereas intrathecal injection of mir-124 had an obvious anti-nociceptive effect in rodents with inflammation or nerve injury (Willemen et al., 2012). Interestingly, a recent study showed an upregulation of mir-500 in the bladder smooth muscle cells in patients with bladder pain syndrome (Sanchez Freire et al., 2010). Currently, whether mir-500 is involved in the neuropathic pain induced by chemotherapeutic drugs or nerve injury, as well as the pathophysiological involvement of mir-500 in the setting of neuropathic pain, remains unclear.
Processing of nociceptive information in the spinal dorsal horn is modulated substantially by the local GABAergic and glycinergic inhibitory circuits. It was reported that intrathecal treatment of the GABA receptor antagonist produced remarkable hypersensitivity to innocuous tactile stimuli (Malan et al., 2002; Gwak et al., 2006) and transgenic mice lacking specific subunits of the GABA receptor developed hyperalgesia and allodynia (Schuler et al., 2001). GABA is synthesized from glutamate by two isoforms of glutamate decarboxylase (GAD), GAD67 and GAD65, which are encoded by two separate genes, Gad1 and Gad2, respectively. GAD67 is distributed throughout the cell, whereas GAD65 localizes preferentially in synaptic terminals (Esclapez et al., 1994). Studies have shown that GAD67 was responsible for >90% of basal GABA synthesis (Asada et al., 1997) and knock-down of Gad1, but not Gad2, resulted in cell-autonomous deficits in axon branching, perisomatic synapse formation, and complexity of the innervation fields (Chattopadhyaya et al., 2007), suggesting a specific role of GAD67 in maintining proper neural function. Furthermore, depletion of Gad1 markedly reduced the amplitude and frequency of miniature IPSCs (mIPSCs) in central neurons (Lau and Murthy, 2012) and chronic constriction of the sciatic nerve significantly decreased the expression of GAD67 (Vaysse et al., 2011). These results implied the potential involvement of GAD67 downregulation in the adaptation of spinal nociceptive circuits in the setting of neuropathic pain. Currently, whether GAD67 contributes to neuropathic pain and the underlying mechanism remained unclear. The present study focused on the mechanism of mir-500-mediated modulation of GAD67 and its functional significance in rodent models with neuropathic pain induced by paclitaxel or L5 ventral root transection (L5-VRT).
Materials and Methods
Animals.
Male Sprague Dawley rats (180–230 g) were obtained from Sun Yat-Sen University Animal Care and Use Committee. To examine the role of mir-500 in neuropathic pain, we create rno-mir-500 knock-out (KO) Sprague Dawley rats by microinjection of TALENs in fertilized eggs (Cyagen Biosciences). Briefly, the rno-mir-500 gene (Gene ID: 100314112; MirBase accession no. MIMAT0005321), which was located on rat chromosome X, was selected as the TALEN target site. TALENs were constructed using the Golden Gate Assembly and confirmed by sequencing. TALEN mRNA, generated by in vitro transcription, was injected into fertilized eggs for KO rat production. The founders, which were missing 13 bases in the both strands compared with the wild-type (WT) DNA sequence, were genotyped by PCR followed by DNA-sequencing analysis. After that, heterozygous F1 were produced by mating F0 with WT rats. Heterozygous F1 rats were determined and mated to produce F2. Homozygote F2 (mir-500−/−) and WT (mir-500+/+) was determined by DNA-sequencing analysis (see Fig. 5A). All animals were housed in a temperature-controlled room (22 ± 1°C) with a 12 h light/dark cycle. All animal experimental procedures were approved by local Animal Care and Use Committee and performed in accordance with the animal care and ethical guidelines of National Institutes of Health.
Reagents and L5-VRT.
Paclitaxel (Taxol; Bristol-Myers Squibb, 6 mg/ml) was diluted with saline (1:3) and injected (8 mg/kg, i.p., cumulative dose of 24 mg/kg) on 3 alternate days (days 1, 4, and 7). The cholesterol-conjugated mir-500 antagomir and antagomir control were obtained commercially from Ribobio for the in vivo experiments. mir-500 mimic for the in vitro experiments was also purchased from Ribobio.
L5-VRT model was performed as described previously (Li et al., 2002; Zang et al., 2015). Briefly, surgery was performed on rats under deep anesthesia after intraperitoneal administration of sodium pentobarbital (50 mg/kg). An L5 hemilaminectomy was performed to expose the L5 nerve root. The ventral root was gently pulled out with fine forceps and transected 2–3 mm proximal to the DRG and a small portion (2 mm) of the root was dissected. In the sham group, an identical operation was performed for exposure of the L5 ventral root, but the nerve was not transected.
Intrathecal application of miRNA modulators and behavioral test.
Intrathecal injection of mir-500 antagomir (0.5 nmol/d) was performed according to previously described methods (Li et al., 2015). In brief, a polyethylene-10 catheter was inserted into the rat's subarachnoid space through L5 to L6 spinal cord segmental intervertebral space and the tip of the catheter was located at the L5 spinal segmental level. mir-500 antagomir was intrathecal injected 30 min before the first dose of paclitaxel or L5-VRT surgery and continued for 5 d.
The withdrawal threshold of the hindpaw was determined by applying mechanical stimuli to the plantar surface of the hindpaw using Von Frey hairs. The 50% withdrawal threshold was defined as the lowest force that produced five or more responses of 10 stimuli (Liu et al., 2010).
Thermal hyperalgesia was tested using a plantar test (7370; Ugo Basile Plantar Test Apparatus) as described previously (Hargreaves et al., 1988). Briefly, a radiant heat source beneath a glass floor was aimed at the plantar surface of the hindpaw. Three measurements of hindpaw withdrawal latency were taken for each hindpaw in each test session. The hindpaw was tested alternately with >5 min intervals between the consecutive tests. Three measurements of latency per side were averaged as the result of each test. A 20 s cutoff was set to prevent the tissue damage. The experimenter who conducted the behavioral tests was blinded to all treatments.
Application of miRNA modulators in glioma cell line C6.
Malignant glioma cell line C6 was cultured in regular medium (DMEM with 10% fetal bovine serum and 1% antibiotics). Then, 20 pmol of mir-500 mimic and scrambled control miRNA (Ribo) were transfected into C6 cells in 24-well plates using Lipofectamine 2000 (Invitrogen) according to the manufacturer's instructions. The expression of GAD67 was examined using the Western blot method 1 d after transfection.
RNA isolation and microarray experiments.
Total RNA was extracted from spinal dorsal horn tissue using TRIzol reagent (Invitrogen). RNA labeling and hybridization on miRNA microarray chips were conducted as described previously (Liu et al., 2004a). Briefly, the purified RNA using the mirVANA miRNA Isolation Kit (Applied Biosystems) was labeled with fluorescein and hybridization was performed on miRNA microarray chips (miRNA microarray V4.0; CapitalBio) containing 1965 probes in triplicate. These probes corresponded to 988 human, 627 mice, and 350 rat miRNA genes designed based on the miRBase Release 12.0. Three independent RNA samples on different time points from spinal dorsal horn tissue in the rats treated with paclitaxel or vehicle were hybridized with miRNA microarrays separately. Hybridization intensity values from individual samples were filtered and normalized to per-chip mean values. We considered candidate miRNAs with a signal >500 as positive.
Western blot.
Western blot analyses were performed as described previously (Huang et al., 2014; Li et al., 2015). Briefly, lysates of spinal dorsal horn tissues or cells were prepared and separated by gel electrophoresis (SDS-PAGE) and transferred onto a polyvinylidene fluoride membrane. The blots were placed in blocking buffer for 1 h at room temperature and incubated with primary antibody against GAD67 (1:1000; Abcam) overnight at 4°C. The blots were then incubated with horseradish peroxidase-conjugated IgG. Electrochemiluminescence (Pierce) was used to detect the immune complex. The band was quantified with computer-assisted imaging analysis system (ImageJ).
Plasmid construction and luciferase reporter assay.
A luciferase reporter assay was performed to test the binding of mir-500 to its target gene Gad1. A 1233 bp segment of rat Gad1 3′-UTR containing a presumed mir-500 complementary site (seed sequence, GTGCAT) was amplified by PCR using rat cDNA as a template. The PCR products were inserted into the psi-CHECK2 plasmid (Applied Biosystems) and efficient insertion was confirmed by sequencing. To test the binding specificity, the seed sequence of mir-500 was mutated from GTGCAT to CACGTA (see Fig. 2H). For the luciferase reporter assays, 0.5 μg of firefly luciferase reporter plasmid and 20 pmol of mir-500 mimic and scrambled control miRNA (Ribo) were transfected into 293T cells in 24-well plates using Lipofectamine 2000 (Invitrogen) according to the manufacturer's instruction. The b-gal vector was used as a transfection control. One day after transfection, the cells were assayed using the Dual-Luciferase Reporter Assay System (Promega). Briefly, cells were lysed with 35 μl of passive lysis buffer per well. The culture plate was incubated for 15 min with shaking at 400 rpm at room temperature. To record luminescence, 100 μl of Luciferase Assay Reagent II was added to 30 μl of lysate for first measurement (firefly luciferase) and 100 μl of Stop & Glo Reagent was added for the second measurement (Renilla luciferase). The Renilla/firefly luciferase ratio was used for statistical analysis. Three to four replicate wells were used in each experiment
FISH and immunohistochemistry.
Animals were anesthetized with sodium pentobarbital (50 mg/kg, i.p.) and perfused transcardially with 4% paraformaldehyde in 0.1 m DEPC-PBS. The lumbar enlargement of spinal cord was removed, postfixed with the same fixative for 24 h at 4°C, and paraffin embedded. FISH analysis was performed according to the manufacturer's guidelines. Briefly, the spinal cord tissues were cut into 10-μm-thick transverse sections and washed with 75% alcohol. The sections were then hybridized at 42°C for 24 h with 5′-digoxigenin-labeled mir-500 oligonucleotide probe (1:800, EXQON). Sections were washed and incubated with blocking buffer for 60 min at room temperature. Sections were incubated at 4°C overnight with anti-digoxigenin-POD polyclonal antibody (1:500; Roche) and GAD67 antibody (1:200; Santa Cruz Biotechnology). Sections were then incubated with TSA Plus Fluorescein (1:250; PerkinElmer) and Alexa Fluor 555-conjugated secondary antibody (1:500; Thermo Fisher) at 37°C for 30 min. The stained sections were then examined with a Zeiss confocal microscope equipped with a digital camera (Zeiss). The quantification of immunofluorescence staining was performed by counting the number of mir-500 and GAD67-immunoreactive positive cells in each section. The proportion of mir-500-expressing GABAergic neurons was determined as the number of mir-500- and GAD67-positive cells divided by the total number of GAD67-positive cells. In each rat, every fourth section was picked from a series of consecutive sections and four sections were counted. An average proportion of mir-500-expressing GABAergic neurons were obtained.
For GAD67 staining, sections (16 μm) were cut and processed for immunohistochemistry with primary antibody for GAD67 (1:200; Santa Cruz Biotechnology). After incubation overnight at 4°C, the sections were incubated with fluorescein isothiocyanate-conjugated secondary antibodies for 1 h at room temperature. The stained sections were then examined with a Leica fluorescence microscope and images were captured with a Leica DFC350 FX camera. For quantification of GAD67 immunostaining, the immunoreactive-positive area was analyzed with a Leica Qwin V3 image system.
PCR amplification of miRNAs and mRNAs.
Stem-loop qRT-PCR assays using TaqMan miRNA probes (Applied Biosystems) to quantify mature miRNAs were performed, as previously reported (Chen et al., 2005; Li et al., 2009). Real-time PCR was performed using a TaqMan PCR kit and the Applied Biosystems 7500 Sequence Detection System with the following conditions: 95°C for 10 min, 60 cycles of 95°C, 15 s, and 60°C 1 min. All reactions, including no-template controls, were run in triplicate. The relative expression ratio of miRNA in the rat's spinal tissues was quantified using the 2−δδCT method.
To detect the mRNA of GAD67, real-time qRT-PCR was performed using SYBR Green qPCR SuperMix (Invitrogen) and the ABI PRISM7500 Sequence Detection System. The reactions were set up based on the manufacturer's protocol. PCR conditions were as follows: incubation at 95°C for 3 min, followed by 40 cycles of thermal cycling (10 s at 95°C, 20 s at 58°C, and 10 s at 72°C). All primers used are listed in Table 1.
Specific primer sequences
Spinal cord slice preparation and patch-clamp recording.
Spinal cord slices were prepared and whole-cell patch-clamp recording was performed as described previously with minor adjustments (Moore et al., 2002; Huang et al., 2014). Briefly, male Sprague Dawley rats (age 25–35 d) were deeply anesthetized with sodium pentobarbital (50 mg/kg, i.p.) and the spinal cord was quickly removed into ice-cold incubation solution consisting of the following (in mm): 95 NaCl, 1.8 KCl, 1.2 KH2PO4, 0.5 CaCl2, 7 MgSO4, 26 NaHCO3, 15 glucose, and 50 sucrose, oxygenated with 95% O2/5% CO2, pH 7.4, at a measured osmolarity of 310–320 mOsm/L. The dura mater and the ventral roots were removed. Parasagittal slices with attached dorsal roots (10–15 mm long) were cut at 400–600 μm thickness using a vibrating microslicer (DTK-1000; Dosaka) and kept in the incubation solution at 33°C for at least 60 min. After incubation, an individual slice with attached dorsal root was transferred to a recording chamber and continually perfused with oxygenated ACSF solution at room temperature. The dorsal root was gently led into the suction electrode. Lamina I-II neurons were visualized using a 60× water-immersion objective on an upright infrared Nikon microscope. EPC-10 amplifier and the PULSE program (HEKA Electronics) were used with a pipette (4–6 MΩ) with internal solution containing the following (in mm): 120 CsSO4, 20 TEA-Cl, 2 MgCl2, 2 Na2ATP, 10 EGTA, 0.5 NaGTP, biocytin 5, 10 HEPES, pH 7.28, with CsOH at a measured osmolarity of 300 mOsm/L. Stimulation was applied to the attached dorsal roots at 0.1 ms pulse width via a suction electrode with an isolated current stimulator (A360; World Precision Instruments). The neurons were voltage clamped at a holding potential of 0 mV, thus minimizing the contribution of NMDA and AMPA/kainite receptor-mediated events (Yoshimura and Nishi, 1993; Baccei and Fitzgerald, 2004). After incubation with 1 μm strychnine, IPSCs were evoked by stimulating the attached dorsal root via the suction electrode. At the end of all experiments, bicuculline (10 μm) was bath applied and completely blocked IPSCs, confirming that the recorded currents were mainly mediated by GABAA receptors. For GABAergic mIPSCs, slices were incubated with antagonists of glycinergic and glutamatergic ionotropic transmission [1 μm strychnine, 50 μm d-(−)-2-amino-5-phosphonopentanoic acid, and 10 μm 6-cyano-7-nitroquinoxaline-2,3-dione] and were recorded in the presence of lidocaine (2 mm) to block action-potential-dependent effects. In all experiments, only one neuron per spinal cord slice was used.
Statistical analysis.
All data are expressed as means ± SD and were analyzed using SPSS version 13.0 software. Western blot, qRT-PCR, and electrophysiological data were analyzed by two-way ANOVA followed by Tukey's post hoc test. For behavioral analysis, one-way or two-way ANOVA with repeated-measures followed by Tukey's post hoc test was performed. The criterion for statistical significance was p < 0.05. The sample size was based on our and others' experience in painful behavior studies to provide power analysis.
Results
Decreased GABAergic synaptic function and GAD67 expression in superficial dorsal horn is induced by paclitaxel or L5-VRT
We first examined the change of GABAergic inhibitory synaptic connections in superficial spinal dorsal horn of the rats with persistent hyperalgesia and allodynia induced by paclitaxel. Paclitaxel (3 × 8 mg/kg; cumulative dose, 24 mg/kg) induced chronic marked thermal hyperalgesia (Fig. 1A) and mechanical allodynia (Fig. 1B) on days 4, 6, and 10 d, but not day 1, compared with the vehicle group. Next, we found that the slope of the input–output curve for GABAergic IPSCs was decreased in superficial dorsal horn neurons at day 4 and 10, but not day 1, after paclitaxel injection (Fig. 1C). Notably, the evoked IPSCs were completely blocked by bath application of bicuculline (Fig. 1D). Importantly, the reduced frequency of mIPSCs (Fig. 1E) was observed in the superficial dorsal horn neurons at days 4 and 10 after injection of paclitaxel. These results further confirmed that paclitaxel treatment reduced the presynaptic GABA releases in the dorsal horn and mediated the chronic pain. The study also showed that LF-VRT, another rodent model of chronic neuropathic pain (Li et al., 2002; Zang et al., 2015), significantly induced the persistent thermal hyperalgesia (Fig. 1F) and mechanical allodynia (Fig. 1G). Furthermore, similar to the results in the paclitaxel-treated rats, the slope of the input–output curve for GABAergic IPSCs (Fig. 1H) and the frequency of mIPSCs (Fig. 1I) were significantly decreased on day 5 and 15, but not day 1, after L5-VRT compared with the sham group. These results indicated a decreased presynaptic GABA release onto the superficial dorsal horn neurons and persistent pain-related behavior induced by paclitaxel or VRT.
Persistent pain decreases GABAergic synaptic function by inhibiting presynaptic GABA release. A, B, Time course of changes in withdrawal latency and withdrawal threshold after intraperitoneal paclitaxel injection. **p < 0.01 versus the vehicle group; n = 12 in each group. C, Representative traces of GABA IPSCs evoked by various stimulation intensities in spinal superficial dorsal horn neurons. *p < 0.05, **p < 0.01 on d 4 after paclitaxel group versus the corresponding vehicle group. #p < 0.05, ##p < 0.01 on d 10 after paclitaxel group versus corresponding vehicle group. n = 21 in each group. D, IPSCs were abolished by bicuculline. E, Representative traces of mIPSCs in neurons from the four groups (n = 60 cells for each group). F, G, Time course of changes in withdrawal latency and withdrawal threshold after L5-VRT treatment. **p < 0.01 versus the sham group. n = 12 in each group. H, Representative traces of GABA IPSCs evoked by various stimulation intensities in spinal superficial dorsal horn neurons. **p < 0.01 d 4 after L5-VRT group versus corresponding vehicle group. ##p < 0.01 d 10 after L5-VRT group versus corresponding vehicle group n = 21 in each group). I, Representative traces of mIPSCs in neurons from the L5-VRT or sham groups. n = 49 cells for each group.
Depletion of GAD67 led to a 93% reduction in GABA levels, which suggested the critical role of GAD67 in the synthesis of GABA neurotransmitter (Asada et al., 1997; Martin, 2000). We then determined the level of GAD67 protein and mRNA in the spinal dorsal neurons after the development of persistent pain. Although the expression of GAD67 protein was significantly decreased at days 4 and 10 (Fig. 2A), the GAD67 mRNA was not significantly changed at any time point (Fig. 2B) after paclitaxel injection. Similarly, the levels of GAD67 protein were also decreased at day 15 (Fig. 2C) and the mRNA level was not changed after L5-VRT surgery (Fig. 2D). Immunohistochemical analysis further revealed that GAD67 was downregulated throughout the superficial dorsal horn on day 10 after paclitaxel treatment or on day 15 after VRT surgery (Fig. 2E). These results suggested that a potential posttranscriptional mechanism was involved in the GAD67 downregulation in neuropathic pain induced by paclitaxel or L5-VRT.
The GAD67 protein, but not mRNA, is decreased by persistent pain. A, Representative blots and histogram showing downregulation of GAD67 induced by paclitaxel treatment. **p < 0.01 versus corresponding vehicle group. n = 6 in each group. B, Paclitaxel treatment did not affect the mRNA level of GAD67. n = 6 in each group. C, Western blots and group data of GAD67 protein in L5-VRT rats. **p < 0.01 versus corresponding sham group. n = 6 in each group. D, Analysis of GAD67 mRNA level after L5-VRT. n = 6 in each group. E, Immunohistochemical analysis revealing decreased GAD67 levels in the superficial dorsal horn after paclitaxel treatment or L5-VRT. n = 6 in each group.
Increased mir-500 regulates GAD expression through the predicted sites
Noncoding RNAs, including miRNAs, repress gene expression at the posttranscriptional level through binding, usually at the 3′-UTR of a target mRNA, thereby providing a dynamic brake on the protein synthesis (Kim et al., 2009). To identify the potential involvement of miRNAs in the pathogenesis of neuropathic pain, we collected the tissue of L4–L6 spinal dorsal horn at days 4 and 10 after paclitaxel treatment and performed genome-wide miRNA screening to detect the change in miRNAs. Because we were interested in the mechanism underlying GAD67 reduction in the spinal dorsal horn, two conditions, the potential binding site from databases of TargetScan and miRanda analysis and at least 2.5-fold upregulation of miRNAs, were set as the criteria to screen the target miRNAs. The results showed that 35 miRNAs of the total 350 were increased with a >2.5-fold change on days 4 and 10 after paclitaxel treatment (Fig. 3A) and 4 miRNAs, including mir-500, mir-188, mir-181, and mir-331, were related to GAD67 expression from the analysis of these 2 criteria (Fig. 3B). qRT-PCR analyses further showed that mir-500 was significantly increased on days 4 and 10, but not day 1, after paclitaxel treatment (Fig. 3C). However, Mir-188 and mir-331 only had a significant increase on day 10 and mir-181 did not show a significantly change after paclitaxel treatment (Fig. 3C). Increased expression of mir-500 was also observed in the spinal dorsal horn at days 5 and 15, but not day 1, after L5-VRT in the rats (Fig. 3D).
Increased mir-500 is induced by persistent pain directly regulated by GAD67 expression. A, Heat maps of miRNAs found to be significantly upregulated via microarray analysis in the spinal dorsal horn of rats on days 4 and 10 after paclitaxel treatment compared with the vehicle group. B, Database analysis showed that four miRNAs, including mir-500, mir-188, mir-181, and mir-331, are increased relative to the Gad1 gene. C, qRT-PCR analysis validated the upregulation of mir-500, mir-188, and mir-331, but not mir-181, after paclitaxel treatment. *p < 0.05, **p < 0.01 versus corresponding vehicle group. n = 6 in each group. D, Expression level of mir-500 in the spinal dorsal horn after L5-VRT. **p < 0.01 versus corresponding sham group. n = 6 in each group. E, FISH showed that mir-500 is coexpressed on GAD67-positive cells in the spinal dorsal horn. F, Representative blots and histogram showing downregulation of GAD67 induced by mir-500 mimic in a C6 cell. **p < 0.01 versus corresponding scramble group. n = 6 in each group. G, Luciferase activity of reporter constructs containing either the 3′-UTR of Gad1 or mutated Gad1 site 1, site 2, and site 3 or mutated all of Gad1 three sites, was cotransfected with mir500 mimic. ##p < 0.01 versus scramble group, *p < 0.05, **p < 0.01 versus mir500 mimic group. n = 6. H, Schematic representation of 3′-UTR of Gad1 gene demonstrating multiple mutated complementary binding sites of mir-500. I, Schematic representation of the specific mir-500 target sequence within the 3′ UTR of Gad1. The seed sequence and mutated sequence are indicated by bold letters.
To verify whether mir-500 is actually involved in GAD67 expression, we first determined the expression of mir-500 and GAD67 using the FISH method. In all examined sections, 768 cells were GAD67-positive and 576 were mir-500-positive. All (100%) mir-500-positive cells were colocalized with the GAD67-positive cells in the spinal dorsal horn (Fig. 3E). To explore the functional association between mir-500 and GAD67, we transfected the mir-500 mimic into the C6 cells, which expressed a detectable amount of endogenous GAD67. The results showed that the mir-500 mimic significantly inhibited the GAD67 expression in C6 cells 1 d after transfection (Fig. 3F). Next, we constructed a luciferase reporter vector (psi-CHECK2 vector) in which the 3′-UTR of Gad1 was fused to the luciferase coding sequence. We then transfected human embryonic kidney 293 (HEK-293) cells with luciferase reporter vector and mir-500 mimic. The results showed that mir-500 mimic at dose of 20 pmol or 50 pmol, but not 5 pmol, significantly reduced the luciferase signal compared with the scrambled control miRNA (Fig. 3G), suggesting that mir-500 directly targeted the 3′-UTR region of GAD67 mRNA to modulate its expression. Further silico analysis showed that mir-500 had three predicted binding sites in the 3′-UTR region of GAD67 mRNA, so we then generated mutant constructs for each of the sites individually (mir-500-1m, mir-500-2m, mir-500-3m) and a triple mutant with all three binding sites mutated (mir-500-m) in the 3′-UTR region of Gad1 (Fig. 3H,I). We cotransfected the mutant constructs with the corresponding mir-500 mimic into HEK-293 cells and assessed the relative extent of inhibition. Mutation of the first predicted site in the 3′-UTR region of Gad1 (mir-500-1 m) significantly abolished the suppression of luciferase signal induced by mir-500 mimic and mutation in the third site (mir-500-3m) partially reversed the suppression of luciferase activity by the mir-500 mimic, whereas mutation in the second site (mir-500-1m) did not change mir-500 mimic-induced suppression (Fig. 3G). These results suggested that 635–659 in the 3′-UTR of Gad1 was the critical site for mir-500 to regulate the expression of GAD67.
mir-500 contributes to neuropathic pain via regulation of GAD67
To assess the involvement of mir-500 in the induction of persistent pain, we blocked mir-500 function using a mir-500 antagomir. Antagomirs are synthetic miRNA inhibitors that are chemically modified for enhanced stability and cholesterol conjugated for high-efficiency delivery (Liang et al., 2014; Zhou et al., 2016). The results showed that intrathecal injection of mir-500 antagomir (0.5 nmol/d for 5 consecutive d) alone did not change the withdrawal latency to thermal stimuli (Fig. 4A) or the withdrawal threshold to mechanical stimuli in naive rats (Fig. 4B). However, compared with antagomir control, intrathecal injection of mir-500 antagomir (0.5 nmol/d for 5 consecutive d) significantly elevated the withdrawal latency (Fig. 4A) and withdrawal threshold (Fig. 4B) and increased the frequency of mIPSCs of superficial dorsal horn neurons in the rats treated with paclitaxel (Fig. 4C). qRT-PCR analyses showed that mir-500 antagomir significantly decreased the mir-500 upregulation on days 4 and 10 after paclitaxel treatment (Fig. 4D). Western blots also showed that mir-500 antagomir significantly increased the expression of GAD67 in the spinal dorsal horn on day 10 after initial treatment with paclitaxel (Fig. 4E). Similarly, mir-500 antagomir also inhibited thermal hyperalgesia (Fig. 4F) and mechanical allodynia (Fig. 4G) induced by L5-VRT and increased the frequency of mIPSCs of superficial dorsal horn neurons after L5-VRT (Fig. 4H). Notably, intrathecal injection of antagomir failed to affect the GAD67 expression, frequency of mIPSCs, and pain behavior in the control rats.
mir-500 modulates GABAergeic synapses and contributed to behavioral hypersensitivity. A, B, Intrathecal injection of mir-500 antagomir for 5 consecutive d increased the withdrawal latency to thermal stimuli and withdrawal threshold to mechanical stimuli. **p < 0.01 versus antagomir group; ##p < 0.01 versus corresponding antagomir control + paclitaxel group. n = 12 in each group. C, Intrathecal injection of mir-500 antagomir increased the frequency of mIPSCs after paclitaxel treatment. **p < 0.01 versus antagomir group; ##p < 0.01 versus corresponding antagomir control + paclitaxel group. n = 51 cells for each group. D, mir-500 antagomir significantly decreased mir-500 upregulation on days 4 and 10 after paclitaxel treatment. **p < 0.01 versus antagomir control group. n = 6 for each group. E, Representative blots and histogram showing that mir-500 antagomir inhibited paclitaxel-induced downregulation of GAD67. **p < 0.01 versus antagomir group; ##p < 0.01 versus corresponding antagomir control + paclitaxel group. n = 12 in each group. F, G, Intrathecal injection of mir-500 antagomir for 5 consecutive d attenuated thermal hyperalgesia and mechanical allodynia induced by L5-VRT. **p < 0.01 versus antagomir group; #p < 0.05, ##p < 0.01 versus corresponding antagomir control + L5-VRT group. n = 12 in each group. H, Intrathecal injection of mir-500 antagomir increased the frequency of mIPSCs in dorsal horn neurons in rats after L5-VRT. **p < 0.01 versus antagomir group; ##p < 0.01 versus corresponding antagomir control + VRT group. n = 60 cells for each group.
KO of mir-500 increases mIPSC frequency and prevents behavioral hypersensitivity induced by paclitaxel or L5-VRT
To confirm the role of mir-500 in persistent pain, we created a mir-500 KO rat using TALENs. The newborn F2 rats were genotyped by PCR followed by DNA-sequencing analysis to identify homozygotes F2 (Fig. 5A). We used homozygote F2 rats (mir-500−/−) to investigate the potential mir-500-mediated regulation of mIPSC frequency in the dorsal horn. The results showed that the frequency of mIPSCs in the KO rats was significantly increased compared with that in WT rats after paclitaxel treatment (Fig. 5B). Moreover, compared with the WT group, the KO rats exhibited a significant remission in thermal hyperalgesia (Fig. 5C) and mechanical allodynia (Fig. 5D) after paclitaxel treatment. Western blot analysis showed that GAD67 expression exhibited a significant recovery in the KO rats on day 10 after paclitaxel treatments (Fig. 5E). Similarly, the mir-500 KO rats exhibited a significant recovery of GAD67 expression (Fig. 5F) and higher withdrawal latency (Fig. 5G) and withdrawal threshold (Fig. 5H) after the L5-VRT surgery. These results further confirmed that mir-500 plays an important role in the neuropathic pain induced by paclitaxel or L5-VRT.
The KO of mir-500 increased mIPSC frequency and ameliorated the behavioral hypersensitivity induced by paclitaxel or L5-VRT. A, DNA-sequencing analysis showing that the F2 rats used for the experiment are homozygote rats for KO mir-500. The deleted bases are marked by blue. B, KO of mir-500 increased the frequency of mIPSC in dorsal horn neurons. **p < 0.01 versus corresponding WT group; ##p < 0.01 versus corresponding WT + paclitaxel group. n = 52 cells for each group. C, D, Withdrawal latency and withdrawal threshold were significantly increased in KO rats after paclitaxel treatment compared with WT rats also treated with paclitaxel. **p < 0.01 versus corresponding WT group; ##p < 0.01 versus WT + paclitaxel group. n = 12 rats for each group. E, F, Representative blots and histogram showing that the KO of mir-500 inhibited downregulation of GAD67 induced by paclitaxel or L5-VRT. **p < 0.01 versus WT group; ##p < 0.01 versus corresponding WT + paclitaxel or WT + L5-VRT group. n = 12 in each group. G, H, KO of mir-500 increased the withdrawal latency to thermal stimuli and withdrawal threshold to mechanical stimuli in rats after L5-VRT compared with those in WT rats after L5-VRT. **p < 0.01 versus corresponding WT group; #p < 0.05, ##p < 0.01 versus corresponding WT + paclitaxel or WT + L5-VRT group. n = 12 rats for each group.
Discussion
In a rat model of neuropathic pain, we first showed that the chemotherapeutic drug paclitaxel or L5-VRT substantially enhanced mir-500 expression and consequently impaired inhibitory GABAergic transmission in spinal dorsal horn neurons, contributing to the development of persistent behavioral hypersensitivity. These results were supported by evidence that rats with mir-500 KO or treated with mir-500 antagomir demonstrated the rescue of GABA synaptic function in the spinal dorsal horn neurons and the attenuated pain behavior. In addition, in vivo and in vitro studies further showed that 635–659 in 3′-UTR of gad1 might be the critical binding site for mir-500 to regulate the expression of GAD67. To our knowledge, this is the first study to investigate the functional significance of mir-500 and the underlying molecular mechanisms in rodent models of neuropathic pain, providing critical insights into how mir-500 participates in the regulation of GAD67 expression to contribute to the induction of neuropathic pain.
Role of GAD67 in the spinal dorsal horn in neuropathic pain
GAD67 is a critical enzyme for the synthesis of cytoplasmic GABA and control GABA vesicular filling and synaptic homeostasis (Lau and Murthy, 2012). Impaired synthesis and release of GABA may cause the loss of GABAergic inhibition, leading to the hyperexcitation of neurons. For instance, peripheral injury induces a loss of GABA inhibition in spinal neurons (Moore et al., 2002) and enhancing GABA synaptic inhibition alleviates injury-induced pain (Knabl et al., 2008). Our results suggested that the chemotherapeutic drug paclitaxel or L5-VRT substantially downregulates GAD67 expression, leading to the impairment of GABA synaptic inhibition in the spinal dorsal horn and sensitized pain behavior. This was consistent with the previous reports that viral delivery of Gad1 reduces neuropathic pain (Liu et al., 2004b; Liu et al., 2008). In the present study, the downregulation of GAD67 was coincidently observed with the impairment of GABA synaptic inhibition in the dorsal horn neurons of the modeled rats, which was consistent with previous studies showing that cytosolic concentration of the GABA transmitter is a critical factor in its vesicular loading (Hartman et al., 2006) and that lack of GAD67 significantly decreases the GABA level (Asada et al., 1997). Previous studies also showed that homozygous deletion of Gad1 (Gad1−/−) significantly reduces mIPSC amplitude and frequency in cultured hippocampal neurons from newborn mice pups (Lau and Murthy, 2012), whereas the frequency, but not amplitude, of mIPSCs was slightly, though insignificantly, decreased in Gad1+/− neurons (Lau and Murthy, 2012). These results implied that adaptation of mIPSCs was largely dependent on the content of GAD67. In the present study, substantial reduction of GAD67 was observed in the dorsal horn neurons, which might be responsible for the decreased mIPSC frequency and attenuated GABAergic input–output response in rodent models of neuropathic pain.
Furthermore, studies have shown that paclitaxel treatment and direct nerve damage may induce a number of neuroinflammatory responses, including increases of proinflammatory cytokines in the dorsal horn. These cytokines, including TNF-α and IL-6, suppress GABA-mediated IPSCs via TNFR1 and the IL6R (Kawasaki et al., 2008; Zhang et al., 2013). However, whether these proinflammatory cytokines suppressed GAD67 expression, thus leading to GABAergic suppression in dorsal horn and behavioral sensitization, needs to be further studied in the setting of neuropathic pain. The mechanism underlying GAD67 reduction was complex among the different tissues. For example, various factors such as transcription factors Zif268 (Yanagawa et al., 1997), histone acetylation (Tremolizzo et al., 2002), and DNA methylation (Kundakovic et al., 2007) participated in GAD67 expression via modification of the Gad1 promoter. The present study provided a novel theory of the mechanism underlying GAD67 downregulation in the spinal dorsal horn that mir-500 may affect the function of GAD67 mRNA and, consequently, mediate the reduction of GAD67 in spinal dorsal horn neurons after paclitaxel treatment or L5-VRT. Furthermore, neuropathic pain generally involves the altered expression of many genes (Lacroix-Fralish et al., 2007) and it is likely that other genes in addition to Gad1 are also the targets of persistent pain-induced chromatin remodeling. GAD65 is another key rate-limiting-rate and is targeted preferentially to presynaptic terminals of central neurons for GABA synthesis (Tian et al., 1999; Patel et al., 2006). Studies have showed that GAD65 in different regions could participate in pain regulation by decreasing inhibition among neurons through reducing the levels of cellular GABA (Gwak and Hulsebosch, 2011; Zhang et al., 2011; Lorenzo et al., 2014). For example, suppression of Gad2 transcription through histone hypoacetylation in brainstem nucleus raphe magnus results in the impairment of GABA synaptic function and sensitized pain behavior (Zhang et al., 2011). Our data showed an important role of GAD67 in the spinal dorsal horn in pain modulation, whereas this did not completely preclude a potential involvement of GAD65. We note that the detailed mechanism by which GAD65 might play a role in pain modulation needs further study.
Role of mir-500-mediated targeting of GAD67 in neuropathic pain
Several studies have revealed the dysregulation of miRNA expression in the spinal cord under neuropathic pain conditions (Genda et al., 2013; Elramah et al., 2014), whereas the intracellular processes largely remained elusive. Intrathecal injection of mir-124 had an obvious antinociceptive effect in rodents with inflammation or nerve injury (Willemen et al., 2012). Using a miRNA microarray, the present study showed that four miRNAs, including mir-500, mir-188, mir-181, and mir-331, were significantly increased in the spinal dorsal horn after paclitaxel treatment (>2.5-fold change) compared with the vehicle group. qRT-PCR analysis further confirmed that mir-500 exhibited the most significant increase and the pattern was consistent with paclitaxel-induced behavioral hypersensitivity. Similarly, an increase of mir-500 was also observed in the dorsal horn in rats after L5-VRT. These results were supported by the observation that mir-500 incurred a significant increase in the bladder smooth muscle cells of patients with bladder pain syndrome (Sanchez Freire et al., 2010). All these studies suggested that mir-500 might be a critical molecule to mediating pain syndromes. Although our present results and those of others (Arai et al., 2013) suggested that miR-188 might be involved in hyperalgesia, the function of miR-188 was not investigated in the current study, mainly because mir-500 had a much greater increase in spinal dorsal horn after paclitaxel treatment. Notably, we provided the morphological evidence to demonstrate the distribution of mir-500 in GABAergic neurons in the dorsal horn and also identified the machinery for the direct binding of mir-500 in the critical site (635–659) in the GAD67 mRNA sequence. These results strongly implied a potential role of mir-500 in the modulation of GAD67 expression in the dorsal horn in rodent models of neuropathic pain. In the present study, we first reported that inhibition of mir-500 by KO or antagomir remarkably attenuated the GAD67 reduction and suppression of GABAergic transmission in the dorsal horn, as well as the behavioral hypersensitivity in the rats after paclitaxel treatment or L5-VRT. It is well known that multiple mechanisms appear to be responsible for neuropathic pain, including ectopic activity in primary afferents, induction of central sensitization, and loss of inhibition (GABAergic neurons and/or glycinergic) in the spinal dorsal horn (Sivilotti and Woolf, 1994; Woolf and Mannion, 1999). In the present study, we only observed the modulation of GABAergic IPSCs in superficial dorsal horn neurons in rodent models of neuropathic pain. Therefore, it is reasonable that mir-500 antagomir treatment completely restored GABA-mediated IPSC frequency, yet the behavioral responses did not return to the baseline level. Together, the present results show for the first time a mir-500-involved mechanism to regulate GABAergic transmission in that dorsal horn and contribute to the induction of neuropathic pain in a rodent model.
Footnotes
This work was supported by the National Natural Science Foundation of China (Grants 81271474, 81300966, and 81171469) and the Fundamental Research Funds for the Central Universities (Grant 15ykjco4b).
The authors declare no competing financial interests.
- Correspondence should be addressed to Wenjun Xin, Zhongshan Medical School, Guangdong Province Key Laboratory of Brain Function and Disease, Sun Yat-Sen University, 74 Zhongshan Rd. 2, Guangzhou, China. xinwj{at}mail.sysu.edu.cn