术后慢性疼痛是指由手术引起、继发于术后急性疼痛且持续时间超过3个月的慢性疼痛[1]。其病因复杂，机制未明，尚缺乏科学的干预手段和有效治疗方案，对动物机体康复不利，导致医疗资源浪费，还严重影响生活质量，产生不良后果[2]。为明析术后慢性疼痛的机制和致病因素，实现动物“疼痛最小化”，兽医工作者进行了大量研究。术后疼痛的实质是神经递质和受体等物质参与的多维神经调节的结果，涉及多种物质变化[3]，但目前尚不明确究竟何种物质在其发生、传递和调制的过程中起决定作用。严重的术后急性疼痛是术后慢性疼痛的重要预测因素[4-5]，因此，甄别引发术后急性到慢性疼痛转化这一复杂过程的关键物质，对术后慢性疼痛机制的深入研究和治疗靶点的选择尤为重要。
近年来，代谢组学技术的应用促进了疼痛相关物质研究的发展[6-8]，因此，本试验拟用GC-MS代谢组学技术探讨大鼠术后急性疼痛和慢性坐骨神经紧缩损伤疼痛模型(CCI)的相关代谢变化，以期筛选出术后急性向慢性疼痛转化的关键物质及其相关代谢通路，为术后慢性疼痛的研究奠定基础。
1 材料与方法 1.1 仪器与试剂Agilent GC-MS气质联用仪(7890B-5977A)；Agilent DB-5MS色谱柱(30 m×0.25 mm×0.25 μm)；电子测痛仪；L-2-氯苯丙氨酸(上海恒创生物)；HPLC级氯仿、甲醇和正己烷(CNW，德国)；超纯水由Milli-Q超纯水仪制备，其余试剂均为市售分析纯。
1.2 方法1.2.1 实验动物分组及模型构建 SPF级SD大鼠65只，体重(250±10)g，南方医科大学实验动物中心提供，饲养于华南农业大学实验动物中心。大鼠随机分为空白对照组(n=5)、假手术组(n=30)和手术组(n=30)。参照Bennett等[9]的方法制作大鼠CCI模型，具体方法：大鼠经腹腔注射水合氯醛麻醉后剔除右侧大腿被毛，消毒备皮。平行股骨外侧作长约1 cm切口，拨开皮肤、筋膜和肌肉暴露坐骨神经主干。挑起坐骨神经主干并结扎，共4道，间距2 mm，结扎力度以小腿肌肉轻微颤动为宜，关闭切口；假手术组暴露坐骨神经但不结扎，其余操作同手术组；空白对照组不作任何处理。
1.2.2 机械缩足反射阈值测量 Von Frey电子测痛仪测定大鼠机械缩足反射阈值(mechanical withdrawal threshold，MWT)，大鼠于测痛装置内适应环境30 min，将Von Frey刺激大鼠术侧足底，力度由小到大，直至出现躲避反射即为其MWT值，重复测定5次，结果取平均值。
1.2.3 样本采集与预处理 术前和术后第1、3、5、7、14、21天，大鼠经乙醚麻醉后，断头法处死，取下丘脑及胸椎T3~T5段脊髓经PBS冲洗后，液氮中保存。称取样品于离心管中，加入内标L-2-氯苯丙氨酸20 μL和600 μL甲醇水溶液，-80 ℃放置2 min，研磨；加入120 μL氯仿。冰水浴超声提取，4 ℃静置。离心10 min(12 000 r·min-1, 4 ℃)，吸取400 μL上清液转入衍生瓶中，挥干样本。加入甲氧胺盐酸盐吡啶溶液80 μL进行肟化反应，随后加入80 μL BSTFA衍生试剂和20 μL的正己烷。70 ℃孵育60 min，室温冷却后用于GC-MS检测。每10个样品中插入1个由所有样本的提取液等体积混合制备的质控样品(QC)。
1.2.4 气相色谱-质谱检测及数据分析 GC-MS检测条件参考文献[10]的方法并结合本试验适当更改，质谱的离子源温度为230 ℃，质量扫描范围为50~500 m·z-1；气相色谱进样口温度为260 ℃, 升温程序见表 1。GC-MS的原始数据由Anal-ysisBaseFileConverter和MS-DIAL软件进行预处理，将手术组或假手术组在相同取样时间的样品数据归为同一组，即手术组(D1、D3、D5、D7、D14、D21组)，假手术组(S1、S3、S5、S7、S14、S21组)；将数据导入SIMCA-P软件中进行多元统计分析，找到组间差异代谢物，并将筛选出的差异代谢物以KEGG数据库为背景，经MBRole进行代谢通路富集分析。
|
表 1 气相色谱升温程序 Table 1 Temperature program of GC |
大鼠机械缩足反射阈值结果见图 1，结果显示，大鼠MWT术前基础值比较无显著差异。模型建立后，手术组大鼠MWT较术前显著降低，随后逐渐恢复，直至术后第21天MWT值仍低于术前基础值，差异显著(P < 0.05)，提示CCI模型建模成功。
|
图 1 大鼠机械缩足反射阈值 Fig. 1 The mechanical withdrawal threshold of rat |
质控样品的TIC叠加图及单个代表性样品TIC图见图 2，结果表明，采用的分析方法具有良好的重复性和稳定性，分析程序合理。
|
A. 质控样品TIC图；B. 手术组大鼠样品TIC图 A. TIC of quality control samples; B. TIC of surgical group sample 图 2 样品的总离子流色谱图 Fig. 2 Total ions chromatogram |
样本PCA模型得分图如图 3所示，QC样本紧密分布，表明试验具有较好的稳定性和重复性。为更直观观察组间差异，进一步建立每两分析组间共11组的PCA、PLS-DA、OPLS-DA模型。PCA得分图表明，组间分离趋势明显；PLS-DA模型解释能力参数R2Y(0.842~0.995)，预测能力参数Q2 (0.584~0.937)；OPLS-DA模型参数R2Y(0.842~ 0.991)，Q2(0.474~0.805)；在建立的OPLS-DA模型下进行200次响应排序检验考察模型的质量，各组模型预测能力参数Q2均在-0.925~-0.156；除D21/S21组的PLS-DA模型参数(R2X=0.491，R2Y=0.842，Q2=0.181)外，模型稳定性、预测能力均符合分析要求，模型未出现过拟合。以手术组D3和D5组两组样本的分析为例，见图 4。
|
横坐标表示第1主成分，即PC1，用t[1]表示；纵坐标表示第2主成分，即PC2，用t[2]表示 X-axis show the first principal component, t[1], Y-axis show the second principal component, t[2] 图 3 所有分析样本PCA得分图 Fig. 3 PCA score plots of all sample group |
|
A、B、C分别代表样本的PCA、PLS-DA、OPLS-DA得分图；方形代表D3组的5个样本，三角形代表D5组的5个样本；得分图的横坐标表示第1主成分，即PC1，用t[1]表示；纵坐标表示第2主成分，即PC2，用t[2]表示；D是OPLS-DA模型的响应排序检验，纵坐标表示R2和Q2的取值 A, B, C show the sample score plots of PCA, PLS-DA, OPLS-DA, square mean 5 samples in the D3 group and triangles mean 5 samples in the D5 group of surgical group; X-axis show the first principal component, t[1], Y-axis show the second principal component, t[2]. D mean the permutation test for OPLS-DA model, Y-axis mean the value of R2 and Q2 图 4 两组样本得分图 Fig. 4 Sample score plots of two group |
经多元统计分析、T检验和变异倍数分析(VIP＞1，P＜0.05，∣log2FC∣≥2)，从224种代谢物中筛选出35种代谢物作为候选差异代谢物，其中，24个显著下调，9个显著上调，见表 2。为表征候选差异代谢物间的聚落关系，对代谢物的定量信息进行层次聚类，见图 5；在候选差异代谢物中，N-乙酰天冬氨酸、泛酸、天冬氨酸、葡萄糖、β-丙氨酸、3-羟基丁酸与疼痛密切相关，可能是术后急性疼痛向慢性疼痛转化的潜在生物标志物，为表征代谢物在不同时间的变化趋势，对6种代谢物的平均表达量绘制箱线图，见图 6。
|
|
表 2 主要差异代谢物及其相关代谢通路 Table 2 The main differential metabolites and related metabolic pathways |
|
图 5 各组分析样本差异代谢物热图 Fig. 5 Heatmap of differential metabolites in each sample group |
|
A. N-乙酰天冬氨酸; B. β-丙氨酸; C. 3-羟基丁酸; D. 天冬氨酸; E. 葡萄糖; F. 泛酸 A. N-acetylaspartic acid; B. Beta-alanine; C. 3-hydroxybutyric acid; D. Aspartate; E. Glucose; F. Pantothenic acid 图 6 代谢物平均表达量箱线图 Fig. 6 Box plots of average metabolites expression |
将差异代谢物映射到KEGG数据库中，通过MBRole通路分析功能进行通路富集分析，结果35种候选差异代谢物参与了45个代谢途径。代谢通路的P值表示该代谢通路富集的显著性，以10为底取对数绘制代谢通路富集图，P值越小，-lg(P-value)越大；图中红线和蓝线分别示意P值为0.01和0.05(图 7)。
|
A.D 1组vs D 3组；B. D 3组vs D 5组；C. D 5组vs D 7组；D. D 7组vs D 14组；E. D14组vs D21组 A.Group D 1 vs group D 3; B. Group D 3 vs group D 5; C. Group D 5 vs group D 7; D. Group D 7 vs group D 14; E. Group D 14 vs group D 21 图 7 代谢通路富集图 Fig. 7 Enrichment map of metabolic pathways |
本试验成功建立了大鼠术后急性疼痛和慢性坐骨神经紧缩损伤疼痛模型，并应用GC-MS代谢组学技术比较手术组与对照组、手术组D1、D3、D5、D7、D14、D21各组间的差异代谢物，鉴定出包括N-乙酰天冬氨酸(NAA)、泛酸、天冬氨酸、β-丙氨酸、3-羟基丁酸、葡萄糖等35种差异代谢物及45条潜在代谢通路可能参与了术后急性向慢性疼痛的转化。
N-乙酰天冬氨酸是衡量神经细胞破坏程度的标志物[11]，浓度降低提示神经元丢失、结构破坏或功能异常[12-13]。应用核磁共振波谱技术(MRS)，发现慢性疼痛患者经颅直流电刺激治疗后，疼痛程度减轻，且脑代谢物中NAA浓度较治疗前显著升高[14-15]。进一步研究发现NAA的浓度与疼痛的强度呈负相关[16]，本试验中，NAA浓度在术后D1 vs D3组和D14 vs D21组中表达下调，与这些研究结果相同，其平均表达量与MWT值变化趋势一致，说明NAA在急性向慢性疼痛转化的过程中发挥了作用。
泛酸以辅酶A的生物活性形式挥发其功能，参与含碳架物质的能量代谢[17]，其表达失调表明机体能量失衡。神经损伤后释放大量ATP，作用于P2X3受体参与伤害性信息的传递和痛觉敏化[18-19]。本试验中泛酸浓度下调，其代谢紊乱引起的能量失衡可能导致信号传导异常或细胞内外离子交换障碍，从而在术后急性疼痛向慢性疼痛转化过程发挥作用。辅酶A也是机体应激反应重要观测指标皮质醇的必需辅助因子[20]，伤害性信号促使下丘脑-垂体-肾上腺皮质轴兴奋增强，皮质醇等应激激素水平升高，引起术后应激反应，进而导致糖代谢紊乱[21]。葡萄糖是脑细胞能量代谢的最主要的直接来源，本试验中葡萄糖浓度下调，可能是手术伤害性刺激引起神经元的活性增强、代谢旺盛，从而参与疼痛的形成。
β-丙氨酸为不参与蛋白质合成的氨基酸，与主要的抑制性神经递质甘氨酸和γ-氨基丁酸结构相似，能够激活甘氨酸和γ-氨基丁酸受体，进而影响突触的可塑性[22-23]。伤害感受器可塑性变化在急性疼痛向慢性疼痛转化中具有重要作用[24]。试验结果发现β-丙氨酸代谢紊乱，说明其可能通过调控神经系统的可塑性变化从而参与了术后急性向慢性疼痛转变的过程。
天冬氨酸是重要的兴奋性氨基酸类神经递质，参与维持神经系统的结构与功能，浓度过高会介导神经元的死亡。本试验天冬氨酸浓度升高，提示手术引起的伤害性刺激传入脊髓背角使神经突触和胶质细胞释放大量的天冬氨酸[25]，天冬氨酸等兴奋性氨基酸过度释放会导致NMDA受体上调，使突触的可塑性出现长时程易化，造成中枢敏化[26]，从而参与了术后慢性疼痛的形成。
3-羟基丁酸具有多种生物活性，可以为神经胶质细胞提供能量，并通过抗氧化、线粒体保护和增加脑源性神经营养因子的表达等方式对神经起到保护作用[27]，同时还可以抑制神经胶质细胞的凋亡。胡华辉等[28]发现3-羟基丁酸在急性脊髓损伤模型大鼠脊髓中浓度降低，与本试验中3-羟基丁酸在手术组D7 vs D14组中表达下调相似，说明其参与了疼痛形成过程中神经保护与修复的过程。3-羟基丁酸还可以作为表观遗传修饰调控因子发挥功能，预示其可能在疼痛的表观遗传调控中发挥功能[29-30]。
4 结论使用GC-MS代谢组学技术对大鼠术后急性疼痛和慢性坐骨神经紧缩损伤疼痛模型的下丘脑和脊髓代谢物进行初步研究，筛选出N-乙酰天冬氨酸、泛酸、天冬氨酸、β-丙氨酸、3-羟基丁酸、葡萄糖等35种物质可能与术后急性疼痛向慢性疼痛转变密切相关。但由于试验条件和试验设计的限制，这些差异代谢物在术后疼痛中的作用及其具体机制仍有待进一步深入研究。
| [1] | SCHUG S A, LAVAND'HOMME P, BARKE A, et al. The IASP classification of chronic pain for ICD-11[J]. Pain, 2019, 160(1): 45–52. DOI: 10.1097/j.pain.0000000000001413 |
| [2] |
远立国, 熊惠军, 贾坤, 等. 动物疼痛识别与评估研究进展[J]. 畜牧兽医学报, 2014, 45(10): 1592–1599.
YUAN L G, XIONG H J, JIA K, et al. Advances in animal identification and assessment of pain[J]. Acta Veterinaria et Zootechnica Sinica, 2014, 45(10): 1592–1599. DOI: 10.11843/j.issn.0366-6964.2014.10.003 (in Chinese) |
| [3] |
徐涛, 陈画菡, 谌颜, 等. P物质与犬术后急性疼痛强度关系[J]. 中国兽医学报, 2019, 39(2): 303–306.
XU T, CHEN H H, SHEN Y, et al. Relationship of substance P and acute pain intensity in dogs[J]. Chinese Journal of Veterinary Science, 2019, 39(2): 303–306. (in Chinese) |
| [4] | CHAPMAN C R, VIERCK C J. The transition of acute postoperative pain to chronic pain: an integrative overview of research on mechanisms[J]. J Pain, 2017, 18(4): 351. e1–359. e38. |
| [5] | BUVANENDRAN A, VALLE C J D, KROIN J S, et al. Acute postoperative pain is an independent predictor of chronic postsurgical pain following total knee arthroplasty at 6 months: a prospective cohort study[J]. Reg Anesth Pain Med, 2019, 44(3): 287–296. |
| [6] | PATTI G J, YANES O, SHRIVER L P, et al. Metabolomics implicates altered sphingolipids in chronic pain of neuropathic origin[J]. Nat Chem Biol, 2012, 8(3): 232–234. DOI: 10.1038/nchembio.767 |
| [7] | FINCO G, LOCCI E, MURA P, et al. Can urine metabolomics be helpful in differentiating neuropathic and nociceptive pain? A proof-of-concept study[J]. PLoS One, 2016, 11(3): e0150476. DOI: 10.1371/journal.pone.0150476 |
| [8] | WANG H T, LIU W, CAI Y H, et al. Glutaminase 1 is a potential biomarker for chronic post-surgical pain in the rat dorsal spinal cord using differential proteomics[J]. Amino Acids, 2016, 48(2): 337–348. DOI: 10.1007/s00726-015-2085-z |
| [9] | BENNETT G J, XIE Y K. A peripheral mononeuropathy in rat that produces disorders of pain sensation like those seen in man[J]. Pain, 1988, 33(1): 87–107. DOI: 10.1016/0304-3959(88)90209-6 |
| [10] |
杨冰玉, 李晓莉, 姚遥, 等. 基于GC-MS的慢性髓系白血病代谢轮廓分析[J]. 中国实验血液学杂志, 2017, 25(6): 1585–1591.
YANG Y B, LI X L, YAO Y, et al. Metabolomics of chronic myelogenous leukemia based on GC-MS[J]. Journal of Experimental Hematology, 2017, 25(6): 1585–1591. (in Chinese) |
| [11] |
章容, 戴武, 张继, 等. 初诊2型糖尿病患者胰岛素强化治疗前后海马区NAA/Cr、Cho/Cr的变化[J]. 安徽医科大学学报, 2015, 50(4): 515–518.
ZHANG R, DAI W, ZHANG J, et al. Changes of Cho/Cr and NAA/Cr in hippocampus, before and after intensive insulin therapy in newly diagnosed type 2 diabetic patients[J]. Acta Universitatis Medicinalis Anhui, 2015, 50(4): 515–518. (in Chinese) |
| [12] |
曹琳玉, 丁虹琬, 姚开云, 等. N-乙酰天冬氨酸在中枢神经系统中的研究进展[J]. 药学学报, 2020, 55(8): 1754–1761.
CAO L Y, DING H W, YAO K Y, et al. Progress in research of N-acetylaspartate in central nervous system[J]. Acta Pharmaceutica Sinica, 2020, 55(8): 1754–1761. (in Chinese) |
| [13] | KAMEDA T, FUKUI S, TOMINAGA R, et al. Brain metabolite changes in the anterior cingulate cortex of chronic low back pain patients and correlations between metabolites and psychological state[J]. Clin J Pain, 2018, 34(7): 657–663. DOI: 10.1097/AJP.0000000000000583 |
| [14] | CHANG L D, MUNSAKA S M, KRAFT-TERRY S, et al. Magnetic resonance spectroscopy to assess neuroInflammation and neuropathic pain[J]. J Neuroimmune Pharmacol, 2013, 8(3): 576–593. DOI: 10.1007/s11481-013-9460-x |
| [15] | AUVICHAYAPAT P, KEERATITANONT K, JANYACHAREON T, et al. The effects of transcranial direct current stimulation on metabolite changes at the anterior cingulate cortex in neuropathic pain: a pilot study[J]. J Pain Res, 2018, 11: 2301–2309. DOI: 10.2147/JPR.S172920 |
| [16] | WIDERSTRÖM-NOGA E, GOVIND V, ADCOCK J P, et al. Subacute pain after traumatic brain injury is associated with lower insular N-acetylaspartate concentrations[J]. J Neurotrauma, 2016, 33(14): 1380–1389. DOI: 10.1089/neu.2015.4098 |
| [17] | GHEITA A A, GHEITA T A, KENAWY S A. The potential role of B5:a stitch in time and switch in cytokine[J]. Phytother Res, 2020, 34(2): 306–314. DOI: 10.1002/ptr.6537 |
| [18] | XIANG X E, WANG S S, SHAO F B, et al. Electroacupuncture stimulation alleviates CFA-Induced inflammatory pain via suppressing P2X3 expression[J]. Int J Mol Sci, 2019, 20(13): 3248. DOI: 10.3390/ijms20133248 |
| [19] | BELE T, FABBRETTI E. The scaffold protein calcium/calmodulin-dependent serine protein kinase controls ATP release in sensory ganglia upon P2X3 receptor activation and is part of an ATP keeper complex[J]. J Neurochem, 2016, 138(4): 587–597. DOI: 10.1111/jnc.13680 |
| [20] |
张宏如, 陶嘉磊, 卢圣锋, 等. 从针刺治疗应激性疾病论针刺与应激的相关性[J]. 针刺研究, 2016, 41(6): 560–564.
ZHANG H R, TAO J L, LU S F, et al. From treatment of stress disease by acupuncture to correlation between acupuncture and stress[J]. Acupuncture Research, 2016, 41(6): 560–564. (in Chinese) |
| [21] |
邢丽娇. 手术伤害性刺激模型的18F-FDG PET脑功能成像研究[D]. 北京: 北京协和医学院, 2016.
XING L J. Brain functional imaging study of surgical induced nociception with 18F-FDG PET[D]. Beijing: Peking Union Medical College, 2016. (in Chinese) |
| [22] | SEINO Y, OHASHI N, KOHNO T. The endogenous agonist, β-alanine, activates glycine receptors in rat spinal dorsal neurons[J]. Biochem Biophys Res Commun, 2018, 500(4): 897–901. DOI: 10.1016/j.bbrc.2018.04.183 |
| [23] |
刘钰, 何永恒. 疼痛机制的分子生物学研究进展[J]. 中华中医药学刊, 2017, 35(2): 373–377.
LIU Y, HE Y H. Review on recent studies in molecular biology of pain mechanism[J]. Chinese Archives of Traditional Chinese Medicine, 2017, 35(2): 373–377. (in Chinese) |
| [24] | LUO C, KUNER T, KUNER R. Synaptic plasticity in pathological pain[J]. Trends Neurosci, 2014, 37(6): 343–355. DOI: 10.1016/j.tins.2014.04.002 |
| [25] |
张美玉, 焦玥, 刘洋, 等. 川芎嗪对神经病理性疼痛模型大鼠的镇痛作用及对脑中枢敏化区氨基酸类神经递质含量的影响[J]. 中国药理学与毒理学杂志, 2018, 32(12): 921–929.
ZHANG M Y, JIAO Y, LIU Y, et al. Effect of ligustrazine on analgesic effect and amino acid transmitters in central sensitizations of brain regions in neuropathic pain model rats[J]. Chinese Journal of Pharmacology and Toxicology, 2018, 32(12): 921–929. DOI: 10.3867/j.issn.1000-3002.2018.12.003 (in Chinese) |
| [26] |
崔东, 李泽华, 宋学军. 慢性疼痛的脊髓机制[J]. 中国疼痛医学杂志, 2017, 23(9): 641–647.
CUI D, LI Z H, SONG X J. Spinal mechanisms of chronic pain[J]. Chinese Journal of Pain Medicine, 2017, 23(9): 641–647. DOI: 10.3969/j.issn.1006-9852.2017.09.001 (in Chinese) |
| [27] | WANG L L, YANG Z Q, FAN F, et al. PHBHHx facilitated the residence, survival and stemness maintain of transplanted neural stem cells in traumatic brain injury rats[J]. Biomacromolecules, 2019, 20(9): 3294–3302. DOI: 10.1021/acs.biomac.9b00408 |