MicroRNA And Circular RNA Expression In Affected Skin Of Patients With Postherpetic Neuralgia.

Mechanisms of postherpetic neuralgia (PHN) are still not clear. Transcripts such as microRNA (miRNA) and circular RNA (circRNA) in the affected skin may take part in the initiation and development of this neuropathic pain; however, their expression profiles in skins of PHN patients have not been reported. The PHN affected skin and the mirror skin were collected and subjected to miRNA and circRNA microarray, and expression profiles were comparatively analyzed. There were 317 differently expressed miRNAs in PHN affected skin compared with mirror skin (fold change ≥2.0), and 13 of them showed fold change >10 in the PHN skin. Only one circRNA, hsa_circRNA_405463 showed fold change >2 in PHN skin, however, 31 circRNAs with fold change ≥1.5. To evaluate functions of differential miRNAs, their target mRNAs were predicted and bioinformatics analyses including gene ontology, Kyoto Encyclopedia of Genes and Genomes pathway were conducted. Target mRNAs significantly (P<0.05) enriched in 85 pathways, such as FoxO, AMPK, MAPK and pathway. These data reported for the first time that miRNA and circRNA differentially expressed in the PHN skin and these transcripts with abnormal expression could be potential targets to treat PHN.

Introduction

Postherpetic neuralgia (PHN) is a typical neuropathic pain, which lasts more than 1 month,1 and profoundly affects the quality of life.2 The mechanisms of PHN are complicated, including the peripheral pathology3–5 and plastic changes in the central nervous system.4–8 Elucidating molecular changes during neuropathic pain is fundamental for the development of mechanism-oriented treatments.9 Understanding the molecular changes in PHN skin may help develop skin-targeting external medications or treatments.

MicroRNAs (miRNAs) are small non-coding molecules, which post-transcriptionally modulate pathological processes including neuropathic pain.10,11 Existing studies have shown that miRNAs are involved in the development of neuropathic pain and are expected to become clinical biomarkers and intervention targets.12,13 miRNA expression differences were detected in serum between two kinds of patients: herpes zoster and PHN.14 In patients with peripheral neuropathy, abnormal expression of miR-21-5p, miR-146a, miR-132-3p and miR-155 were reported in the skin, blood leukocytes and sural nerves.15,16 Patients with fibromyalgia syndrome showed abnormal expression of miR-let-7d in skin and leukocytes, which may be associated with lesions in peripheral skin neurofibromas.17

Circular RNAs (circRNAs) are non-coding RNAs which interact with miRNAs.18 circRNAs can act as miRNA sponges19–21 and regulate gene expression through a circRNA-miRNA-mRNA pattern.22,23 To date, there are few reports on the expression change and mechanism of circRNA in pain situation.24–26 We previously found that circRNA expression profile changed in the dorsal spinal cord of CCI rats.27 It is reported that intrathecally silence circHIPK3 with shRNA alleviated neuropathic pain in diabetic rats.24

It seems feasible and effective to interfere with the expression of miRNAs in skin. Animal studies suggested that interfering with the expression of miRNAs in the skin may affect the prognosis of diseases. Skin is an important lesion site and potential drug delivery location3 in PHN; however, the expression profiles of miRNA and circRNA in PHN skin have not been reported. In this study, the differential expression of miRNA and circRNA were detected with microarrays, and differential miRNAs’ target mRNAs were subjected to gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses to evaluate their possible functions and mechanisms in PHN situation.

Methods

Participants

This study was conducted in accordance with the Declaration of Helsinki. The present study was approved by the Medical Ethics Committee of the Affiliated Hospital of Zunyi Medical University and registered at the Chinese Clinical Trial Registry (ChiCTR1800017821) and all patients provided written informed consent. Patients were recruited from the Department of Pain Medicine of the local hospital from September 2018 to March 2019. All PHN patients in this study were diagnosed by two or more chief physicians. Pain intensity was evaluated with an 11-point numeric rating scale (NRS-11, values from 0 to 10). All of the patients recruited claimed intense pain (NRS scores ≥5) and pain durations were more than 1 month after the HZ healed. For easy sampling and to ensure that the skin locations were identical and comparable, only PHN lesion in skin of the thoracic segment was enrolled. Five PHN patients were enrolled in the microarray experiments and described in Table 1.

