Skip to main content
Download PDF

Chronic constriction injury-induced changes in circular RNA expression profiling of the dorsal root ganglion in a rat model of neuropathic pain

BMC Neuroscience volume 23, Article number: 64 (2022) Cite this article

Abstract

Background

The pathogenesis of neuropathic pain (NP) has not been fully elucidated. Gene changes in dorsal root ganglia (DRG) may contribute to the development of NP. Circular RNAs (circRNAs) are a class of endogenous noncoding RNAs that form covalently closed loop structures and are crucial for genetic and epigenetic regulation. However, little is known about circRNA changes in DRG neurons after peripheral nerve injury.

Methods

A sciatic nerve chronic constriction injury (CCI) model was established to induce neuropathic pain. We performed genome-wide circRNA analysis of four paired dorsal root ganglion (DRG) samples (L4–L5) from CCI and negative control (NC) rats using next-generation sequencing technology. The differentially expressed circRNAs (DEcircRNAs) were identified by differential expression analysis, and the expression profile of circRNAs was validated by quantitative PCR. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses were performed to predict the function of DEcircRNAs.

Results

A total of 374 DEcircRNAs were identified between CCI and NC rats using circRNA high-throughput sequencing. Among them, 290 were upregulated and 84 were downregulated in the CCI group. The expression levels of nine DEcircRNAs were validated by qPCR. Functional annotation analysis showed that the DEcircRNAs were mainly enriched in pathways and functions, including ‘dopaminergic synapse,’ ‘renin secretion,’ ‘mitogen-activated protein kinase signaling pathway,’ and ‘neurogenesis.’ Competing endogenous RNA analysis showed that the top 50 circRNAs exhibited interactions with four pain-related microRNAs (miRNAs). Circ:chr2:33950934–33955969 was the largest node in the circRNA–miRNA interaction network.

Conclusions

Peripheral nerve injury-induced neuropathic pain led to changes in the comprehensive expression profile of circRNAs in the DRG of rats. DEcircRNAs may advance our understanding of the molecular mechanisms underlying neuropathic pain.

Peer Review reports

Background

Neuropathic pain develops after injury of the somatosensory system [1, 2]. It is characterized by spontaneous continuous or paroxysmal pain and stimuli-evoked pain [1, 3] and has emerged as one of the most challenging pain syndromes, affecting 7–10% of the population [1]. Moreover, the incidence of neuropathic pain is increasing with the aging of the global population [1]. Despite recent improvements, the management of neuropathic pain remains unsatisfactory, mainly due to an insufficient understanding of the molecular mechanisms underlying the syndrome [2, 4]. Therefore, there is a need to reveal the underlying mechanisms to benefit the clinical management of neuropathic pain.

The dorsal root ganglion (DRG) contains a large proportion of the body’s sensory neurons, which are critical for transducing sensory information from the periphery to the central nervous system [5, 6]. Accumulating evidence has highlighted that the development of neuropathic pain is associated with changes in DRG neurons, such as the release of cytokines and dysregulation of ion channels [5]. In addition to changes in the expression of individual genes, recent studies have highlighted gene expression alterations in DRGs associated with neuropathic pain. After examining gene expression in the DRG of normal rats and peripheral axotomy-operated rats, Xiao et al. [7] identified significant differences in the expression of 173 genes, including neuropeptides, ion channels, and synaptic vesicle proteins. Reinhold et al. [8] used fluorescent neuronal tracers to confirm the differential expression patterns between damaged DRG neurons and adjacent neurons [7]. Furthermore, universal alterations of microRNA (miRNA) expression in DRG have been demonstrated in neuropathic pain. Li et al. [9] identified a total of 114 differentially expressed miRNAs in the DRG of a neuropathic pain rat model [8]. In another study, the expression of 42 miRNAs was significantly altered in the DRG of a streptozotocin-induced diabetic neuropathic pain mouse model [9].

Circular RNAs (circRNAs) are a large family of covalently closed RNA molecules with regulatory roles in gene transcription, miRNA function, and protein production [10,11,12]. They are expressed in diverse cell types and preferentially expressed in neural tissues [12]. The expression levels of circRNAs are dynamically modulated in neurons where are essential for synaptic plasticity and neuronal function [13]. In a recent study, 469 circRNAs were found to be differentially expressed in the spinal dorsal horn of chronic constriction injury (CCI) rats [14]. Wang et al. found that circHIPK3, which is involved in neuroinflammation, was associated with NP grade in type 2 diabetic patients. Silencing circHIPK3 can alleviate pain in diabetic NP rats [15].

