Back to Journals » Journal of Pain Research » Volume 10

Chronic constriction injury of sciatic nerve changes circular RNA expression in rat spinal dorsal horn

Authors Cao S, Deng W, Li Y, Qin B, Zhang L, Yu S, Xie P, Xiao Z, Yu T

Received 13 April 2017

Accepted for publication 17 June 2017

Published 17 July 2017 Volume 2017:10 Pages 1687—1696

DOI https://doi.org/10.2147/JPR.S139592

Checked for plagiarism Yes

Review by Single anonymous peer review

Peer reviewer comments 3

Editor who approved publication: Dr Katherine Hanlon



Song Cao,1,2,* Wenwen Deng,3,* Ying Li,1 Bangyong Qin,1 Lin Zhang,2 Shouyang Yu,2 Peng Xie,2 Zhi Xiao,4 Tian Yu2

1Department of Pain Medicine, Affiliated Hospital of Zunyi Medical University, 2Guizhou Key Laboratory of Anesthesia and Organ Protection, Zunyi Medical University, 3Department of Cardiology, Affiliated Hospital of Zunyi Medical University, 4Research Center for Medicine and Biology, Zunyi Medical University, Zunyi, Guizhou, China

*These authors contributed equally to this work

Background: Mechanisms of neuropathic pain are still largely unknown. Molecular changes in spinal dorsal horn may contribute to the initiation and development of neuropathic pain. Circular RNAs (circRNAs) have been identified as microRNA sponges and involved in various biological processes, but whether their expression profile changes in neuropathic pain condition is not reported.
Methods: To test whether neuropathic pain influences circRNA expression, we developed a sciatic chronic constriction injury (CCI) model in rats. The CCI ipsilateral spinal dorsal horns of lumbar enlargement segments (L3–L5) were collected, and the total RNA was extracted and subjected to Arraystar Rat circRNA Microarray. Quantitative real-time polymerase chain reaction (qPCR) was used to confirm the circRNA expression profile. To estimate functions of differential circRNAs, bioinformatics analyses including gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes Pathway analyses were performed for the top 100 circRNAs and circRNA–microRNA networks were constructed for the top 10 circRNAs.
Results: circRNA microarrays showed that 469 circRNAs were differentially expressed between CCI and sham-operated rats (fold change ≥2). In all, 363 of them were significantly upregulated, and the other 106 were downregulated in the CCI group. Three of them (circRNA_013779, circRNA_008008, and circRNA_003724) overexpressed >10 times after CCI insult. Expression levels of eight circRNAs were verified using qPCR. GO analysis revealed that thousands of predicted target genes were involved in the biological processes, cellular component, and molecular function; in addition, dozens of these genes were enriched in the Hippo signaling pathway, MAPK signaling pathway, and so on. Competing endogenous RNAs analysis showed that circRNA_008008 and circRNA_013779 are the two largest nodes in the circRNA–microRNA interaction network of the top 10 circRNAs.
Conclusion: CCI resulted in a comprehensive expression profile of circRNAs in the spinal dorsal horn in rats. CircRNAs in the dorsal horn could be helpful to reveal molecular mechanisms of neuropathic pain.

Keywords: neuropathic pain, circular RNA, circRNA–microRNA interaction, microarray

Introduction

Millions of people are suffering from neuropathic pain worldwide, which is still a major health and economic burden.1,2 The exact molecular mechanisms of neuropathic pain are largely unknown and the elucidation of its molecular mechanisms will be beneficial for the development of mechanism-oriented treatment.3

The spinal dorsal horn harbors secondary neurons that receive nociceptive stimuli. It is reported that the dorsal horn plays an important role in the generation of neuropathic pain and the occurrence of central sensitization of pain.47 Dorsal horn is also the location where inhibition of the transmission of nociceptive stimuli takes place because the descending inhibitory neuronal terminals end here;8,9 therefore, it is presumed that regulatory molecules play a role in the dorsal horn of the spinal cord in neuropathic pain, because neuroplastic changes in the peripheral (such as the dorsal root ganglion) and central nervous system (such as the spinal dorsal horn and brain) of the pain processing network take part in the development of pathological pain.10,11