Table 1

Demographic And Clinical Variables Of Five Patients In Microarray Experiments

No.	Age (Year)	Gender	Location Of Lesion	Pain Duration (Month)	NRS Score
1	65	M	Right T2-4	24.0	7
2	74	M	Right T5-7	1.0	9
3	64	M	Left T6-8	1.0	6
4	77	F	Right T4-6	1.0	7
5	63	F	Left T3-5	12.0	8

Abbreviations: M, male; F, female; T, thoracic; NRS, numerical rating scale.

Skin Collection

For each patient, two circular skins on the back were sampled. Patients took a prone position, local anesthetics (1% lidocaine) was injected at the lesion site and the mirror site (symmetrical site of the opposite side) before the skin was taken with biopsy punchers (with 5 mm diameter and circular cross section). The bottom of the skin samples reaches the dermis, with a thickness of ~2 mm.

RNA Isolation, Purification And Hybridization

Total RNA from each sample was isolated and purified as we reported.28 miRNA labeling and array hybridization was according to Exiqon’s manual. After stopping the labeling procedure, the Hy3™-labeled samples were hybridized on the miRCURYTM LNA Array (v.19.0, Exiqon, USA) according to array manual. Then, the slides were scanned using the Axon GenePix 4000B microarray scanner (Axon Instruments, Foster City, USA).

For circRNA microarray, the RNA preparation and microarray hybridization were performed as we reported27 based on Arraystar’s standard protocols.29 Total RNA from each sample was treated with Rnase R (Epicentre Inc.) to remove linear RNAs and to enrich circRNA. Then, the enriched circRNAs were amplified and transcribed into fluorescent cRNA using a random priming method (Arraystar Super RNA Labeling Kit; Arraystar, USA). The labeled cRNAs were hybridized onto the Arraystar Human circRNA Array (v.2.0, 8×15K, Arraystar). At last, the arrays were scanned by the Agilent Scanner G2505C (Agilent, USA).

Microarray Data Analysis

For miRNAs, scanned images were imported into GenePix Pro 6.0 software (Axon) for grid alignment and data extraction. Replicated miRNAs were averaged and miRNAs that intensities ≥30 in all samples were chosen for calculating normalization factor. Expressed data were normalized using the median normalization. For circRNAs, quantile normalization of raw data and subsequent data processing were performed using the R software limma package. After quantile normalization of the raw data, low-intensity filtering was performed, and the circRNAs that at least 5 of 10 samples have flags in “P” or “M” (“All Targets Value”) were retained for further analyses. After normalization, significant differentially expressed miRNAs and circRNAs between two groups were identified through fold change and P-value. MiRNAs and circRNAs with fold changes ≥2 and P-value ≤0.05 were selected as differentially expressed genes. Hierarchical clustering was used to display the distinguishable expression pattern between PHN and control samples.

Real-Time Quantitative PCR Validation

Five differential miRNAs were randomly selected and subjected to real-time quantitative PCR (RT-qPCR). RNAiso Plus (Takara, Japan) was used to isolate total RNA from skins according to manufacturer’s instructions. Total RNA was prepared as standard template for reverse transcription with the Mir-X miRNA First-Strand Synthesis Kit (Clontech Laboratories, Japan). The miRNA-specific 5ʹ primers were as follows: hsa-miR-3664-3p, 5ʹ-CGTCTCAGGAGTAAAGACAGAGTT-3ʹ; hsa-miR-4714-3p, 5ʹ-CCAACCTAGGTGGTCAGAGTT-3ʹ; hsa-miR-16-5p, 5ʹ-CCGTAGCAGCACGTAAATATTGGCG-3ʹ; hsa-miR-20a-5p, 5ʹ-CGGCTAAAGTGCTTATAGTGCAGGTAG-3ʹ; hsa-miR-505-5p, 5ʹ-ACGGGAGCCAGGAAGTATTGATGT-3ʹ; hsa-let-7a-5p, 5ʹ-CGGCGTGAGGTAGTAGGTTGTATAGTT-3ʹ. RT-qPCR was performed on a CFX Connect Real-Time system (Bio-Rad, USA) in a 25 μL tube containing 12.5 μL of TB Green™ Premix Ex Taq™ II (Takara), and 0.5 μL of mRQ 3ʹ Primer according to manufacturer’s instructions. U6 was used as an endogenous reference. The 2−ΔΔCt method was used to calculate the relative expression of miRNAs.

