脊髓损伤中细胞类型特异性的m1A动态变化与小胶质细胞表型转化及神经元代谢适应相关

《PLOS Computational Biology》：Cell-type-specific m1A dynamics are associated with microglial phenotypic transformation and neuronal metabolic adaptation during spinal cord injury

　　

m1A（N1-methyladenosine，N1-甲基腺苷）是一种重要的表观遗传机制，可调控包括脊髓损伤（SCI）在内的多种疾病的发生与发展。为探究SCI后m1A水平的整体变化，研究人员分析了SCI样本的转录组测序数据，并基于m1A调控因子水平赋予m1A评分。本研究中，m1A评分是根据m1A调控基因（writers/erasers/readers）表达量推算的替代指标，并不直接测量RNA m1A修饰水平。结果显示，SCI后第1天m1A评分升高，随后下降，至第3天和第7天低于基线。进一步分析揭示，小胶质细胞与神经元是m1A评分变化最显著的两类细胞：小胶质细胞中m1A评分各时间点均降低，神经元中则均升高。通过拟时序分析与功能富集分析，m1A评分可能关联小胶质细胞表型转换与神经元能量代谢，并在体内外实验中获验证。总之，本研究揭示了SCI后bulk与单细胞水平m1A特征性变化，提示其与神经元功能潜在联系，并为探索m1A相关调控因子作为SCI治疗靶点提供理论依据。

研究背景与立项依据：脊髓损伤（spinal cord injury，SCI）是一类导致感觉、运动及自主神经功能丧失的严重神经系统疾病，全球患者基数大且急性期的药物或手术疗效缺乏证据支持。神经元再生能力弱、原发与继发氧化应激损伤，以及小胶质细胞由稳态向炎性表型极化加剧神经炎症，是阻碍修复的核心病理环节。RNA表观修饰中m⁶A研究较多，而m¹A（N1-methyladenosine，N1-甲基腺苷）在SCI中的动态变化与细胞类型特异性功能长期空白。m¹A由写入蛋白（writers，如TRMT10C、TRMT6-TRMT61A）、擦除蛋白（erasers，如ALKBH3、ALKBH1、FTO）与阅读蛋白（readers，YTH家族）共同调控，既往证据多集中于肿瘤与少数神经退行病，直接在SCI组织内测量RNA m¹A修饰难度较大，因此研究人员以调控因子表达构建m¹A评分为替代指标，在bulk RNA-seq、scRNA-seq与snRNA-seq（单细胞细胞核测序）多层级上解析其动力学。

主要技术方法概览：研究人员整合已发表的成年小鼠SCI单细胞资源（GSE162610、GSE172167，含Sham及1、3、7 dpi样本），以Seurat完成质控与批次校正；用UCell/GSVA将10个m¹A调控因子（Ythdc1、Ythdf1/2/3、Trmt10c、Trmt61a、Trmt6、Alkbh3、Alkbh1、Fto）打包为基因集计算细胞m¹A评分；辅以pySCENIC推断转录因子活性、Monocle2拟时序追踪小胶质细胞状态转换、SuperCell构建metacell降低稀疏噪声；并行开展体内小鼠T9撞击SCI模型与体外HT22神经元OGD/R（氧糖剥夺/复氧）验证，以qPCR检测呼吸链亚基（ND1、COX1、ATP6）及ATP定量。

研究结果

2.1 SCI后bulk水平m¹A评分动态及其与能量代谢通路显著相关：bulk RNA-seq显示除Ythdf3外，各调控因子与m¹A评分在1 dpi上升、3/7 dpi回落低于基线；GSVA/GSEA表明m¹A评分正相关线粒体电子传递链、ATP代谢与氧化磷酸化，提示亚急性期评分走低或对应能量代谢紊乱。

2.2 SCI后单细胞水平m¹A评分的细胞异质性：15类细胞聚类后，神经元m¹A评分伤后持续上行，小胶质细胞、巨噬细胞、周细胞、OPC等持续下行，中性粒/少突胶质/室管膜细胞呈先升后降；证明m¹A调控具明显细胞类型偏好。

2.3 髓系细胞转录网络失调与m¹A评分密切相关：髓系中仅小胶质细胞在数量增多的同时评分显著下降；pyscenic显示Etv6等正关转录因子活性减弱、writer表达受抑；Monocle2拟时序重叠小胶质从homeostatic→inflammatory/interferon→dividing→homeostatic转换，migrating亚型持续高驻；体内qPCR以Tmem119/Cxcl2/Pclaf/Adam8/Isg15标记验证五亚群时序切换与评分下行共现。

2.4 SCI后神经元m¹A评分显著差异：因神经元胞体大、常规scRNA-seq捕获偏倚，研究人员用snRNA-seq补位，发现神经元各时间点m¹A评分单向攀升，其余细胞变化不显著，说明伤后神经元转录态与高m¹A调控活性绑定。

2.5 m¹A评分与伤后神经元能量代谢相关：神经元差异基因GO显示在氧化磷酸化、ATP代谢、PI3K信号抑制背景下，metacell级别UCell评分相关性分析指出m¹A评分负调氧化磷酸化/线粒体ATP合成/有氧呼吸，正调mRNA代谢与神经元黏附，并提炼出候选关键基因集。

2.6 脊髓损伤后神经元能量代谢显著下降：snRNA-seq见呼吸链亚基（ND1/COX1/ATP6）转录上调，HT22 OGD 12h后复氧0.5h亚基mRNA继续升高、24h仍高于对照；但ATP检测（OGD 1h）显示缺血期ATP骤降、复氧0.5h回升、24h回基线，解释为单位时间内呼吸链亚基代偿性高表达对抗低氧、而实时能量产出仍受抑。

讨论与结论翻译：研究人员在讨论中指出，SCI微环境具炎症、星形胶质增生、神经元变性等多重动态，bulk与单细胞互补必要；m¹A评分作为调控基因表达代理而非直接修饰量度，须谨慎引申。结果发现bulk层1 dpi短暂升高后走低，单细胞层小胶质评分持续降、神经元持续升，前者共现于homeostatic→inflammatory→dividing→homeostatic表型轨线，后者负关联氧化磷酸化并伴缺血期ATP下降与呼吸链亚基代偿。小胶质m¹A研究极少，神经元中m¹A与OGD应激转录已有外围证据但因果未明。作者自陈局限：HT22非原代脊髓神经元、体内外均未敲除/过表达m¹A酶故不能断言因果，仅提供关联假设。结论段原意浓缩为：本研究在bulk、scRNA-seq及snRNA-seq三层揭示SCI后m¹A特征时空变化，锁定小胶质细胞表型转换与神经元能量代谢两类最显著生物学过程，支持后续以TRMT6/TRMT61A、ALKBH3等m¹A调控因子为干涉节点的机制研究，论文发表于《PLOS Computational Biology》。