Animals
Male C57BL/6 mice, 8 weeks old, were used in this study unless otherwise noted. Mice were housed under 12-hour light/dark cycles in a pathogen-free room with clean bedding and free access to food and water. Cage and bedding changes were performed each week. All animal studies were performed following protocols approved by the Animal Experimentation Ethics Committee, the Chinese University of Hong Kong or Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Science.
Spared Nerve Injury (Sni) Surgery
Mice were anesthetized by gas anesthesia (isoflurane 1–2%). For induction of SNI neuropathy, the left sciatic nerve was exposed at the level of the trifurcation into the sural, tibial, and common peroneal nerves. The tibial and common peroneal nerves were tightly ligated and severed, leaving the sural nerve intact. For the sham group, the left sciatic nerve was exposed as in the SNI procedure but was not further manipulated.
Behavioral Assessment
For the von Frey test, mice were habituated on an elevated platform with a mesh floor for 30 min. The plantar surface of the hind paw was stimulated with a series of calibrated von Frey filaments (NC12775-99, North coast medical, Morgan Hill, CA, USA). The 50% paw withdrawal thresholds were calculated using the Up-Down method.
For the pinprick test, the plantar hind paw was stimulated by stroking using a 25G safety pin. (Ratings: 0 = no response, 1 = paw withdrawal, 2 = flicking of the paw, and 3 = licking of the paw).
For the chemogenetics experiment, behavioral tests were performed at 30–120 minutes after intraperitoneal (i.p.) administration of clozapine-n-oxide (CNO) (3 mg/kg, BrainVTA, Wuhan, Hubei, China) or saline to virus-infected mice.
Drug Administration
R-HNK (SML1873, Sigma-Aldrich, St Louis, MO, USA) or saline was administered via i.p. injection (10 mg/kg) twice daily for 3 consecutive days. For continuous Bdnf delivery (14 days), recombinant Bdnf protein (0.5 µg/90 µL, PeproTech, 450-02, USA) or saline (0.9%, 90 µL, vehicle) was loaded into ALZET osmotic pumps (1003D, RWD, Shenzhen, Guangdong, China) implanted in PrL.
Transcriptome Analysis
Mice were divided into three groups and treated accordingly (n = 3 mice/group), namely Sham/Saline (Sham), SNI/Saline (SNI), and SNI/R-HNK (R-HNK). On the 14th day after surgery, PrL tissues were collected and snap-frozen in liquid nitrogen. Total RNA was isolated from the PrL using TRIzol Reagent (15596026, Thermo Fisher Scientific, Pittsburgh, PA, USA) according to the manufacturer’s instructions. RNA was purified with RNeasy Micro Kit 50 (74004, Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Library construction and sequencing services were provided by Novogene using the Illumina Novaseq 6000 system (Novogene, Tianjin, China) following a standard protocol. The sequencing data were then pre-processed to remove low-quality reads and adaptors, hisat2 (version 2.0.5) to align to mouse genome (version GRCm38), and feature counts (version 2.0.3) to quantify gene counts. Differentially expressed genes (DEGs) were identified using the DESeq2 R-package with default settings. Pathways and gene ontology (GO) annotations were conducted with Enrichr (http://amp.pharm.mssm.edu/Enrichr/) using DEGs between sham and SNI groups. The clustering of similar GO terms under the Leiden algorithm was done using the enrichment analysis visualizer provided by Enrichr Appyters. (https://github.com/MaayanLab/appyter-catalog/tree/main/appyters/Enrichment_Analysis_Visualizer). For gene set enrichment analysis (GSEA), the entire gene list was ranked according to the values of \(-log10\left(pvalues\right) sign\left(log2Foldchange\right)\) of each gene. fGSEA R-package was used to analyze GO biological processes. The log2Foldchange and p-values were generated with DESeq2 R-package. The enrichment plot was drawn with fGSEA. The volcano plots, scatter plots, bar charts, bubble plots, and ridge plots were drawn with the ggplot2 R-package.
Stereotaxic Injections
Stereotaxic injection of recombinant adeno-associated virus (rAAV) was performed according to a standard stereotaxic surgery procedure using an automated nanoliter syringe (Nanoject II, Drummond Scientific, Broomall, PA) with the following parameters: a total volume of 300 nL or 250 nL at 40 nL/minute for injection and hold for 5 minutes before withdrawal or the micropipette. The rAAVs used in this study were listed in the Supplementary file. S1. For PrL transduction, viral particles were injected at anteroposterior (AP): +1.8–1.9 mm; medial-lateral (ML): -0.5 mm; dorsal-ventral (DV): -2.1 mm. rAAV was infused into the ventrolateral PAG (vlPAG) at AP: +4.6 mm, ML: -0.8 mm, and DV: -3 mm.
Fiber Photometry
Fiber photometry was performed using the Multi-Channel Fiber Photometry Device (Inper, Hangzhou, Zhejiang, China). We calibrated the power of the LED lamp (λ = 465 nm) that excites GCaMP to 30 µW and the power of the LED lamp (λ = 405 nm) that records the baseline to 10 µW. We recorded CaMKII neurons in PrL of mice injected with rAAV9-CaMKIIa-GCaMP6m-WPRE-hGH-pA. To record vlPAG-projecting neurons in PrL, rAAV9retro-CaMKIIa-Cre-EGFP-WPRE-pA and rAAV9-EF1a-DIO-GCaMP6f-WPRE-hGH-pA were injected into PrL and vlPAG respectively. Mice were allowed to recover for 7 days, followed by pain model development and drug treatment (2 ⋅ 3 injections). Mice were habituated to the behavioral assessment environment three days before the recording day (postoperative day 14). On the day of the test, mice were connected to the photometric system and acclimated in a chamber (as in the mechanical threshold test) on a wire mesh for 30 minutes. The GCaMP signals were recorded in an 8-minute window divided into three slots: resting state: 3 minutes; stimuli sections (with 1 g of von Frey filament every 40 s on the left hind paw): 4 minutes, and then another resting state: 1 minute.
Whole-cell Patch Clamp Recordings
For the patch clamp study, 4-week-old mice were used for pain model development and drug treatment. Brain slices were collected for recording at 7 weeks of age. Brains were quickly harvested and transferred to ice-cold dissection solution equilibrated with 95% O2 and 5% CO2. Coronal sections (300 µm) were cut with Leica VT1200S (Leica Biosystems, Buffalo Grove, IL, USA) and incubated at 30℃ for 30 minutes in artificial cerebrospinal fluid (aCSF, 126 mM NaCl, 2.5 mM KCl, 125 mM NaH2PO4, 26 mM NaHCO3, 10 mM D-Glucose, 0.5 mM CaCl2, and 10 mM MgSO4) equilibrated with 95% O2 and 5% CO2. Layer 5 pyramidal neurons in PrL were identified under 40 X magnification. The recording electrode was filled with intracellular solution (124 mM NaCl, 2.5 mM KCl, 1.2 mM NaH2PO4, 24 mM NaHCO3, 5 mM HEPES, 12.5 mM D-Glucose, 2 mM CaCl2, 1.5 mM MgSO4, and 1.5 mM QX-314). The patch-pipettes have a resistance of 6–9 mΩ. We placed the tip of the glass electrode at 1/3 − 2/3 of the target neurons. For the whole-cell recording, the voltage was held at -70 mV to record the sEPSC of pyramidal neurons, and at the voltage was held at + 10 mV for the sIPSC. The total recording time for each cell was 3 minutes and data from the median 1 minute was used for the analysis.
The Ratio Of Ampar/nmdar Eepscs
Acute brain slices were prepared as above and placed in aCSF containing 50 µM bicuculine. To measure AMPAR/NMDAR eEPSC ratio, the recording electrode containing the intracellular solution was attached to the target neuron in PrL, while the stimulation electrode was placed 200–250 µm away from the target neuron. The eEPSCs were then evoked with increasing currents ranging from 0-160 µA with an increment of 20 µA at each step. The AMPAR- and NMDAR-eEPSCs were recorded with a holding potential of -70 mV and + 40 mV, respectively.
Dendritic Spines Analysis
In each group, 2–4 segments per CaMKII neuron, including 10–16 neurons, from 5 animals were used to reconstruct z-stacks into 3D models for analysis. To quantify dendritic spine density and morphology, Z-stacked images were converted to maximum intensity projections and analyzed using Imaris 7 software (Bitplane, Concord, MA, USA). Images of dendritic fragments were collected from secondary to tertiary dendrites. The Filament Tracer module with default settings and the auto-path mode were used for tracing the dendritic fragments. If necessary, we manually corrected the dendritic fragments and spines detection. We calculated dendritic spine density by counting the number of spines per 10 µm of dendritic length. Co-localization of spines with GluA1 puncta was 3D-reconstructed and analyzed using the “surface” and “spot” functions of Imaris. Dendritic spines were classified into mushroom spines, thin spines, stubby spines, and filopodia spines using Imaris X Tension, and the density of each spine type was calculated. Masked surface channels were created for the dendritic shaft, filament, and each spine class. Colocalized puncta in the dendritic spines were subtracted from the total co-localized puncta in each filament to determine the co-localized voxels within the spines. The co-localized volume of the individual spine was determined by co-localizing the masked spine channel with the masked filament channel of interest.
Statistical Analysis
Data were checked for normality using the Shapiro-Wilk test. Tests of significance were conducted using one-way ANOVA, followed by Bonferroni's post hoc tests. All statistics were performed using GraphPad Prism (version 9, GraphPad Software, San Diego, CA, USA) with statistical significance set at alpha = 0.05. Data were presented as mean ± standard error of the mean (mean ± SEM).