Zhang et al. found that cytoplasmic circAnks1a can act as an endogenous miR-324-3p sponge, resulting in increased VEGFB expression, which was positively correlated with pain behavior induced by nerve injury [16]. In a CCI model, Zhang et al. confirmed that circ_0005075 can target miR-151a-3p, contributing to the pain behavior of rats by inducing NOTCH2 expression [17]. In another study, Cai et al. observed that ciRS-7, sponged to miR-135a-5p, contributed to the development of NP by upregulating autophagy and inflammation [18]. All of these study results indicated that alterations in the expression of circRNAs may play roles in neuropathic pain [14]. However, knowledge about the altered expression of circRNAs in the DRG of rats with neuropathic pain is limited.

To gain a global view of aberrant circRNA expression in the DRG of CCI rats, we established the CCI model and evaluated circRNA expression levels in the DRG using next-generation sequencing technology. Differentially expressed circRNAs (DEcircRNAs) between CCI and NC rats were identified and validated by quantitative PCR (qPCR). Functional annotation was implemented to explore the biological function of DEcircRNAs.

Materials and methods

Animals

Adult Sprague Dawley rats (male, 200–250 g) were acquired from the Experimental Animal Center of Fudan University (Shanghai, China). Age-matched rats were blinded and randomly assigned to CCI or negative control (NC) groups (n = 7/group). Rats were housed in a pathogen-free condition with a standard 12-h light/dark cycle and free access to rodent food and water. The animal experimental protocols (No. 2019-079) were approved by the Institutional Animal Care and Use Committee of Fudan University and complied with the ethical guidelines of the International Association for the Study of Pain to ensure minimal animal use and discomfort. All rats were acclimatized to laboratory conditions before the tests. All behavioral experiments were carried out between 9:00 a.m. and 12:00 p.m.

Experimental design

Our research goal was to examine the changes in circRNA expression in the DRG after surgery. All rats were randomized to the NC group and CCI surgery group, and mechanical allodynia and thermal hyperalgesia were assessed after surgery. All animals were euthanized by an overdose of sodium pentobarbital in 14 days after surgery. L4–L5 DRGs of rats ipsilateral to the nerve surgery were harvested for subsequent total RNA extraction.

Rat model of neuropathic pain

The CCI rat model was produced under sodium pentobarbital anesthesia (40 mg/kg, intraperitoneal injection) based on a procedure established by Bennett and Xie [19]. The sciatic nerves on the left side were exposed, and then ligatures were tied loosely around the nerve with approximately 1 mm space between knots. Rats in the NC group were subjected to sciatic nerve exposure without ligation.

Behavioral test

Paw withdrawal mechanical threshold (PWMT) test

According to previous studies [2, 20, 21], PWMT in response to mechanical stimuli was carried out with an electronic von Frey paw pressure anesthesiometer (IITC Life Science, Woodland Hills, CA, USA), which consists of a hand-held force transducer with a fixed 200 μm diameter tungsten wire tip. The experimenter applied the electronic von Frey probe tip perpendicularly to the plantar surface of the hind paw as described previously. Brisk withdrawal or flinching evoked by applications was considered a positive response. The test was performed 1 day before and 1, 3, 5, 7, 9, and 14 days after surgery. The three measurements of the mechanical threshold were averaged.

Paw withdrawal thermal latency (PWTL) tests

According to the protocol described previously [2, 20, 21], PWTL in response to radiant heat was carried out using a heat-mediated pain stimulator (Model 336; IITC Life Science). The time from onset of application to hind paw withdrawal was measured, and each rat was tested three times in 10-min intervals. The three measurements of thermal latency were averaged.

The test was also performed 1 day before and 1, 3, 5, 7, 9, and 14 days after surgery.

RNA isolation and quality control

The rats were anesthetized with sodium pentobarbital (40 mg/kg, intraperitoneal injection) and sacrificed by cervical dislocation at the end of the behavioral test. L4–L5 DRGs were excised from a single rat in each group (CCI and NC) and rapidly frozen at − 80 °C. Total RNA of each DRG sample was isolated using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s instructions. The quality and quantity of total RNA samples were measured with an ND-1000 NanoDrop spectrophotometer (Nano-Drop Technologies, Wilmington, DE, USA). The integrity of each RNA sample was assessed by denaturing agarose gel electrophoresis.

Next-generation RNA sequencing analysis