MiRNA Prediction

Target mRNAs of the top 100 differential miRNAs were predicted with Arraystar’s home-made miRNA target prediction software based on the TargetScan (v.7.1, )30 and miRDB (v.5.0, ) prediction algorithm. Only the mRNAs predicted by two databases were included in the GO and pathway analyses.

Bioinformatics Analysis

GO and pathway analyses were performed to explore the function of differential miRNAs. The mRNAs of top 100 differential miRNAs were subjected to GO annotation in terms of Biological Processes, Cellular Components and Molecular Functions. Pathways defined by KEGG, Biocarta and Reactome () were identified by DAVID ().

Statistical Analysis

The results were reported as mean±SD. Statistically significant differences of miRNA expression between PHN and control group were evaluated with paired t-tests with GraphPad Prism (v.7.0, GraphPad Software, USA). P<0.05 was considered to be statistically significant.

Results

miRNA Profiles

Ten samples from five PHN patients were subjected to miRNA microarray (Figure 1). The miRCURY™ LNA Array (v7.0, Exiqon) detected 2080 human miRNAs (). Overall, 317 of them were calculated as differential miRNAs (fold change ≥2.0, P<0.05), 67 were significantly up-regulated, and the other 250 were down-regulated in PHN group (). The miRNA, hsa-miR-4772-5p, down-regulated 21 times in PHN skin. Thirteen of these differentially expressed miRNAs showed expression fold change >10. The top 10 up- and down-regulated miRNAs are listed in Table 2.

Figure 1

miRNA expression comparison between the PHN skin and normal skin. (A) Scatter plot shows the miRNA expression differences between two groups. The miRNAs above the top green line and below the bottom green line indicated differential miRNAs (fold change ≥2.0 or ≤ −2.0, respectively). (B) The volcano plot shows the distribution of differential miRNAs (indicated with red blocks). Horizontal green line indicated the threshold of P-value, above which indicates P- values <0.05; left and right vertical green line indicated fold change threshold of −2.0 and 2.0, respectively, and between which indicates genes with fold changes <2.0. (C) Hierarchical clustering of top 25 up- and down-regulated miRNAs between groups.

Table 2

Differential miRNAs Between PHN Skin And Mirror Skin (n=5)

MiRNA Name	Fold Change	Intensity (PHN)	Intensity (Control)	P-value
Top 10 up-regulated in PHN
hsa-miR-4491	11.5	498	40	0.018
hsa-miR-502-5p	9.6	304.7	29	0.012
hsa-miR-4528	9.3	116.7	11.5	0.026
hsa-miR-4721	8.5	29.5	3.2	0.042
hsa-miR-760	7.3	171.9	21.8	0.024
hsa-miR-495-3p	6.8	150.7	21	0.012
hsa-miR-382-5p	6.6	1765.4	244.9	0.041
hsa-miR-4506	6.1	295.3	45	0.003
hsa-miR-1258	6.1	1839.7	278.9	0.026
hsa-miR-330-5p	5.9	45.1	7.1	0.018
Top 10 down-regulated in PHN
hsa-miR-4772-5p	−21.0	1.2	33.6	0.011
hsa-miR-2682-5p	−16.7	18.8	293.4	0.005
hsa-miR-3678-3p	−16.2	13.7	205.7	9.3E−07
hsa-miR-3678-5p	−15.8	3.4	51	0.004
hsa-miR-5579-3p	−15.5	97.9	1417.6	0.042
hsa-miR-3664-3p	−15.0	265.4	3719.7	0.015
hsa-miR-4692	−13.2	15.4	241.6	0.041
hsa-miR-4680-3p	−12.5	31.9	373.8	0.003
hsa-miR-3187-3p	−11.9	2.3	31.6	0.033
hsa-miR-518e-3p	−11.7	24.2	266.2	0.000

