Single-cell analysis of dorsal root ganglia reveals metalloproteinase signaling in satellite glial cells and pain

Satellite glial cells (SGCs) are among the most abundant non-neuronal cells in dorsal root ganglia (DRGs) and closely envelop sensory neurons that detect painful stimuli. However, little is still known about their homeostatic activities and their contribution to pain. Using single-cell RNA sequencing (scRNA-seq), we were able to obtain a unique transcriptional profile for SGCs. We found enriched expression of the tissue inhibitor metalloproteinase 3 (TIMP3) and other metalloproteinases in SGCs. Small interfering RNA and neutralizing antibody experiments revealed that TIMP3 modulates somatosensory stimuli. TIMP3 expression decreased after paclitaxel treatment, and its rescue by delivery of a recombinant TIMP3 protein reversed and prevented paclitaxel-induced pain. We also established that paclitaxel directly impacts metalloproteinase signaling in cultured SGCs, which may be used to identify potential new treatments for pain. Therefore, our results reveal a metalloproteinase signaling pathway in SGCs for proper processing of somatosensory stimuli and potential discovery of novel pain treatments.


Introduction
Sensory neurons in dorsal root ganglia (DRGs) can detect and differentiate diverse sensations, including pain. Pain is a physiological reaction to hazardous stimuli that threaten homeostasis and survival. However, its abnormal processing -due to damage to the nervous system, cancer, diabetes, infection, autoimmune diseases, or chemotherapy -can result in pathological pain conditions (Basbaum et al., 2009;Scholz and Woolf, 2002). These states are often linked to debilitating chronic pain with symptoms like spontaneous pain, hyperalgesia and allodynia. Though sensory neurons and their related pathways have been viewed as main contributors to chronic pain, current clinical treatments targeting neurons have low analgesic efficacy and frequently have adverse effects Grosser et al., 2017). It's now evident that neurons are not the only ones driving various clinical pain states; non-neuronal cells in the nervous system, such as macrophages, T cells, and glial cells, also play important roles (Ji et al., 2013;Ji et al., 2016;Scholz and Woolf, 2007). Among these cells, satellite glial cells (SGCs) that envelop the soma of peripheral sensory neurons, have been overlooked but are gaining attention as novel players in pain and therapeutic targets (Hanani and Spray, 2020).
Satellite glial cells (SGCs) are abundant in peripheral ganglia, such as DRGs, and have multiple roles: providing metabolic and trophic support to sensory neurons, buffering extracellular ions and neurotransmitters, and enabling neuromodulation (Hanani, 2005;Hanani and Spray, 2020). SGCs have characteristics which make them ideal partners with sensory neurons and immune cells: they respond quickly to changes in neuron activity, they express neurotransmitter receptors and GPCRs, and they can respond to the release of mediators such as ATP, TNF, and IL-1β (Berta et al., 2012;Gu et al., 2010;Kim et al., 2016;Suadicani et al., 2010). These communications maintain neuronal homeostasis (Hanani, 2005;Pannese, 2010) and enable optimized responses to external stimuli. However, excessive or inappropriate communication might lead to excessive inflammatory and pain responses after injury (Grace et al., 2014;Ji et al., 2014). Indeed, injury can not only sensitize neurons, but also activate SGCs in sensory ganglia, resulting in exaggerated pain sensitivity and chronic pain Ji et al., 2013). SGC mediators have been shown to powerfully modulate neuronal functions in acute and chronic pain, and to reduce the effectiveness of opioid analgesia (Berta et al., 2012). Therefore, a better understanding of SGC biology and functions could lead to new therapeutic options for pain treatment.
Here, we have taken advantage of advances in single cell technologies, which allow us to investigate the transcriptomes of nearly all cell populations in the nervous system (Shapiro et al., 2013), to better understand the biology and functions of SGCs in DRGs. These technologies have already provided valuable insights into the transcriptomes of neuronal and non-neuronal cells in DRGs, but most of them have focused on the heterogeneity of peripheral sensory neurons and neuronal regeneration (Avraham et al., 2020;Renthal et al., 2020;Sharma et al., 2020). What has been largely overlooked in these studies is the heterogeneity of SGCs and the identification of new protein signaling pathways related to pain. Using single-cell RNA sequencing (scRNAseq), we identified two major subpopulations of SGCs, one of which was enriched in genes associated with metalloproteinase signaling, including the tissue inhibitor metalloproteinase 3 (TIMP3). Changes in TIMP3 expression and function in vivo influenced modulation of somatosensory stimuli and chemotherapy-induced pain, and in vitro modulation of metalloproteinase signaling in SGCs can be used to screen and identify novel drugs for pain. Our findings reveal a previously unknown molecular heterogeneity of SGCs and metalloproteinase signaling for the maintenance of homeostasis in DRG tissues and the discovery of new pain treatments.