Total RNA from four matched DRG samples of CCI and NC rats was treated with the Epicenter Ribo-Zero rRNA Removal Kit (Illumina, San Diego, CA, USA) and RNAse R (Epicenter, Madison, WI, USA) to remove ribosomal and linear RNA, respectively. A strand-specific sequencing library was prepared by using the NEBNext Ultra Directional RNA Library Prep Kit (Illumina). Then, RNA sequencing (RNA-seq) was performed on the Illumina HiSeq 2500 system following the manufacturer’s recommendations. All sequencing procedures and analyses were performed by RiboBio (Guangzhou, China). High-throughput sequencing (HTSeq) software was used to count the numbers of sequence reads mapped to each circRNA [22]. TopHat mapping and TopHat-Fusion mapping reads were used by RPM (Reads Per Million mapped reads) to calculate the expression of individual circRNA [23, 24]. The total number of reads spanning back-spliced junctions was used as an absolute measure of circRNA abundance. ANNOVAR was used to analyse the functional elements from which the circRNAs were derived [25]. DEcircRNAs between the CCI and NC groups were identified using fold change (FC) filtering and Q value. The Q value was acquired by adjusting the p value using the Benjamini‒Hochberg method. We defined the statistical criteria for selecting aberrantly expressed circRNAs using a q-value of < 0.001 with an FC of ≥ 2 or ≤ 0.5. The cluster profile of DEcircRNAs was revealed by two-way hierarchical clustering analysis based on their expression values.

Functional annotation

The possible functions of circRNAs were assumed to be related to their host genes. Kyoto Encyclopedia of Genes and Genomes (KEGG) [26,27,28] and Gene Ontology (GO) [29] analyses were performed for DEcircRNAs to reveal neuropathic pain-associated biological pathways and functions. Fisher’s exact test was used to compute the statistical significance of the enrichment. The selection criterion for significantly enriched terms was set at FDR < 0.05.

qPCR

cDNA samples were prepared by reverse transcribing total RNAs using the PrimeScript RT Reagent Kit (TaKaRa, Tokyo, Japan). DEcircRNAs that were expressed in all samples and showed significant expression differences were selected for qPCR analysis. The expression levels of the selected DEcircRNAs were validated by qPCR assays using SYBR Green qPCR Master Mix (ABI, Los. Angeles, CA, USA) and specific primer sequences (Additional file 1: Table S1). The expression levels of circRNAs were normalized to GAPDH expression and analysed by the 2−ΔΔCt method.

Statistical analyses