RT-qPCR Validation

To validate the microarray data, RT-qPCR was employed to detect miRNA expression in PHN skin and normal mirror skin in another cohort of PHN patients (n=12). Five miRNAs were randomly selected. Expression levels detected by the microarray and RT-qPCR are presented in Figure 2A. The expression trends detected by the two methods were consistent, which suggested high reliability of the microarray results.

Figure 2

Validation of five differentially expressed miRNAs by RT-qPCR and gene ontology enrichment analysis. (A) Validation of five differentially expressed miRNAs by RT-qPCR. Data were from 6 to 12 patients. Expression levels were normalized to U6 expression. Bars represent the mean±SD. Paired t-tests. *P<0.05; **P<0.01 compared with control group. (B–D) Enriched gene ontology (GO) terms correspond to the 4077 target mRNAs of the top 100 miRNAs. (B) Top 10 significantly (Sig) enriched biological processes (BP) terms. (C) Top 10 enriched cellular component (CC) terms. (D) Top 10 enriched molecular function (MF) terms.

Bioinformatics Analyses Of Target mRNA Of Differential miRNAs

A total of 4077 target mRNAs of the top 100 miRNAs were predicted by miRDB and TargetScan. These mRNAs were subjected to GO enrichment and the top 20 significantly enriched (P<0.05) GO terms are listed in Figure 2B–D. These mRNAs enriched in Biological Processes, Cellular Component and Molecular Function terms, such as Protein binding (GO: 0005515, with 542 genes) etc.

KEGG pathway analysis showed that the 4077 target mRNAs significantly (P<0.05) enriched in 85 pathways (), and the top 20 pathways are listed in Figure 3. The FoxO (45 mRNAs), ErbB (31 mRNAs), MAPK (78 mRNAs) and AMPK (39 mRNAs) pathway were among the top 10 (Figure 3).

Figure 3

Top 20 KEGG pathways significantly enriched by target mRNAs of top 100 miRNAs. (A) Top 20 pathways significantly enriched (P<0.05) by the 4077 target mRNAs of top 100 miRNAs. (B) The mRNA numbers enriched in the top 20 pathways. Selection counts refer to the count of mRNAs’ entities directly associated with the listed pathway.

Overview Of circRNA Profiles

Ten samples from five PHN patients were subjected to circRNA microarray. The Arraystar Human CircRNA Array v.2.0 detected 12,257 circRNAs (). Hierarchical clustering and scatter plot visualization showed the circRNAs expression levels were different (Figure 4), although only the hsa_circRNA_405463 show differential expression fold change ≥2.0 between the PHN and control group. There were 23 circRNAs were up-regulated and 8 were down-regulated in PHN skin (fold change ≥1.5), and if the threshold was set as 1.2, up to 260 circRNAs were up-regulated and 154 were down-regulated in PHN skin (fold change ≥1.2, ). The top 10 up- and down-regulated circRNAs are listed in Table 3.

Figure 4

circRNA expression comparison between the PHN skin and normal skin. (A) Scatter plot shows the circRNA expression differences between two groups. The circRNAs above the top green line and below the bottom green line indicated differential circRNAs (fold change ≥1.2 or ≤ −1.2, respectively). (B) The volcano plot shows the distribution of differential circRNAs (indicated with red blocks). Horizontal green line indicated the threshold of P-value, above which indicates P- values <0.05; left and right vertical green line indicated fold change threshold of −1.2 and 1.2, respectively, and between which indicates genes with fold changes <1.2. (C) Hierarchical clustering of top 25 up- and down-regulated circRNAs between groups.

Table 3

Differential circRNAs Between PHN Skin And Mirror Skin (n=5)

