GPNMB+疾病相关小胶质细胞状态在帕金森病中的整合转录组鉴定及SCFA–HDAC轴评估

《Neuroscience Informatics》：Integrated Transcriptomic Analysis Identifies a GPNMB+ Disease-Associated Microglial State in Parkinson's Disease and Evaluates the SCFA–HDAC Axis

　　

小胶质细胞活化诱导的神经炎症是帕金森病（PD）的关键过程，但相关细胞状态及肠道来源短链脂肪酸（SCFAs）的调控能力尚未完全阐明。本研究将黑质批量转录组分析与中脑单核转录组分析相结合，以考察PD相关小胶质细胞活化并评估所提出的SCFA–HDAC通路。通路水平分析显示多巴胺能和突触通路下调，同时神经炎症通路上调。单核分析将该信号定位于一个糖蛋白非转移性黑色素瘤蛋白B（GPNMB）阳性疾病相关小胶质细胞（DAM）样活化亚簇，该亚簇在本队列PD供者中持续富集。通过扰动连接性分析和分子对接，本计算替代性评估不支持HDAC抑制对该活化程序的一致逆转，核心靶点GPNMB的调控具有药物特异性而非一致性。上述发现鉴定了一个与PD神经炎症相关的候选GPNMB阳性活化小胶质细胞表型，并表明现有证据对所提出的SCFA–HDAC轴的支持有限。

