Wu Shaogen, Marie Lutz Brianna, Miao Xuerong, Liang Lingli, Mo Kai, Chang Yun-Juan, Du Peicheng, Soteropoulos Patricia, Tian Bin, Kaufman Andrew G, Bekker Alex, Hu Yali, Tao Yuan-Xiang
Department of Obstetrics and Gynecology, Nanjing Drum Tower Hospital, Nanjing University Medical School, Jiangsu, China Department of Anesthesiology, New Jersey Medical School, Rutgers, The State University of New Jersey, Newark, NJ, USA.
Department of Anesthesiology, New Jersey Medical School, Rutgers, The State University of New Jersey, Newark, NJ, USA.
Mol Pain. 2016 Mar 11;12. doi: 10.1177/1744806916629048. Print 2016.
Peripheral nerve injury leads to changes in gene expression in primary sensory neurons of the injured dorsal root ganglia. These changes are believed to be involved in neuropathic pain genesis. Previously, these changes have been identified using gene microarrays or next generation RNA sequencing with poly-A tail selection, but these approaches cannot provide a more thorough analysis of gene expression alterations after nerve injury.
The present study chose to eliminate mRNA poly-A tail selection and perform strand-specific next generation RNA sequencing to analyze whole transcriptomes in the injured dorsal root ganglia following spinal nerve ligation. Quantitative real-time reverse transcriptase polymerase chain reaction assay was carried out to verify the changes of some differentially expressed RNAs in the injured dorsal root ganglia after spinal nerve ligation.
Our results showed that more than 50 million (M) paired mapped sequences with strand information were yielded in each group (51.87 M-56.12 M in sham vs. 51.08 M-57.99 M in spinal nerve ligation). Six days after spinal nerve ligation, expression levels of 11,163 out of a total of 27,463 identified genes in the injured dorsal root ganglia significantly changed, of which 52.14% were upregulated and 47.86% downregulated. The largest transcriptional changes were observed in protein-coding genes (91.5%) followed by noncoding RNAs. Within 944 differentially expressed noncoding RNAs, the most significant changes were seen in long interspersed noncoding RNAs followed by antisense RNAs, processed transcripts, and pseudogenes. We observed a notable proportion of reads aligning to intronic regions in both groups (44.0% in sham vs. 49.6% in spinal nerve ligation). Using quantitative real-time polymerase chain reaction, we confirmed consistent differential expression of selected genes including Kcna2, Oprm1 as well as lncRNAs Gm21781 and 4732491K20Rik following spinal nerve ligation.
Our findings suggest that next generation RNA sequencing can be used as a promising approach to analyze the changes of whole transcriptomes in dorsal root ganglia following nerve injury and to possibly identify new targets for prevention and treatment of neuropathic pain.
周围神经损伤会导致受损背根神经节的初级感觉神经元基因表达发生变化。这些变化被认为与神经性疼痛的发生有关。此前,这些变化已通过基因微阵列或采用多聚腺苷酸尾选择的新一代RNA测序技术得以确定,但这些方法无法对神经损伤后的基因表达改变进行更全面的分析。
本研究选择去除mRNA多聚腺苷酸尾选择,并进行链特异性新一代RNA测序,以分析脊神经结扎后受损背根神经节中的全转录组。开展定量实时逆转录聚合酶链反应试验,以验证脊神经结扎后受损背根神经节中一些差异表达RNA的变化。
我们的结果显示,每组均产生了超过5000万(M)条带有链信息的配对比对序列(假手术组为5187万 - 5612万条,脊神经结扎组为5108万 - 5799万条)。脊神经结扎6天后,在受损背根神经节中总共27463个已鉴定基因中,11163个基因的表达水平发生了显著变化,其中52.14%上调,47.86%下调。在蛋白质编码基因中观察到最大的转录变化(91.5%),其次是非编码RNA。在944个差异表达的非编码RNA中,长散在非编码RNA的变化最为显著,其次是反义RNA、加工转录本和假基因。我们观察到两组中均有相当比例的 reads 比对到内含子区域(假手术组为44.0%,脊神经结扎组为49.6%)。使用定量实时聚合酶链反应,我们证实了脊神经结扎后所选基因(包括Kcna2、Oprm1以及长链非编码RNA Gm21781和4732491K20Rik)的一致差异表达。
我们的研究结果表明,新一代RNA测序可作为一种有前景的方法,用于分析神经损伤后背根神经节全转录组的变化,并可能为神经性疼痛的预防和治疗确定新的靶点。