Values were expressed as group mean and standard errors (mean ± SEM) and analysed using GraphPad Prism 8 (GraphPad Software, San Diego, CA, USA), Microsoft Excel (Microsoft, Washington, DC, USA), and R software version 3.5.3 (http://www.r-project.org/). Q-value of < 0.001 with an FC of ≥ 2 or ≤ 0.5 was used to identify significant differences in circRNA expression. According to previous studies [30, 31], two-way ANOVA followed by Bonferroni’s post-hoc tests was used to analyse the behavioral data. p < 0.05 was considered statistically significant, *p < 0.05, **p < 0.01, compared with NC groups.

Results

Development of mechanical and thermal hypersensitivity after CCI

To reveal circRNA expression changes induced by neuropathic pain, rat CCI models were established. Mechanical allodynia and thermal hyperalgesia were used to evaluate the models. As revealed by PWT and PWL tests, PWT and PWL were significantly decreased in the CCI group from postoperative day (POD) 3 to POD14 compared with the NC group (Fig. 1).

Fig. 1
figure 1

Chronic constriction injury-induced nociceptive behavior of rats. Mechanical allodynia (a) and thermal hyperalgesia (b) at each time point after CCI surgery. n = 7/group. All data are represented as the mean ± SEM, two-way ANOVA tests followed by Bonferroni’s post-hoc tests, **p < 0.01, CCI versus negative control (NC)

Full size image

Overview of circRNA profiles

RNA-seq was performed to profile circRNA expression in four paired DRG samples. A total of 13,172 circRNAs were detected in the L4–L5 DRG samples from the CCI and NC groups (Additional file 1: Table S2). Among these, 7514 circRNAs contained at least three unique back-spliced reads. According to their host gene location, the analysis of chromosome distribution showed that these 7514 circRNAs were derived from all chromosomes, including the mitochondrial genome. Chromosomes 1 to 3 produced more than 3000 circRNAs, and several hundred circRNAs were distributed on other chromosomes. Only 10 circRNAs were detected on chromosome Y, and 4 circRNAs were detected on the mitochondrial genome (Fig. 2a). According to their location on the genome, the host genes of these circRNAs were derived from exonic regions (62.0%), intergenic regions (24.0%), and intronic regions (6.6%) (Fig. 2b).

Fig. 2
figure 2

Chromosome distribution of circular RNAs. a Chromosome distribution of circular RNAs (circRNAs). b Genome distribution of circRNA sequence fragments

Full size image

CCI induces genome-wide circRNA differential expression in DRGs

Following differential expression analysis, 374 DEcircRNAs were identified between the CCI and NC groups by using FC filtering (FC ≥ 2 or ≤ 0.5) and FDR (Q < 0.001). Among the DEcircRNAs, 84 were downregulated and 290 were upregulated in the CCI group. A scatter plot was used to show variation in circRNA expression between these two groups of DRG samples (Fig. 3a). Moreover, a volcano plot was drawn to visualize the DEcircRNAs with statistical significance (Fig. 3b). Hierarchical clustering analysis was performed to group samples based on the expression levels of DEcircRNAs across samples. As shown by the heatmap, the CCI and NC samples showed distinguishable circRNA expression profiles. Venn diagram was made according to average of RPM in two groups and indicated the number of overlapping and distinct circRNAs between NC and CCI rats. (Fig. 3c, d, Additional file 2: Table S2 and Additional file 3: Table S3). The top 50 upregulated circRNAs are listed in Table 1, and the downregulated circRNAs are listed in Table 2.

Fig. 3
figure 3

Differentially expressed circRNAs in dorsal root ganglions between the negative control and chronic constriction injury model. a Scatter plots showing circular RNA (circRNA) expression differences in the L4–L5 dorsal root ganglions (DRGs) of chronic constriction injury (CCI) rats with respect to the negative control (NC) group. n = 4/group. The circRNAs above the top red line or below the bottom red line indicated differential circRNAs (fold change ≥ 2 or ≤ 0.5, and FDR < 0.001) respectively. b Volcano plot indicating differentially expressed circRNAs (DEcircRNAs) in the L4–L5 DRGs of CCI models. Red dots represent circRNAs with significant expression, whereas black dots represent circRNAs with no significant difference. c Heatmap of DEcircRNAs shows hierarchical clustering of DEcircRNAs of rats in the CCI group compared to the NC group. In the clustering analysis, upregulated and downregulated genes are shown in red and green, respectively. d The Venn diagram indicates the number of overlapping and distinct circRNAs between NC and CCI rats

Full size image
Table 1 Top 50 upregulated circRNAs between NC and CCI rats
Full size table
Table 2 Top 50 upregulated circRNAs between NC and CCI rats
Full size table

Functional annotation of DEcircRNAs

To predict the functions and pathways related to the DEcircRNAs identified above, KEGG and GO analyses were performed for the host genes of DEcircRNAs. According to our results, the most significantly enriched KEGG pathways included ‘dopaminergic synapse,’ ‘renin secretion,’ ‘long-term depression,’ ‘metabolic pathways,’ ‘mitogen-activated protein kinase (MAPK) signaling pathway,’ and ‘cyclic guanosine monophosphate-protein kinase G signaling pathway.’ Next, we performed GO analysis of the linear counterparts that produced DEcircRNAs. Our data showed that the most significantly enriched GO terms included ‘generation of neurons,’ ‘neurogenesis,’ ‘GTPase regulator activity,’ ‘negative regulation of metabolic process,’ and ‘cellular process.’ The top 30 significantly enriched KEGG and GO terms are shown in Fig. 4a, b. The results suggest that these signaling pathway and molecular function-related to DEcircRNAs may be involved in the pathogenesis and development of chronic pain.

Fig. 4
figure 4

Differentially expressed circRNA functional analysis between NC and CCI rats. a Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enriched by DEcircRNAs. The vertical axis is the pathway category. The horizontal axis is the −log10 (p value) of the pathway and the number of differentially expressed circRNAs (DEcircRNAs). Different colors are used to distinguish cellular processes, environmental information processing, human diseases, metabolism and organismal systems. b The most enriched Gene Ontology (GO) terms of DEcircRNAs. The vertical coordinate is the enriched GO terms, and the horizontal coordinate is the rich factor. Rich factor value represents the enrichment degree. The size of the points represents the number of DEcircRNAs in the term. The colors of the points correspond to different p value ranges. c CircRNA–microRNA competing endogenous RNA network. The red nodes represent the upregulated circRNAs, and the green nodes represent the downregulated circRNAs. The size of each node represents the functional connectivity of each circRNA. The line color indicates potential interactions between circRNAs and microRNAs, which were predicted by using TargetScan and miRanda. d CircRNA expression changes validated by quantitative PCR (qPCR). qPCR analysis of differences in the expression levels of DEcircRNAs in the dorsal root ganglion between the CCI model and NC groups. The results were calculated by normalizing to GAPDH in the same sample with the 2−ΔΔCt method. Changes in relative levels of circRNAs expressed as fold of the controls. All values are expressed as the mean ± SEM, Student’s t test, *p < 0.05 (n = 3)

Full size image

CircRNA–miRNA competing endogenous RNA network

There is emerging evidence that miRNAs play an important role in the development of chronic pain. In our previous study, we found that miRNA-146a-5p, a key molecule in the negative regulation of innate immunity, has dysregulated expression in the DRG of rats with chronic pain [32]. In addition, the aberrant expression of other miRNAs in DRG has also been reported, such as miR-155-5p [33], miR-21 [34], and let-7a-5p [35]. Several circRNAs have been recently identified that exert important regulatory effects as miRNA sponges. Potential connections between circRNAs and microRNAs involved in chronic pain were predicted by using TargetScan and miRanda, and the interaction results were displayed by Cytoscape. Four chronic pain-related miRNAs exhibited an interaction with the top 50 circRNAs according to the size of p value (Fig. 4c and Additional file 4: Table S4).

Validation of circRNA expression using qPCR

To validate the screened circRNA candidates, qPCR analysis was performed to amplify the back-spliced junctions of nine typical DEcircRNAs (chr1:166208849–166259335, chr2:26590643–26650618, chr3:112752502–112762808, chr4:132075419–132113888, chr7:64125692–64131186, chr10:8331904–8371302, chr13:84161876–84168811, chr15:51936366–51944899, and chr20:50725400–50764570). As detected by circRNA HTS, these circRNAs were expressed in all CCI and NC samples and showed significant expression differences between groups. In addition, all of these circRNAs were predicted to have at least one binding site with miR-146a-5p. According to our qPCR results, with the exception of chr10:8331904–8371302, chr13:84161876–84168811, and chr20:50725400–50764570, the expression levels of all remaining DEcircRNAs were consistent with those measured by sequencing analysis (Fig. 4d).

Discussion

CircRNAs are mainly alternative splice variants of pre-mRNA [36]. Alterations in universal circRNA expression have been identified in various diseases, including cancer [37], Alzheimer’s disease [38], and type 2 diabetes [39]. In our study, we revealed circRNA expression changes at the genome-wide level by comparing circRNA expression in the DRGs of CCI and NC rats.

The CCI is one of the most frequently used models of neuropathic pain [40]. As shown by our pain behavior test, PWMT and PWTL were significantly decreased in the CCI group compared with the NC group, indicating that the CCI rat model was successfully established. CircRNA samples were prepared from the DRGs of CCI and NC rats and were subjected to circRNA sequencing. Differential expression analysis identified upregulated expression changes in 290 circRNAs and downregulated expression changes in 84 circRNAs of DRG neurons of CCI induced neuropathic pain rats compared with the control group. The number of downregulated DEcircRNAs was great less than that of the upregulated ones. We listed top 50 upregulated DEcircRNAs and top 50 downregulated DEcircRNAs respectively. Among these DEcircRNAs, circ:chr13:14728984–14734632 in CCI rats was upregulated more than tenfold and the circ:chr20:37502198–37560959 in DRG of CCI rats was downregulated more than 12 fold compared with the rats in control group. In order to confirm the reliability of our sequencing data, we picked nine circRNAs with high expression level and significant difference in the expression between CCI and NC groups. Expression patterns of these circRNAs was performed by qPCR and the consistent trends were compared between qPCR results and the circRNAs HTSeq.

As we known, circRNAs could expressed dynamically and exert their function through multiple mechanisms [41]. However, the vast majority of circRNAs are generated from the sequences of their host genes, and some circRNAs also play roles in transcriptional or posttranscriptional regulation of their host genes [36]. Therefore, the functions of circRNAs may be relevant to their host genes. In our study, we first conducted the GO analysis and KEGG pathway analysis for the host genes of the DEcircRNAs to predict the function of them. According to our results, a wide variety of KEGG and GO terms were identified, indicating a broad range of functional modulations in the DRG of neuropathic pain rats. The GO analysis indicated that the DEcircRNAs were presumably involved in development of neurons and nervous tissue, the regulation of GTPase, metabolic and cellular process.

The top five significant functional terms conducted by KEGG pathway annotation included ‘dopaminergic synapse,’ ‘renin secretion,’ ‘MAPK signaling pathway,’ and ‘neurogenesis.’ It could be reasonable, as many of these processes were considered contributing to the development and maintenance of chronic neuropathic pain. Dopaminergic functions are commonly modulated in neuropathic pain, and administration of dopaminergic agents may relieve neuropathic pain [42]. Renin is a component of the renin–angiotensin–aldosterone system, and increased levels of renin are associated with neuropathic pain [43]. Increased phosphorylation of MAPK in the DRG is an early event in neuropathic pain following peripheral nerve injury, and the extracellular signal-related kinase/MAPK pathway is considered a novel therapeutic target for neuropathic pain [44]. Dysregulation of neurogenesis in the hippocampus is also suggested to be involved in neuropathic pain [45]. Therefore, we hypothesized that the biological functions of DRG are widely modulated in neuropathic pain.

In our and others’ previous research, miRNAs were proved to play an important role in development of neuropathic pain. miR-146a and miR-155, taking a crucial part in inflammation, immunity as well as tumor, were highly expressed in DRG neurons of chronic neuropathic pain rats [32, 33]. miR-21 may serve as an endogenous ligand, acting on TLR8 in endosomes and lysosomes, and involve in the maintenance of neuropathic pain [34]. Recent evidence revealed that the essential function of circRNAs is serving as miRNAs sponges. We predicted the circRNAs/miRNA interaction using TargetScan and miRanda and then constructed ceRNA network. According to our results, we deduced that these circRNAs may bind pain related miRNAs and inhibit their biological functions. Our results implied that it is essential and worthwhile to further study potential functions of these pain related DEcircRNAs in chronic neuropathic pain. Gao et al. analyzed the circRNA expression in the dorsal horn of the spinal cord using microarrays. They only used fold changes ≥ 2 as selection criterion to select significant DEcircRNAs and found expression change of 469 circRNAs in CCI rats compared with sham-operated rats. Similar to our KEGG pathway annotation results, they also reported the MAPK signaling pathway was related to neuropathic pain [14]. However, we first identified circular RNA expression profiles of the dorsal root ganglion in a rat model of neuropathic pain. Our results may help to elucidate peripheral mechanism of neuropathic pain.

The main strength of this study was that circRNA expression in DRGs was compared between CCI and NC rats, which provided a global picture of circRNA expression alterations in DRG. The limitation of our study was that the CCI model was the only neuropathic model included in our study. Thus, additional neuropathic models are required in future studies to validate our results. Further functional studies in vivo should also be performed to elucidate the roles of DEcircRNAs in neuropathic pain.

Conclusions

We developed a CCI model and identified 374 significant DEcircRNAs in the DRGs of CCI and NC rats. These DEcircRNAs may be involved in various biological pathways and functions, which are potentially modulated in DRGs by neuropathic pain. Further bioinformatics analyses suggested that these DEcircRNAs could be sponged with miRNAs, which targeted mRNAs associated with chronic pain. Our study may provide new insights into the molecular mechanisms underlying neuropathic pain.

Availability of data materials

There are no data, software, databases, or applications/tools available apart from those reported in the present study. All data are provided in the manuscript and supplementary data section.

Abbreviations

cGMP:

Cyclic guanosine monophosphate

ERK:

Extracellular signal-regulated protein kinase

MAPK:

Mitogen-activated protein kinase

PKG:

Protein kinase G

References

  1. Colloca L, Ludman T, Bouhassira D, Baron R, Dickenson AH, Yarnitsky D, Freeman R, Truini A, Attal N, Finnerup NB, et al. Neuropathic pain. Nat Rev Dis Primers. 2017;3:17002.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Wang Z, Liu F, Wei M, Qiu Y, Ma C, Shen L, Huang Y. Chronic constriction injury-induced microRNA-146a-5p alleviates neuropathic pain through suppression of IRAK1/TRAF6 signaling pathway. J Neuroinflam. 2018;15(1):179.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Bouhassira D, Attal N. Translational neuropathic pain research: a clinical perspective. Neuroscience. 2016;338:27–35.

    Article  CAS  PubMed  Google Scholar 

  4. Baron R, Binder A, Wasner G. Neuropathic pain: diagnosis, pathophysiological mechanisms, and treatment. Lancet Neurol. 2010;9(8):807–19.

    Article  PubMed  Google Scholar 

  5. Krames ES. The role of the dorsal root ganglion in the development of neuropathic pain. Pain Med. 2014;15(10):1669–85.

    Article  PubMed  Google Scholar 

  6. Liem L, van Dongen E, Huygen FJ, Staats P, Kramer J. The dorsal root ganglion as a therapeutic target for chronic pain. Reg Anesth Pain Med. 2016;41(4):511–9.

    Article  PubMed  Google Scholar 

  7. Reinhold AK, Batti L, Bilbao D, Buness A, Rittner HL, Heppenstall PA. Differential transcriptional profiling of damaged and intact adjacent dorsal root ganglia neurons in neuropathic pain. PLoS ONE. 2015;10(4): e0123342.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Li S, Yu B, Wang S, Gu Y, Yao D, Wang Y, Qian T, Ding F, Gu X. Identification and functional analysis of novel micro-RNAs in rat dorsal root ganglia after sciatic nerve resection. J Neurosci Res. 2012;90(4):791–801.

    Article  CAS  PubMed  Google Scholar 

  9. Gong Q, Lu Z, Huang Q, Ruan L, Chen J, Liang Y, Wang H, Yue Y, Feng S. Altered microRNAs expression profiling in mice with diabetic neuropathic pain. Biochem Biophys Res Commun. 2015;456(2):615–20.

    Article  CAS  PubMed  Google Scholar 

  10. Salzman J. Circular RNA expression: its potential regulation and function. Trends Genet. 2016;32(5):309–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, Maier L, Mackowiak SD, Gregersen LH, Munschauer M, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495(7441):333–8.

    Article  CAS  PubMed  Google Scholar 

  12. Floris G, Zhang L, Follesa P, Sun T. Regulatory role of circular RNAs and neurological disorders. Mol Neurobiol. 2017;54(7):5156–65.

    Article  CAS  PubMed  Google Scholar 

  13. Hanan M, Soreq H, Kadener S. CircRNAs in the brain. RNA Biol. 2017;14(8):1028–34.

    Article  PubMed  Google Scholar 

  14. Cao S, Deng W, Li Y, Qin B, Zhang L, Yu S, Xie P, Xiao Z, Yu T. Chronic constriction injury of sciatic nerve changes circular RNA expression in rat spinal dorsal horn. J Pain Res. 2017;10:1687–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Wang L, Luo T, Bao Z, Li Y, Bu W. Intrathecal circHIPK3 shRNA alleviates neuropathic pain in diabetic rats. Biochem Biophys Res Commun. 2018;505(3):644–50.

    Article  CAS  PubMed  Google Scholar 

  16. Zhang SB, Lin SY, Liu M, Liu CC, Ding HH, Sun Y, Ma C, Guo RX, Lv YY, Wu SL, et al. CircAnks1a in the spinal cord regulates hypersensitivity in a rodent model of neuropathic pain. Nat Commun. 2019;10(1):4119.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Zhang Y, Gao T, Li X, Wen CC, Yan XT, Peng C, Xiao Y. Circ_0005075 targeting miR-151a-3p promotes neuropathic pain in CCI rats via inducing NOTCH2 expression. Gene. 2021;767: 145079.

    Article  CAS  PubMed  Google Scholar 

  18. Cai W, Zhang Y, Su Z. ciRS-7 targeting miR-135a-5p promotes neuropathic pain in CCI rats via inflammation and autophagy. Gene. 2020;736: 144386.

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  20. Liu F, Shen X, Su S, Cui H, Fang Y, Wang T, Zhang L, Huang Y, Ma C. Fcγ receptor I-coupled signaling in peripheral nociceptors mediates joint pain in a rat model of rheumatoid arthritis. Arthritis Rheumatol. 2020;72(10):1668–78.

    Article  CAS  PubMed  Google Scholar 

  21. Liu F, Wang Z, Qiu Y, Wei M, Li C, Xie Y, Shen L, Huang Y, Ma C. Suppression of MyD88-dependent signaling alleviates neuropathic pain induced by peripheral nerve injury in the rat. J Neuroinflammation. 2017;14(1):70.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9.

    Article  CAS  PubMed  Google Scholar 

  23. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Kim D, Salzberg SL. TopHat-Fusion: an algorithm for discovery of novel fusion transcripts. Genome Biol. 2011;12(8):R72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16): e164.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51.

    Article  CAS  PubMed  Google Scholar 

  27. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Tiwari V, Anderson M, Yang F, Tiwari V, Zheng Q, He S-Q, Zhang T, Shu B, Chen X, Grenald SA, et al. Peripherally acting μ-opioid receptor agonists attenuate ongoing pain-associated behavior and spontaneous neuronal activity after nerve injury in rats. Anesthesiology. 2018;128(6):1220–36.

    Article  CAS  PubMed  Google Scholar 

  31. Tiwari V, He S-Q, Huang Q, Liang L, Yang F, Chen Z, Tiwari V, Fujita W, Devi LA, Dong X, et al. Activation of µ-δ opioid receptor heteromers inhibits neuropathic pain behavior in rodents. Pain. 2020;161(4):842–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Wang Z, Liu F, Wei M, Qiu Y, Ma C, Shen L, Huang YJJ. Chronic constriction injury-induced microRNA-146a-5p alleviates neuropathic pain through suppression of IRAK1/TRAF6 signaling pathway. J Neuroinflamm. 2018;15(1):179.

    Article  Google Scholar 

  33. Dou L, Lin H, Wang K, Zhu G, Zou X, Chang E, Zhu Y. Long non-coding RNA CCAT1 modulates neuropathic pain progression through sponging miR-155. Oncotarget. 2017;8(52):89949–57.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Zhang ZJ, Guo JS, Li SS, Wu XB, Cao DL, Jiang BC, Jing PB, Bai XQ, Li CH, Wu ZH, et al. TLR8 and its endogenous ligand miR-21 contribute to neuropathic pain in murine DRG. J Exp Med. 2018;215(12):3019–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Motti D, Lerch JK, Danzi MC, Gans JH, Kuo F, Slepak TI, Bixby JL, Lemmon VP. Identification of miRNAs involved in DRG neurite outgrowth and their putative targets. FEBS Lett. 2017;591(14):2091–105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Barrett SP, Salzman J. Circular RNAs: analysis, expression and potential functions. Development. 2016;143(11):1838–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Patop IL, Kadener S. circRNAs in cancer. Curr Opin Genet Dev. 2018;48:121–7.

    Article  CAS  PubMed  Google Scholar 

  38. Lukiw WJ. Circular RNA (circRNA) in Alzheimer’s disease (AD). Front Genet. 2013;4:307.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Zhao Z, Li X, Jian D, Hao P, Rao L, Li M. Hsa_circ_0054633 in peripheral blood can be used as a diagnostic biomarker of pre-diabetes and type 2 diabetes mellitus. Acta Diabetol. 2017;54(3):237–45.

    Article  CAS  PubMed  Google Scholar 

  40. Jaggi AS, Jain V, Singh N. Animal models of neuropathic pain. Fundam Clin Pharmacol. 2011;25(1):1–28.

    Article  CAS  PubMed  Google Scholar 

  41. Weng Y, Wu J, Li L, Shao J, Li Z, Deng M, Zou W. Circular RNA expression profile in the spinal cord of morphine tolerated rats and screen of putative key circRNAs. Mol Brain. 2019;12(1):79.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Hagelberg N, Jaaskelainen SK, Martikainen IK, Mansikka H, Forssell H, Scheinin H, Hietala J, Pertovaara A. Striatal dopamine D2 receptors in modulation of pain in humans: a review. Eur J Pharmacol. 2004;500(1–3):187–92.

    Article  CAS  PubMed  Google Scholar 

  43. Muthuraman A, Kaur P. Renin–angiotensin–aldosterone system: a current drug target for the management of neuropathic pain. Curr Drug Targets. 2016;17(2):178–95.

    Article  CAS  PubMed  Google Scholar 

  44. Ma W, Quirion R. The ERK/MAPK pathway, as a target for the treatment of neuropathic pain. Expert Opin Ther Targets. 2005;9(4):699–713.

    Article  CAS  PubMed  Google Scholar 

  45. Grilli M. Chronic pain and adult hippocampal neurogenesis: translational implications from preclinical studies. J Pain Res. 2017;10:2281–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by grants from the National Natural Science Foundation of China (Grant Nos. 81801114, 81801376, 81901125, and 81901130] and the China Postdoctoral Science Foundation (Grant No. 2018M631389). The funding bodies did not have any influence on the design of the study; collection, analysis, and interpretation of the data; or writing of the manuscript.

