Single-nuclei RNA Sequencing Reveals Distinct Transcriptomic Signatures of Rat Dorsal Root Ganglia in a Chronic Discogenic Low Back Pain Model

Single-nuclei RNA Sequencing Reveals Distinct Transcriptomic Signatures of Rat Dorsal Root Ganglia in a Chronic Discogenic Low Back Pain Model | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results Single-nuclei RNA Sequencing Reveals Distinct Transcriptomic Signatures of Rat Dorsal Root Ganglia in a Chronic Discogenic Low Back Pain Model View ORCID Profile Sydney M. Caparaso , View ORCID Profile Ishwarya Sankaranarayanan , David J. Lillyman , View ORCID Profile Theodore J. Price , View ORCID Profile Rebecca A. Wachs doi: https://doi.org/10.1101/2025.02.19.639130 Sydney M. Caparaso a Department of Biological Systems Engineering, University of Nebraska-Lincoln , Lincoln, USA ; Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Sydney M. Caparaso Ishwarya Sankaranarayanan b Department of Neuroscience, Center for Advanced Pain Studies, University of Texas at Dallas , Dallas, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Ishwarya Sankaranarayanan David J. Lillyman a Department of Biological Systems Engineering, University of Nebraska-Lincoln , Lincoln, USA ; Find this author on Google Scholar Find this author on PubMed Search for this author on this site Theodore J. Price b Department of Neuroscience, Center for Advanced Pain Studies, University of Texas at Dallas , Dallas, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Theodore J. Price Rebecca A. Wachs a Department of Biological Systems Engineering, University of Nebraska-Lincoln , Lincoln, USA ; Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Rebecca A. Wachs For correspondence: rebecca.wachs{at}asu.edu Abstract Full Text Info/History Metrics Preview PDF Abstract Chronic low back pain (LBP), often correlated with intervertebral disc degeneration, is a leading source of disability worldwide yet remains poorly understood. Current treatments often fail to provide sustained relief, highlighting the need to better understand the mechanisms driving discogenic LBP. During disc degeneration, the extracellular matrix degrades, allowing nociceptive nerve fibers to innervate previously aneural disc regions. Persistent mechanical and inflammatory stimulation of nociceptors can induce plastic changes within dorsal root ganglia (DRG) neurons, characterized by altered gene expression, enhanced excitability, and lowered activation thresholds. Although these transcriptional changes have been described in other pain states, including osteoarthritis, they remain underexplored in discogenic LBP. To address this gap, this study represents the first application of comprehensive single-nuclei RNA sequencing of DRG neurons in a rat model of chronic discogenic LBP. Eighteen distinct DRG subpopulations were identified and mapped to existing mouse and cross-species atlases revealing strong similarities in neuronal populations with the mouse. Differential expression analysis revealed increased expression of pain-associated genes, including Scn9a and Piezo2 , and neuroinflammatory mediators such as Fstl1 and Ngfr , in LBP animals. Axial hypersensitivity, measured using grip strength, significantly correlated with increased expression of Scn9a, Fstl1, and Ngfr, which suggests their role in maintaining axial hypersensitivity in this model. These findings establish a relationship between DRG transcriptomic changes and axial hypersensitivity in a discogenic LBP model, identifying potential molecular targets for non-opioid treatments and advancing understanding of discogenic LBP mechanisms. 1. Introduction Chronic low back pain (LBP) affects over 500 million people worldwide, resulting in significant disability and reduced quality of life [ 4 , 13 , 27 , 31 ]. Disc degeneration is highly associated with chronic LBP [ 11 , 15 , 42 , 44 , 62 , 70 , 73 ]. In up to 42% of chronic LBP patients, the pain is attributed to the disc, termed discogenic LBP [ 28 , 45 , 47 ]. However, existing treatments— including nonsteroidal anti-inflammatory drugs (NSAIDs) [ 60 , 71 ], opioids [ 32 ], and surgery [ 8 ]— often fail to provide long-term relief [ 18 , 43 , 59 ], highlighting the need to develop a comprehensive understanding of the mechanisms driving discogenic LBP. Patients with discogenic LBP exhibit biological and anatomical changes within the disc. During disc degeneration, the extracellular matrix (ECM) breaks down, disrupting its ability to maintain hydration and native neuroinhibitory properties [ 33 , 36 ]. This breakdown allows nociceptive nerve fibers to innervate deep into the disc through axonal sprouting. Innervation has been observed in previously aneural regions including the nucleus pulposus (NP) core and the inner annulus fibrosus (AF) [ 7 , 20 , 24 ]. In addition, this mechanical breakdown of the disc is often observed with a catabolic microenvironment [ 10 , 34 , 74 ]. Quantitative sensory testing (QST) studies in LBP patients have shown decreased pain thresholds and heightened responses to mechanical [ 12 , 53 ] and cold stimuli [ 12 , 52 ] suggesting a process of nociceptor plasticity, or transcriptional changes to peripheral and central neurons with an implication on the dorsal root ganglia (DRG) as the primary site for sensory input. Importantly, studies of genetic trait loci in human DRGs have observed links between chronic pain conditions, such as temporomandibular disorder, and pathways involving nociceptor sensitization and immune responses [ 58 ]. However, there is a gap in understanding how transcriptional changes in nociceptors contribute to discogenic LBP. Animal models have replicated several features of chronic discogenic LBP, including ECM breakdown, inflammation, and nociceptor innervation into the disc [ 2 , 40 , 49 , 55 ]. These animal models also recapitulate some subsets of pain-like behaviors like those exhibited in human patients [ 22 , 25 , 35 , 37 , 39 , 40 , 41 , 50 , 56 ]. Evidence from preclinical models of other pain states, such as knee osteoarthritis (OA), suggest that DRG neurons likely undergo plastic changes observed through the upregulation of ion channels and inflammatory markers [ 54 , 76 ]. Importantly, pain-related ion channel targeting has alleviated pain-like behavior in OA models; specifically, Nav1.7 [ 30 ], Trpv4 [ 68 ], and Piezo2 [ 54 ] have shown efficacy in mitigating OA progression and reducing pain-like response. Therefore, characterization of nociceptor plasticity in the context of chronic discogenic LBP could provide valuable insights into the molecular drivers of this disease. Therefore, the goal of this work was to characterize DRGs in our discogenic rat LP model to understand nociceptive plastic changes in response to injury and correlation with pain-like behavior. To our knowledge, this study represents the first application of comprehensive single-nuclei RNA sequencing (snRNAseq) of rat DRG neurons. Herein, a high-resolution transcriptomic atlas of DRG subpopulations was generated to identify distinct nociceptor subtypes and their gene expression profiles in response to disc injury. In addition, our findings advance our understanding of DRG-specific changes associated with painful disc degeneration and highlight potential gene targets that correlate with pain-like behaviors. 2. Materials and methods 2.1. Study Design and Animals A comprehensive dataset of dorsal root ganglia (DRGs) from 12 healthy adult female and male Sprague Dawley rats (Envigo) was deeply sequenced to establish a robust anchor dataset. This dataset was critical for accurately annotating rat neuronal subpopulations for subsequent pain studies. Herein, this study is referred to as the “anchor study,” as shown in ( Fig. 1 ). For this anchor study, 6 rats (n=3 male and n=3 female) were harvested from both the University of Nebraska-Lincoln and the University of Texas at Dallas by two unique experimenters for a total of 12 rats. All procedures were conducted under approved Institutional Animal Care and Use Committee (IACUC) protocols at respective universities. This anchor dataset initially identified sixteen distinct clusters, providing unique DRG subpopulation annotations and establishing a foundational rat DRG atlas. After establishing the anchor dataset for snRNAseq of rat DRGs, a “chronic LBP study” examined transcriptional changes in DRG neurons in a model of chronic discogenic LBP as shown in ( Fig. 1 ). For this study, 12 adult female Sprague Dawley rats (Envigo) from a larger cohort were randomly selected in (n=6 injured) and (n=6 sham) groups for sequencing. These samples were mapped onto the anchor dataset, enabling confident label transfer and precise classification of DRG subpopulations. Download figure Open in new tab Figure 1. Schematic overview of experimental design for defining transcriptomic clusters in rat DRGs from chronic discogenic LBP model Top (Cohort 1) : SnRNA seq of T13-L1 DRG was performed from uninjured control rats, and the analysis was conducted to create a high-resolution anchor dataset. T13-L1 bilateral DRGs from healthy male and female Sprague Dawley rats were harvested from UNL and UTDallas experimenters (n=12; 3 male and female per site at UNL and UTDallas). This dataset serves as a reference map for DRG subpopulations in healthy rats. Top (Cohort 2) : Adult female Sprague Dawley rats underwent surgical intervention to induce mechanical disc degeneration at the L5-L6 intervertebral disc level, followed by a 15-week period to allow for the development of painful disc degeneration. T13-L1 bilateral DRGs were harvested from sham and injured rats (n=6 per group). Labels defined from anchor study were transferred to the chronic LBP study. Bottom : Schematic of nuclei isolation for snRNAseq and clustering. Nuclei were isolated from DRGs, mechanically dissociated, and filtered to ensure single-nuclei suspension. Nuclei suspensions were then barcoded and sequenced using the Parse Bioscience. The data were analyzed for clustering visualization, revealing distinct clusters in the DRG. subpopulations based on gene expression profiles. Clusters were annotated to identify DRG subpopulations according to prior published human and mouse DRG datasets [ 6 , 61 ]. Dorsal root ganglia (DRG); low back pain (LBP); single-nuclei RNA sequencing (snRNAseq); University of Nebraska-Lincoln (UNL); University of Texas at Dallas (UTDallas). 2.2. Anchor study 2.2.1. Humane euthanasia and DRG tissue harvest Animals were humanely euthanized using CO 2 inhalation as the primary method, followed by bilateral pneumothorax as a secondary method. DRGs from both sexes were carefully harvested based on methods from our lab [ 23 ]. In short, the spinous and transverse processes were removed, the spinal cord was carefully cut out, and the DRGs were extracted. During harvest, the connective tissue, spinal, and peripheral nerve processes were trimmed from the DRG explants. Bilateral T13-L1 DRGs were combined per animal for snRNAseq and flash frozen in liquid nitrogen before further processing. These DRG levels were selected as L1 DRGs are a primary contributor to innervation in the L5-L6 disc in rats [ 49 , 55 , 65 ]. 2.2.1.1. DRG nuclei isolation from adult Sprague Dawley rats Prior to handling DRG samples, all surfaces and tools were cleaned with 70% ethanol and an RNase decontamination solution (RNAzap, Thermofisher). Nuclei were isolated based on methods from the manufacturer’s protocol (Parse Evercode V2 kit) and literature [ 6 ]. In short, T13-L1 bilateral DRGs per animal were chopped into <1mm small pieces using Bonn scissors (14084-08, FST) in a 1.5 mL DNA LoBind tube (Eppendorf, 0030108418) with 500 µL of homogenization buffer followed by mechanically dissociating using the Dounce homogenizer with 15 pestle strokes ( Fig. 1 ). The buffer was freshly prepared by combining 0.25 M sucrose, 150 mM KCl, 5 mM MgCl₂, 1 M Tris buffer, pH 8.0, 0.1 mM DTT (ThermoFisher, D1532), cOmplete protease inhibitor (Sigma, 11836153001), 0.1% Triton X-100 (Sigma, T8787), and 0.2 U/µL RNAsin (Promega, N2515). Then, the tissue was filtered through a 70 µm cell strainer (Corning, 352360) into a 50 mL conical tube. The nuclei suspension was centrifuged at 500 x g for 5 minutes at 4°C, and the supernatant was discarded. The nuclei pellet was washed once in wash buffer containing 1% RNase-free BSA, 0.2 U/µL RNase inhibitor, 1X DPBS without calcium and magnesium (ThermoFisher, 14190144) and centrifuged at 500 x g for 5 minutes. The nuclei were counted using 0.4% Trypan Blue and hemocytometer to determine the number of cells. (ThermoFisher, 15-250-061) and hemocytometer. Nuclei were fixed and processed using the manufacturer’s instructions using the Evercode Nuclei Fixation v2 kit (Parse Biosciences, ECF2003). The libraries were prepared according to the manufacturers’ instructions (Parse Evercode V2). Libraries were sequenced using Nextseq2000 at the University of Texas at Dallas Genome Core. 2.2.1.2. Bioinformatics analysis The sequenced libraries were processed using Parse Biosciences Evercode Software Suite and adapted from prior methods [ 6 ]. Initial quality control steps included filtering out nuclei with low gene counts (10%), or suspected doublets to retain only high-quality nuclei, were processed using the Seurat pipeline. The resulting expression profiles were normalized and scaled to account for differences in sequencing depth. Dimensionality reduction was performed using clustering with Uniform Manifold Approximation and Projection (UMAP) for clustering and visualization of DRG subpopulations. To annotate the anchor study dataset, a prior published single-cell cross-species atlas from Bhuiyan et al. [ 6 ] was used to transfer levels using Conos. This established a detailed anchor set for rat DRG populations which was subsequently used as a reference to assign labels to a dataset containing both sham and injured animals. Differential expression analysis was conducted using the DEseq2 package within the labeled clusters between sham and injured animals. Data visualization was conducted in R. Data can be accessed at https://paincenter.utdallas.edu/sensoryomics/ . 2.3. Chronic low back pain study 2.3.1. Animal acclimation and husbandry All animal procedures followed the National Institutes of Health guidelines and the Public Health Service Policy on Humane Care and Use of Laboratory Animals and were approved by the IACUC at the University of Nebraska-Lincoln. Adult female Sprague Dawley rats (Envigo) were housed under a reverse 12-hour light/dark cycle with ad libitum access to food and water. Before the start of the study, each animal was acclimated to experimenters for a minimum of five minutes per day for eight days. 2.3.2. Surgical procedure to elicit chronic low back pain During surgery, rats were anesthetized using 2-3% isoflurane in oxygen and given a single subcutaneous injection of Buprenorphine SR (0.75 mg/kg) for post-operative pain management. The L5-L6 disc was located using the iliac crest as a landmark, and tissue dissection was performed according to previously established methods in our lab [ 40 ]. The abdomen was shaved and disinfected with betadine and iodine. A midline incision was made along the ventral abdomen at the level of the iliac crest, extending through the skin, subcutaneous tissue, and fascia to expose the retroperitoneal space while avoiding major blood vessels. The L5-L6 disc was bilaterally punctured using a 3 mm dissecting needle (Roboz, RS-6066) and swept in a 90° arc six times in the transverse plane, following an established protocol [ 40 ]. The abdominal incision was closed using a continuous subcuticular suture pattern. Animals were observed 1 hour after surgery and then every 12 hours for 3 days to monitor well-being and wound healing. Post-operative behavioral assessments began after a two-week recovery period. By the end of the study two sham animals were euthanized due to surgical complications or illness, making the final animal numbers per group: sham = 10, injured = 12. These animals made up the final pool from which animals were randomly selected for snRNAseq. 2.3.3. Animal health pain-like behavior assessments Prior to all behavioral testing, animals were acclimated to equipment for at least five minutes per animal per apparatus. Animals were randomly assigned to groups and experimenters were blinded to all animal assignments. Animals were housed under a reverse light cycle and all behavioral data was collected with red light to reduce stress. Health and well-being were monitored bi-weekly by measuring weight, and any weight loss exceeding 10% from the previous measurement prompted veterinary consultation. Pain-like behaviors were evaluated bi-weekly post-operatively using two evoked assays (grip strength and pressure algometry). Prior to the start of the study, two baseline collections per animal were averaged for each assay. 2.3.3.1. Grip strength Axial hypersensitivity was quantified using grip strength (San Diego Instruments, 2325-0063) following methods previously described by our lab [ 40 ]. Animals underwent acclimation to the testing environment for at least 5 minutes before each test. In short, animals were held by the hindquarters which allowed them to grip a mesh metal wire that was connected to a grip strength force sensor. After the animal gripped the bar, the animal hold was then shifted to the base of the tail and the animal was gently and uniformly pulled back from the bar over the course of 1-1.5 seconds until it released the mesh. Three trials were recorded, and the average maximum force (N) was used. Grip strength values were normalized to baseline measurements for each animal. 2.3.3.2. Pressure algometry Pressure sensitivity at the L5-L6 region was assessed using an electronic von Frey aesthesiometer (IITC, 2391) using a blunt probe following methods previously described by our lab [ 40 ]. In short, prior to testing, each animal was acclimated to the testing environment for at least 5 minutes. During testing, the animal was carefully wrapped in a clean cotton t-shirt and loosely restrained by one experimenter, while a second experiment applied the blunt probe to the skin over the L5-L6 area. The pressure gradually increased until the animal demonstrated a pain-like response, such as avoidance behavior or vocalization. Two measurements were taken for each animal, and the average was recorded. All pressure thresholds were normalized to the baseline measurements for each animal. 2.3.4. Motion segment processing Tissue was processed following methods previously described by our lab [ 40 ]. In short, the L5-L6 motion segments were carefully isolated from each animal, cleaned, and fixed in 3 mL of 4% paraformaldehyde (PFA, Sigma Aldrich, 441244-1KG) for 24 hours at room temperature, with continuous agitation at 180 rpm on an orbital shaker. After fixation, the motion segments were rinsed three times for 15 minutes each with 1X PBS and subsequently placed in 3 mL of Immunocal (Fisher Scientific, NC9044643) to decalcify for 18 hours at room temperature with constant agitation on an orbital shaker. Once decalcified, the motion segments were washed three times with 1X PBS for 15 minutes each and then soaked in a 30% sucrose solution (wt/v) prepared in 1X PBS for 24 hours at 4°C. The segments were cryoembedded in Optimal Cutting Temperature Compound (Scigen 4586), frozen at −80°C, and sectioned in the sagittal plane into 40 μM slices. One injured animal was excluded from analysis due to issues with tissue harvesting, making the final animal numbers for subsequent motion segment histological analyses and behavioral analyses: sham = 10, injured = 11. 2.3.5. Immunohistochemistry Immunohistochemistry (IHC) was performed to assess the degree of nerve innervation in the disc. A peptidergic C-fiber marker for calcitonin gene-related peptide (CGRP) was selected due to prior demonstrated efficacy of staining nerve fibers in the disc [ 2 , 40 , 48 , 57 ]. Three 40 μM thick sections of the L5-L6 motion segments from each animal were post-fixed in 4% paraformaldehyde (PFA) for 15 minutes, followed by two 5-minute washes in 1X PBS and a final 5-minute wash in 1X PBST. The sections were then blocked for 1 hour in a blocking buffer consisting of 1X PBST, 3% goat serum (Sigma-Aldrich, G9023), and 0.01% Tween-20 (Fisher Scientific, BP337-100). Following blocking, the sections were then incubated overnight (16 hours) with a primary antibody diluted in blocking buffer, using 1:100 mouse anti-CGRP (Abcam, Ab81887). After incubation, the slides were washed three times for 15 minutes each with PBST, then incubated for 2 hours with a secondary antibody 1:500 goat anti-mouse AF488 (Abcam, Ab150117) diluted in blocking buffer. Following secondary incubation, the sections were washed three additional times with PBST for 15 minutes each, followed by a 10-minute incubation with 1:1000 DAPI (ThermoFisher, D1306) in 1X PBS for nuclear staining. Finally, the sections were mounted using Prolong Gold antifade reagent (ThermoFisher, P36930). Slides were imaged using a Zeiss LSM 800 confocal microscope (Carl Zeiss Microscopy) with a 10X objective. Whole motion segments were captured using the tiling feature and z-stacks to visualize the entire section thickness. Both brightfield and fluorescent images were taken with the following excitation/emission settings: 353/465 nm for DAPI and 493/517 nm for AF488. 2.3.6. Image analysis Immunohistochemistry images were first processed using ZEN Blue (Carl Zeiss Microscopy), where Z-stack projections were generated and saved as .czi files. These files were then imported into QuPath-0.4.3 (Pete Bankhead) for further analysis and annotation of disc regions. The wand and brush tools were used to create annotations on each section, representing distinct disc regions, including the dorsal ligament, dorsal annulus fibrosus (AF), nucleus pulposus (NP), the inner two-thirds of the ventral AF, the outer one-third of the ventral AF, granulation tissue, and ventral ligament. After completing the annotations, the ‘Pixel Classifier’ tool was applied to threshold CGRP-positive staining within each region. The percentage of positive staining was calculated by normalizing the CGRP-positive area to the total area of the respective disc region. 2.3.7. Statistical analysis Behavioral data are presented as mean ± 95% confidence interval (CI). IHC data are presented as mean ± standard deviation. All data were analyzed using GraphPad Prism V9. Normality was assessed using a Shapiro-Wilk test to determine subsequent analyses. Nerve IHC were analyzed using a Kruskal-Wallis test with Dunn’s post hoc. Behavioral data were analyzed using a Two-Way ANOVA with Dunnett’s post hoc test. Results were considered statistically significant when p≤0.05. Differential gene expression analysis was performed between sham and injured animals, looking at genes averaged across all neuronal clusters and within specific clusters. Pseudobulked count matrices were generated by summing raw counts per gene across all neuronal nuclei within each sample before applying DESeq2 normalization and modeling. Preprocessing included quality filtering and normalization to account for variations in sequencing depth. Genes were considered significantly upregulated in injured animals with a Log Fold Change (LogFC) threshold of 1.35, and p-adjusted ≤0.05. Differential expression analysis between sham and injured animals within each defined cluster was conducted using a non-parametric Wilcoxon rank-sum test. Comparisons between select genes across neuronal clusters were conducted using a non-parametric Mann-Whitney U test. Spearman correlations were calculated between gene expression levels and pain-like behavior metrics (grip strength and pressure algometry) at week 15 as well as CGRP-positive nerve staining in the disc. 3. Results 3.1. Anchor study 3.1.1. snRNAseq and clustering The anchor set consisted of 13,346 nuclei, which were clustered into 16 distinct clusters ( Fig. 2A ). High-resolution snRNA-seq metrics demonstrated robust sequencing quality, with a median of 1,830 transcripts per cell, a median of 1,155 genes per cell, and a mean of 97,967 reads per cell ( Fig. 2B-C ). The average mitochondrial RNA content across samples remained below 10%, suggesting minimal cell damage or lysis during processing. These metrics indicate the data had sufficient sequencing depth and quality for downstream analysis. Neuron-specific genes were examined, and strong expressions of both Snap25 , and Rbfox3 (NeuN) were observed across neuronal clusters, with lower expression in Schwann cell cluster 2, which expresses Schwann cell markers ( Fig. 2D, E ). Dot plot analyses revealed the percentage of cells expressing various genes across clusters, highlighting specific expression patterns. Ntrk3 was highly expressed in clusters 4 and 13, Trpm8 was predominantly expressed in cluster 15, and Il31ra showed strong expression in cluster 14, along with their corresponding average expression levels ( Fig. 2F ). UMAP projections highlighted spatial distribution and distinct expression patterns of genes such as Calca, Il31ra, Ntrk2, and Gfra2 in specific cluster populations ( Fig. 2G ). The clusters were confidently annotated using our previously published dataset, which includes genes from mouse species [ 6 ]. Download figure Open in new tab Figure 2. DRG sequencing and clustering analysis in the anchor study to establish subpopulation labels (A) UMAP plot of predicted clusters for T13-L1 bilateral DRGs harvested from male and female Sprague Dawley rats (n=6 male and female). Quality control metrics for RNA sequencing in the anchor dataset: (B) Volin Plot showing total number of molecules detected within a cell (nCount_RNA) and (C) Volin Plot showing number of genes detected in each cell (nFeature_RNA). (D) Volin Plot showing the expression levels of the neuronal marker Snap25 and (E) Rbfox3 (neuronal nuclei) across clusters. (F) Dot plots showing the percentage of cells expressing key genes across clusters. UMAP showing gene expression markers (G). Uniform Manifold Approximation and Projection; (UMAP). N=6 rats per group. Using Conos label transfer, eighteen neuronal subtypes were identified across 13,346 nuclei ( Fig. 3A ). These annotations were based on the harmonized atlas, of which 15 neuronal subtypes have been described as having at least one known function in rodent models [ 6 ]. UMAP projections revealed distinct clustering of nuclei, with populations showing similar proportions in both mouse and rat DRG datasets ( Fig. 3B-C ). The subpopulations include A-fiber subtypes such as Pvalb-expressing proprioceptors, Ntrk3high+Ntrk2 Aβ-rapid-adapting (RA) low-threshold mechanoreceptors (LTMRs), Ntrk3high+S100a16 Aβ-slow-adapting (SA) LTMRs, Ntrk3low+Ntrk2 Aδ-LTMRs, Calca+Bmpr1b Aδ-high-threshold mechanoreceptors (HTMRs), and Calca+Smr2 Aδ-HTMRs. The C-fiber subtypes include Calca+Sstr2, Calca+Adra2a, and Calca+Dcn nociceptors; Calca+Oprk1, a nociceptor expressing the opioid receptor kappa 1; Trpm8 cold nociceptors; and Th C-LTMRs. Additional C-fiber subtypes include Mrgprd nociceptors, Mrgpra3+Mrgprb4 C-LTMRs, Mrgpra3+Trpv1 pruriceptors, and Sst pruriceptors. Notably, Rxfp1, a rare subtype previously described to express high levels of Trpv1 [ 6 ], was also identified in our anchor dataset. Lastly, a small representation of an Atf3 cluster was also identified, characterized by the expression of this neuronal injury-associated gene. While the clusters represent shared cell types across species, specific clusters exhibited unique distributions, as shown in the pie chart, reflecting species-specific transcription patterns ( Fig. 3C-D ). Notably, the proportion of Th C-LTMRs was higher in mice (13.37%) compared to rats (8.68%), whereas Pvalb-expressing proprioceptors were more prevalent in rats (7.96%) than in mice (4.06%), but these differences may also be due to the level of DRGs assessed. In addition, the proportion of Ntrk3high+Ntrk2 Aβ-RA LTMRs was higher in rats (4.88%) compared to mice (2.81%). Dot plot analyses comparing the anchor dataset and the reference dataset demonstrated the expression profiles of key genes across clusters in both species ( Fig. 3D-E ). Trpm8 , associated with cold sensing, was predominantly expressed in specific clusters in both mouse and rat, showing similar gene expression patterns across species. Pvalb, a marker of proprioceptive neurons, exhibited consistent cluster enrichment patterns alongside genes such as Esrng , Vsnl1 , and Tafa2 in both species. Similarly, the peptidergic cluster population Calca+Bmpr1b displayed unique clustering patterns, characterized by distinct expression of genes including Tafa1 , Kcnq5 , and Csmd3 . These findings highlight both conserved and species-specific transcriptional features in neuronal subtypes. Next, using a mouse reference dataset, Bhuiyan et al. [ 6 ] and https://paincenter.utdallas.edu/sensoryomics/sensoryomics-dashboard/ [ 63 ], we identified 49 rat-specific genes distinguishing rat-specific transcriptional features. These genes were present in this rat dataset but were absent from the mouse reference dataset ( Fig. 3F ). These genes included Abcg3l1 ( ATP-binding cassette, subfamily G) , Dennd2b (DENN domain containing 2B) , Dsp (Desmoplakin) a gene involved in cell adhesion , Kcne2 (Potassium Channel, Voltage-Gated Subfamily E Regulatory Beta Subunit 2) , Lhfpl6 (Tetraspan subfamily member 6) , Nrdc (Nardilysin), a metallopeptidase that regulates growth factors and inflammatory response, and Tlcd4 (TLC Domain-Containing 4), which were not detected in the mouse dataset [ 6 ]. Download figure Open in new tab Figure 3. Species-Specific and Conserved Neuronal Subtypes in DRG ( A ) UMAP showing 18 neuronal subtypes after Conos label transfer with labelled clusters of neuronal subtypes. ( B-C) Species-specific differences in neuronal subtype proportions between mouse and rat. ( D-E ) Dot plot analyses demonstrated conserved gene expression profiles across mouse and rat for clusters. ( F ) Comparative analyses identified 49 rat-specific genes absent in mouse datasets, including Abcg3l1, Dennd2b, Dsp , and Kcne2 , which are associated with cell adhesion, structural integrity, and inflammatory responses. Six rats (n=3 male and n=3 female) were harvested from both the University of Nebraska-Lincoln and the University of Texas at Dallas by two unique experimenters for a total of 12 rats. Uniform Manifold Approximation and Projection; (UMAP). 3.2. Chronic LBP study 3.2.1. snRNAseq and clustering Libraries were generated from the bilateral T13-L1 dorsal root ganglia (DRG) of a chronic low back pain cohort, which included 12 female rats—6 from the injured group and six from the sham group. A total of 39,058 nuclei were recovered across both conditions. Nuclei had a median of 1,162 transcripts, and the sequencing depth was an average of 16,395 reads per nuclei. Conos label transfer was applied to the chronic lower back pain cohort using the previously annotated anchor set. UMAP projections again identified 18 neuronal subpopulations. The cluster subtypes showed broad similarities between the injured and sham groups ( Fig. 4A-B ). Cluster proportions between sham and injured rats were compared ( Fig. 4C ). While most clusters were consistently represented in both groups, we observed slight increases in the proportions of injury-associated clusters in the injured animals. For example, the Ntrk3low+Ntrk2 (Aδ-LTMRs) cluster showed a cell proportion of 3.2% in the injured group compared to 2.3% in the sham group, reflecting injury-induced transcriptional changes. Download figure Open in new tab Figure 4. snRNAseq atlas comparing sham and injured DRG subpopulations in a chronic LBP model ( A ) UMAP projection showing 18 distinct clusters identified from the anchor dataset; labels were transferred to the chronic LBP dataset and consistent across sham and injured animals. ( B ) Clusters annotated using the Bhuiyan et al. [ 6 ] dataset for both anchor and chronic LBP datasets, showing 18 defined clusters. ( C ) Relative proportions of each DRG subpopulation in sham and injured animals. ( D ) All differentially expressed genes by cluster. N=6 animals per group. Uniform Manifold Approximation and Projection; (UMAP). Differential gene expression analysis was performed between sham and injured animals across all neuronal clusters as well as within specific clusters. DEGs were identified between sham and injured animals across clusters ( Fig. 4D ). This approach inherently dilutes cluster-specific signals, as variations in gene expression within specific subpopulations may be masked by the larger, heterogeneous population. The highest proportion of DEGs was observed in the Ntrk3low+Ntrk2 (Aδ-LTMRs), Calca+Bmpr1b (Aδ-HTMRs), and Pvalb (A-fiber proprioceptors) populations. In the injured group, several pain-associated genes showed significant upregulation compared to sham log2FC of 0.44 or higher. For instance, within the Ntrk3low+Ntrk2 population (Aδ-LTMRs) cluster, Scn9a and Piezo2 , key ion channels involved in pain sensitivity, were significantly upregulated in injured animals. Specifically, Scn9a exhibited a log2 fold change of 0.58 while Piezo2 showed a log2 fold change of 0.44. The Ntrk3low+Ntrk2 population (Aδ-LTMRs) are known to innervate hair follicles and detect hair deflection [ 1 ] but have no known role in disc innervation. Additionally, Ngfr (nerve growth factor receptor), a marker linked to neurotrophic signaling in pain pathways, displayed significantly increased expression within this cluster in injured animals with a log2FC of 0.66. Other genes associated with disc degeneration and cellular senescence were also significantly elevated in the injured group in this cluster. Importantly, Spock2 (sparc/osteonectin, cwcv and kazal-like domains proteoglycan 2), related to extracellular matrix integrity and implicated in disc degeneration and Fstl1 , related to nucleus pulposus cell senescence and intervertebral disc degeneration were significantly increased in the Ntrk3low+Ntrk2 population (Aδ-LTMRs) ( Fig 4D ). These data suggest injured animals may exhibit alterations in genes associated with pain signaling and disc degeneration mechanisms. The highest proportion of differentially expressed genes were found in the Ntrk3low+Ntrk2 (Aδ-LTMRs), Calca+Bmpr1b (Aδ-HTMRs), and Pvalb (A-fiber proprioceptors) populations ( Fig. 4D ). Consistent DRG clusters across sham and injured animals suggest that injury does not alter overall DRG subpopulations but affects gene expression within existing neuronal subtypes. Differential expression analysis was conducted across all neuronal clusters to identify transcriptional changes following injury. A total of 75 representative DEGs with an adjusted p-value (p-adj) < 0.05 are shown ( Fig. 5A ), selected from a larger set of significant DEGs identified across all neuronal clusters. These genes were categorized into functional groups, including neuronal cytoskeleton and axonal transport, ion channels and transporters, chaperones and heat shock proteins, signal transduction and synaptic proteins, metabolism and enzymes, extracellular matrix and adhesion, and other. Expression patterns of selected genes of interest are shown across neuronal clusters ( Fig. 5B–E ). Ngfr expression was broadly increased in injured animals, particularly in Pvalb (A-fiber proprioceptors) and Ntrk3high+S100a16 (Aβ-SA LTMRs) clusters ( Fig. 5B ). Scn9a expression was elevated following injury, particularly in Trpm8 (cold nociceptors) and Mrgpr3+Mrgprb4 (C-LTMRs) clusters ( Fig. 5C ). Fstl1 , an extracellular matrix-related gene, exhibited cluster-specific injury-induced upregulation, with the highest expression observed in Pvalb (A-fiber proprioceptors) and Mrgpr3+Mrgprb4 (C-LTMRs) clusters ( Fig. 5D ). Similarly, Spock2 showed distinct cluster-dependent changes, with the most pronounced increases in Pvalb (A-fiber proprioceptors) and Ntrk3low+Ntrk2 (Aδ-LTMRs) clusters ( Fig. 5E ). These data suggest injury-induced transcriptional adaptations across neuronal subtypes, highlighting potential mechanisms contributing to altered neuronal function and extracellular matrix remodeling in response to injury. Download figure Open in new tab Figure 5. Differential expression of key genes across neuronal clusters in sham and injured animals (A) Gene classification table grouping differentially expressed genes (p-adj <0.05) by functional category. Showing 75 DEGs from a larger pool. (B-E) Violin plots of Ngfr (B) , Scn9a (C) , Fstl1 (D) , and Spock2 (E) expression across neuronal clusters, comparing sham (black) and injured (red) animals. Expression levels are log-transformed (log2) for visualization. N=6 animals. P-adj value across clusters per gene displayed on each violin plot. 3.2.2. Differentially expressed pain-related genes correlate with pain-like behavior Animals developed progressive axial hypersensitivity, emerging around week 8 and remaining statistically significant compared to sham controls through week 15, with injured animals exhibiting grip strength reduced to ≤82% of baseline ( Fig. 6A and B ). Significant correlations were identified between grip strength and expression levels of Scn9a (p=0.0102), Ngfr (p=0.0275), Fslt1 (p=0.0025), and Spock2 (p=0.0205) across all clusters, suggesting these genes may play a role in the maintenance of axial hypersensitivity. A trend was also observed between grip strength and Cdkn2b (p=0.067) ( Fig. 6C ) across all clusters. Download figure Open in new tab Figure 6. Spearman correlations of differentially expressed genes with axial hypersensitivity in sham and injured animals measured by grip strength (A) Schematic of the grip strength assay to evaluate axial hypersensitivity. (B) Normalized grip strength at week 15 post-injury, showing significantly lower grip strength in injured animals indicating axial hypersensitivity. (C) Spearman correlations between differentially expressed genes relevant to pain-like behavior and disc degeneration ( Scn9a, Ngfr, Fstl1, Spock2, and Cdkn2b ) and grip strength. Significant correlations were found for Scn9a , Ngfr , Fstl1 , and Spock2 , suggesting these genes may contribute to persistent axial hypersensitivity. A trend toward significance was observed for Cdkn2b (p=0.067). N=6 animals per group. Black symbols represent sham, red symbols represent injured animals. *=p≤0.05. The arrow indicates that higher grip strength values correspond to less pain-like behavior, reflecting a value closer to the baseline (1). Significant differences in pressure algometry were observed between sham and injured animals from weeks 8 to 15, with injured animals exhibiting pressure thresholds reduced to ≤81% of baseline ( Fig 7A-B ); however, no significant correlations or trends were observed between pressure algometry and differentially expressed genes across all clusters ( Fig. 7C ). Together, these data suggest that axial hypersensitivity is a prominent pain phenotype in this discogenic LBP model with correlations to pain-related genes in the DRG. In addition, localized touch sensitivity, measured through pressure algometry, showed a significant decrease in injured animals suggesting a pain phenotype present in our animals; however, no significant correlations with DRG changes were observed with pressure algometry. Download figure Open in new tab Figure 7. Spearman correlations of differentially expressed genes of interest between sham and injured animals with pain-like behavior measured by pressure algometry (A) Schematic of the pressure algometry assay to evaluate pain-like behavior. (B) Normalized pressure algometry at week 15 post-injury, showing significantly lower pressure algometry in injured animals indicating pain-like behavior. (C) Spearman correlations between differentially expressed genes relevant to pain-like behavior and disc degeneration ( Scn9a, Ngfr, Fslt1, Spock2, and Cdkn2b ) and pressure algometry. No significance was observed, suggesting that gene expression changes does not correlate with pain-like behavior as measured by pressure algometry in this dataset. N=6 animals per group. Black symbols represent sham, red symbols represent injured animals. *=p ≤ 0.05. The arrow indicates that higher pressure algometry values correspond to less pain-like behavior, reflecting a value closer to the baseline (1). 3.2.3. Differentially expressed pain-related genes do not correlate with CGRP nerves in the disc Injured discs showed a significant increase in CGRP+ nerve staining compared to sham discs, indicating heightened peptidergic nociceptor innervation post-injury ( Fig. 7A-B ). Specifically, the dorsal annulus fibrosus (AF) exhibited a 5.96-fold increase, the nucleus pulposus (NP) showed a 2.08-fold increase, and the ventral AF demonstrated a 2.79-fold increase in CGRP+ staining ( Fig. 8B ). However, no significant correlations were observed between CGRP staining across the disc and the differentially expressed genes across all clusters ( Fig. 8C ), suggesting that while nociceptor density increases with injury, gene expression profiles in the DRG may not be associated with nociceptive nerve density in the disc. Download figure Open in new tab Figure 8. Spearman correlations of differentially expressed genes of interest between sham and injured animals relative to CGRP-positive nerve staining in the disc (A) Diagram of a sham disc segmented into regions of interest for analysis, labeled as follows: a = dorsal ligament, b = dorsal annulus fibrosus (AF), c = nucleus pulposus (NP), d = ventral annulus fibrosus inner 2/3, e = ventral annulus fibrosus outer 1/3, f = granulation, g = ventral ligament. Dorsal and ventral directions of the disc are shown. (B) Quantification of CGRP-positive area in the disc at week 15, segmented by disc region. Data show significantly increased CGRP+ staining in the dorsal AF, NP, and ventral AF regions in injured discs, indicating increased peptidergic nociceptors are present post-injury. (C) Spearman correlations between pain- and degeneration-related genes ( Scn9a, Ngfr, Fstl1, Spock2, Cdkn2b ) and CGRP-positive area revealed no significant correlations, suggesting gene expression changes do not align with increased CGRP-positive innervation in this dataset. N=6 animals per group. Black symbols represent sham, red symbols represent injured animals. Scale bar = 1 mm. 4. Discussion Discogenic LBP is a complex and persistent condition with limited clinically effective long-term treatments. A key feature of this condition is the abnormal growth of nociceptor axons into previously aneural disc regions [ 7 , 20 ]. Repeated stimulation of nociceptors by noxious mechanical and inflammatory stimuli can elicit plastic changes, as observed in various pain states [ 26 ]. These plastic changes include increased expression and presentation of ion channels, enhanced excitability, and a lowered activation threshold [ 19 , 26 , 29 ]. Despite the prevalence of discogenic LBP, the lack of effective treatments could potentially be attributed to an incomplete understanding of the underlying mechanisms. This study addresses this gap by generating the first comprehensive transcriptional map of DRG neurons in a rat model of discogenic LBP using snRNAseq. These data highlight correlations between changes in nociceptor gene expression, pain-like behaviors, and nerve innervation in the disc. Importantly, these findings reveal significant upregulation of key pain-sensing ion channels, neuroinflammatory markers, and disc degeneration markers in DRGs of LBP animals, suggesting that mechanical hypersensitivity and inflammatory signaling are central to the pain phenotype in this model. Building on prior literature, the following sections will explore specific genes of interest, such as Scn9a , Piezo2 , Fstl1 , Ngfr , and Spock2 and their potential role in this discogenic LBP model. Additionally, the implications of specific sensory neuron subtypes will be addressed, along with the limitations of this study and directions for future research. In this study, Scn9a , which encodes the Nav1.7 sodium channel, was significantly upregulated in DRGs of injured animals (p-adj≤0.05) across all clusters and significantly correlated with reduced grip strength, a measure of axial hypersensitivity (p=0.0102). This suggests that Scn9a may play a role in pain-like behavior in this discogenic associated LBP model. Previous research has demonstrated that upregulation of Scn9a lowers activation thresholds in nociceptors, facilitating pain transmission in various chronic conditions, including neuropathic pain and OA [ 17 ]. In both an experimental model of neuropathic pain and in human neuropathic patients, DRG Scn9a upregulation was linked to enhanced pain sensitivity [ 38 ], highlighting its relevance across multiple pain conditions. Taken together, these data highlight the key role that Scn9a may play in maintaining hypersensitivity across pain states and specifically in discogenic LBP. Our findings also demonstrated that Piezo2 , a key neuronal mechanosensitive ion channel [ 14 , 66 , 69 ], was significantly upregulated in DRGs of injured animals (p-adj≤0.05) across all clusters. This suggests a component of mechanical hypersensitivity in this model of discogenic LBP. However, Piezo2 did not significantly correlate with axial hypersensitivity as measured by grip strength (p=0.1267), indicating it may not directly drive the observed pain-like behaviors in this model. Prior studies have shown that Piezo2 transduces mechanical signals into neuronal activation and contributes to pain sensitivity, with a well-documented role in mechanical allodynia in preclinical pain models [ 51 , 54 ]. These prior studies suggest the role of Piezo2 as key driver in conditions characterized by mechanical pain, where nociceptor excitability is increased in response to mechanical stimuli. Piezo2 has an established role in various pain models [ 51 , 54 ], and showed a log2 fold change of 0.44 in our discogenic LBP model. Although no correlations were observed between overall Piezo2 expression and pain-like behaviors in this model, Piezo2 expression was increased in pain-sensing nerve subclasses, particularly in Ntrk3low+Ntrk2 Aδ-LTMRs. Therefore, Piezo2 may remain a target for further examination in discogenic LBP models. Our findings demonstrated significant upregulation of Fstl1 (Follistatin-like protein 1) (p-adj≤0.05) and Ngfr (p75NTR) (p-adj≤0.05) in DRGs of injured animals across all clusters, suggesting a role for neuroinflammation in this model. Notably, Fstl1 expression strongly correlated with axial hypersensitivity (p=0.0025), suggesting its potential contribution to pain-like behavior in this model. Fstl1 has been shown to be elevated in inflamed and degenerative disc tissues in a rabbit model of intervertebral disc degeneration [ 75 ]. Fstl1 promotes inflammation and nucleus pulposus cell senescence through the TLR4/NF-κB pathway, exacerbating disc degeneration [ 75 ], which suggests a key role in chronic LBP. Additionally, Ngfr , which encodes the p75 nerve growth factor receptor, plays a critical role in neurotrophic signaling and nociceptor plasticity [ 67 ]. NGF signaling through p75NTR and TrkA receptors has been shown to support nerve innervation and promote nociceptor sensitization in animal models of OA [ 3 ]. Increased Ngfr expression in our injured animals across all clusters suggests that neuroinflammation may also contribute to the persistent pain observed in our animals. These data highlight the potential role of inflammatory signaling in sustained nociceptor plasticity across various chronic pain states. Our findings revealed a significant upregulation of Spock2 (p≤0.05), a gene associated with extracellular matrix integrity and remodeling, across all clusters in injured DRGs of this discogenic LBP model ( Fig. 4D ). Further, Spock2 was significantly correlated with axial hypersensitivity in this model ( Fig. 5C ) (p=0.0205). In this model, the increased expression of Spock2 may reflect active disc degeneration and its effects on nociceptor plasticity within the DRGs. Spock2 has been identified as a potential biomarker for LBP patients. Notably, genome-wide association studies from the UK Biobank identified Spock2 as a genetic risk factor for low back pain, linking this gene to the pathophysiology of disc degeneration and pain in human patients [ 21 ]. Together, these findings suggest Spock2 may serve as both a biomarker of disc degeneration and discogenic LBP. This study highlights DRG subtype-specific changes in Ntrk3low+Ntrk2 Aδ-LTMRs and Calca+Sstr2 peptidergic C-fibers, which exhibited increased proportions in injured animals (3.3% and 2.3% expression in injured and sham animals, respectively). Ntrk3low+Ntrk2 Aδ-LTMRs are associated with hair follicle innervation and detect stimuli on the skin through Ntrk2 (brain derived neurotrophic factor receptor) and Ntrk3 (NT-3 receptor) [ 1 ]. A study by Dhandapani et al. [ 16 ], showed that mice lacking this Aδ-LTMRs of Ntrk2-positive neurons exhibit reduced touch sensitivity and fail to respond to mechanical stimulation after neuropathic injury, suggesting that these Aδ-fibers may contribute to the development of allodynia. These fibers showed significant upregulation of pain-relevant genes such as Scn9a , Piezo2 , and Ngfr , implicating these genes in heightened sensitivity observed in this chronic LBP model. Although their presence in painful intervertebral discs is unknown, their involvement in allodynia suggests a potential role in discogenic pain. Similarly, Calca+Sstr2 C-fibers (8.5% and 7.1% expression in injured and sham animals, respectively), showed elevated expression of Spock2 and Fstl1 , suggesting these subtypes play a key role in driving pain and structural degeneration in discogenic LBP. In the Calca+Smr2 C-fiber population, genes such as Ngfr , Calm1 , Kcna2 , and Calca were differentially expressed in injured animals, suggesting these genes may have a role nociceptor plasticity in this DRG neuron subtype. Notably, Calca , which encodes CGRP—a key mediator of pain and inflammation [ 72 ]— showed increased expression. This was further supported by elevated CGRP+ nerve staining in the disc, indicating heightened nociceptor innervation post-injury, particularly within the dorsal annulus fibrosus. There are some limitations to our study. First, our findings are based on a rat model. Although our LBP model recapitulates many aspects of human LBP, it does not fully capture the complexity of human DRG structures or nociceptor diversity [ 5 , 64 ]. Importantly, when annotating the DRG dataset for this study it was observed that the rat DRGs were more closely related to mice than humans. During analysis, 18 neuronal subtypes were identified across mouse and rat DRG datasets, using Conos label transfer. This study highlights both conserved and species-specific transcriptional profiles. Notably, Ntrk3high+Ntrk2 Aβ-RA LTMRs were more prevalent in rats, while Th C-LTMRs were more prevalent in mice, suggesting interspecies differences in sensory neuron subpopulations. The rat dataset revealed 49 unique genes absent in mice, including Dsp (cell adhesion and structural integrity) and Nrdc (regulation of growth factors and inflammatory response), which may contribute to species-specific injury and repair mechanisms. These differences underscore the need to consider species-specific neuronal distributions and gene expression patterns in pain studies. An additional limitation is that analyzed gene expression was analyzed at a single time point (15 weeks post-injury), capturing chronic stage molecular profiles but potentially missing earlier dynamic changes in DRG plasticity. Longitudinal studies capturing gene expression at multiple time points would provide insights into the progression of nociceptor changes and could reveal potential intervention points [ 76 ]. In our study, a small percentage of Atf3 cells (0.13%) were observed. Studies of various pain models have examined significant changes in Atf3 , a marker of nerve injury, in timepoints up to 28 days [ 9 , 46 ]. However, due to our extended timepoint of 15 weeks post-injury, the transient expression changes in Atf3 or other acute markers of neuronal injury and stress may be missed. Lastly, no interventional studies were conducted to examine the potential role of genes that showed altered gene expression in pain in this LBP model. This can be assessed in future studies, in particular for Nav1.7. Along those lines, the findings highlight the potential of DRG-targeted interventions. Another opportunity would be a nociceptor-specific Piezo2 knockout model like that used in experimental OA demonstrated significant reductions in pain-like behaviors and joint damage [ 54 ]. These results suggest that similar strategies could be adapted for discogenic LBP, with a focus on targeting DRG-specific mechanisms to alleviate pain and improve clinical outcomes. Since these models are more amenable to using rats, genetic tools for targeting sensory neurons in this species would be useful. To conclude, this study provides insights into the molecular changes driving chronic discogenic LBP by offering the first comprehensive transcriptional map of DRG neurons using snRNAseq in a rat model. This work establishes a foundational rat atlas of DRG subpopulations, identifying 18 unique clusters and emphasizing the importance of nociceptor subtype-specific changes in discogenic LBP. Key genes, including Scn9a , Ngfr , Fstl1 , and Spock2 , were significantly correlated with axial hypersensitivity. This work establishes a foundation for understanding the molecular mechanisms underlying chronic discogenic LBP, and highlights opportunities for developing targeted, non-opioid therapies to enhance therapeutic translation and patient outcomes. Conflict of interest SMC, IS, RAW, TJP have no conflicts of interest to disclose. Acknowledgments The authors thank Kayla Ney and Evie Barnett for their support with pressure algometry and disc histology and immunohistochemistry, respectively. The authors thank the Genome Core at the University of Texas at Dallas for technical support. This work was funded by the National Institutes of Health (R01AR080926) and supported by the NIH grant (NS065926) to the Center for Advanced Pain Studies at the University of Texas at Dallas. Sequencing data is available at GEO at accession number GSE289659. References [1]. ↵ Abraira VE , Ginty DD . The sensory neurons of touch . Neuron 2013 ; 79 ( 4 ): 618 – 639 . OpenUrl CrossRef PubMed Web of Science [2]. ↵ Aoki Y , Takahashi Y , Ohtori S , Moriya H , Takahashi K . Distribution and immunocytochemical characterization of dorsal root ganglion neurons innervating the lumbar intervertebral disc in rats: a review . Life Sci 2004 ; 74 ( 21 ): 2627 – 2642 . OpenUrl CrossRef PubMed Web of Science [3]. ↵ Ashraf S , Bouhana KS , Pheneger J , Andrews SW , Walsh DA . Selective inhibition of tropomyosin-receptor-kinase A (TrkA) reduces pain and joint damage in two rat models of inflammatory arthritis . Arthritis Res Ther 2016 ; 18 ( 1 ): 97 . OpenUrl CrossRef PubMed [4]. ↵ Balague F , Mannion AF , Pellise F , Cedraschi C . Non-specific low back pain . Lancet 2012 ; 379 ( 9814 ): 482 - 491 . OpenUrl CrossRef PubMed Web of Science [5]. ↵ Bhuiyan SA , Xu M , Yang L , Semizoglou E , Bhatia P , Pantaleo KI , Tochitsky I , Jain A , Erdogan B , Blair S , Cat V , Mwirigi JM , Sankaranarayanan I , Tavares-Ferreira D , Green U , McIlvried LA , Copits BA , Bertels Z , Del Rosario JS , Widman AJ , Slivicki RA , Yi J , Sharif-Naeini R , Woolf CJ , Lennerz JK , Whited JL , Price TJ , Renthal W . Harmonized cross-species cell atlases of trigeminal and dorsal root ganglia . Science Advances 2024 ; 10 ( 25 ): eadj9173 . OpenUrl CrossRef PubMed [6]. ↵ Bhuiyan SA , Xu M , Yang L , Semizoglou E , Bhatia P , Pantaleo KI , Tochitsky I , Jain A , Erdogan B , Blair S , Cat V , Mwirigi JM , Sankaranarayanan I , Tavares-Ferreira D , Green U , McIlvried LA , Copits BA , Bertels Z , Del Rosario JS , Widman AJ , Slivicki RA , Yi J , Woolf CJ , Lennerz JK , Whited JL , Price TJ , Gereau RWt , Renthal W. Harmonized cross-species cell atlases of trigeminal and dorsal root ganglia . bioRxiv 2023 . [7]. ↵ Binch AL , Cole AA , Breakwell LM , Michael AL , Chiverton N , Creemers LB , Cross AK , Le Maitre CL . Nerves are more abundant than blood vessels in the degenerate human intervertebral disc . Arthritis Res Ther 2015 ; 17 : 370 . [8]. ↵ Brantigan JW , Neidre A , Toohey JS . The Lumbar I/F Cage for posterior lumbar interbody fusion with the variable screw placement system: 10-year results of a Food and Drug Administration clinical trial . Spine J 2004 ; 4 ( 6 ): 681 – 688 . OpenUrl CrossRef PubMed [9]. ↵ Braz JM , Basbaum AI . Differential ATF3 expression in dorsal root ganglion neurons reveals the profile of primary afferents engaged by diverse noxious chemical stimuli . Pain 2010 ; 150 ( 2 ): 290 – 301 . OpenUrl CrossRef PubMed Web of Science [10]. ↵ Burke JG , Watson RW , McCormack D , Dowling FE , Walsh MG , Fitzpatrick JM . Intervertebral discs which cause low back pain secrete high levels of proinflammatory mediators . J Bone Joint Surg Br 2002 ; 84 ( 2 ): 196 – 201 . OpenUrl CrossRef PubMed [11]. ↵ Cheung KM , Karppinen J , Chan D , Ho DW , Song YQ , Sham P , Cheah KS , Leong JC , Luk KD . Prevalence and pattern of lumbar magnetic resonance imaging changes in a population study of one thousand forty-three individuals . Spine (Phila Pa 1976) 2009 ; 34 ( 9 ): 934 – 940 . OpenUrl CrossRef [12]. ↵ Clauw DJ , Williams D , Lauerman W , Dahlman M , Aslami A , Nachemson AL , Kobrine AI , Wiesel SW . Pain sensitivity as a correlate of clinical status in individuals with chronic low back pain . Spine (Phila Pa 1976) 1999 ; 24 ( 19 ): 2035 – 2041 . OpenUrl CrossRef [13]. ↵ Collaborators GBDO . Global, regional, and national burden of osteoarthritis, 1990-2020 and projections to 2050: a systematic analysis for the Global Burden of Disease Study 2021 . Lancet Rheumatol 2023 ; 5 ( 9 ): e508 – e522 . OpenUrl CrossRef PubMed [14]. ↵ Coste B , Mathur J , Schmidt M , Earley TJ , Ranade S , Petrus MJ , Dubin AE , Patapoutian A . Piezo1 and Piezo2 are essential components of distinct mechanically activated cation channels . Science 2010 ; 330 ( 6000 ): 55 – 60 . OpenUrl Abstract / FREE Full Text [15]. ↵ de Schepper EI , Damen J , van Meurs JB , Ginai AZ , Popham M , Hofman A , Koes BW , Bierma-Zeinstra SM . The association between lumbar disc degeneration and low back pain: the influence of age, gender, and individual radiographic features . Spine (Phila Pa 1976) 2010 ; 35 ( 5 ): 531 – 536 . OpenUrl CrossRef [16]. ↵ Dhandapani R , Arokiaraj CM , Taberner FJ , Pacifico P , Raja S , Nocchi L , Portulano C , Franciosa F , Maffei M , Hussain AF , de Castro Reis F , Reymond L , Perlas E , Garcovich S , Barth S , Johnsson K , Lechner SG , Heppenstall PA . Control of mechanical pain hypersensitivity in mice through ligand-targeted photoablation of TrkB-positive sensory neurons . Nat Commun 2018 ; 9 ( 1 ): 1640 . OpenUrl CrossRef PubMed [17]. ↵ Dib-Hajj SD , Yang Y , Black JA , Waxman SG . The Na(V)1.7 sodium channel: from molecule to man . Nat Rev Neurosci 2013 ; 14 ( 1 ): 49 – 62 . OpenUrl CrossRef PubMed [18]. ↵ Fairbank J , Frost H , Wilson-MacDonald J , Yu LM , Barker K , Collins R , Spine Stabilisation Trial G. Randomised controlled trial to compare surgical stabilisation of the lumbar spine with an intensive rehabilitation programme for patients with chronic low back pain: the MRC spine stabilisation trial . BMJ 2005 ; 330 ( 7502 ): 1233 . OpenUrl Abstract / FREE Full Text [19]. ↵ Fingleton C , Smart K , Moloney N , Fullen BM , Doody C . Pain sensitization in people with knee osteoarthritis: a systematic review and meta-analysis . Osteoarthritis Cartilage 2015 ; 23 ( 7 ): 1043 – 1056 . OpenUrl CrossRef PubMed [20]. ↵ Freemont AJ , Peacock TE , Goupille P , Hoyland JA , O’Brien J , Jayson MI . Nerve ingrowth into diseased intervertebral disc in chronic back pain . Lancet 1997 ; 350 ( 9072 ): 178 – 181 . OpenUrl CrossRef PubMed Web of Science [21]. ↵ Freidin MB , Tsepilov YA , Palmer M , Karssen LC , Suri P , Aulchenko YS , Williams FMK , Group CMW . Insight into the genetic architecture of back pain and its risk factors from a study of 509,000 individuals . Pain 2019 ; 160 ( 6 ): 1361 – 1373 . OpenUrl CrossRef PubMed [22]. ↵ Fukui D , Kawakami M , Yoshida M , Nakao S , Matsuoka T , Yamada H . Gait abnormality due to spinal instability after lumbar facetectomy in the rat . Eur Spine J 2015 ; 24 ( 9 ): 2085 – 2094 . OpenUrl CrossRef PubMed [23]. ↵ Gane BM , Caparaso SM , Lee FS , Redwine AL , Wachs RA. Harvesting, Embedding, and Culturing Dorsal Root Ganglia in Multi-compartment Devices to Study Peripheral Neuronal Features . J Vis Exp 2024 ( 208 ). [24]. ↵ Garcia-Cosamalon J , del Valle ME , Calavia MG , Garcia-Suarez O , Lopez-Muniz A , Otero J , Vega JA . Intervertebral disc, sensory nerves and neurotrophins: who is who in discogenic pain? J Anat 2010 ; 217 ( 1 ): 1 – 15 . OpenUrl CrossRef PubMed [25]. ↵ Glaeser JD , Tawackoli W , Ju DG , Yang JH , Kanim LE , Salehi K , Yu V , Saidara E , Vit JP , Khnkoyan Z , NaPier Z , Stone LS , Bae HW , Sheyn D . Optimization of a rat lumbar IVD degeneration model for low back pain . JOR Spine 2020 ; 3 ( 2 ): e1092 . OpenUrl CrossRef [26]. ↵ Gold MS , Gebhart GF . Nociceptor sensitization in pain pathogenesis . Nat Med 2010 ; 16 ( 11 ): 1248 – 1257 . OpenUrl CrossRef PubMed [27]. ↵ Gore M , Sadosky A , Stacey BR , Tai KS , Leslie D . The burden of chronic low back pain: clinical comorbidities, treatment patterns, and health care costs in usual care settings . Spine (Phila Pa 1976) 2012 ; 37 ( 11 ): E668 – 677 . OpenUrl CrossRef [28]. ↵ Hadjipavlou AG , Tzermiadianos MN , Bogduk N , Zindrick MR . The pathophysiology of disc degeneration: a critical review . J Bone Joint Surg Br 2008 ; 90 ( 10 ): 1261 – 1270 . OpenUrl CrossRef PubMed [29]. ↵ He BH , Christin M , Mouchbahani-Constance S , Davidova A , Sharif-Naeini R . Mechanosensitive ion channels in articular nociceptors drive mechanical allodynia in osteoarthritis . Osteoarthritis Cartilage 2017 ; 25 ( 12 ): 2091 – 2099 . OpenUrl CrossRef PubMed [30]. ↵ Hestehave S , Allen HN , Gomez K , Duran P , Calderon-Rivera A , Loya-Lopez S , Rodriguez-Palma EJ , Khanna R . Small molecule targeting Na V 1.7 via inhibition of CRMP2-Ubc9 interaction reduces pain-related outcomes in a rodent osteoarthritic model . Pain 2025 ; 166 ( 1 ): 99 – 111 . OpenUrl CrossRef PubMed [31]. ↵ Hoy D , Bain C , Williams G , March L , Brooks P , Blyth F , Woolf A , Vos T , Buchbinder R . A systematic review of the global prevalence of low back pain . Arthritis Rheum 2012 ; 64 ( 6 ): 2028 – 2037 . OpenUrl CrossRef PubMed Web of Science [32]. ↵ Ivanova JI , Birnbaum HG , Schiller M , Kantor E , Johnstone BM , Swindle RW . Real-world practice patterns, health-care utilization, and costs in patients with low back pain: the long road to guideline-concordant care . Spine J 2011 ; 11 ( 7 ): 622 – 632 . OpenUrl CrossRef PubMed Web of Science [33]. ↵ Johnson WE , Caterson B , Eisenstein SM , Hynds DL , Snow DM , Roberts S . Human intervertebral disc aggrecan inhibits nerve growth in vitro . Arthritis Rheum 2002 ; 46 ( 10 ): 2658 – 2664 . OpenUrl CrossRef PubMed Web of Science [34]. ↵ Kepler CK , Markova DZ , Dibra F , Yadla S , Vaccaro AR , Risbud MV , Albert TJ , Anderson DG . Expression and relationship of proinflammatory chemokine RANTES/CCL5 and cytokine IL-1beta in painful human intervertebral discs . Spine (Phila Pa 1976) 2013 ; 38 ( 11 ): 873 – 880 . OpenUrl CrossRef [35]. ↵ Lai A , Ho L , Evashwick-Rogler TW , Watanabe H , Salandra J , Winkelstein BA , Laudier D , Hecht AC , Pasinetti GM , Iatridis JC . Dietary polyphenols as a safe and novel intervention for modulating pain associated with intervertebral disc degeneration in an in-vivo rat model . PLoS One 2019 ; 14 ( 10 ): e0223435 . OpenUrl CrossRef PubMed [36]. ↵ Lama P , Le Maitre CL , Harding IJ , Dolan P , Adams MA . Nerves and blood vessels in degenerated intervertebral discs are confined to physically disrupted tissue . J Anat 2018 ; 233 ( 1 ): 86 – 97 . OpenUrl CrossRef PubMed [37]. ↵ Leimer EM , Gayoso MG , Jing L , Tang SY , Gupta MC , Setton LA . Behavioral Compensations and Neuronal Remodeling in a Rodent Model of Chronic Intervertebral Disc Degeneration . Sci Rep 2019 ; 9 ( 1 ): 3759 . OpenUrl CrossRef PubMed [38]. ↵ Li Y , North RY , Rhines LD , Tatsui CE , Rao G , Edwards DD , Cassidy RM , Harrison DS , Johansson CA , Zhang H , Dougherty PM . DRG Voltage-Gated Sodium Channel 1.7 Is Upregulated in Paclitaxel-Induced Neuropathy in Rats and in Humans with Neuropathic Pain . J Neurosci 2018 ; 38 ( 5 ): 1124 – 1136 . OpenUrl Abstract / FREE Full Text [39]. ↵ Li Z , Liu H , Yang H , Wang J , Wang H , Zhang K , Ding W , Zheng Z . Both expression of cytokines and posterior annulus fibrosus rupture are essential for pain behavior changes induced by degenerative intervertebral disc: An experimental study in rats . J Orthop Res 2014 ; 32 ( 2 ): 262 – 272 . OpenUrl CrossRef PubMed [40]. ↵ Lillyman DJ , Lee FS , Barnett EC , Miller TJ , Alvaro ML , Drvol HC , Wachs RA . Axial hypersensitivity is associated with aberrant nerve sprouting in a novel model of disc degeneration in female Sprague Dawley rats . JOR Spine 2022 ; 5 ( 3 ): e1212 . OpenUrl CrossRef PubMed [41]. ↵ Liu Q , Wang X , Hua Y , Kong G , Wu X , Huang Z , Huang Z , Liu J , Yang Z , Zhu Q . Estrogen Deficiency Exacerbates Intervertebral Disc Degeneration Induced by Spinal Instability in Rats . Spine (Phila Pa 1976) 2019 ; 44 ( 9 ): E510 – E519 . OpenUrl CrossRef [42]. ↵ Livshits G , Popham M , Malkin I , Sambrook PN , Macgregor AJ , Spector T , Williams FM . Lumbar disc degeneration and genetic factors are the main risk factors for low back pain in women: the UK Twin Spine Study . Ann Rheum Dis 2011 ; 70 ( 10 ): 1740 – 1745 . OpenUrl Abstract / FREE Full Text [43]. ↵ Lu Y , Guzman JZ , Purmessur D , Iatridis JC , Hecht AC , Qureshi SA , Cho SK . Nonoperative management of discogenic back pain: a systematic review . Spine (Phila Pa 1976) 2014 ; 39 ( 16 ): 1314 – 1324 . OpenUrl CrossRef [44]. ↵ Luoma K , Riihimaki H , Luukkonen R , Raininko R , Viikari-Juntura E , Lamminen A . Low back pain in relation to lumbar disc degeneration . Spine (Phila Pa 1976) 2000 ; 25 ( 4 ): 487 – 492 . OpenUrl CrossRef [45]. ↵ Maher C , Underwood M , Buchbinder R . Non-specific low back pain . Lancet 2017 ; 389 ( 10070 ): 736 – 747 . OpenUrl CrossRef PubMed [46]. ↵ Matsuura Y , Ohtori S , Iwakura N , Suzuki T , Kuniyoshi K , Takahashi K . Expression of activating transcription factor 3 (ATF3) in uninjured dorsal root ganglion neurons in a lower trunk avulsion pain model in rats . Eur Spine J 2013 ; 22 ( 8 ): 1794 – 1799 . OpenUrl CrossRef PubMed [47]. ↵ Modic MT , Ross JS . Magnetic resonance imaging in the evaluation of low back pain . Orthop Clin North Am 1991 ; 22 ( 2 ): 283 – 301 . OpenUrl CrossRef PubMed [48]. ↵ Mohd Isa IL , Abbah SA , Kilcoyne M , Sakai D , Dockery P , Finn DP , Pandit A . Implantation of hyaluronic acid hydrogel prevents the pain phenotype in a rat model of intervertebral disc injury . Sci Adv 2018 ; 4 ( 4 ): eaaq0597 . OpenUrl FREE Full Text [49]. ↵ Morinaga T , Takahashi K , Yamagata M , Chiba T , Tanaka K , Takahashi Y , Nakamura S , Suseki K , Moriya H . Sensory innervation to the anterior portion of lumbar intervertebral disc . Spine (Phila Pa 1976) 1996 ; 21 ( 16 ): 1848 – 1851 . OpenUrl CrossRef [50]. ↵ Muralidharan A , Park TSW , Mackie JT , Gimenez LGS , Kuo A , Nicholson JR , Corradini L , Smith MT . Establishment and Characterization of a Novel Rat Model of Mechanical Low Back Pain Using Behavioral, Pharmacologic and Histologic Methods . Front Pharmacol 2017 ; 8 : 493 . [51]. ↵ Murthy SE , Loud MC , Daou I , Marshall KL , Schwaller F , Kuhnemund J , Francisco AG , Keenan WT , Dubin AE , Lewin GR , Patapoutian A . The mechanosensitive ion channel Piezo2 mediates sensitivity to mechanical pain in mice . Sci Transl Med 2018 ; 10 ( 462 ). [52]. ↵ Nouwen A , Cloutier C , Kappas A , Warbrick T , Sheffield D . Effects of focusing and distraction on cold pressor-induced pain in chronic back pain patients and control subjects . J Pain 2006 ; 7 ( 1 ): 62 – 71 . OpenUrl CrossRef PubMed Web of Science [53]. ↵ O’Neill S , Manniche C , Graven-Nielsen T , Arendt-Nielsen L . Association between a composite score of pain sensitivity and clinical parameters in low-back pain . Clin J Pain 2014 ; 30 ( 10 ): 831 – 838 . OpenUrl CrossRef PubMed [54]. ↵ Obeidat AM , Wood MJ , Adamczyk NS , Ishihara S , Li J , Wang L , Ren D , Bennett DA , Miller RJ , Malfait AM , Miller RE . Piezo2 expressing nociceptors mediate mechanical sensitization in experimental osteoarthritis . Nat Commun 2023 ; 14 ( 1 ): 2479 . OpenUrl CrossRef PubMed [55]. ↵ Ohtori S , Takahashi K , Moriya H . Existence of brain-derived neurotrophic factor and vanilloid receptor subtype 1 immunoreactive sensory DRG neurons innervating L5/6 intervertebral discs in rats . J Orthop Sci 2003 ; 8 ( 1 ): 84 – 87 . OpenUrl CrossRef PubMed [56]. ↵ Oichi T , Taniguchi Y , Soma K , Chang SH , Yano F , Tanaka S , Saito T . A Mouse Intervertebral Disc Degeneration Model by Surgically Induced Instability . Spine (Phila Pa 1976) 2018 ; 43 ( 10 ): E557 – E564 . OpenUrl CrossRef [57]. ↵ Ozawa T , Ohtori S , Inoue G , Aoki Y , Moriya H , Takahashi K . The degenerated lumbar intervertebral disc is innervated primarily by peptide-containing sensory nerve fibers in humans . Spine (Phila Pa 1976) 2006 ; 31 ( 21 ): 2418 – 2422 . OpenUrl CrossRef [58]. ↵ Parisien M , Khoury S , Chabot-Dore AJ , Sotocinal SG , Slade GD , Smith SB , Fillingim RB , Ohrbach R , Greenspan JD , Maixner W , Mogil JS , Belfer I , Diatchenko L . Effect of Human Genetic Variability on Gene Expression in Dorsal Root Ganglia and Association with Pain Phenotypes . Cell Rep 2017 ; 19 ( 9 ): 1940 – 1952 . OpenUrl CrossRef PubMed [59]. ↵ Parker LM , Murrell SE , Boden SD , Horton WC . The outcome of posterolateral fusion in highly selected patients with discogenic low back pain . Spine (Phila Pa 1976) 1996 ; 21 ( 16 ): 1909 – 1916 ; discussion 1916-1907. OpenUrl CrossRef [60]. ↵ Qaseem A , McLean RM , O’Gurek D , Batur P , Lin K , Kansagara DL , Clinical Guidelines Committee of the American College of P, Commission on Health of the P, Science of the American Academy of Family P , Cooney TG , Forciea MA , Crandall CJ , Fitterman N , Hicks LA , Horwitch C , Maroto M , McLean RM , Mustafa RA , Tufte J , Vijan S , Williams JW , Jr. Nonpharmacologic and Pharmacologic Management of Acute Pain From Non-Low Back, Musculoskeletal Injuries in Adults: A Clinical Guideline From the American College of Physicians and American Academy of Family Physicians . Ann Intern Med 2020 ; 173 ( 9 ): 739 – 748 . OpenUrl CrossRef PubMed [61]. ↵ Qi L , Iskols M , Shi D , Reddy P , Walker C , Lezgiyeva K , Voisin T , Pawlak M , Kuchroo VK , Chiu IM , Ginty DD , Sharma N . A mouse DRG genetic toolkit reveals morphological and physiological diversity of somatosensory neuron subtypes . Cell 2024 ; 187 ( 6 ): 1508 – 1526 e1516. OpenUrl CrossRef PubMed [62]. ↵ Raastad J , Reiman M , Coeytaux R , Ledbetter L , Goode AP . The association between lumbar spine radiographic features and low back pain: a systematic review and meta-analysis . Semin Arthritis Rheum 2015 ; 44 ( 5 ): 571 – 585 . OpenUrl CrossRef PubMed [63]. ↵ Ray PR , Shiers S , Caruso JP , Tavares-Ferreira D , Sankaranarayanan I , Uhelski ML , Li Y , North RY , Tatsui C , Dussor G , Burton MD , Dougherty PM , Price TJ . RNA profiling of human dorsal root ganglia reveals sex differences in mechanisms promoting neuropathic pain . Brain 2023 ; 146 ( 2 ): 749 – 766 . OpenUrl CrossRef PubMed [64]. ↵ Renthal W , Chamessian A , Curatolo M , Davidson S , Burton M , Dib-Hajj S , Dougherty PM , Ebert AD , Gereau RWt , Ghetti A , Gold MS , Hoben G , Menichella DM , Mercier P , Ray WZ , Salvemini D , Seal RP , Waxman S , Woolf CJ , Stucky CL , Price TJ . Human cells and networks of pain: Transforming pain target identification and therapeutic development . Neuron 2021 ; 109 ( 9 ): 1426 – 1429 . OpenUrl CrossRef PubMed [65]. ↵ Sameda H , Takahashi Y , Takahashi K , Chiba T , Ohtori S , Moriya H . Dorsal root ganglion neurones with dichotomising afferent fibres to both the lumbar disc and the groin skin. A possible neuronal mechanism underlying referred groin pain in lower lumbar disc diseases . J Bone Joint Surg Br 2003 ; 85 ( 4 ): 600 – 603 . OpenUrl CrossRef [66]. ↵ Shin KC , Park HJ , Kim JG , Lee IH , Cho H , Park C , Sung TS , Koh SD , Park SW , Bae YM . The Piezo2 ion channel is mechanically activated by low-threshold positive pressure . Sci Rep 2019 ; 9 ( 1 ): 6446 . OpenUrl CrossRef PubMed [67]. ↵ Sobue G , Yamamoto M , Doyu M , Li M , Yasuda T , Mitsuma T . Expression of mRNAs for neurotrophins (NGF, BDNF, and NT-3) and their receptors (p75NGFR, trk, trkB, and trkC) in human peripheral neuropathies . Neurochem Res 1998 ; 23 ( 6 ): 821 – 829 . OpenUrl CrossRef PubMed Web of Science [68]. ↵ Sun H , Sun Z , Xu X , Lv Z , Li J , Wu R , Fei Y , Tan G , Liu Z , Liu Y , Shi D . Blocking TRPV4 Ameliorates Osteoarthritis by Inhibiting M1 Macrophage Polarization via the ROS/NLRP3 Signaling Pathway . Antioxidants (Basel) 2022 ; 11 ( 12 ). [69]. ↵ Szczot M , Nickolls AR , Lam RM , Chesler AT . The Form and Function of PIEZO2 . Annu Rev Biochem 2021 ; 90 : 507 – 534 . OpenUrl CrossRef PubMed [70]. ↵ Teraguchi M , Yoshimura N , Hashizume H , Muraki S , Yamada H , Minamide A , Oka H , Ishimoto Y , Nagata K , Kagotani R , Takiguchi N , Akune T , Kawaguchi H , Nakamura K , Yoshida M . Prevalence and distribution of intervertebral disc degeneration over the entire spine in a population-based cohort: the Wakayama Spine Study . Osteoarthritis Cartilage 2014 ; 22 ( 1 ): 104 – 110 . OpenUrl CrossRef PubMed [71]. ↵ Vane JR , Botting RM . Anti-inflammatory drugs and their mechanism of action . Inflamm Res 1998 ; 47 Suppl 2 : S78 – 87 . OpenUrl CrossRef PubMed Web of Science [72]. ↵ Walsh DA , Mapp PI , Kelly S . Calcitonin gene-related peptide in the joint: contributions to pain and inflammation . Br J Clin Pharmacol 2015 ; 80 ( 5 ): 965 – 978 . OpenUrl CrossRef PubMed [73]. ↵ Wang Y , Videman T , Battie MC . ISSLS prize winner: Lumbar vertebral endplate lesions: associations with disc degeneration and back pain history . Spine (Phila Pa 1976) 2012 ; 37 ( 17 ): 1490 – 1496 . OpenUrl CrossRef [74]. ↵ Weiler C , Nerlich AG , Bachmeier BE , Boos N . Expression and distribution of tumor necrosis factor alpha in human lumbar intervertebral discs: a study in surgical specimen and autopsy controls . Spine (Phila Pa 1976) 2005 ; 30 ( 1 ): 44 – 53 ; discussion 54. OpenUrl CrossRef [75]. ↵ Yan X , Ding JY , Zhang RJ , Zhang HQ , Kang L , Jia CY , Liu XY , Shen CL . FSTL1 Accelerates Nucleus Pulposus Cell Senescence and Intervertebral Disc Degeneration Through TLR4/NF-kappaB Pathway . Inflammation 2024 ; 47 ( 4 ): 1229 – 1247 . OpenUrl CrossRef PubMed [76]. ↵ Zaki S , Smith MM , Little CB . Pathology-pain relationships in different osteoarthritis animal model phenotypes: it matters what you measure, when you measure, and how you got there . Osteoarthritis Cartilage 2021 ; 29 ( 10 ): 1448 – 1461 . OpenUrl CrossRef PubMed View the discussion thread. Back to top Previous Next Posted February 24, 2025. Download PDF Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Single-nuclei RNA Sequencing Reveals Distinct Transcriptomic Signatures of Rat Dorsal Root Ganglia in a Chronic Discogenic Low Back Pain Model Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Single-nuclei RNA Sequencing Reveals Distinct Transcriptomic Signatures of Rat Dorsal Root Ganglia in a Chronic Discogenic Low Back Pain Model Sydney M. Caparaso , Ishwarya Sankaranarayanan , David J. Lillyman , Theodore J. Price , Rebecca A. Wachs bioRxiv 2025.02.19.639130; doi: https://doi.org/10.1101/2025.02.19.639130 Share This Article: Copy Citation Tools Single-nuclei RNA Sequencing Reveals Distinct Transcriptomic Signatures of Rat Dorsal Root Ganglia in a Chronic Discogenic Low Back Pain Model Sydney M. Caparaso , Ishwarya Sankaranarayanan , David J. Lillyman , Theodore J. Price , Rebecca A. Wachs bioRxiv 2025.02.19.639130; doi: https://doi.org/10.1101/2025.02.19.639130 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Neuroscience Subject Areas All Articles Animal Behavior and Cognition (7624) Biochemistry (17651) Bioengineering (13873) Bioinformatics (41887) Biophysics (21424) Cancer Biology (18566) Cell Biology (25465) Clinical Trials (138) Developmental Biology (13365) Ecology (19871) Epidemiology (2067) Evolutionary Biology (24293) Genetics (15591) Genomics (22478) Immunology (17715) Microbiology (40331) Molecular Biology (17150) Neuroscience (88492) Paleontology (666) Pathology (2828) Pharmacology and Toxicology (4817) Physiology (7635) Plant Biology (15114) Scientific Communication and Education (2044) Synthetic Biology (4286) Systems Biology (9817) Zoology (2268)