Animal and human studies
All animal experiments conformed to the guidelines recommended by the National Institute of Health, the International Association for the Study of Pain, and the National Centre for the Replacement, Refinement, and Reduction of Animals in Research ARRIVE guidelines. These experiments were approved by the Institutional Animal Care and Use Committee at the University of Cincinnati (UC). Experiments using mice began at 8-10 weeks of age for both male and female CD1 mice (Cat# 022, Charles River Laboratories, Wilmington, MA). RNAscope experiments used Plp1-tdTomato mice, which were generated by crossing the inducible Plip1-Cre/ERT mice (Cat# 005975) with the Rosa26-LoxP-STOP-LoxP-tdTomato mice (Cat# 007909) from the Jackson Laboratory (Bar Harbor, ME). Recombination was induced by oral delivery of a 200 µl tamoxifen/corn oil solution (20 mg/ml) for three consecutive days one week prior to RNAscope. Mice were housed four per cage at 22 ± 0.5 • C under a controlled 14/10 h light/dark cycle with free access to food and water. Studies involving human DRG tissues were performed on deidentified biospecimens and approved by the UC Institutional Review Boards.

DRG tissue dissociation
Adult mice were anesthetized with isoflurane and DRGs were pooled from three male mice from the cervical (C1-C8), thoracic (T1-T13), and lumbar (L1-L6) segments. The DRGs were dissociated according to an adapted protocol (Malin et al., 2007). DRGs were bilaterally removed and placed in ice-cold Hank's balanced salt solution (HBSS, Thermo Fisher Scientific, Waltham, MA) supplemented with 1% pen/strep and 1% HEPES, then finely dissected to remove all nerve roots. The DRGs were incubated in papain in HBSS (60 U/ml,cat# P3125,Millipore Sigma,St. Louis,MO) for 20 min at 37 • C, then in collagenase in HBSS (1 mg/ml, cat# C6885, Millipore Sigma) for another 20 min at 37 • C. They were triturated with 1 ml and then 200 µl pipettes, passed through a 40µm filter, and centrifuged through 15% bovine serum albumin (BSA) solution to remove cell doublets and debris. DRG cells were resuspended in PBS-BSA 0.01% and a small portion counterstained with trypan blue to count viable cells with a hemocytometer. Cells were diluted at 1000 cells/µl and considered as passing quality control if >90% were viable and virtually no cell debris was visible.

Single-cell RNA library preparation, sequencing, and analyses
Single-cell RNA sequencing was conducted by the Gene Expression Core at Cincinnati Children's Hospital using the 10x Chromium Single Cell Kit (v3, 10x Genomics, Pleasanton, CA). Approximately 15,000 single cells were loaded onto the microfluidic chip. Downstream processes included reverse transcription, cDNA synthesis/amplification, and library preparation according to the manufacturer's instructions. The library was sequenced on a NovaSeq 6000 System (Illumina, San Diego, CA). An expression matrix was obtained using 10X's CellRanger pipeline, imported into R, and analyzed with Seurat V3 package (Stuart et al., 2019). Quality control filters included removing cells with a total number of expressed genes outside the range 200-7,000 and more than 10% of reads mapping to mitochondrial genes. A total of 10,172 cells were selected for principal component analysis and visualized by Uniform Manifold Approximation and Projection (UMAP). Clusters were identified at various resolutions using the Louvain modularity optimization as implemented in Seurat V3 package. Optimal clustering resolution was obtained by exploring gene expressions of well-known cellular marker genes (Suppl. Table 1). Differentially expressed genes, defined by having a fold change greater than 1.25 and an adjusted pvalue below 0.05, were identified in Seurat V3 package. Subsequently, the data was analyzed for enrichment of GO terms using Enrichr (Kuleshov et al., 2016) and a background gene list including all genes detected in our sequencing (Suppl. Table 1).

Rnascope (fluorescent in situ hybridization)
Mice were anesthetized with isoflurane and transcardially perfused with 20 ml of PBS followed by 20 ml of 4% paraformaldehyde in PBS at 4 • C. The lumbar L3-L5 DRG were then removed, post-fixed in the same fixative for 2 h, and placed in 30% sucrose solution overnight at 4 • C. Sections of the tissue were cut in a cryostat (12 µm) and mounted onto Superfrost Plus slides (Thermo Fisher Scientific). The slides were dried for 15 min at room temperature and stored at − 80 • C until processing. RNA was detected using the RNAscope Multiplex Fluorescent V2 Assay (Advanced Cell Diagnostics, Newark, CA) with the omission of the initial on-slide fixation step. Probes used included Mm-Gja1-C4 (Cat# 486191-C4), Mm-Entpd2-C2 (Cat# 479511), Mm-Timp3-C2 (Cat# 425321-C2), Mm-Mmp14 (Cat# 432001) and Mm-Adam17 (Cat# 418011). After probe incubation and fluorescence signal amplification, slides were incubated in 1:5000 DAPI (Cat# D1306, Thermo Fisher Scientific) in 0.1 M PB for 1 min before washing, air-drying, and coverslipping with Prolong Gold Antifade mounting medium (Thermo Fisher Scientific). Images were acquired using a Keyence BZ-X800 microscope and analyzed using NIH Image J open source software (Schindelin et al., 2012).

Gene knockdown by small interfering RNA (siRNA)
A total of 3 µg of siRNA targeting Timp3 or a non-targeting control siRNA were diluted in a 10 µl solution containing 2.62 µl of in vivo-jetPEI (Cat# 201-10G, Polyplus, New York, NY), and 5% glucose, and injected intrathecally (i.t.) twice (with an interval of 6 h between injections) on day 0. On day 2 following siRNA administration, we assessed the effect of Timp3 knockdown in mice using mechanical and thermal (heat and cold) sensitivity assays, as well as a motor performance test (see section 2.7 Mouse model and behavioral assays). To initially assess the specificity of Timp3 siRNA over other TIMPs in vitro, we transfected SGC cultures (see section 2.9 Satellite glial cell culture) on day 10 with 25 nM of Timp3 or control siRNA for 48 h. Transfection was achieved by mixing 10 µM of siRNA with 250 μl OptiMEM media (Cat# 31985062, Thermo Fisher Scientific) containing 7.5 μl TransIT-X2 (Cat# MIR 6010, Mirus Bio, Madison, WI). To then assess the in vivo specificity and involvement of various downstream pathways associated with mechanical hypersensitivity in Timp3 knockdown mice, we injected an exogenous recombinant TIMP3 (100 ng/site, i.t.), a general endogenous tissue inhibitor of MMPs (TIMP-1, 4 pmol/site, i.t.), MMP2 and MMP14 inhibitors (10 µg/site, i.t.), TACE/ADAM17 inhibitor (TAPI-2, 1 µg/site, i.t.), and a neutralizing antibody for TNF-α (5 µg/site, i.t.).

Mouse model of neuropathic pain
To produce an animal model of chemotherapy-induced neuropathic pain (CINP), mice were treated with paclitaxel (cat# T7402, Millipore Sigma), as previously described Tonello et al., 2020). In brief, a 6 mg/mL stock of paclitaxel was diluted with Cremophor EL and 95% dehydrated ethanol (1:1 ratio), and administered intraperitoneally at a dosage of 2 mg/kg in saline, every other day for a total of 4 injections. Control animals received an equivalent volume of a vehicle containing proportional amounts of Cremophor EL and 95% dehydrated ethanol, also diluted in saline. Signs of peripheral neuropathy similar to those found in patients have been validated by numerous investigators in this non-tumor-bearing animal model, including the time-dependent development of mechanical and cold allodynia (Flatters et al., 2017;Li et al., 2015). All experiments in naïve and paclitaxeltreated mice were performed by blinded observers, and no more than three assays were performed each testing day.

Behavioral assays
Mechanical sensitivity (von Frey assay) -Paw withdrawal thresholds to mechanical stimuli were evaluated using a calibrated von Frey monofilament series (Stoelting Co., Wood Dale, IL). Mice were first acclimated in individual, clear Plexiglas boxes on an elevated wire mesh platform to access the hind paws' plantar surface. A series of von Frey monofilaments (0.02, 0.07, 0.16, 0.4, 0.6, 1.0, and 1.4 g) was applied perpendicularly to the hindpaw's plantar surface. The test began with the 0.6 g hair. A positive response was a clear paw withdrawal or shaking. Whenever a positive response occurred, the next lower hair was applied, and whenever a negative response occurred, the next higher hair was applied. The testing consisted of 6 up-and-down stimuli and the pattern of response was converted to a 50% von Frey threshold (Chaplan et al., 1994). The anti-allodynic effect of each drug was quantified with respect to the area under the threshold-time curve (trapezoidal method) over the 24 h post-injection testing period. Data are reported as a percentage of the maximum possible anti-allodynia, calculated as a ratio of actual anti-allodynia to a hypothetical situation where the drug brought withdrawal thresholds to their original baseline at all post-injection time points (Sorge et al., 2015).
Heat sensitivity (Hargreaves assay) -The Plantar Test apparatus (IITC Life Science, Woodland Hills, CA) was used to assess paw withdrawal latency in response to heat stimulation, according to the Hargreaves method (Hargreaves et al., 1988). Animals were first acclimated in plastic observation boxes for 20 min. Then, the mid-plantar surface of the mouse hindpaw was exposed to a radiant heat source through a glass floor until paw withdrawal. The intensity of the heat was adjusted to produce a baseline of about 15 s in naïve mice, with a maximum cutoff of 20 s. Each mouse was tested three times, with at least two minutes between recordings. The average of all recordings was used for statistical analysis.
Cold sensitivity (Dry ice assay) -Mice were acclimated to the Plantar Test apparatus (IITC Life Science) for 20 min before cold stimulation with dry ice, according to a previously described protocol (Brenner et al., 2012). A 10 ml syringe, sectioned above the Leur lock and tightly packed with finely crushed dry ice, was pressed firmly on the bottom of the tempered glass directly below the hindpaw. The paw withdrawal latency of the mice was measured with a stopwatch. A standard baseline of 20 s was observed in naïve mice, with a maximum cutoff of 30 s. All mice were tested three times, with at least 5 min between recordings. The average of all recordings was used for statistical analysis.
Motor performance (Rotarod assay) -Motor performance was assessed using a Rotarod Treadmill from IITC Life Science. Mice received two consecutive days of training before testing. The latency to fall was measured with an accelerated rotation speed from 5 to 35 rpm over 3 min. All mice were tested three times with at least 10 min between recordings and the average of all recordings was used for statistical analysis.

Satellite glial cell cultures
Freshly dissociated DRG tissues were used to prepare SGC cultures. Mice were anesthetized with isoflurane and lumbar, thoracic and cervical DRG tissues were aseptically removed. DRG tissues were then incubated in papain (60 U/ml, Cat# P3125, Millipore Sigma) for 30 min at 37 • C followed by collagenase (1 mg/ml, Cat# C6885, Millipore Sigma) for another 30 min at 37 • C. DRG cells were dissociated using a 1 ml and then a 200 μl pipette tip, filtered through a 40 µm and then a 10 µm cell strainer and cultured in DMEM low glucose media (Cat# 12320032, Thermo Fisher Scientific) with 10% fetal bovine serum, 1% pen/strep and 1% fungizone (Cat# 15290018, Thermo Fisher Scientific). Cultures were maintained at 37 • C with 5% carbon dioxide for around 10 days before use. Using this method, we were able to achieve a purity greater than 90% for SGCs. This was assessed by colocalizing glutamine synthetase (GS) with DAPI staining (data not shown).

Cell culture viability
Cell cultures were assessed for viability 24 h after exposure to vehicle or paclitaxel using the PrestoBlue Cell Viability Reagent (Cat# A13261, Thermo Fisher Scientific). Briefly, 50 µl was added to the culture media for 30 min, and 200 µl of media was collected for absorbance quantification. Absorbance was measured at 600 nm using the EnVision plate reader (PerkinElmer, Waltham, MA).

Immunofluorescence
DRG tissues and SGC cultures -DRG tissues were dissected and fixed with 4% paraformaldehyde (PFA) for 2 h. Tissues were then cryoprotected in 30% sucrose and sectioned at 12 µm. Cell cultures were fixed with 4% PFA for 15 min and were processed immediately for immunofluorescence. Tissues and cultures were washed in PBS, blocked with 1% BSA, 0.2% Triton X-100 in PBS (BSA solution), and incubated in the BSA solution overnight at 4 • C with the following primary antibodies: TIMP3 Ab (1:500, Cat# MAB973, R&D Systems), FABP7 Ab (1:1000, Cat# PA524949, Thermo Fisher Scientific) or GFAP Ab (1:1000, Cat# MAB360, Millipore Sigma). Tissues and cultures were then incubated with appropriate secondary antibodies conjugated to Alexa Fluor 546 (1:1000, Thermo Fisher Scientific) for 1 h at room temperature in BSA solution, followed by DAPI (Thermo Fisher Scientific) in 0.1 M PB for 1 min at room temperature. Finally, tissues and cultures were air-dried and coverslipped with Prolong Gold Antifade mounting medium (Thermo Fisher Scientific). Images were acquired using a Keyence BZ-X800 microscope, and images were analyzed using NIH Image J open source software (Schindelin et al., 2012).
Skin tissues -Biopsies from the plantar surface of the hind paws of adult mice were dissected, fixed with 4% PFA for 2 h, and cryoprotected in 30% sucrose overnight. The tissues were then sliced into 40-μm sections and collected in PBS solution for free-floating staining. Sections were blocked for 1 h at room temperature in BSA solution and incubated with primary antibodies against the pan-neuronal marker PGP9.5 (1:500; Cat# Z511601-2; Agilent, Santa Clara, CA) overnight at 4 • C. Following this, the sections were incubated with the secondary antibody Alexa Fluor 488 (1:1000; Thermo Fisher Scientific) for 1 h at room temperature, and DAPI was used for counterstaining. A Keyence BZ-X800 microscope was used to acquire images from multiple sections of skin tissues, which were selected at random. Intradermal neuronal fibers (IENFs) -nerve fibers that crossed the dermal/epidermal junction into the epidermis -were analyzed using NIH Image J open source software (Schindelin et al., 2012). The density of the IENFs was determined as the total number of fibers/length of epidermis (IENFs/cm), as previously described .

Western blotting
Western blot was used to validate the tissue-specific knockdown of TIMP3 in both DRG and spinal cord tissues. The tissues were quickly dissected from adult mice, homogenized in RIPA lysis buffer (Millipore Sigma, Cat# 20-188) containing MS-SAFE protease and phosphate inhibitors (Millipore Sigma, Cat# MSSAFE). Protein concentration was measured with Qubit protein assay kit (Thermo Fisher Scientific, Cat# Q33211) and 20 µg of protein was loaded for each lane. Proteins were separated using NuPAGE™ 4-12% Bis-Tris Protein Gel (Thermo Fisher Scientific, Cat# NP0323BOX) and transferred to PVDF membranes. Blocking was done with 2% BSA followed by overnight incubation of TIMP3 antibody (R&D Systems, Cat# MAB973, 1:500) and GAPDH antibody (Cell Signaling, Cat# 3683S, 1:1000) at 4 • C. After that, the membranes were incubated with appropriate HRP-conjugated secondary antibody (Cell Signaling, 1:2000). Blots were detected using chemiluminescence method (SuperSignal™ West Dura Extended Duration Substrate, Thermo Fisher Scientific, Cat# 34075) and scanned with iBrightTM FL1000 Imaging System (Thermo Fisher Scientific). Bands were quantified using NIH software Image J (Schindelin et al., 2012) and relative values of each protein were normalized by GAPDH protein levels.

Quantitative Real-Time RT-PCR (qPCR)
Dissected DRG tissues and cell cultures were stored in TRIzol reagent (Cat# 15596018, Thermo Fisher Scientific). Total RNA was extracted, quantified, assessed for quality, and converted into cDNA according to the RT-PCR method. Primers were obtained from PrimerBank (Wang et al., 2012) and are listed in Suppl. Table 3. qPCR was performed on a QuantStudio 3 Real-Time PCR System (Thermo Fisher Scientific) using PowerUp SYBR Green Master Mix (Cat# A25741, Thermo Fisher Scientific). All samples were normalized by Gapdh expression. Relative transcriptional expression ratios per condition were calculated using the Pfaffl method (Pfaffl, 2001).

Enzyme linked immunosorbent assay (ELISA)
The ELISA kit for mouse TIMP3 was purchased from MyBioSource (Cat# MBS703363) and used according to the manufacturer's protocol. Cell cultures were harvested 24 h after incubation with paclitaxel (300 nM) alone or in combination with pioglitazone (10 µM), and homogenized in RIPA lysis buffer (Cat# 89900, Thermo Fisher Scientific) containing protease and phosphatase inhibitors. Initial protein concentrations were measured using Qubit TM (Thermo Fisher Scientific).

Statistical analysis
Prism 9.0 (GraphPad, San Diego, CA) was used for statistical analysis. Sample sizes were designed to generate 80% power at two-sided P < 0.05. Specific statistical tests used are specified in each figure legend. Generally, an unpaired, two-tailed Student's t-test was used for analyses between two groups, whereas one-way or two-way ANOVA with a specific multiple-comparisons test was used for multiple groups and/or conditions. Data are presented as mean ± standard error of mean (SEM), with values derived from independent biological replicates. Technical replicates and data transformations are described in the specific methods (e.g. the description of data transformation for anti-allodynic effect is in the behavioral assays section). Adobe Illustrator 25.0 (Adobe, San Jose, CA) was used for illustrations and figure organization.

Single-cell RNA sequencing can identify the major non-neuronal cell types present in DRG tissues
Research in the past has primarily focused on the diversity and function of DRG neurons Sharma et al., 2020;Usoskin et al., 2015), yet the diversity and function of non-neuronal cells surrounding these neurons has received less attention. To characterize these nonneuronal cells, we performed single cell RNA sequencing of cervical, thoracic, and lumbar DRGs pooled from three naïve mice, using the Chromium Single Cell Gene Expression Solution (10x Genomics) and unbiased uniform manifold approximation and projection (UMAP) (Becht et al., 2018) analysis (Fig. 1A). A total of 10,172 single DRG cells were considered (see filtering details in Methods and Suppl. Fig. 1) and organized into 12 distinct clusters (c00-c11, Fig. 1B), with a similar number of unique molecular identifiers (UMIs) and an average of around 2000 detected genes for each cluster (Fig. 1C). Each cluster was classified as a cell subtype based on its expression of unique cell markers ( Fig. 1D and Suppl. Table 1), which were compared to previous studies that have described canonical markers and functions for individual somatosensory neuron subtypes and characterized non-neuronal cells in the nervous system (Suppl. Fig. 2A and 2B). This cell-type classification is consistent with previous RNA-seq of DRG cells and nuclei (Avraham et al., 2020;Renthal et al., 2020;Zeisel et al., 2018). Mapping of well-known marker genes onto the UMAP plots revealed clusters of both neuronal cells (c00 and c10, expression of Tubb3) and non-neuronal cells (c01-c09 and c11, expression of B2m), including macrophages (c07, expression of Aif1) and glial cells (c01-c03, expression of Sox10) (Fig. 1E).

Characterization of molecular and functional heterogeneity of SGCs
The goal of this study was to obtain a comprehensive transcriptome of SGCs, which tightly envelop DRG neuron cell bodies ( Fig. 2A) and regulate their environment and function. SGCs represent the largest subtype of glial cells in the DRG, with two distinct clusters ( Fig. 1B and 1C). Although SGCs are derived from Schwann cells and share some functions with astrocytes (Hanani and Spray, 2020), both clusters are transcriptionally distinct from Schwann cells and astrocytes (Fig. 2B, Suppl. Fig. 3A and Suppl. Table 2). As previously reported (Avraham et al., 2020), transcripts for Glul/GS -a historically used SGC marker -were found in nearly all DRG cells (data not shown). Closer examination of transcripts in each SGC cluster showed that other historically defined markers, such as Slc1a3, Kcnj10, and Gja1, are enriched in SGC I cluster (Fig. 2C), compared to SGC II cluster, which is characterized by exclusive   expression of the enzyme ectonucleoside triphosphate diphosphohydrolase 2 (Entpd2, Fig. 2D). In situ hybridization of Gja1 and Entpd2 confirmed two distinct clusters of SGCs surrounding the same DRG neurons (Fig. 2D). These SGCs, though surrounding the same neurons, presented different biological processes (Fig. 2E). SGC I cluster was enriched in cholesterol and fatty acid process-related genes (e.g., Fabp7, fatty acid binding protein 7, Suppl. Fig. 3B), while SGC II cluster was enriched in genes involved in cellular response to transforming growth factor beta (TGFß) response (e.g., Tgfbr3, TGFß receptor III, Suppl. Fig. 3B) and cell adhesion (e.g., matrix metalloproteinases Mmp2, Mmp14 and disintegrin and metalloproteases Adam17 and Adam23, Suppl. Fig. 3B). These SGC biological processes differed from those associated with Schwann cells and astrocytes (Suppl. Fig. 3C), highlighting unique transcriptional and functional characteristics of the two SGC subtypes.

Confirmation of TIMP3 expression and metalloprotease signaling in SGCs
The top biological process in the SGCs II cluster is associated with extracellular matrix organization (Fig. 2E). Matrix metalloproteinases (MMPs) and their inhibitors, tissue inhibitors of metalloproteinases (TIMPs), are involved in tissue homeostasis and inflammation by regulating extracellular matrix turnover and proteolytic processes (Khokha et al., 2013). Evidence shows expression of MMP2, MMP9 and TIMP1 in peripheral sensory neurons and SGCs, as well as their importance in pain (Kawasaki et al., 2008;Knight et al., 2019;Tonello et al., 2019). Genes enriched in SGCs I have been explored for their role in pain (Gja1 and Kcnj10) and regeneration (Fabp7), while genes in SGCs II are less known, including Mmp2, Mmp14, and Timp3 (Suppl. Fig. 3B). TIMP3 is unique among TIMPs due to its inhibition of a broad range of substrates (Qi et al., 2003), including disintegrin and metalloproteinases (ADAMs) such as Adam17 (Fig. 3A), also enriched in SGCs II (Suppl. Fig. 3B). Confirming this, fluorescent in situ hybridization showed Mmp14, Timp3, and Adam17 transcripts to be highly expressed in SGCs ( Fig. 3B and Suppl. Fig. 4A), and Timp3 co-expressed with Mmp14 and Adam17 in cells surrounding DRG neurons (Fig. 3C). Furthermore, in situ hybridization images from the Allen Mouse Spinal Cord Atlas (https://mousesp inal.brain-map.org/), which include DRG tissues, corroborate our data, showing a restricted Timp3 staining pattern in DRG tissues compared to Gja1 and Tubb3 staining patterns detected in both DRG and spinal cord tissues (Suppl. Fig. 4B). Of note, both mouse and human DRG tissues contained TIMP3 mRNA (Fig. 3D), and immunofluorescence showed characteristic SGC circular staining patterns for TIMP3 protein around human DRG neurons (Fig. 3E). Thus, TIMP3 is expressed in SGCs across species, and may control DRG tissue homeostasis by regulating MMP14, ADAM17, and their downstream signaling.

TIMP3 expression modulates response to somatosensory stimuli in naïve mice
TIMPs are secreted proteins known to control tissue homeostasis and function (Khokha et al., 2013), but the role of TIMP3 in DRG tissue remains unknown. To address this, we used a previously described siRNA approach  and injected a control siRNA or Timp3 siRNA intrathecally two days prior to behavioral, pharmacological, and biochemical testing (Fig. 4A). We confirmed the tissue-specificity of this approach and showed that the expression levels of TIMP3 protein were significantly decreased in DRG, but not in spinal cord tissues, of mice treated with Timp3 siRNA compared to those treated with control siRNA (Fig. 4B and Suppl. Fig. 5. Furthermore, we demonstrated through qPCR that Timp3 siRNA specifically and significantly reduced Timp3 mRNA expression levels in SGC cultures compared to other TIMPs (Suppl. Fig. 5B). Since DRGs contain the cell bodies of peripheral sensory neurons transmitting somatosensory information and SGCs can modulate this information , we tested the innate responses to various somatosensory stimuli in mice treated with Timp3 siRNA and control siRNA. We observed that Timp3 knockdown did not impair motor coordination (Fig. 4C), but caused the development of mechanical allodynia for up to 14 days (Fig. 4D). To a lesser extent, these mice also developed heat hyperalgesia at day 4 and cold allodynia at days 4 and 7 after injection (Fig. 4D).
To validate our siRNA approach, we tested whether mechanical allodynia in Timp3 knockdown mice could be rescued by treatment with a TIMP3 recombinant protein (rTIMP3). We found that rTIMP3 (100 ng/ site, i.t.) significantly reversed the mechanical allodynia in mice previously injected with Timp3 siRNA (Fig. 4E). It is well-reported that TIMP3 mainly functions as an endogenous inhibitor of MMPs and ADAM17 (Qi et al., 2003). Therefore, we tested the effects of general MMP inhibitors (TIMP1 recombinant protein, rTIMP1, 4 pmol/site) and specific inhibitors of MMP2 (MMP2 inh., 10 µg/site) and MMP14 (MMP14 inh., 10 µg/site), which all significantly reversed the mechanical allodynia in mice previously injected with Timp3 siRNA (Fig. 4E). TIMP3 is known to inhibit tumor necrosis factor (TNF)-a converting enzyme (ADAM17), and TNF-a has a central role in modulating DRG and neuronal homeostasis and somatosensory information (Liu et al., 2002;Mohammed et al., 2004;Qi et al., 2003;Segond von Banchet et al., 2009). In line with this, we found that inhibition of ADAM17 (TAPI-2, 1 µg/site) and TNF-a (neutralizing antibody, 10 µg/ site) significantly reversed the mechanical allodynia induced by Timp3 siRNA (Fig. 4E). Importantly, we have demonstrated that rTIMP3, rTIMP1, and TAPI-2 do not significantly alter the mechanical threshold in mice injected with a control siRNA (Suppl. Fig. 6). This finding rules out the possibility that the reversal of Timp3 siRNA-induced allodynia is due to the indiscriminate analgesic efficacy of these inhibitors.
To corroborate the extracellular nature of TIMP3 and our siRNA data, we intrathecally injected mice with a neutralizing antibody (Ab) against TIMP3 (10 µg/site in PBS) or with an IgG control (10 µg/site in PBS). Compared to the control treatment, treatment with TIMP3 Ab induced a significant mechanical allodynia in male and female mice, lasting up to 3 and 6 h, respectively (Fig. 4F). However, this treatment did not induce heat hyperalgesia (Suppl. Fig. 7A) or cold allodynia (Suppl. Fig. 7B), likely due to the transitory inhibition of TIMP3 by the Ab treatment compared to the long-lasting inhibition of the siRNA approach. Together, these data suggest that TIMP3 expression and metalloproteinase signaling (e.g. MMP14 and ADAM17 enriched in SGCs) modulate the responses to mechanical and thermal stimuli.

TIMP3 expression controls the development of chemotherapyinduced neuropathic pain (CINP)
Metalloproteinases are emerging as major players in neuropathic pain (Ji et al., 2009). We and other groups have reported an increase in expression levels of various MMPs and ADAM17 activity in an animal model of chemotherapy-induced neuropathic pain (Cirrincione et al., 2020;Tonello et al., 2019;Wang et al., 2021;Wu et al., 2015). This condition is often refractory to current treatments. We hypothesized that TIMP3, an endogenous inhibitor for various MMPs and ADAM17, may be regulated and, therefore, offers a novel therapeutic approach for CINP. To test this hypothesis, we used a previously reported CINP mouse model (Li et al., 2015;Tonello et al., 2019). In this model, paclitaxel (PAX, 2 mg/kg) is intraperitoneally injected every other day (days 0, 2, 4, and 6 with a final cumulative dose of 8 mg/kg) to mimic the treatment protocols used in clinic. Neuropathic pain with a similar phenotype to that in patients, including development of mechanical and cold allodynia that peak at 14 days after the first PAX treatment (Fig. 5A), has been validated. We found a significant decrease of TIMP3 expression levels in DRG tissues of mice treated with PAX compared to those treated with a vehicle control ( Fig. 5B and C).
Since our data demonstrated that treatment with recombinant TIMP3 (rTIMP3) reversed mechanical allodynia in Timp3 knockdown mice (Fig. 4E), we reasoned that similar treatment would reverse paclitaxel-induced allodynia. Indeed, a single intrathecal delivery of rTIMP3 (100 ng/site) significantly reversed both mechanical and cold allodynia ( Fig. 5D and F). The analgesic effects of rTIMP3 lasted for few hours and were dose-dependent (3 -100 ng/site) for the reverse of mechanical allodynia (Fig. 5E). To investigate whether treatments with rTIMP3 at the time of PAX injections can prevent the development of allodynia, repeated intrathecal administrations of rTIMP3 (100 ng/site, delivered every other day for 4 days at the same time as PAX injections) were conducted. Remarkably, this significantly prevented the development of both mechanical and cold allodynia in male mice for up to 14 days ( Fig. 5H and I). Similarly, paclitaxel-induced mechanical allodynia was also prevented in female mice (Suppl. Fig. 8A, B).
In line with its role in regulating metalloproteinase signaling, these administrations of rTIMP3 significantly reversed the paclitaxel-induced increases of Mmp14 and Adam17 mRNAs in DRG tissues (Suppl. Fig. 9A). Moreover, pro-inflammatory cytokine Tnf expression levels were significantly decreased, whereas the anti-inflammatory cytokine Tgfb1 expression levels were significantly increased in DRG tissues of mice administered rTIMP3 (Suppl. Fig. 9A). However, there were no changes in the expression levels of the antioxidant gene Sod1 in DRG tissues of mice administered rTIMP3 (Suppl. Fig. 9A). Paclitaxel-treatment is known to reduce intra-epidermal nerve fiber (IENF) density in the plantar surface of the paw (Siau et al., 2006), probably due to Data are expressed as mean ± SEM and statistically analyzed by two-tailed t-test (B, C, E) and Two-way ANOVA followed by Sidak's post hoc test (D, F). *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001. mitochondrial damage, oxidative stress and high bioenergetic demand of these peripheral nerve (Bennett et al., 2014). Consistent with our previous report , paclitaxel treatment reduced the IENF density. However, the administration of rTIMP3 did not reverse the loss of IENF density induced by paclitaxel (Suppl. Fig. 9B). Together, these data indicate that rTIMP3 functions as an analgesic in a preclinical model of CINP by regulating metalloproteinase signaling and decreasing neuroinflammation in DRG tissues.

Establishment of a potential SGC discovery platform to identify novel treatments for neuropathic pain
Systemic treatments, including paclitaxel, tend to accumulate in DRG tissues due to their rich vascularization with unique endothelial cells that allow drug penetration (Jimenez-Andrade et al., 2008). To test the hypothesis that paclitaxel can directly target SGCs and alter TIMP3 signaling in DRG tissues, we generated SGC cultures and incubated them with a vehicle control solution or PAX at a concentration of 300 nM (Fig. 6A) -similar to concentrations previously reported in rodent DRG tissues (Cavaletti et al., 2000). We found that PAX significantly increased the expression of the glial fibrillary acid protein (GFAP) in SGC cultures (Fig. 6B), indicating that SGCs are directly activated by paclitaxel exposure without impairing their viability (Fig. 6C). We also sought to test whether TIMP3 signaling is altered in these SGC cultures. In line with our observation in DRG tissues from mice treated with PAX ( Fig. 5B and C), we found that SGC cultures exposed to PAX showed significant decreases in Timp3 mRNA, as well as significant increases in Mmp2, Mmp14, and Adam17 mRNAs (Fig. 6D). Based on these prevent paclitaxel-induced mechanical and cold allodynia (n = 5-6). BL = baseline. Data expressed as mean ± SEM, statistically analyzed by two-tailed t-test (C), two-way ANOVA followed by Sidak's post hoc test (D, F, H, I), and one-way ANOVA followed by Tukey's post hoc test (E). *P < 0.05, **P < 0.01, ***P < 0.001. observations, we hypothesize that drugs that can normalize TIMP3 signaling in SGC cultures exposed to PAX may provide new therapeutic strategies for chemotherapy-induced neuropathic pain.
To test this hypothesis, we used the synthetic prosaptide (Tx14) and the anti-diabetic medication pioglitazone (PGZ). Tx14 is reported to have glioprotective properties and to prevent paclitaxel toxicity (Campana et al., 1998;Meyer et al., 2013). PGZ, in addition to acting on peroxisome proliferator-activated receptor γ, has been demonstrated to rescue the expression of TIMP3 in skeletal muscle of diabetic patients (Tripathy et al., 2013). We found that these drugs did not reverse the effects of PAX on canonical SGC genes, but consistently normalized the transcriptional regulation of Mmp14, and Adam17 ( Fig. 6E and Suppl. Fig. 10A). Although we surprisingly observed transcriptional downregulation of Timp3 in cultures treated with PAX and PGZ, the TIMP3 protein levels increased in these cultures compared to those treated with only PAX (Suppl. Fig. 10B). This suggests a discrepancy between transcript and protein regulation, but it is also consistent with the protein rescue effect previously reported in the literature (Tripathy et al., 2013). More importantly, administration of Tx14 (10 µg/site, i.t.) or PGZ (100 µg/site, i.t.) every other day for 4 days, concomitantly with paclitaxel ( Fig. 6F), significantly prevented the development of paclitaxel-induced mechanical allodynia in mice ( Fig. 6G and H). We also demonstrated that PGZ prevented paclitaxel-induced cold allodynia (Suppl. Fig. 10C). Together, these data suggest a potential cell discovery platform based on the regulation of TIMP3 signaling in SGC cultures exposed to paclitaxel to identify new or repurposed FDA-approved therapeutic drugs.

Discussion
Dorsal root ganglia are complex tissues containing neuronal and nonneuronal cells . However, current analgesic drugs mostly target DRG neurons and are associated with side effects, addiction and limited efficacy . It is well-known that DRG neurons have unique molecular profiles and are tuned to distinct stimuli, including painful stimuli, yet accumulating evidence suggests nonneuronal cells in DRGs also play a role in neuronal homeostasis and processing of these stimuli. Recently, scRNA-seq revealed diverse nonneuronal cells surrounding sensory neurons (Renthal et al., 2020;Zeisel et al., 2018). A similar strategy uncovered unique molecular profiles in SGCs in DRGs and their roles in development and neuronal regeneration (Avraham et al., 2020;Tasdemir-Yilmaz et al., 2021). SGCs are also known for their activation and function in preclinical models of pain (Hanani and Spray, 2020;Ji et al., 2013). To further characterize the molecular profile of SGCs, we conducted scRNA-seq. We identified two distinct SGC populations and a mechanism by which TIMP3 controls homeostasis and pathological pain through SGC metalloproteinase signaling. This regulation of metalloproteinase signaling in cultured SGCs can offer a novel, reliable platform for non-neuronal and nonaddictive pain therapies.
There is a growing body of literature implicating SGCs in neuronal homeostasis and pathological pain (Hanani and Spray, 2020). However, past studies on the role of SGCs in pain have considered them a homogeneous population and focused on a few molecules (Ji et al., 2013). Recently, single-cell sequencing has revealed a more unbiased view of the heterogeneity of glial cells, including SGCs (Avraham et al., 2021;Mapps et al., 2022). In agreement with this, we have identified at least two major SGC subsets in DRGs of naïve mice. We confirmed these distinct subsets through in situ fluorescent hybridization and observed a mosaic distribution around the same DRG neurons, suggesting neuronal-subtype independent distribution of these SGCs. Both subsets (SGCs I and II) may affect neuronal functions and participate in homeostasis and various pathological pain conditions. It is well-documented that SGCs can influence DRG neurons by releasing cytokines and altering expression of proteins such as GLAST1 (encoded by the gene Slc1a3), KIR4.1 (encoded by the gene Kcnj10), and CX43 (encoded by the gene Gja1) Retamal et al., 2017;Silva et al., 2017;Takeda et al., 2011). Interestingly, we found an enrichment of these genes in SGCs I, suggesting a major role of this SGC subset in pain. SGCs I also had an enrichment of genes related to fatty acid synthesis and metabolism, which have been linked to neuronal regeneration (Avraham et al., 2020).
Which genes are enriched in SGCs II subset? What are their functions? Functional analysis of SGCs II revealed a potential role in glial cell differentiation, which may support previous reports suggesting SGCs as a population of developmentally arrested Schwann cells that can proliferate after peripheral nerve injury (Donegan et al., 2013;George et al., 2018). However, our data indicated that both SGCs I and SGCs II are transcriptionally distinct from Schwann cells. Whether SGCs proliferate after injury remains a matter of debate (Jager et al., 2020). We found enrichment of genes associated with extracellular matrix organization, including various matrix metalloproteinases and disintegrin and metalloproteinase proteins. These transcripts include Adam10, Adam17, Mmp2, Mmp14, and the tissue inhibitor of metalloproteinase Timp3. TIMP3 is unique for its inhibition of several membrane-bound molecules with sheddase functions, such as MMP14, MMP3, and ADAM17 (which is responsible for the maturation of TNF) (Mohammed et al., 2004;Qi et al., 2003). Due to its plethoric signaling and functions, TIMP3 has become increasingly recognized for its role in homeostasis, neurodegeneration, and inflammation (Mohammed et al., 2004;Qi et al., 2003;Rosenberg, 2009).
We confirmed strong expression of TIMP3 in SGCs surrounding DRG neurons and its colocalization with MMP14 and ADAM17, suggesting a role of TIMP3 in somatosensation by controlling these metalloproteases. To determine whether TIMP3 is involved in somatosensation, we used local knockdown of Timp3 in DRG tissue of naïve mice using siRNA. These mice developed somatosensory hypersensitivities, including robust and long-lasting mechanical allodynia that persisted for up to 14 days. We previously reported that siRNA knockdown lasts about 5-7 days (Berta et al., 2012;Berta et al., 2014;Lee et al., 2018), therefore the mechanical allodynia lasted beyond the knockdown period, suggesting major homeostatic changes in the DRG tissues. This may be due to the involvement of TIMP3 in regulating multiple downstream target proteins, including metalloproteinases and inflammatory mediators. Indeed, pharmacological blockade of various metalloproteases, ADAM17, and extracellular TNF significantly reversed the mechanical allodynia caused by knockdown of TIMP3. TNF is a key mediator of inflammation and pain . TNF levels are increased in DRGs in animal models of inflammatory and neuropathic pain (Leung and Cahill, 2010). Increased expression and release of TNF are observed in SGCs of mouse models of acute herpetic neuralgia and paclitaxelinduced neuropathic pain (Silva et al., 2017;Wu et al., 2015). Exogenous delivery of TNF in naïve mice elicited mechanical and thermal hypersensitivities, and direct application of TNF to cultured DRG neurons increased their excitability (Jin and Gereau, 2006). Therefore, we suggest TIMP3 in SGCs as a key physiological regulator of neuronal and pain hypersensitivities by controlling metalloproteinase signaling and inflammatory mediators such as TNF.
To assess the role of TIMP3 in pathological pain, we used a paclitaxel-induced mouse model of CINP. This model was used for several reasons. Firstly, paclitaxel increases GFAP expression in SGCs, suggesting their involvement in this type of pain (Warwick and Hanani, 2013). Secondly, it has been reported that paclitaxel increases expression levels of proteins related to metalloprotease signaling, such as MMP2, MMP9, ADAM17 and TNF Wu et al., 2015). Lastly, there is currently no FDA-approved therapy to treat chronic intractable cancer pain (CINP) in cancer patients, which may contribute to the reported misuse of opioids by these patients (Kwon et al., 2013;Reyes-Gibby et al., 2016). Here, we suggest TIMP3 as a novel therapeutic target for CINP, showing that rescuing decreased TIMP3 expression levels in DRGs with an exogenous recombinant TIMP3 protein can reverse and prevent paclitaxel-induced mechanical and cold allodynia in male mice.
Previous studies have revealed sex differences in immune and glial cell signaling that underlies mechanical hypersensitivity (Mogil, 2020). Our results showed that TIMP3 inhibition with a neutralizing antibody induced mechanical allodynia in both male and female naïve mice, implying that SGC signaling is independent of sex. In addition, our lab's previous data has shown that astrocytic signaling, which is similar to that observed in SGCs and responsible for mechanical allodynia, is also independent of sex. For example, TNF induced similar levels of mechanical hypersensitivity in both naïve males and females (Chen et al., 2018). However, sex-dependent glial signaling has been mostly observed in mice with pathological pain (Sorge et al., 2015;Tansley et al., 2022;Taves et al., 2016), and may not be evident in naïve mice. To address this, we tested and demonstrated that recombinant TIMP3 treatment can prevent paclitaxel-induced mechanical allodynia not only in male, but also in female mice. Overall, SGC signaling seems independent of sex and TIMP3 treatment may prevent neuropathic pain in both male and female mice.
Mechanistically, we have shown that this treatment normalizes TIMP3 signaling by reversing the increased expression of Mmp14 and Adam17 mRNA induced by paclitaxel. Additionally, it potentially reduces neuroinflammation in DRGs by decreasing the transcriptional expression of the pro-inflammatory cytokine Tnf, while increasing the transcriptional expression of the anti-inflammatory mediator Tgfb. Further studies are ongoing to elucidate how TIMP3 regulates neuroinflammation in DRG tissues. In addition to its well-known role in the maturation of pro-inflammatory cytokines (Mohammed et al., 2004), it may also regulate the infiltration of immune cells in DRG through its action on the vascular endothelial growth factor receptor-2 (Qi et al., 2003). Surprisingly, TIMP3 treatment did not rescue the loss of IENF density. Although the loss of IENF density is a common characteristic of chemotherapy-induced neuropathies (Bennett et al., 2014), its contribution to pain remains unclear. Interestingly, pre-clinical and clinical studies have found that paclitaxel-induced pain often precedes deficits in sensory nerve function (Boyette-Davis et al., 2015;Ma et al., 2018), and pain from a high dose of paclitaxel can occur independently from functional alterations of sensory nerves (Flatters and Bennett, 2006). Targeting neuroinflammation in DRGs to treat CIPN is supported by multiple preclinical studies, which have demonstrated that depletion of macrophages, blockade of TNF, and promotion of anti-inflammatory cytokines IL-10 and TGFβ all reduce paclitaxel-induced pain (Krukowski et al., 2016;Tonello et al., 2020;Zhang et al., 2016). Nonsteroidal anti-inflammatory drugs are used to treat cancer patients, but their analgesic effects against CIPN require further clinical studies (Desforges et al., 2022).
Chemotherapeutic treatments tend to accumulate in DRG tissues due to their rich and fenestrated vascularization (Jimenez-Andrade et al., 2008). Few studies, however, have investigated the direct interaction of paclitaxel with glia cells (Liu et al., 2019;Wu et al., 2015). Our culture studies show that paclitaxel increases the reactivity of SGCs, which is consistent with the transcriptional changes in TIMP3 signaling observed in the DRGs of mice treated with paclitaxel. This culture may be used to screen different drugs foe CINP. As a proof-of-concept, we tested the ability of the glioprotective peptide Tx14 (Campana et al., 1998;Meyer et al., 2013) and antidiabetic medication pioglitazone to reverse the transcriptional changes in TIMP3 signaling in cultured SGCs treated with paclitaxel. We found that transcriptional changes in Mmp14 and Adam17 transcripts were increased in paclitaxel-treated culture, but normalized by both drugs Tx14 and pioglitazone. More importantly, both drugs were able to prevent CINP in mice. We hypothesize that pioglitazone prevents CINP partially through the increase of TIMP3 expression in DRGs, as it can rescue the expression of TIMP3 in skeletal muscle of diabetic patients (Tripathy et al., 2013). However, the mechanisms underlying the regulation of TIMP3 signaling in SGC by pioglitazone and its analgesic effects warrant further investigation. It is also important to consider that pioglitazone can also alleviate neuropathic pain through spinal mechanisms involving peroxisome proliferator-activated receptor γ (Griggs et al., 2016;Griggs et al., 2015). However, a recent study has proposed that activation of PPARγ by pioglitazone may mitigate lupus nephritis symptoms in mice by suppressing miR-21-5p and promoting the increase of TIMP3 (Liu and Zhang, 2021). Interestingly, we observed different regulation of TIMP3 mRNA and protein in SGC cultures treated with pioglitazone. The correspondence between mRNA and protein expression levels is notoriously poor, in part due to different time courses and post-transcriptional regulation (Vogel and Marcotte, 2012). This may also explain the increased level of GFAP protein, but decreased mRNA levels, in SGC cultures treated with paclitaxel. Previously, we showed that voltagegated sodium channel β2 subunit protein levels were up-regulated in DRG tissues by nerve injury despite no mRNA regulation (Pertin et al., 2005).
This study fully characterized the transcriptional profile of SGCs and identified new therapeutic targets for pain treatment. Despite these strengths, limitations exist. First, the number of genes for most dorsal root ganglion (DRG) neurons was lower than expected in the single-cell RNAseq, preventing us from analyzing their heterogeneity or interactions with other non-neuronal cells. This is likely due to the technical limitation of 10X Genomic microfluidic channels, which are optimized for cells smaller than DRG neurons. Glia and immune cells, which are smaller than DRG neurons, have a number of genes consistent with other single-cell RNAseq studies Zeisel et al., 2018). Single-nucleus RNAseq can be used for better transcriptional analysis of various DRG neuronal populations (Renthal et al., 2020). Second, our study focused on metalloproteinase signaling because of its involvement in tissue homeostasis, inflammation and pain (Ji et al., 2009). We propose targeting TIMP3 to alleviate pain since it is expressed in both mouse and human SGCs and is key in regulating various MMPs and pro-inflammatory and algesic mediators, such as TNF. For this purpose, we delivered the recombinant TIMP3 protein to the DRG tissues via intrathecal injection, a commonly used delivery route in preclinical studies . However, we did not confirm whether the TIMP3 protein reached the DRG tissues, and intrathecal injections may not be practical in a clinical setting. Further studies are needed to address these limitations. It is also important to note that the therapeutic potential of targeting metalloproteinase signaling has been hampered by lack of specificity and functional compensation among various MMPs (Khokha et al., 2013;Vandenbroucke and Libert, 2014). Future studies should investigate other potential therapeutic targets such as ion channels and G-protein-coupled receptors. Finally, cultured SGCs may not fully recapitulate their heterogeneity and state in vivo (Avraham et al., 2021;Jager et al., 2022). Nevertheless, we showed that paclitaxel induces similar transcriptional changes in SGC cultures and mouse DRGs, and provided a proof-of-concept for their use to identify new analgesics for CINP.
In summary, pain studies have mainly focused on neurons, yet only a few new analgesics have been introduced in the clinic over the past decades. Satellite glial cells are now seen as important players in neuronal homeostasis and disease. Here, we present the transcriptional profile of SGCs and explain a previously unknown mechanism of TIMP3 signaling, which modulates nociception and neuropathic pain. The therapeutic potential of TIMP3 signaling is supported by its expression in human SGCs, role in pain modulation in mice, and development of a potential cellular discovery platform to find new analgesics or repurpose existing FDA-approved drugs.

Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Data availability
Data will be made available on request.