研究背景方面，帕金森病（PD）为全球第二常见的神经退行性疾病，以黑质致密部多巴胺能神经元丢失及路易体病理形成为特征，运动和Non-motor表现共存构成严重疾病负担。近二十年来，神经炎症逐渐被认为不仅是伴随过程，而是黑质致密部细胞丢失的关键病理机制，其中中枢神经系统固有免疫效应细胞小胶质细胞受到广泛关注。慢性活化的小胶质细胞不仅分泌过量促炎因子，还与神经元损伤形成自我强化循环，但PD背景下慢性炎症发生的细胞状态及其介导机制尚未完全确定。在机制层面，异常积聚的α-突触核蛋白作为上游激活因素，通过模式识别受体激活小胶质细胞核因子κB（NF-κB）及NLRP3炎症小体信号通路，该过程已在人原代小胶质细胞中得到验证。全基因组关联研究（GWAS）提供了PD的遗传线索，其中7号染色体风险位点指向GPNMB基因；功能研究显示GPNMB通过与α-突触核蛋白相互作用增加疾病易感性，降低其表达可减少神经元对致病性α-突触核蛋白的摄取；蛋白水平上GPNMB在PD患者黑质中选择性上调，并在溶酶体应激时进一步升高。上述分子观察与小胶质细胞异质性理论相契合，疾病相关小胶质细胞（DAM）假说提示小胶质细胞由稳态向活化状态转变，伴随脂质代谢和吞噬相关基因改变，人样本中GPNMB与APOE、SPP1、TREM2等标记基因共同构成DAM特征。与此同时，肠道菌群及其衍生物短链脂肪酸（SCFAs）对神经炎症的潜在调控作用受到关注，荟萃分析显示PD患者肠道产SCFA细菌显著减少，粪便乙酸、丙酸和丁酸水平降低，但SCFAs在PD相关神经炎症中的作用存在持续争议，部分研究支持其保护多巴胺能神经元和小胶质细胞，另有研究报道其加重α-突触核蛋白病理和小胶质细胞炎症，且SCFAs通过组蛋白去乙酰化酶（HDAC）抑制发挥作用的确切过程尚未系统研究。既往研究集中于SCFAs对神经炎症的总体影响，未针对明确的小胶质细胞活化状态严格评估该机制，造成微生物代谢物与小胶质细胞靶点证据整合的空白。本研究针对两个具体空白：其一，尽管PD存在广泛小胶质细胞活化的证据，人中脑单细胞分辨率下疾病相关小胶质细胞的精确转录状态仍未完全界定；其二，尽管SCFAs被提出通过HDAC抑制调控神经炎症，该机制尚未针对分子定义的PD小胶质细胞活化特征进行评估。为此，本研究整合PD黑质批量转录组与中脑单核转录组，并采用扰动连接性分析和分子对接评估所提出的SCFA–HDAC通路，对SCFA功能采取中立而非预设的立场。研究主要贡献包括：分子对接证实SCFAs与HDAC8催化口袋的结构可结合性，但预测亲和力一致偏弱，为所提出的SCFA–HDAC轴提供了计算边界；扰动连接性分析显示HDAC抑制不能一致逆转该小胶质细胞活化程序，核心标记GPNMB的调控具有药物特异性而非一致性；利用单核转录组学在PD人中脑鉴定并量化了一个GPNMB⁺ DAM样小胶质细胞亚群，将组织水平神经炎症信号与离散的、分子定义的细胞状态相联系。研究人员开展的研究为整合分析，所得出主要结论为：PD中存在一个与神经炎症相关的候选GPNMB⁺活化小胶质细胞表型，而现有证据对所提出的SCFA–HDAC轴支持有限。其重要意义在于将组织水平神经炎症特征与可量化的分子定义小胶质细胞状态相连接，为PD神经炎症研究提供候选细胞与分子靶点，同时明确从间接计算证据推断代谢干预机制的方法学限度。该论文发表在《Neuroscience Informatics》。技术方法方面，本研究数据均来自公开数据库，包括PD黑质批量转录组数据集GSE7621（16例PD/9例对照）、PD中脑单核转录组数据集GSE157783（6例对照和5例特发性PD病例，源自Smaji?等）、LINCS L1000药物扰动数据集GSE92742和GSE70138中的HDAC抑制剂扰动特征，以及PDB 3F0R（人HDAC8晶体结构，保留催化Zn²⁺）；采用预排序基因集富集分析（GSEA）评估黑质通路水平改变，采用整合单细胞流程完成单核数据的质控、标准化、降维、供者感知Harmony整合、聚类和基于标记的注释，并以供者水平作为生物学重复进行比例检验和差异表达分析；以LINCS HDAC抑制剂转录扰动特征为机制替代，计算连接性评分并与随机零分布比较；分子对接评估乙酸、丙酸、丁酸及丙戊酸在HDAC8活性位点内的结合模式。研究结果方面，3.1 帕金森病黑质通路特征：多巴胺能下调伴随神经炎症通路系统性上调。通过黑质批量转录组通路水平分析，GSEA共发现275条通路显著富集，多巴胺能和突触通路下调，神经递质转运标准化富集评分（NES）为-2.29、多巴胺转运为-2.28、神经递质水平调控为-2.25，与黑质多巴胺能神经元选择性丢失的病理特征吻合；相反，67条神经炎症和先天免疫应答通路在PD组均上调，Toll样受体信号NES为+2.03、经NF-κB的肿瘤坏死因子α信号为+1.50，补体活化和白细胞介素6相关信号同时上调。结果提示PD黑质呈现多巴胺能神经元功能降低与组织水平炎症升高共存的转录组模式，但经多重检验校正后无单个基因保持显著差异调控，提示该数据集疾病信号在通路水平较单基因水平更稳健；由于批量转录组无法解析炎症信号的细胞来源，遂采用单核分析进行细胞状态定位。3.2 单核转录组将神经炎症定位至小胶质细胞：GPNMB⁺ DAM样活化状态在帕金森病中特异性扩增。中脑单核分析显示炎症相关信号与小胶质细胞状态相关，质控后共分析40,927个细胞核，经降维聚类鉴定出中脑主要细胞类型，包括小胶质细胞、少突胶质细胞、星形胶质细胞、神经元、少突胶质前体细胞、内皮细胞和周细胞，各细胞类型由特异性标记基因定义，小胶质细胞以CSF1R、P2RY12和CX3CR1高表达为代表。按此注释，PD队列小胶质细胞核比例较高（P = 0.0087），但由于该比例以发生显著神经元丢失的组织总核数为分母，其相对增加部分可能反映细胞组成变化而非小胶质细胞绝对增殖。小胶质细胞再聚类揭示以GPNMB、APOE和CD9表达为标志的DAM样活化表型，该亚群占总小胶质细胞的比例在PD队列中约高1.73倍（P = 0.0043），且本队列内全部5例PD供者的每供者比例均高于全部6例对照供者；鉴于样本量小（n = 5 vs 6），该供者水平分离应解释为队列内观察，需在更大的人中脑数据集中独立重复后方可推广。供者水平假体批量（pseudobulk）差异表达分析显示，核心标记GPNMB在PD患者小胶质细胞中显著上调，log2倍数变化为+1.30（P = 0.004），是唯一达到供者水平统计学显著的基因；APOE、TREM2等疾病标记基因及IL1B、TNFAIP3、NFKBIA等NF-κB轴基因呈一致上调趋势，P2RY12、CX3CR1等稳态标记基因呈一致下调。需注意单核测序无法捕获翻译后释放至细胞外的分泌型炎症趋化因子如CCL2和CXCL8，因此本研究的炎症表征基于细胞内活化标记转录本而非活化小胶质细胞完整分泌输出，该GPNMB⁺亚簇的实际炎症活性可能被低估。鉴于GPNMB作为PD风险相关基因已被验证，其在活化小胶质细胞亚群中的一致升高支持将其作为PD相关小胶质细胞状态的候选标记；由于单核测序对促炎趋化因子捕获不完全且供者数量有限，该发现应视为需独立验证的转录组状态标记。3.3 SCFA–HDAC通路对DAM活化程序的评估：无证据表明HDAC抑制逆转DAM特征。对所提出的SCFA–HDAC通路的评估未显示HDAC抑制能一致逆转疾病相关小胶质细胞活化程序的证据，在较小和较大药物扰动数据集中，HDAC抑制剂特征拮抗DAM样特征的能力与随机零分布无差异；基因表达方向检查进一步显示，HDAC抑制剂扰动后APOE、CD9、CD68、SPP1和TREM2等活化相关基因上调多于下调，核心标记GPNMB呈药物特异性和分歧性调控，伏立诺他（vorinostat）相关下调而曲古抑菌素A（trichostatin A）相关上调。值得注意的是，尽管丙戊酸在较小扰动数据集中对GPNMB呈下调趋势，该效应在较大数据集中未获统计学验证。上述结果提示现有间接扰动证据不支持SCFAs通过HDAC抑制逆转疾病相关小胶质细胞活化程序。3.4 SCFA与HDAC8结合的评估：结构可结合但结合弱。分子对接提示SCFA分子可占据HDAC8催化口袋但预测结合较弱，丁酸对接于催化锌离子附近并与过渡态稳定位点的Tyr306相互作用，SCFA相关HDAC抑制剂丙戊酸也占据催化腔靠近锌离子处；乙酸、丙酸、丁酸和丙戊酸的羧基距催化锌离子2.2–3.2 ?，预测结合能一致偏弱（-2.9至-4.7 kcal/mol）。由于这些差异处于对接函数典型评分误差范围内，且小羧酸片段接近可靠对接的尺寸下限，上述数值仅解释为与SCFAs已知毫摩尔级效力相符的弱而非特异性结合，不作为配体间定量亲和力排序。上述对接结果仅支持SCFAs与靶酶的结构可结合性，不等同于对疾病相关小胶质细胞活化程序的功能性逆转，与连接性分析的阴性结论相呼应。讨论部分总结如下：通过整合黑质批量与中脑单核转录组，本研究提示PD组织水平神经炎症信号与小胶质细胞重编程密切相关；在Smaji?等同一队列报道的泛胶质活化基础上，本研究进一步将弥漫性小胶质细胞活化解析为离散、可量化的GPNMB⁺ DAM样亚群，为解释PD病理的神经免疫组分提供锚点。批量转录组水平多巴胺能通路下调与神经炎症通路上调并存，提示黑质细胞组成与组织生理的双重改变；单核分辨率进一步提示炎症应答富集于小胶质细胞而非均匀分布于所有细胞类型，与小胶质细胞状态应解释为连续转录程序而非刚性分类的观点一致。该小胶质细胞亚群中GPNMB升高与PD患者黑质GPNMB选择性升高的报道一致，支持转录水平小胶质细胞活化与蛋白水平疾病应激反应之间的可能联系；但须承认本研究描述的GPNMB⁺ DAM样状态完全基于单一公开数据集的转录组水平定义，该计算特征与人组织直接病理验证之间差距显著，该转录状态是否对应PD黑质中形态学或功能上独特的小胶质细胞群体、小胶质细胞GPNMB蛋白丰度是否随疾病严重程度变化，均需在独立死后队列中经免疫组化或流式细胞术验证，因此本研究发现最好解释为关于疾病相关小胶质细胞状态的转录组假说而非已确认的病理实体。本研究与既往单细胞图谱工作的关系可调和且值得细究：既往人黑质单核图谱提示PD遗传风险主要富集于少突胶质细胞和多巴胺能神经元，小胶质细胞富集有限；本研究发现不与上述遗传水平结论矛盾，而是从疾病状态动力学角度予以补充，携带遗传风险的细胞类型与疾病进展中显示活化的细胞类型未必重合，DAM样小胶质细胞活化可能代表对黑质微环境退变的应答而非原发驱动，这与近期更大样本单细胞图谱报道的细胞类型特异性应答相呼应；Smaji?等描述了广泛胶质活化和PD特异性神经元状态，本研究再分析将部分活化定位至可测量的GPNMB⁺小胶质细胞亚群，提高了表征的细胞精度。另一关键发现是从中立计算视角评估了SCFAs通过HDAC抑制调控该过程的假说：药物扰动数据连接性分析显示HDAC抑制作为扰动类别不能逆转疾病相关小胶质细胞活化程序，且HDAC抑制剂扰动后多个活化标记上调多于下调、GPNMB调控具药物特异性；这些分歧效应可能反映抑制剂类别差异，因为脂肪酸羧酸盐在机制上不同于异羟肟酸和苯甲酰胺类HDAC抑制剂。分子对接强化了同样的审慎解释：尽管SCFAs可占据HDAC8催化口袋并接近催化锌离子，其弱预测结合亲和力与弱、毫摩尔级抑制剂而非强直接调节剂相符。综上，这些分析不足以支持SCFAs通过简单HDAC抑制机制逆转PD相关小胶质细胞活化的观点，与PD神经炎症表观遗传调控的更广泛复杂性一致。一个重要警示是连接性分析所用扰动特征主要来源于非神经细胞系（包括癌症来源和上皮细胞模型），而非原代小胶质细胞或中枢神经系统来源细胞；小胶质细胞在转录和表观遗传上与外周免疫细胞及非神经细胞系不同，具有由TGF-β和CSF1等CNS特异性信号维持的独特稳态特征，其对HDAC抑制的调控应答可能与非神经模型显著不同；这种跨组织推断引入无法通过计算校正的方向性偏倚，因此阴性连接性结果应解释为现有非神经替代数据中证据缺失，而非HDAC抑制在小胶质细胞中无效的直接证据。分子对接分析的另一考量是仅选择HDAC8作为代表靶点：选择HDAC8是因为具有分辨良好的晶体结构（PDB 3F0R）和清晰催化Zn²⁺位点，适合小羧酸配体对接；但I类HDAC包括HDAC1、HDAC2、HDAC3和HDAC8四个亚型，它们共享保守催化机制和锌螯合活性位点，但在调节亚基关联、组织表达和底物选择性上不同，HDAC1、HDAC2和HDAC3主要存在于核共阻遏复合物中且在免疫和神经细胞中表达高于HDAC8，因此SCFA介导的HDAC抑制若与小胶质细胞相关，可能优先通过HDAC1或HDAC3而非HDAC8发挥作用；本研究对接结果不应推广至整个I类家族，未来计算或生化工作应评估全部四个亚型的结合，以HDAC8为代表引入的结构不确定性限制了对接结论的范围。本研究评估结果凸显了关于SCFAs在PD中作用方向的长期争论：既往动物研究报道丁酸对多巴胺能损伤的保护作用，另有研究观察到SCFA混合物引起小胶质细胞活化和α-突触核蛋白病理，产SCFA细菌与疾病风险的孟德尔随机化分析联系亦不一致；本研究在细胞状态水平对SCFAs作用的批判性评估与SCFAs未必有益、其在PD中的影响很大程度上取决于浓度、暴露时间和特定细胞状态的观念相吻合；细菌产生的丁酸影响表观遗传改变和患者症状的发现从另一角度表明SCFAs的主要影响间接通过精细的表观遗传和代谢系统而非直接调控特定活化程序实现；SCFAs作为介导肠道菌群与小胶质细胞活动关系的信使，其调控机制分析需考虑全部作用模式而非依赖单一抗炎理论。研究局限包括：所有数据来自公开转录组和结构数据库，结论仍为关联性或预测性而非实验因果性；单核数据集供者数量有限（n = 5 PD vs 6对照），GPNMB⁺亚簇观察到的完全供者水平分离虽内部一致，但鉴于队列规模小应谨慎解释，需在更大的人中脑单核数据集中独立重复后方可视为PD的稳健判别特征；单核测序无法捕获CCL2和CXCL8等分泌型炎症趋化因子，GPNMB⁺亚簇炎症活性可能被低估；连接性分析所用扰动特征主要来自非神经细胞系，引入跨组织偏倚，限制其作为小胶质细胞HDAC调控机制替代的有效性；分子对接仅针对HDAC8，未捕获I类HDAC潜在亚型差异，结果不应推广至HDAC1、HDAC2或HDAC3；分子对接仅提供配体–口袋相容性的计算预测，不应解释为直接功能证据。未来研究应直接检验SCFAs是否在原代或iPSC来源小胶质细胞中调控GPNMB⁺小胶质细胞活化程序，评估SCFA在全部I类HDAC亚型中的结合，并在独立死后PD队列中验证GPNMB蛋白表达以确定其作为生物标志物或治疗靶点的价值。总体而言，本研究鉴定了一个GPNMB⁺疾病相关活化小胶质细胞状态作为PD神经炎症的候选特征，并界定了所提出的SCFA–HDAC通路当前的证据限度。研究结论翻译如下：通过整合PD患者黑质批量转录组与中脑单核转录组数据，本研究鉴定了一个与PD相关的GPNMB⁺小胶质细胞活化状态的富集。该发现将组织水平神经炎症特征与可量化的小胶质细胞状态相联系，为未来PD神经炎症研究提供了候选细胞和分子靶点。采用药物扰动连接性分析和分子对接的评估提示，目前不支持SCFAs通过HDAC通路作为该活化程序的一致抑制剂。关键靶点GPNMB的调控似乎具有药物特异性而非普遍性。上述结果既阐明了GPNMB⁺活化状态的发现价值，也明确了从间接计算证据推断代谢干预机制的方法学限度。