MicroRNAs are noncoding posttranscriptional molecules and important biomarkers, which modulate numerous pathophysiological processes including neuropathic pain.12 For example, microRNA-7a,13 microRNA-1,14,15 and microRNA-30b16 in the dorsal root ganglion have been shown to be involved in the induction of neuropathic pain. Although one study reported that chronic constriction injury (CCI) did not lead to relevant changes in spinal miRNA expression levels in rats, Genda et al17 detected >100 differentially expressed microRNAs in the dorsal horn of CCI rats. Another study reported that microRNA-203 in the dorsal horn decreased as many as 10 fold in CCI rats.18

Circular RNAs (circRNAs) are another class of noncoding RNAs which interact with microRNAs.19 CircRNAs act as miRNA sponges2022 and regulate gene expression through a circRNA–microRNA–mRNA pathway.23,24 Accordingly, we believed that altered expression of circRNAs may accompany the development and progress of neuropathic pain. To date, studies have been performed probing the promising role of circRNAs as biomarkers of diseases such as Alzheimer’s disease,25 coronary artery disease,26 tumor,27,28 and other diseases. However, whether circRNAs show different expression profile in neuropathic pain states has not been reported.

To explore the expression of circRNAs in neuropathic pain condition, we developed the widely used CCI neuropathic pain model and detected the differentially expressed circRNAs between sham-operated and neuropathic pain rats with circRNA microarrays. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were carried out according to the microarray results and bioinformatics predictions and the circRNA–microRNA–mRNA interaction mode.

Methods

Animals

The experimental protocols were approved by the Experimental Animal Care and Use Committee of Zunyi Medical University. All experiments were conducted in accordance with the ethical guidelines of the International Association for the Study of Pain. All efforts were made to minimize animal suffering and to reduce animal use. Sprague Dawley rats (male, 7–8 weeks, 250–300 g) were housed under approved conditions with 12/12 hour light/dark cycles. Rats were free to eat and to drink water.

Induction of neuropathic pain

Rats underwent CCI surgeries to induce neuropathic pain on the left hind limb. Animals were randomly allocated to sham or CCI group. Neuropathic pain was induced with CCI as previously described by Bennett and Xie29 with a 5-0 chromic gut suture.30 In brief, after the intraperitoneal injection of sodium pentobarbital (35 mg/kg), the left sciatic nerve was loosely ligated with four sutures distant by 1–1.5 mm at the upstream of the nerve trifurcation of tibial, sural, and the common peroneal nerve. The muscle and skin were closed with sutures. Rats of the sham group underwent the same anesthesia and surgical procedures but the sciatic nerves were not ligated. Besides baseline tests, mechanical pain intensity was monitored at 2, 6, 10, and 14 days post-surgery.

Nociceptive behavioral testing

Mechanical hypersensitivity was determined by a male experimenter using an electronic von Frey plantar aesthesiometer (IITC, Wood Dale, IL, USA). Before tests, rats were given 10 minutes to habituate to the test environment. A rigid tip was applied against the mid-plantar surface of the left hind paw. The paw withdrawal threshold was automatically recorded by the device and the cutoff was set at 50 g. The rigid tip was presented perpendicularly to the plantar surface, and brisk withdrawal or paw flinching were considered as positive responses; the digital number presented on the monitor was recorded as the paw mechanical withdrawal threshold (MWT). Three successive stimuli were applied. MWT for individual animals was represented by the mean values (calculated from the three successive stimuli).

Ipsilateral dorsal horn samples preparation

After MWT tests on the 14th day post-surgery, rats were deeply anesthetized with isoflurane and were rapidly decapitated. The lumbar enlargement segments (L3–L5) of spinal cords were transversely sectioned and hemi-dissected along the midline. Only the dorsal half of the lumbar enlargement at the operated side of both sham-operated rats and CCI rats was collected and stored at −80°C.

RNA isolation, purification, and hybridization

One sample for each group (sham and CCI) was collected by pooling six ipsilateral L3–L5 dorsal spinal cords. RNA preparation and microarray hybridization were performed based on Arraystar’s standard protocols (Arraystar Inc, Rockville, MD, USA).31 The RNA amount of each sample was quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA) and treated with Rnase R (Epicentre Inc, Madison, WI, USA) to remove linear RNAs and to enrich circRNAs. Then, the enriched circRNAs were amplified and transcribed into fluorescent cRNA using a random priming method (Arraystar Super RNA Labeling Kit; Arraystar Inc). The labeled cRNAs were hybridized onto the Arraystar Rat circRNA Array (8×15K, Arraystar Inc). At last, the arrays were scanned by an Agilent Scanner G2505C (Agilent Technologies, Santa Clara, CA, USA).