Author information

Author notes

  1. Wanxia Xiong, Min Wei, and Li Zhang contributed equally to this work

Authors and Affiliations

  1. Department of Anesthesiology, Zhongshan Hospital, Fudan University, Shanghai, China

    Wanxia Xiong, Jie Wang & Zhiyao Wang

  2. Center for Anesthesiology, Beijing Anzhen Hospital, Capital Medical University, Beijing, China

    Min Wei

  3. Department of Anesthesiology, Beijing Friendship Hospital, Capital Medical University, Beijing, China

    Li Zhang

  4. Department of Human Anatomy, Histology and Embryology, Institute of Basic Medical Sciences Chinese Academy of Medical Sciences, School of Basic Medicine Peking Union Medical College, Beijing, China

    Fan Liu

Authors
  1. Wanxia Xiong
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Min Wei
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Li Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Jie Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Fan Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Zhiyao Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

FL and ZW conceived of the study and participated in its design and coordination to draft the manuscript. ZW, MW and LZ performed the CCI and behavioral assessment of pain. JW, MW, LZ, ZW and WX acquired, analyzed and interpreted the data. All authors read and approved the final manuscript. All authors reviewed the manuscript.

Corresponding authors

Correspondence to Fan Liu or Zhiyao Wang.

Ethics declarations

Ethics approval and consent to participate

