《Frontiers in Molecular Neuroscience》：Peripheral exosomal miRNAs as candidate molecular indicators associated with acupuncture at ST36 in an LPS-induced neuroinflammation model
编辑推荐：
本研究旨在通过经验贝叶斯（empirical Bayes）重分析方法，在脂多糖（LPS）诱导的神经炎症模型中，鉴定经足三里（ST36）针刺治疗后呈现类对照表达模式的外周循环外泌体microRNA（miRNA）候选分子。研究人员对来自C57BL/6雄性小鼠的两套
本研究旨在通过经验贝叶斯（empirical Bayes）重分析方法，在脂多糖（LPS）诱导的神经炎症模型中，鉴定经足三里（ST36）针刺治疗后呈现类对照表达模式的外周循环外泌体microRNA（miRNA）候选分子。研究人员对来自C57BL/6雄性小鼠的两套独立Affymetrix GeneChip miRNA 4.0微阵列数据集，采用limma经验贝叶斯框架并辅以Benjamini-Hochberg错误发现率（FDR）校正进行重新分析。经检测调用（detection-call）过滤后，共鉴定出286个共同表达的探针。类对照模式候选分子需同时满足三项标准：LPS效应的FDR < 0.05；LPS+ST36与CON比较的FDR > 0.05；且logFC方向与LPS效应方向相反。随后进行了靶基因预测与通路富集分析。研究共鉴定出8个miRNA候选分子：其中5个被LPS上调并在ST36针刺后呈类对照模式（let-7d-5p、miR-103-3p、miR-107-3p、miR-451a、miR-7047-5p），3个被LPS下调并在ST36针刺后呈类对照模式（miR-5128、miR-16-2-3p、miR-101a-5p）。上述候选分子与原始基于倍数变化分析所鉴定的候选分子无重叠。对1,596个预测靶基因的功能富集分析显示，丝裂原活化蛋白激酶（MAPK）、磷脂酰肌醇3-激酶/蛋白激酶B（PI3K-Akt）及Wnt信号通路显著富集，基因本体（Gene Ontology）分析则突出显示其富集于神经元及突触相关过程。基于limma的重新分析揭示了一组被基于阈值方法所遗漏的ST36相关外泌体miRNA，为深入研究针刺在神经炎症中的外周效应提供了探索性候选分子指标。
论文解读：基于limma重分析鉴定ST36针刺相关外周外泌体miRNA候选分子指标
一、研究背景与科学问题
针刺的分子机制研究长期依赖脑组织分析，该方法具有侵入性且难以转化为临床监测手段。循环外泌体microRNA（miRNA）作为外周血中稳定存在、可反映全身炎症状态且无需组织采样的分子指标，为神经炎症相关中枢神经系统疾病的监测提供了极具吸引力的替代方案。然而，外泌体miRNA与针刺治疗效果之间的关联尚缺乏系统表征。足三里（ST36）针刺已被证实可调节神经-免疫通路，包括迷走-肾上腺抗炎轴，并在多种疾病中发挥系统性炎症调节作用。尽管已有研究报道针刺刺激与miRNA变化存在关联，但ST36针刺在神经炎症中外周miRNA特征性改变仍未得到完整定义。此外，有学者指出电针的抗炎机制涉及多系统、多水平、多靶点的调节，并非仅局限于迷走-肾上腺轴，这进一步凸显了开展全面分子表征的必要性。
研究人员前期曾在LPS诱导的抑郁模型中报道了ST36针刺治疗后血清外泌体miRNA的改变，为外周miRNA谱可能反映针刺诱导的免疫调节提供了初步证据。然而，该研究候选分子的鉴定依赖于倍数变化阈值和非配对t检验，缺乏跨独立数据集的正式统计验证，结果的稳健性受限。limma包中实现的经验贝叶斯框架通过将各探针的方差估计向共同先验收缩，可显著提升小样本微阵列研究的统计功效，这一优势对于原始数据每组仅3只小鼠（n = 3）的实验设计尤为重要。
二、研究内容与主要结论
研究人员对同一批微阵列数据应用limma重分析，旨在鉴定ST36治疗后呈现类对照表达模式的外泌体miRNA候选分子，并表征其预测功能靶点。研究共鉴定出8个类对照模式候选分子，其中5个被LPS上调并在ST36后呈类对照模式（let-7d-5p、miR-103-3p、miR-107-3p、miR-451a、miR-7047-5p），3个被LPS下调并在ST36后呈类对照模式（miR-5128、miR-16-2-3p、miR-101a-5p）。这些候选分子与原始基于倍数变化方法鉴定的结果完全无重叠，凸显了统计方法学在小样本miRNA研究中的关键作用。预测靶基因的功能富集分析显示，PI3K-Akt、MAPK及Wnt等信号通路显著富集，基因本体（GO）分析则突出富集于前脑发育、突触组装、神经元分化调控及突触结构或活性调控等神经元与突触相关过程。该研究的重要意义在于，通过统计重分析揭示了被传统阈值方法遗漏的ST36相关外泌体miRNA，为针刺外周效应的深入研究提供了探索性候选分子指标，并为外周血样本在针刺治疗反应监测中的转化应用奠定了基础。
三、关键技术方法
研究人员主要采用了以下关键技术方法：其一，数据来源于两套独立的Affymetrix GeneChip miRNA 4.0微阵列数据集，均来自C57BL/6雄性小鼠血清外泌体样本，第一套包含CON（对照）、ST36（仅ST36针刺）和LPS+ST36（LPS诱导炎症并接受ST36治疗）三组，每组3只；第二套包含CON和LPS两组，每组3只。两套数据集在同一实验室采用相同芯片平台、试剂批次和小鼠品系生成，但作为独立微阵列批次处理。其二，数据预处理采用Affymetrix Power Tools（APT）软件中的miRNA水平RMA+DABG算法进行背景扣除、标准化和探针集汇总，并以DABG p < 0.05作为检测标准进行探针过滤。其三，差异表达分析采用limma经验贝叶斯框架，配合Benjamini-Hochberg错误发现率（FDR）校正，使用无截距设计矩阵（?0 + group）配合contrasts.fit()和eBayes()（robust = TRUE, trend = FALSE）实现。其四，类对照模式候选分子定义采用三项标准同时满足的保守策略。其五，靶基因预测采用miRNAtap包（v1.46.0）的几何平均秩聚合方法，整合miRDB、TargetScan、PicTar、DIANA-microT和miRanda五个数据库，要求至少两个来源共识；通路富集分析采用clusterProfiler包（v4.20.0）进行KEGG通路和GO生物学过程富集。此外，还通过multiMiR包（v1.34.0）获取实验验证的miRNA-靶基因互作进行敏感性分析。
四、研究结果
4.1 差异表达分析结果
limma重分析在LPS与CON比较中，于352个表达探针中鉴定出123个FDR显著探针。在第一套数据集中，LPS+ST36与CON比较在423个表达探针中产生49个FDR显著探针，而ST36与CON比较未产生任何FDR显著探针。研究人员指出，鉴于每组仅3只小鼠的有限统计功效，该零结果不应被解读为ST36对健康对照miRNA表达无影响的证据，缺乏证据不等于证据缺失。
4.2 类对照模式候选分子
将三项标准的类对照模式定义应用于286个共同探针，共获得8个候选分子。5个被LPS上调并在ST36后呈类对照模式：mmu-let-7d-5p（LPS logFC = +2.69）、mmu-miR-103-3p（+1.77）、mmu-miR-107-3p（+1.34）、mmu-miR-451a（+1.34）和mmu-miR-7047-5p（+1.30）。3个被LPS下调并在ST36后呈类对照模式：mmu-miR-5128（LPS logFC = ?1.02）、mmu-miR-16-2-3p（?0.96）和mmu-miR-101a-5p（?0.71）。这8个候选分子与原始倍数变化分析鉴定的miRNA完全无重叠，表明统计方法可显著影响小样本miRNA研究的候选分子选择。
4.3 功能富集结果
靶基因预测共获得1,596个独特预测靶基因。KEGG通路富集鉴定出57条显著通路，最强富集见于PI3K-Akt信号通路（p.adj = 0.0004）、MAPK信号通路（p.adj = 0.0011）、Hippo信号通路（p.adj = 0.0016）和Wnt信号通路（p.adj = 0.0071）。轴突导向亦显著富集（p.adj = 0.0132）。GO分析揭示1,193个显著生物学过程条目，最强富集见于前脑发育、突触组装、神经元分化调控和突触结构或活性调控（均p.adj < 10??）。敏感性分析中，限制于实验验证的miRNA-靶基因互作（8,357个独特验证靶基因）后，57条预测靶基因KEGG通路中有55条仍显著富集，表明富集结果对限制于验证互作具有较好的稳健性。
五、讨论与结论
讨论部分指出，本研究所鉴定的8个类对照模式候选分子完全未被原始倍数变化阈值方法所发现，凸显了统计方法学在小样本miRNA研究中的实际重要性。limma经验贝叶斯框架通过将各探针方差估计向共同先验收缩来提高灵敏度，稳健经验贝叶斯收缩进一步限制了离群探针对该先验的影响，这在每组样本量较小导致个体方差估计不可靠时尤为关键。在候选分子中，let-7d-5p和miR-103-3p是公认的神经炎症信号调节因子，这些miRNA在ST36治疗后呈现的类对照模式与针刺调节神经-免疫通路的既有文献一致。MAPK和PI3K-Akt信号通路的富集与ST36刺激在LPS诱导神经炎症背景下与神经可塑性相关基因网络相联系的候选生物信息学推断通路相容，但该推断仍属假设生成性质，尚未经实验证实。GO条目在突触和神经元过程的汇聚同样与神经炎症机制一致，但不能独立确立该机制。值得注意的是，ST36在健康对照中未产生FDR显著miRNA变化，与LPS+ST36与CON比较中观察到的49个FDR显著探针形成对比，该对比与候选miRNA具有炎症背景依赖性的生物学活性相容，但不能确证。循环外泌体miRNA的转化意义在于其可及性，外周血采样在临床环境中可行，外泌体miRNA的稳定性使其适合作为监测工具。但需注意物种间保守性问题，mmu-miR-7047-5p和mmu-miR-5128为小鼠特异性或人类直系同源物表征有限，而let-7d-5p和miR-103-3p在哺乳动物中广泛保守。研究局限性包括：小样本量限制统计功效且无法表征个体变异；分析属相关性探索性质，未进行体外或体内功能验证；第一套数据集缺乏LPS单独组需跨数据集整合，残余批次效应无法完全排除；仅使用雄性小鼠，结果可能不适用于雌性；缺乏独立复制数据集，候选miRNA未经RT-qPCR等正交方法进行技术验证。
研究结论为：基于limma的外泌体miRNA微阵列数据重分析鉴定出8个ST36相关类对照模式候选分子，这些候选分子未被原始倍数变化方法检测到。其预测靶基因汇聚于神经炎症和神经可塑性相关通路，该模式与ST36针刺在所用LPS诱导神经炎症模型中的机制性联系相容，但并非证明。这些发现凸显了统计方法学在小样本miRNA研究中的重要性，并确定了值得进一步验证的探索性候选分子指标，为针刺研究中外周生物标志物的开发奠定了基础。