Microarray data analysis

Agilent Feature Extraction software (version 11.0; Agilent Technologies) was used to analyze acquired array images. Quantile normalization and subsequent data processing were performed using the R software Limma package (http://www.bioconductor.org/packages/release/bioc/html/limma.html). Differentially expressed circRNAs between two samples were identified through fold change filtering. CircRNAs having fold changes ≥2 were selected as of significantly differential expression. Hierarchical clustering was performed to show the distinguishable circRNAs expression pattern between sham and CCI samples.

Real-time polymerase chain reaction (PCR) validation

Quantitative real-time PCR (qPCR) was used to confirm the circRNA expression profiles obtained from the microarray data. Total RNA aforementioned was reversely transcribed into cDNA using a circRNA qRT-PCR kit (GeneSeed, Guangzhou, China) with TransScript RT enzyme (M-MLV) and reverse transcription primer according to the manufacturer’s protocol. The relative gene expression was determined using a CFX Connect Real-Time system (Bio-Rad Laboratories Inc, Hercules, CA, USA). All samples were normalized to the signal generated from GAPDH. Data are shown as fold change (2−ΔΔCt). Each sample was tested in triplicate. Primer sequences are listed in Table S1.

MicroRNA prediction

The circRNA–microRNA interaction was predicted with Arraystar’s home-made microRNA target prediction software based on the TargetScan32 and miRanda33 prediction algorithm.

MicroRNA target genes prediction and bioinformatics analysis

To further investigate the functional roles of microRNA, putative targets of miRNAs were predicted by the TargetScan software. GO analysis was performed to explore the functional roles of target genes in terms of biological processes, cellular components, and molecular functions. Biological pathways defined by KEGG, Biocarta, and Reactome (http://www.genome.jp/kegg/) were identified by Database for Annotation, Visualization and Integrated Discovery (https://david.ncifcrf.gov/).

CircRNA–microRNA interaction network

To further elucidate correlations between circRNAs and microRNA, potential circRNA–microRNA interaction analyses were conducted and maps were drawn with Cytoscape. The size of each node represents the number of putative microRNA functionally connected to each circRNA.

Statistical analysis

Results are reported as mean ± standard deviation. Statistically significant differences between the sham and CCI group were estimated by the repeated measures two-way analysis of variance using GraphPad Prism 7.0 (GraphPad Software Inc, La Jolla, CA, USA). P<0.05 was considered as being statistically significant.

Results

Development of mechanical hypersensitivity following CCI surgery

Compared with the sham group, rats in the CCI group showed significantly higher mechanical hypersensitivity on the ipsilateral hind paw at all of the four tested time points after CCI surgery (all P<0.001; Figure 1).

Figure 1 CCI-induced mechanical hypersensitivity.

Notes: Changes in MWT were assessed with the electronic von Frey filament on ipsilateral hind paws of CCI and sham-operated rats. n=6 in both groups. All data are represented as mean ± SD. Statistical analyses consisted of repeated measures two-way ANOVA tests. ***P<0.001 between sham and CCI groups.

Abbreviations: CCI, chronic constriction injury; SD, standard deviation; ANOVA, analysis of variance; MWT, mechanical withdrawal thresholds.

Overview of circRNA profiles

The Arraystar Rat circRNA Microarray detected 12,770 rat circRNAs (Supplementary materials). Hierarchical clustering and scatter plot visualization showed the circRNAs expression levels were different (Figure 2). Overall, 469 circRNAs were listed as differential circRNAs (fold change ≥2) between the CCI and sham group. Up to 363 circRNAs were significantly upregulated and 106 were downregulated in the CCI dorsal horn compared with that of sham-operated rats (Supplementary materials). The top 20 upregulated and downregulated circRNAs are listed in Table 1.

Figure 2 CircRNA expression between the CCI and sham-operated rats.

Notes: (A) Hierarchical clustering shows differential circRNA expression profiles between two groups. The upregulated circRNAs in the CCI group are shown in red, and the downregulated circRNAs are indicated in green. (B) Box plots show the distribution of circRNAs for the two compared samples. The distributions were nearly the same after normalization. (C) Scatter plots assess the circRNA expression differences between the two compared groups. The circRNAs above the top green line and below the bottom green line indicated differential circRNAs (fold change ≥2 or ≤−2, respectively).

Abbreviations: CCI, chronic constriction injury; circRNA, circular RNA.

Table 1 Top 20 upregulated and downregulated circRNAs between sham-operated and CCI rats (sham vs CCI)

Abbreviations: circRNAs, circular RNAs; CCI, chronic constriction injury.

qPCR validation

To validate the microarray data, qPCR was employed to detect circRNA expression. A total of eight circRNAs (four upregulated and four downregulated in the CCI group) were selected. Expression levels detected by the microarray and qPCR are presented in Figure 3. The expression trends detected by the two methods were consistent with each other, which demonstrated the high reliability of the microarray results.

Figure 3 Validation of eight differentially expressed circRNAs by qPCR.

Notes: qPCR results were calculated with 2−ΔΔCt method and normalized to GAPDH expression in the sham and CCI group. Bars represent the mean ± SD; positive fold change means overexpression in the CCI group. Each qPCR assay was performed three times.

Abbreviations: CCI, chronic constriction injury; SD, standard deviation; circRNA, circular RNA; qPCR, quantitative real-time polymerase chain reaction.

MicroRNA prediction and bioinformatics analyses

MicroRNA prediction was carried out with Arraystar’s home-made miRNA target prediction software based on miRanda and TargetScan. Target genes of the 505 predicted miRNAs of the top 100 differential circRNAs were predicted. We used miRDB (version 5) to forecast miRNAs’ target genes. Results of the GO enrichment analysis of the top 100 circRNAs with identified target genes are shown in Figure 4. GO analysis revealed thousands of target genes involved in the biological processes, such as cellular component organization (GO 0048856) etc, the cellular component, such as nucleoplasm (GO 0012505) etc, and molecular function, such as binding (GO 0005488) etc. KEGG analysis showed that there were 112 pathways related to the 100 circRNAs, including Hippo signaling pathway, MAPK signaling pathway, endocytosis and longevity regulating pathway etc (Supplementary materials). The top 20 high enrichment score pathways are listed in Figure 5.

Figure 4 GO enrichment analysis. GO enrichment corresponds to the top 100 circRNAs.

Abbreviations: circRNAs, circular RNAs; GO, gene ontology.

Figure 5 Top 20 KEGG pathways with the highest enrichment score.

Notes: Pathways correspond to the top 100 circRNAs. The vertical axis is the pathway category and the horizontal axis is the –LgP of the pathway; LgP is the logarithm of the P-value, and P<0.05 was considered significant.

Abbreviations: circRNAs, circular RNAs; KEGG, Kyoto Encyclopedia of Genes and Genomes.

CircRNA–microRNA competing endogenous RNAs (ceRNA) network

Potential connections between the top 10 circRNAs and microRNAs were predicted by using TargetScan and miRanda and the interactions were displayed with Cytoscape. As shown in Figure 6, seven of the top 10 circRNAs exhibited interactions with microRNAs. CircRNA_001059 and circRNA_000167 were the two largest nodes in the network (with eight predicted microRNAs).

Figure 6 CircRNA–microRNA (ceRNA) network.

Notes: The circRNA–microRNA ceRNA network was drawn with the Cytoscape software. The size of each green node represents the functional connectivity of each circRNA. Rno_circRNA_008008 and rno_circRNA_013779 are the two largest nodes in the network.

Abbreviations: circRNAs, circular RNAs; ceRNA, competing endogenous RNAs.

Discussion

In the present study, we analyzed the circRNA expression in the dorsal horn of the spinal cord, which is an important location for pain regulation, synaptic plasticity, and neuropathic pain treatment.5 Comparison between sham-operated and CCI rats revealed expression change of 469 circRNAs. This study indicates that CCI of the peripheral nerve induces a lot of change of circRNAs in the rat spinal dorsal horn.

Despite the potential importance of circRNAs reported in several types of cancer,34 heart diseases,3538 and other diseases, there is no report on the functional roles of circRNAs in the regulation of physical or pathological pain. In this study, in a rat neuropathic pain model, we found 363 circRNAs significantly upregulated and 106 circRNAs significantly downregulated in the spinal dorsal horn compared with normal ones. Three circRNAs (rno_circRNA_013779, rno_circRNA_013779, and rno_circRNA_013779) in CCI rats were upregulated more than 10 fold and the fold change of circRNA_013779 even reached 28. Interestingly, the number of downregulated circRNAs was less than that of the upregulated ones (363 vs 106); in addition, the magnitudes in the CCI group were much smaller. The top downregulated circRNA, rno_circRNA_007512, downregulated only 2.8 times in the CCI spinal dorsal horn. To evaluate the reliability of rat circRNA microarrays, we picked eight circRNAs of relatively big difference in the expression between two groups and with high expression levels in both groups. Expression patterns of these circRNAs were validated by qPCR. Consistent trends were observed between the microarrays and qPCR results (Figure 2).

Recent evidence revealed that circRNAs can function as microRNA sponges22,39 and regulate parent gene expression to affect disease outcomes.37,39 Despite the lack of literature reporting the exact functions of these differential circRNAs, we investigated the potential targets of the top 100 differential circRNAs (all of them were upregulated in the CCI group). In total, by using the miRanda software, 100 circRNAs were identified as the potential microRNA sponge of 505 microRNAs. GO analysis and KEGG pathway annotation were conducted to investigate the functions of related microRNAs.31 GO enrichment analysis revealed that thousands of target genes were involved in the GO term: biological processes, cellular component, and molecular function, which indicates that these differential circRNAs may influence genes of varied biological function. KEGG analysis showed that there were 112 pathways related to the 100 circRNAs, the top three were reported to be related to neuropathic pain: the Hippo signaling pathway,40 MAPK signaling pathway,4143 and endocytosis.44,45

The circRNA–microRNA interaction network analysis was conducted for the top 10 circRNAs in this study. After filtering, we found that seven circRNAs may interact with microRNAs, all of them with more than three predicted target microRNAs. Two of them, circRNA_008008 and circRNA_013779, targeted eight microRNAs. CircRNAs were believed to negatively regulate miRNAs, and contribute substantially to the ceRNA network. It has been reported that ciRS-7, as a circular miR-7 inhibitor, harbors >60 conventional miR-7 binding sites, which is far more than any known linear sponge.46 Sex determining region Y was another identified miRNA sponge and functioned as a miR-138 sponge.47 According to our results, we hypothesized that some of these circRNAs may act as an inhibitor of miRNAs by acting as microRNA sponges and binding them. Our results implied that it is worthwhile to further investigate potential biological functions of these dysregulated circRNAs in the development of nerve constriction-induced neuropathic pain.

In the present study, although we pooled six dorsal horn samples in both groups, repeated microarray detections (ie, three repeats) and a set P threshold will be better to detect and decrease the number of differential circRNAs. In addition, comprehensive studies on circRNA expression at different time points and different nerve tissues (such as the peripheral sciatic nerve and dorsal root ganglion) in different neuropathic pain models are warranted. Literature has indicated that males and females respond differently in chronic pain conditions. Because the results of this study involve male rats,48,49 we must realize that sex difference may influence the outcome of differential circRNAs. Functional studies of specific circRNAs in pain medicine are needed to clarify circRNAs roles in pain regulation.

Conclusion

Neuropathic pain rats displayed a circRNA expression change. Bioinformatics analyses showed that these circRNAs could be miRNA sponges and these microRNAs targeted mRNAs of various biological functions. CircRNA studies will be helpful to reveal molecular mechanisms of neuropathic pain.

Acknowledgments

The authors thank Dr Yi Zhang and Dr Wenjing Zhou for edits and revision of this manuscript. The study was supported by the National Natural Science Foundation of China (grant number 81660201 to Song Cao).

Disclosure

The authors report no conflicts of interest in this work.

References

1.

Van HO, Austin SK, Khan RA, Smith BH, Torrance N. Neuropathic pain in the general population: a systematic review of epidemiological studies. Pain. 2014;155(4):654–662.

2.

Loeser JD. Economic implications of pain management. Acta Anaesthesiol Scand. 1999;43(9):957–959.

3.

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

4.

Ji RR, Chamessian A, Zhang YQ. Pain regulation by non-neuronal cells and inflammation. Science. 2016;354(6312):572–577.

5.

Ji RR, Xu ZZ, Gao YJ. Emerging targets in neuroinflammation-driven chronic pain. Nat Rev Drug Discov. 2014;13(7):533–548.

6.

West SJ, Bannister K, Dickenson AH, Bennett DL. Circuitry and plasticity of the dorsal horn – toward a better understanding of neuropathic pain. Neuroscience. 2015;300:254.

7.

Todd AJ. Neuronal circuitry for pain processing in the dorsal horn. Nat Rev Neurosci. 2010;11(12):823–836.

8.

Pertovaara A. Noradrenergic pain modulation. Prog Neurobiol. 2006;80(2):53–83.

9.

Bannister K, Bee LA, Dickenson AH. Preclinical and early clinical investigations related to monoaminergic pain modulation. Neurotherapeutics. 2009;6(4):703–712.

10.

Hu R, Li L, Li D, et al. Downregulation of Cdh1 signalling in spinal dorsal horn contributes to the maintenance of mechanical allodynia after nerve injury in rats. Mol Pain. 2016;12.

11.

Gangadharan V, Kuner R. Pain hypersensitivity mechanisms at a glance. Dis Model Mech. 2013;6(4):889.

12.

Andersen HH, Duroux M, Gazerani P. MicroRNAs as modulators and biomarkers of inflammatory and neuropathic pain conditions. Neurobiol Dis. 2014;71(11):159–168.

13.

Sakai A, Saitow F, Miyake N, Miyake K, Shimada T, Suzuki H. miR-7a alleviates the maintenance of neuropathic pain through regulation of neuronal excitability. Brain. 2013;136(9):2738–2750.

14.

Kusuda R, Cadetti F, Ravanelli MI, et al. Differential expression of microRNAs in mouse pain models. Mol Pain. 2011;7(1):1–8.

15.

Neumann E, Hermanns H, Barthel F, Werdehausen R, Brandenburger T. Expression changes of microRNA-1 and its targets connexin 43 and brain-derived neurotrophic factor in the peripheral nervous system of chronic neuropathic rats. Mol Pain. 2015;11(1):1–9.

16.

Shao J, Cao J, Wang J, et al. MicroRNA-30b regulates expression of the sodium channel Nav1.7 in nerve injury-induced neuropathic pain in the rat. Mol Pain. 2016;12.

17.

Genda Y, Arai M, Ishikawa M, Tanaka S, Okabe T, Sakamoto A. microRNA changes in the dorsal horn of the spinal cord of rats with chronic constriction injury: a TaqMan® low density array study. Int J Mol Med. 2013;31(1):129–137.

18.

Li H, Huang Y, Ma C, Yu X, Zhang Z, Shen L. MiR-203 involves in neuropathic pain development and represses Rap1a expression in nerve growth factor differentiated neuronal PC12 cells. Clin J Pain. 2015;31(1):36–43.

19.

Wilusz JE, Sharp PA. A circuitous route to noncoding RNA. Science. 2013;340(6131):440–441.

20.

Hansen TB, Jensen TI, Clausen BH, et al. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495(7441):384–388.

21.

Zheng Q, Bao C, Guo W, et al. Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs. Nat Commun. 2016;7:11215.

22.

Thomson DW, Dinger ME. Endogenous microRNA sponges: evidence and controversy. Nat Rev Genet. 2016;17(5):272.

23.

Memczak S, Jens M, Elefsinioti A, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495(7441):333–338.

24.

Wang Y, Yu X, Luo S, Han H. Comprehensive circular RNA profiling reveals that circular RNA100783 is involved in chronic CD28-associated CD8(+)T cell ageing. Immun Ageing. 2015;12(1):1–10.

25.

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

26.

Zhao Z, Li X, Gao C, et al. Peripheral blood circular RNA hsa_circ_0124644 can be used as a diagnostic biomarker of coronary artery disease. Sci Rep. 2017;7:39918.

27.

Li P, Chen S, Chen H, et al. Using circular RNA as a novel type of biomarker in the screening of gastric cancer. Clin Chim Acta. 2015;444:132.

28.

Weng W, Wei Q, Toden S, et al. Circular RNA ciRS-7 – a promising prognostic biomarker and a potential therapeutic target in colorectal cancer. Clin Cancer Res. Epub 2017 Feb 7.

29.

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.

30.

Parent AJ, Tétreault P, Roux M, et al. Descending nociceptive inhibition is modulated in a time-dependent manner in a double-hit model of chronic/tonic pain. Neuroscience. 2015;315(3):70–78.

31.

Su H, Lin F, Xia D, Shen L. Profiling and bioinformatics analyses reveal differential circular RNA expression in radioresistant esophageal cancer cells. J Transl Med. 2016;14(1):225.

32.

Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS. MicroRNA targets in Drosophila. Genome Biol. 2003;5(1):R1.

33.

Pasquinelli AE. MicroRNAs and their targets: recognition, regulation and an emerging reciprocal relationship. Nat Rev Genet. 2012;13(4):271–282.

34.

Li P, Chen H, Chen S, et al. Circular RNA 0000096 affects cell growth and migration in gastric cancer. Br J Cancer. 2017;116(5):626–633.

35.

Vausort M, Salgadosomoza A, Zhang L, et al. Myocardial infarction-associated circular RNA predicting left ventricular dysfunction. J Am Coll Cardiol. 2016;68(11):1247–1248.

36.

Du WW. Foxo3 circular RNA promotes cardiac senescence by modulating multiple factors associated with stress and senescence responses. Eur Heart J. 2017;38(18):1402–1412.

37.

Wang K, Long B, Liu F, et al. A circular RNA protects the heart from pathological hypertrophy and heart failure by targeting miR-223. Eur Heart J. 2016;37(33):2602.

38.

Tang CM, Zhang M, Huang L, et al. CircRNA_000203 enhances the expression of fibrosis-associated genes by derepressing targets of miR-26b-5p, Col1a2 and CTGF, in cardiac fibroblasts. Sci Rep. 2017;7:40342.

39.

Zheng J, Liu X, Xue Y, et al. TTBK2 circular RNA promotes glioma malignancy by regulating miR-217/HNF1β/Derlin-1 pathway. J Hematol Oncol. 2017;10(1):52.

40.

Li N, Lim G, Chen L, et al. Spinal expression of Hippo signaling components YAP and TAZ following peripheral nerve injury in rats. Brain Res. 2013;1535:137–147.

41.

Huang Q, Mao XF, Wu HY, et al. Cynandione A attenuates neuropathic pain through p38β MAPK-mediated spinal microglial expression of β-endorphin. Brain Behav Immun. 2017;62:64–77.

42.

Ji RR, Suter MR. p38 MAPK, microglial signaling, and neuropathic pain. Mol Pain. 2007;3:33.

43.

Chen Y, Li G, Huang LY. p38 MAPK mediates glial P2X7R-neuronal P2Y1R inhibitory control of P2X3R expression in dorsal root ganglion neurons. Mol Pain. 2015;11:68.

44.

Chen W, Lu Z. Upregulated TLR3 promotes neuropathic pain by regulating autophagy in rat with L5 spinal nerve ligation model. Neurochem Res. 2017;42(2):634–643.

45.

Altmann C, Hardt S, Fischer C, et al. Progranulin overexpression in sensory neurons attenuates neuropathic pain in mice: role of autophagy. Neurobiol Dis. 2016;96:294–311.

46.

Zhao ZJ, Shen J. Circular RNA participates in the carcinogenesis and the malignant behavior of cancer. RNA Biol. 2017;14(5):514–521.

47.

Li J, Yang J, Zhou P, et al. Circular RNAs in cancer: novel insights into origins, properties, functions and implications. Am J Cancer Res. 2015;5(2):472–480.

48.

Sorge RE, Mapplebeck JC, Rosen S, et al. Different immune cells mediate mechanical pain hypersensitivity in male and female mice. 2015;18(8):1081–1083.

49.

Mapplebeck JC, Beggs S, Salter MW. Sex differences in pain: a tale of two immune cells. Pain. 2016;157(Suppl 1):S2–S6.

Creative Commons License © 2017 The Author(s). This work is published and licensed by Dove Medical Press Limited. The full terms of this license are available at https://www.dovepress.com/terms.php and incorporate the Creative Commons Attribution - Non Commercial (unported, v3.0) License. By accessing the work you hereby accept the Terms. Non-commercial uses of the work are permitted without any further permission from Dove Medical Press Limited, provided the work is properly attributed. For permission for commercial use of this work, please see paragraphs 4.2 and 5 of our Terms.