CircRNA Name	Fold Change	Intensity (PHN)	Intensity (Control)	P-value
Top 10 up-regulated in PHN
hsa_circRNA_101797	1.9	136.9	57.9	0.018
hsa_circRNA_058819	1.6	216.1	109.5	0.045
hsa_circRNA_063013	1.6	114.1	66.4	0.010
hsa_circRNA_404988	1.6	60.3	40.5	0.026
hsa_circRNA_407176	1.6	317.5	158.4	0.009
hsa_circRNA_406848	1.5	70.1	42.8	0.040
hsa_circRNA_001379	1.5	536.7	271.5	0.008
hsa_circRNA_406106	1.5	213.7	116.3	0.036
hsa_circRNA_404432	1.5	392.4	209.9	0.003
hsa_circRNA_102538	1.5	183.1	106.7	0.014
Top 10 down-regulated in PHN
hsa_circRNA_405463	−2.0	1809.7	2057.5	0.007
hsa_circRNA_077409	−1.8	128.2	190.6	0.037
hsa_circRNA_011772	−1.8	82.6	122.9	0.029
hsa_circRNA_101731	−1.8	160.5	235.2	0.007
hsa_circRNA_104137	−1.8	1126.4	1502.8	0.002
hsa_circRNA_104664	−1.7	263.1	347.6	0.026
hsa_circRNA_025988	−1.7	51.5	78.2	0.004
hsa_circRNA_101130	−1.6	161.6	204.4	0.029
hsa_circRNA_101341	−1.6	101.9	142.1	0.002
hsa_circRNA_026457	−1.6	3393.7	3519.8	0.022

Discussion

Skin is the lesion site and the location where hyperalgesia and allodynia appear in PHN patients. For the first time, the expression profiles of miRNAs and circRNAs in the PHN skin were detected with microarrays. Comprehensive expression changes of transcripts, especially miRNAs, were found between PHN and normal skin. This indicates that PHN skin and these differentially expressed molecules can be targets for PHN treatment.

miRNAs can regulate biological processes by inhibiting the function of their target genes. In this study, as many as 317 miRNAs were found differentially expressed in PHN skin, and 13 of them showed expression fold change >10, indicating that these differential miRNAs may influence their target genes’ expression as well as PHN pathologies. To evaluate functions of these miRNAs, their target mRNAs were predicted and KEGG pathway results indicated that they may influence the PHN pathology via AMPK, MAPK, ErbB, FoxO pathway, etc. It has been extensively reported that AMPK,31,32 MAPK13,33,34 and ErbB35,36 pathway involved in chronic and neuropathic pain in animal studies. Although FoxO’s function has not been confirmed in neuropathic pain model, bioinformatics predicted that FoxO pathway in the dorsal spinal cord27 or dorsal root ganglion37 was related to neuropathic pain.

It seems feasible and effective to interfere with the expression of miRNAs in skin. Animal studies suggested that interfering with the miRNA expression in the skin affects the prognosis of diseases. For example, interfering NR1 subunit of the NMDA receptor by intradermally injecting lentiviral vector-encoded miRNA-based shRNA reduced inflammatory pain in rats.38 Intradermally injecting miR-129 and/or miR-335 agomirs into the wound edges accelerated wound healing in rats with diabetes.39 Some miRNAs have been intradermally given and tested in clinical trials.40

Exploration of circRNAs’ roles in pain situation is still in infancy. Only a few studies focused on the expression change and mechanism of circRNA in pain situation.24–27 This study added evidence and suggested that PHN induced expression changes of circRNA in peripheral tissue. Further studies with gain- and/or loss-of-function experiments of these miRNAs and circRNAs will help understand their functions in PHN situation.

It is hard to tell whether the differential miRNAs and circRNAs are from the nerve endings of DRG neurons or cells that make up the skin. It is also not certain that these transcripts are present in cells because neurons and other cells may release miRNAs41 and circRNAs42 to the extracellular space via extracellular vesicles, such as exosomes. In situ hybridization experiments may help explain this question. Nevertheless, interventions targeting these differential miRNAs or circRNAs at the affected skin may help to alleviate PHN.

Conclusions

MiRNAs and circRNAs differentially expressed in the skin of PHN patients. Target genes of differential miRNAs enriched into pathways related to pathological pain regulation. These differential miRNAs and circRNAs as well as these predicted pathways could be potential targets to treat PHN.