《Frontiers in Neurology》：Application of weighted gene co-expression network analysis to explore the targets of treadmill training in rats with spinal cord injury
编辑推荐：
背景：脊髓损伤（SCI）指脊髓任何损伤导致功能丧失。此类损伤常接受后续跑台训练（TM）治疗，本研究旨在揭示促进TM介导SCI后康复的机制。方法：识别差异表达基因（DEGs），并进行加权基因共表达网络分析（WGCNA）以识别TM相关模块。随后开展共表达网络构建、
背景：脊髓损伤（SCI）指脊髓任何损伤导致功能丧失。此类损伤常接受后续跑台训练（TM）治疗，本研究旨在揭示促进TM介导SCI后康复的机制。方法：识别差异表达基因（DEGs），并进行加权基因共表达网络分析（WGCNA）以识别TM相关模块。随后开展共表达网络构建、功能注释、蛋白质–蛋白质相互作用（PPI）网络及转录因子（TF）-miRNA-靶基因调控网络分析。研究人员建立SCI大鼠模型并使用逆转录定量聚合酶链反应（RT-qPCR）验证关键RNA。使用双荧光素酶报告基因实验检测miR-135a-5p与JAK2和STAT1之间预测的相互作用。结果：在SCI 8d + TM vs. SCI 8d和SCI 14d + TM vs. SCI 14d比较组中分别筛选出640和877个差异表达基因，并将其合并为DEGs，进一步缩小至四个与谷氨酸能突触和核糖体生物发生通路相关的TM相关模块。PPI网络评估确定v-Akt murine thymoma viral oncogene homolog 1（AKT1）、Janus kinase 2（JAK2）、tumor necrosis factor（TNF）和signal transducer and activator of transcription 1（STAT1）为枢纽节点。JAK2和STAT1也显示在模块b中存在相互作用。TF-miRNA-靶基因调控网络揭示STAT1和JAK2均为miR-135a的靶标。RT-qPCR结果证实TM后AKT1和STAT1上调，JAK2、TNF和miR-135a下调。双荧光素酶报告基因实验揭示miR-135a-5p与JAK2和STAT1的直接结合。结论：这些发现提示TM通过调控AKT1、JAK2、TNF、STAT1和miR-135a促进SCI恢复。
脊髓损伤（SCI，指脊髓任何损伤导致暂时或永久功能丧失）可由创伤或非创伤因素引起，临床表现从疼痛、麻木到瘫痪和失禁不等，并可伴压疮、肌萎缩、呼吸问题和感染等并发症；全球发病率估计为236至4,187/百万。运动可改善慢性SCI患者肌力与体能，TM可增强不完全性SCI后下肢运动功能和步态恢复，动物模型中适量TM也促进功能恢复、神经肌肉再生和骨骼肌可塑性；运动还可通过调控靶向PTEN/mTOR等通路的miRNA影响损伤脊髓可塑性，ERK1/2信号及BDNF/KCC2等亦参与TM效应。然而，TM介导SCI大鼠恢复的潜在机制仍 largely unknown。2015年Baligand et al.在NYU-MASCIS中度挫伤SCI大鼠模型中进行TM，并对比目鱼肌进行表达谱分析，识别TGFβ/SMAD3和BMP通路为SCI及TM后肌肉重建的关键驱动因素，但未深入探索TM影响SCI大鼠所涉及的基因。为此，研究人员利用其GSE45550微阵列数据集开展系统生物信息学分析，并用RT-qPCR验证关键标志物，以揭示TM介导SCI后康复的机制，为理解SCI后创伤后医疗的深远影响提供依据。研究人员得出结论：TM通过调控AKT1、JAK2、TNF、STAT1和miR-135a促进SCI恢复，并提示谷氨酸能突触和核糖体生物发生通路参与其中。
主要技术方法（≤250字）：使用GEO来源的GSE45550数据集，含24个大鼠比目鱼肌样本（SCI 8d、SCI 14d、SCI 8d + TM、SCI 14d + TM各n=6，GPL1355 Affymetrix Rat Genome 230 2.0 Array）；以RMA预处理、limma筛选DEGs（p<0.05，未校正，无fold-change阈值）；以WGCNA识别TM相关模块；以DAVID进行GO注释，以pre-ranked GSEA评估谷氨酸能突触和核糖体生物发生通路；以STRING和Cytoscape构建PPI网络并用MCODE识别模块，结合CTD与GeneCards交集筛选关键基因；以WebGestal ORA构建TF-miRNA-靶基因网络；建立20只成年雄性Sprague Dawley大鼠SCI模型（T9钳夹损伤，分组含SCI 8d、SCI 14d、SCI 8d + TM、SCI 14d + TM及未损伤/未损伤+TM），用RT-qPCR验证关键RNA；用双荧光素酶报告基因实验验证miR-135a-5p与JAK2、STAT1 3′UTR（3′非翻译区）结合。
研究结果：
Differential expression analysis（差异表达分析）：在SCI 8d + TM vs. SCI 8d和SCI 14d + TM vs. SCI 14d中分别获得640个（428上调、212下调）和877个（439上调、438下调）DEGs，并将两组DEGs并集用于后续WGCNA。
WGCNA：确定软阈值power=16，构建无标度网络并将DEGs划分为30个模块，识别7个与TM相关且|r|>0.8的模块，重点为MElightcyan1、MEbisque4、MEdarkseagreen4和MEgreen。
Co-expression network construction and functional annotation（共表达网络构建与功能注释）：对四个模块节点构建共表达网络；GO注释显示MElightcyan1与突触传递和谷氨酸通路相关，MEbisque4与肌动蛋白丝聚合负调控相关，MEdarkseagreen4与钙离子跨膜转运相关，MEgreen与DNA模板转录负调控相关。合并节点得492个基因（308上调、190下调）；上调基因主要与离子跨膜转运和谷氨酸能突触相关，下调基因主要与神经管闭合和真核生物核糖体生物发生相关；pre-ranked GSEA显示谷氨酸能突触基因集在第8天富集（NES=1.67，adjusted p=0.0108），第14天谷氨酸能突触及两个时间点核糖体生物发生未显著富集。
PPI network construction and key gene identification（PPI网络构建与关键基因识别）：PPI网络含249个节点和476条边；AKT1度最高（degree=35），JAK2（degree=18）、TNF（degree=14）、STAT1（degree=14）为高度节点；识别模块a（score=7.714，8节点27边）和模块b（score=5.419，8节点19边），JAK2与STAT1在模块b中相互作用；CTD、GeneCards与数据集交集得到AKT1、JAK2、TNF、STAT1四个关键基因；14天SCI + TM组中JAK2和TNF表达显著降低，AKT1和STAT1显著升高（p<0.05）。
Construction of the TF-miRNA-target regulatory network（TF-miRNA-靶基因调控网络构建）：识别14个TFs和7个miRNAs，产生310个TF-miRNA-靶基因对，网络显示STAT1和JAK2均为miR-135a的靶标。
RT-qPCR：验证显示AKT1和STAT1在SCI 8d + TM和SCI 14d + TM组显著上调，JAK2、TNF和miR-135a较相应SCI 8d和SCI 14d组显著下调（p<0.05）。
Dual-luciferase reporter assays（双荧光素酶报告基因实验）：miR-135a-5p mimic显著降低JAK2和STAT1野生型3′UTR报告基因活性（p<0.05），突变型无显著变化，支持miR-135a-5p与JAK2和STAT1 3′UTR预测位点直接结合。
讨论部分总结：原文将结果与其他研究相联系，指出TM可经PI3K/Akt通路上调神经营养因子、抑制损伤组织凋亡并促进运动功能恢复，Akt通路在脊髓中经内皮NO通路发挥保护作用，AKT1与肌母细胞分化相关；TNFα可诱导少突胶质细胞和神经元凋亡，并参与炎性痛和脊髓突触可塑性，etanercept可减轻SCI组织损伤、促进髓鞘再生和运动功能；抑制STAT1可减少组织损伤和运动障碍，JAK2/STAT1在模块b中相互作用提示其参与SCI病理结局；miR-135a降低与GATA-3升高相关，GATA-3和miR-135a可能支持运动诱导的免疫抑制，Wang et al.报道miR-135a-5p通过靶向SP1和ROCK促进SCI功能恢复，而本研究提示STAT1和JAK2亦为miR-135a靶标。核糖体生物发生影响细胞翻译能力和肌肉生长维持，谷氨酸能神经元可支配肌纤维并恢复运动，因此原文提示TM可能通过谷氨酸能突触和核糖体生物发生通路促进SCI康复。治疗视角上，AKT1相关信号激活可能增强神经元存活和轴突再生，抑制过度JAK2/STAT信号可能减轻神经炎症和神经元凋亡，TNF可被etanercept等中和生物制剂靶向，但证据主要限于临床前，仍需实验验证。局限性包括：未校正P<0.05筛选DEGs增加假阳性风险；发现数据集和RT-qPCR样本来自比目鱼肌，不能单独确立脊髓内分子修复；验证主要限于RNA表达和双荧光素酶报告基因实验，蛋白丰度、磷酸化、细胞特异性定位及AKT1、JAK2、TNF、STAT1、miR-135a-5p的功能效应需敲低、过表达或药理学干预研究；仅评估第8天和第14天及雄性大鼠，行为学结果未与分子测量相关；未整合蛋白质组、磷酸化蛋白质组、代谢组或单细胞数据。
研究结论（翻译）：总之，本研究识别AKT1、JAK2、TNF、STAT1和miR-135a为与SCI后跑台训练反应相关的候选分子因素。综合来看，这些发现还提示TM可能通过谷氨酸能突触和核糖体生物发生通路的变化发挥SCI康复作用。
