全血miRNA候选优先排序用于纵向疼痛变化提示早期类风湿关节炎中突触和神经元组织通路参与

《Molecular Pain》：Whole-blood miRNA candidate prioritization for longitudinal pain change suggests synaptic and neuronal organization pathway involvement in early rheumatoid arthritis

　　

疼痛是慢性炎症性疾病的主要临床结局，然而个体内疼痛变化的分子相关性仍不清楚，且必须与炎症变化区分开来。研究人员利用RA-MAP/TACERA早期类风湿关节炎队列的纵向全血多组学数据，优先筛选其6个月变化与视觉模拟量表疼痛变化（ΔVAS）共变的微RNA（miRNA）候选。在151名参与者和1,881个miRNA特征中，研究人员应用了一种稳定性导向的delta组学工作流程，结合子采样弹性网络选择频率、自举方向过滤以及针对ΔCRP和ΔESR调整的选择后秩残差偏相关。这确定了29个方向稳定的miRNA候选；在ΔCRP/ΔESR调整后，16个在事后BH-FDR<0.05下在这些29个候选中保留了符号一致的支撑。匹配的miRNA-mRNA delta分析未发现任何单个miRNA-mRNA配对能通过配对水平的BH-FDR校正。在主要的p<0.05名义配对选择阈值下，涉及117个靶基因的123个反向配对产生了探索性基因本体生物学过程富集，涉及突触和神经元组织；然而，当配对选择阈值收紧至p<0.01时，没有GO术语保持BH显著。探索性血浆慢性颈痛比较在39个符号一致检测匹配中产生22个（56.4%；精确95% CI，39.6%–72.2%），与50%的随机期望无统计学差异，且没有匹配通过BH-FDR校正。这些发现概述了一种稳定性导向、炎症调整的方法，用于优先筛选与疼痛变化相关的血液miRNA候选。这16个miRNA子集应被解释为产生假设的候选池，而非经过验证的生物标志物组合，并为未来的纵向验证和功能研究提供了聚焦基础。

慢性炎症性疾病中，疼痛是核心临床结局，但主观量表无法区分疼痛变化源于生物学缓解、中枢适应还是心理因素。血液因易于反复采样，且包含免疫、血管和神经免疫信号，成为寻找分子相关性的理想样本。既往研究多为横断面或基线预测，难以揭示个体内疼痛变化。微RNA（miRNA）在血液中稳定存在，可调控多个靶mRNA，是高维转录变化的紧凑代表。然而，在类风湿关节炎（RA）背景下，循环miRNA和疼痛均受炎症影响，需调整炎症标志物。RA-MAP/TACERA早期RA纵向队列提供了基线和6个月随访的全血多组学数据，使研究人员能够开发一种稳定性导向、炎症调整的delta（Δ，变化量）组学方法，优先排序与疼痛变化相关的miRNA候选。该研究发表在《Molecular Pain》。

研究人员利用151名参与者和1,881个miRNA特征，设计四阶段工作流程。主要技术方法包括：子采样弹性网络（subsampling elastic-net）选择频率、自举（bootstrap）方向过滤、留一法（jackknife）敏感性分析、针对ΔCRP和ΔESR的秩残差偏相关、经验贝叶斯收缩（ashr）、multiMiR数据库验证靶标的匹配miRNA-mRNA反向共变分析、基因本体（GO）生物学过程和KEGG通路富集，以及外部慢性颈痛血浆数据集的符号一致性比较。样本来源为RA-MAP/TACERA队列（全血）和一个外部9例慢性颈痛血浆数据集。

结果如下。

Phase 1: miRNA-only候选优先排序。通过全数据集Spearman相关、子采样弹性网络选择频率≥30%且名义p<0.05、以及自举95%置信区间不跨零，从1,881个特征中筛选出29个方向稳定的候选miRNA。最强的候选是hsa-mir-6852（ρ=0.319，选择频率96%）。所有29个候选在留一法删除任意单个参与者后均保持方向不变。

Phase 2: 炎症标志物调整的敏感性分析。对ΔCRP和ΔESR进行秩残差偏相关调整后，29个候选全部保留方向，其中20个保持名义显著，16个在29个候选内的事后BH-FDR<0.05下保留支持，构成16-miRNA核心子集。临床活动性调整（加入ΔTJC（压痛关节计数）和ΔSJC（肿胀关节计数））后，FDR保留候选降至7个；血细胞计数调整后保留16个，其中15个与核心子集重叠。经验贝叶斯收缩显示，16个核心候选的效应量大幅衰减（中位相对衰减99.05%），但方向全部保留，且无候选达到局部错误符号率（LFSR）<0.05。基于16个候选符号加权的集合评分与ΔVAS相关（ρ=0.650；置换p=9.999×10^-5），在置换检验中显著。

Phase 3: 靶标水平和通路水平的探索性生物学背景。在匹配的151名参与者miRNA-mRNA delta数据中，4,967个数据库验证的配对无一个通过配对水平BH-FDR校正。在p<0.05名义阈值下，123个反向配对涉及117个靶基因，GO生物学过程富集显示突触组织、细胞连接组织、神经元突出形态发生和轴突发生等术语显著（如突触组织GO:0050808，BH调整p=0.00318）。然而，当配对选择阈值收紧至p<0.01时，仅24个反向配对，无GO术语保持BH显著。KEGG通路无显著结果。

Phase 4: 探索性外部疼痛变化比较。外部慢性颈痛血浆数据集中，16个核心miRNA全部匹配到39个检测列。其中22个（56.4%）与Phase 1方向一致，但精确95% CI包含50%，与随机无差异；无匹配通过BH-FDR校正。

讨论部分总结：研究的主要贡献是提出一种稳定性导向、炎症调整的delta多组学工作流程，而非验证具体生物标志物。16-miRNA核心子集应视为产生假设的候选池，而非验证面板。由于无miRNA-mRNA配对通过FDR校正，且外部比较无显著一致性，这些结果强调需要在更大的独立纵向队列中进行评估。未来方向包括独立发现/验证队列、针对性定量PCR或数字PCR检测，以及miRNA模拟物或抑制剂的功能研究。研究结论翻译如下：总之，这项工作概述了一种稳定性导向、炎症调整的delta多组学方法，用于优先排序与疼痛变化相关的全血miRNA候选。16-miRNA核心子集应被解释为产生假设的候选池，而非经过验证的生物标志物组合。由于没有单个miRNA-mRNA配对通过配对水平FDR校正，且外部比较既未显示统计学上可检测的集合符号一致性（高于随机）也无任何FDR显著的匹配，主要意义是一个优先排序的miRNA候选集，需要在更大的独立delta队列中评估。