All animal procedures performed in this study were reviewed and approved by the Institutional Animal Care and Use Committee of Fudan University and were conducted in accordance with the guidelines of the International Association for the Study of Pain. The study is reported in accordance with ARRIVE guidelines.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

List of DNA primer sequences designed for qPCR.

Additional file 2: Table S2.

List of total circRNAs detected in the DRG of NC and CCI rats. RNA-seq was performed to profile circRNA expression in four paired DRG samples. A total of 13,172 circRNAs were detected in the L4–L5 DRG samples from the CCI and NC groups (4 mice per group).

Additional file 3: Table S3.

List of DEcircRNAs were identified between the CCI and NC groups. Following differential expression analysis, 374 DEcircRNAs were identified between the CCI and NC groups (4 mice per group) by using FC filtering (FC ≥ 2 or ≤ 0.5) and FDR (Q <  0.001).

Additional file 4: Table S4.

List of circRNAs and miRNA predicted by using miRanda and TargetScan. Potential connections between circRNAs and microRNAs were predicted by using TargetScan and miRanda, and the interaction results were displayed by Cytoscape. All the circRNAs and miRNAs were listed.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xiong, W., Wei, M., Zhang, L. et al. Chronic constriction injury-induced changes in circular RNA expression profiling of the dorsal root ganglion in a rat model of neuropathic pain. BMC Neurosci 23, 64 (2022). https://doi.org/10.1186/s12868-022-00745-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12868-022-00745-5

Keywords

BMC Neuroscience

ISSN: 1471-2202

Contact us